Magnetic, Structural, and Electronic Properties of CrOCl with the PBE Functional

CrOCl is a van der Waals-layered insulator with an antiferromagnetic ground state, making it a promising platform for exfoliation and the exploration of low-dimensional magnetism. An accurate ab initio description is therefore essential. Previous density-functional studies have shown that DFT+$U$ calculations may erroneously favor ferromagnetic order depending on the choice of parametrization, an issue that cannot be remedied by simply adjusting the value of $U$. Here, we demonstrate that an explicit Hubbard correction is unnecessary: the PBE functional correctly reproduces the AFM ground state while simultaneously improving the description of structural properties. Moreover, PBE provides a reliable account of the electronic structure. These findings clarify the role of correlation effects in CrOCl and identify PBE as a robust starting point for future ab initio studies of CrOCl-based materials.

💡 Research Summary

This paper presents a comprehensive first‑principles investigation of the van‑der‑Waals layered insulator chromium oxochloride (CrOCl), focusing on its magnetic ground state, structural parameters, and electronic structure. CrOCl is experimentally known to adopt an orthorhombic Pmmn lattice at room temperature and to undergo a transition to a monoclinic P2₁/m antiferromagnetic (AFM) superstructure below the Néel temperature (Tₙ ≈ 13.5 K). Earlier density‑functional theory (DFT) studies have commonly employed the DFT + U approach with the Perdew‑Burke‑Ernzerhof (PBE) exchange‑correlation functional. However, many of these works reported a ferromagnetic (FM) ground state when using the Dudarev parametrization of U, in clear disagreement with experiment. Some authors suggested that using an unpolarized PBE functional, or adding a Hund’s J term via the Lichtenstein scheme, could recover the AFM order, but these remedies introduce additional parameters and ambiguity.

The authors set out to test whether the plain, spin‑polarized PBE functional—without any Hubbard U or J corrections—can correctly describe CrOCl. They perform systematic calculations using VASP with PAW potentials, a 650 eV plane‑wave cutoff, and dense k‑point meshes (19 × 21 × 13 for the primitive cell). Both collinear and non‑collinear spin configurations are examined, and spin‑orbit coupling (SOC) is included in selected runs. Van‑der‑Waals interactions are accounted for with the D3(BJ) dispersion correction. For validation, all‑electron LAPW calculations are carried out with Wien2k, and magnetic exchange constants J_ij are extracted via maximally localized Wannier functions and the Green‑function magnetic‑force theorem (TB2J code).

Key findings are as follows:

1. Magnetic Energetics – At the experimental lattice constants, plain PBE predicts the AFM configuration to be lower in energy than FM by ~23 meV per unit cell, in excellent agreement with experiment. LDA yields a similar energy difference (~29 meV). The meta‑GGA SCAN functional also favors AFM, though the margin is smaller (~4 meV). In contrast, the Dudarev‑type PBE + U (U = 3 eV) incorrectly stabilizes FM, while the Lichtenstein‑type PBE + U + J (U = 3 eV, J = 1.5 eV) restores the AFM order. All screened hybrid functionals tested (HSE06, PBE0, B3LYP, etc.) strongly favor FM, indicating that the increased exact‑exchange fraction over‑stabilizes the ferromagnetic alignment. Inclusion of SOC changes the AFM‑FM energy difference by less than 1 meV, confirming that SOC is a secondary effect.

2. Structural Optimization – When the lattice is fully relaxed with a D3 vdW correction, plain PBE reproduces the low‑temperature experimental cell volume within 0.8 % and yields lattice constants a and b within 0.2 % of the measured values; the c‑axis is slightly overestimated by 0.8 %. The monoclinic angle α is modestly exaggerated (≈90.22° versus 90.07° experimentally) but remains realistic. By contrast, PBE + U systematically overestimates the volume (by 2–4 % depending on U) and produces longer Cr–Cr distances, while PBE + U + J improves the angle but still deviates more than plain PBE. LDA underestimates the volume considerably. Thus, the superior structural performance of PBE is not a consequence of a fixed geometry but stems from its balanced treatment of exchange‑correlation and vdW forces.

3. Electronic Structure – The density of states (DOS) computed with PBE shows a band gap of ~1.05 eV, consistent with the limited optical absorption data (onset just below 1.5 eV) measured at room temperature. PBE + U and PBE + U + J enlarge the gap to 2.2–2.5 eV, while SCAN yields a gap comparable to PBE. Hybrid functionals predict much larger gaps (≈3.5–4 eV) and display strong hybridization between Cr d‑states and Cl/O p‑states, which correlates with their tendency to favor FM. Projected DOS analysis reveals that in the AFM ground state the occupied t₂g manifold (dₓ²₋ᵧ², d_zx, d_yz) is fully spin‑up polarized, corresponding to a high‑spin d³ configuration, while the e_g orbitals remain empty. This simple crystal‑field picture explains why strong on‑site correlation corrections are unnecessary for CrOCl.

4. Magnetic Exchange – Wannier‑based extraction of exchange constants confirms that the dominant nearest‑neighbor Cr–Cr interaction is antiferromagnetic, consistent with the experimentally observed Néel temperature. The magnitude of J extracted from PBE matches the energy difference between AFM and FM configurations, further supporting the adequacy of PBE for describing magnetic interactions in this material.

Overall, the study demonstrates that CrOCl is a rare example of a layered magnetic insulator whose essential physics—magnetic ordering, lattice geometry, and insulating gap—are captured accurately by the semi‑local PBE functional when combined with a modest vdW correction. The need for adjustable Hubbard U or hybrid exchange is eliminated, simplifying future computational explorations of CrOCl‑based heterostructures, strain engineering, and low‑dimensional magnetism. The authors recommend using plain PBE + D3 as the default first‑principles approach for CrOCl and related transition‑metal oxyhalides, reserving more expensive many‑body techniques (GW, DMFT) for quantitative refinements of the band gap or excited‑state properties.