Symmetry crossover and excitation thresholds at the neutral-ionic transition of the modified Hubbard model

Abstract

Exact ground states, charge densities, and excitation energies are found using valence bond methods for N-site modified Hubbard models with uniform spacing. At the neutral-ionic transition (NIT), the ground state has a symmetry crossover in $4n$ and $4n+2$ rings with periodic and antiperiodic boundary conditions, respectively. Large site energies $\ensuremath{\Delta}$ stabilize a paired state of the half-filled chain, while large U stabilizes a covalent state. Finite-transfer integrals t shift the NIT to the covalent side of $U\ensuremath{-}2\ensuremath{\Delta}.$ Exact results to $N=16$ in the full basis and to $N=22$ in a restricted basis for large U and $\ensuremath{\Delta}$ are extrapolated to obtain the crossover and charge density of extended chains. The modified Hubbard model has a continuous NIT between a diamagnetic band insulator on the paired side and a paramagnetic Mott insulator on the covalent side. The singlet-triplet (ST), singlet-singlet (SS), and charge gaps for finite N indicate that the ST and SS gaps close at the NIT with increasing U, and that the charge gap vanishes only there. Finite-N excitations constrain all singularities to $\ifmmode\pm\else\textpm\fi{}0.1t$ of the symmetry crossover. The NIT is interpreted as a localized ground state (GS) with finite gaps on the paired side and an extended GS with vanishing ST and SS gaps on the covalent side. The charge gap and charge stiffness indicate a metallic GS at the transition that, however, is unconditionally unstable to dimerization. Finite $\ensuremath{\Delta}$ breaks electron-hole $(e\ensuremath{-}h)$ symmetry, but the modified Hubbard model has an extended $e\ensuremath{-}h$ symmetry, and a strong mixing of spin and charge excitations is limited to a few t's about the NIT. Exact finite-size results complement other approaches to valence or ferroelectric transitions in organic charge-transfer salts or in inorganic oxides, and to electron-vibration coupling and structural instabilities in one-dimensional systems.