Emergent magnetic and charge ordered phases in freestanding ultrathin ce{LaVO3}
Pith reviewed 2026-05-22 03:45 UTC · model grok-4.3
The pith
Freestanding ultrathin LaVO3 films show thickness-dependent magnetic transitions from stripe antiferromagnetism in monolayers to bulk-like C-type antiferromagnetism, with surface charge transfer creating ferromagnetic states beyond four van
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the intrinsic magnetic order in freestanding LaVO3 films depends on thickness, recovering the bulk C-AFM configuration only for thicker films, while polar-catastrophe-driven surface doping beyond four layers produces mixed stripe-AFM and ferromagnetic surface states, and that hole doping of the monolayer drives a ferromagnetic state accompanied by striped charge order that changes from a simple stripe pattern at 0.5 holes per formula unit to a 3:1 stripe pattern at 0.25 holes per formula unit.
What carries the argument
Polar catastrophe driven charge transfer that dopes surface layers and thereby stabilizes new magnetic and charge-ordered states.
If this is right
- Monolayer LaVO3 adopts stripe antiferromagnetic order.
- Films thicker than four layers develop ferromagnetic surface states due to charge transfer while interior layers remain bulk-like C-AFM.
- Hole doping of the monolayer stabilizes ferromagnetism.
- Charge ordering appears as a striped pattern at 0.5 holes per formula unit and a 3:1 stripe pattern at 0.25 holes per formula unit.
Where Pith is reading between the lines
- Varying film thickness could be used to engineer surface-dominated magnetic states in oxide membranes.
- The doping-induced charge-ordering patterns suggest that similar superstructures might appear under gate-controlled doping in related vanadates.
- The results imply that freestanding films of other polar perovskites could exhibit analogous surface reconstruction effects.
Load-bearing premise
Standard density-functional-theory approximations are accurate enough to identify the correct magnetic and charge-ordered ground states in these freestanding films.
What would settle it
Experimental measurement of magnetic order in a freestanding monolayer LaVO3 film that fails to show stripe antiferromagnetism, or in a five-layer film that fails to show ferromagnetic surface states, would falsify the predicted sequence.
Figures
read the original abstract
Transition metal oxide perovskites are an ideal platform for exploring the interplay between spin, orbital, charge and lattice degrees of freedom. Among them, \ce{LaVO3} has been extensively studied in heterostructures and superlattices, where exotic phases have been reported. Motivated by the advances in freestanding oxide membranes, we investigate the intrinsic properties of freestanding ultrathin \ce{LaVO3} films using density functional theory. Our calculations reveal a sequence of magnetic phase transitions with thickness, starting from stripe-AFM in monolayer until the bulk like C-AFM is recovered. Beyond four layers, polar catastrophe driven charge transfer dopes the surface layers giving rise to stripe-AFM and ferromagnetic surface states while the central layers remain bulk like. We further explore this fact by studying charge doped monolayer, discovering that hole doping drives the system into ferromagnetic state. Doping also induced charge ordering in the system. A striped charge ordering pattern is observed at 0.5 h/fu, while a 3:1 stripe pattern emerges at 0.25 h/fu, indicating that the periodicity of the superstructure changes with doping concentration.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses density-functional theory to study freestanding ultrathin LaVO3 films. It reports a thickness-driven sequence of magnetic ground states, beginning with stripe antiferromagnetic order in the monolayer and recovering bulk-like C-type antiferromagnetic order for thicker films. Beyond approximately four layers, polar-catastrophe-induced charge transfer is claimed to dope the surface layers, stabilizing stripe-AFM and ferromagnetic surface states while the interior remains bulk-like. Separate calculations on charge-doped monolayers show that hole doping drives the system ferromagnetic and induces charge ordering, with a striped pattern at 0.5 holes per formula unit and a 3:1 stripe pattern at 0.25 holes per formula unit.
Significance. If the reported phase sequence and doping-induced orders prove robust, the work would illustrate how reduced dimensionality and polar electrostatics can generate emergent magnetic and charge-ordered states in LaVO3 membranes, complementing existing heterostructure studies and offering testable predictions for experiments on freestanding oxide films.
major comments (2)
- [Results on thickness dependence and polar catastrophe] The central thickness-dependent magnetic crossover (monolayer stripe-AFM to bulk-like C-AFM beyond four layers) and the surface states induced by polar-catastrophe charge transfer rest on total-energy comparisons performed within a single DFT+U setup. No scan of the Hubbard U parameter, no comparison of different exchange-correlation functionals, and no reported energy differences with error bars are provided, even though the relative stability of C-AFM, stripe-AFM, and FM states in V^{3+} (d^{2}) Mott insulators is known to shift by tens of meV per formula unit when U is varied by 1–2 eV. This directly affects the reliability of the reported four-layer threshold and the assignment of surface magnetic states.
- [Charge-doped monolayer calculations] The charge-ordering patterns reported for the doped monolayer (striped order at 0.5 h/fu and 3:1 stripe at 0.25 h/fu) are obtained from the same unbenchmarked DFT+U framework. Without a demonstration that these patterns survive changes in U or functional, or without comparison to higher-level methods such as DMFT, the specific periodicities and the claim that doping drives ferromagnetism remain sensitive to the chosen approximation.
minor comments (2)
- [Abstract and Computational Methods] The abstract and methods section should explicitly state the exchange-correlation functional, the value of U (and J), the k-point mesh, and the slab vacuum thickness used for the freestanding films.
- [Figures showing spin and charge densities] Figure captions for the magnetic and charge-order configurations should include the computed energy differences relative to the lowest state so that the stability margins are immediately visible.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review of our manuscript. We address each major comment in detail below, indicating the revisions we will implement to improve the robustness of our findings.
read point-by-point responses
-
Referee: [Results on thickness dependence and polar catastrophe] The central thickness-dependent magnetic crossover (monolayer stripe-AFM to bulk-like C-AFM beyond four layers) and the surface states induced by polar-catastrophe charge transfer rest on total-energy comparisons performed within a single DFT+U setup. No scan of the Hubbard U parameter, no comparison of different exchange-correlation functionals, and no reported energy differences with error bars are provided, even though the relative stability of C-AFM, stripe-AFM, and FM states in V^{3+} (d^{2}) Mott insulators is known to shift by tens of meV per formula unit when U is varied by 1–2 eV. This directly affects the reliability of the reported four-layer threshold and the assignment of surface magnetic states.
Authors: We agree that the sensitivity of the magnetic ground states to the Hubbard U parameter requires explicit demonstration. Our original calculations employed the PBE+U functional with U = 4 eV, a value consistent with prior literature on LaVO3. To address this point, we have performed additional total-energy comparisons for monolayer, trilayer, and five-layer films using U = 3 eV and U = 5 eV. The sequence of magnetic phases and the recovery of bulk-like C-AFM order beyond approximately four layers remains qualitatively unchanged, although the precise layer threshold shifts by at most one layer. We will add a supplementary figure showing the U-dependent energy differences (including statistical uncertainties from k-point convergence) and a short discussion of why PBE+U is appropriate given existing benchmarks on bulk LaVO3. Hybrid-functional calculations remain computationally prohibitive for the supercell sizes employed, but we will reference relevant literature comparisons. revision: yes
-
Referee: [Charge-doped monolayer calculations] The charge-ordering patterns reported for the doped monolayer (striped order at 0.5 h/fu and 3:1 stripe at 0.25 h/fu) are obtained from the same unbenchmarked DFT+U framework. Without a demonstration that these patterns survive changes in U or functional, or without comparison to higher-level methods such as DMFT, the specific periodicities and the claim that doping drives ferromagnetism remain sensitive to the chosen approximation.
Authors: We concur that the doping-induced ferromagnetic state and charge-ordering patterns must be tested for robustness. We have carried out additional calculations for the 0.25 and 0.5 hole-doped monolayers at U = 3 eV and U = 5 eV. Both the ferromagnetic ground state and the striped (0.5 h/fu) and 3:1 (0.25 h/fu) charge-ordering patterns persist, with energy differences varying by less than 15 meV per formula unit. These results, together with a brief analysis of how superstructure periodicity evolves with doping, will be incorporated into the revised manuscript. A comprehensive DMFT study of the doped freestanding films lies outside the scope of the present work owing to the large supercells required; we will instead cite existing DMFT literature on doped vanadates to contextualize the DFT+U findings. revision: yes
Circularity Check
No circularity: DFT total-energy results are independent of inputs by construction
full rationale
The paper reports direct outputs from standard DFT+U total-energy comparisons for different magnetic configurations and doping levels in LaVO3 slabs. These include thickness-dependent transitions from stripe-AFM to C-AFM and doping-induced FM plus charge-order patterns. No step reduces a claimed prediction to a fitted parameter, self-citation, or ansatz that is defined in terms of the target result. The methodology is self-contained against external benchmarks (experimental bulk properties or higher-level methods), with no load-bearing self-citation chain or renaming of known results. This is the normal honest case for a computational materials paper.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Density functional theory with standard approximations suffices to determine the magnetic and charge-ordered ground states of LaVO3.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
An effective U of 3 eV is considered, consistent with prior studies on vanadates... We compute the isotropic exchange interactions through energy mapping of various spin-ordered phases.
-
IndisputableMonolith/Foundation/DimensionForcing.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Our calculations reveal a sequence of magnetic phase transitions with thickness, starting from stripe-AFM in monolayer until the bulk like C-AFM is recovered.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
J. G. Bednorz and K. A. M¨ uller, Possible highTc super- conductivity in the Ba−La−Cu−O system, Z. Phys. B - Condensed Matter64, 189 (1986)
work page 1986
-
[2]
S. Jin, T. H. Tiefel, M. McCormack, R. Fastnacht, R. Ramesh, and L. Chen, Thousandfold change in resis- tivity in magnetoresistive La−Ca−Mn−O films, Science 264, 413 (1994)
work page 1994
- [3]
-
[4]
A. Ohtomo and H. Hwang, A high-mobility electron gas at the LaAlO 3/SrTiO3 heterointerface, Nature427, 423 (2004)
work page 2004
-
[5]
D. Lu, D. J. Baek, S. S. Hong, L. F. Kourkoutis, Y. Hikita, and H. Y. Hwang, Synthesis of freestanding single-crystal perovskite films and heterostructures by 8 etching of sacrificial water-soluble layers, Nat. Mater.15, 1255 (2016)
work page 2016
-
[6]
F. M. Chiabrera, S. Yun, Y. Li, R. T. Dahm, H. Zhang, C. K. Kirchert, D. V. Christensen, F. Trier, T. S. Jes- persen, and N. Pryds, Freestanding perovskite oxide films: Synthesis, challenges, and properties, Ann. Phys. 534, 2200084 (2022)
work page 2022
- [7]
-
[8]
S. S. Hong, M. Gu, M. Verma, V. Harbola, B. Y. Wang, D. Lu, A. Vailionis, Y. Hikita, R. Pentcheva, J. M. Rondinelli, and H. Y. Hwang, Extreme tensile strain states in La 0.7Ca0.3MnO3 membranes, Science368, 71 (2020)
work page 2020
-
[9]
D. Ji, S. Cai, T. R. Paudel, H. Sun, C. Zhang, L. Han, Y. Wei, Y. Zang, M. Gu, Y. Zhang, W. Gao, H. Huyan, W. Guo, D. Wu, Z. Gu, E. Y. Tsymbal, P. Wang, Y. Nie, and X. Pan, Freestanding crystalline oxide perovskites down to the monolayer limit, Nature570, 87 (2019)
work page 2019
-
[10]
R. Xu, J. Huang, E. S. Barnard, S. S. Hong, P. Singh, E. K. Wong, T. Jansen, V. Harbola, J. Xiao, B. Y. Wang, S. Crossley, D. Lu, S. Liu, and H. Y. Hwang, Strain- induced room-temperature ferroelectricity in SrTiO 3 membranes, Nat. Commun.11, 3141 (2020)
work page 2020
- [11]
-
[12]
A. Mahajan, D. Johnston, D. Torgeson, and F. Borsa, Magnetic properties of LaVO 3, Phys. Rev. B46, 10966 (1992)
work page 1992
-
[13]
H. C. Nguyen and J. B. Goodenough, Magnetic studies of some orthovanadates, Phys. Rev. B52, 324 (1995)
work page 1995
-
[14]
S. Miyasaka, Y. Okimoto, M. Iwama, and Y. Tokura, Spin-orbital phase diagram of perovskite-typeRVO 3 (R=rare-earth ion or Y), Phys. Rev. B68, 100406 (2003)
work page 2003
-
[15]
Z. Fang and N. Nagaosa, Quantum versus Jahn-Teller orbital physics in YVO 3 and LaVO 3, Phys. Rev. Lett. 93, 176404 (2004)
work page 2004
-
[16]
L. D. Tung, A. Ivanov, J. Schefer, M. R. Lees, G. Bal- akrishnan, and D. M. Paul, Spin, orbital ordering, and magnetic dynamics of LaVO 3: Magnetization, heat ca- pacity, and neutron scattering studies, Phys. Rev. B78, 054416 (2008)
work page 2008
- [17]
-
[18]
G. R. Blake, T. T. M. Palstra, Y. Ren, A. A. Nugroho, and A. A. Menovsky, Transition between orbital order- ings in YVO 3, Phys. Rev. Lett.87, 245501 (2001)
work page 2001
-
[19]
G. R. Blake, T. T. M. Palstra, Y. Ren, A. A. Nu- groho, and A. A. Menovsky, Neutron diffraction, x-ray diffraction, and specific heat studies of orbital ordering in YVO3, Phys. Rev. B65, 174112 (2002)
work page 2002
- [20]
-
[21]
D. A. Mazurenko, A. A. Nugroho, T. T. M. Palstra, and P. H. M. van Loosdrecht, Dynamics of spin and orbital phase transitions in YVO3, Phys. Rev. Lett.101, 245702 (2008)
work page 2008
-
[22]
N. E. Massa, C. Piamonteze, H. C. N. Tolentino, J. A. Alonso, M. J. Mart´ ınez-Lope, and M. T. Casais, Phonon activity and intermediate glassy phase of YVO 3, Phys. Rev. B69, 054111 (2004)
work page 2004
-
[23]
Y. Tao, D. L. Abernathy, T. Chen, T. Yildirim, J. Yan, J. Zhou, J. B. Goodenough, and D. Louca, Lattice and magnetic dynamics in the YVO 3 Mott insulator stud- ied by neutron scattering and first-principles calculations, Phys. Rev. B105, 094412 (2022)
work page 2022
- [24]
-
[25]
N. Wadehra, R. Tomar, R. M. Varma, R. Gopal, Y. Singh, S. Dattagupta, and S. Chakraverty, Planar hall effect and anisotropic magnetoresistance in polar-polar interface of LaVO 3-KTaO3 with strong spin-orbit cou- pling, Nat. Commun.11, 874 (2020)
work page 2020
- [26]
-
[27]
Y. Liu, Z. Liu, M. Zhang, Y. Sun, H. Tian, and Y. Xie, Superconductivity in epitaxially grown LaVO 3/KTaO3 (111) heterostructures, Chin. Phys. B32, 037305 (2023)
work page 2023
- [28]
- [29]
-
[30]
U. L¨ uders, W. Sheets, A. David, W. Prellier, and R. Fr´ esard, Room-temperature magnetism in LaVO3/SrVO3 superlattices by geometrically confined doping, Phys. Rev. B80, 241102 (2009)
work page 2009
-
[31]
P. Hohenberg and W. Kohn, Inhomogeneous electron gas, Phys. Rev.136, B864 (1964)
work page 1964
-
[32]
W. Kohn and L. J. Sham, Self-consistent equations in- cluding exchange and correlation effects, Phys. Rev.140, A1133 (1965)
work page 1965
-
[33]
G. Kresse and J. Hafner, Ab initio molecular dynamics for open-shell transition metals, Phys. Rev. B48, 13115 (1993)
work page 1993
-
[34]
G. Kresse and J. Furthm¨ uller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B54, 11169 (1996)
work page 1996
-
[35]
P. E. Bl¨ ochl, Projector augmented-wave method, Phys. Rev. B50, 17953 (1994)
work page 1994
-
[36]
J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)
work page 1996
-
[37]
S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys, and A. P. Sutton, Electron-energy-loss spec- tra and the structural stability of nickel oxide: An LSDA+U study, Phys. Rev. B57, 1505 (1998)
work page 1998
-
[38]
J. Varignon, N. C. Bristowe, E. Bousquet, and P. Ghosez, Coupling and electrical control of structural, orbital and magnetic orders in perovskites, Sci. Rep.5, 15364 (2015)
work page 2015
-
[39]
Q. Dai, U. Eckern, and U. Schwingenschl¨ ogl, Effects of oxygen vacancies on the electronic structure of the (LaVO3)6/SrVO3 superlattice: a computational study, 9 New J. Phys.20, 073011 (2018)
work page 2018
-
[40]
R. F. L. Evans, W. J. Fan, P. Chureemart, T. A. Ostler, M. O. A. Ellis, and R. W. Chantrell, Atomistic spin model simulations of magnetic nanomaterials, J. Phys.: Condens. Matter26, 103202 (2014)
work page 2014
- [41]
-
[42]
S. Y. Park, A. Kumar, and K. M. Rabe, Charge-order- induced ferroelectricity in LaVO 3/SrVO3 superlattices, Phys. Rev. Lett.118, 087602 (2017)
work page 2017
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.