Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2407.01936 · Physical Review Letters

Many-Body Effects on High-Harmonic Generation in Hubbard Ladders

Yuta Murakami, Thomas Hansen, Shintaro Takayoshi, Lars Bojer Madsen, Philipp Werner

Published Mar 7, 2025• 7 citations
arXivHTMLPDF

Abstract

We unveil multifaceted effects of spin-charge couplings on the high-harmonic generation (HHG) in Mott insulators by analyzing the two-leg ladder Hubbard model, where spin dynamics activated by the interchain hopping t_{y} drastically modifies the HHG features. When the two chains are decoupled (t_{y}=0), HHG originates from the dynamics of coherent doublon-holon pairs because of spin-charge separation. With increasing t_{y}, the doublon-holon pairs lose their coherence due to their interchain hopping and resultant nonlocal spin strings. Furthermore, the HHG signal from spin polarons-charges dressed by spin clouds-leads to an additional plateau in the HHG spectrum. For large t_{y}, we identify unconventional HHG processes involving three elementary excitations-two polarons and one magnon. Our study is an essential step toward a microscopic understanding of nontrivial many-body effects on HHG in correlated materials beyond conventional semiconductors.

DOI
10.1103/physrevlett.134.096504

Related papers

High-harmonic generation in one-dimensional Mott insulators

Yuta Murakami, Shintaro Takayoshi, Akihisa Koga +1 more

We study high-harmonic generation (HHG) in the one-dimensional Hubbard model in order to understand its relation to elementary excitations as well as the similarities and differences to semiconductors. The simulations are based on the infinite time-evolving block decimation (iTEBD) method and exact diagonalization. We clarify that the HHG originates from the doublon-holon recombination, and the scaling of the cutoff frequency is consistent with a linear dependence on the external field. We demonstrate that the subcycle features of the HHG can be reasonably described by a phenomenological three step model for a doublon-holon pair. We argue that the HHG in the one-dimensional Mott insulator is closely related to the dispersion of the doublon-holon pair with respect to its relative momentum, which is not necessarily captured by the single-particle spectrum due to the many-body nature of the elementary excitations. For the comparison to semiconductors, we introduce effective models obtained from the Schrieffer-Wolff transformation, i.e., a strong-coupling expansion, which allows us to disentangle the different processes involved in the Hubbard model: intraband dynamics of doublons and holons, interband dipole excitations, and spin exchanges. These demonstrate the formal similarity of the Mott system to the semiconductor models in the dipole gauge, and reveal that the spin dynamics, which does not directly affect the charge dynamics, can reduce the HHG intensity. We also show that the long-range component of the intraband dipole moment has a substantial effect on the HHG intensity, while the correlated hopping terms for the doublons and holons essentially determine the shape of the HHG spectrum. A new numerical method to evaluate single-particle spectra within the iTEBD method is also introduced.

arXiv:2010.13255•Physical review. B./Physical review. B•Jan 8, 2021•64 citations

Tracking doublon-holon dynamics in high-harmonic generation from Mott insulators

Tao-Yuan Du, Hui-Ru Li, Bo Li +1 more

High-harmonic generation (HHG) in strongly correlated Mott insulators is investigated using exact diagonalization and time-dependent density-matrix propagation of a laser-driven one-dimensional Hubbard chain. By projecting onto equilibrium Hubbard bands, we use the doublon population and its dynamics as a diagnostic to analyze intraband (spin-wave-like) and interband (doublon-holon creation) excitation channels. A filling-dependent crossover emerges: Bloch-like intraband response at dilute filling, mixed dynamics at intermediate filling, and interband-dominated HHG with plateau and cutoff near half filling. In the considered parameter range, increasing interaction strength $U$ strongly suppresses interband contributions through the enlarged Mott gap and correlation-induced localization. Intra- and interband current decomposition reveals opposing flows below the Mott gap (Δ_{\mathrm{Mott}}) and selective dephasing suppression of interband coherence, enhancing net doublon accumulation. Time-frequency analysis uncovers the filling-dependent features of quantum trajectories, manifesting in distinct below-Δ_{\mathrm{Mott}} emission. This doublon-based analysis provides a transparent link between equilibrium spin-charge separation and nonequilibrium strong-field response, and clarifies how dephasing modifies interband coherence and doublon accumulation.

