Search

The development of high-intensity ultrashort laser pulses unlocks the potential of pump-probe spectroscopy in sub-femtosecond timescale. Notably, subcycle pump pulses can generate electronic states unreachable by conventional multicycle pulses, leading to a phenomenon that we refer to as subcycle-pulse engineering. In this study, we employ the time-dependent density-matrix renormalization group method to unveil the transient absorption spectra of superconducting and metallic states in nearly half-filled one-dimensional and two-dimensional Hubbard models excited by an ultrashort halfcycle pulse, which can induce a current with inversion-symmetry breaking. In a superconducting state realized in attractive on-site interactions, we find the transient activation of absorptions at energies corresponding to the amplitude modes of superconducting and charge-density-wave states. On the other hand, in a metallic state realized in the two-dimensional model with repulsive on-site interactions, we obtain another type of absorption enhancements, which are distributed broadly in spin excitation energies. These findings indicate that superconducting and metallic states are sensitive to an ultrashort halfcycle pulse, leading to the transient activations of optical absorptions with their respective mechanisms., Comment: 7 pages, 4 figures

In periodically driven (Floquet) systems, evolution typically results in an infinite-temperature thermal state due to continuous energy absorption over time. However, before reaching thermal equilibrium, such systems may transiently pass through a meta-stable state known as a prethermal state. This prethermal state can exhibit phenomena not commonly observed in equilibrium, such as discrete time crystals (DTCs), making it an intriguing platform for exploring out-of-equilibrium dynamics. Here, we investigate the relaxation dynamics of initially prepared product states under periodic driving in a kicked Ising model using the IBM Quantum Heron processor, comprising 133 superconducting qubits arranged on a heavy-hexagonal lattice, over up to $100$ time steps. We identify the presence of a prethermal regime characterised by magnetisation measurements oscillating at twice the period of the Floquet cycle and demonstrate its robustness against perturbations to the transverse field. Our results provide evidence supporting the realisation of a period-doubling DTC in a two-dimensional system. Moreover, we discover that the longitudinal field induces additional amplitude modulations in the magnetisation with a period incommensurate with the driving period, leading to the emergence of discrete time quasicrystals (DTQCs). These observations are further validated through comparison with tensor-network and state-vector simulations. Our findings not only enhance our understanding of clean DTCs in two dimensions but also highlight the utility of digital quantum computers for simulating the dynamics of quantum many-body systems, addressing challenges faced by state-of-the-art classical simulations., Comment: 8 pages, 4 figures + Supplemental Material (17 pages, 20 figures)

Recent advancements in pulse laser technology have facilitated the exploration of non-equilibrium spectroscopy of electronic states in the presence of strong electric fields across a broad range of photon energies. The Keldysh crossover serves as an indicator that distinguishes between excitations resulting from photon absorption triggered by near-infrared multicycle pulses and those arising from quantum tunneling induced by terahertz pulses. Using time-dependent density-matrix renormalization group, we investigate the emergence of the Keldysh crossover in a one-dimensional (1D) Mott insulator. We find that the Drude weight is proportional to photo-doped doublon density when a pump pulse induces photon absorption. In contrast, the Drude weight is suppressed when a terahertz pulse introduces doublons and holons via quantum tunneling. The suppressed Drude weight accompanies glassy dynamics with suppressed diffusion, which is a consequence of strong correlations and exhibits finite polarization decaying slowly after pulse irradiation. In the quantum tunneling region, entanglement entropy slowly grows logarithmically. These contrasting behaviors between the photon-absorption and quantum tunneling regions are a manifestation of the Keldysh crossover in 1D Mott insulators and provide a novel methodology for designing the localization and symmetry of electronic states called subcycle-pulse engineering., Comment: 7 pages, 5 figures

