다중스케일 모델링으로 CBRAM 셀 전류‑전압 특성 예측

초록

본 논문은 금속/산화물/금속 구조를 갖는 전도성 브리징 랜덤 액세스 메모리(CBRAM) 셀의 I‑V 특성을 다중스케일 시뮬레이션 프레임워크로 예측한다. 유한요소법(FEM) 모델에 필요한 물성값을 ab initio 계산과 머신러닝 기반 경험적 모델에서 추출해 피팅 파라미터를 최소화하고, 전기·열·이온 전달을 동시에 고려한다. Ag/a‑SiO₂/Pt 구조에 적용해 실험 데이터와 높은 일치를 보였으며, 얇은 나노필라멘트에서의 줄열 효과도 정량적으로 분석하였다.

상세 분석

이 연구는 CBRAM 셀의 전기·화학·열적 거동을 원자 수준에서부터 매크로 스케일까지 일관되게 연결하는 다중스케일 시뮬레이션 체인을 구축한다. 첫 단계에서는 밀도범함수이론(DFT) 기반의 ab initio 분자동역학(AIMD)과 네들드 엘라스틱 밴드(NEB) 계산을 이용해 Ag 이온의 확산 장벽, 산화·환원 반응의 전위차, 그리고 금속 필라멘트의 전자·포논 전도도를 얻는다. 여기서 얻은 원자‑레벨 파라미터는 반경 1 nm 이하의 나노필라멘트 형성·소멸 메커니즘을 정확히 기술하는 데 필수적이다.

다음으로, 이러한 파라미터를 머신러닝으로 보강한 경험적 포스필드에 입력해 장시간·대규모 시뮬레이션을 수행한다. 특히, Ag‑a‑SiO₂ 계면에서의 전하 전달과 이온 이동을 설명하기 위해 부틀러‑볼머식 전기화학 반응 속도식을 도입하고, 전기장에 의한 이온 이동을 Nernst‑Planck 방정식으로 모델링한다. 전자 전송은 비균일 전위 분포를 고려한 터널링 전류 모델을 사용해 필라멘트가 전극 간격을 좁히는 순간 급격히 증가하는 전류를 재현한다.

핵심은 유한요소법(FEM) 기반의 연속체 모델에 원자‑레벨 파라미터를 직접 매핑함으로써, 전기·열·기계적 연동을 완전 자가일관적으로 해결한다는 점이다. 열 전달은 전자·포논 전도도와 필라멘트의 저항 발열을 포함해 푸아송 방정식으로 풀며, 필라멘트 직경이 5 nm 이하로 감소할 경우 10 µA 수준 전류에서도 온도가 수십 켈빈 상승한다는 결과를 도출한다. 이러한 줄열 현상은 RESET 단계에서 필라멘트 용해를 촉진해 전류‑전압 히스테리시스에 결정적 영향을 미친다.

마지막으로, Ag/a‑SiO₂/Pt 구조에 대한 실험 I‑V 곡선과 시뮬레이션 결과를 비교했을 때, 전압 설정값(V_set)과 전류 제한값(I_cc)을 제외하고는 거의 파라미터 조정 없이 높은 일치를 보였다. 이는 전통적인 KMC·연속체 모델에서 요구되는 10여 개 이상의 피팅 파라미터를 대폭 감소시킨 것으로, 새로운 소재·구조 설계 시 사전 예측 능력을 크게 향상시킨다.