VASP-Maschinenlernkraftfeld Berechnet Silizium-Phononenspektrum
Vienna Ab initio Simulation Package (VASP:Das VASP-Handbuch – VASP Wiki ) ist ein Computerprogramm zur Materialmodellierung auf atomarer Ebene von Grund auf, wie etwa Berechnungen der elektronischen Struktur und quantenmechanische Molekulardynamik.
PhonopieWillkommen bei Phonopy — Phonopy v.2.37.1) ist ein Python-Toolkit zum Berechnen von Phononenbandstrukturen, thermischen Eigenschaften, Gruppengeschwindigkeiten und anderen Phononen-bezogenen Größen auf harmonischer und quasiharmonischer Ebene.
Tutorial-Inhalte
In diesem Tutorial verwenden wir ein automatisiertes Skript, um den Berechnungsprozess des Kraftfeld-Phononenspektrums des maschinellen Lernens zu demonstrieren. In diesem Tutorial lernen Sie den grundlegenden Prozess der Berechnung des Phononenspektrums mithilfe von Kraftfeldern des maschinellen Lernens kennen:
- Vorbereiten der Störungsdifferenz-Superzellenstruktur
- Legen Sie die Kraftfeldparameter für maschinelles Lernen fest und berechnen Sie die Gesamtenergie aller Strukturen
- Berechnen Sie die Kraftkonstantenmatrix
- Berechnen Sie das Phononenspektrum aus der Kraftkonstantenmatrix
Eingabedatei
Die Eingabedatei enthält
├── POSCAR-unitcell
├── clean.sh
├── run.sh
└── run_vasp.sh
├── pt
│ ├── INCAR
│ ├── KPOINTS
│ ├── ML_FF
│ └── POTCAR(需自行上传)
POSCAR-Einheitszelle
Si #硅结构
5.38930000000000
0.0000000000000000 0.5071343999939496 0.5071343999939496
0.5071343999939496 0.0000000000000000 0.5071343999939496
0.5071343999939496 0.5071343999939496 0.0000000000000000
2
Direct
0.8750000000000000 0.8750000000000000 0.8750000000000000
0.1250000000000000 0.1250000000000000 0.1250000000000000
sauber.sh
#!/bin/bash
rm -r *.yaml band.pdf band.yaml FORCE_SETS vasp poscar *out SPOSCAR
# 删除不必要文件
ausführen.sh
#!/bin/bash
rm -r vasp poscar
#准备微扰差分超胞结构
##########################
phonopy -d --dim 2 2 2 --pa auto -c POSCAR-unitcell
mkdir poscar
mv POSCAR-unitcell pp
mv POSCAR-* poscar/
##########################
#计算所有结构的总能
##########################
mkdir vasp
cd vasp
Pnum=$(ls -l ../poscar/ -IR | grep "^-" | wc | awk -F ' ' '{print $1}')
cp ../run_vasp.sh .
t_head="for i in {1.."
t_tail="}"
sed -i "3c ${t_head}${Pnum}${t_tail}" run_vasp.sh #生成vasp计算脚本
./run_vasp.sh
##########################
cd ../
mv pp POSCAR-unitcell
#计算力常数矩阵
##########################
phonopy -f vasp/*/vasprun.xml > pfcout
##########################
#根据力常数矩阵计算声子谱
##########################
phonopy-load --band "0.0 0.0 0.0 0.5 0.0 0.5 0.625 0.25 0.625, 0.375 0.375 0.75 0.0 0.0 0.0 0.5 0.5 0.5" --band-labels "Γ X K K Γ L" -p -s
##########################
run_vasp.sh
#!/bin/bash
#计算 vasp 流程自动脚本
for i in {1..8}
do
rm -r $i
mkdir $i
cd $i
ii=$(printf "%03d" $i)
cp ../../poscar/POSCAR-${ii} POSCAR
cp ../../pt/* .
mpirun -n 1 vasp_std
cd ../
done
Dateien im Ordner pt
INCAR-Standardstatikberechnung
ISTART = 1 (若有波函数、读取波函数)
ISPIN = 1 (非极化计算)
Static Calculation
ISMEAR = 0 (高斯占据)
SIGMA = 0.05 (高斯展宽)
NELM = 60 (最大电子步)
EDIFF = 1E-08 (SCF 收敛精度)
Einstellungen für das Kraftfeld für maschinelles Lernen
ML_LMLFF = T (使用机器学习力场)
ML_ISTART = 2 (仅预测模式)
KPOINTS (wird nicht verwendet, wenn maschinelles Lernen aktiviert ist, legen Sie einfach einen fest)
K-Spacing Value to Generate K-Mesh: 0.040
0
Gamma
4 4 4
0.0 0.0 0.0
POTCAR
Die Pseudopotentialkombination der entsprechenden Elemente des Systems, hier das Pseudopotential von Si
ML_FF
Die Kraftfeldparameterdatei für maschinelles Lernen, die dem vortrainierten System entspricht. Den genauen Trainingsablauf finden Sie unter Si's Kraftfeld für maschinelles LernenLernprogramm.
Erste Schritte
1. Klonen Sie den Container
Suchen Sie das Arbeitsverzeichnis des Tutorials und klonen Sie den Container
2. Container aufstellen
选择 4090 — 按量付费 — vasp 6.3.0-cuda11.8 — 工作空间

加载完毕后打开工作空间

打开终端

上传准备好的硅赝势,可以使用官网例子里的赝势 POTCAR 放到目录中

3. Installieren Sie die Phonopie-Umgebung
conda install -c conda-forge phonopy
Geben Sie dann y ein und drücken Sie die Eingabetaste, um der Installation zuzustimmen

4. Führen Sie das Skript aus
chmod 777 *.sh
./run.sh
5. Betrachten Sie das Phononenspektrum
Sehen Sie sich die Datei band.pdf an
