Atomistic study of finite temperature properties in ferroelectric BiAlO$_3$

B. K. Mani; Chandan Kumar Vishwakarma; Indranil Mal; I. Ponomareva; Mohd Zeeshan

Atomistic study of finite temperature properties in ferroelectric BiAlO$_3$

T0 review · 3 major / 0 minor · reviewed 2026-07-11 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read First-principles Monte Carlo simulations of BiAlO3 predict a cubic-to-rhombohedral ferroelectric transition at 1160 K with 81 μC/cm² spontaneous polarization and strong ferroelectric–antiferrodistortive coupling under temperature, field, an

desk verdict Solid first finite-T map of BiAlO3 with standard effective-Hamiltonian MC; useful for the subfield, but the same LDA parameterization that overestimates TC by ~370 K undercuts the quantitative stress/field diagrams. read the letter →

arxiv 2607.03842 v1 pith:5VMCE43N submitted 2026-07-04 cond-mat.mtrl-sci

Indranil Mal , Chandan Kumar Vishwakarma , Mohd Zeeshan , I. Ponomareva , B. K. Mani This is my paper · ORCID

classification cond-mat.mtrl-sci PACS 77.80.-e77.84.-s64.70.K

keywords BiAlO3lead-freeferroelectricseffectiveHamiltonianMonteCarlosimulationrhombohedralR3cantiferrodistortivecouplinghydrostaticpressurestress-inducedphasetransitions

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper builds a first-principles effective Hamiltonian for the lead-free perovskite BiAlO3 and uses Monte Carlo simulations to predict its finite-temperature ferroelectric behavior. The model shows a single transition from cubic paraelectric Pm-3m to rhombohedral ferroelectric R3c at 1160 K, with spontaneous polarization of about 81 μC/cm² along [111] at low temperature. Ferroelectric and oxygen-octahedra tilt (antiferrodistortive) modes stay tightly coupled across temperature and electric field. Hydrostatic pressure lowers both polarization and the Curie temperature, while uniaxial and biaxial stresses open additional intermediate phases. The work supplies the missing finite-temperature atomistic picture for a material whose experimental polarization and transition temperature have been hard to reconcile with zero-kelvin calculations, and it maps how external loads can retune those properties for device use.

What carries the argument

First-principles effective Hamiltonian whose degrees of freedom are the ferroelectric local mode u, the antiferrodistortive oxygen-octahedra rotation ω (R25), and homogeneous strain η, with all interaction parameters fixed by LDA-DFT; the energy is sampled by single-flip Metropolis Monte Carlo annealing on 20×20×20 supercells.

What would settle it

A high-quality bulk single-crystal or carefully prepared ceramic measurement of the Curie temperature and saturated room-temperature polarization of BiAlO3 that either confirms values near 1160 K and 50–80 μC/cm² or shows that the true bulk values remain near the lower experimental range (~800 K and <30 μC/cm²).

Watch

Extended reading notes

Core claim

An effective Hamiltonian parameterized from density-functional theory and sampled by Monte Carlo annealing predicts that BiAlO3 undergoes a single first-order-like transition from cubic Pm-3m to rhombohedral R3c at TC = 1160 K, develops a spontaneous polarization of ~81 μC/cm² along [111], and maintains strong coupling between the polar local mode and the R25 antiferrodistortive tilts under temperature, electric field, hydrostatic pressure, and uniaxial/biaxial stress.

Load-bearing premise

The truncated effective Hamiltonian (only polar Γ modes, R25 tilts, and homogeneous strain, with parameters fixed by zero-temperature LDA) remains quantitatively accurate at finite temperature and under stress, even though it already overestimates the experimental Curie temperature by several hundred kelvin.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