The recent experimental observations of loop current in Sr$_{2}$IrO$_{4}$, YBa$_{2}$Cu$_{3}$O$_{7}$, and Sr$_{14}$Cu$_{24}$O$_{41}$ have inspired a theoretical study that broadly redefines loop current as a manifestation of quantum liquid crystals. Using the density-matrix renormalization group method, we investigate the emergence of spin loop-current (sLC) textures in carrier-doped (i) excitonic insulators, (ii) orbital-selective Mott insulators, and (iii) two-dimensional Mott insulators, modeled by a two-orbital Hubbard model on a ladder lattice in (i) and (ii) and a single-orbital Hubbard model on a square lattice in (iii). Calculating the spatial distribution of spin current around a bond to which a pinning field is applied, we find conditions for longer-ranged sLC correlations. In system (i), when using the model parameters employed to describe the excitonic condensation, we find that a sLC texture appears near half filling, associated with an excitonic condensation in a spin channel. In system (ii), using typical sets of model parameters for BaFe$_{2}$Se$_{3}$, we find that a sLC texture appears at electron fillings where a block-type antiferromagnetism develops. In system (iii), introducing a next-nearest-neighbor hopping $t'\sim -0.25$ (in units of the nearest-neighbor hopping) suggested for high-$T_\text{c}$ cuprates, we find that an axial-sLC texture emerges at hole-carrier density $\delta=0.125$, where the charge stripe simultaneously appears., Comment: 15 pages, 12 figures

Photoinduced nonequilibrium states in the Mott insulators reflect the fundamental nature of competition between itinerancy and localization of the charge degrees of freedom. The spin degrees of freedom will also contribute to the competition in a different manner depending on lattice geometry. We investigate pulse-excited optical responses of a half-filled two-leg Hubbard ladder and compare them with those of a one-dimensional extended Hubbard chain. Calculating the time-dependent optical conductivity, we find that strong mono-cycle pulse inducing quantum tunneling gives rise to anomalous suppression of photo-induced in-gap weight, leading to negative weight. This is in contrast to finite positive weight in the Hubbard chain. The origin of this anomalous behavior in the two-leg ladder is attributed to photoinduced localized exciton that reflects strong spin-singlet dimer correlation in the ground state., Comment: 6 pages, 4 figures, to appear in Phys. Rev. B

Owning to their high controllability, laser pulses have contributed greatly to our understanding of strongly correlated electron systems. However, typical multicycle pulses do not control the symmetry of systems that plays an important role in the emergence of novel quantum phases. Here, we demonstrate that subcycle pulses whose oscillation is less than one period within a pulse envelope can control inversion and time-reversal symmetries in the electronic states of the one-dimensional extended Hubbard model. Using an ultrashort subcycle pulse, one can generate a steady electric current (SEC) in a photoexcited state due to an Aharonov-Bohm flux instantaneously introduced through the phase of an electric field. Consequently, time-reversal symmetry is broken. In contrast, a broad subcycle pulse does not induce SEC but instead generates electric polarization, thus breaking inversion symmetry. Both symmetry breakings in a photoexcited state can be monitored by second harmonic generation. These findings provide a new methodology for designing the symmetries of electronic states and open up a new field of subcycle-pulse engineering., Comment: 8 pages, 7 figures

The photoexcitation of a Mott insulator on a square lattice weakens the intensity of both single- and two-magnon excitations as observed in time-resolved resonant-inelastic X-ray scattering and time-resolved Raman scattering, respectively. However, the spectral changes in the low-energy regions below the magnons have not yet been clearly understood. To uncover the nature of the photoinduced low-energy magnetic excitations of the Mott insulator, we numerically investigate the transient magnetic dynamics in a photoexcited half-filled Hubbard model on a square lattice. After turning off a pump pulse tuned for an absorption edge, new magnetic signals clearly emerge well below the magnon energy in both single- and two-magnon excitations. We find that the low-energy excitations are predominantly created via excitonic states at the absorption edge. These exciton-assisted magnetic excitations may provide a possible explanation for the low-energy spectral weight in a recent time-resolved two-magnon Raman scattering experiment for insulating YBa$_2$Cu$_3$O$_{6.1}$., Comment: 7 pages, 6 figures, 2 pages supplemental note with 3 figures

