THz 광학으로 보는 샤이바 상태 유도 무갭 초전도 현상

초록

본 논문은 자기불순물에 의해 형성된 샤이바(Y‑SR) 결합 상태가 전자쌍을 약화시켜 전통적인 s‑파 초전도체를 무갭(gapless) 상태로 전이시킬 수 있음을, 완전 자기일관적인 수치 계산을 통해 입증한다. 무갭 상태에서는 단일입자 에너지 갭이 사라지지만 초전도 순서 매개변수 Δ₀는 남아 있다. 결과적으로 저주파(ℏω < 2|Δ₀|) 영역에서 선형 THz 전도도가 유한한 흡수를 보이며, 비선형 응답은 힉스 모드에 의해 지배되어 여전히 코히런트한다. 이는 단일입자 소산과 집단적 응집체 동역학 사이의 근본적인 차이를 드러내며, THz 분광법이 무갭 초전도체를 탐지하는 새로운 수단이 될 수 있음을 시사한다.

상세 분석

이 연구는 s‑wave 초전도체에 무작위로 배치된 자기불순물(밀도 n_i, 교환 상수 J)을 포함한 전자‑스핀 상호작용 모델을 시작점으로 삼는다. Nambu‑스핀 공간에서의 전자 연산자를 이용해 해밀토니안을 구성하고, Dyson 방정식과 무작위 위상 근사(RPA)를 통해 자기불순물에 의한 셀프‑에너지 Σ(ω)를 도출한다. 핵심은 Shiba가 제시한 복소화된 주파수 ˜ω와 갭 ˜Δ₀를 이용한 자기일관적 정규화 방정식(식 9‑10)이다. 여기서 γ_s는 교환 상호작용에 의한 준입자 이완률을, η는 단일 샤이바 결합 상태의 에너지 분열을 나타낸다.

복소 방정식의 분기점과 브랜치 컷을 처리하기 위해 저자들은 식 12의 6차 다항식 형태로 변환하였다. 이 다항식은 여섯 개의 복소 근을 제공하지만, 물리적으로 의미 있는 해는 retarded Green 함수의 조건(Im y ≥ 0)과 무불순물 한계(y ≈ x)에서 연속성을 유지하는 근 하나로 제한된다. 이러한 선택 기준은 수치적 안정성을 보장하고, 전체 파라미터 공간(온도 T, 불순물 농도 n_i)에서 수렴을 확보한다.

자기일관적 루프는 (i) ˜ω/˜Δ₀를 구하고, (ii) 이를 식 10에 대입해 ˜Δ₀와 ˜ω를 얻으며, (iii) 식 13의 비정상 그린 함수를 통해 anomalous Green 함수 F_k를 계산한다. 최종적으로 식 14의 gap 방정식을 이용해 실제 초전도 순서 매개변수 Δ₀를 결정한다. 이 과정은 전자 스펙트럼, 전도도, 비선형 응답 등 다양한 물리량을 일관되게 산출한다.

결과적으로, 불순물 농도가 임계값을 초과하면 η·Δ₀ 중심의 Y‑SR 밴드가 넓어져 단일입자 밀도 상태 ρ(ω)가 갭 내부까지 확장된다. Δ₀는 여전히 유한하지만, 단일입자 에너지 갭이 사라지는 무갭 초전도 상태가 형성된다. 선형 THz 전도도 σ₁(ω)는 전통적인 2|Δ₀| 임계 이하에서도 유한한 흡수를 보이며, 이는 ρ(ω) ∝ Im G(ω) ∝ γ_s와 직접 연결된다. 반면 비선형 THz 펄스(3ω, 5ω 등)에서 관측되는 힉스 모드 진동은 여전히 ℏω_H ≈ 2|Δ₀|에 집중되고, 감쇠가 거의 없으며, 이는 무갭 상태에서도 초전도 응집체의 위상 일관성이 유지됨을 의미한다.

이러한 차이는 실험적으로 THz 시간분해 광학(THz‑TDS)과 비선형 펌프‑프로브 기법을 통해 구분 가능하며, 특히 무접점(non‑contact) 방식으로 전자 스펙트럼과 집단 모드 동역학을 동시에 탐지할 수 있는 새로운 검출 스킴을 제공한다. 또한, Shiba 밴드가 토폴로지적 밴드 구조를 형성할 경우 마요라나 모드와의 결합을 통한 양자계산 응용 가능성도 시사한다.