## Discovery of an isostructural phase transition within orthorhombic phase field of CaTiO3

Saurabh Tripathi1, Anil Kumar2, Sharmila Shirodkar3, Masatomo Yashima4, U. V. Waghmare2, Dhananjai Pandey1

1 School of Materials Science and Technology, Indian Institute of Technology (Banaras Hindu University), Varanasi-221005, India

2 Theoretical  Sciences  Unit,  Jawaharlal  Nehru  Centre  for  Advanced  Scientific  Research,  Bangalore 560064, India

3 Department  of  Chemistry  and  Materials  Science,  Graduate  School  of  Science  and  Engineering,  Tokyo Institute of Technology, 2-12-1-W4-17, O-okayama, Meguro-ku, Tokyo, 152-8551, Japan

Earth's lower mantle extending from 670 to 2,990 km deep is predominantly composed of a perovskite-type (Mg,Fe)SiO 3 phase 1,2 .    The perovskite phase undergoes a structural phase transition to a post-perovskite phase responsible for D'' layer seismic discontinuity 2,3 at  about  2690  km  depth  in  the  lowermost  region  of  the  lower  mantle. However, structural basis of other seismic discontinuities occurring in the upper region of the  lower  mantle  (700  km  to  1,200  km  deep)  remains  unexplained 4-7 ,  as  no  apparent change in the crystal symmetry of the orthorhombic perovskite phase has been reported . 5 We present here unambiguous evidence for a non-apparent isostructural phase transition 8 in  the  stable  orthorhombic  perovskite  phase  of  CaTiO   which  may  have  relevance  to 3 phase  transitions  in  the  perovskite  phase  of  (Mg,Fe)SiO   also,  as  both  the  compounds 3 have similar structure, tolerance factor and thermochemical properties 9-11 . Our results are based  on  the  analysis  of  neutron  powder  diffraction  patterns  using  Rietveld  and  mode crystallography  techniques  and  are  supported  by  density  functional  and  Landau  theory calculations.  The  present  results  on  CaTiO   would  encourage  search  for  isostructural 3 phase  transition  in  the  perovskite  phase  of  (Mg,Fe)SiO   that  may  provide  clue  to  the 3 unexplained geophysical phenomena in the upper part of the earth's lower mantle.

The  family  of  oxide  perovskites  with  the  general  formula AB O3 is  named  after  the mineral name of CaTiO3 9-10 where A and B are larger and smaller cations, respectively. CaTiO3, whose structure and the tolerance factor are similar to the perovskite phase of the earth's lower mantle material (Mg,Fe)SiO 3 10-11 , undergoes structural phase transitions at easily achievable temperatures 12 , as compared to (Mg,Fe)SiO3 that requires extremely high temperatures and pressures . CaTiO 2 3 has therefore been used as a model system to mimic the phase transition behaviour of (Mg,Fe)SiO .  The stable high temperature cubic 3 phase (space  group Pm3m)  of CaTiO3 transforms to a tetragonal  phase (space group I 4/ mcm ) at T ~1636 K  which  can  be visualized  to occur through condensation of one of the components of the R 4+ mode at q = ‰ ‰ ‰ of the cubic Brillouin zone culminating in anti-phase rotated  oxygen octahedra in the a a c 0 0 -tilt system 12-13 . This tetragonal phase on further cooling below T = 1512 K  transforms to an orthorhombic phase in the Pnma space group  with  a b a  tilt system -+ -14-15 due to the condensation of  M 3+ mode 12-13 at q= ‰ ‰ 0  involving  in-phase rotation and  another component of the R 4+ mode causing additional  anti-phase  rotation  of  the  octahedra.  This  orthorhombic  phase  is  widely accepted  as  the  stable  phase  of  CaTiO 3 12-13 under  ambient  conditions  and  of  MgSiO 3 under  the  lower  mantle  conditions .        From  an  analysis  of  the  temperature  dependent 4 neutron diffraction patterns of CaTiO  within the stability field of the orthorhombic Pnma 3 phase, we  present evidence for discontinuous changes in the unit cell volume, pseudomonoclinic cell distortion angle and the in-phase as well as anti-phase tilt angles in between 1000-1069K indicative of a first order isostructural phase transition that does not involve any change of space group and/or Wyckoff site symmetries. We show that

this  transition  is  phonon  driven  as  the  amplitudes  of  the  R 4+ and  M 3+ phonon  modes, obtained from amplimode analysis of neutron powder diffraction data and first principles calculations, exhibit anomalous non-analytic behavior at the transition temperature. We also  present  results  of  Landau  theory  to  show  how  a  synergy  between  the  two  order parameters emerging at higher temperature phase transitions involving in-phase and antiphase tilts of oxygen octahedra can lead to such an isostructural phase transition.

