Surface lone-pair polarization probed by quantum-geometric transport in tellurium

Technical Summary: Surface Lone-Pair Polarization Probed by Quantum-Geometric Transport in Tellurium

Problem Statement
Stereochemically active lone pairs are ubiquitous microscopic sources of polarity in molecules and solids, yet their collective behavior in crystals is often obscured by symmetry or confined to surfaces. In trigonal tellurium (Te), a paradigmatic example, the bulk crystal structure consists of helical chains where interchain interactions are mediated by lone pairs. While the bulk possesses $C_3$ rotational symmetry that enforces the exact cancellation of in-plane dipoles (resulting in zero net polarization), the loss of coordination at surfaces in films or flakes breaks this equivalence. This allows for the emergence of an in-plane surface polarization. Recent transport experiments in Te flakes revealed a giant anisotropy in linear resistivity ($\rho_{xx}$) and a large longitudinal second-harmonic response ($\chi_{xxx}$) with distinct angular and gate-voltage dependencies. The fundamental question addressed is identifying the microscopic degree of freedom probed by these linear and nonlinear transport signals, particularly given that standard Berry-curvature dipole mechanisms typically predict Hall-like responses rather than the observed dominant longitudinal effects.

Methodology
The authors employ a multi-faceted approach combining semiclassical transport theory, quantum-geometric formalism, and first-principles calculations:

1. Theoretical Framework: The study utilizes a semi-classical wavepacket dynamics approach within Boltzmann transport theory. The authors decompose the crystal potential $V(r)$ into inversion-even ($V_{even}$) and inversion-odd ($V_{odd}$) components. The inversion-odd component, $V_{odd}$, models the surface lone-pair polarization field $E_0$.
2. Hamiltonian Construction: A multi-band $k \cdot p$ Hamiltonian is constructed near the $H$ point of the Brillouin zone. Using Löwdin downfolding, the authors project the full Hamiltonian onto the two uppermost valence bands ($H_4$ and $H_5$). Crucially, they include virtual coupling to remote bands ($H_6, H_6'$) to generate the necessary quantum-geometric structure (Christoffel symbols) in the low-energy manifold.
3. Quantum-Geometric Corrections: The inversion-odd potential induces a gauge-invariant positional shift of the wavepacket center of mass and a second-order correction to the wavepacket energy. These effects generate quantum-geometric (QG) corrections to the wavepacket velocity, which are linear and quadratic in the effective polar field $E_0$.
4. Transport Calculations: The authors derive expressions for linear and nonlinear conductivities. They relate the measured harmonic voltages (first and second harmonic) to the QG corrections, specifically identifying how the linear response probes the second moment of the polarization field ($\langle E_0^2 \rangle$) and the nonlinear response probes the first moment ($\langle E_0 \rangle$).
5. Microscopic Modeling: A coarse-grained lattice model is developed to describe the thermodynamics of the surface lone-pair polarization, treating the local dipole weights as variables subject to surface-induced symmetry breaking.
6. First-Principles Validation: Density Functional Theory (DFT) calculations are performed on Te slabs to visualize the local polar slopes (gradients of the Hartree potential) and verify the scaling of the surface polarization with flake thickness.

Key Contributions and Results

• Identification of the Mechanism: The paper demonstrates that the observed transport anomalies in Te are driven by a net surface lone-pair polarization texture ($E_0$). This texture acts as an inversion-odd component of the crystal potential, producing QG corrections to the semi-classical velocity.
• Distinct Transport Signatures:
  • Linear Response: The surface polarization generates an inversion-even geometric correction to the linear conductivity, $\delta\sigma_{xx}^{geo} \propto E_0^2$. This correction governs the first-harmonic resistance, reproducing the experimentally observed $\sin^2\theta$ angular dependence.
  • Nonlinear Response: The polarization generates a dominant longitudinal nonlinear conductivity, $\chi_{xxx}^{geo} \propto E_0$. This term dictates the second-harmonic response, reproducing the observed $\sin^3\theta$ angular dependence. The sign of $\chi_{xxx}$ depends on the chirality of the Te sample (related to the direction of $E_0$).
• Scaling Laws: The theory predicts specific scaling behaviors with respect to the chemical potential (controlled by gate voltage). The linear geometric correction scales as $\mu^{-5/2}$, and the nonlinear response scales as $\mu^{-7/2}$. These predictions show close agreement with experimental gating data.
• Thickness Dependence: The model predicts that the average surface polarization scales as $1/d$ (where $d$ is the flake thickness) due to the exponential decay of the polarization field into the bulk. This is confirmed by DFT calculations and experimental data showing a $1/d$ dependence in the measured nonlinear signal.
• Thermodynamic Connection: The paper establishes a direct link between transport signals and the thermodynamics of the surface polarization. The second-harmonic voltage tracks the first thermal moment (net polarization $\langle P_x \rangle_T$), while the linear anisotropy tracks the second thermal moment (fluctuating polar strength $\langle P_x^2 \rangle_T$).

Significance and Claims
The paper claims to provide a microscopic route to understanding and engineering polarization-driven, quantum-geometric electronic devices based on tellurium allotropes. By demonstrating that rectified voltages in Te flakes are directly proportional to the surface lone-pair polarization, the work suggests a mechanism for designing controllable, radio-frequency quantum-geometric rectifiers. The authors emphasize that this framework extends beyond the interpretation of existing data, predicting the existence of higher-order nonlinear responses and offering a general approach for engineering QG rectifiers in materials with frustrated polar structures or surface-confined dipole textures. The work highlights that quantum-geometric transport serves as a dual probe for both the fluctuating polar strength and the net inversion-odd surface order.