에너지의 변형률 미분으로 셀 힘과 선형 탄성 응답을 계산하기
원자 힘이 원자 좌표에 대한 에너지 기울기라면, 셀을 늘리고 비트는 변형률에 대한 기울기는 응력입니다. 작은 변형 주변의 곡률은 탄성계수를 줍니다.
원자 스케일 총에너지에서 거시적 탄성상수와 기계적 안정성을 어떻게 얻을까요?
이번 장의 답응력은 셀 부피당 에너지의 변형률 미분이고 탄성텐서는 응력의 선형 기울기입니다. 대칭에 맞는 독립 변형을 ±로 적용하고 내부 원자를 완화해 수렴된 응답을 fit하며 안정조건을 확인합니다.
01에너지 미분으로서의 stress tensor
작은 변형률 $\varepsilon_{ij}$에 대한 총에너지 변화는 응력텐서와 연결됩니다. 코드마다 인장·압축 부호 convention이 다를 수 있으므로 출력 정의를 확인해야 합니다.
외압 0의 셀 최적화에서는 원자 힘과 응력이 모두 허용치 아래여야 합니다. 원자만 완화하고 셀을 고정하면 constrained 구조의 응력이 남으며, 이것이 의도한 epitaxial strain 모델인지 실수인지 구분해야 합니다.
02탄성상수를 얻는 두 경로
strain-stress 방법은 작은 ±변형을 가하고 계산된 응력을 선형 fit합니다. strain-energy 방법은 에너지의 2차 곡률을 fit합니다. 응력은 한 계산에서 여러 성분을 주지만 수치 잡음에 민감하고, 에너지는 안정적이지만 2차 작은 차이를 정확히 맞춰야 합니다.
- 평형 셀과 엄격한 stress 수렴을 먼저 확보
- 약 ±0.5~1.5%의 여러 strain으로 선형·2차 영역 검사
- clamped-ion과 relaxed-ion 탄성을 구분
- 양·음 strain을 함께 써 비선형 편향 제거
03Born 기계적 안정조건
안정한 결정에서는 임의의 작은 strain에 대한 2차 에너지 변화가 양수여야 하므로 탄성 행렬이 양의 정부호여야 합니다. 입방정에서는 $C_{11}-C_{12}>0$, $C_{11}+2C_{12}>0$, $C_{44}>0$ 같은 조건으로 단순화됩니다.
기계적 안정은 dynamical stability와 다릅니다. 장파장 균일 strain에는 안정해도 유한 q의 phonon 모드가 불안정할 수 있습니다. 탄성과 phonon은 서로 보완하는 검사입니다.
04이 글의 용어 사전
본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.
- strain tensor
- 길이와 각도의 작은 상대 변화를 나타내는 대칭 2차 텐서입니다.
- stress tensor
- 에너지를 strain으로 미분한 단위면적당 내부 힘의 텐서입니다.
- 탄성텐서elastic tensor
- 선형 영역에서 stress와 strain을 잇는 4차 텐서입니다.
- clamped-ion 탄성
- strain을 가한 뒤 내부 원자좌표를 고정한 고주파적 탄성 응답입니다.
- relaxed-ion 탄성
- 각 strain에서 내부 좌표를 완화해 얻는 정적 탄성 응답입니다.
- Born 안정조건
- 작은 균일 strain에 대한 에너지 2차 형식이 양수가 되도록 하는 탄성상수 조건입니다.
05참고자료와 더 읽을거리
- Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
- Blöchl, Projector augmented-wave method (1994) — 평면파 효율성과 핵 주변의 전전자 정보를 잇는 PAW 방법의 원 논문입니다.
- Perdew, Burke & Ernzerhof, Generalized Gradient Approximation Made Simple (1996) — 재료 계산에서 널리 쓰이는 PBE-GGA 교환-상관 범함수의 원 논문입니다.