In  CaTiO3 12 , as  also  in  MgSiO 3 16 , the  hhh  reflections  are  doublet  and  hence  the equivalent  elementary  perovskite  cell  can  be  taken  to  be  pseudomonoclinic  (for  more details, see the structure refinement section of the supplementary file). The pseudomonoclinic cell parameters of CaTiO3 at room temperature show the relationship 17 : a pm =cpm&lt;bpm, α = =90 γ 0 ≠β . The temperature variation of the elementary perovskite cell parameters  a pm =  c pm ,  b pm and  the  monoclinic  angle β pm of  CaTiO ,    as  obtained  from 3 Rietveld refinements,  shown in Fig 1 reveals that the relationship a pm =cpm&lt;bpm  changes to  a pm =cpm&gt;bpm somewhere between 1000 to 1069K. Further, the monoclinic angle γ pm changes discontinuously in between 1000-1069K (see inset 'a' in Fig. 1). The unit cell volume also shows a discontinuous change in the same temperature range (see inset 'b' in Fig. 1). These discontinuous changes in the unit cell parameters occur  within the stability field of the orthorhombic phase in the Pnma space group, below the cubic to tetragonal and tetragonal to orthorhombic phase transitions at 1636 and 1512K, respectively, and therefore  represent  an  isostructural  phase  transition  (ISPT),  which  is  quite  rare  in nature 9,18,19 .

