QuantumATK Forum
QuantumATK => General Questions and Answers => Topic started by: fangbiaolee on January 12, 2023, 11:51
-
(https://pic1.zhimg.com/80/v2-8b8f6a164105c094a75d5c24d5254858_r.jpg)
Phase diagram of KPF6 (potassium hexafluorophosphate)
(https://pic1.zhimg.com/80/v2-e7575f51346336684cc9bae3e8a78b16_r.jpg) (https://pic1.zhimg.com/80/v2-d0fb7eeae85f7b4ef669a6ce7f90541e_r.jpg) KPF6 crystal
(https://picx.zhimg.com/80/v2-6506e255c18b79423325045eb8b48e4b_r.jpg) (https://picx.zhimg.com/80/v2-ff638614a99f46edbf02a8ffbded1a7d_r.jpg) The script
(https://pic4.zhimg.com/80/v2-433589d7cb799467d38216f34f1f338d_r.jpg)
The trajectory of KPF6
Simulation conditions:
1000 atoms(5x5x5);
time_step=1*femtoSecond,
reservoir_temperature=350*Kelvin,
reservoir_pressure=[1000, 1000, 1000, 1000, 1000, 1000]*bar,
thermostat_timescale=100*femtoSecond,
barostat_timescale=1000*femtoSecond,
initial_velocity=initial_velocity,
heating_rate=0*Kelvin/picoSecond,
compression_rate=0*bar/femtoSecond,
coupling_mask=[[True, True, True], [True, True, True], [True, True, True]],
chain_length=3,
1. Why does the simulation box eventually become so distorted? According to the phase diagram, KPF6 crystal may undergo phase transformation. Is it necessary to define the simulation box as a triclinic system?
2. The chemical bond of PF6 has broken, but it is obvious that the chemical bond of PF6 will not break at this temperature.
3. Are these results due to the inaccuracy of LCAO? However, the AIMD based on LCAO should predict the physical laws, even though the value may have some deviation.
I have been troubled by this problem for a long time and can't solve it with all kinds of methods. If someone can help me solve it, I will thank him very much.
-
I suppose a lot of unexpected things happen at 1,000x atmospheric pressure. Maybe first make sure the calculation is stable at standard pressure.