arXiv:2605.01388•Physical review. B./Physical review. B•Jun 16, 2026•0 citations

Doping effects in high-harmonic generation from correlated systems

Thomas Hansen, Lars Bojer Madsen

Using the one-dimensional Hubbard model, which is commonly used for describing, e.g., high-${T}_{c}$ superconducting cuprates, we study high-harmonic generation (HHG) from doped, correlated materials. Doping is modeled by changing the number of electrons in the lattice from the conventional half-filling case. For relatively small Hubbard $U$, i.e., small electron-electron correlation, we find little to no effect of doping on the dynamics and the HHG spectra. For increasing $U$ the degree of doping has a marked effect on the dynamics and spectra. We explain these findings through the quasiparticle-based doublon-holon picture. The dynamics are separated into two types, first, doublon and holon movement, and, second, doublon-holon pair creation and annihilation. Doping results in all configurations containing doublons or holons. Those quasiparticles can move at no extra cost in energy regardless of the correlation level. This motion at no energy cost increases the high-harmonic gain for low- and medium-harmonic orders. We discuss that in the high-$U$ limit, antiferromagnetic ordering becomes increasingly unlikely with increasing doping rates and explain an associated drop in the high-order harmonics relative to the case of half-filling.

arXiv:2209.04311•Physical review. B./Physical review. B•Dec 20, 2022•21 citations

Spectral function of Mott-insulating Hubbard ladders: From fractionalized excitations to coherent quasiparticles

Chun Yang, Adrian E. Feiguin

We study the spectral function of two-leg Hubbard ladders with the time-dependent density matrix renormalization group method (tDMRG). The high-resolution spectrum displays features of spin-charge separation and a scattering continuum of excitations with coherent bands of bound states ``leaking'' from it. As the interleg hopping is increased, the continuum in the bonding channel moves to higher energies and spinon and holon branches merge into a single coherent quasiparticle band. Simultaneously, the spectrum undergoes a crossover from a regime with two minima at incommensurate values of ${k}_{x}$ (a Mott insulator), to one with a single minimum at ${k}_{x}=\ensuremath{\pi}$ (a band insulator). We identify the presence of a continuum of scattering states consisting of a triplon and a polaron. We analyze the processes leading to quasiparticle formation by studying the time evolution of charge and spin degrees of freedom in real space after the hole is created. At short times, incoherent holons and spinons are emitted but after a characteristic time $\ensuremath{\tau}$ charge and spin form polarons that propagate coherently.

arXiv:1903.08267•Physical review. B./Physical review. B•Jun 10, 2019•14 citations

Excitonic effects on high-harmonic generation in Mott insulators

Mina Udono, K. Sugimoto, Tatsuya Kaneko +1 more

To study excitonic effects on high-harmonic generation (HHG) in Mott insulators, we investigate pumped nonequilibrium dynamics in the one-dimensional extended Hubbard model. By employing time-dependent calculations based on the exact diagonalization and infinite time-evolving block decimation methods, we find the strong enhancement of the HHG intensity around the exciton energy. The subcycle analysis in the sub-Mott-gap regime shows that the intensity region of the time-resolved spectrum around the exciton energy splits into two levels and oscillates following the driving electric field. This excitonic dynamics is qualitatively different from the dynamics of free doublon and holon but favorably contributes to HHG in the Mott insulator.

arXiv:2205.06390•Physical review. B./Physical review. B•Jun 10, 2022•29 citations

Correlation effects in low dimensional quantum systems

Chun Yang