In the a b -+ a  tilt system of the Pnma phase of CaTiO , the tilt angles for the anti-phase -3 rotations (a ) of the oxygen octahedra about the pseudocubic  [100]  and [001]  axes are -p p equal and are equivalent to a single tilt through an angle φ about one of the two diad axes parallel to the [101]  and [101]  directions, where the subscript p refers to the ideal high p p symmetry elementary perovskite cell.  The  in-phase  octahedral  rotation  (b )  through  an + angle ζ is about the [010]  axis. Both the tilt angles, p φ and ζ , can be determined from the coordinates of  the apical oxygen O , which changes from 0.25, 0, 0.75 in the cubic phase I to 0.25+u, v, 0.75-w in the orthorhombic structure, using the relationships given in Ref 2021 :    It  is  evident  from Fig 2 (a, b), which depicts the temperature variation of  the antiphase  ( φ )  and  in-phase  ( ζ )  tilt  angles,  that  both  the  angles  undergo  a  discontinuous change  in between 1000-1069K similar to the  unit cell parameters shown in Fig 1. Since the tilt angles are the macroscopic order parameters for the AFD phase transition, their discontinuous change in between 1000-1069K in Fig 2 (a, b) confirms the existence of a first order AFD transition of isostructural type. To understand the likely role of phonons in the AFD isostructural transition, we also determined the amplitudes of the R 4+ and M3+ modes of the cubic perovskite structure 22-23 from the refined positional coordinates using the Amplimodes suite 24-25 .  We find that the amplitudes of  R 4+ and  M 3+ modes change discontinuously at the isostructural phase transition (see Fig 2 (a, b). This suggests the role of phonons in this transition, as confirmed further using first principles calculations.

For a system with more than one unstable mode, there could be competition as well as coupling  between  the  modes  to  drive  the  system  to  the  most  stable  structure.  The amplitude  of  the  unstable  mode,  which  is  a  measure  of  deformation  of  the  high

symmetric  high  temperature  structure  that  yields  the  ground  state  structure,  provides important  information  about  the  structural  phase  transitions.  For  a  given  structure,  the amplitude of unstable mode shows a continuous behaviour as a function of temperature or  external  field,  but  shows  a  non-analytic  behaviour  close  to  a  phase  transition.  The unstable modes at M and R points in the Brillouin zone of cubic CaTiO3 'freeze-in' or distort the structure  leading  to  the  orthorhombic  phase  at  low  temperatures.  Our calculations of the cubic structure show that the M and R modes are unstable in CaTiO , 3 in  agreement  with  the  findings  of  previous  workers 26-27 ,  and  their  frequency  decreases monotonically with temperature as shown in Fig.3 (a).  To understand the behaviour of M and R mode amplitudes and associated phase transitions with temperature, as observed in our  experiments,  we  allow  the  system  to  relax  to  minimum  energy  by  freezing  the unstable zone boundary phonon modes using a 2 x 2 x 2 supercell.  The energy profile as a function of atomic displacements along the soft eigen modes provides information about the mode amplitude (of atomic displacements) in the minimum energy structure. We determined energy profile as a function of atomic displacements for the experimental structural  parameters  at  different  T  (see  Fig.  2  (a,  b)  of  Extended  data  figure  file).  We find that energy gain by freezing of the M 3+ and R4+ modes decreases with the increase of temperature (a similar trend can be seen in the increase of phonon frequencies for these modes  with  temperature  in  Fig  3(a)).  The  amplitudes  of  M 3+ and  R 4+ modes  decrease continuously  with  increasing  temperature  but  show  a  non-analytical  behaviour  with  a discontinuous  change  in  between  1000  to  1069K  (see  Fig.  3(a)),  consistent  with  our experimental observations shown in Fig. 2(a, b). Above the isostructural phase transition

temperature, the temperature dependence of the two amplitudes are continuous within the stability limit of the Pnma phase.

In  complement  to  first-principles  demonstration  of  the  discontinuity  in  R 4+ and  M 3+ phonon modes across the isostructural transition occurring somewhere in between 1000 to 1069K (as a result of observed discontinuity in volume), we now use phenomenological  Landau  theory  and  show  that  this  isostructural  phase  transition originates  from  the  synergistic  coupling  between  order  parameters  emerging  at  higher temperature  transitions.  From  the  Landau  free  energy  of  higher  temperature  phase transitions (see Supplementary file for complete derivation), we obtain free energy as a function of x and y, the two order parameters related to the R 4+ mode (along (0,a,a)) and M3+ mode (along (a,0,0)), respectively 28 .

$$& F = e _ { 1 } x + \frac { a _ { 1 } } { 2 } ( \frac { T } { T _ { \text{C} } } + p ) x ^ { 2 } + \frac { b _ { 1 } } { 3 } x ^ { 3 } + \frac { c _ { 1 } } { 4 } x ^ { 4 } \\ & + e _ { 2 } y + \frac { a _ { 2 } } { 2 } ( \frac { T } { T _ { \text{H} } } + p ) y ^ { 2 } + \frac { b _ { 2 } } { 3 } y ^ { 3 } + \frac { c _ { 2 } } { 4 } y ^ { 4 } \\ & + d _ { 1 } x ^ { 2 } y + d _ { 2 } x y ^ { 2 } + d _ { 3 } x y + d _ { 4 } x ^ { 2 } y ^ { 2 },$$

where {e1, a , b , c } are the Landau coefficients for the order parameter (x) associated 1 1 1 with the R4+ mode (along (0,a,a)), and {e , a , b , c } are the Landau coefficients of the 2 2 2 2 order parameter (y) associated with M 3+ mode [along (a,0,0)], and d 's (i=1 to 4) are the i coupling constants. Renormalization of the coefficients of the quadratic and linear terms (refer  to  Supplementary  file)  gives  the  effective  transition  temperatures  of  T C  and  T H. We  fit  them  to  1000  K  to  reproduce  the  experimentally  observed  isostructural  phase transition (refer to Supplementary file). The values of parameters were fixed at e = 0.3, 1

a1 =1, b = 2, c = 1, T 1 1 C= 1000 K and p= 0, and e2= 0.1, a2 = 1, b = 2, c = 1, T 2 2 H= 1000 K and  q=  0,  such  that  the  value  of  free  energy  functions  of  the  decoupled x and y order parameters is 0 at T C and TH, respectively. The coupling coefficients d , d , d , d  are 0.3, 1 2 3 4 0.1,  0.4  and  0.1  respectively.  Equilibrium  values  of x and y are  obtained  through minimization of free energy at each temperature 29 and are depicted in Fig. 3 (b), which shows  a  clear  discontinuity  in  both  the x and y order  parameters  near  1000K,  in agreement  with  experimental  observations.  Thus,  the  isostructural  phase  transition  is  a result of nonlinear coupling between two distinct three-component order parameters that have  already  broken  the  symmetry,  and  the  richness  in  the  symmetry  equivalent structures possible.

We  believe  that  the  discovery  of  the  isostructural  phase  transition  in  CaTiO   will 3 encourage  search  for  a  similar  transition  within  the  stability  field  of  the  orthorhombic perovskite phase of MgSiO3/(Mg,Fe)SiO3 that may hold clue to the hitherto unexplained geophysical  phenomena   in  the    upper  part  of  the  earth's  lower  mantle  consisting  of 5 perovskite phase of  (Mg,Fe)SiO  predominantly. 3

Methods: The details of sample preparation, neutron data collection, Rietveld refinements and Amplimode analysis, first principles calculations and Landau theory are given in the supplementary files.

## References

- 1.    Peyronneau,  J.  &amp;  Poirer,  J.  P.  Electrical  conductivity  of  the  Earth's  lower  mantle. Nature 342 , 537-539, 1989.
- 2. Oganov, A. R. &amp; Ono, S. Theoretical and experimental evidence for a post-perovskite phase of MgSiO3 in Earth's D" layer. Nature 430 , 445-448 (2004).

- 3.  Murakami,  M.,  Hirose,  K.,  Kawamura,  K.,  Sata,  N.  &amp;  Ohishi,  Y.  Post-Perovskite Phase Transition in MgSiO3. Science 304 , 855-858 (2004).
- 4.  Stixrude,  L.  &amp;  Cohen,  R.  E.  Stability  of  orthorhombic  MgSiO   perovskite  in  the 3 Earth's lower mantle. Nature 364 , 613-616 (1996).
- 5.  Revenaugh,  J.  &amp;  Jordan,  T.  H.  Mantle  layering  from  ScS  reverberations  2.  The transition zone. J. Geophysical Research. 96 , 763-779 (1991).
- 6. Lay, T. Structure of the core-mantle transition zone: a chemical and thermal boundary layer. EOS 70 , 49, 54-55, 58-59 (1989).
- 7. Kawakatsu, H. &amp; Niu, F. Seismic evidence for a 920-km discontinuity in the mantle. Nature 371 , 301-305 (1994).
- 8. Singh, A.,  Pandey, V., Kotnala, R. K. &amp; Pandey, D. Direct Evidence for Multiferroic Magnetoelectric  Coupling  in  0.9BiFeO -0.1BaTiO3 3 . Phys.  Rev.  Lett. 101 , 247602-1247601-4 (2008).
- 9. Galasso, F. S. Structure, properties, and preparation of perovskite-type compounds. Pergamon press 1969.
- 10. Salje, E. Characteristics of perovskite-related materials. Phil. Trans. R. Soc. Lond. A 328 , 409-416 (1989).
- 11. Gautason, B. &amp; Muehlenbachs, K. Oxygen diffusion in perovskite: implications for electrical conductivity in the lower mantle. Science 260 , 518-521 (1993).
- 12.  Ali, R. &amp; Yashima, M. J. Space group and crystal structure of the Perovskite CaTiO3 from 296K to 1720K. Solid State Chem. 178 , 2867-2872 (2005).
- 13.  Yashima,  M.  &amp;  Ali,  R.  Structural  phase  transition  and  octahedral  tilting  in  the calcium titanate perovskite CaTiO . Solid. State. Ionics. 3 180 , 120-126 (2009).

- 14.    Kennedy , B.  J., Howard,  C.  J.  and  Chakoumakos,  B.  C.  Phase  transitions  in perovskite  at  elevated  temperatures  -  a  powder  neutron  diffraction  study.  J.  Phys.: Condens. Matter 11 , 1479-1488 (1999).
- 15. Glazer, A.M. Simple Ways of Determining Perovskite Structures. Acta Crystallogr. Sec. A 31, 756-762 (1975).
- 16.  Meade,  C.,  Kwang  Mao,  H.  &amp;  Hu,  J.  High-Temperature  Phase  Transition  and Dissociation of (Mg,Fe)SiO3 Perovskite at Lower Mantle Pressures. Science 268, 17431745 (1995).
- 17. Ball , C. J., Begg, B. D., Cookson, D.J., Thorogood, G. J. &amp; Vance, E. R. Structures in the System CaTiO3/SrTiO3. J. Solid. State. Chem. 139 , 238-247, (1998).
- 18. Lee, S., Pirogov, A., Kang, M., Jang, K.-H., Yonemura, M., Kamiyama, T., Cheong, S.-W.,  Gozzo  F.,  Shin  N.,  Kimura  H.,  Noda  Y.,  &amp;  Park,  J,-G.  Giant  magneto-elastic coupling in multiferroic hexagonal manganites, Nature 451 , 805-808 (2008).
- 19. Bhattacharjee, S., Taji, K., Moriyoshi, C., Kuroiwa, Y., &amp; Pandey, D. Temperature induced isostructural phase  transition, associated large negative volume expansion, and the  existence  of  a  critical  point  in  the  phase  diagram  of  the  multiferroic  (1-x)BiFeO 3 xPbTiO3 solid solution system, Phys. Rev. B 84 , 104116 (1-11) 2011
- 20  Ranjan,  R.,  Pandey,  D.,  Siruguri,  V.,  Krishna,  P.  S.  R.  &amp;  Paranjape,  S.  K.  Novel structural features and phase transition behaviour of (Sr 1-x Cax)TiO3: I. Neutron diffraction study. J. Phys.: Condens. Matter 11 , 2233-2246 (1999).
- 21. Ranjan, R. Ph.D. Thesis of B.H.U., Varanasi, India.
- 22.  Darlington,  C.  N.  W.  Normal-mode  analysis  of  the  structures  of  perovskites  with tilted octahedra. Acta. Cryst. A 58 , 66-71 (2002).

- 23.  Darlington,  C.  N.  W.  Normal-mode  analysis  of  the  structures  of  perovskites  with tilted octahedra. Erratum. Acta. Cryst. A 58 , 299-300 (2002).
- 24. Aroyo, M. I., Kirov, A., Capillas, C., Perez-Mato, J. M. and Wondratschek. H. Bilbao Crystallographic  Server  II:  Representations  of  crystallographic  point  groups  and  space groups.  Acta Cryst. A 62 , 115-128 (2006).
- 25. Aroyo,  M.  I.,  Perez-Mato,  J.  M.,  Capillas,  C.  Kroumova,  E.,  Ivantchev,  S., Madariaga,  G.,  Kirov,  A.  &amp;  Wondratschek,  H.  Bilbao  Crystallographic  Server  I: Databases  and  crystallographic  computing  programs.  Zeitschrift  fuer  Kristallographie, 221 , 1, 15-27 (2006).
- 26. Cockayne, E. &amp; Burton, B. P. Phonons and static dielectric constant in CaTiO  from 3 first principles. Phys. Rev. B 62 , 3735-3743 (2000).
- 27.  Tripathi,  S.,  Kumar,  A.,  Waghmare,  U.  V.  &amp;  Pandey,  D.  Effect  of  NaNbO 3 substitution on the quantum paraelectric behavior of CaTiO . Phys. Rev. B 3 81 , 212101222108 (2010).

28. iso.byu.edu. Stokes, Harold.  T., Hatch,  Dorian.  M.,  &amp;  Campbell,  Branton.  J. ISOTROPY Software Suite.

- 29. Wolfram  Research, Inc., Mathematica, Version 7.0, Champaign, IL (2008). Registered  Name:  Waghmare,  Umesh.  V,  Organization:  Jawarharlal  Nehru  Centre  for Advanced Scientific Research, Licence Number: L3283-6940.

Supplementary Information is linked to the online version of the paper at www.nature.com/nature.

Author  Contributions: D. P. conceived  the  problem  and  applied  the  concept  of isostructural phase transition, supervised the analysis of the neutron diffraction patterns

as also the theoretical calculations and prepared the final version of the manuscript. ST analysed  the  neutron  diffraction  patterns  and  prepared  the  first  draft  of  the  paper.  MY provided the neutron diffraction data  collected on a sample prepared in  his laboratory. The details of sample preparation are given in Ref. No. 12. The data used in Ref. No. 12 was reanalyzed by ST leading to the discovery of the isostructural phase transition in the orthorhombic phase of CaTiO3. AK and UW carried out first principle calculations while SS and UW carried out the Landau theory calculations. AK, SS and UW prepared the first draft of the theoretical portion of the manuscript.

## Figure Legends

Fig. 1 Variation of equivalent elementary perovskite cell parameters a pm= c pm , b pm with temperature.  Inset shows variation of (a) monoclinic angle β pm and (b) volume of CaTiO3 with temperature. Error bars are displayed in the figures.

Fig. 2 Variation of (a) inphase tilt angle ( ζ ) and amplitude of M 3+ (A) primary mode, (b) antiphase tilt ( φ ) and amplitude of R 4+ (A) primary mode of Pnma phase of CaTiO3 with temperature. Error bars for tilt angles are displayed in the figures.

Fig. 3 (a) Variation of frequency ( ω ) and amplitude (A) of soft R and M primary modes for  cubic  structure  of  CaTiO   with  temperature,  (b)  'x'  and  'y'  order  parameters  with 3 temperature 29 . Note that since the values of 'x' and 'y' are negative, here we have plotted the absolute values.

Fig. 1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0335v1-figure-944cdb97-5192-4d0b-9420-6701051d031f.jpg)

