pith. sign in

arxiv: 2605.22415 · v1 · pith:WYGPLE75new · submitted 2026-05-21 · ❄️ cond-mat.str-el · cond-mat.mes-hall· cond-mat.mtrl-sci

Emergent magnetic and charge ordered phases in freestanding ultrathin ce{LaVO3}

Pith reviewed 2026-05-22 03:45 UTC · model grok-4.3

classification ❄️ cond-mat.str-el cond-mat.mes-hallcond-mat.mtrl-sci
keywords LaVO3ultrathin filmsmagnetic phase transitionscharge orderingdensity functional theorypolar catastrophefreestanding membranessurface doping
0
0 comments X

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.

This paper uses density functional theory to study freestanding ultrathin LaVO3 films and finds that their magnetic ground state changes systematically with thickness. Monolayer films adopt a stripe antiferromagnetic order that evolves toward the bulk C-type antiferromagnetic state as layers are added. In films thicker than four layers, a polar catastrophe drives charge transfer that dopes the surface layers, producing stripe antiferromagnetic and ferromagnetic surface states while the central layers retain bulk-like order. The same calculations on charge-doped monolayers show that hole doping stabilizes ferromagnetism and induces charge ordering whose periodicity depends on the doping level.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2605.22415 by Ashutosh Anand, Mukul Kabir.

