Effective theory of quantum phases in the dipolar planar rotor chain

Abstract: In this work, we develop a theoretical description of the collective behavior of interacting dipolar planar rotors by using time independent perturbation theory and a small angle quadratic approximation. The ground state properties for both the ordered and disordered quantum phases of the system are directly calculated and analyzed. Time-independent perturbation theory is shown to be appropriate for the disordered phase. For the ordered phase, we construct a quadratic approximation based on the stable equilibrium configurations of the dipolar ordering; we show that the inclusion of the quartic terms from the expansion of the potential energy are essential to correct the shift in the energy spectrum due to quantization ambiguities. Numerical techniques such as Exact Diagonalization and Density Matrix Renormalization Group are used for the benchmark the quality of both approximations.

• The paper develops an analytical framework capturing quantum phase transitions by balancing rotational kinetic energy and dipolar interactions.
• It employs perturbative analysis and harmonic approximations to derive closed-form expressions for ground state energy, angular momentum variance, and polarization.
• The results offer actionable insights for designing quantum simulations and experiments with dipolar molecular chains, confirmed by numerical methods.

Effective Theory of Quantum Phases in the Dipolar Planar Rotor Chain

Introduction and Motivation

The study develops a controlled analytical framework for understanding the quantum phases realized in chains of interacting dipolar planar rotors. These systems, where each constituent possesses a rotational degree of freedom and interacts via nearest-neighbor dipole-dipole couplings, have direct implications in condensed matter physics, quantum simulation, and quantum computation with molecular arrays. The competition between rotational kinetic energy (parameterized by the rotational constant $B$) and the dipolar interaction (scaling as $\mu^2$) generates a quantum phase transition between disordered and ordered ferroelectric phases.

The universal relationship between $B$ and $\mu$ controlling the system's phase structure is encapsulated in a phase diagram (Figure 1), which demarcates the ordered, disordered, and critical regions for experimentally relevant endofullerene chain realizations.

Figure 1: Phase diagram for endofullerene chains showing ordered (blue), disordered (red), and critical (white) phases as a function of molecular rotational constant $B$ and dipole moment $\mu$. The black curve indicates critical boundary.

Model Formulation

The Hamiltonian describes $N$ planar rotors coupled by nearest-neighbor dipole-dipole interactions and possessing kinetic energy from free rotations; the total Hamiltonian is

$$H = \sum_{i} L_i^2 + g \sum_{\langle ij \rangle} \left[ \sin \varphi_i \sin \varphi_j - 2 \cos \varphi_i \cos \varphi_j \right]$$

where $g \propto \mu^2 / r^3$ encapsulates the effective interaction strength. The model features a $Z_2$ symmetry (invariance under $\mu^2$0 for all $\mu^2$1), leading to a quantum phase transition at a critical $\mu^2$2.

The spatial structure is a linear co-planar chain (Figure 2).

Figure 2: The $\mu^2$3 planar rotor chain arranged linearly in a co-planar geometry.

Analytical Treatments

Distinct analytical frameworks are constructed for each phase.

Disordered Phase ($\mu^2$4)

Treating $\mu^2$5 as a perturbation, ground state properties are computed to second-order via standard time-independent perturbation theory in the angular momentum basis. For the ground state energy, variance of total angular momentum, polarization, and correlation functions, closed forms are derived, e.g.

$\mu^2$6

Polarization vanishes in this regime. This approach is quantitatively validated via exact diagonalization for small $\mu^2$7 (Figure 3).

Figure 3: Ground state energy and variance of total angular momentum for $\mu^2$8 planar rotors versus $\mu^2$9: agreement between ED and analytical results is manifest in both phases.

Ordered Phase ($B$0)

The system is expanded around one of the two symmetry-broken ordered configurations ($B$1 or $B$2), and a quadratic (harmonic) approximation is used for small fluctuations. After Fourier transformation into normal modes, the Hamiltonian becomes a sum of decoupled QHOs (normal mode phonons). Quantization is implemented via second quantization, yielding the spectrum and observables in terms of the normal mode frequencies.

The ground state energy is, up to leading order and in the thermodynamic limit,

$B$3

where $B$4 is a complete elliptic integral. Corrections to the spectrum from quartic (and higher) terms are also analyzed, showing that only the quartic term provides a finite correction in the large-$B$5 limit: a constant shift of $B$6 per degree of freedom, which is required for agreement with exact diagonalization (Figure 4).

Figure 4: Absolute errors between ED and analytical results reveal a spectrum shift for $B$7 traceable to quantization ambiguities on compact manifolds.

Systematic comparison with DMRG for large $B$8 further verifies the harmonic theory, including its limitations in the critical domain (Figures 5, 6, and 7). As illustrated in Figure 5, finite-size convergence and the analytic thermodynamic limit are established for both strongly ordered and disordered regimes.

Figure 5: Chemical potential as function of chain size $B$9 for various $\mu$0; analytic predictions closely track DMRG except near criticality where correlations grow.

Figure 6: Analytical vs. DMRG for several observables (energy, angular momentum variance, polarization, orientational correlation) as functions of $\mu$1 for $\mu$2.

Figure 7: Relative discrepancies between analytic and DMRG calculations across all observables, outside the critical regime, are minimal.

Quantization Ambiguity and Higher-Order Corrections

A key theoretical insight is the necessity to consistently handle quantization on compact manifolds ($\mu$3) rather than $\mu$4. Neglecting this affects spectrum and observable calculations in the ordered phase. This is resolved by incorporating quartic terms in the potential expansion, verified explicitly in the exactly-solvable two-rotor case via mapping to Mathieu equations. The shift is demonstrated to be a universal correction per mode in the large-$\mu$5 regime.

Implications and Outlook

This analytic framework yields quantitatively robust predictions for ground state energy, angular momentum, polarization, and correlation functions across the full range of physical parameters, excepting a narrow critical window near the quantum phase transition. The methods presented circumvent limitations of numerical approaches (ED, DMRG, PIMC, neural-network wavefunctions) regarding system size and computational cost.

The results directly inform the phase structure of realistic molecular systems—e.g., endofullerene and water-in-nanotube chains—enabling rapid evaluation of phase diagrams across molecular parameter space. For quantum simulation contexts, analytic access to excitation spectra and correlations can guide design and interpretation in quantum device architectures based on molecular arrays.

A natural extension is the rigorous global development of a $\mu$6 field theory to describe the ordered phase, which could account for fluctuation-mediated corrections and provide analytic access even closer to the critical regime. The presented techniques are also adaptable for computing observables such as entanglement entropy and linear response functions.

Conclusion

A controlled and detailed analytic treatment for the dipolar planar rotor chain problem is achieved, with explicit capture of quantum phase structure, numerical accuracy confirmed against state-of-the-art methods, and illumination of foundational aspects such as quantization on manifolds. The approach provides actionable insight for both theoretical studies of quantum phase transitions in low-dimensional rotor systems and practical guidance for experimental platforms exploiting collective rotational degrees of freedom.