This is the first finite-temperature atomistic treatment of BiAlO3. They build a Zhong–Vanderbilt–Rabe effective Hamiltonian (Γ polar + R25 AFD + homogeneous strain) from LDA phonons and energy landscapes, then run Metropolis MC on a 20^{3} supercell. Ground state is correctly R3c, FE and AFD stay locked along [111], and they deliver the previously missing T-dependence of PS, hysteresis, coercive field, strain, plus hydrostatic and uniaxial/biaxial stress phase sequences. That fills a real gap; prior work stopped at 0 K DFT or fragmentary experiment.

What they do well is transparent and standard. Phonon analysis, double-well maps, and the full parameter table are clear. They openly note the TC overestimate (1160 K vs experimental lower bound >793 K) and the large PS discrepancy with measured values, attributing both to LDA binding energy and missing defects. FE–AFD coupling under field is cleanly shown, and the stress-induced multi-phase sequences (cubic–T–R under tensile uniaxial, opposite under biaxial) are new predictions worth having on the table.

The soft spot is exactly the one the stress-test flags, and it is real but not fatal. Every coupling coefficient is frozen at zero-K LDA; nothing is re-tuned or cross-checked against finite-T data. So the same model that already misses TC by hundreds of kelvin is used to claim quantitative field and stress diagrams. That weakens the “accurately predict” language in the abstract, but the qualitative trends (pressure suppresses TC, stress splits the transition) are still useful. Missing error bars and no released code are minor annoyances, not structural flaws. Circularity is low: parameters come from independent DFT, observables are genuine MC output.

This is for people who already work on lead-free perovskites or effective-Hamiltonian modeling. A serious referee should see it; the method is established, the gap is genuine, and the limitations are stated. I would cite the phase sequences and the FE–AFD locking if I needed BiAlO3 numbers, with the usual LDA-caveat. Engage.

Referee Report

3 major / 0 minor

Summary. The manuscript constructs a first-principles effective Hamiltonian for BiAlO3 that retains the polar Γ local mode, the R25 antiferrodistortive (AFD) mode and homogeneous strain (Eq. 1), with all coefficients obtained from LDA-DFT (Table I). Single-flip Metropolis Monte Carlo annealing on a 20×20×20 supercell is then used to map the finite-temperature phase diagram, electric-field hysteresis and the response to hydrostatic pressure and uniaxial/biaxial stress. The simulations recover a single first-order-like cubic (Pm-3m) to rhombohedral (R3c) transition at TC = 1160 K, a low-temperature spontaneous polarization of ~81 µC cm-2 along [111], strong FE–AFD coupling that persists under field, and a sequence of stress-induced multi-phase transitions (Fig. 5).

Significance. A transferable, first-principles-based atomistic model for BiAlO3 fills a genuine gap: previous DFT work was limited to T = 0 K and experimental polarization values scatter by nearly an order of magnitude. The predicted FE–AFD coupling, remanent strain and stress-tunable phase sequence are concrete, falsifiable results that can guide thin-film and high-pressure experiments on this lead-free candidate. The parameterization follows the established Zhong–Vanderbilt–Rabe protocol and is fully tabulated, so the work is in principle reproducible.

major comments (3)

Abstract and §III claim that the model “accurately predict[s]” finite-temperature properties, yet TC = 1160 K exceeds the experimental lower bound (>793 K) by ~370 K and the room-temperature PS (~58 µC cm-2 under field) is several times larger than measured values (9–29 µC cm-2). The authors attribute the discrepancy to overestimated binding energy and missing defects, but never quantify the sensitivity of the stress- and field-dependent diagrams (Figs. 3–5) to this systematic error. A short discussion or re-scaling test that shows which qualitative features survive a ~30 % reduction in the energy scale is needed before the diagrams can be presented as predictive.
§II and Fig. 1 establish that the global energy minimum is obtained only when Γ and R25 are frozen together, yet the subsequent Monte Carlo runs never report the free-energy barrier or the order of the transition beyond a visual inspection of the order-parameter jump. Because the multi-phase sequences under uniaxial/biaxial stress (Fig. 5b,c) rest on the relative stability of these modes, a quantitative measure of the FE–AFD coupling strength (e.g., the Gxxxx, Gxxyy, Gxyxy terms of Table I evaluated along the transition path) would strengthen the central claim that the coupling remains strong over a broad temperature and field range.
The electric-field and stress results (Figs. 3–5) are obtained with a single 20×20×20 supercell and no finite-size or statistical-error analysis. Near the first-order-like transition and under stress, domain-wall or nucleation effects can shift the apparent coercive fields and the TC1/TC2 boundaries by tens of kelvin. At least a brief check with a second supercell size (or block-averaging error bars on the order parameters) is required to confirm that the reported multi-phase sequences are not finite-size artifacts.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: DFT-derived effective-Hamiltonian parameters feed independent Monte Carlo predictions of finite-T order parameters, TC, hysteresis and stress sequences.

