Dear all,
Hi! I hope you are fine. I have a question?!
Why do we have residual stress when simulation box is relaxed under NVT?
Thanks
Nili
The simulation box can’t be relaxed if you only did NVT runs. NVT forces the simulation box to have constant volume which whatever you’re simulating may not like. If you use NPT with the aniso keyword, you should have negligible residual stresses (and thus a relaxed state).
Will
Dear professor Pisani,
Thanks for your help. But, when I use NPT in this relaxation, my Strain Energy Density in Compression (not in tension) isn’t minimum at first step which is weird.
As my result shows, the minimum value appears in third step. How can I fix it?
-24654.8 -24663.4 -24669.5 -24663.7 -24666.5 -24661 -24646.5 -24632.4 -24614.6 -24591.4 -24554 -24518.9 -24477.4 -24421.4 -24366.4 -24304.7 -24241.6 -24167.6 -24092.6 -24009.1 -23927.2 -23840.7 -23748.2 -23766.6 -23895 -23890.3 -23848.2 -23783.2 -23729.6 -23653.8 -23561.1
Best Regards,
Nili
You don’t need to fix it. The difference is only -0.06% from the initial value. That’s small enough to not worry about it.
If this data is from a pure NPT relaxation (no loads), then your average stresses (and the elastic strain energy) will trend toward zero (which is what you want). Negligible stresses are near-zero stresses. Eventually your stresses will look like typical MD pressure fluctuations.
Hope this helps,
Will