광격자

1. 개요[편집]

광격자 Optical Lattice
퍼텐셜	$V(x)=V_0\sin^2(kx)$, 격자 간격 $d=\lambda/2$
에너지 단위	반동에너지 $E_R=\hbar^2k^2/2m$, 깊이 $s=V_0/E_R$
고유값 문제	마티외 방정식 ($q=s/4$)
깊은 격자	현장 진동수 $\hbar\omega=2\sqrt{V_0E_R}$, 터널링 $J\propto s^{3/4}e^{-2\sqrt s}$
유효 모형	보스-허바드 / 페르미-허바드
대표 실험	초유체-모트 절연체 전이 (Greiner 외, 2002)
응용	양자 시뮬레이터 · 광격자 시계 · 양자기체 현미경

결정을 이해하려고 결정을 계산하는 대신, 결정을 하나 만들어서 직접 물어보는 쪽으로 가면 어떨까. 단 이온은 빼고, 불순물도 빼고, 격자 상수는 광학 파장 규모로 키워서.

광격자(optical lattice)는 마주 오는 레이저 빔이 만드는 정상파의 주기적 세기 분포를 쌍극자 퍼텐셜로 삼아, 중성원자를 인공 결정 격자에 가두는 계다. 1차원에서 퍼텐셜은

$$V(x) = V_0\sin^2(kx), \qquad k=\frac{2\pi}{\lambda}$$

이고 격자 간격은 빔 파장의 절반 $d=\lambda/2$ 다. 실제 고체의 격자 상수보다 세 자릿수쯤 크고, 이온도 결함도 포논도 없으며, 무엇보다 깊이 $V_0$ 와 원자 간 상호작용을 실시간으로 돌릴 수 있다. 응집물질 이론이 수십 년간 종이 위에서만 다루던 모형 해밀토니안이, 여기서는 실험 장치의 노브가 된다.

원리는 광집게와 같은 쌍극자 힘이다. 레이저가 원자에 유도하는 쌍극자 모멘트가 다시 그 빔의 세기 기울기를 느끼며, 퍼텐셜은 대략 세기를 이탈맞춤(detuning)으로 나눈 값에 비례한다. 적색 이탈맞춤이면 원자가 세기 최대점(정상파의 배)으로, 청색이면 최소점(마디)으로 모인다. 원자를 여기까지 데려오는 것은 레이저 냉각과 증발 냉각의 몫이고, 이 문서는 그 다음부터를 다룬다.

2. 반동에너지 — 이 바닥의 자연 단위[편집]

광격자 문헌을 열면 온도도 깊이도 터널링도 전부 반동에너지

$$E_R = \frac{\hbar^2k^2}{2m}$$

단위로 적혀 있다. 광자 하나를 흡수한 원자가 얻는 운동에너지이며, 격자 문제의 유일한 자연스러운 에너지 척도다. 무차원 깊이 $s\equiv V_0/E_R$ 하나로 밴드 구조가 통째로 결정되기 때문에, 원자 종류와 파장을 바꿔도 $s$ 만 같으면 밴드 구조는 같다. 전형적인 실험 값은 $s\sim 5\!-\!30$ 이고, $s\gtrsim 15$ 면 대개 “깊은 격자”로 취급한다.

3. 마티외 방정식과 블로흐 밴드[편집]

주기 퍼텐셜이므로 블로흐 정리가 그대로 적용된다. 정상상태 슈뢰딩거 방정식

$$-\frac{\hbar^2}{2m}\psi'' + V_0\sin^2(kx)\,\psi = E\psi$$

에 $z=kx$ 를 넣고 $\sin^2 z=(1-\cos 2z)/2$ 로 풀어 쓰면

$$\frac{d^2\psi}{dz^2} + \left[\frac{E}{E_R}-\frac{s}{2}+\frac{s}{2}\cos 2z\right]\psi = 0$$

이 되어, 정확히 마티외 방정식 $\psi''+(a-2q\cos 2z)\psi=0$ 의 꼴이다. 대응은

$$a = \frac{E}{E_R}-\frac{s}{2}, \qquad |q| = \frac{s}{4}$$

이며 $q$ 의 부호는 원점을 $\pi/2$ 옮기면 뒤집히는 관례 문제다. 마티외 방정식의 특성값 곡선이 곧 광격자의 밴드 가장자리이고, 인스-스트럿 선도의 불안정 혀가 밴드갭에 대응한다. 응집물질에서 배운 밴드 구조와 수학물리에서 배운 마티외 함수가 여기서 같은 그림의 앞뒤라는 사실이 이 계의 교육적 매력이다.