Figure 1
Figure 1. Figure 1: (a) √ 2× √ 2 supercell of the orthorhombic structure. (b) Phase diagram of the J1−J2 Heisenberg model, depending on the frustration angle θ = tan−1 (J2/J1) the ground state can be FM, Neel-AFM or stripe-AFM. The gray region in the phase diagram corresponds to the spin-liquid regime. a Hubbard-like on-site Coulomb repulsion U, using the rotationally invariant scheme [37]. An effective U of 3 eV is considere… view at source ↗
Figure 2
Figure 2. Figure 2: (a) Representative crystal structure to illustrate [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: (a) Schematic evolution of magnetic ground state in LaVO [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: (a) Layer resolved total vanadium charge across [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (a) Stripe and (b) Checkerboard charge order [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Relative energies of different charge ordering pat [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Relative energies of different charge ordering pat [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

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)
  1. [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.
  2. [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)
  1. [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.
  2. [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

2 responses · 0 unresolved

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
  1. 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

  2. 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

0 steps flagged

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

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the adequacy of DFT for this strongly correlated system and on the applicability of the polar-catastrophe picture to freestanding films; no explicit free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption Density functional theory with standard approximations suffices to determine the magnetic and charge-ordered ground states of LaVO3.
    Invoked implicitly throughout the computational results; standard but unverified for freestanding ultrathin films.

pith-pipeline@v0.9.0 · 5742 in / 1149 out tokens · 44773 ms · 2026-05-22T03:45:40.222350+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

42 extracted references · 42 canonical work pages

  1. [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)

  2. [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)

  3. [3]

    Fiebig, T

    M. Fiebig, T. Lottermoser, D. Meier, and M. Trassin, The evolution of multiferroics, Nat. Rev. Mater.1, 1 (2016)

  4. [4]

    Ohtomo and H

    A. Ohtomo and H. Hwang, A high-mobility electron gas at the LaAlO 3/SrTiO3 heterointerface, Nature427, 423 (2004)

  5. [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)

  6. [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)

  7. [7]

    Xiang, P

    L. Xiang, P. Cai, Q. Li, Q. Shi, T. Miao, Y. Bai, Y. Yu, F. Lan, S. Guo, G. Chen, W. Wang, L. Yin, Y. Zhang, and J. Shen, Boosting ferromagnetism in free- standing electronically phase separated manganite thin films, Phys. Rev. Mater.8, 054417 (2024)

  8. [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)

  9. [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)

  10. [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)

  11. [11]

    Bordet, C

    P. Bordet, C. Chaillout, M. Marezio, Q. Huang, A. San- toro, S. Cheong, H. Takagi, C. Oglesby, and B. Batlogg, Structural aspects of the crystallographic-magnetic tran- sition in LaVO3 around 140 K, J. Solid State Chem.106, 253 (1993)

  12. [12]

    Mahajan, D

    A. Mahajan, D. Johnston, D. Torgeson, and F. Borsa, Magnetic properties of LaVO 3, Phys. Rev. B46, 10966 (1992)

  13. [13]

    H. C. Nguyen and J. B. Goodenough, Magnetic studies of some orthovanadates, Phys. Rev. B52, 324 (1995)

  14. [14]

    Miyasaka, Y

    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)

  15. [15]

    Fang and N

    Z. Fang and N. Nagaosa, Quantum versus Jahn-Teller orbital physics in YVO 3 and LaVO 3, Phys. Rev. Lett. 93, 176404 (2004)

  16. [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)

  17. [17]

    Kawano, H

    H. Kawano, H. Yoshizawa, and Y. Ueda, Magnetic be- havior of a mott-insulator YVO 3, J. Phys. Soc. Jpn.63, 2857 (1994)

  18. [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)

  19. [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)

  20. [20]

    Ulrich, G

    C. Ulrich, G. Khaliullin, J. Sirker, M. Reehuis, M. Ohl, S. Miyasaka, Y. Tokura, and B. Keimer, Magnetic neu- tron scattering study of YVO 3: Evidence for an orbital peierls state, Phys. Rev. Lett.91, 257202 (2003)

  21. [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)

  22. [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)

  23. [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)

  24. [24]

    Hotta, T

    Y. Hotta, T. Susaki, and H. Hwang, Polar discontinuity doping of the LaVO 3/SrTiO3 interface, Phys. Rev. Lett. 99, 236805 (2007)

  25. [25]

    Wadehra, R

    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)

  26. [26]

    Halder, M

    S. Halder, M. Garg, N. S. Mehta, A. Kumari, R. Sharma, T. Das, S. Chakraverty, and G. Sheet, Unconventional su- perconductivity at LaVO3/SrTiO3 interfaces, ACS Appl. Electron. Mater.4, 5859 (2022)

  27. [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)

  28. [28]

    Tomar, S

    R. Tomar, S. Kakkar, C. Bera, and S. Chakraverty, Anisotropic magnetoresistance and planar Hall effect in (001) and (111) LaVO 3/SrTiO3 heterostructures, Phys. Rev. B103, 115407 (2021)

  29. [29]

    Sheets, B

    W. Sheets, B. Mercey, and W. Prellier, Effect of charge modulation in (LaVO 3)m(SrVO3)n superlattices on the insulator-metal transition, Applied Physics Letters91 (2007)

  30. [30]

    L¨ uders, W

    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)

  31. [31]

    Hohenberg and W

    P. Hohenberg and W. Kohn, Inhomogeneous electron gas, Phys. Rev.136, B864 (1964)

  32. [32]

    Kohn and L

    W. Kohn and L. J. Sham, Self-consistent equations in- cluding exchange and correlation effects, Phys. Rev.140, A1133 (1965)

  33. [33]

    Kresse and J

    G. Kresse and J. Hafner, Ab initio molecular dynamics for open-shell transition metals, Phys. Rev. B48, 13115 (1993)

  34. [34]

    Kresse and J

    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)

  35. [35]

    P. E. Bl¨ ochl, Projector augmented-wave method, Phys. Rev. B50, 17953 (1994)

  36. [36]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)

  37. [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)

  38. [38]

    Varignon, N

    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)

  39. [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)

  40. [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)

  41. [41]

    Arima, Y

    T. Arima, Y. Tokura, and J. B. Torrance, Variation of op- tical gaps in perovskite-type 3d transition-metal oxides, Phys. Rev. B48, 17006 (1993)

  42. [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)