This dissertation discussed four low-dimensional systems: Hubbard ladder and its 2D extension, anisotropic Hubbard ladder, one-dimensional two-orbital t-J model, and Kondo model. All of them are of strongly correlated. Doped cuprates, as strongly correlated systems, have attracted a lot of attention in the condensed matter physics community, not just from the purely theoretical perspective, but also due to the novel applications beyond the conventional semiconductor paradigm. They display complicated phase diagram, defying the conventional rigid band theory, shaking our common sense and driving us to wrestle with many-body physics. With the aid of advanced experimental techniques such as ARPES, its various exotic electronic structure, such as waterfalls, hole pockets, Fermi arcs and kinks have been unveiled and remain open to interpretation. To attack such complicated problems, people usually build various models to mimic the CuO2 planes, which is the common feature of various cuprate compounds and is believed to be mainly responsible for their complicated properties, especially the high temperature superconductivity. In Chapter 3, we simulated the spectral properties of the 2- leg single-band Hubbard ladder and its 2-dimensional extension at and away from half-lling using the numerical tool: the combination of Cluster Perturbation Theory (CPT) with the time-dependent density matrix renormalization group method (tDMRG) as a solver. The main features of the spectrum can be described with a mean-eld dispersion, with kinks and pseudogap traced back to scattering between spin and charge degrees of freedom. As a bonus, we nd that the spectrum of hubbard ladder involves spin-polaron, the bound state of spinon and holon. Interacting one-dimensional electron systems generically behave dierent from higher dimensional systems in the sense that its low-energy excitations are described by Luttinger liquid theory as collective bosonic modes carrying spin and charge separately, the so-called spinons and holons, rather than the coherent quasi-particles in Fermi liquid theory. This phenomenon is known as spin-charge separation, an intrinsic property of 1D system. It can be thought of as a particular case of electron fractionalization. When multiple 1D chains get coupled, how the picture changes? In Chapter 4, we study the spectral function of the anisotropic 2-leg Hubbard ladder, one of the simplest multiple-chain systems. The spectrum can be understood as the scattering resultants of spinon and holons: As the inter-chain hopping is increased, the continuum in the bonding channel moves to higher energies and spinon and holon branches merge into a single coherent quasi-particle band. Simultaneously, the spectrum undergoes a crossover from a regime with two minima at incommensurate values of kx (a Mott insulator), to one with a single minimum at kx = (a band insulator). In addition, the spectrum contains a continuum of scattering states consisting of a triplon and a polaron, accompanied by the formation of bound states with higher spin 3 2 , ferromagnetic \spin bags" in which the hole can move freely. To analyze the formation of quasi-particle, we study the time-evolution of charge and spin degrees of freedom in real space after a hole is injected, which indicates that incoherent holons and spinons can emerge into polarons, that can propagate coherently, after a characterstic time . Multiorbital systems can involve Mott and Hund physics, and are rich in various phenomena, including orbital selective Mott transitions, spin-orbital separations, and pairing density wave due to the complex interplay between spin, charge and orbital degrees of freedom. Therefore, it can capture many unique features in real strongly correlated systems, such as iron pnictides and ruthenates, and thus attract a lot of investigations. In Chapter 5, we comprehensively study a one-dimensional two-orbital model at and below quarter lling, which realizes a number of unconventional phases, by means of DMRG and some analytical calculations. In particular, we nd that the system supports an extensive density wave in which excitons quasicondense with nite center-of-mass momentum. In this phase, excitons behave as hardcore bosons without charge order. In addition, excitons can pair to form biexcitons in a state that is close to a charge density wave instability. Finally, a genuine orbital-selective paired state (metallic plus spin-gapped superconductor) can be realized when pairing dominates over the interorbital repulsion. Kondo problem, as one of the most studied and best understood correlated systems, has been studied by various techniques, but the internal structure of its ground state has two deeply con icting pictures: At strong coupling limit, the system is essentially composed of a scattering center, which is a tightly bound singlet formed by the impurity and a localized electron at r0, and the remaining conduction electrons that are decoupled to the impurity and remain almost \free". At weak enough couplings, the singlet extends to the entire volume and the impurity will couple mostly to one electron at the Fermi level. In chapter 6, we provide a coherent picture to unveil the puzzle in the context of quantum information theory by deconstructing the Kondo singlet, both in real and energy space, and measuring the spin correlations between the impurity and the invididual free electron eigenstate based on two dierent measures. As a bonus, we can identify \projected natural orbitals" that can be used as a \good" approximate wavefunctions for simulations, and a quasi-order parameter that determines the crossover from weak to strong coupling regimes and furthermore oers the means to measure the Kondo screening length K. This screening length might be used to characterize the \size" of the Kondo cloud.