수치적으로는 마티외 함수를 부를 필요도 없다. 평면파 기저 $\psi_q(x)=\sum_{n=-n_{\max}}^{n_{\max}} c_n\,e^{i(q+2nk)x}$ 로 전개하면 해밀토니안이 대각선에 $(q+2nk)^2\hbar^2/2m + V_0/2$, 부대각선에 $-V_0/4$ 인 삼중대각 실대칭 행렬이 되고, $n_{\max}\sim10$ 이면 낮은 밴드는 기계 정밀도로 수렴한다. 물리 계산 숙제로 15분이면 끝나는 문제인데 실험 전체의 해석 기반이 된다.

4. 깊은 격자 — 조화 근사, 터널링, 현장 상호작용[편집]

$s\gg1$ 이면 각 우물 바닥을 조화 퍼텐셜로 근사할 수 있다. $V_0\sin^2(kx)\approx V_0k^2x^2$ 에서

$$\hbar\omega = 2\sqrt{V_0 E_R} = 2\sqrt{s}\,E_R$$

를 얻는다. 밴드갭이 $\sqrt s$ 로 커지고, 밴드폭은 지수적으로 줄어든다. 후자가 결정적이다. 최저 밴드의 폭 $4J$ 에서 정의되는 이웃 간 터널링 진폭은 깊은 격자에서

$$\frac{J}{E_R} \simeq \frac{4}{\sqrt\pi}\,s^{3/4}\,e^{-2\sqrt{s}}$$

라는 점근형을 갖는다. $s$ 를 10에서 20으로 올리는 것만으로 $J$ 가 한 자릿수 넘게 떨어지므로, 레이저 세기 하나로 「금속 ↔ 절연체」 사이를 왕복할 수 있다. 실제 고체에서는 상상도 못 할 조작이다.

한편 원자끼리는 s-파 산란으로 접촉 상호작용을 하며, 같은 우물에 두 원자가 있을 때 드는 에너지가 현장 상호작용

$$U = \frac{4\pi\hbar^2 a_s}{m}\int |w(\mathbf r)|^4\,d^3r$$

다. $a_s$ 는 s-파 산란길이, $w$ 는 최저 밴드의 와니에 함수다. 깊은 격자에서 $w$ 를 가우스로 근사하면 $U$ 가 $s^{3/4}$ 로 완만히 커지는 반면 $J$ 는 지수적으로 죽으므로, $U/J$ 는 격자 깊이만으로 몇 자릿수를 훑는다. $a_s$ 자체를 페슈바흐 공명으로 따로 조절할 수도 있으니 노브가 두 개다.

여기서 와니에 함수가 왜 중요한지가 분명해진다. 블로흐 함수는 퍼져 있어 「어느 우물에 있는가」를 말할 수 없고, 국소화된 와니에 기저로 갈아타야 $J$ 와 $U$ 라는 격자 모형의 계수가 정의된다. 위상 선택을 제대로 해야 최대 국소화된 와니에 함수를 얻는다는 사실은 밴드 구조 계산 쪽 이야기와 완전히 같다.

5. 보스-허바드 모형과 초유체-모트 절연체 전이[편집]

최저 밴드만 남기면 유효 해밀토니안이 보스-허바드 모형이 된다.

$$\hat H = -J\sum_{\langle i,j\rangle}\hat b_i^\dagger \hat b_j + \frac{U}{2}\sum_i \hat n_i(\hat n_i-1) - \mu\sum_i \hat n_i$$

이 모형을 광격자에서 실현할 수 있다는 제안이 1998년에 나왔고(약슈·촐러 등), 2002년 그라이너 등이 $^{87}$Rb 응축체를 3차원 격자에 얹어 초유체-모트 절연체 전이를 관측했다. 실험의 신호가 우아하다. 격자를 끄고 자유 팽창시킨 뒤 흡수 영상을 찍으면 위상 결맞음이 있는 초유체 쪽에서는 브래그 조건을 만족하는 뚜렷한 간섭 봉우리가 나타나고, 격자를 깊게 해 모트 절연체로 넘어가면 봉우리가 사라져 뭉개진 분포가 된다. 다시 얕게 하면 봉우리가 되살아난다 — 가역적이므로 결어긋남으로 인한 파괴가 아니라 진짜 상전이라는 증거다.

두 극한의 대비는 명확하다.

	초유체 ($U/J$ 작음)	모트 절연체 ($U/J$ 큼)
바닥상태	모든 자리에 걸친 결맞은 상태	자리마다 정수 개 고정
자리당 원자수 요동	큼 (푸아송 꼴)	억제됨
위상 결맞음	있음	없음
들뜸	무갭 (음파 모드)	유갭 ($\sim U$)
압축률	유한	0

