Dear all
I am trying to calculate phase transition of CsSnBr3 from pm3m 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 xaxis and enthalpy at yaxis, 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 xaxis and total energy at yaxis, 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
Phase transition under pressure?
Phase transition under pressure?
 Attachments

 CsSnBrcubic.in
 (1.03 KiB) Downloaded 10 times

 CsSnBr3tet.in
 (920 Bytes) Downloaded 9 times
Re: Phase transition under pressure?
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
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
Re: Phase transition under pressure?
Dear Eric
When i was varying lattice parameters of the noncubic phase, i was varying in somehow to keep a/c ratio constant..
I have tried optimization at hydrostatic pressures also, The noncubic phase was turning to almost the cubic one, means it gives me nearly the same results...
Regards
Nawzad
When i was varying lattice parameters of the noncubic phase, i was varying in somehow to keep a/c ratio constant..
I have tried optimization at hydrostatic pressures also, The noncubic phase was turning to almost the cubic one, means it gives me nearly the same results...
Regards
Nawzad
Re: Phase transition under pressure?
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
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.
Re: Phase transition under pressure?
Dear Eric
check that
Thanks for your help always
Best regards
Nawzad
check that
Thanks for your help always
Best regards
Nawzad