Fig. 2

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0335v1-figure-43694b23-77a9-4e32-abe7-c04d86af97d9.jpg)

-1

c

Fig. 3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0335v1-figure-7bb0a00e-171b-4c00-a89c-021cc9dba0ba.jpg)

## Extented data

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0335v1-figure-dbe56747-6224-4f0f-aa0c-e57a27b146ab.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0335v1-figure-76d25ed7-2104-4b30-b173-02544a15d752.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0335v1-figure-081c2889-409b-496e-ab34-e665f55be5b3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0335v1-figure-c9b52f79-7920-451c-a275-ccc585583267.jpg)

Extended  data  Fig.  1  (a,  b,  c,  and  d)  depicts  the  observed,  calculated  and  difference profiles obtained after final cycle of refinement for T= 296K, 868K, 1069K and 1423K using Pnma space group.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0335v1-figure-bab5afe6-f08e-4281-afb1-7e364fa7cdfe.jpg)

Extended data Fig. 2 (a, b) Energy profile as a function of atomic displacements for the experimental structural parameters at different Temperature.

## Supplementary Information

## Structure refinement

The Rietveld refinements using the neutron powder diffraction patterns of CaTiO 3 were performed by FULLPROF  suite . The complete Fullprof suite and documentation 1 is available at the web . The peak shape was assumed to be pseudo-Voigt function with 2 asymmetry.  The  background  of  each  profile  was  approximated  by  a  6  coefficient polynomial. Except for the occupancy parameters of the atoms, which were kept at their nominal composition, all other parameters, i.e., scale factor, zero correction, background, half  width  parameters  along  with  mixing  parameters,  lattice  parameters,  positional coordinates, and thermal parameters, were refined.

