OS-Nav: Quasi-1D Transition Metal Oxide

• OS-Nav is a transition metal oxide that exhibits pronounced quasi-1D behavior in electronic transport and magnetism due to distinct orbital ordering of symmetric and antisymmetric t₂g states.
• Its low-energy model, based on a three-orbital Hubbard Hamiltonian, reveals specific crystal field splittings and significant anisotropy in hopping amplitudes along chain versus transverse directions.
• The magnetic exchange interactions arise from a competition between double exchange and superexchange, explaining the strong ferromagnetic coupling along chains and complex interchain dynamics.

OS-Nav (Orthorhombic Sodium Vanadate, NaV₂O₄) is a transition metal oxide characterized by its double chains of edge-sharing VO₆ octahedra and pronounced quasi-one-dimensional (quasi-1D) behavior in both electronic transport and magnetism. Despite a crystallographically three-dimensional arrangement, NaV₂O₄ exhibits a strong anisotropy in its conductivity and magnetic exchange interactions, arising primarily from an orbital ordering pattern allowed by its orthorhombic $Pnma$ symmetry. Central to this behavior is the division of the V $t_{2g}$ states into symmetric (delocalized, metallic) and antisymmetric (localized) combinations under the mirror reflection $y \rightarrow -y$ (Pchelkina et al., 2012).

1. Low-Energy Hamiltonian and Crystal Field Splittings

The electronic structure of NaV₂O₄ is described by a three-orbital (t₂g) Hubbard Hamiltonian in the basis of local crystal-field eigenstates at each vanadium site: $$H = \sum_{\langle ij\rangle}\sum_{m m'} t_{ij}^{m m'}\,c_{i m}^\dagger c_{j m'} + \sum_{i}\sum_{m}\Delta_{m} n_{i m} + H_{\rm int}$$ where $c_{i m}^\dagger$ is the creation operator for an electron in the $m$th t₂g Wannier orbital at site $i$, $\Delta_m$ are crystal-field splittings, $t_{ij}^{m m'}$ are hopping integrals, and $H_{\rm int}$ denotes on-site Coulomb interactions parameterized by $\mathcal{U}\approx3.15$–$3.2$ eV and Hund’s coupling $\mathcal{J}\approx0.63$ eV.

Distinct crystal-field splitting patterns are observed at the two non-equivalent vanadium sites (V1, V2), with the energetically lowest state for both being symmetric (S1), followed by an antisymmetric state (A), and then a higher-energy symmetric state (S2): $\begin{array}{c|ccc} & \Delta_{\rm S1} & \Delta_{\rm A} & \Delta_{\rm S2} \ \hline \mathrm{V1} & 0 & 46 & 204 \ \mathrm{V2} & 0 & 10 & 264 \ \end{array}$ Typical nearest-neighbor hopping amplitudes for the $(\mathrm{S1}, \mathrm{A}, \mathrm{S2})$ basis are:

