실린더 구조를 활용한 메타표면에서 BIC(연속체 속박 상태)를 결합하여 2차 고조파를 증대시키는 연구 결과를 재현했다. 이 기술은 초표면 분야에서 주목받고 있으며, 실린더 배열만으로도 기존 구조 대비 2차 고조파 강도를 10배까지 높일 수 있음을 시뮬레이션을 통해 확인했다.

BIC는 특정 구조 내에서 빛이 빠져나가지 못하도록 가두는 현상을 의미한다. COMSOL 모델링 시 핵심은 광장이 완벽하게 국소화되는 기하학적 파라미터를 찾는 것이다. 아래는 구조 파라미터 설정 예시다:

period = 820e-9  # 주기 설정: 운동량 정합 조건 충족 필수
radius = 0.28 * period  # 반경 계수: 미세 조정 필요
height = 350e-9  # 높이: 모드 결합에 영향
material = 'Si'  # 실리콘: χ² 값은 낮지만 구조적 이점 활용

여기서 radius의 0.28 계수는 매우 민감하다. 0.01만 변경해도 Q값이 급감할 수 있다. 실수로 정수를 사용하면 시뮬레이션이 수렴하지 않을 위험이 크다.

모드 여기에는 파수 스캔 방식을 적용했다. 자동화 스크립트로 k-공간을 탐색하며 BIC 지점을 찾는다:

% 파수 스캔 파라미터
k_min = 0.4 * (2 * pi / period);
k_max = 0.6 * (2 * pi / period);
steps = 50; % 30단계 초과 시 시스템 불안정 가능성
for k = linspace(k_min, k_max, steps)
    update_model_parameters(k); % 최적화 함수 사용
    solve_eigenfrequency();
    if abs(imag(freq)) < 1e-3 % Q값이 최대가 되는 지점 식별
        break;
    end
end

1e-3 임계값은 신중히 선택해야 한다. 너무 크면 BIC 지점을 놓치고, 너무 작으면 계산 시간이 과도하게 증가한다. 과도한 설정으로 인해 서버가 과부하에 걸린 사례도 있었다.

2차 고조파 증대의 핵심은 전기장 국소화다. 시뮬레이션 데이터에서 E장 성분을 추출하여 분석한다:

# 전기장 증대 인자 추출
E_norm = sqrt(Ex**2 + Ey**2 + Ez**2)
enhancement = max(E_norm)**4  # 비선형 효과는 4제곱 의존성
print(f"전기장 증대: {enhancement:.1e}배")

4제곱을 사용하는 이유는 2차 고조파 강도가 기저 주파수 장의 4제곱에 비례하기 때문이다. 이 관계를 간과하면 결과 해석에 오류가 발생할 수 있다.

실전 팁: 시뮬레이션에서 자동 메시 생성은 피해야 한다. 수동으로 경계층 메시를 설정하여 표면 장 변화를 정확히 포착하면서도 계산 부담을 최소화한다:

경계층 수: 5
성장률: 1.2
최소 요소 크기: 2nm

자동 메시를 사용하면 전기장 분포가 비현실적으로 나타나 해석이 어려워질 수 있다. 예를 들어, 튜닝이 잘못된 경우 전기장 패턴이 추상화처럼 보여 연구실에서 논의가 길어진 적이 있다.

재현 결과는 이론 예측과 일치했다. 특히 1100nm 파장에서 2차 고조파 신호가 급격히 증가하는 현상을 관찰했다. BIC 덕분에 실리콘과 같은 약한 비선형 재료도 강력한 비선형 효과를 발휘할 수 있음을 입증했다.