현재 나노스케일 기술 노드에서 열 안정성 장벽을 높이기 위해 이중 계면 자기 터널 접합(MTJ)이 개발되었습니다. 강자성체/중금속 구조 사용으로 인해 이중 계면 MTJ에는 드잘로신스키-모리야 상호작용(DMI)이 필연적으로 존재합니다. 기존 연구들은 DMI가 단일 계면 스핀전달 토크(STT) MTJ에 해로운 영향을 미친다고 보고했습니다. 본 연구에서는 이중 계면 STT-MTJ에서 DMI의 해로운 영향이 거의 제거될 수 있음을 입증했습니다. 이 결론은 루더만-키텔-카스야-요시다(RKKY) 상호작용이 DMI를 억제하는 효과 덕분입니다. 상세 메커니즘을 이론적 모델과 마이크로마그네틱 시뮬레이션 결과로 분석했습니다. 본 연구는 이중 계면 STT-MTJ의 복합 자유층에서 DMI를 적절히 제어하는 것의 중요성을 강조합니다.
서론
자기 랜덤 액세스 메모리(MRAM)는 낮은 전력 소비, 고밀도, 빠른 접근 속도, 준무한 내구성, CMOS 기술과의 좋은 호환성 덕분에 차세대 비휘발성 메모리의 가장 유망한 후보입니다. MRAM의 기본 소자는 자기 터널 접합(MTJ)으로, 두 강자성층(고정층과 자유층) 사이에 터널 장벽이 샌드위치 구조로 되어 있습니다. 수직 이방성 발전 덕분에 MTJ 특징 크기가 40nm 아래, 심지어 1×nm까지 축소되었습니다. 그러나 서브-40nm MTJ의 과제는 적절한 열 안정성 장벽 E = μ₀MₛHₖV/2를 유지하는 것입니다. 이 식에서 V는 자유층 부피이므로 MTJ 축소에 따라 E가 감소하여 데이터 보존 시간이 줄어듭니다. 이 문제를 극복하기 위해 이중 계면 MTJ가 제안되어 서브-40nm 기술 노드에서 충분히 높은 E를 달성했습니다. 복합 자유층으로 두 개의 결합된 강자성층을 사용하면 이중 계면 MTJ에서 등가 부피(V)가 증가하여 열 안정성 장벽이 향상되고, 감쇠 상수가 감소하여 낮은 스위칭 전류를 유지합니다.
이중 계면 MTJ에서 강자성체/중금속(FM/HM) 구조는 성능 최적화에 중요한 역할을 합니다. 한편으로는 FM/HM 구조가 스핀-궤도 결합(SOC)을 증가시켜 수직 이방성을 유도합니다. 다른 한편으로는 중금속이 복합 자유층의 두 강자성층 사이의 스페이서 역할을 하여 루더만-키텔-카스야-요시다(RKKY) 상호작용을 제공하며, 이는 두 층의 자화를 강자성적으로 결합시켜 마치 동일한 층처럼 거동하게 합니다. 또한 최근 연구들은 중금속의 강한 SOC와 강자성체의 원자 스핀이 결합하여 비대칭 교환 결합인 드잘로신스키-모리야 상호작용(DMI)을 형성함을 보였습니다. 따라서 FM/HM 구조를 가진 이중 계면 MTJ에서는 DMI가 자연스럽게 유도됩니다. DMI는 카이랄 자기 텍스처(스핀 나선, 스커미온, 네엘형 도메인 벽)를 선호하며 자화 역학에 극적인 영향을 미칩니다. 이중 계면 MTJ에서는 두 개의 FM/HM 계면과 추가적인 RKKY 상호작용을 함께 고려해야 하므로 DMI의 역할이 더 복잡해집니다. 따라서 이중 계면 MTJ에서 DMI의 영향을 밝히는 것이 중요합니다.
본 논문에서는 처음으로 DMI와 RKKY 상호작용 하에서 이중 계면 MTJ의 스위칭 공정을 연구합니다. 이중 계면 MTJ는 스핀전달 토크(STT)로 스위칭되며, 이는 MRAM 데이터 쓰기의 주류 방식입니다. 최근 DMI가 STT 스위칭에 해로운 영향을 미친다는 보고가 있었습니다. 본 연구 결과는 이중 계면 MTJ에서 RKKY 상호작용이 DMI의 해로운 영향을 억제하여 빠른 스위칭과 더 균일한 역학을 가져옴을 보여줍니다. 본 연구는 이중 계면 STT-MTJ가 부정적 계면 효과에 대해 강건함을 입증합니다.
실험 방법
본 연구에서 다룬 소자는 그림 1a에 도시된 바와 같이 FM/HM/FM 구조를 복합 자유층으로 가집니다. HM 층 두께는 유도된 RKKY 상호작용이 두 FM 층을 강자성적으로 결합시키도록 적절히 조절됩니다. 한 FM 층은 자석적으로 더 부드러워 FL1(자유층 1)로, 다른 층은 더 단단하여 FL2(자유층 2)로 표기합니다. 이중 계면 MTJ에 전류를 가해 복합 자유층의 자화를 스위칭하며 STT를 생성합니다. 본 연구에서는 기준층에서 FL1로 전달되는 STT만 고려하고 FL1과 FL2 사이의 다른 토크는 무시합니다. 이 단순화된 모델은 선행 연구들과 일치합니다. DMI는 두 FM/HM 및 HM/FM 계면에서 유도되며, 서로 다른 카이랄리티로 인해 반대 부호를 가집니다.
a 연구 소자 개략도(기타 층 생략). b 시간 의존적 mz(단위 자화의 수직 성분) 대표 결과. Case A: σ=1×10⁻³ J/m², D₁=D₂=0(FL1 빨강, FL2 파랑). Case B: σ=1×10⁻³ J/m², D₁=1 mJ/m², D₂=-1 mJ/m²(FL1 주황, FL2 청록). Case C: σ=1×10⁻⁴ J/m², D₁=D₂=0(FL1 초록, FL2 검정)
이중 계면 MTJ에서 FL1과 FL2의 자화 역학은 마이크로마그네틱 시뮬레이션으로 연구했습니다. 단위 자화 벡터의 시간 진화는 다음 란다우-리프시츠-길버트(LLG) 방정식으로 지배됩니다. 오픈소스 마이크로마그네틱 시뮬레이션 도구인 OOMMF 패키지를 사용하여 소자 구조를 모델링하고 LLG 방정식을 풀어 자화 역학을 분석했습니다.
RKKY 에너지는 자기 모멘트 쌍 mᵢ와 mⱼ 사이에 Eᵢⱼ = σ[1 - mᵢ·mⱼ]/Δᵢⱼ로 정의됩니다. 여기서 σ는 두 표면 사이의 쌍선형 표면 교환 계수입니다. FL1과 FL2가 강자성적으로 결합하면 σ > 0으로 RKKY 상호작용은 mᵢ를 mⱼ와 평행하게 만들려는 경향이 있습니다. DMI는 큰 SOC를 가진 제3의 원자를 통해 인접 원자 스핀 S₁과 S₂에 작용하며, 해밀토니안은 H_DM = -D₁₂·(S₁×S₂)로 표현됩니다. 따라서 DMI는 S₁과 S₂ 사이의 균일성을 저하시키며 RKKY 상호작용과 경쟁합니다.
결과 및 고찰
먼저 시간 의존적 mz(단위 자화의 수직 성분) 대표 시뮬레이션 결과를 그림 1b에 보입니다. RKKY 상호작용이 충분히 강하면(예: σ=1×10⁻³ J/m², Case A와 B), FL1과 FL2가 결합되어 DMI 유무와 관계없이 자화 역학이 거의 동일합니다. DMI 도입이 자화 스위칭 과정을 왜곡시키는 것도 확인됩니다(Case B). 이는 선행 결과와 일치하며 DMI의 비대칭 교환 때문으로 보입니다. RKKY 상호작용이 충분히 강하지 않으면 FL1과 FL2의 자화 역학이 이상적으로 결합되지 않아 유의미한 차이가 관찰됩니다(Case C). 아래 시뮬레이션 결과는 별도 언급 없이는 충분히 강한 RKKY 상호작용 하에서 얻었습니다.
이후 다양한 RKKY 상호작용 하에서 스위칭 속도를 연구했습니다. 스위칭 속도는 mz가 0에 도달하는 시간(스위칭 시간)으로 반영됩니다. D₁과 D₂는 각각 양수와 음수 값으로 설정했습니다. 그림 2에 결과가 나와 있습니다. DMI가 없으면 스위칭 시간은 RKKY 상호작용 강화에 따라 증가하며, 이는 선행 결과와 일치합니다. 더 강한 RKKY 상호작용은 FL1과 FL2의 자화 역학을 더 응집되게 하여 복합 자유층의 등가 이방성을 효과적으로 증가시키기 때문입니다. 그러나 DMI가 존재하면 스위칭 시간의 RKKY 강도 의존성이 더 혼돈스러워집니다. 이 혼돈은 주로 FL1과 FL2 사이의 이방성 불일치에 기인하며 추후 설명하겠습니다. 이 결과들은 이중 계면 MTJ의 스위칭 거동에 DMI의 영향이 무시할 수 없음을 입증합니다.
RKKY 강도에 따른 스위칭 시간(σ는 로그 스케일). D₁과 D₂는 같은 값이나 부호가 반대
다음으로 DMI 효과를 더 자세히 연구했습니다. 그림 3은 DMI 강도에 따른 스위칭 시간을 보여줍니다. 그림 3a에서는 D₁과 D₂를 의도적으로 같은 양수로 설정했습니다(현실에서는 부호가 반대). 이는 시뮬레이션 모델 검증을 위한 가상 사례입니다. 물리 이론 관점에서 두 양수 DMI의 해로운 효과는 강자성적으로 결합된 RKKY 상호작용 하에서 누적됩니다. 따라서 D₁과 D₂ 증가에 따라 스위칭 시간이 증가할 것으로 예상되며, 이는 선행 연구와 일치합니다. 그림 3a 결과와 잘 일치하여 시뮬레이션 모델의 타당성이 검증되었습니다. 반면 그림 3b에서처럼 D₁과 D₂가 반대 부호면 DMI 해로운 효과가 완화되어 스위칭 시간 변동이 그림 3a보다 훨씬 작습니다. 그림 3b 곡선이 정확히 단조롭지는 않은데, 국소적 변동은 나중에 설명하겠습니다. 주목할 점은 D₁과 D₂ 크기를 적절히 조절하면 두 계면의 DMI 효과가 상쇄될 수 있다는 것입니다(그림 3c). 이는 카이랄리티 이론으로 설명 가능합니다.
DMI 강도에 따른 스위칭 시간. a D₁, D₂ 같은 양수. b D₁, D₂ 같은 크기이나 부호 반대. c t₁D₁+t₂D₂=0 조건. d 두께나 이방성 상수 변경하면서 t₁D₁+t₂D₂=0 유지. 파란선: t₁=2nm 변경. 빨간선: t₁=1.5nm 변경. 삼각형 데이터: σ=3×10⁻³ J/m². 원 데이터: σ=1×10⁻³ J/m²
DMI 에너지는 E_DM = t ∬ D[mₓ(∂mz/∂x) - mz(∂mₓ/∂x) + mᵧ(∂mz/∂y) - mz(∂mᵧ/∂y)] d²r = tDε_DM로 표현됩니다. 여기서 D는 연속 DMI 상수, t는 강자성층 두께입니다. 앞서 언급했듯 충분히 강한 RKKY 상호작용 하에서 FL1과 FL2의 자화 역학은 거의 동일하므로 같은 ε_DM을 얻습니다. 그러면 FL1과 FL2의 총 DMI 에너지는 E_tot = (t₁D₁ + t₂D₂)ε_DM으로 계산됩니다. 따라서 D₁/D₂ = -t₂/t₁로 설정하면 큰 σ에서 두 층의 DMI 효과를 완전히 상쇄할 수 있으며, 이는 그림 3c와 일치합니다. 이 결론은 그림 3d의 추가 결과로 더욱 검증되었습니다.
복합 자유층의 등가 DMI 크기(D_eq)는 D_eq = (t₁D₁ + t₂D₂)/(t₁ + t₂)로 표현되며, 이는 이중 계면 MTJ에서 DMI 효과를 정량적으로 분석하는 데 사용됩니다. 이 식의 유효성을 검증하기 위해 그림 4a에 두 그룹의 시뮬레이션 결과를 보입니다. 두 곡선은 같은 D_eq지만 다른 {D₁, D₂} 쌍으로 얻었습니다. 두 곡선 사이에 약간의 차이가 있지만 전체 경향은 유사하며 DMI가 STT 스위칭에 해로운 영향을 미침을 검증합니다. 이 차이는 FL1과 FL2가 다른 이방성 상수를 가져 자화 역학의 국소적 불확실 진동을 일으키기 때문입니다(그림 4c). 같은 현상이 그림 2와 3b에서도 관찰됩니다. 반면 그림 4b,d는 FL1과 FL2의 이방성 상수를 같게 설정한 이상적 경우로, 두 곡선이 잘 일치하여 위 D_eq 표현이 이중 계면 MTJ의 등가 DMI 효과를 잘 설명함을 나타냅니다.
a,b D_eq에 따른 스위칭 시간. 각 D_eq는 D_eq=(t₁D₁+t₂D₂)/(t₁+t₂)에 따라 두 개의 서로 다른 {D₁,D₂} 쌍으로 얻음. 빨간 곡선: D₁ 변화, D₂=1 mJ/m² 고정. 파란 곡선: D₁, D₂ 항상 같은 값. σ=1×10⁻² J/m². a 다른 파라미터는 표 1 설정. b Kᵤ₁=Kᵤ₂=0.7 mJ/m³로 이상적 경우. c,d 각각 a,b에 해당하는 시간 의존적 mz 대표 결과
마지막으로 자화 역학의 시간 진화를 더 자세히 분석했습니다. 그림 5는 자화 스위칭 중 시간 의존적 에너지를 보여줍니다. FL1과 FL2의 DMI 에너지는 D₁, D₂의 부호와 크기에 따라 누적되거나 상쇄됩니다. 이 경향은 위 이론 모델과 잘 일치합니다. 또한 RKKY 에너지는 낮은 값으로 유지되어 FL1과 FL2의 자기 모멘트가 동기적으로 구동됨을 확인합니다. RKKY와 DMI 필드 분포는 그림 6에 나와 있으며, RKKY 필드는 경우에 따라 다른 역할을 합니다. 첫째, 0이 아닌 DMI 경우(Case 2, 3)에서 RKKY 필드가 0 DMI 경우(Case 1)보다 훨씬 강합니다. DMI 존재 시 자기 텍스처의 추가적 비균일성을 극복해야 하기 때문으로 이해됩니다. 둘째, D₁과 D₂가 반대 부호면 RKKY 필드가 FL1과 FL2 모두에서 DMI 필드에 저항합니다(Case 2). 결과적으로 DMI가 약화되어 자화 역학이 더 균일해집니다. 반대로 D₁과 D₂가 같은 부호면 RKKY 필드가 한 강자성층에서는 DMI 필드에 저항하지만 다른 층에서는 돕습니다(Case 3). 따라서 전체 DMI 필드가 여전히 자화 역학에 일정한 영향을 미쳐, D₁과 D₂가 같은 부호면 DMI가 상쇄될 수 없음을 검증합니다.
시간 진화에 따른 DMI 및 RKKY 에너지. a D₁=1.5 mJ/m², D₂=-1 mJ/m²(DMI 효과 상쇄). b D₁=D₂=1 mJ/m²(DMI 효과 누적). c D₁=1 mJ/m², D₂=-1 mJ/m²(DMI 효과 완화되나 상쇄 안됨). d D₁=D₂=0DMI 및 RKKY 필드 공간 분포. 각 경우의 대표적 시간 순간 결과. Case 1: D₁=D₂=0. Case 2: D₁=1.5 mJ/m², D₂=-1 mJ/m²(DMI 효과 상쇄). Case 3: D₁=D₂=1 mJ/m²(DMI 효과 누적)
그림 7은 자화 스위칭 중 FL1과 FL2의 마이크로마그네틱 구성을 보여줍니다. 모든 경우에서 도메인 벽이 나타나지만, 특정 시간 순간에 다른 특징을 관찰할 수 있습니다. DMI가 비균일 자기 텍스처를 선호하는 것으로 알려져 있지만, 그림 7에서 DMI 효과가 상쇄되면 DMI 존재 하에서도 균일한 자화가 형성됩니다(mz=-0.5 시점, Case 2). 이는 위 이론 모델을 다시 검증합니다. 또한 D₁과 D₂가 같은 부호면 자화 역학이 더 비균일해집니다(Case 3에서 도메인 벽이 항상 나타남). 더 작은 MTJ로 시뮬레이션한 결과도 보입니다(그림 7 마지막 두 행). Case 2(DMI 상쇄)와 Case 3(DMI 미상쇄) 간 마이크로마그네틱 구성 차이가 더 두드러집니다.
자화 스위칭 중 마이크로마그네틱 구성. Case 1-3은 그림 6과 같은 파라미터로 구성
결론
본 연구는 이중 계면 STT-MTJ에서 DMI 효과를 포괄적으로 연구했습니다. 알려진 바와 같이 이중 계면 MTJ는 열 안정성 장벽 향상을 위해 개발되었습니다. 본 연구에서는 이중 계면 MTJ의 또 다른 장점, 즉 DMI의 해로운 효과를 억제할 수 있음을 입증했습니다. 두 강자성층의 DMI는 적절한 값과 반대 부호로 구성되면 억제되거나 상쇄될 수 있으며, 이는 이중 계면 STT-MTJ 구조에서 자연스럽게 만족됩니다. 결론을 설명하기 위한 이론적 모델을 제안하고, 자화 역학에서 DMI가 수행하는 역할을 밝히기 위해 마이크로마그네틱 시뮬레이션 결과를 논의했습니다. 본 연구는 이중 계면 STT-MTJ에서 DMI를 최소화하는 실현 가능한 접근법을 제공합니다.