• Along the chains: $t_{1,1'}^{\rm S1,S1}\simeq130$ meV, $t_{1,1'}^{\rm A,A}\simeq-18$ meV, $t_{1,1'}^{\rm S2,S2}\simeq66$ meV.
• Between chains: $t_{1,5}^{\rm S1,S1}\simeq7.7$ meV, $t_{1,8}^{\rm S1,S1}\simeq10.8$ meV, $t_{1,2}^{\rm S1,S2}\simeq2.0$ meV.

2. Orbital Ordering and Its Symmetry Consequences

The $Pnma$ symmetry of NaV₂O₄ enforces a mirror operation $y\rightarrow-y$ at each V site. Diagonalization of the local crystal field yields three t₂g orbitals per site: $\phi_{S1}^i,\;\phi_{A}^i,\;\phi_{S2}^i$ with $\phi_A^i$ being antisymmetric (localized) and $\phi_{S1}^i, \phi_{S2}^i$ symmetric (delocalized) under the mirror operation.

At V1, explicit real-harmonic composition is: $$\begin{aligned} \phi_{A}^1 &\simeq 0.94\,|xy\rangle - 0.33\,|yz\rangle \ \phi_{S1}^1 &\simeq 0.84\,|3z^2-r^2\rangle + 0.19\,|zx\rangle + 0.51\,|x^2-y^2\rangle \ \end{aligned}$$ The antisymmetric $\phi_A$ state carries one localized electron per V, whereas the symmetric orbitals are partially filled, forming an itinerant, quasi-1D conduction band.

3. Band Structure and Directional Dispersion

Restricting to the low-energy, symmetric sector (orbital ordering $\cal O1$), the effective model can be cast as a two-band tight-binding (TB) Hamiltonian for $\phi_{S1}$ states on the V1 and V2 sublattices. The dominant hopping is along the $b$-axis (chain direction), with $t_b\approx130$ meV, while transverse hoppings are much weaker ($t_\perp\approx8$–$11$ meV).

The band Hamiltonian: $$H_{\rm band}(k_b,k_\perp) = \begin{pmatrix} \varepsilon_S + 2t_b\cos k_b & 2t_\perp\cos\frac{k_\perp}{2} \ 2t_\perp\cos\frac{k_\perp}{2} & \varepsilon_S + 2t_b\cos k_b \ \end{pmatrix}$$ with band energies

$$E_{\pm}(k_b,k_\perp) = \varepsilon_S + 2t_b\cos k_b \pm 2t_\perp\cos\frac{k_\perp}{2}$$

The bandwidth along $b$ ($4|t_b|\sim0.5$ eV) is much larger than that perpendicular to $b$ ($\sim0.04$ eV), manifesting quasi-1D metallicity.

4. Magnetic Exchange Interactions and Double Exchange

Magnetic couplings are described by a Heisenberg model: $$\hat{\cal H}_S = -\sum_{i>j}J_{ij}\,\mathbf{e}_i \cdot \mathbf{e}_j$$ with anisotropic exchange constants. Calculated NN exchange values (in meV): $\begin{array}{c|cc} \mathrm{bond} & J_{ij}^{\mathrm{along}\;b} & J_{ij}^{\mathrm{perp}} \ \hline \text{1--1'} & 30.1\ (\mathrm{V1}\!-\!\mathrm{V1}) & 7.7\ (\mathrm{V1}\!-\!\mathrm{V2}) \ \text{5--5'} & 30.0\ (\mathrm{V2}\!-\!\mathrm{V2}) & 3.8\ (\mathrm{V2}\!-\!\mathrm{V2}) \ \end{array}$ All NN couplings are ferromagnetic, with maxima along the chain. The total exchange decomposes into double-exchange (DE) and superexchange (SE) contributions: $J_{ij} = J_{ij}^{\rm DE} + J_{ij}^{\rm SE}$ For bond 1–1', $J^{\rm DE}_{1,1'}\simeq45.5$ meV, $J^{\rm SE}_{1,1'}\simeq-18.8$ meV, reflecting the half-metallic and correlated nature. Correlations that decrease the spin splitting $\Delta_{\rm ex}$ enhance $|J^{\rm SE}|$, promoting antiferromagnetic (AFM) interchain couplings and possibly stabilizing helical or $\uparrow\downarrow\uparrow\downarrow$ spin arrangements as observed experimentally.

5. Electrical Resistivity Anisotropy and Transport

In the relaxation-time approximation, the conductivity tensor $\sigma_{\alpha\beta}$ is expressed as: $$\sigma_{\alpha\beta} = \frac{e^2\tau}{V}\sum_l\int\frac{d\mathbf{k}}{(2\pi)^3} v_l^\alpha(\mathbf{k})\,v_l^\beta(\mathbf{k})\, \left(-\frac{\partial f}{\partial\epsilon}\right)_{\epsilon_l(\mathbf{k})}$$ with $v_l^\alpha = \partial\epsilon_l/\partial k_\alpha$. For a constant relaxation time,

$$\frac{\rho_\perp}{\rho_\parallel} \simeq \frac{\langle v_x^2\rangle}{\langle v_b^2\rangle} \approx 20\text{–}35$$

In the FM half-metallic state with $\cal O1$: $$\frac{\rho_{xx}}{\rho_{yy}} \simeq 34,\quad \frac{\rho_{zz}}{\rho_{yy}} \simeq 19$$ These values are consistent with experimental $\rho_\perp/\rho_\parallel > 20$. Antiferromagnetic ordering increases the resistivity anisotropy (up to $\sim210$ in one transverse direction) but the primary cause remains the orientation and delocalization of the symmetric t₂g orbitals.

6. Correlation Effects and Experimental Consequences

Although the three-dimensional V–V network might suggest isotropic behavior, the symmetry-allowed orbital ordering into one antisymmetric (localized) and two symmetric (delocalized) t₂g states crucially determines the electronic and magnetic structure (Pchelkina et al., 2012). The result is a single, strongly dispersive half-metallic band along the chain axis, imparting both the quasi-1D character of transport and pronounced anisotropy in magnetic exchanges.

Theoretical calculations align well with transport and magnetic measurements on NaV₂O₄, providing compelling evidence that the observed quasi-1D conducting and magnetic properties derive from emergent electronic structure rather than crystal geometry alone. Correlation effects shift the balance between double exchange and superexchange, and are instrumental in explaining experimentally observed AFM phases and modulated spin patterns. A plausible implication is that the interplay of crystal symmetry, orbital ordering, and electronic correlations in OS-Nav serves as a model for other transition-metal oxides exhibiting emergent low-dimensional physics.