Amplitudes  of  phonon  modes  were  calculated  using  AMPLIMODES  package  from Bilbao Crystallographic Server 3,4 .

The asymmetric unit of the orthorhombic phase of CaTiO3 in the Pnma space group consists of 4 atoms, viz Ca, Ti, O1 (apical oxygen) and O2 (planar oxygen) as

| Ca   |   1 | 4c   | 0.0+/-u   | 0.25    | 0+/-w    |
|------|-----|------|-----------|---------|----------|
| Ti   |   1 | 4b   | 0.0       | 0.0     | 0.5      |
| O    |   1 | 8d   | 0.25+/-u  | 0.0+/-v | 0.75+/-w |
| O    |   2 | 4c   | 0.50+/-u  | 0.75    | 0+/-w    |

There are 7  refineable  atomic  xyz  coordinates.  The  coordinates  of  all  the  atoms  in  the asymmetric unit at various temperatures from 296K to 1423K were obtained by Rietveld refinements. Fig. 1 (a, b, c, and d) of the Extended data figure file, depicts the observed, calculated and difference profiles obtained after final cycle of refinement for T= 296K,

868K, 1069K and 1423K using Pnma space group. Evidently, all the fits are very good as can be seen from the near flatness of the difference profile.

In  order  to  explain  the  isostructural  phase  transition  we  have  expressed  orthorhombic structure  (supercell)  of  CaTiO   with  Pnma  space  group  having  a b a   tilt  system  into 3 -+ -pseudomonoclinic elementary perovskite cell (subcell).  The relationships  between the orthorhombic cell parameters A  , B 0 0 and C0 and the elementary pseudocubic  perovskite cell parameters a pm, b pm , and c pm  are : (a pm)  = (c 2 pm)   = (A /2)  + (C /2) ,  A 2 0 2 0 2 0 2 =  a pm 2 + cpm 2 -2a pmcpmcos β , and C 0 2 =  a pm 2 + c pm 2 +2apmcpmcos β , and B0=2bpm, where β pm is the  angle  between  the  a pm  and  c pm  axes.  The  other  two  angles  are α pm= γ pm=90 0 ≠ β pm making the elementary cell pseudomonoclinic. For  A =C , the elementary perovskite cell 0 0 axes are orthogonal to each other. The hhh type reflections are doublet and singlet for the monoclinic  and  orthogonal  cells,  respectively.  In  CaTiO ,  as  also  in  MgSiO ,  the  hhh 3 3 reflections are doublet and hence the equivalent elementary perovskite cell can be taken to be pseudomonoclinic.

## Mode Crystallography:

Mode  crystallographic  analysis  was  done  using  AMPLIMODE  suite  at  Bilbao Crystallographic  Server 3,4 , using  Rietveld  refined  coordinates.  AMPLIMODES  is  a computer program which calculates the  amplitude of  symmetry adapted modes for the low symmetry distorted structure with respect to a high symmetry structure. These modes are given in relative units using atomic displacements of the atoms in the asymmetric unit of the distorted phase.

