도 1은 본 발명에 따른 검사 방법을 개략적으로 도시한 순서도이다. 이와 같 은 방법으로 주기 구조물의 구조 및 내부 구성 성분을 비파괴적으로 검사할 수 있다. 특히, 주기 구조물에 자연적으로 형성되는 산화막이나 의도적으로 형성된 표면 코팅층 등을 포함한 미세 구조를 정밀하게 검사할 수 있다.
주기 구조물의 측정
도 2 및 도 3은 주기 구조물의 검사 장치를 개략적으로 도시한 모식도이다. 도 2를 참조하면, 검사 장치는 광원(100), 검출기(110) 및 프로세서(120)를 포함한다. 기판(130)에 검사 대상인 주기 구조물(200)이 위치할 경우, 광원(100)은 주기 구조물(200)에 특정 파장 또는 여러 파장을 갖는 빛을 조사한다. 주기 구조물(200)에 입사한 빛 중 일부는 투과하고, 일부는 반사한다. 반사된 빛은 검출기(110)에서 검출되고, 프로세서(120)에서는 검출기(110)에서 측정된 반사파의 반사율이 산출된다. 투과된 빛 또한 검출기(110)에서 검출되고, 프로세서(120)에서는 검출기(110)에서 측정된 투과파의 투과율이 산출된다.
상기 검사 장치는 도 3에 도시된 바와 같이, 편광기(140)를 더 구비할 수 있다. 이 경우, 광원(100)에서 발생한 빛은 편광기(140)를 거쳐 TE 모드의 빛 또는 TM 모드의 빛으로 편광되어 주기 구조물(200)에 입사된다.
주기 구조물(200)에 빛을 입사시키면, 입사한 빛은 반사되는 빛과 투과되는 빛으로 나뉘어진다. 본 발명에서는 빛의 반사 및 투과에 있어서 가장 근본적인 2개의 편광 상태, 즉 TE 모드 및 TM 모드의 반사율 및 투과율을 산출하여 주기 구조물의 비파괴 검사를 수행한다.
예를 들어, 주기 구조물에 빛을 입사하여 측정된 반사율 및 투과율과 관계된 물리량은 TE 모드 전기장의 입사파에 대한 반사파 및 투과파의 진폭 또는 위상과 관계된 물리량 및 TM 모드 자기장의 입사파에 대한 반사파 및 투과파의 진폭 또는 위상과 관련된 물리량의 조합으로 이해될 수 있다.
가상 주기 구조물
도 4는 가상 주기 구조물의 기하학적 구조를 도시한 사시도이고, 도 5는 도 4의 가상 주기 구조물의 단면을 복수의 층으로 나눈 상태를 도시한 것이다.
상기 주기 구조물(200a)을 예를 들어 1차원, 2차원, 또는 3차원의 주기구조의 형태를 갖는 반도체 소자라고 가정하면, 상기 가상 주기 구조물(200a)은 각 층마다 실리콘 등의 해당 물질로 이루어진 물질부분(n
1 부터 n
L까지)과 공기층 등 입사부의 물질부분
의 두 가지 물질이 수평적으로 주기적인 구조를 갖는다.
그러나 반도체 공정 등의 실제 환경에서는, 주기 구조물을 완벽한 진공 상태에서 제조하지 못하기 때문에 실제 주기 구조물의 표면에는 공기 또는 수분과의 접촉에 의해 산화막이 형성된다. 또한, 공정단계에서 주기 구조물의 표면에 의도적인 코팅층을 형성하거나, 주기 구조물 표면에 거칠기층(roughness layer)이 존재하여 실제의 기하학적 형태를 추정하는데 도 4의 주기 구조물은 한계가 있다.
도 6은 본 발명에 따른 가상 주기 구조물(200b)의 기하학적 구조를 도시한 것이다. 이러한 가상 주기 구조물(200b)은 표면에 산화막 등의 표면층(210)이 형성되어 있으며, 그 단면을 복수의 층으로 분할하면 도 7에 도시한 바와 같이 적어도 3 개의 물질이 주기적으로 반복되는 구조를 갖게 된다. 상기 가상 주기 구조물(200b)은 골(groove) 영역에 해당하는 제3물질을 사이에 두고 반복적인 주기로 형성된 마루(ridge) 영역을 포함하며, 이 마루 영역은 제1물질로 구성된 중심부와 제2물질로 구성되며 상기 제2물질은 중심부의 외면에 형성되는 표면층을 포함한다. 도 7에서
는 층
l 에서 마루(ridge) 영역(제1물질), 골(groove) 영역(제3물질), 표면층 영역(제2물질)을 각각 나타낸다.
상기 표면층(제2물질)은 산화막 또는 코팅층일 수 있으며, 경우에 따라서는 주기 구조물 표면의 거칠기층일 수도 있다. 상기 골 영역에 해당하는 제3물질은 기상, 액상, 또는 고상일 수 있다.
예를 들어, 상기 가상 주기 구조물(200b)을 반도체 소자로 가정하면, 직사각형 단면의 구조의 복수 층(1 내지 L)에서, 최상부층(제1층)을 제외하고, 중간의 각 층마다 실리콘 등의 물질로 이루어진 제1물질, 산화막 또는 코팅층으로서 존재하는 표면의 제2물질 및 공기층이나 액상 또는 고상형태의 제3물질이 수평적으로 주기적인 반복 구조를 가지게 된다.
이와 같이 산화막이나 코팅층, 또는 표면의 거칠기층 등의 표면층(210)을 고려하여 가상 주기 구조물(200b)을 설정하게 되면 가상 구조물의 반사율 및 투과율 을 실제 주기 구조물과 더욱 근사하게 산출할 수 있다. 따라서, 나노 수준의 미세 주기 구조물의 구조 및 성분을 매우 정확하게 측정할 수 있다. 구체적으로는 주기 구조의 기하학적인 외관 형태 및 내부 구성 성분, 그 아래에 존재하는 박막 구조들의 두께까지 비교 분석할 수 있다.
가상 주기 구조물의 물리량 산출
본 발명에서는 가상 주기 구조물의 물리량을 산출하기 위해서 가상 주기 구조물의 한 형태를 가정하여 영차 구조물(zero-th order periodic structure)로 정의하고, 이 영차 구조물을 섭동영역에서 기하학적 또는 물리적 변화를 가한 섭동된 주기 구조물(perturbed periodic structure)을 얻는다.
영차 구조물에 대한 빛의 반사율 및 투과율, 섭동영역에서의 영차 구조물의 그린함수(Green's function)를 계산하고, 리프만-슈빙거 적분방정식 (Lippmann-Schwinger Equation)을 반복법(iteration)으로 풀어서, 이들 결과로부터 섭동된 주기 구조물의 반사율, 투과율, 또는 반사율 및 투과율에 관계된 물리량을 구한다.
영차 구조물에 대한 빛의 반사율 및 투과율은 맥스웰 방정식에 주기조건을 더한 계산을 통하여 추정할 수 있다. 이를 위하여, 도 7에 도시된 바와 같이 가상 주기 구조물(200b)을 직사각형 단면 형태의 복수의 층(1 내지 L)으로 나눈다.
이어서, 가상 주기 구조물의 분할된 층에서의 유전율 함수를 푸리에(Fourier) 급수로 전개한다. 가상 주기 구조물(200b)에 빛이 입사한다고 가정할 때, 입사파, 반사파, 투과파의 전자기파를 평면파의 합으로 두고, 각 층에서의 맥 스웰 방정식의 해를 고유함수 모드의 합으로 전개하여 그 전개계수들을 경계치 조건을 사용하여 결정한다.
상기 전개계수를 이용하여 TE 모드 전기장의 입사파에 대한 반사파 및 투과파의 진폭 또는 위상, 및 TM 모드 자기장의 입사파에 대한 반사파 및 입사파의 진폭 또는 위상을 산출하고, 영차 구조물의 반사율
및 투과율
을 계산한다.
이때, 가상 주기 구조물(200b)에서 분할된 층의 수와 푸리에 급수 전개의 항의 수를 크게 하면 보다 정밀한 해를 얻을 수 있다.
상기 가상 주기 구조물의 각 층에서의 푸리에(Fourier) 급수 전개계수는, 마루 영역을 구성하는 제1물질 및 제2물질과 골영역을 구성하는 제3물질 각각의 복소 굴절률, 주기 Λ에 대해서 분할된 층 l 의 상기 제1물질 및 제2물질이 차지하는 영역의 비율(f l ), 상기 제2물질만이 차지하는 비율(2δ l /Λ) 및 분할된 층의 첫 번째 층의 중심에 대해서 층 l 의 중심이 x축 방향으로 얼마나 떨어져 있는지를 나타내는 변수(t l )를 사용하여 나타낼 수 있으며, 이에 대해서는 후술한다.
섭동 영역에서의 영차 구조물의 그린함수도 상기 영차 구조물의 투과율 또는 반사율을 계산하는 방법을 적용하여 구할 수 있다.
리프만-슈빙거 적분방정식을 통하여 섭동 영역의 각 층에서 TE 모드 전기장과 TM 모드 자기장의 결합파 전개계수를 계산한 다음, 이들을 사용하여 섭동된 구조물에 의한 반사율
및 투과율
을 각각 계산한다.
본 발명에서는 입사 매질(영역 I)과 기판(영역 II) 사이의 주기 구조물을
층으로 나누는데 일부 또는 전체의 층이 균일한 물질로 이루어진 것을 포함한다. 섭동 영역 전체의 층 모두가 각 층 대로 균일한 물질로 이루어진 경우는 다층 Ga 구조, 모든 층이 하나의 균일한 물질로 이루어진 경우는 단층 Ga 구조, 섭동 영역을 기판 위의 공기층으로 두는 경우는 G0 구조로 정의하고, 각 층이 균일한 물질은 아니면서 모든 층이 주어진 층 하나의 반복으로 이루어진 경우는 단층 G1 구조로 정의할 수 있다.
본 발명에서는 후술하는 바와 같이 특히 표면층이 존재하는 주기 구조물을 고려하여 상기 G0, Ga, G1 구조에 대한 구체적인 계산 결과를 제공한다.
측정값과
산출값의
비교
주기 구조물의 실제 측정된 반사율 및 투과율 또는 그와 관계된 물리량과 가상 주기 구조물의 산출된 반사율 및 투과율 또는 그와 관계된 물리량을 비교하여, 일정 오차 범위 이내에서 동일할 경우, 실제 주기 구조물의 구조가 가상 주기 구조물(200b)의 구조와 일치한다고 판단할 수 있다.
측정된 반사율 및 투과율과 산출된 반사율 및 투과율을 비교함에 있어서, 컴퓨터 등의 별도 장치를 이용할 수 있으며 측정값과 산출값의 비교를 위한 소프트웨어를 이용할 수도 있을 것이다. 이와 같은 방법으로 주기 구조물(200)의 실제 구조를 정밀하게 검사할 수 있다.
주기 구조물의 측정 실험에서는 반사율(R)과 투과율(T)에 관련된 다양한 물 리량을 측정할 수 있으며, 측정된 물리량과 관련된 추정값을 수학적인 방법에 따라 계산된 반사율(R)과 투과율(T)을 사용하여 얻을 수 있다.
실시예
본 발명의 일실시예로서, 도 8에서와 같이 표면층이 형성된 가상 주기 구조물을 설정하였다. 이 주기 구조물은 60nm의 반복적인 주기를 가지며, 실리콘 결정으로 이루어진 마루 영역의 두께가 45nm이고, 마루 영역 표면에 산화물이 5nm 두께로 형성되어 있으며, 마루의 기울기는 81°인 경우로 설정하였다.
이러한 가상 주기 구조물의 기울기, 주기와 마루의 폭 및 높이는 임의의 값으로 얼마든지 변경 가능하다. 또한, 그 기하학적인 형태도 임의의 모양을 가질 수 있다.
이와 같은 가상 주기 구조물에 대하여 본 발명에 따른 방법을 이용하여 실제 주기 구조물의 측정값과 비교 가능한 물리량들을 산출하였다.
도 9 및 도 10은 TE 모드 및 TM 모드에서 도 8에서와 같이 산화막을 고려한 가상 주기 구조물과 산화막을 고려하지 않은 가상 주기 구조물의 반사율을 계산한 결과를 도시한 그래프들이다. 즉, 가상 주기 구조물이 산화막을 고려한, 즉 산화막이 도포된 실리콘으로 이루어진 경우와 산화막을 고려하지 않은, 즉 순수 실리콘으로 이루어진 경우를 가정하여 입사 파장에 대한 TE 모드와 TM 모드에서의 반사율을 계산한 것이다.
도시된 결과에 따르면, 순수 실리콘의 반사율과 산화막이 도포된 실리콘의 반사율의 계산값에 차이가 상당히 존재함을 확인할 수 있다. 즉, 산화막 등의 표면층의 존재를 가정했는지 유무에 따라 가상 주기 구조물이 실제 주기 구조물에 근접하는 정도가 크게 달라진다. 따라서, 산화막 등의 표면층이 형성되어 있는 실제 주기 구조물의 구조를 측정하는 경우, 가상 주기 구조물을 순수 실리콘으로만 이루어진 것으로 가정하여 계산하는 경우에는 절대로 실제 주기구조물과 일치하는 계산값을 구해낼 수가 없으며 산화막이 도포된 실리콘으로 가정하여 계산하는 경우에만 실제 주기 구조물과의 일치하는 정도를 증가시킬 수 있음이 자명하다.
도 11 및 도 12는 도 9 및 도 10의 계산 결과로부터 일반적인 타원편광분석기(ellipsometry)에서 측정하는 물리량들을 추출한 결과를 도시한 그래프이다. 이 그래프들은 본 발명에 따라 산출된 결과로부터 다양한 물리량들을 얻을 수 있음을 보여주고 있으며, 구체적으로는 TM모드(Rp)와 TE모드(Rs)의 반사율의 비의 크기와 위상을 일반적인 타원편광분석법에서 측정하는 물리량들 (ψ, Δ)로 계산한 결과를 도시하고 있다. 이 그래프들의 결과에서도, 순수 실리콘으로 이루어진 가상 주기 구조물과 산화막이 도포된 실리콘으로 이루어진 가상 주기 구조물의 분석 결과에 큰 차이가 있음을 확인할 수 있다. 이와 같이, 반사율 또는 투과율로부터 다양한 물리량을 추출하여 실제 주기 구조물의 구조와 비교할 경우, 산화막의 존재를 고려하여 가상 주기 구조물의 구조를 가정할 경우에 보다 정확한 구조를 예측할 수 있다.
이하에서는 도 6 및 도 7에 도시된 가상 주기 구조물(200b)의 반사율 및 투 과율을 계산하는 과정에 대해 상술한다.
I.
리프만
-
슈빙거
방정식:
TE
모드
주기 구조물의
축 방향 주기성에 의하여 맥스웰(Maxwell) 방정식의 TE 모드 해는,
식 (1):
식 (2):
로 나타낼 수 있다. 여기서
이고, 플로퀘(Floquet) 조건에 의해
,
이고,
는 진공에서의 입사파의 파장이며,
는 입사영역의 굴절율이며,
는 입사각이다. 또한,
는 진공의 유전율이며,
는 진공의 투자율이다.
축 방향으로 주기
를 갖는 주기가 1차원인 주기구조의 유전율 함수
를 푸리에(Fourier) 급수로 다음과 같이 전개할 수 있다.
식 (3):
맥스웰 방정식으로부터 식 (1)의
와 식 (2)의
는 다음과 같이 결합되어 있다.
식 (4):
식 (5):
이며, K
x는
원소가
인 대각행렬이며,
은
원소가
인 유전율 조화성분으로 만들어지는 행렬이다.
이하, 결합파 기저(coupled-wave basis)를 나타내는 첨자의 범위를
에서
까지로 제한할 것이다.
식 (4)에 있는 두 1차 미분방정식을 결합하면,
식 (6):
이 얻어진다.
영차 주기구조물에 대한 방정식을 첨자 [0]을 사용하여 나타내면 유사한 식,
식 (7):
을 얻으며, 이에 대응하는 그린 함수는,
식 (8):
로 나타낼 수 있다.
표준적인 방법에 따라 식 (6), 식 (7), 식 (8)을 결합하면 하기의 리프만-슈빙거(Lippamann-Schwinger) 방정식이 얻어진다.
식 (9):
여기서,
이다. 식 (9)에서 섭동 영역 밖에서는
이므로 적분영역은 섭동 영역으로 제한된다.
이 절에서는 영차 주기구조에서의 관련 양을 나타내는 첨자 [0]을 생략한다.
A. 주기구조 영역의 층 나누기
산화막, 코팅막 또는 표면층이 있는 주기구조를 공통의 주기
를 갖는 평행한 직사각형 모양의 층이 쌓인 것으로 근사한다. 이렇게 함으로써 유전율 함수
의
의존성은 층을 나타내는 지수
에게 전가된다. 그러면 주어진 층
에서의 유전율 함수는
로 나타나고 다음과 같이 표현될 수 있다.
식 (10):
식 (10)에서의 전개계수
은 각 층에서 다음과 같이 주어진다.
식 (11):
식 (12):
식 (13):
식 (14):
식 (15):
식 (16):
여기서,
는 층
에서 마루(ridge) 영역, 골(groove) 영역, 표면층 영역에서 각각의 복소수 굴절율을 나타내며,
은 주기
에 대해서 층
의 마루 영역과 표면층 영역이 차지하는 영역의 비율이며,
는 표면층 영역만이 차지하는 비율이며,
(단,
)은 제1층의 중심에 대해서 층
의 중심이
축 방향으로 얼마나 떨어져 있는지를 나타내는 변수이다.
상기 층 나누기에서 일부 또는 전체의 층이 균일한 물질로 이루어진 것을 포 함한다. 전체의 층 모두가 각 층 대로 균일한 물질로 이루어진 경우는 상술한 바와 같이, 다층 Ga 구조이며, 모든 층이 하나의 균일한 물질로 이루어진 경우는 단층 Ga 구조이다. 또한 각 층이 균일한 물질은 아니면서 모든 층이 주어진 층 하나의 반복으로 이루어진 경우는 단층 G1 구조이다. 특히, 도 6의 영역 II가 공기층일 수도 있다.
영역 Ⅰ과 Ⅱ에서 전기장과 자기장의 해는 플로퀘(Floquet) 조건에 의해 다음과 같이 나타낼 수 있다.
식 (17):
식 (18):
식 (19):
식 (20):
여기서,
식 (21):
식 (22):
식 (23):
식 (24):
식 (25):
이며,
과
는 입사 매질과 기판의 복소수 굴절율이다. 식 (21)의 우변의 첫째 항은 입사파를 나타낸다.
주기구조영역(
부터
까지)에서의 전기장과 자기장을 푸리에(Fourier) 급수 전개에 의해 다음과 같이 결합파 기저로 전개할 수 있다.
식 (26):
식 (27):
식 (7)을 층
에 적용하면 식 (26)의
(
부터
까지)들이 만족해야 할 방정식,
식 (28):
를 얻으며, 여기서,
식 (29):
식 (28)의 해
와
를 다음과 같이 적을 수 있다.
식 (30):
식 (31):
여기서,
은
의
개의 고유벡터들로 만들어지는 정방행렬이며,
은
의
개의 고유값의 양의 제곱근을 원소로 갖는 대각행렬이며,
이다.
과
은 경계조건의 적용에 의해 정해져야 하는 열벡터이다.
은 층
과 층
의 경계면의
좌표값이다.
식 (21), 식 (22), 식 (23), 식 (24)에 있는 결합파 기저의 성분들은 다음과 같이 행렬형태로 적을 수 있다.
식 (32):
식 (33):
식 (34):
식 (35):
여기서,
와
은
번째 대각원소가 각각
과
인 대각행렬이며,
는 모두
개의 성분을 갖는 열벡터이며,
,
,
,
이다.
(
부터
까지)에서 전기장과 자기장이 만족해야 할 경계조건으로부터,
식 (36):
을 얻는다. 여기서,
식 (37):
이고,
식 (38):
이다.
이 대칭행렬이므로
의 고유벡터들은 서로 수직이 된다. 따라서 규격화된 고유벡터 행렬
을 사용하면 수치계산의 효율을 높이기 위해
대신
를 사용할 수 있다.
식 (36)에서의 순환식은 궁극적으로 다음 형태의 식을 준다.
식 (39):
식 (40):
을 얻는다. 이
과
을 식 (36)에 대입하면
(
부터
까지)에 대한 순환식,
식 (41):
을 얻는다. 초기화 값
을 사용하여
에서 시작하여 이 순환식을 되풀이하여 적용하면 모든 층에서의
값들이 얻어진다.
이제 식 (40)의
(
부터
까지)을 얻을 수 있는 알고리듬이 구해졌으 므로 남은 문제는
들을 모두 구해야 하는 것이다.
을 알기 때문에 입사 매질로부터 층
로 전환행렬법(transfer matrix method)을 적용함으로써
을 구할 수 있다.
이를 위해 식 (36)의 순환식을 다음과 같이 변형시킨다.
식 (42):
여기서,
식 (43):
이다. 초기화값
을 가지고 식 (42)를 되풀이하여 적용하면 모든
값들이 정해진다. 이로써 앞에서 구한
과 식 (40)에 의해
이 결정되어 식 (30)의
가 구해진다. 이 절에서 첨자 [0]을 생략한 것을 생각하면, 식 (9)에 나오는 좌변의
가 결정된 것이다. 다음 절에서 식 (9)의 피적분 함수에 나오는
을 구한다.
식 (44):
로 된다. 이 절에서도 II절에서와 같이 영차 주기구조물과 관련된 양을 나타내는 첨자 [0]을 생략한다. 주어진 층
에 대해서
은 더 이상
의 함수가 아니다. 식 (44)의 해는 다음과 같이 나타낼 수 있다.
식 (45):
여기서,
과
은 III절에서와 같이
의 고유벡터들과 고유값들로부터 만들어지는 정방행렬이며, 주어진 층
에서
에 무관하고
에 의존하는 정방행렬
과
은
와
가 만족해야 하는 경계조건에 의해 결정된다.
,
,
의 첨자
에서 첫 번째
은 관찰점
가 층
내부에 위치함을 뜻하고 둘째
은 원천점
이 층
내부에 위치함을 뜻한다. 영역 I은 첨자
을, 영역 II는 첨자
을 부여하면 이들 영역에서의 해는,
식 (46):
로 나타낼 수 있다. 여기서,
,
,
,
은 앞에서 정의한 것과 같으며,
과
은 영행렬이다(앞 절에서
는 영이 아닌 열벡터이다.).
A. 기판(
)으로부터
원천층(
)으로
전환행렬법의
적용
과
가
(
부터
까지)에서 만족해야 하는 경계조건으로부터 식 (36)과 비슷한 형태의 순환식,
식 (47):
을 얻는다. 여기서
과
는 식 (37)과 (38)에서 정의된 대로이다. 식 (47)을 되풀이하여 적용하면 궁극적으로,
식 (48):
식 (49):
가 얻어진다. 식 (49)를 식 (47)에 대입하면 식 (41)과 완전히 같은 순환식,
식 (50):
이 얻어진다. 초기화 값
도 같으므로 같은 계산을 되풀이하지 않고 이미 II절에서 구한
의 값을 그대로 사용할 수 있다.
과
을 연결짓는 식 (49)와
을 구하는 알고리듬을 주는 식 (50)이 얻었으므로
인
을 구하는 문제에서 남은 일은
을 구하는 것이다. 첫 단계로
을
으로 표현한다. 이를 위해 식 (47)을 다음과 같이 변형시킨다.
식 (51):
이 식은 식 (49)와 함께 다음의 순환식을 준다.
식 (52):
이 식을 반복하여 적용하면
부터
까지의 모든
에 대하여,
식 (53):
을 얻는다. 따라서
이 일단 구해지면
과
이 모두 구해진다.
를 구하기 전에
인
을 구하는 문제로 돌아간다. 이는
과
을 구하는 문제이다.
B.
입사매질(
)로부터
원천층(
)으로
전환행렬법의
적용
과
가
(
부터
까지)에서 만족해야 하는 경계조건으로부터,
식 (54):
을 얻는다. 여기서
과 N
l ±는 앞에서 정의한 대로이다. 식 (54)를 반복하여 적용하면,
식 (55):
을 얻고 이 식에서
을 소거하면 다음과 같이
을
으로 표현할 수 있다.
식 (56):
식 (56)을 사용하여 식 (54)에서
과
을 소거하면 순환식,
식 (57):
을 얻는다. 초기화값
을 가지고
에서 출발하여 이 순환식을 되풀이하여 적용하면 모든
이 구해진다.
를 알고리듬에 의해 구해지는
과 아직 미지의 값
으로 표현했으므로 남은 문제는
을 구하는 것이다. 이를 위해
을
으로 표현한다. 식 (47)에 주어진 것과 같은 형태의 식,
식 (58):
이
일때
에서 경계조건의 적용에 의해 얻어진다. 이 식과 식 (56)으로부터
에 대한 순환식,
식 (59):
을 얻고 이것으로부터
(
부터
까지)를
으로 표현한 식,
식 (60):
을 얻는다.
이 구해지면
과
(
부터
까지)들이 구해진다.
이제
과
을 구하자.
에서
이 연속이고 식 (44)의 우변의 델타함수에 의해
가 불연속이므로,
식 (61):
를 얻는다.
을 사용하여 식 (61)을 다음과 같이 표현할 수 있다.
식 (62):
식 (63):
식 (62), 식 (63), 식 (53), 식 (60)에 의해 식 (45)의 그린 함수가 완전히 정해진다.
IV
. 수치해석:
TE
모드
상술한 바와 같이, G0는 섭동 영역이 단층 균일매질이며, Ga는 주어진 층이 균일매질인데 단층과 다층으로 나눌 수 있다. G1은 주어진 층이 균일매질이 아니며 역시 단층과 다층으로 나눌 수 있다. G1 다층은 가장 일반적인 경우로서 위에서 언 급한 모든 경우를 포함한다. 따라서, 이하 G1 다층에 대해 분석한다.
A.
리프만
-
슈빙거
적분방정식의 이산화
영차(zeroth order) 주기 구조물의 물리적 성질 및 기하학적 구조와 입사파의 정보에 의해 정해진 두 양
및
과 섭동 퍼텐셜
를 입력변수로 하여, 섭동된 주기구조물의 물리적 성질, 기하학적 구조, 그리고 동일한 입사파 정보에 의해 결정되어야 할 미지수
를 리프만-슈빙거 방정식의 풀이법인 반복법(iteration)을 사용하여 구할 것이다.
층 (
부터
까지) 내부에 위치하는
를
로 보내고
층 내부에 위치하는
를
로 보내고 식 (53)과 (60)을 사용하여 층
에서 섭동 함수
이 상수값
를 갖는다고 가정하면 식 (9)가 다음과 같이 됨을 볼 수 있다.
식 (64):
식 (65):
식 (64)에서
일 경우 우변의
은 어떤 항도 발생시키지 않으며
일 우변의
은 어떤 항도 발생시키지 않음에 주의해야 한다.
정밀도를 높이기 위해서 내삽 공식,
식 (66):
식 (67):
을 사용하여 식 (64)와 식 (65)에 나오는 적분을 해석적으로 수행하면,
부터
까지에 대해 식 (64)는,
식 (68):
식 (69):
식 (70):
로 표현된다.
또한, 식 (65)는,
식 (71):
로 표현된다.
식 (68), 식 (69), 식 (70), 식 (71)을 모아서 다음과 같이 하나의 확장된 행렬형태의 1차 방정식계로 나타낼 수 있다.
식 (72):
여기서
와
는
개의 층성분
과
을 각각 갖는 열벡 터이며, 각 층성분
과
은 모두
개의 결합파 기저성분을 가진 열벡터이다.
는
개의 층성분을 갖는 정방행렬이며 각각의 층성분은 아래와 같이 주어지며 모두
개의 결합파 기저성분을 갖는 정방행렬이다.
식 (73):
식 (74):
식 (75):
식 (76):
B. 특수한 사례
특수한 경우로,
(i) 각 층이 균일 매질이면 (다층 Ga 구조),
식 (77):
로 되고, 이로 인해
이 대각행렬이 되고, 그 결과 식 (73) 및 식 (74)의 행렬요소
과
이 모두 대각행렬이 되어 더욱 간단해진 표현식을 얻을 수 있으며,
(ii) 각 층의 영차 유전율 함수가 모두 동일할 경우(단층 G1 경우), 즉
일 때,
식 (78):
식 (79):
에 의해 간단해진 표현의 행렬요소
과
을 얻을 수 있으며,
(iii) 더 나아가 각 층이 균일하면(단층 Ga 구조),
식 (80):
식 (81):
에 의해 더욱 간단한 표현의 행렬요소를 얻을 수 있으며,
(iv) 섭동영역을 기판 바로 위의 입사매질로 취할 경우(G0 구조), 단층 Ga의 구조에 더하여
으로 둘 수 있고 행렬요소도 훨씬 간단해진다.
결합파 기저 공간에서의 행렬과 벡터의 곱의 계산속도를 높이기 위해 다음의 행렬,
식 (82):
(여기서
는 K×K 단위행렬이며
는
원소가
로 주어지는 K×K 정방행렬이며
와
은 크기가
인 영 열벡터와 영 행벡터임)에 의한 상사변환에 의해 섭동행렬을 블록대각화 시킬 수 있다.
식 (72)의 행렬 방정식을 반복법으로 풀기 위해 초기화 값으로
=
으로 시작한다. 이것을 식 (72)에 대입하면,
식 (83):
을 얻는다. 주어진 오차한계 내에서 식 (83)이 성립하지 않으면 새로운 값
를 식 (83)의 우변의 값으로 둔다.
식 (84):
이 값을 다시 식 (72)에 대입하여,
식 (85):
을 얻는다. 이 식의 좌우변의 값을 비교하여 만족스러운 값을 얻지 못하면 새로운
를 식 (85)의 우변의 값으로 둔다.
식 (86):
이
를 가지고 식 (72)가 성립하는지 시험한다. 만족할 만한
를 얻을 때까지 이 과정을 되풀이한다.
최종적인 계산목표는 결합파 기저의 차수
과 관련된 반사율
들과 투과율
들이다. 반사율 열벡터
는 다음과 같이 추출해낼 수 있다.
식 (87):
여기서
는 영차 반사율 행렬이다. 주차수(principal order) 반사 회절파의 반사율을 얻고자 하면
을 추출하면 된다.
식 (88):
으로 주어진다. 여기서
는 영차 투과율 행렬이며,
는 주차수(principal order) 투과 회절파의 투과율이다.
V.
리프만
-
슈빙거
방정식:
TM
모드
주기 구조물의
축 방향 주기성에 의하여 맥스웰(Maxwell) 방정식의 TE 모드 해
와
는,
식 (89):
식 (90):
로 나타낼 수 있다. 유전율 함수
의 역수를 푸리에(Fourier) 급수로 다음과 같이 전개할 수 있다.
식 (91):
*맥스웰 방정식에 의해
와
는 결합된 1차 미분방정식 2개를 만족하는데
를 소거하면
에 관한 1개의 2차 미분방정식,
식 (92):
이 얻어진다. 여기서,
식 (93):
이며
는
원소가
인 정방행렬이다. 동일 주기의 영차 구조물에 대해서 첨자 [0]을 사용하면 비슷한 식을 얻는다.
식 (94):
대응하는 그린 함수는,
식 (95):
로 주어진다. 표준적인 방법에 따라 식 (92), 식 (94), 식 (95)을 결합하면 다음의 TM 모드에 대한 리프만-슈빙거 방정식,
식 (96):
이 절에서는 영차 주기구조에서의 관련 양을 나타내는 첨자 [0]을 생략한다. 영역 Ⅰ과 Ⅱ에서의 전기장과 자기장의 해는 플로퀘(Floquet) 조건에 의해 결합파 기저로 전개한 형태로 나타낼 수 있다.
식 (97):
식 (98):
식 (99):
식 (100):
여기서,
식 (101):
식 (102):
식 (103):
식 (104):
이다. 식 (101)의 우변의 첫째 항은 입사파를 나타내는 항이다. 비슷하게, 주기구조영역(
부터
까지)에서의 전기장과 자기장을 푸리에(Fourier) 급수 전개에 의해 다음과 같이 결합파 기저로 전개할 수 있다.
식 (105):
식 (106):
식 (92)를 층
에 적용하면 식 (105)의
(
부터
까지)들이 만족해야 할 방정식,
식 (107):
를 얻으며, 여기서,
은 II절에서 정의한 대로이며
에 무관한
은 다음 과 같이 정의된다.
식 (108):
식 (107)의 해는 다음과 같이 주어진다.
식 (109):
식 (110):
여기서
은
의
개의 고유벡터들로 만들어지는 정방행렬이며,
은
의
개의 고유값의 양의 제곱근을 대각원소로 갖는 대각행렬이며,
이며,
과
은 경계조건의 적용에 의해 정해져야 하는
개의 성분을 갖는 열벡터이다.
식 (101), 식 (102), 식 (103), 식 (104)에 있는 결합파 기저의 성분들은 다음과 같이 행렬형태로 적을 수 있다.
식 (111):
식 (112):
식 (113):
식 (114):
여기서,
,
,
,
,
,
,
,
은 II절에서 정의한 대로이며,
,
이다.
(
부터
까지)에서 전기장과 자기장이 만족해야 할 경계조건을 적용하면 식 (36)에서 주어진 순환식이 나온다. 계산의 나머지 부분은 TE 모드에서와 정확히 같은 방법으로 행해질 수 있다.
에 나타나는
의 역행렬과 관련한 논평 하나로 이 절을 마무리한다.
비록
과
이 대칭행렬이더라도 그 곱
은 일반적으로 대칭행렬이 아니다. 따라서
의 고유벡터들
(
부터
까지)은 서로 수직이 아니다. 즉, m≠n에 대해
이다. 그러므로 이 경우 적절한 규격화에 의해
은 단위행렬
로 표현될 수 없다. 그 대신
은 대각행렬이 될 수 있기 때문에
로 규격화시킬 수 있다. 이 경우
대신
를 사용할 수 있다.
식 (115):
로 된다. 이 절에서도 II절과 VI절에서와 같이 영차 주기구조물과 관련된 양을 나타내는 첨자 [0]을 생략한다. 주어진 층
에 대해서
과
은 더 이상
의 함수가 아니다.
과
이
의 고유벡터와 고유치에 의해 만들어진 행렬이라는 것을 제외하고는 식 (115)의 해가 식 (45)의 TE 모드의 해와 완전히 같은 형태이다. 이때 같은 표현의
와
를 가진다. 이 절에서의
과
은 물론 VI절에서 구한 것과 완전히 같다. 영역 I로 표시되는 입사매질과 영역 II로 표시되는 기판에서의 해는 식 (46)에 정의된 그대로이다.
VIII
. 수치해석:
TM
모드
수치계산의 정밀도를 높이기 위하여 식 (66)과 비슷한
에 대한 내삽 공식을 사용하여 적분을 해석적으로 수행하면 식 (97)의 리프만-슈빙거 적분방정식은 다음과 같이 이산화된다.
식 (116):
식 (117):
식 (118):
식 (119):
식 (116), 식 (117), 식 (118), 식 (119)를 모아서 다음과 같이 하나의 확장된 행렬형태의 1차 방정식계로 나타낼 수 있다.
식 (120):
여기서
는
개의 층성분
과
을 각각 갖는 열벡터이며 각 층성분
과
이 모두
개의 결합파 기저 성분을 가진 열벡터이다.
는 식 (73), 식 (74), 식 (75), 식 (76)에서 정의되었으며,
은
개의 층성분을 갖는 정방행렬이며 각각의 층성분은 다음과 같이 주어지며 모두
개의 결합파 기저 성분을 갖는 정방행렬이다.
식 (121):
식 (122):
식 (123):
식 (124):
TE 모드에서와 비슷하게 다층 Ga 구조, 단층 G1 구조, 단층 Ga 구조, G0 구조에 대한 대응되는 식을 얻을 수 있다.
식 (120)의 행렬 방정식은 초기화 값
을 가지고 TE 모드에서와 같이 반복법으로 푼다.
반사율 열벡터
는 다음과 같이 추출해낼 수 있다.
식 (125):
여기서
는 영차 반사율 행렬이다. 주차수(principal order) 반사 회절파의 반사율을 얻고자 하면
을 추출하면 된다.
식 (126):
으로 주어진다. 여기서
는 영차 투과율 행렬이며,
는 주차수(principal order) 투과 회절파의 투과율이다.
이렇게 하여 계산한 TE 모드와 TM 모드에 대한 각각의 반사율 및 투과율을 실제 측정한 반사율 및 투과율과 비교함으로써, 여러 주기 구조물, 예를 들어 홀로그래픽 격자, 표면부조 및 다층 격자구조, 평면 유전 또는 흡수 홀로그래픽 격자, 임의 단면 유전체 및 흡수 표면부조 격자, 2차원 표면부조 격자 또는 비등방성 격자구조의 비파괴적 분석에 적용될 수 있다. 한편, 상기 주기 구조물은 상술한 예에 한정되지 않음은 물론이다.
이상에서 설명한 본 발명의 실시예는 단지 예시의 목적을 위해 개시된 것이고, 본 발명에 대한 통상의 지식을 가지는 당업자라면 본 발명의 사상과 범위 안에서 다양한 수정, 변경, 부가가 가능할 것이다.