자기장 하 블록 전자의 평균 밀도와 2차 비선형 응답

초록

이 논문은 3차원 다밴드 결정에서 약한 균일 자기장이 전자 평균 밀도에 미치는 1차와 2차 효과를 이론적으로 계산한다. 절연체에서는 선형 응답이 잘 알려진 Streda 공식과 일치하지만, 금속에서는 페르미면의 궤도 자기모멘트가 추가된다. 2차 응답에서는 퀀텀 메트릭 텐서가 기하학적(의사)자기모멘트를 생성함을 보이며, 이 결과가 부피·압력 변화와 연결된다. 전 과정은 게이지 불변이며, 밴드 내·외 전이와 완화 과정을 모두 포함한다.

상세 분석

본 연구는 Bloch 전자들의 평균 전자밀도 n을 자기장 B에 대한 멱급수 n=n^{(0)}+σ^{(1)}{a}B{a}+σ^{(2)}{ab}B{a}B_{b}+… 로 전개하고, 각각의 계수를 미세하게 유도한다. 선형계수 σ^{(1)}{a}에 대해서는, 화학퍼텐셜이 밴드갭 안에 있을 때는 Streda 공식 ∂n/∂B=−(e)^{-1}σ{H}와 정확히 일치한다는 것을 증명한다. 그러나 화학퍼텐셜이 밴드 내부에 위치한 금속에서는, 페르미함수의 미분 f’{n}가 비제로가 되므로, 페르미면에 존재하는 궤도 자기모멘트 m{n}이 추가 항 ∝∑{n,k}f’{n}m_{n}을 만든다(식 3). 이는 기존의 Streda 식에 대한 중요한 수정이며, 금속에서의 전자밀도 변화를 정확히 기술한다.

2차 응답 σ^{(2)}{ab}는 훨씬 복잡한 구조를 가진다. 저자들은 전자-전자 상호작용을 무시한 비상호작용 모델을 가정하고, Bloch 함수 |u{n\mathbf{k}}⟩의 k-공간 미분을 이용해 전자-자기장 상호작용을 정규화한다(q→0 한계에서 위치 연산자의 특이성을 회피). 이 과정에서 나타나는 핵심 물리량이 바로 퀀텀 메트릭 텐서 g_{ab}^{n}=Re⟨∂{a}u{n}|(1-P_{n})|∂{b}u{n}⟩이다. 이 텐서는 Bloch 상태 사이의 거리(복소 사영 평면에서의 회전)를 측정하며, 식 (4)와 (33)에서 보듯이 B^{2}에 비례하는 의사 자기모멘트 m^{(2)}{ab}∝e^{2}g{ab}v_{a}v_{b}를 생성한다. 즉, 실제 스핀·궤도 각운동량이 아닌, 파동함수의 기하학적 회전에 의해 발생하는 “기하학적 자기모멘트”가 2차 밀도 변화를 주도한다.

또한, 저자들은 σ^{(2)}{ab}를 밴드 인덱스 조합에 따라 11+21, 12+22, … 등 6개의 그룹으로 분류하고(식 31), 각 항이 intraband, interband, diamagnetic 등 물리적 프로세스를 어떻게 반영하는지 상세히 분석한다. 특히, intraband 항 중에서 g{ab}^{n}에 의존하는 항이 가장 큰 기여를 하며, 이는 식 (32)–(34)에서 제시된 “기하학적 응답” 형태로 나타난다.

부피·압력과의 연결 고리도 흥미롭다. 평균 밀도 n이 변하면 부피 V도 변해야 하므로 V의 B-전개식(5,6)을 도입하고, n와 V 사이의 관계 n=N/V를 이용해 V^{(1)}/V^{(0)}=−n^{(1)}/n^{(0)} 등을 도출한다. TRS가 깨진 경우, 선형 부피 변화 v^{(1)}_{a}가 Streda와 동일한 계수를 가지며, 이는 “자기 부피 효과”로 해석된다. 압력 역시 동일한 구조를 갖는 전개식으로 표현되어, 압력-자기장 상관관계가 양자 홀 효과와 직접 연결됨을 보여준다.

마지막으로, 저자들은 2밴드 모델 h_{0}=ℏv(k_{x}σ_{y}−k_{y}σ_{x})를 이용해 구체적인 계산을 수행한다. 이 모델에서 FS는 원형이며, f’’_{n}가 FS 근처에 국한되므로, 식 (32)의 기하학적 항이 intraband 응답을 지배한다. 계산 결과 n^{(2)}/n^{(0)}≈10^{-4}% 수준의 밀도 변화를 예측했으며, 보다 정밀한 수치는 실제 결정 구조를 포함한 대규모 수치 시뮬레이션이 필요함을 언급한다.

전체적으로, 이 논문은 Bloch 전자들의 평균 밀도에 대한 1차·2차 자기장 응답을 완전히 게이지 불변적으로 정식화하고, 특히 퀀텀 메트릭이 비선형 응답에 미치는 기하학적 역할을 최초로 명시함으로써, 전자학, 토폴로지 물질, 그리고 고체의 열역학적 특성(부피·압력) 사이의 새로운 연결 고리를 제시한다.