[QE-users] Pressure/Stress in vc-relax calculation
Inna Martha Nangoi
inna_nangoi at yahoo.com.br
Fri Nov 8 14:51:03 CET 2019
Dear all,
I'm proceeding with a vc-relax calculation with conv_thr = 10-6 ; energy and force default values. The pressure goes to zero but at the last structure the pressure increases to ~ -8 kbar. Since the pressure is negative one must think that the cell is expanding, but what it is observed (checking the cell parameters) is that the cell is contracting.Why is this happening?
Thank you in advance,
Dr. Inna NangoiGrupo de Físico-Química de Sólidos e Interfaces (GFQSI)http://www.ufjf.br/gfqsi/Universidade Federal de Juiz de Fora - UFJF
INPUT: / &SYSTEM ibrav = 0, celldm(1) = 1.8897265, nat = 408, ntyp = 2, ecutwfc = 40 , ecutrho = 400 , / &ELECTRONS conv_thr = 1.D-6 , / &IONS / &CELL /CELL_PARAMETERS cubic 16.467196779 0.003155548 0.004141842 0.003137055 16.453425288 -0.000761157 0.004171643 -0.000775474 16.446571996ATOMIC_SPECIES O 16.00000 O.pbe-n-kjpaw_psl.1.0.0.UPF H 1.01000 H.pbe-kjpaw_psl.1.0.0.UPFATOMIC_POSITIONS angstromO 9.321137481 14.266560377 9.327818241H 9.682105309 14.676153052 8.481823531H 10.049781829 14.333756337 10.015860139O 11.166283791 6.044399971 2.975926238.. . K_POINTS gamma
OUTPUT:
End of BFGS Geometry Optimization
Final enthalpy = -5995.1370960771 RyBegin final coordinates new unit-cell volume = 28794.22023 a.u.^3 ( 4266.86316 Ang^3 ) density = 0.95375 g/cm^3
CELL_PARAMETERS (alat= 1.88972650) 16.229014053 0.002449560 0.005745845 0.002429308 16.219588894 0.000479082 0.005772101 0.000461059 16.209758644
ATOMIC_POSITIONS (angstrom)O 9.180431407 14.143089383 9.191719262H 9.561522937 14.567241497 8.352842759H 9.908913197 14.170313644 9.891402550O 11.026600724 5.996092777 2.933422632H 11.877618746 6.446014321 2.622523705H 10.719968316 6.437111173 3.792164500O 5.141235135 15.134646177 1.987199104H 4.933636936 15.342422258 1.017207703H 5.853001846 14.414601411 1.987420436 . . . H 13.338033020 5.155154854 13.251553812O 11.079195321 9.183121298 6.146092164H 10.379755422 9.912994593 6.143398011H 11.912644022 9.559263869 5.711046728End final coordinates .
. . Total force = 0.008044 Total SCF correction = 0.000070
Computing stress (Cartesian axis) and pressure
negative rho (up, down): 9.616E-06 0.000E+00 total stress (Ry/bohr**3) (kbar) P= -8.32 -0.00005717 -0.00000010 -0.00000025 -8.41 -0.01 -0.04 -0.00000010 -0.00005654 0.00000006 -0.01 -8.32 0.01 -0.00000025 0.00000006 -0.00005604 -0.04 0.01 -8.24
total stress (Ry/bohr**3) (kbar) P= -12.00 total stress (Ry/bohr**3) (kbar) P= -10.74 total stress (Ry/bohr**3) (kbar) P= -7.11 total stress (Ry/bohr**3) (kbar) P= -2.92 total stress (Ry/bohr**3) (kbar) P= 0.65 total stress (Ry/bohr**3) (kbar) P= 0.70 total stress (Ry/bohr**3) (kbar) P= 0.30 total stress (Ry/bohr**3) (kbar) P= 0.08 total stress (Ry/bohr**3) (kbar) P= -0.19 total stress (Ry/bohr**3) (kbar) P= -0.43 total stress (Ry/bohr**3) (kbar) P= -0.58 total stress (Ry/bohr**3) (kbar) P= -0.40 total stress (Ry/bohr**3) (kbar) P= -0.06 total stress (Ry/bohr**3) (kbar) P= 0.25 total stress (Ry/bohr**3) (kbar) P= 0.42 total stress (Ry/bohr**3) (kbar) P= 0.43 total stress (Ry/bohr**3) (kbar) P= 0.26 total stress (Ry/bohr**3) (kbar) P= -0.03 total stress (Ry/bohr**3) (kbar) P= -0.29 total stress (Ry/bohr**3) (kbar) P= -0.53 total stress (Ry/bohr**3) (kbar) P= -0.37 total stress (Ry/bohr**3) (kbar) P= -0.19 total stress (Ry/bohr**3) (kbar) P= 0.18 total stress (Ry/bohr**3) (kbar) P= 0.36 total stress (Ry/bohr**3) (kbar) P= 0.50 total stress (Ry/bohr**3) (kbar) P= 0.33 total stress (Ry/bohr**3) (kbar) P= 0.06 total stress (Ry/bohr**3) (kbar) P= -0.09 total stress (Ry/bohr**3) (kbar) P= -0.34 total stress (Ry/bohr**3) (kbar) P= -0.28 total stress (Ry/bohr**3) (kbar) P= -0.10 total stress (Ry/bohr**3) (kbar) P= 0.03 total stress (Ry/bohr**3) (kbar) P= 0.25 total stress (Ry/bohr**3) (kbar) P= 0.25 total stress (Ry/bohr**3) (kbar) P= -8.32
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20191108/8278ee49/attachment.html>
More information about the users
mailing list