FDTD 튜토리얼: 실리콘 Mach–Zehnder 변조기

1. 개요: 무엇을 시뮬레이션할 것인가

Mach–Zehnder modulator (MZM)는 전기적 제어 신호를 광학 진폭(또는 세기) 변화로 변환하는 간섭형 광자 부품입니다. 빛은 두 도파관 팔로 분기되고, 상대적인 위상차를 축적한 뒤 다시 결합됩니다. 출력 전력은 두 팔 사이의 위상 관계에 의존합니다: 두 팔이 동위상으로 재결합하면 출력은 밝고, 역위상으로 재결합하면 출력은 억제됩니다.

실리콘 포토닉스에서 MZM은 고속 데이터 통신, 코히런트 광학, 광 상호연결, 집적 마이크로파 포토닉스에 널리 사용됩니다. 실제 장치는 일반적으로 플랫폼에서 가능한 전기광학 메커니즘을 사용하여 위상 제어를 구현합니다 (실리콘의 경우 이는 종종 PN 접합, MOS 커패시터를 통한 캐리어 분산 또는 열광학 조정을 위한 히터입니다). 구동 메커니즘과 무관하게, 장치의 광학 코어는 동일합니다: 분기기, 두 개의 유도 팔, 그리고 재결합기입니다.

이 튜토리얼에서는 대표적인 실리콘 Mach–Zehnder 변조기 구조의 FDTD 시뮬레이션을 실행합니다. 여기서의 목표는 완전한 전기 공동 시뮬레이션을 구축하는 것이 아니라, 시간 영역 장 전파를 사용하여 유도된 에너지가 어떻게 분기되고, 두 팔을 따라 이동하며, 출력에서 다시 결합하는지에 대한 직관을 기르는 것입니다. 전력 밀도 스냅샷을 확인하고 검출기 시간 추적을 해석할 것이며, 여기에는 팔에 위치한 검출기들과 두 팔이 다시 합쳐지는 더 하류에 위치한 검출기 사이의 전파 지연이 포함됩니다.

수치적 관점에서 이 예제는 의도적으로 “무거운” 예제입니다: 시뮬레이션 도메인이 비교적 크고 디스크에 상당한 양의 데이터를 덤프합니다. 따라서 이는 전체 FDTD 스택의 좋은 스트레스 테스트가 됩니다: 기하 정의, 소스 주입, CPML 흡수 경계, 시간 스텝 진행, 검출기 추출, 그리고 출력 파이프라인 (스냅샷과 검출기 뷰어)을 포함합니다.

2. 새 시뮬레이션 만들기

새 시뮬레이션 버튼을 클릭하여 New simulation 창을 엽니다. 이 창에서 FDTD examples 카테고리를 선택한 다음 Silicon Mach–Zehnder modulator 데모를 더블클릭합니다 (목록의 맨 아래에 있습니다). 메시지가 나타나면 시뮬레이션을 로컬 디스크에 저장합니다.

이 예제는 실행 중 많은 파일을 기록합니다 (스냅샷, 검출기 출력, 중간 데이터). 따라서 네트워크 드라이브나 가상화/클라우드 동기화 디렉터리(예: OneDrive)에 저장하지 마십시오. 로컬 SSD가 이상적입니다. 느린 저장소에 저장하면 실행 시간이 극적으로 증가하거나, 출력 기록 중 인터페이스가 느려질 수 있습니다.

3. 메인 창에서 방향 잡기

장치는 Device structure 탭의 3D 보기에 나타납니다 (??). 두 개의 간섭계 팔(racetrack 형태의 도파관 경로)과 구조를 따라 배치된 색이 다른 검출기 객체(보라색, 녹색, 빨간색)를 식별할 수 있어야 합니다. 정확한 색은 중요하지 않지만, 기하 구조에서 신호가 어디서 모니터링되는지 쉽게 볼 수 있게 해 줍니다.

이 예제에서 주로 사용할 것은 다음과 같습니다:

• FDTD 계산을 시작하는 Run simulation (▶) 또는 F9.
• 스냅샷과 검출기 디렉터리/파일을 여는 Output 탭.
• 시간을 따라 전력 밀도가 어떻게 변화하는지 관찰하는 Snapshot viewer.
• 각 검출기의 power.csv 시간 추적을 플롯하는 Detector power viewer.

