Phase transition under pressure?

Total energy, geometry optimization, DFT+U, spin....

Moderator: bguster

Locked
User avatar
new_986
Posts: 77
Joined: Wed Aug 04, 2010 9:27 pm

Phase transition under pressure?

Post by new_986 » Mon May 25, 2020 12:12 pm

Dear all
I am trying to calculate phase transition of CsSnBr3 from pm-3m cubic to p4/mbm tetragonal under pressure
the cubic contains 5 atoms per unit cell while tetragonal contains 10 atoms per unit cell
I have varied lattice parameters of the unit cell of both and calculated energy then I have calculated enthalpy (E+pv)
When I plot between P at x-axis and enthalpy at y-axis, the graph of both phases should cross each other at a certain pressure (transition pressure) BUT they don't
When I plot between volume at x-axis and total energy at y-axis, the graph of both CROSS each other at a certain pressure
I have some questions please:
1. Can anybody explain to me what is happening?
2. I have considered z value of tetragonal = 2, because it contains 10 atoms twice more than number of atoms of tetragonal (10). I want tot know if this is correct??
3. I have used eV,GPa and A^3 units for energy, pressure and unit cell volume to calculate enthalpy, Is this correct?

Thanks in advance

Nawzad
Attachments
CsSnBrcubic.in
(1.03 KiB) Downloaded 194 times
CsSnBr3tet.in
(920 Bytes) Downloaded 174 times
Capture.JPG
Capture1.JPG

ebousquet
Posts: 469
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: Phase transition under pressure?

Post by ebousquet » Thu Jun 04, 2020 4:03 pm

Dear Nawzad,
If you constrain the cell with acell, then you compute the internal energy at constrained strain and not at hydrostatic pressure. This is important mostly when the system is not cubic since the constrain along a,b and c will not represent the same pressure. I guess this can be your problem.
Best wishes,
Eric

User avatar
new_986
Posts: 77
Joined: Wed Aug 04, 2010 9:27 pm

Re: Phase transition under pressure?

Post by new_986 » Sun Jun 07, 2020 11:51 am

Dear Eric
When i was varying lattice parameters of the non-cubic phase, i was varying in somehow to keep a/c ratio constant..
I have tried optimization at hydrostatic pressures also, The non-cubic phase was turning to almost the cubic one, means it gives me nearly the same results...

Regards

Nawzad

ebousquet
Posts: 469
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: Phase transition under pressure?

Post by ebousquet » Mon Jun 08, 2020 9:27 am

Dear Nawzad,
In your case, I would check the phonon dispersions of the cubic 5 atom cell at different pressure to check whether you have an instability or not. If no instability is present, it will mean that the cubic phase is indeed stable in your calculation. If an instability appears (reported as negative frequencies in Abinit) then the eigenvector of the unstable phonon will give you the distortions that allow to lower the energy of the crystal.
Best wishes,
Eric
Last edited by ebousquet on Sat Jun 20, 2020 9:22 am, edited 2 times in total.

User avatar
new_986
Posts: 77
Joined: Wed Aug 04, 2010 9:27 pm

Re: Phase transition under pressure?

Post by new_986 » Thu Jun 11, 2020 12:13 pm

Dear Eric

check that
Thanks for your help always

Best regards

Nawzad

Locked