﻿

## Phase transition under pressure?

structure determination,...

Moderators: jzwanzig, jolafc

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

### Phase transition under pressure?

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
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?

Attachments
CsSnBrcubic.in
CsSnBr3tet.in

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

### Re: Phase transition under pressure?

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

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

### Re: Phase transition under pressure?

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

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

### Re: Phase transition under pressure?

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.

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

Dear Eric

check that