4. 시뮬레이션 실행

F9를 누르거나 파란색 재생 버튼(▶)을 클릭하여 실행을 시작합니다. 이 시뮬레이션은 비교적 크며, 격자가 크고 스냅샷 및 검출기 데이터를 디스크에 덤프하기 때문에 완료하는 데 상당한 시간이 걸립니다.

성능이 괜찮은 노트북에서는 일반적인 실행 시간이 약 10–15분입니다. 실행이 진행되는 동안 Output 탭을 열어 파일이 나타나는 것을 볼 수 있습니다. 특히 snapshots/ 디렉터리가 보이고, 시뮬레이션이 진행됨에 따라 검출기 출력(Detector 0, Detector 1, Detector 2)이 나타납니다.

5. 전력 밀도 스냅샷 보기

Output 탭에서 snapshots/ 폴더를 엽니다 (??) 그러면 snapshot viewer가 실행됩니다. 대표적인 snapshot-viewer 창은 ??에 나와 있습니다. 파란색 플러스 버튼을 사용하여 powerdensity.csv를 뷰어에 추가합니다. 그 다음 슬라이더를 사용하여 시간을 따라 이동할 수 있습니다.

이 시뮬레이션을 보는 가장 유용한 방법은 장치 내부를 전파하는 유도 에너지를 따라가는 것입니다. 초기 시간에는 소스가 입력 도파관으로 주입되고 분기기 영역에 도달하는 것을 볼 수 있습니다. 시간이 진행됨에 따라 전력은 두 팔로 나뉩니다. 이후 스냅샷에서는 전력이 racetrack 팔을 따라 계속 전파되는 것이 보입니다. 더 나중에는 두 팔이 다시 결합하고, 에너지가 출력 섹션으로 라우팅되는 것을 볼 수 있습니다.

대표적인 스냅샷은 ??– ??에 나와 있습니다. 장치 매개변수를 바꾸지 않더라도, 이 스냅샷들은 이것이 간섭 구조라는 점을 분명하게 보여줍니다: 에너지가 분기되고, 물리적으로 분리된 두 경로를 따라가며, 그다음 다시 합쳐집니다.

6. 검출기 출력 보기

시뮬레이션이 실행되는 동안(그리고 완료 후에도) Output 탭에 검출기 출력이 나타납니다 (??). 각 검출기 항목은 작은 CCD 스타일 아이콘처럼 보입니다. 검출기를 열려면 이를 더블클릭한 다음, 해당 검출기 디렉터리 안의 power.csv를 더블클릭하여 detector power viewer를 엽니다.

이 예제는 세 개의 검출기를 사용합니다. 두 개의 검출기는 두 간섭계 팔을 모니터링하고(그래서 유도 에너지가 각 팔 모니터에 도달하는 시간을 비교할 수 있음), 세 번째 검출기는 두 팔이 재결합하는 최종 출력 단계의 더 하류에 배치됩니다. 세 개의 검출기 플롯은 ??, ??, 그리고 ??에 나와 있습니다.

이 실행에서는 두 팔 검출기에서 전력이 대략 \(3 \times 10^{-13}\,\mathrm{s}\) 부근에서 상승하기 시작하는 반면, 더 하류에 있는 검출기는 (대략 그 두 배 정도의 시간에) 더 늦게 시작하는 것을 볼 수 있습니다. 이는 두 팔이 이동하고 재결합한 뒤의 전파 경로에 더 하류에 배치되었기 때문입니다. 이것은 구조를 통한 전파 지연을 직접적이고 직관적으로 보여줍니다: 팔 모니터가 먼저 응답하고, 최종 단계 모니터는 유도 에너지가 물리적으로 더 멀리 이동해야 하므로 나중에 응답합니다.

추적이 완전히 매끄럽지는 않은데, 이는 우리가 강한 시간 평균을 하지 않기 때문입니다; 대신 모니터를 통과하는 시간 의존 유도 성분을 볼 수 있습니다. 이러한 플롯을 “순간적인” 전력으로 해석하는 것은 보통 바람직하지 않으며; 이 튜토리얼에서 더 유용한 해석은 도달 시간과 유도 에너지가 계속 전파되고 간섭함에 따라 신호 외피가 어떻게 진화하는가입니다.