•Jan 1, 2018•0 citations
We unveil multifaceted effects of spin-charge couplings on the high-harmonic generation (HHG) in Mott insulators by analyzing the two-leg ladder Hubbard model, where spin dynamics activated by the interchain hopping t_{y} drastically modifies the HHG features. When the two chains are decoupled (t_{y}=0), HHG originates from the dynamics of coherent doublon-holon pairs because of spin-charge separation. With increasing t_{y}, the doublon-holon pairs lose their coherence due to their interchain hopping and resultant nonlocal spin strings. Furthermore, the HHG signal from spin polarons-charges dressed by spin clouds-leads to an additional plateau in the HHG spectrum. For large t_{y}, we identify unconventional HHG processes involving three elementary excitations-two polarons and one magnon. Our study is an essential step toward a microscopic understanding of nontrivial many-body effects on HHG in correlated materials beyond conventional semiconductors.
We study high-harmonic generation (HHG) in the one-dimensional Hubbard model in order to understand its relation to elementary excitations as well as the similarities and differences to semiconductors. The simulations are based on the infinite time-evolving block decimation (iTEBD) method and exact diagonalization. We clarify that the HHG originates from the doublon-holon recombination, and the scaling of the cutoff frequency is consistent with a linear dependence on the external field. We demonstrate that the subcycle features of the HHG can be reasonably described by a phenomenological three step model for a doublon-holon pair. We argue that the HHG in the one-dimensional Mott insulator is closely related to the dispersion of the doublon-holon pair with respect to its relative momentum, which is not necessarily captured by the single-particle spectrum due to the many-body nature of the elementary excitations. For the comparison to semiconductors, we introduce effective models obtained from the Schrieffer-Wolff transformation, i.e., a strong-coupling expansion, which allows us to disentangle the different processes involved in the Hubbard model: intraband dynamics of doublons and holons, interband dipole excitations, and spin exchanges. These demonstrate the formal similarity of the Mott system to the semiconductor models in the dipole gauge, and reveal that the spin dynamics, which does not directly affect the charge dynamics, can reduce the HHG intensity. We also show that the long-range component of the intraband dipole moment has a substantial effect on the HHG intensity, while the correlated hopping terms for the doublons and holons essentially determine the shape of the HHG spectrum. A new numerical method to evaluate single-particle spectra within the iTEBD method is also introduced.
High-harmonic generation (HHG) in strongly correlated Mott insulators is investigated using exact diagonalization and time-dependent density-matrix propagation of a laser-driven one-dimensional Hubbard chain. By projecting onto equilibrium Hubbard bands, we use the doublon population and its dynamics as a diagnostic to analyze intraband (spin-wave-like) and interband (doublon-holon creation) excitation channels. A filling-dependent crossover emerges: Bloch-like intraband response at dilute filling, mixed dynamics at intermediate filling, and interband-dominated HHG with plateau and cutoff near half filling. In the considered parameter range, increasing interaction strength UUU strongly suppresses interband contributions through the enlarged Mott gap and correlation-induced localization. Intra- and interband current decomposition reveals opposing flows below the Mott gap (Δ_{\mathrm{Mott}}) and selective dephasing suppression of interband coherence, enhancing net doublon accumulation. Time-frequency analysis uncovers the filling-dependent features of quantum trajectories, manifesting in distinct below-Δ_{\mathrm{Mott}} emission. This doublon-based analysis provides a transparent link between equilibrium spin-charge separation and nonequilibrium strong-field response, and clarifies how dephasing modifies interband coherence and doublon accumulation.
Using the one-dimensional Hubbard model, which is commonly used for describing, e.g., high-Tc{T}_{c}Tc​ superconducting cuprates, we study high-harmonic generation (HHG) from doped, correlated materials. Doping is modeled by changing the number of electrons in the lattice from the conventional half-filling case. For relatively small Hubbard UUU, i.e., small electron-electron correlation, we find little to no effect of doping on the dynamics and the HHG spectra. For increasing UUU the degree of doping has a marked effect on the dynamics and spectra. We explain these findings through the quasiparticle-based doublon-holon picture. The dynamics are separated into two types, first, doublon and holon movement, and, second, doublon-holon pair creation and annihilation. Doping results in all configurations containing doublons or holons. Those quasiparticles can move at no extra cost in energy regardless of the correlation level. This motion at no energy cost increases the high-harmonic gain for low- and medium-harmonic orders. We discuss that in the high-UUU limit, antiferromagnetic ordering becomes increasingly unlikely with increasing doping rates and explain an associated drop in the high-order harmonics relative to the case of half-filling.
We study the spectral function of two-leg Hubbard ladders with the time-dependent density matrix renormalization group method (tDMRG). The high-resolution spectrum displays features of spin-charge separation and a scattering continuum of excitations with coherent bands of bound states ``leaking'' from it. As the interleg hopping is increased, the continuum in the bonding channel moves to higher energies and spinon and holon branches merge into a single coherent quasiparticle band. Simultaneously, the spectrum undergoes a crossover from a regime with two minima at incommensurate values of kx{k}_{x}kx​ (a Mott insulator), to one with a single minimum at ${k}_{x}=\ensuremath{\pi}$ (a band insulator). We identify the presence of a continuum of scattering states consisting of a triplon and a polaron. We analyze the processes leading to quasiparticle formation by studying the time evolution of charge and spin degrees of freedom in real space after the hole is created. At short times, incoherent holons and spinons are emitted but after a characteristic time $\ensuremath{\tau}$ charge and spin form polarons that propagate coherently.
To study excitonic effects on high-harmonic generation (HHG) in Mott insulators, we investigate pumped nonequilibrium dynamics in the one-dimensional extended Hubbard model. By employing time-dependent calculations based on the exact diagonalization and infinite time-evolving block decimation methods, we find the strong enhancement of the HHG intensity around the exciton energy. The subcycle analysis in the sub-Mott-gap regime shows that the intensity region of the time-resolved spectrum around the exciton energy splits into two levels and oscillates following the driving electric field. This excitonic dynamics is qualitatively different from the dynamics of free doublon and holon but favorably contributes to HHG in the Mott insulator.
This dissertation discussed four low-dimensional systems: Hubbard ladder and its 2D extension, anisotropic Hubbard ladder, one-dimensional two-orbital t-J model, and Kondo model. All of them are of strongly correlated. Doped cuprates, as strongly correlated systems, have attracted a lot of attention in the condensed matter physics community, not just from the purely theoretical perspective, but also due to the novel applications beyond the conventional semiconductor paradigm. They display complicated phase diagram, defying the conventional rigid band theory, shaking our common sense and driving us to wrestle with many-body physics. With the aid of advanced experimental techniques such as ARPES, its various exotic electronic structure, such as waterfalls, hole pockets, Fermi arcs and kinks have been unveiled and remain open to interpretation. To attack such complicated problems, people usually build various models to mimic the CuO2 planes, which is the common feature of various cuprate compounds and is believed to be mainly responsible for their complicated properties, especially the high temperature superconductivity. In Chapter 3, we simulated the spectral properties of the 2- leg single-band Hubbard ladder and its 2-dimensional extension at and away from half-lling using the numerical tool: the combination of Cluster Perturbation Theory (CPT) with the time-dependent density matrix renormalization group method (tDMRG) as a solver. The main features of the spectrum can be described with a mean-eld dispersion, with kinks and pseudogap traced back to scattering between spin and charge degrees of freedom. As a bonus, we nd that the spectrum of hubbard ladder involves spin-polaron, the bound state of spinon and holon. Interacting one-dimensional electron systems generically behave dierent from higher dimensional systems in the sense that its low-energy excitations are described by Luttinger liquid theory as collective bosonic modes carrying spin and charge separately, the so-called spinons and holons, rather than the coherent quasi-particles in Fermi liquid theory. This phenomenon is known as spin-charge separation, an intrinsic property of 1D system. It can be thought of as a particular case of electron fractionalization. When multiple 1D chains get coupled, how the picture changes? In Chapter 4, we study the spectral function of the anisotropic 2-leg Hubbard ladder, one of the simplest multiple-chain systems. The spectrum can be understood as the scattering resultants of spinon and holons: As the inter-chain hopping is increased, the continuum in the bonding channel moves to higher energies and spinon and holon branches merge into a single coherent quasi-particle band. Simultaneously, the spectrum undergoes a crossover from a regime with two minima at incommensurate values of kx (a Mott insulator), to one with a single minimum at kx = (a band insulator). In addition, the spectrum contains a continuum of scattering states consisting of a triplon and a polaron, accompanied by the formation of bound states with higher spin 3 2 , ferromagnetic \spin bags" in which the hole can move freely. To analyze the formation of quasi-particle, we study the time-evolution of charge and spin degrees of freedom in real space after a hole is injected, which indicates that incoherent holons and spinons can emerge into polarons, that can propagate coherently, after a characterstic time . Multiorbital systems can involve Mott and Hund physics, and are rich in various phenomena, including orbital selective Mott transitions, spin-orbital separations, and pairing density wave due to the complex interplay between spin, charge and orbital degrees of freedom. Therefore, it can capture many unique features in real strongly correlated systems, such as iron pnictides and ruthenates, and thus attract a lot of investigations. In Chapter 5, we comprehensively study a one-dimensional two-orbital model at and below quarter lling, which realizes a number of unconventional phases, by means of DMRG and some analytical calculations. In particular, we nd that the system supports an extensive density wave in which excitons quasicondense with nite center-of-mass momentum. In this phase, excitons behave as hardcore bosons without charge order. In addition, excitons can pair to form biexcitons in a state that is close to a charge density wave instability. Finally, a genuine orbital-selective paired state (metallic plus spin-gapped superconductor) can be realized when pairing dominates over the interorbital repulsion. Kondo problem, as one of the most studied and best understood correlated systems, has been studied by various techniques, but the internal structure of its ground state has two deeply con icting pictures: At strong coupling limit, the system is essentially composed of a scattering center, which is a tightly bound singlet formed by the impurity and a localized electron at r0, and the remaining conduction electrons that are decoupled to the impurity and remain almost \free". At weak enough couplings, the singlet extends to the entire volume and the impurity will couple mostly to one electron at the Fermi level. In chapter 6, we provide a coherent picture to unveil the puzzle in the context of quantum information theory by deconstructing the Kondo singlet, both in real and energy space, and measuring the spin correlations between the impurity and the invididual free electron eigenstate based on two dierent measures. As a bonus, we can identify \projected natural orbitals" that can be used as a \good" approximate wavefunctions for simulations, and a quasi-order parameter that determines the crossover from weak to strong coupling regimes and furthermore oers the means to measure the Kondo screening length K. This screening length might be used to characterize the \size" of the Kondo cloud.