은·금 나노구조를 위한 다상 구아시 포텐셜 파라미터화

초록

본 연구는 밀도범함수이론(DFT)으로 구한 은(Ag)과 금(Au)의 FCC, BCC, SC, 다이아몬드형 및 다이머 구조의 결합에너지 곡선을 이용해 구아시(Gupta) 포텐셜을 다중상태에 걸쳐 재파라미터화하였다. 파라미터의 전이성을 검증하기 위해 55원자 Ag 클러스터에 대해 2 ns 분자동역학(MD) 시뮬레이션을 수행했으며, 초기 FCC 구조가 실온에서 icosahedral 구조로 변환되는 것을 확인하였다. 또한, 새롭게 얻은 Ag 파라미터를 사용해 500 ± 50 K의 증착 온도에서 1000원자까지의 나노클러스터를 원자별 증착 방식으로 성장시켰고, 이 온도 구간이 가장 낮은 에너지 구조를 생성함을 제시한다.

상세 분석

이 논문은 전통적인 구아시 포텐셜이 특정 결정구조에만 최적화된 한계를 극복하고자, 여러 결정상과 저차원 구조에 대한 DFT 기반 에너지 데이터를 동시에 활용한 파라미터화 전략을 제시한다. 먼저, VASP와 같은 고성능 DFT 코드를 이용해 Ag와 Au의 FCC, BCC, SC, 가상의 다이아몬드형, 그리고 두 원자 다이머에 대해 전자밀도와 총 에너지를 광범위한 결합거리 범위(≈2.0–4.0 Å)에서 계산하였다. 각 구조별 에너지-거리 곡선은 포텐셜 파라미터인 A, ξ, p, q, r₀를 비선형 최소제곱법으로 동시에 피팅함으로써, 단일 구조에 국한되지 않는 전역 최적화를 달성하였다. 특히, 다이머와 다이아몬드형 구조는 표면 원자 비율이 높아 표면 에너지와 다체 상호작용을 강하게 반영하므로, 파라미터가 나노스케일 클러스터와 촉매 표면에 적용될 때의 전이성을 검증하는 중요한 기준이 된다.

피팅 결과, 기존 문헌에 보고된 Gupta 파라미터와 비교했을 때, 특히 p와 q 값이 약 10 % 정도 조정되어 전자구름의 비선형 응집 효과를 더 정확히 재현한다는 점이 눈에 띈다. 이러한 조정은 Ag와 Au 모두에서 FCC와 BCC 사이의 에너지 차이를 보다 현실적으로 재현하며, SC와 다이아몬드형 같은 메타안정 구조에서도 과도한 에너지 오차를 방지한다.

전이성을 검증하기 위해 55원자 Ag 클러스터(icosahedral 구조가 가장 안정적인 것으로 알려짐)를 초기 FCC 격자 형태로 배치하고, 300 K에서 2 ns(2 ps 타임스텝) MD 시뮬레이션을 수행하였다. 시뮬레이션 결과, 클러스터는 초기 FCC 배열이 급격히 재배열되어 icosahedral 대칭을 획득함을 확인했으며, 이는 새로운 파라미터가 표면 원자 재배열 에너지 장벽을 정확히 포착함을 의미한다.

또한, Ag 원자들을 1 Å 간격으로 순차적으로 증착하면서 온도(300–700 K)를 변화시킨 원자별 성장 시뮬레이션을 수행하였다. 에너지 최소화와 구조 분석 결과, 500 ± 50 K 구간에서 형성된 클러스터는 평균 결합거리와 전자밀도 분포가 가장 균일하고, 결함 밀도가 최소이며, icosahedral 또는 다중-다면체 형태를 선호한다. 이는 실험적 증착 및 어닐링 공정에서 최적 온도 범위가 약 500 K임을 이론적으로 뒷받침한다.

전반적으로, 다중상태 DFT 데이터를 기반으로 한 Gupta 파라미터화는 전통적인 단일상태 피팅 대비 전이성(transferability)이 크게 향상되었으며, 나노클러스터 성장, 촉매 표면 재구성, 그리고 금속 나노와이어 설계 등 다양한 나노스케일 응용에 바로 활용 가능함을 보여준다.