Numerical renormalization group (NRG) is formulated for nonequilibrium steady-state by converting finite-lattice many-body eigenstates into scattering states. Extension of the full-density-matrix NRG for a biased Anderson impurity model, simplified by formulating with the original orbital basis as the Hamiltonian, enables detailed studies of the sub-Kondo spectral evolution in the zero-temperature limit, confirming the double-resonance structure at bias of the Kondo energy scale $T_K$. The distribution shows distinct multi-scale spectral features at energy $Ï$ below the Kondo scale ($Ï\lesssim T_K$) and near the bias ($Ï\gtrsim V$), leading to the nonequilibrium temperature $T_{\rm loc}$ local to the Kondo dot scaling as $k_BT_{\rm loc}\approx V$ for $V\gg T_K$. The current-voltage relation in the low-temperature limit ($T\ll T_K$) deviates from the unitary limit as the bias exceeds the Kondo scale ($V/2\gtrsim T_K$) and reaches the current saturation regime.
논문 ID : 2503.14400제목 : Nonequilibrium Statistics of Biased Kondo Resonance저자 : Jong E. Han (State University of New York at Buffalo)분류 : cond-mat.str-el (강상관 전자계)발표 시간 : 2025년 10월 14일논문 링크 : https://arxiv.org/abs/2503.14400 본 논문은 유한 격자점 다체 고유상태를 산란상태로 변환함으로써 비평형 정상상태에 대한 수치 재정규화군(NRG) 방법을 제시한다. 원시 궤도 기저를 해밀토니안 표현으로 사용하여 전체 밀도행렬 NRG를 편향된 Anderson 불순물 모델로 확장함으로써, 영온도 극한에서 부-콘도 스펙트럼 진화를 상세히 연구할 수 있게 되었다. 콘도 에너지 스케일 T K T_K T K 편향 하에서의 이중 공명 구조를 확인했다. 분포는 콘도 에너지 스케일 이하(ω ≲ T K \omega \lesssim T_K ω ≲ T K )와 편향 근처(ω ≳ V \omega \gtrsim V ω ≳ V )에서 독특한 다중 스케일 스펙트럼 특징을 나타내며, 이는 콘도 양자점의 비평형 온도 T l o c T_{\rm loc} T loc 이 V ≫ T K V \gg T_K V ≫ T K 일 때 k B T l o c ≈ V k_BT_{\rm loc} \approx V k B T loc ≈ V 로 스케일링되도록 한다. 낮은 온도 극한(T ≪ T K T \ll T_K T ≪ T K )에서 편향이 콘도 에너지 스케일을 초과할 때(V / 2 ≳ T K V/2 \gtrsim T_K V /2 ≳ T K ), 전류-전압 관계는 단일성 극한에서 벗어나 전류 포화 영역에 도달한다.
핵심 과제 : 편향된 양자계의 비평형 통계역학 문제를 푸는 것은 지난 수십 년간 거대한 도전 과제였으며, 특히 강상관 양자 불순물 모델의 경우 더욱 그렇다.중요성 :콘도 문제는 강상관 비평형 기술의 시험대로서 양자 다체계의 비평형 거동을 이해하는 데 필수적이다 양자 불순물 모델은 동역학 평균장 이론을 통해 응축물질 극한을 푸는 데 핵심적인 역할을 한다 기존 방법의 한계 :비평형 불순물 해석기의 이론적 엄밀성과 수치 효율성이 만족스러운 수준에 도달하지 못했다 계산의 어려움으로 인해 편향 하에서 콘도 단중항이 어떻게 이중 피크 부-콘도 영역으로 분해되는지에 대한 이해가 제한된다 모델의 스펙트럼 특징에 비해 스펙트럼의 비평형 통계에 대한 이해가 부족하다 연구 동기 :성숙한 NRG 방법을 기반으로 접근 가능한 정상상태 비평형 확장 구축 편향 하에서 콘도 공명 진화에 대한 명확한 이해 제공 복잡한 불순물 상호작용에 적용 가능한 일반적 방법 개발 방법론적 혁신 : 유한 격자점 고유상태를 산란상태로 변환하여 정상상태 비평형 문제를 처리하는 새로운 비평형 NRG 방법 제시기술적 단순화 : 원시 궤도 기저 표현을 사용하여 해밀토니안을 표현함으로써 구현을 단순화하고 일반화된 문제로의 확장에 더 적합하게 함물리적 발견 :V / 2 ≈ T K V/2 \approx T_K V /2 ≈ T K 에서의 콘도 피크 분열 거동 확인다중 스케일 분포 특징 및 점유수 역전 현상 발견 비평형 온도의 스케일링 관계 T l o c ≈ V T_{\rm loc} \approx V T loc ≈ V 규명 수송 특성 : 낮은 온도 극한에서의 전류 포화 현상 검증, T K ≲ V / 2 ≲ 3 T K T_K \lesssim V/2 \lesssim 3T_K T K ≲ V /2 ≲ 3 T K 영역에서 확인해밀토니안은 다음과 같이 정의된다:
H ^ = ∑ σ ϵ d n d σ + U ( n d ↑ − 1 2 ) ( n d ↓ − 1 2 ) + ∑ α k σ ϵ k c α k σ † c α k σ − t Ω ∑ α k σ ( d σ † c α k σ + c α k σ † d σ ) \hat{H} = \sum_\sigma \epsilon_d n_{d\sigma} + U\left(n_{d\uparrow} - \frac{1}{2}\right)\left(n_{d\downarrow} - \frac{1}{2}\right) + \sum_{\alpha k\sigma} \epsilon_k c^\dagger_{\alpha k\sigma} c_{\alpha k\sigma} - \frac{t}{\sqrt{\Omega}} \sum_{\alpha k\sigma} (d^\dagger_\sigma c_{\alpha k\sigma} + c^\dagger_{\alpha k\sigma} d_\sigma) H ^ = ∑ σ ϵ d n d σ + U ( n d ↑ − 2 1 ) ( n d ↓ − 2 1 ) + ∑ α kσ ϵ k c α kσ † c α kσ − Ω t ∑ α kσ ( d σ † c α kσ + c α kσ † d σ )
여기서 α = L , R \alpha = L,R α = L , R 은 소스 및 드레인 저장소를 나타내고, U U U 는 쿨롱 반발 매개변수, ϵ d \epsilon_d ϵ d 는 불순물 에너지 준위이다.
정상상태 비평형 밀도행렬은 다음과 같이 구성된다:
ρ ^ n e q = 1 Z n e q e − β ( H ^ − V Y ^ ) \hat{\rho}_{\rm neq} = \frac{1}{Z_{\rm neq}} e^{-\beta(\hat{H} - V\hat{Y})} ρ ^ neq = Z neq 1 e − β ( H ^ − V Y ^ )
여기서 Hershfield의 Y 연산자는 다음과 같이 정의된다:
Y ^ = 1 2 ∑ k σ ( ψ L k σ † ψ L k σ − ψ R k σ † ψ R k σ ) \hat{Y} = \frac{1}{2}\sum_{k\sigma}(\psi^\dagger_{Lk\sigma}\psi_{Lk\sigma} - \psi^\dagger_{Rk\sigma}\psi_{Rk\sigma}) Y ^ = 2 1 ∑ kσ ( ψ L kσ † ψ L kσ − ψ R kσ † ψ R kσ )
핵심 혁신은 다음 관계식을 통해 Y 연산자를 구성하는 것이다:
Y ^ = i η − L + i η Y ^ 0 \hat{Y} = \frac{i\eta}{-\mathcal{L} + i\eta}\hat{Y}_0 Y ^ = − L + i η i η Y ^ 0
여기서 L \mathcal{L} L 은 리우빌리안 연산자이고, Y ^ 0 \hat{Y}_0 Y ^ 0 는 불순물이 분리되었을 때의 초기 Y 연산자이다. 에너지 고유상태 표현에서:
⟨ α ∣ Y ^ ∣ β ⟩ = i η − E α + E β + i η ⟨ α ∣ Y ^ 0 ∣ β ⟩ \langle\alpha|\hat{Y}|\beta\rangle = \frac{i\eta}{-E_\alpha + E_\beta + i\eta}\langle\alpha|\hat{Y}_0|\beta\rangle ⟨ α ∣ Y ^ ∣ β ⟩ = − E α + E β + i η i η ⟨ α ∣ Y ^ 0 ∣ β ⟩
Wilson 체인(WC) 시뮬레이션 저장소, 체인 길이 N W C = 100 N_{\rm WC} = 100 N WC = 100 , 재정규화 매개변수 Λ = 1.2 \Lambda = 1.2 Λ = 1.2 환경 궤도 처리: 편향 하에서 환경 궤도의 자유 에너지 변화 고려 밀도행렬 가중치는 Wilson 체인 길이 n = M n = M n = M 에서 피크에 도달하며, D M ≈ V / 2 D_M \approx V/2 D M ≈ V /2 를 만족한다 유지 상태 수 N k e p t = 2000 N_{\rm kept} = 2000 N kept = 2000 z 평균화는 8개의 z 값을 사용하여 Wilson 체인의 이산성 감소 각 NRG 반복의 힐베르트 공간 차원은 16 × N k e p t 16 \times N_{\rm kept} 16 × N kept 강 콘도 영역 : U / Γ = 10 U/\Gamma = 10 U /Γ = 10 , Γ = 0.1 \Gamma = 0.1 Γ = 0.1 , T K = 0.0044 T_K = 0.0044 T K = 0.0044 입자-정공 대칭 : ϵ d = 0 \epsilon_d = 0 ϵ d = 0 입자-정공 비대칭 : ϵ d = U / 4 \epsilon_d = U/4 ϵ d = U /4 온도 : T = 0.1 T K T = 0.1T_K T = 0.1 T K (낮은 온도 극한)재정규화 매개변수 : Λ = 4 \Lambda = 4 Λ = 4 스펙트럼 확대 매개변수 : δ = 0.2 × max ( T K , T , V / 2 ) \delta = 0.2 \times \max(T_K, T, V/2) δ = 0.2 × max ( T K , T , V /2 ) 무한소 매개변수 : η = 0.1 × max ( T K , T , V / 2 ) \eta = 0.1 \times \max(T_K, T, V/2) η = 0.1 × max ( T K , T , V /2 ) 스펙트럼 함수 : ρ Σ ( ω ) \rho_\Sigma(\omega) ρ Σ ( ω ) (자기에너지 기법을 사용하여 계산)분포 함수 : f ( ω ) = G d < ( ω ) / [ 2 π i ρ ( ω ) ] f(\omega) = G^<_d(\omega)/[2\pi i\rho(\omega)] f ( ω ) = G d < ( ω ) / [ 2 πi ρ ( ω )] 비평형 온도 : T l o c 2 = 6 π 2 ∫ − ∞ ∞ ω [ f ( ω ) − Θ ( − ω ) ] d ω T_{\rm loc}^2 = \frac{6}{\pi^2}\int_{-\infty}^{\infty} \omega[f(\omega) - \Theta(-\omega)]d\omega T loc 2 = π 2 6 ∫ − ∞ ∞ ω [ f ( ω ) − Θ ( − ω )] d ω 전류 : 그래픽 기법을 통해 계산비상호작용 공명 에너지 준위 모델에서 방법의 정확성 검증 작은 편향에서 전류는 I = G 0 V I = G_0 V I = G 0 V 를 따르며, 여기서 G 0 = 1 / ( 2 π ) G_0 = 1/(2\pi) G 0 = 1/ ( 2 π ) 는 단일성 양자 전도도이다 큰 편향에서 정확한 관계식 I = ( Γ / π ) tan − 1 ( V / 2 Γ ) I = (\Gamma/\pi)\tan^{-1}(V/2\Gamma) I = ( Γ/ π ) tan − 1 ( V /2Γ ) 를 만족한다 영 편향 : 스펙트럼 함수는 ρ Σ ( 0 ) ≈ 1 / ( π Γ ) \rho_\Sigma(0) \approx 1/(\pi\Gamma) ρ Σ ( 0 ) ≈ 1/ ( π Γ ) 에서 단일성 극한을 만족한다임계 편향 : V / 2 ≈ T K V/2 \approx T_K V /2 ≈ T K 일 때, 콘도 공명은 피크 분열을 나타내기 시작한다높은 편향 : 피크 위치는 ± V / 2 \pm V/2 ± V /2 에 정렬되지만, 선폭은 T K T_K T K 보다 훨씬 크므로 약한 상관으로 간주되어야 한다낮은 편향 : 분포 함수는 페르미-디랙 함수에 가깝다중간 편향 : ∣ ω ∣ ≳ V |ω| \gtrsim V ∣ ω ∣ ≳ V 에서 크고 넓은 피크가 발달한다점유수 역전 : 분열된 콘도 피크 사이(− V / 2 ≲ ω ≲ V / 2 -V/2 \lesssim \omega \lesssim V/2 − V /2 ≲ ω ≲ V /2 )에서 f ′ ( 0 ) < 0 f'(0) < 0 f ′ ( 0 ) < 0 이 나타난다핵심 발견: 높은 편향 극한에서
T l o c ≈ V for V / 2 ≫ T K T_{\rm loc} \approx V \quad \text{for} \quad V/2 \gg T_K T loc ≈ V for V /2 ≫ T K
이는 국소 전자 온도가 편향에 정비례함을 나타낸다.
낮은 온도 극한 (T = 0.1 T K T = 0.1T_K T = 0.1 T K ):작은 편향에서 단일성 극한 G 0 = ( 2 π ) − 1 G_0 = (2\pi)^{-1} G 0 = ( 2 π ) − 1 을 만족한다 V / 2 ≳ T K V/2 \gtrsim T_K V /2 ≳ T K 일 때 단일성 극한에서 벗어나기 시작한다T K ≲ V / 2 ≲ 3 T K T_K \lesssim V/2 \lesssim 3T_K T K ≲ V /2 ≲ 3 T K 영역에서 전류 포화가 나타난다높은 온도 (T = T K T = T_K T = T K ):전류 포화 현상이 나타나지 않는다 전류 포화가 강상관 효과에서 비롯됨을 증명한다 논문은 콘도 문제를 다루기 위한 다양한 비평형 기법을 체계적으로 검토한다:
해석적 이론 : Bethe ansatz 방법그래픽 방법 : Landauer 공식, 섭동론수치 방법 : 실시간 그래픽 시뮬레이션, 허시간 시뮬레이션, 텐서 네트워크 방법재정규화군 : 실시간 RG, 섭동 RG, 함수 RGSNRG 방법 : Anders가 제시한 정상상태 NRG는 비상호작용 산란상태를 사용하여 Y ^ 0 \hat{Y}_0 Y ^ 0 를 구성한다본 방법 : 원시 해밀토니안 기저를 사용하여 더 직접적인 구현을 실현하고 복잡한 불순물 상호작용으로의 확장에 적합하다방법 검증 : Hershfield 산란상태 이론과 NRG 방법을 성공적으로 결합하여 비평형 강상관계에 대한 신뢰할 수 있는 수치 도구를 제공한다물리적 그림 : V / 2 ≈ T K V/2 \approx T_K V /2 ≈ T K 에서의 콘도 피크 분열을 확인하고, 부-콘도 영역의 복잡한 다중 스케일 구조를 규명했다비평형 특징 : 점유수 역전과 비평형 온도의 선형 스케일링 관계를 발견하여 강상관 비평형 상태를 이해하기 위한 새로운 통찰력을 제공한다수송 현상 : 전류 포화 현상을 검증하여 이론 예측을 확인했다계산 비용 : 두 개의 Wilson 체인을 처리해야 하므로 힐베르트 공간이 빠르게 증가하여 유지 상태 수를 제한한다매개변수 의존성 : 방법은 무한소 매개변수 η \eta η 의 선택에 민감하며 신중한 조정이 필요하다온도 제한 : 주로 낮은 온도 영역 T ≪ T K T \ll T_K T ≪ T K 에 적용 가능하다확장 적용 : 복잡한 불순물 상호작용 모델에 적용격자 이론 : 동역학 평균장 이론을 통해 비평형 격자 모델에 적용실험 비교 : 양자점 수송 실험과의 상세한 비교이론적 엄밀성 : Hershfield의 엄격한 이론 틀을 기반으로 하며 수학적 유도가 명확하다방법론적 혁신 : 유한 격자점 고유상태를 산란상태로 변환하는 영리한 방법으로 산란상태의 명시적 구성의 복잡성을 회피한다수치적 신뢰성 : 비상호작용 극한을 통해 검증하여 방법의 정확성을 보장한다물리적 통찰력 : 콘도 공명의 비평형 진화에 대한 깊이 있는 이해를 제공한다계산 복잡도 : 평형 상태 NRG에 비해 계산 비용이 현저히 증가한다매개변수 조정 : η \eta η 매개변수의 선택이 경험적이며 결과의 보편성에 영향을 미칠 수 있다실험 검증 : 실제 양자점 실험과의 직접적인 비교가 부족하다방법론적 기여 : 비평형 강상관계에 대한 새로운 수치 도구 제공이론적 가치 : 콘도 효과의 비평형 거동에 대한 이해 심화응용 전망 : 더 복잡한 불순물 모델 및 격자계로 확장 가능양자점 수송 실험의 이론적 분석 강상관 전자계의 비평형 연구 동역학 평균장 이론의 불순물 해석기 분자 접합 전자공학의 이론적 모델링 논문은 해당 분야의 중요한 문헌을 인용하고 있으며, 다음을 포함한다:
콘도의 원본 연구 및 Hewson의 고전 교과서 Hershfield의 비평형 정상상태 이론 Wilson의 수치 재정규화군 방법 Anders의 정상상태 NRG 방법 다양한 비평형 양자 불순물 이론 및 수치 방법 본 논문은 비평형 양자 다체 이론 분야에서 중요한 기여를 하였으며, 강상관 비평형계를 연구하기 위한 엄격하면서도 실용적인 수치 방법을 제공한다. 특히 콘도 효과의 비평형 거동을 이해하는 데 있어 중요한 가치를 지닌다.