full rationale

The derivation chain is the standard first-principles effective-Hamiltonian workflow. Phonon instabilities and energy landscapes (Fig. 1) are obtained from independent zero-temperature LDA-DFT; the resulting coefficients (Table I) completely fix Eq. (1). Subsequent single-flip Metropolis Monte Carlo annealing and field/stress sweeps then generate the finite-temperature observables (TC = 1160 K, PS(T), hysteresis loops, pressure/stress phase sequences). These observables are not among the fitted inputs, nor are they forced by construction from the zero-K energy surfaces. Self-citations ([29–31] and earlier Mani/Ponomareva applications) merely point to the established methodological framework; none of them supplies a uniqueness theorem or an ansatz that closes the loop on the BiAlO3 results. Discrepancies with experiment (overestimated TC, larger PS) are openly noted and attributed to missing defects and LDA binding-energy error, confirming that the finite-T numbers are genuine (if imperfect) predictions rather than re-statements of the inputs. Score 1 reflects only the routine methodological self-citation, which is not load-bearing.

Assumptions & free parameters 2 free parameters · 3 assumptions · 0 invented entities

The central finite-T predictions rest on (i) the standard effective-Hamiltonian energy functional whose functional form is taken from the literature, (ii) a large set of numerical coefficients extracted once and for all from LDA-DFT, and (iii) the usual Monte Carlo sampling assumptions. No new physical entities are postulated; the free parameters are the DFT-derived coefficients themselves.

free parameters (2)

Full set of FE, AFD, elastic and coupling coefficients (κ2, α, γ, j1–j7, B11, B12, B44, B1xx …, Gxxxx …, Z*, ε∞, etc.) = see Table I (atomic units)
Approximately 30 numerical values listed in Table I; each is obtained by fitting the effective-Hamiltonian energy surface to a series of LDA-DFT total-energy calculations on frozen-phonon and strained supercells. The finite-T observables inherit whatever systematic bias is present in those DFT numbers.
Cubic reference lattice constant a0 = 7.03225 a.u.
Fixed at the LDA value 7.03225 a.u.; all subsequent strain and polarization conversions depend on this choice.

assumptions (3)

domain assumption The low-energy physics of BiAlO3 is completely captured by the local polar mode u, the R25 AFD mode ω, and homogeneous strain η; higher-order or zone-boundary modes may be neglected.
Stated in Sec. II after the phonon analysis; the energy functional (Eq. 1) is written solely in these variables.
domain assumption LDA exchange-correlation functional yields sufficiently accurate energy surfaces and Born effective charges for quantitative finite-T predictions.
All parameters in Table I are computed with LDA; the known LDA underestimation of lattice constants and overestimation of soft-mode frequencies is accepted without correction.
domain assumption Metropolis Monte Carlo sampling of a 20 imes20 imes20 supercell with the stated annealing schedule converges to the thermodynamic equilibrium state.
Methodology paragraph; finite-size and equilibration errors are not quantified.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Atomistic study of finite temperature properties in ferroelectric BiAlO$_3$." pith.science (2026). https://pith.science/paper/5VMCE43N

@misc{pith2026260703842,
  author       = {Pith},
  title        = {Pith review of: Atomistic study of finite temperature properties in ferroelectric BiAlO$_3$},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5VMCE43N}},
  note         = {Machine review of arXiv:2607.03842}
}

