Dielectric response of modified Hubbard models with neutral-ionic and   Peierls transitions
  The dipole P(F) of systems with periodic boundary conditions (PBC) in a static electric field F is applied to one-dimensional Peierls-Hubbard models for organic charge-transfer (CT) salts. Exact results for P(F) are obtained for finite systems of N = 14 and 16 sites that are almost converged to infinite chains in deformable lattices subject to a Peierls transition. The electronic polarizability per site, \alpha_{el} = (\partial P/\partial F)_0, of rigid stacks with alternating transfer integrals t(1 +/- \delta) diverges at the neutral-ionic transition for \delta = 0 but remains finite for \delta > 0 in dimerized chains. The Peierls or dimerization mode couples to charge fluctuations along the stack and results in large vibrational contributions, \alpha_{vib}, that are related to \partial P/\partial \delta and that peak sharply at the Peierls transition. The extension of P(F) to correlated electronic states yields the dielectric response \kappa of models with neutral-ionic or Peierls transitions, where \kappa peaks >100 are found with parameters used previously for variable ionicity \rho and vibrational spectra of CT salts. The calculated \kappa accounts for the dielectric response of CT salts based on substituted TTFs (tetrathiafulvalene) and substituted CAs (chloranil). The role of lattice stiffness appears clearly in models: soft systems have a Peierls instability at small \rho and continuous crossover to large \rho, while stiff stacks such as TTF-CA have a first-order transition with discontinuous \rho that is both a neutral-ionic and Peierls transition. The transitions are associated with tuning the electronic ground state of insulators via temperature or pressure in experiments, or via model parameters in calculations. 
Zoltan G. Soos (1), Sharon A. Bewick (1), Andrea Peri (2), Anna   Painelli (2) ((1) Department of Chemistry, Princeton Unive...