$n=1$ 채움에서 전이점의 평균장 추정은 $(U/J)_c = z\left(2n+1+2\sqrt{n(n+1)}\right)$ 이며 배위수 $z=6$ 인 3차원 정육면체 격자에서 약 35가 나온다. 그런데 양자 몬테카를로 계산은 약 29를 준다. 평균장이 요동을 과소평가해 전이를 늦게 보는 전형적인 양상이고, 1차원에서는 아예 전이의 성격 자체가 달라져(BKT 부류) 평균장이 정성적으로도 틀린다. 광격자가 「이론 검증대」로 쓰이는 이유가 여기 있다.

6. 무엇을 계산하는 도구인가 — 양자 시뮬레이터[편집]

파인만이 1982년에 던진 제안이 실물이 된 것이 이 계다. 격자 위 강상관 다체 문제는 힐베르트 공간이 지수적으로 커서 고전 컴퓨터로 정직하게 풀 수 없다.

• 정확 대각화: 자리 수 수십 개가 한계.
• DMRG·텐서 네트워크: 1차원에서는 사실상 정확하지만 2차원 이상에서 얽힘 면적 법칙이 불리해진다.
• 양자 몬테카를로: 보손계는 웜 알고리즘 등으로 매우 강력하지만, 페르미온과 좌절 스핀계는 부호 문제로 막힌다.

마지막 항목이 양자 시뮬레이터로서의 광격자의 존재 이유다. 페르미-허바드 모형의 도핑 영역은 고온초전도의 최소 모형으로 지목되면서도 고전 계산이 부호 문제로 접근 못 하는 대표 사례인데, $^{40}$K·$^6$Li 페르미 원자를 격자에 넣으면 자연이 대신 그 계를 시간 발전시켜 준다. 반강자성 스핀 상관 관측까지는 도달했고, 진짜 관심사인 저온 d-파 짝짓기 영역은 온도(엔트로피)를 더 낮추는 것이 남은 과제다. 냉원자 쪽에서 「온도」보다 「엔트로피 per 입자」를 말하는 이유가 이것이다.

주의할 점도 분명하다. 광격자가 구현하는 것은 이상화된 모형 해밀토니안이지 실제 물질이 아니다. 장거리 쿨롱, 실제 포논, 다중 궤도, 무질서가 전부 빠져 있다. 그것이 장점(깨끗한 모형 검증)이자 한계(구리 산화물의 실제 $T_c$ 는 안 나온다)다. 그리고 검증(V&V)의 관점에서 골치 아픈 문제가 하나 더 있다 — 양자 시뮬레이터의 답을 무엇으로 검산할 것인가. 고전 계산이 안 되는 영역을 풀려고 만든 장치인데, 바로 그 영역에서 정답을 모른다. 현재의 답은 고전 계산이 가능한 극한에서 교차 검증하고, 서로 다른 실험군의 장치를 비교하는 것이다.

7. 광격자 시계와 양자기체 현미경[편집]

광격자 시계는 완전히 다른 방향의 응용이다. 원자를 격자에 가두면 램-디케 영역에 들어가 반동 이동과 1차 도플러 이동이 사라지고, 중성원자를 수천 개 동시에 분광할 수 있어 신호가 커진다. 문제는 격자 레이저 자체가 두 시계 준위를 서로 다르게 밀어낸다는 것인데, **두 준위의 동적 편극률이 같아지는 「마법 파장」**을 고르면 그 차이가 1차에서 상쇄된다. $^{87}$Sr 의 경우 약 813 nm 가 그 파장이다. 이 아이디어로 광격자 시계는 세슘 마이크로파 표준을 여러 자릿수 앞질렀고, 지금은 상대론적 중력 적색편이를 실험실 책상 높이 차이에서 재는 수준까지 왔다.

양자기체 현미경은 반대로 공간 분해능 쪽 돌파구였다. 고해상도 대물렌즈로 격자 자리 하나하나를 분해해 형광 영상을 찍으면 원자 배치를 자리 단위로 직접 읽는다. 2009~2010년경 보손계에서 처음 구현됐고 이후 페르미온으로 확장됐다. 이전까지 냉원자 관측량은 자유 팽창 후의 운동량 분포(즉 파수공간 정보)가 거의 전부였는데, 이 장치가 실공간 상관함수를 직접 측정 가능하게 만들었다. 스핀-스핀 상관, 개수 요동, 심지어 얽힘 엔트로피 관련 양까지 측정 대상이 됐다. 다만 한 번 찍으면 원자가 데워져 날아가므로 매 실험이 파괴 측정이고, 통계는 반복 실행으로 쌓는다.

8. 수치적으로 다룰 때[편집]

