We define, compute and analyze the nonequilibrium differential optical conductivity of the one-dimensional extended Hubbard model at half-filling after applying a pump pulse, using the time-dependent density matrix renormalization group method. The melting of the Mott insulator is accompanied by a suppression of the local magnetic moment and ensuing photogeneration of doublon-holon pairs. The differential optical conductivity reveals $(i)$ mid-gap states related to parity-forbidden optical states, and $(ii)$ strong renormalization and hybridization of the excitonic resonance and the absorption band, yielding a Fano resonance. We offer evidence and interpret such a resonance as a signature of nonequilibrium optical excitations resembling excitonic strings, (bi)excitons, and unbound doublon-holon pairs, depending on the magnitude of the intersite Coulomb repulsion. We discuss our results in the context of pump and probe spectroscopy experiments on organic Mott insulators.
We study the photoinduced breakdown of a two-orbital Mott insulator and resulting metallic state. Using time-dependent density matrix renormalization group, we scrutinize the real-time dynamics of the half-filled two-orbital Hubbard model interacting with a resonant radiation field pulse. The breakdown, caused by production of doublon-holon pairs, is enhanced by Hunds exchange, which dynamically activates large orbital fluctuations. The melting of the Mott insulator is accompanied by a high to low spin transition with a concomitant reduction of antiferromagnetic spin fluctuations. Most notably, the overall time response is driven by the photogeneration of excitons with orbital character that are stabilized by Hunds coupling. These unconventional Hund excitons correspond to bound spin-singlet orbital-triplet doublon-holon pairs. We study exciton properties such as bandwidth, binding potential, and size within a semiclassical approach. The photometallic state results from a coexistence of Hund excitons and doublon-holon plasma.
The generalization of matrix product states (MPS) to continuous systems, as proposed in the breakthrough paper [F. Verstraete, J.I. Cirac, Phys. Rev. Lett. 104, 190405(2010)], provides a powerful variational ansatz for the ground state of strongly interacting quantum field theories in one spatial dimension. A continuous MPS (cMPS) approximation to the ground state can be obtained by simulating an Euclidean time evolution. In this Letter we propose a cMPS optimization algorithm based instead on energy minimization by gradient methods, and demonstrate its performance by applying it to the Lieb Liniger model (an integrable model of an interacting bosonic field) directly in the thermodynamic limit. We observe a very significant computational speed-up, of more than two orders of magnitude, with respect to simulating an Euclidean time evolution. As a result, much larger cMPS bond dimension D can be reached (e.g. D = 256 with moderate computational resources) thus helping unlock the full potential of the cMPS representation for ground state studies.
The Lieb-Liniger model describes one-dimensional bosons interacting through a repulsive contact potential. In this work, we introduce an extended version of this model by replacing the contact potential with a decaying exponential. Using the recently developed continuous matrix product states techniques, we explore the ground state phase diagram of this model by examining the superfluid and density correlation functions. At weak coupling superfluidity governs the ground state, in a similar way as in the Lieb-Liniger model. However, at strong coupling quasi-crystal and super-Tonks-Girardeau regimes are also found, which are not present in the original Lieb-Liniger case. Therefore the presence of the exponentially-decaying potential leads to a superfluid/super-Tonks-Girardeau/quasi-crystal crossover, when tuning the coupling strength from weak to strong interactions. This corresponds to a Luttinger liquid parameter in the range $K in (0, infty)$; in contrast with the Lieb-Liniger model, where $K in [1, infty)$, and the screened long-range potential, where $K in (0, 1]$.
Using the time-dependent density-matrix renormalization group (tDMRG), we study the time evolution of electron wave packets in one-dimensional (1D) metal-superconductor heterostructures. The results show Andreev reflection at the interface, as expected. By combining these results with the well-known single-spin-species electron-hole transformation in the Hubbard model, we predict an analogous spin Andreev reflection in metal-Mott insulator heterostructures. This effect is numerically confirmed using 1D tDMRG, but it is expected to be present also in higher dimensions, as well as in more general Hamiltonians. We present an intuitive picture of the spin reflection, analogous to that of Andreev reflection at metal-superconductors interfaces. This allows us to discuss a novel antiferromagnetic proximity effect. Possible experimental realizations are discussed.
We report a quantum phase transition between orbital-selective Mott states, with different localized orbitals, in a Hunds metals model. Using the density matrix renormalization group, the phase diagram is constructed varying the electronic density and Hubbard $U$, at robust Hunds coupling. We demonstrate that this transition is preempted by charge fluctuations and the emergence of free spinless fermions, as opposed to the magnetically-driven Mott transition. The Luttinger correlation exponent is shown to have a universal value in the strong-coupling phase, whereas it is interaction dependent at intermediate couplings. At weak coupling we find a second transition from a normal metal to the intermediate-coupling phase.
Here we study the electronic properties of cuprate/manganite interfaces. By means of atomic resolution electron microscopy and spectroscopy, we produce a subnanometer scale map of the transition metal oxidation state profile across the interface between the high $T_c$ superconductor YBa$_2$Cu$_3$O$_{7-delta}$ and the colossal magnetoresistance compound (La,Ca)MnO$_3$. A net transfer of electrons from manganite to cuprate with a peculiar non-monotonic charge profile is observed. Model calculations rationalize the profile in terms of the competition between standard charge transfer tendencies (due to band mismatch), strong chemical bonding effects across the interface, and Cu substitution into the Mn lattice, with different characteristic length scales.
We investigate the real-time dynamics of photoexcited electronic instabilities in a charge-transfer system model, using the time-dependent density matrix renormalization group method. The model of choice was the quarter-filled one-dimensional extended Peierls-Hubbard Hamiltonian interacting with classical few-cycle electromagnetic radiation. The results show that only one electronic instability drives the main features of the photogenerated time-dependent behavior. Indeed, the photoresponse of the system shows a large enhancement of the $4k_F$ (bond and charge) instability whereas the $2k_F$ state remains largely unaffected. This conclusion holds regardless of the nature of the optical excitations and whether the system is perturbed resonantly or not. Our results suggest potential applications of charge-transfer systems with slow phononic dynamics as optoelectronic switching devices.
The orbital-selective Mott phase (OSMP) of multiorbital Hubbard models has been extensively analyzed before using static and dynamical mean-field approximations. In parallel, the properties of Block states (antiferromagnetically coupled ferromagnetic spin clusters) in Fe-based superconductors have also been much discussed. The present effort uses numerically exact techniques in one-dimensional systems to report the observation of Block states within the OSMP regime, connecting two seemingly independent areas of research, and providing analogies with the physics of Double-Exchange models.
Using time-dependent density-matrix renormalization group, we study the time evolution of electronic wave packets in the one-dimensional extended Hubbard model with on-site and nearest neighbor repulsion, U and V, respectively. As expected, the wave packets separate into spin-only and charge-only excitations (spin-charge separation). Charge and spin velocities exhibit non-monotonic dependence on V. For small and intermediate values of V, both velocities increase with V. However, the charge velocity exhibits a stronger dependence than that of the spin, leading to a more pronounced spin-charge separation. Charge fractionalization, on the other hand, is weakly affected by V. The results are explained in terms of Luttinger liquid theory in the weak-coupling limit, and an effective model in the strong-coupling regime.