In order to work with the AMPLIMODE analysis, we need to know the parent [or high symmetry (HS)] structure, the [low symmetry structure (LS)] and the transformation matrix for the two phases, and the symmetry adapted primary and secondary modes that transform  the  HS  phase  into  the  LS  phase.  In  the  case  of  CaTiO ,  the  high  symmetry 3 structure is cubic with space group Pm 3 m , where A site cation, B site cation and  O go to 1b,  1a and 3d  Wyckoff sites, with   fractional coordinates  0.5 0.5 0.5, 0.0 0.0 0.0 and 0.5 0.0 0.0, respectively.

The  transformation  matrices  relating  the  unit  cell  of  the  low  symmetry  room temperature phase in the Pnma space group with the unit cell of the high symmetry cubic phase of CaTiO3 is given below:

$$[ \begin{matrix} 1 & 0 & 1 ] \, [ & 1 / 2 ] \\ [ & 0 & 2 & 0 ] \, [ & 0 ] \\ [ & - 1 & 0 & 1 ] \, [ & 1 / 2 ] \end{matrix} ]$$

The  Wyckoff  sites  of  the  cubic  phase  get  split  via  condensation  of  various symmetry  adapted  modes  to  give  the  Pnma  space  group  for  the  room  temperature structures of CaTiO .  It has been shown that seven modes are required to arrive at the 3 Pnma  structure  from  the  Pm 3 m  structure 5,6 .  The  irreps  corresponding  to  these  seven modes  are:  octahedral  tilting  modes  R4+  (anti-phase  tilt)  and  M3+  (in-phase  tilt), octahedral distortion modes X5+, R5+ and M2+, and A-site cation displacement modes R5+ and X5+ .  Using the package 'Symmodes' at the Bilbao Crystallographic server , we 7 8 find that the displacive phase  transition from  cubic  space  group  Pm 3 m  to  the orthorhombic space group Pnma requires the splitting of the cubic Wyckoff sites through these vibrational modes in the following manner:

- (i) Ca at 1b site of Pm3m goes to 4c site via R5+(1) and X5+(1) symmodes with coordinates as (0.00+u, 0.25, 0.00-w),
- (ii) O at 3d site goes to O1 at 8d (0.25+u, 0.00+v, 0.75-w) and O2 at 4c (0.50+u, 0.75, 0.00-w) sites via R4+(1), R5+(1), X5+(1), M2+(1), and M3+(1) symmodes and
- (iii) Ti  at  1a  site  remains  undisplaced  from  its  position  and  occupies  4b  site  with coordinates (0.00, 0.00, 0.50),

The  above  representations  of  irreps  are  obtained  from  Bilbao  crystallographic  server. Thus  this process gives us amplitudes  of phonon  modes  using  Rietveld  refined coordinates of low symmetry structure.

## First principle calculations

Density functional theory calculations were performed using ABINIT code  and 9 local  density  approximation  (LDA)  for  the  exchange-correlation  functional.  We  used norm conserving pseudopotentials with an energy cutoff 40 Ha and 6 × 6 × 6 MonkhorstPack mesh  for the sampling of the cubic Brillouin zone. We  approximate  the experimental lattice  constants  of  cubic  structure  for  the  different  values  of  temperature (T) from the volume of orthorhombic unit cells and performed calculations using these approximated cubic lattice constants. The phonon frequencies are the eigen values of the dynamical matrix D:

$$\mathbf D _ { i \alpha j \beta } = ( 1 / \sqrt { m _ { i } } m _ { j } ) \mathbf x \left ( \partial ^ { 2 } \mathbf E / \partial \mathbf x _ { i \alpha } \partial \mathbf x _ { j \beta } \right )$$

where m  is the mass of atom i, x i i α is the displacement of atom i along cartesian direction α and  E  is  total  energy  of  the  system.  The  phonon  dispersion  relations  and  soft  mode eigen vectors are obtained from the linear response calculations in a cubic unit cell. The

phonon frequencies are proportional to the curvature of energy surface as a function of small displacement of atoms from their equilibrium positions (see Eq. 1). The unstable modes, which correspond to the negative curvature of the energy surface, are a signature of  the  structural  phase  transitions.  The  presence  of  unstable  phonon  modes  shows  that system  can  lower  its  energy  from  such  displacement  patterns,  which  also  lowers  the symmetry  of  structure  leading  to  the  structural  phase  transitions.  The  more  is  the curvature  of  energy  surface  for  a  unstable  mode,  larger  is  the  probability  of  system following that displacement pattern.

## Phenomenological theory

R4+ and M3+ are the two modes that lead to the transition of cubic Pm-3m phase to the  Pnma  phase.  R 4+ and  M 3+ modes  correspond  to  tilting  of  octahedra,  and  are  triply degenerate. Let {x 1, x 2, x 3 } and {y1, y 2, y 3 } be the order parameters associated with the R4+ and M3+ modes, respectively. The free energy functional, F(Pm-3m) with Pm-3m as the parent structure and all the six order paramters is given by,