The elucidation of nonequilibrium states in strongly correlated systems holds the key to emergence of novel quantum phases. The nonequilibrium-induced insulator-to-metal transition is particularly interesting since it reflects the fundamental nature of competition between itinerancy and localization of the charge degrees of freedom. We investigate pulse-excited insulator-to-metal transition of the half-filled one-dimensional extended Hubbard model. Calculating the time-dependent optical conductivity with the time-dependent density-matrix renormalization group, we find that strong mono- and half-cycle pulses inducing quantum tunneling strongly suppress spectral weights contributing to the Drude weight $\sigma_\text{D}$, even if we introduce a large number of carriers $\Delta n_\text{d}$. This is in contrast to a metallic behavior of $\sigma_\text{D}\propto \Delta n_\text{d}$ induced by photon absorption and chemical doping. The strong suppression of $\sigma_\text{D}$ in quantum tunneling is a result of the emergence of the Hilbert-space fragmentation, which makes pulse-excited states glassy., Comment: 8 pages, 4 figures + Supplemental Material (11 pages, 5 figures)

The real part of optical conductivity $\text{Re}\sigma(\omega)$ of the Mott insulators has a large amount of information on how spin and charge degrees of freedom interact with each other. By using the time-dependent density-matrix renormalization group, we study $\text{Re}\sigma(\omega)$ of the two-dimensional Hubbard model on a square lattice at half filling. We find an excitonic peak at the Mott-gap edge of $\text{Re}\sigma(\omega)$ not only for the two-dimensional square lattice but also for two- and four-leg ladders. For the square lattice, however, we do not clearly find a gap between an excitonic peak and continuum band, which indicates that a bound state is not well defined. The emergence of an excitonic peak in $\text{Re}\sigma(\omega)$ implies the formation of a spin polaron. Examining the dependence of $\text{Re}\sigma(\omega)$ on the on-site Coulomb interaction and next-nearest neighbor hoppings, we confirm that an excitonic peak is generated from a magnetic effect. Electron scattering due to an electron-phonon interaction is expected to easily suppress an excitonic peak since spectral width of an excitonic peak is very narrow. Introducing a large broadening in $\text{Re}\sigma(\omega)$ by modeling the electron-phonon coupling present in La$_{2}$CuO$_{4}$ and Nd$_{2}$CuO$_{4}$, we obtain $\text{Re}\sigma(\omega)$ comparable with experiments., Comment: 10 pages, 8 figures

Optical measurements in doped Mott insulators have discovered the emergence of spectral weights at mid-infrared (MIR) upon chemical doping and photodoping. MIR weights may have a relation to string-type excitation of spins, which is induced by a doped hole generating misarranged spins with respect to their sublattice. There are two types of string effects: one is an $S^z$ string that is repairable by quantum spin flips and the other is a phase string irreparable by the spin flips. We investigate the effect of $S^{z}$ and phase strings on MIR weights. Calculating the optical conductivity of the single-hole Hubbard model in the strong-coupling regime and the $t$-$J$ model on two-leg ladders by using time-dependent Lanczos and density-matrix renormalization group, we find that phase strings make a crucial effect on the emergence of MIR weights as compared with $S^{z}$ strings. Our findings indicate that a mutual Chern-Simons gauge field acting between spin and charge degrees of freedom, which is the origin of phase strings, is significant for obtaining MIR weights. Conversely, if we remove this gauge field, no phase is picked up by a doped hole. As a result, a spin-polaron accompanied by a local spin distortion emerges and a quasiparticle with a cosine-like energy dispersion is formed in single-particle spectral function. Furthermore, we suggest a Floquet engineering to examine the phase-string effect in cold atoms., Comment: 14 pages, 9 figures

