Orbital and spin ordering phenomena in strongly correlated systems are commonly studied using the local-density approximation + dynamical mean-field theory approach. Typically, however, such simulations are restricted to simplified models (density-de
nsity Coulomb interactions, high symmetry couplings and few-band models). In this work we implement an efficient general hybridization-expansion continuous-time quantum Monte Carlo impurity solver (Krylov approach) which allows us to investigate orbital and spin ordering in a more realistic setting, including interactions that are often neglected (e.g., spin-flip and pair-hopping terms), enlarged basis sets (full d versus eg), low-symmetry distortions, and reaching the very low-temperature (experimental) regime. We use this solver to study ordering phenomena in a selection of exemplary low-symmetry transition-metal oxides: LaMnO3 and rare-earth manganites as well as the perovskites CaVO3 and YTiO3. We show that spin-flip and pair hopping terms do not affect the Kugel-Khomskii orbital-order melting transition in rare-earth manganites, or the suppression of orbital fluctuations driven by crystal field and Coulomb repulsion. For the Mott insulator YTiO3 we find a ferromagnetic transition temperature 50 K, in remarkably good agreement with experiments. For LaMnO3 we show that the classical t2g-spin approximation, commonly adopted for studying manganites, yields indeed an occupied eg orbital in very good agreement with that obtained for the full d 5-orbital Hubbard model, while the spin-spin e_g-t_{2g} correlation function calculated from the full d model is 0.74, very close to the value expected for aligned eg and t2g spins; the eg spectral function matrix is also well reproduced. Finally, we show that the t2g screening reduces the eg-eg Coulomb repulsion by about 10%
We study the mechanism of orbital-order melting observed at temperature T_OO in the series of rare-earth manganites. We find that many-body super-exchange yields a transition-temperature T_KK that decreases with decreasing rare-earth radius, and incr
eases with pressure, opposite to the experimental T_OO. We show that the tetragonal crystal-field splitting reduces T_KK further increasing the discrepancies with experiments. This proves that super-exchange effects, although very efficient, in the light of the experimentally observed trends, play a minor role for the melting of orbital ordering in rare-earth manganites.
The problem of how complex quantum systems eventually come to rest lies at the heart of statistical mechanics. The maximum entropy principle put forward in 1957 by E. T. Jaynes suggests what quantum states one should expect in equilibrium but does no
t hint as to how closed quantum many-body systems dynamically equilibrate. A number of theoretical and numerical studies accumulate evidence that under specific conditions quantum many-body models can relax to a situation that locally or with respect to certain observables appears as if the entire system had relaxed to a maximum entropy state. In this work, we report the experimental observation of the non-equilibrium dynamics of a density wave of ultracold bosonic atoms in an optical lattice in the regime of strong correlations. Using an optical superlattice, we are able to prepare the system in a well-known initial state with high fidelity. We then follow the dynamical evolution of the system in terms of quasi-local densities, currents, and coherences. Numerical studies based on the time-dependent density-matrix renormalization group method are in an excellent quantitative agreement with the experimental data. For very long times, all three local observables show a fast relaxation to equilibrium values compatible with those expected for a global maximum entropy state. We find this relaxation of the quasi-local densities and currents to initially follow a power-law with an exponent being significantly larger than for free or hardcore bosons. For intermediate times the system fulfills the promise of being a dynamical quantum simulator, in that the controlled dynamics runs for longer times than present classical algorithms based on matrix product states can efficiently keep track of.