$$& F ( P m - 3 m ) = g _ { 1 2 } ( T ) ( x _ { 1 } ^ { 2 } + x _ { 2 } ^ { 2 } + x _ { 3 } ^ { 2 } ) + g _ { 1 4 } ( x _ { 1 } ^ { 2 } + x _ { 2 } ^ { 2 } + x _ { 3 } ^ { 2 } ) ^ { 2 } + g _ { 2 4 } ( x _ { 1 } ^ { 4 } + x _ { 2 } ^ { 4 } + x _ { 3 } ^ { 4 } ) \\ & + f _ { 1 2 } ( T ) ( y _ { 1 } ^ { 2 } + y _ { 2 } ^ { 2 } + y _ { 3 } ^ { 2 } ) + f _ { 1 4 } ( y _ { 1 } ^ { 2 } + y _ { 2 } ^ { 2 } + y _ { 3 } ^ { 2 } ) ^ { 2 } + f _ { 2 4 } ( y _ { 1 } ^ { 4 } + y _ { 2 } ^ { 4 } + y _ { 3 } ^ { 4 } ) \\ & + h _ { 1 4 } ( x _ { 1 } ^ { 2 } + x _ { 2 } ^ { 2 } + x _ { 3 } ^ { 2 } ) ( y _ { 1 } ^ { 2 } + y _ { 2 } ^ { 2 } + y _ { 3 } ^ { 2 } ) + h _ { 2 4 } ( x _ { 1 } ^ { 2 } y _ { 1 } ^ { 2 } + x _ { 2 } ^ { 2 } y _ { 2 } ^ { 2 } + x _ { 3 } ^ { 2 } y _ { 3 } ^ { 2 } ).$$

Freezing in of the R 4+ mode along (a,0,0)  gives the  I4/mcm phase. Hence, substituting x1 = xc (constant) and neglecting the constant terms in equation S1 we get,

$$& F ( I 4 / m c m ) = ( g _ { 1 2 } ( T ) + 2 g _ { 1 4 } x _ { c } ^ { \ 2 } ) ( x _ { 2 } ^ { \ 2 } + x _ { 3 } ^ { \ 2 } ) + g _ { 1 4 } ( x _ { 2 } ^ { \ 2 } + x _ { 3 } ^ { \ 2 } ) ^ { 2 } + g _ { 2 4 } ( x _ { 2 } ^ { \ 4 } + x _ { 3 } ^ { \ 4 } ) \\ & + ( f _ { 1 2 } ( T ) + 2 h _ { 1 4 } x _ { c } ^ { \ 2 } ) ( y _ { 1 } ^ { \ 2 } + y _ { 2 } ^ { \ 2 } + y _ { 3 } ^ { \ 2 } ) + f _ { 1 4 } ( y _ { 1 } ^ { \ 2 } + y _ { 2 } ^ { \ 2 } + y _ { 3 } ^ { \ 2 } ) ^ { 2 } + f _ { 2 4 } ( y _ { 1 } ^ { \ 4 } + y _ { 2 } ^ { \ 4 } + y _ { 3 } ^ { \ 4 } ) \\ & + h _ { 1 4 } ( x _ { 2 } ^ { \ 2 } + x _ { 3 } ^ { \ 2 } ) ( y _ { 1 } ^ { \ 2 } + y _ { 2 } ^ { \ 2 } + y _ { 3 } ^ { \ 2 } ) + h _ { 2 4 } ( x _ { c } ^ { \ 2 } y _ { 1 } ^ { \ 2 } + x _ { 2 } ^ { \ 2 } y _ { 2 } ^ { \ 2 } + x _ { 3 } ^ { \ 2 } y _ { 3 } ^ { \ 2 } ).$$

Below 1512 K, the R4+ and M3+ modes freeze into the tetragonal I4/mcm structure along (0,a,a) and (a,0,0) respectively, and give rise to the orthorhombic Pnma phase. Thus, the

temperature  dependence  enters  F  through  g 12 (T)=  (T/1512)-1  and  f 12 (T)=  (T/1512)-1. Since the isostructural phase transition involves freezing in of both R 4+ and M 3+ modes along (0,a,a) and (a,0,0) directions respectively, the order parameters x and y are variable for T &lt;1512 K. Hence substituting x =x = x + x , and y1= y + y  and y =y =0 (x  and y 2 3 o o 2 3 o o are constants), equation S2 reduces to,

$$& a v \text{comma} \mathfrak { y }, \text{equation} \mathfrak { z } \text{in} \colon w, \\ & F ( P n m a ) = [ 4 x _ { o } ( g _ { 1 2 } ( T ) + 2 g _ { 1 4 } x _ { c } ^ { 2 } ) + 4 x _ { o } ^ { 3 } ( 4 g _ { 1 4 } + 2 g _ { 2 4 } ) + 4 h _ { 1 4 } x _ { o } y _ { o } ] x \\ & + [ 2 g _ { 1 2 } ( T ) + 2 g _ { 1 4 } x _ { c } ^ { 2 } + 6 x _ { o } ^ { 2 } ( 2 g _ { 1 4 } + 2 g _ { 2 4 } + h 1 4 y _ { o } ^ { 2 } ) ] x ^ { 2 } \\ & + [ 4 x _ { o } ( g _ { 1 4 } + 2 g _ { 2 4 } ) ] x ^ { 3 } + ( 4 g _ { 1 4 } + 2 g _ { 2 4 } ) x ^ { 4 } \\ & + [ ( f _ { 1 4 } + f _ { 2 4 } ) 4 y _ { o } ^ { 3 } + ( f _ { 1 2 } ( T ) + h _ { 1 4 } x _ { c } ^ { 2 } ) ^ { 2 } y _ { o } + 2 h _ { 2 4 } x _ { c } ^ { 2 } y _ { o } + 4 h _ { 1 4 } x _ { o } ^ { 3 } ] y \\ & + [ f _ { 1 2 } ( T ) + h _ { 1 4 } x _ { c } ^ { 2 } + 6 ( f _ { 1 4 } + f _ { 2 4 } ) y _ { o } ^ { 2 } + 2 h _ { 1 4 } x _ { o } ^ { 2 } ] y ^ { 2 } + [ 4 y _ { o } ( f _ { 1 4 } + f _ { 2 4 } ) ] y ^ { 3 } \\ & + ( f _ { 1 4 } + f _ { 2 4 } ) y ^ { 4 } \\ & + ( 8 h _ { 1 4 } x _ { o } y _ { o } ) x y + ( 4 h _ { 1 4 } y _ { o } ) x ^ { 2 } y + ( 4 h _ { 1 4 } x _ { o } ) x y ^ { 2 } + 2 h _ { 1 2 } x ^ { 2 } y ^ { 2 }$$

