We present a pseudospin formulation for the post-quench dynamics of charge-density-wave (CDW) order in the half-filled spinless Holstein model on a square lattice, assuming spatially homogeneous evolution. This Anderson pseudospin description captures the coherent nonequilibrium dynamics of the coupled electron-lattice system. Numerical simulations reveal three distinct dynamical regimes of the CDW order parameter following a quench-locked oscillations, Landau-damped dynamics, and overdamped relaxation-closely paralleling quench dynamics in BCS superconductors and other electronically driven symmetry-breaking phases. Crucially, however, the presence of dynamical lattice degrees of freedom leads to qualitatively different long-time behavior. In particular, while the oscillation amplitude is reduced in the damped regimes, CDW oscillations do not fully decay but instead persist indefinitely due to feedback from the lattice field. We further show that these persistent oscillations are characterized by a nonequilibrium electronic distribution, which provides an intuitive understanding of both their amplitude and the renormalization of the oscillation frequency relative to the bare Holstein phonon frequency. Our results highlight the essential role of lattice dynamics in nonequilibrium ordered phases and establish a clear distinction between electron-lattice-driven CDW dynamics and their purely electronic counterparts.
We study postquench dynamics of charge density wave (CDW) order in the square-lattice $t\text{\ensuremath{-}}V$ model. The ground state of this system at half-filling is characterized by a checkerboard modulation of particle density. A generalized self-consistent mean-field method, based on the time-dependent variational principle, is employed to describe the dynamical evolution of the CDW states. Assuming a homogeneous CDW order throughout the quench process, the time-dependent mean-field approach is reduced to the Anderson pseudospin method. Quench simulations based on the Bloch equation for pseudospins produce three canonical behaviors of order-parameter dynamics: phase-locked persistent oscillation, Landau-damped oscillation, and dynamical vanishing of the CDW order. We further develop an efficient real-space von Neumann equation method to incorporate dynamical inhomogeneity into simulations of quantum quenches. Our large-scale simulations uncover complex pattern formations in the postquench CDW states, especially in the strong quench regime. The emergent spatial textures are characterized by super density modulations on top of the short-period checkerboard CDW order. Our demonstration of pattern formation in quenched CDW states, described by a simple broken
Z2 symmetry, underscores the importance of dynamical inhomogeneity in quantum quenches of many-body systems with more complex orders.
Charge-density wave phases in quantum materials stem from the complex interplay of electronic and lattice degrees of freedom. Nowadays, various time-resolved spectroscopy techniques allow to actively manipulate such phases and monitor their dynamics in real time. Modeling such nonequilibrium dynamics theoretically is a great challenge and exact methods can usually only treat a small number of atoms and finitely many phonons. We approach the melting of charge-density waves in a Holstein model after a sudden switch-on of the electronic hopping from two perspectives. We prove that in the noninteracting and in the strong-coupling limit, the CDW order parameter on high-dimensional hypercubic lattices obeys a factorization relation for long times, such that its dynamics can be reduced to the one-dimensional case. Secondly, we present numerical results from semiclassical techniques based on the truncated Wigner approximation for two spatial dimensions. A comparison with exact data obtained for a Holstein chain shows that a semiclassical treatment of both the electrons and phonons is required in order to correctly describe the phononic dynamics. This is confirmed, in addition, for a quench in the electron-phonon coupling strength.
We investigate the nonequilibrium dynamics induced by an interaction quench in the semiclassical Holstein model within the Ehrenfest nonadiabatic framework, which describes an isolated hybrid quantum-classical system with strictly conserved total energy. Focusing on the half-filled case, where the equilibrium ground state exhibits commensurate charge-density-wave (CDW) order for any nonzero coupling, we identify three distinct post-quench dynamical regimes as a function of the final electron-phonon coupling: a nonequilibrium metallic state without CDW order, an intermediate regime characterized by slow scale-invariant ordering dynamics, and a frozen CDW state with arrested coarsening and immobile kinks. We analyze the intermediate regime in detail and uncover an unconventional algebraic decay of the kink density,
n∼t−1/3, distinct from both ballistic annihilation and diffusive coarsening in classical dissipative systems. We show that this anomalous exponent arises from the hybrid nature of the dynamics: while the lattice evolves deterministically, the electronic degrees of freedom act as an effective internal bath that induces diffusive kink motion without energy dissipation. An effective reaction-diffusion description, incorporating both annihilation and elastic scattering of kinks, quantitatively accounts for the observed scaling behavior. Our results reveal a distinct coarsening mechanism in isolated hybrid systems, demonstrating how internal quantum dynamics can qualitatively reshape defect kinetics far from equilibrium.
We investigate the out-of-equilibrium dynamics of a photoexcited charge-density-wave (CDW) state in the square-lattice Holstein model, in a setup similar to a pump-probe experiment. At half-filling, the ground state of this system is characterized by a checkerboard modulation of particle densities, accompanied by a concomitant lattice distortion. An efficient real-space dynamics method integrating the von Neumann equation for electron density matrix and Newton equation for classical lattice dynamics is developed to simulate the dynamical evolution of a photoexcited Holstein system. We find that the energy injected by a short pump pulse results in the reduction of the CDW order and the generation of coherent phonons. At strong photoexcitations, while the CDW order is melted in the sense that the time-averaged order parameter vanishes, a dynamical CDW state is self-sustained by the strong coherent phonon oscillations. Our large-scale simulations further uncover a dynamical regime of intermediate fluence where complex pattern formation is induced by the pump pulse. The emergent spatial textures are characterized by super density modulations on top of the short-range checkerboard CDW order. Our demonstration of pattern formation highlights the significance of dynamical inhomogeneity in quantum many-body systems in pump-probe experiments.
We investigate the out-of-equilibrium dynamics of a photo-excited charge-density-wave (CDW) state in the square-lattice Holstein model, in a setup similar to a pump-probe experiment. At half-filling, the ground state of this system is characterized by a checkerboard modulation of particle densities, accompanied by a concomitant lattice distortion. An efficient real-space dynamics method integrating the von~Neumann equation for electron density matrix and Newton equation for classical lattice dynamics is developed to simulate the dynamical evolution of a photo-excited Holstein system. We find that the energy injected by a short pump pulse results in the reduction of the CDW order and the generation of coherent phonons. At strong photoexcitations, while the CDW order is melted in the sense that the time-averaged order parameter vanishes, a dynamical CDW state is self-sustained by the strong coherent phonon oscillations. Our large-scale simulations further uncover a dynamical regime of intermediate fluence where complex pattern formation is induced by the pump pulse. The emergent spatial textures are characterized by super density modulations on top of the short-range checkerboard CDW order. Our demonstration of pattern formation highlights the significance of dynamical inhomogeneity in quantum many-body systems in pump-probe experiments.
The Holstein Hamiltonian describes fermions hopping on a lattice and interacting locally with dispersionless phonon degrees of freedom. In the low-density limit, dressed quasiparticles, polarons and bipolarons, propagate with an effective mass. At higher densities, pairs can condense into a low-temperature superconducting phase and, at or near commensurate filling on a bipartite lattice, to charge density wave (CDW) order. CDW formation breaks a discrete symmetry and hence occurs via a second-order (Ising) transition and therefore at a finite
Tcdw in two dimensions. Quantum Monte Carlo calculations have determined
Tcdw for a variety of geometries, including square, honeycomb, and Lieb lattices. The superconducting transition, on the other hand, in
d=2 is in the Kosterlitz-Thouless universality class and is much less well characterized. In this paper we determine
Tsc for the square lattice for several values of the density $\ensuremath{\rho}$ and phonon frequency ${\ensuremath{\omega}}_{0}$. We find that quasilong-range order sets in at ${T}_{\mathrm{sc}}\ensuremath{\lesssim}t/20$, where
t is the near-neighbor hopping amplitude, consistent with previous rough estimates from simulations which extrapolated to only the temperatures we reach from considerably higher
T. We also show evidence of a discontinuous evolution of the density as the CDW transition is approached at half filling.