• 밴드 구조: 위의 평면파 대각화. $q$ 를 브릴루앙 영역에서 훑으면 끝. 밴드폭에서 $J$, 곡률에서 유효질량을 읽는다.
• 와니에 함수와 $U$: 블로흐 함수의 위상을 실수화하도록 고정한 뒤 푸리에 합. 가우스 근사는 $s\gtrsim15$ 에서만 쓸 만하고, 얕은 격자에서는 꼬리가 두꺼워 $U$ 를 과대평가한다.
• 평균장 상도: 자리별 결맞은 상태 시늉으로 자기일관 풀이. 상도의 얼개는 잡지만 임계점은 위에서 본 대로 어긋난다.
• 1차원 동역학: DMRG/TEBD 로 격자 깊이를 시간에 따라 바꾸는 램프까지 그대로 모사 가능. 「얼마나 천천히 램프해야 단열인가」가 실험 설계에 직결되는 질문이다.
• 평균 규모 기술: 격자가 아주 얕고 상호작용이 약하면 보스-아인슈타인 응축 문서의 그로스-피타옙스키 방정식에 주기 퍼텐셜을 얹는 것으로 충분하다. 깊은 격자에서 이걸 쓰면 안 된다 — 모트 상은 평균장 파동함수로 기술되지 않는다.

9. 여담[편집]

• 광격자를 “인공 결정”이라 부르지만 척도가 완전히 다르다. 격자 상수는 고체의 수백 배, 온도는 나노켈빈, 에너지 척도는 $E_R\sim h\times$ 수 kHz. 모든 것이 느리고 크고 차갑다. 그래서 고체에서는 펨토초 분광이 필요한 동역학을 여기서는 밀리초 단위로 카메라로 찍는다.
• 광결정과 이름이 헷갈리는데 정반대다. 광결정은 유전율을 주기적으로 배열해 빛의 밴드갭을 만들고, 광격자는 빛으로 주기 퍼텐셜을 만들어 물질의 밴드갭을 만든다. 지배 방정식이 같은 꼴이라 계산 코드는 서로 베껴 쓸 수 있다.
• 이 계에서는 「불순물이 없다」가 자랑인 동시에 골칫거리다. 진짜 물질의 흥미로운 현상 중 상당수가 무질서에서 오기 때문에, 냉원자 쪽에서는 스페클 패턴을 일부러 얹어 무질서를 의도적으로 주입한다. 앤더슨 국소화 관측이 그렇게 나왔다.¹
• 초기에는 격자를 “얼마나 깊게 했느냐”를 재는 것 자체가 기술이었다. 요즘 표준은 격자를 짧게 켰다 끄는 카피차-디랙 회절이나 밴드 간 진동수를 측정해 $s$ 를 역산하는 방식이다. 깊이를 모르면 $J$ 도 $U$ 도 모르고, 그러면 논문의 가로축이 통째로 틀린다.²
• 냉원자 논문에서 “온도”보다 “엔트로피”를 강조하는 습관이 처음 보면 이상한데, 격자를 켜고 끄는 과정이 단열이라면 보존되는 것이 엔트로피이지 온도가 아니기 때문이다. 그리고 이 바닥에서 진짜 어려운 일은 냉각이 아니라 엔트로피를 계 밖으로 빼내는 것이다.³

10. 관련 문서[편집]

• 보스-아인슈타인 응축 · 레이저 냉각 · 광집게 · 광결정
• 블로흐 정리 · 와니에 함수 · 밴드 구조 계산 · 마티외 방정식
• 상전이 · 상관함수 · 상태밀도 · 쿠퍼 쌍
• 양자 몬테카를로 · DMRG · 텐서 네트워크 · 힐베르트 공간
• 파울 트랩 · 페닝 트랩 · 언쇼 정리 · 슈뢰딩거 방정식
• 보스-허바드 모형 · 모트 절연체 · 광격자 시계 · 양자 시뮬레이터

11. Footnotes[편집]

1. 이 지점이 은근히 철학적이다. “깨끗하게 만들려고 20년 갈아 넣었는데 이제 일부러 더럽힌다.” 물론 차이는 크다 — 스페클 무질서는 통계 성질을 실험자가 정확히 알고 통제한다. 진짜 물질의 불순물은 아무도 모른다. ↩

2. 그래서 $s$ 교정은 이 바닥의 통과의례다. 두세 가지 독립적인 방법으로 재서 서로 몇 % 안에 들어오는지 보고, 안 맞으면 빔 정렬부터 다시 한다. 검증 및 확인의 정신이 실험실에서 나타나는 형태다. ↩

3. 증발 냉각은 뜨거운 원자를 버려서 엔트로피를 들고 나가게 하는 것이라 원자 수를 대가로 치른다. “원자 100만 개로 시작해서 10만 개 남았다”는 문장은 실패 보고가 아니라 정상적인 실험 일지다. ↩