+

2

2

2

2

14

12

14

14

2

)

4

(

)

4

(

)

8

(

y

x

h

xy

x

h

y

x

y

h

xy

y

x

h

o

o

o

o

+

+

+

## Redefining the Landau coefficients we get,

$$& F = e _ { 1 } x + \frac { a _ { 1 } } { 2 } ( \frac { T } { T _ { c } } + p ) x ^ { 2 } + \frac { b _ { 1 } } { 3 } x ^ { 3 } + \frac { c _ { 1 } } { 4 } x ^ { 4 } \\ & + e _ { 2 } y + \frac { a _ { 2 } } { 2 } ( \frac { T } { T _ { H } } + p ) y ^ { 2 } + \frac { b _ { 2 } } { 3 } y ^ { 3 } + \frac { c _ { 2 } } { 4 } y ^ { 4 } \\ & + d _ { 1 } x ^ { 2 } y + d _ { 2 } x y ^ { 2 } + d _ { 3 } x y + d _ { 4 } x ^ { 2 } y ^ { 2 },$$

Since  the  I4/mcm  to  Pnma  phase  transition  is  relevant  to  equation  S4,  the  critical temperature should be TC = 1512 K. However  all the Landau  coefficients are renormalized, and hence the temperature at which x and y order parameters will freeze into  the  structure  is  no  longer  1512  K.  The  critical  temperature  is  now  replaced  by effective temperature, which is T C= TH= 1000 K. Note that, we fit the T C= TH values to 1000 K to mimic the isostructural phase transition in between 1000-1100K.

Reference :

- 1. J. Rodríguez-Carvajal, Physica B 192 , 55 (1993).

- 2. http://www.ill.eu/sites/fullprof/.
- 3. Aroyo, M. I., Kirov, A., Capillas, C., Perez-Mato, J. M. and Wondratschek. H. Bilbao  Crystallographic  Server  II:  Representations  of  crystallographic  point groups and space groups.  Acta Cryst. A 62 , 115-128 (2006).
- 4. Aroyo,  M.  I.,  Perez-Mato,  J.  M.,  Capillas,  C.  Kroumova,  E.,  Ivantchev,  S., Madariaga, G., Kirov, A. &amp; Wondratschek, H. Bilbao Crystallographic Server I: Databases and crystallographic computing programs. Zeitschrift fuer Kristallographie, 221 , 1, 15-27 (2006).
- 5. C. N. W. Darlington, Acta. Cryst. A 58 , 66(2002(a)).
- 6. C. N. W. Darlington, Acta. Cryst. A 58 , 299(2002(b)).
- 7. K. S. Knight and C. N. W. Darlington, Powder Diffr 20 , 7(2005).
- 8. C. Capillas, E. Kroumova, M. I. Aroyo, J. M. Perez-Mato, H. T. Stokes, and D. M. Hatch, J. Appl. Cryst. 36 , 953(2003).
- 9. http://www.abinit.org/.
Saurabh Tripathi Anil Kumar Sharmila Shirodkar Masatomo Yashima U. V. Waghmare Dhananjai Pandey
Discovery of an isostructural phase transition within orthorhombic phase field of CaTiO3
Earths lower mantle extending from 670 to 2,990 km deep is predominantly
composed of a perovskite-type (Mg,Fe)SiO3 phase1,2. The perovskite phase
undergoes a structural phase transition to a post-perovskite phase responsible
for D" layer seismic discontinuity2,3 at about 2690 km depth in the lowermost
region of the lower mantle. However, structural basis of other seismic
discontinuities occurring in the upper region of the lower mantle (700 km to
1,200 km deep) remains unexplained4-7, as no apparent change in the crystal
symmetry of the orthorhombic perovskite phase has been reported5. We present
here unambiguous evidence for a non-apparent isostructural phase transition8 in
the stable orthorhombic perovskite phase of CaTiO3 which may have relevance to
phase transitions in the perovskite phase of (Mg,Fe)SiO3 also, as both the
compounds have similar structure, tolerance factor and thermochemical
properties9-11. Our results are based on the analysis of neutron powder
diffraction patterns using Rietveld and mode crystallography techniques and are
supported by density functional and Landau theory calculations. The present
results on CaTiO3 would encourage search for isostructural phase transition in
the perovskite phase of (Mg,Fe)SiO3 that may provide clue to the unexplained
geophysical phenomena in the upper part of the earths lower mantle.