Motivated by the recent development of time-resolved resonant-inelastic x-ray scattering (TRRIXS) in photoexcited antiferromagnetic Mott insulators, we numerically investigate momentum-dependent transient spin dynamics in a half-filled Hubbard model on a square lattice. After turning off a pumping photon pulse, the intensity of a dynamical spin structure factor temporally oscillates with frequencies determined by the energy of two magnons in the antiferromagnetic Mott insulator. We find an antiphase behavior in the oscillations between two orthogonal momentum directions, parallel and perpendicular to the electric field of a pump pulse. The phase difference comes from the $B_{1g}$ channel of the two-magnon excitation. Observing the antiphase oscillations will be a big challenge for TRRIXS experiments when their time resolution will be improved by more than an order of magnitude., Comment: 5 pages, 4 figures

We propose an effective model called the "charge model", for the half-filled one-dimensional Hubbard and extended Hubbard models. In this model, spin-charge separation, which has been justified from an infinite on-site repulsion ($U$) in the strict sense, is compatible with charge fluctuations. Our analyses based on the many-body Wannier functions succeeded in determining the optical conductivity spectra in large systems. The obtained spectra reproduce the spectra for the original models well even in the intermediate $U$ region of $U=5-10T$, with $T$ being the nearest-neighbor electron hopping energy. These results indicate that the spin-charge separation works fairly well in this intermediate $U$ region against the usual expectation and that the charge model is an effective model that applies to actual quasi-one-dimensional materials classified as strongly correlated electron systems., Comment: 36 pages, 11 figures

We demonstrate that supervised machine learning (ML) with entanglement spectrum can give useful information for constructing phase diagram in the half-filled one-dimensional extended Hubbard model. Combining ML with infinite-size density-matrix renormalization group, we confirm that bond-order-wave phase remains stable in the thermodynamic limit., Comment: 2 pages, 2 figures

Photoinduced nonequilibrium states can provide new insight into dynamical properties of strongly correlated electron systems. One of the typical and extensively studied systems is the half-filled one-dimensional extended Hubbard model (1DEHM). Here, we propose that the supervised machine learning (ML) can provide useful information for characterizing photoexcited states in 1DEHM. Using entanglement spectra as a training dataset, we construct neural network. Judging from the trained network, we find that bond-spin-density wave (BSDW) order can be enhanced in photoexcited states if the frequency of a driving pulse nearly resonates with gap. We separately calculate the time evolution of local and non-local order parameters and confirm that the correlation functions of BSDW are actually enhanced by photoexcitation as predicted by ML. The successful prediction of BSDW demonstrates the advantage of ML to assist characterizing photoexcited quantum states., Comment: 11 pages, 10 figures

Ultrafast pump-probe spectroscopy with high temporal and spectral resolutions provides new insight into ultrafast nonequilibrium phenomena. We propose that transient interference between pump and probe pulses is realized in pump-probe spectroscopy of band and Mott insulators, which can be observed only after recent developments of ultrafast spectroscopic techniques. A continuum structure in the excitation spectrum of band insulators is found to act as a medium for storing the spectral information of the pump pulse, and the spectrum detected by the probe pulse is interfered with by the medium, generating the transient interference in the energy domain. We also demonstrate the transient interference in the presence of electron correlations in a one-dimensional half-filled Hubbard model. Furthermore, bosons coupled to electrons additively contribute to the interference. Our finding will provide an interpretation of probe-energy-dependent oscillations recently observed in the pump-probe spectrum for a two-dimensional Mott insulator., Comment: 11 pages, 2 figures

We theoretically investigate pump-probe optical responses in the two-dimensional extended Hubbard model describing cuprates by using a time-dependent Lanczos method. At half filling, pumping generates photoinduced absorptions inside the Mott gap. A part of low-energy absorptions is attributed to the independent propagation of photoinduced holons and doublons. The spectral weight just below the Mott gap increases with decreasing the on-site Coulomb interaction $U$. We find that the next-nearest-neighbor Coulomb interaction $V_1$ enhances this $U$ dependence, indicating the presence of biexcitonic contributions formed by two holon-doublon pairs. Photo-pumping in hole-doped systems also induces spectral weights below remnant Mott-gap excitations, being consistent with recent experiments. The induced weights are less sensitive to $V_1$ and may be related to the formation of a biexcitonic state in the presence of hole carriers., Comment: 7pages, 6figures

