스핀 의사밸브 트릴레이어에서 스핀 전송 및 터널 자기저항 분석

초록

본 연구는 FeCr, Fe, Co, NiFe, Ni 전극과 GaSb, InSb, InAs, GaAs, ZnSe 반도체 장벽을 조합한 125가지 PSV 트릴레이어 구조의 스핀 전송을 Slonczewski 모델과 Dresselhaus·Rashba SOC를 포함해 저온에서 Landauer‑Büttiker 단일 채널 방식으로 계산하였다. 최대 TMR 83.60%는 Fe₉₀Cr₁₀/GaSb/Fe₉₀Cr₁₀ 조합에서 얻었으며, Dresselhaus SOC가 Rashba보다 TMR에 더 큰 영향을 미친다. 전극 순서를 바꾸면 Fermi 에너지 차이로 TMR이 변한다.

상세 분석

본 논문은 전통적인 스핀 밸브 구조를 확장해 전극을 Fe₉₀Cr₁₀, Fe, Co, NiFe, Ni 다섯 종류로, 장벽을 III‑V (GaSb, InSb, InAs, GaAs)와 II‑VI (ZnSe) 반도체 다섯 종류로 설정함으로써 총 125개의 조합을 체계적으로 분석한다. 이론적 기반은 Slonczewski가 제시한 스핀 전류 토크 모델이며, 전극 내부의 교환 분열 Δₓₛ와 반도체 내부의 Dresselhaus(β)·Rashba(α) 스핀‑궤도 결합을 명시적으로 포함한다. Hamiltonian은 전극 영역에서 자유 전자와 교환 분열을, 장벽 영역에서는 유효 질량 m*₍b₎와 SOC 항을 각각 포함한 2×2 파울리 행렬 형태로 기술된다.

전극‑장벽 경계에서는 파동함수와 질량 가중 미분 연속성을 적용해 전송 및 반사 계수를 구하고, 이를 통해 전송 확률 Tₛ(k∥,θ) 를 도출한다. 전송 확률은 복소수 파라미터 ρₑₛ, λₑₛ 등으로 표현되며, Dresselhaus와 Rashba가 각각 k_z²·β와 α·sin(2φ) 형태로 전자 스핀을 회전시키는 효과를 만든다. 이때 θ는 전극 간 자기화 방향 차이(평행 θ=0, 반평행 θ=π)를 의미한다.

Landauer‑Büttiker 공식 G = (e²/ħ)∫T(k∥,θ)dk∥ 를 사용해 평행(P) 및 반평행(AP) 전도도를 계산하고, TMR = (G_P−G_AP)/G_AP 로 정의한다. 계산은 0 K 근사에서 단일 채널(k_z 고정)으로 수행했으며, k∥ 범위는 0–1 nm⁻¹, 장벽 두께 a는 1–6 nm로 스캔한다. 전극별 Fermi 파동수 k_F↑, k_F↓와 교환 분열 Δₓₛ는 실험값 및 문헌 보간을 통해 설정했으며, 반도체별 유효 질량, 밴드갭, SOC 상수 β, α는 최신 데이터베이스를 참조했다.

주요 결과는 다음과 같다. (1) TMR 최대값 83.60%는 Fe₉₀Cr₁₀/GaSb/Fe₉₀Cr₁₀ 구조에서 a≈1.92 nm일 때 도출되었으며, 이는 전극의 높은 Δₓₛ와 GaSb의 중간 장벽 높이(0.75 eV) 및 큰 Dresselhaus 상수 β가 결합한 결과이다. (2) Dresselhaus SOC가 Rashba보다 TMR에 미치는 기여도가 현저히 크다; β가 큰 GaSb와 InSb에서 TMR 상승폭이 가장 크게 나타난다. (3) 전극의 자기화 방향과 결정축(