abstract

The lead-free perovskite ferroelectrics captivate researchers with their unique functional properties leading to important technological applications. In the search for a new lead-free perovskite of technological importance, we develop a first-principles based atomistic model to accurately predict the properties of BiAlO$_3$ in experimentally relevant conditions. Consistent with the experimental observations, our simulations predict a rhombohedral ferroelectric ($R3c$) ground state for BiAlO$_3$ facilitated by a structural phase transition from paraelectric (cubic, $Pm\Bar{3}m$) phase. The room-temperature spontaneous polarization and Curie temperature are obtained to be 81 $\mu$C/cm$^2$ (along [111] direction) and 1160 K, respectively. Our simulations reveal strong coupling between ferroelectric and antiferrodistortive modes for a broad spectrum of temperature and electric field. We find that hydrostatic pressure suppresses both spontaneous polarization and Curie temperature, while both uniaxial and biaxial stresses induce multiple phase transitions in BiAlO$_3$.

Figures

Figures reproduced from arXiv: 2607.03842 by the authors.

**Figure 2.** FIG. 2. Temperature evolution of FE order parameter (panel ( [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** FIG. 3. Electric field evolution of polarization (panels (a) [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

Figures from the paper (2 more)

**Figure 4.** Figure 4: FIG. 4. (a) Temperature dependent P-E hysteresis for fields a [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]

**Figure 5.** Figure 5: FIG. 5. Evolution of polarization as a function of temperatu [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

53 extracted references · 3 canonical work pages

[1]

Our simulations reveal strong coupling between ferroelec- tric and antiferrodistortive modes for a broad spectrum of t emperature and electric ﬁeld

direction) and 1160 K, respectively. Our simulations reveal strong coupling between ferroelec- tric and antiferrodistortive modes for a broad spectrum of t emperature and electric ﬁeld. We ﬁnd that hydrostatic pressure suppresses both spontaneous pol arization and Curie temperature, while both uniaxial and biaxial stresses induce multiple phase tr ansitio...

arXiv 2026
[2]

A dense k-mesh of 11 × 11 × 11 and an energy cutoﬀ of 600 eV for plane wave basis were used in the calculations

method to represent the atomic cores and the local- density-approximation (LDA) [35, 36] as the exchange- correlation functional. A dense k-mesh of 11 × 11 × 11 and an energy cutoﬀ of 600 eV for plane wave basis were used in the calculations. All self-consistent-ﬁeld calcula- tions were allowed to converge till the energy diﬀerence was less than 10 − 6 eV...
[3]

In terms of the magnitude of PSat, there is a slight change in the val- ues for the ﬁelds along diﬀerent directions

direction, we observed hysteresis for x-component of polarization, whereas polarization vectors along y− and z-directions remain unaﬀected by the ﬁeld. In terms of the magnitude of PSat, there is a slight change in the val- ues for the ﬁelds along diﬀerent directions. We obtained PSat of 57.8, 59.2, and 60.2 µ C/cm2 for the ﬁelds along [111], [110], and [...
[4]

For this, we simulated the uniaxial and biaxial stresses in the range − 5

direction. For this, we simulated the uniaxial and biaxial stresses in the range − 5. 0 GPa to 5 . 0 GPa in a step of 1 . 0 GPa. The eﬀect of uniaxial stress was simu- lated by applying a single component of the stress tensor (σ ), whereas for biaxial stress, we ﬁxed σ1 = σ2 = σ and σ6 = 0, while allowing all other components to re- lax. We ﬁnd, under the...

2022
[5]

L. Liao, H. J. Fan, B. Yan, Z. Zhang, L. L. Chen, B. S. Li, G. Z. Xing, Z. X. Shen, T. Wu, X. W. Sun, J. Wang, and T. Yu, ACS Nano 3, 700 (2009), publisher: Ameri- can Chemical Society

2009
[6]

H. Xu, F. Sun, E. Li, W. Guo, L. Hua, R. Wang, W. Li, J. Chu, W. Liu, J. Luo, and Z. Sun, Advanced materials , e2414339 (2024)

2024
[7]

R. Guo, L. E. Cross, S.-E. Park, B. Noheda, D. E. Cox, and G. Shirane, Phys. Rev. Lett. 84, 5423 (2000)

2000
[8]

You, W.-Q

Y.-M. You, W.-Q. Liao, D. Zhao, H.-Y. Ye, Y. Zhang, Q. Zhou, X. Niu, J. Wang, P.-F. Li, D.-W. Fu, Z. Wang, S. Gao, K. Yang, J.-M. Liu, J. Li, Y. Yan, and R.-G. Xiong, Science 357, 306 (2017), https://www.science.org/doi/pdf/10.1126/science.aai8535

work page doi:10.1126/science.aai8535 2017

Show all 53 references

[9]

P. M. Rørvik, T. Grande, and M. Einarsrud, Advanced Materials 23, 10.1002/adma.201004676 (2011)

2011 doi
[10]

B. Xu, N. G. Pai, and L. Cross, Materials Letters 34, 157 (1998)

1998
[11]

Ignatiev, Y

A. Ignatiev, Y. Xu, N. Wu, and D. Liu, Materials Science and Engineering: B 56, 191 (1998)

1998
[12]

Scott, Annual Review of Materials Research 41, 229 (2011)

J. Scott, Annual Review of Materials Research 41, 229 (2011)

2011
[13]

A. S. Mischenko, Q. Zhang, J. F. Scott, R. W. What- more, and N. D. Mathur, Science 311, 1270 (2006), https://www.science.org/doi/pdf/10.1126/science.1123811

2006 doi
[14]

T. F. Zhang, X. G. Tang, Q. X. Liu, Y. P. Jiang, X. X. Huang, and Q. F. Zhou, Journal of Physics D: Applied Physics 49, 095302 (2016)

2016
[15]

H. Zhao, J. Miao, L. Zhang, Y. Rong, J. Chen, J. Deng, R. Yu, J. Cao, H. Wang, and X. Xing, Dalton Trans. 45, 1554 (2016)

2016
[16]

H. S. Bhatti, S. T. Hussain, F. A. Khan, and S. Hussain, Applied Surface Science 367, 291 (2016)

2016
[17]

A. K. Yadav, Anita, S. Kumar, A. Panchwanee, V. R. Reddy, P. M. Shirage, S. Biring, and S. Sen, RSC Adv. 7, 39434 (2017)

2017
[18]

Jaouen, A

N. Jaouen, A. C. Dhaussy, J. P. Iti´ e, A. Rogalev, S. Marinel, and Y. Joly, Phys. Rev. B 75, 224115 (2007)

2007
[19]

S´ aghi-Szab´ o, R

G. S´ aghi-Szab´ o, R. E. Cohen, and H. Krakauer, Phys. Rev. Lett. 80, 4321 (1998)

1998
[21]

Chung, W.-J

J.-K. Chung, W.-J. Kim, J. K. Kim, S. S. Kim, and T. K. Song, Integrated Ferroelectrics 87, 25 (2007)

2007
[22]

R. E. Eitel, C. A. Randall, T. R. Shrout, and S.-E. Park, Japanese Journal of Applied Physics 41, 2099 (2002)
[23]

G. P. V. Yu. F. Popov, A. M. Kadomtseva and A. K. Zvezdin, Ferroelectrics 162, 135 (1994)

1994
[24]

Zylberberg, A

J. Zylberberg, A. A. Belik, E. Takayama-Muromachi, and Z.-G. Ye, Chemistry of Materials 19, 6385 (2007)

2007
[25]

Z. Li, H. Liu, Z. Yao, J. Xie, X. Li, C. Diao, A. Ullah, H. Hao, and M. Cao, Ceramics International 45, 22523 (2019)

2019
[26]

A. A. Belik, T. Wuernisha, T. Kamiyama, K. Mori, M. Maie, T. Nagai, Y. Matsui, and E. Takayama- Muromachi, Chemistry of Materials 18, 133 (2006)

2006
[27]

Mangalam, S

R. Mangalam, S. Bhat, A. Iyo, Y. Tanaka, A. Sundaresan, and C. Rao, Solid State Communications 146, 435 (2008)

2008
[28]

J. Y. Son, C. S. Park, and Y.-H. Shin, Applied Physics Letters 92, 222911 (2008)

2008
[29]

Baettig, C

P. Baettig, C. F. Schelle, R. LeSar, U. V. Waghmare, and N. A. Spaldin, Chemistry of Materials 17, 1376 (2005)

2005
[30]

P/suppress lowa´ s-Korus and J

I. P/suppress lowa´ s-Korus and J. Kaczkowski, New Journal of Chemistry 46, 15381 (2022)

2022
[31]

Kaczkowski, Materials Chemistry and Physics 177, 405 (2016)

J. Kaczkowski, Materials Chemistry and Physics 177, 405 (2016)

2016
[32]

Kang, Current Applied Physics 35, 67 (2022)

Y. Kang, Current Applied Physics 35, 67 (2022)

2022
[33]

Zhong, D

W. Zhong, D. Vanderbilt, and K. M. Rabe, Phys. Rev. B 52, 6301 (1995)

1995
[34]

Vanderbilt and W

D. Vanderbilt and W. Zhong, Ferroelectrics 206, 181 (1998)

1998
[35]

B. K. Mani, S. Lisenkov, and I. Ponomareva, Phys. Rev. B 91, 134112 (2015)

2015
[36]

Kresse and J

G. Kresse and J. Furthm¨ uller, Computational Materials Science 6, 15 (1996)

1996
[37]

Kresse and J

G. Kresse and J. Furthm¨ uller, Phys. Rev. B 54, 11169 (1996)

1996
[38]

P. E. Blochl, Phys. Rev. B 50, 17953 (1994)

1994
[39]

D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 (1980)

1980
[40]

J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981)

1981
[41]

B. K. Mani, C.-M. Chang, and I. Ponomareva, Phys. Rev. B 88, 064306 (2013)

2013
[42]

Cisse, W

C. Cisse, W. Zaki, and T. Ben Zineb, International Journal of Plasticity 76, 244 (2016)

2016
[43]

W. He, Q. Li, Q. Yan, N. Luo, Y. Zhang, X. Chu, and D. Shen, Crystals 4, 262 (2014)

2014
[44]

Jang, H.-S

C.-H. Jang, H.-S. Kim, H. Kim, and H.-Y. Cha, Materials 15, 2097 (2022)
[45]

Patel, A

S. Patel, A. Chauhan, A. Chauhan, and R. Vaish, Materials Research Express 2, 035501 (2015)

2015
[46]

Zhang, Z

Y. Zhang, Z. Chen, W. Cao, and Z. Zhang, Applied Physics Letters 111, 172902 (2017)

2017
[47]

X. Zhao, J. Y. Dai, J. Wang, H. L. W. Chan, C. L. Choy, X. M. Wan, and H. S. Luo, Phys. Rev. B 72, 064114 (2005)

2005
[48]

H. Qiao, C. He, Z. Wang, D. Pang, X. Li, Y. Liu, and X. Long, RSC Adv. 7, 32607 (2017)

2017
[49]

R. Liu, C. Gu, X. Tian, F. Zheng, X. Jiang, and Y. Zhang, Ceramics International 48, 22411 (2022)

2022
[50]

Ducharme, V

S. Ducharme, V. M. Fridkin, A. V. Bune, S. P. Palto, L. M. Blinov, N. N. Petukhova, and S. G. Yudin, Phys. Rev. Lett. 84, 175 (2000)

2000
[51]

Kashikar, S

R. Kashikar, S. Lisenkov, and I. Ponomareva, Phys. Rev. B 109, L020101 (2024)

2024
[52]

Mohan and A

Y. Mohan and A. Arockiarajan, AIP Advances 6, 035311 (2016)

2016
[53]

A. L. Narit Funsueb and A. Ngamjarurojana, Phase Transitions 93, 678 (2020)

2020
[54]

McCash, B

K. McCash, B. K. Mani, C.- M. Chang, and I. Ponomareva, Journal of Physics: Condensed Matter 26, 435901 (2014)

2014

Pith tools

Reviewed July 11, 2026 · model on record in the stance chip above.

[1] [1]

Our simulations reveal strong coupling between ferroelec- tric and antiferrodistortive modes for a broad spectrum of t emperature and electric ﬁeld

direction) and 1160 K, respectively. Our simulations reveal strong coupling between ferroelec- tric and antiferrodistortive modes for a broad spectrum of t emperature and electric ﬁeld. We ﬁnd that hydrostatic pressure suppresses both spontaneous pol arization and Curie temperature, while both uniaxial and biaxial stresses induce multiple phase tr ansitio...

arXiv 2026

[2] [2]

A dense k-mesh of 11 × 11 × 11 and an energy cutoﬀ of 600 eV for plane wave basis were used in the calculations

method to represent the atomic cores and the local- density-approximation (LDA) [35, 36] as the exchange- correlation functional. A dense k-mesh of 11 × 11 × 11 and an energy cutoﬀ of 600 eV for plane wave basis were used in the calculations. All self-consistent-ﬁeld calcula- tions were allowed to converge till the energy diﬀerence was less than 10 − 6 eV...

[3] [3]

In terms of the magnitude of PSat, there is a slight change in the val- ues for the ﬁelds along diﬀerent directions

direction, we observed hysteresis for x-component of polarization, whereas polarization vectors along y− and z-directions remain unaﬀected by the ﬁeld. In terms of the magnitude of PSat, there is a slight change in the val- ues for the ﬁelds along diﬀerent directions. We obtained PSat of 57.8, 59.2, and 60.2 µ C/cm2 for the ﬁelds along [111], [110], and [...

[4] [4]

For this, we simulated the uniaxial and biaxial stresses in the range − 5

direction. For this, we simulated the uniaxial and biaxial stresses in the range − 5. 0 GPa to 5 . 0 GPa in a step of 1 . 0 GPa. The eﬀect of uniaxial stress was simu- lated by applying a single component of the stress tensor (σ ), whereas for biaxial stress, we ﬁxed σ1 = σ2 = σ and σ6 = 0, while allowing all other components to re- lax. We ﬁnd, under the...

2022

[5] [5]

L. Liao, H. J. Fan, B. Yan, Z. Zhang, L. L. Chen, B. S. Li, G. Z. Xing, Z. X. Shen, T. Wu, X. W. Sun, J. Wang, and T. Yu, ACS Nano 3, 700 (2009), publisher: Ameri- can Chemical Society

2009

[6] [6]

H. Xu, F. Sun, E. Li, W. Guo, L. Hua, R. Wang, W. Li, J. Chu, W. Liu, J. Luo, and Z. Sun, Advanced materials , e2414339 (2024)

2024

[7] [7]

R. Guo, L. E. Cross, S.-E. Park, B. Noheda, D. E. Cox, and G. Shirane, Phys. Rev. Lett. 84, 5423 (2000)

2000

[8] [8]

You, W.-Q

Y.-M. You, W.-Q. Liao, D. Zhao, H.-Y. Ye, Y. Zhang, Q. Zhou, X. Niu, J. Wang, P.-F. Li, D.-W. Fu, Z. Wang, S. Gao, K. Yang, J.-M. Liu, J. Li, Y. Yan, and R.-G. Xiong, Science 357, 306 (2017), https://www.science.org/doi/pdf/10.1126/science.aai8535

work page doi:10.1126/science.aai8535 2017