We investigate the ground state properties of $\mathrm{Na_2IrO_3}$ based on numerical calculations of the recently proposed {\it ab initio} Hamiltonian represented by Kitaev and extended Heisenberg interactions. To overcome the limitation posed by small tractable system sizes in the exact diagonalization study employed in a previous study (Yamaji $\textit{et al.}$, Phys. Rev. Lett. $\textbf{113}$, 107201 (2014)), we apply two-dimensional density matrix renomalization group, and infinite-size tensor-network method. By calculating at much larger system sizes, we critically test the validity of the exact diagonalization results. The results consistently indicate that the ground state of $\mathrm{Na_2IrO_3}$ is a magnetically ordered state with zigzag configuration in agreement with experimental observations and the previous diagonalization study. Applications of the two independent methods in addition to the exact diagonalization study further uncover a consistent and rich phase diagram near the zigzag phase beyond the accessibility of the exact diagonalization. For example, in the parameter space away from the $\textit{ab initio}$ value of $\mathrm{Na_2IrO_3}$ controlled by the trigonal distortion, we find three phases: (i) an ordered phase with the magnetic moment aligned mutually in $120$ degrees orientation on every third hexagon, (ii) a magnetically ordered phase with a $16$-site unit-cell, and (iii) an ordered phase with presumably incommensurate periodicity of the moment. It suggests that potentially rich magnetic structures may appear in $A_\mathrm{2}\mathrm{IrO_3}$ compounds for $A$ other than Na. The present results also serve to establish the accuracy of the first-principles approach in reproducing the available experimental results thereby further contribute to find a route to realize the Kitaev spin liquid., Comment: 14 pages, 11 figures

We study the Kitaev--Heisenberg model on a triangular lattice by using the two-dimensional density-matrix renormalization group method. Calculating the ground-state energy and spin structure factors, we obtain a ground-state phase diagram of the Kitaev--Heisenberg model. As suggested by previous studies, we find a 120$^\circ$ antiferromagnetic (AFM) phase, a $\mathbb{Z}_2$-vortex crystal phase, a nematic phase, a dual $\mathbb{Z}_2$-vortex crystal phase (the dual counterpart of the $\mathbb{Z}_2$-vortex crystal phase), a $\mathbb{Z}_6$ ferromagnetic phase, and a dual ferromagnetic phase (the dual counterpart of the $\mathbb{Z}_6 $ ferromagnetic phase). Spin correlations discontinuously change at phase boundaries because of first-order phase transitions. We also study the relation among the von Neumann entanglement entropy, entanglement spectrum, and phase transitions of the model. We find that the Schmidt gap closes at phase boundaries and thus the entanglement entropy clearly changes as well. This is different from the Kitaev--Heisenberg model on a honeycomb lattice, where the Schmidt gap and entanglement entropy are not necessarily a good measure of phase transitions., Comment: 10 pages, 12 figures

We study an extended Kitaev-Heisenberg model including additional anisotropic couplings by using two-dimensional density-matrix renormalization group method. Calculating the gound-state energy, entanglement entropy, and spin-spin correlation functions, we make a phase diagram of the extended Kitaev-Heisenberg model around spin-liquid phase. We find a zigzag antiferromagnetic phase, a ferromagnetic phase, a 120-degree antiferromagnetic phase, and two kinds of incommensurate phases around the Kitaev spin-liquid phase. Furthermore, we study the entanglement spectrum of the model and find that entanglement levels in the Kitaev spin-liquid phase are degenerate forming pairs but those in the magnetically ordered phases are non-degenerate. The Schmidt gap defined as the energy difference between the lowest two levels changes at the phase boundary adjacent to the Kitaev spin-liquid phase. However, we find that phase boundaries between magnetically ordered phases do not necessarily agree with the change of the Schmidt gap., Comment: 11 pages, 15 figures