Quasiuniversal Transient Behavior of a Nonequilibrium Mott Insulator Driven by an Electric Field

We use a self-consistent strong-coupling expansion for the self-energy (perturbation theory in the hopping) to describe the nonequilibrium dynamics of strongly correlated lattice fermions. We study the three-dimensional homogeneous Fermi-Hubbard model driven by an external electric field showing that the damping of the ensuing Bloch oscillations depends on the direction of the field and that for a broad range of field strengths a long-lived transient prethermalized state emerges. This long-lived transient regime implies that thermal equilibrium may be out of reach of the time scales accessible in present cold atom experiments but shows that an interesting new quasiuniversal transient state exists in nonequilibrium governed by a thermalized kinetic energy but not a thermalized potential energy. In addition, when the field strength is equal in magnitude to the interaction between atoms, the system undergoes a rapid thermalization, characterized by a different quasiuniversal behavior of the current and spectral function for different values of the hopping. DOI: 10.1103/PhysRevLett.109.260402

Sensory-organ-like response determines the magnetism of zigzag-edged honeycomb nanoribbons

We present an analytical effective theory for the magnetic phase diagram for zigzag-edge terminated honeycomb nanoribbons described by a Hubbard model with an interaction parameter U. We show that the edge magnetic moment varies as ln U and uncover its dependence on the width W of the ribbon. The physics of this owes its origin to the sensory-organ-like response of the nanoribbons, demonstrating that considerations beyond the usual Stoner-Landau theory are necessary to understand the magnetism of these systems. A first-order magnetic transition from an antiparallel orientation of the moments on opposite edges to a parallel orientation occurs upon doping with holes or electrons. The critical doping for this transition is shown to depend inversely on the width of the ribbon. Using variational Monte Carlo calculations, we show that magnetism is robust to fluctuations. Additionally, we show that the magnetic phase diagram is generic to zigzag-edge terminated nanostructures such as nanodots. Furthermore, we perform first-principles modeling to show how such magnetic transitions can be realized in substituted graphene nanoribbons. DOI: 10.1103/PhysRevB.87.085412

Simulation of inhomogeneous distributions of ultracold atoms in an optical lattice via a massively parallel implementation of nonequilibrium strong-coupling perturbation theory

We present a nonequilibrium strong-coupling approach to inhomogeneous systems of ultracold atoms in optical lattices. We demonstrate its application to the Mott-insulating phase of a two-dimensional Fermi-Hubbard model in the presence of a trap potential. Since the theory is formulated self-consistently, the numerical implementation relies on a massively parallel evaluation of the self-energy and the Green's function at each lattice site, employing thousands of CPUs. While the computation of the self-energy is straightforward to parallelize, the evaluation of the Green's function requires the inversion of a large sparse 10(d) x 10(d) matrix, with d > 6. As a crucial ingredient, our solution heavily relies on the smallness of the hopping as compared to the interaction strength and yields a widely scalable realization of a rapidly converging iterative algorithm which evaluates all elements of the Green's function. Results are validated by comparing with the homogeneous case via the local-density approximation. These calculations also show that the local-density approximation is valid in nonequilibrium setups without mass transport.

NaOsO3: A high Neel temperature 5d oxide

The origin of a high Neel temperature in a 5d oxide, NaOsO3, has been analyzed within the mean-field limit of a multiband Hubbard model and compared with the analogous 4d oxide, SrTcO3. Our analysis shows that there are a lot of similarities in both of these oxides on the dependence of the effective exchange interaction strength (J(0)) on the electron-electron interaction strength ( U). However, the relevant value of U in each system puts them in different portions of the parameter space. Although the Neel temperature for NaOsO3 is less than that for SrTcO3, our results suggest that there could be examples among other 5d oxides that have a higher Neel temperature. We have also examined the stability of the G-type antiferromagnetic state found in NaOsO3 as a function of electron doping within GGA + U calculations and find a robust G-type antiferromagnetic metallic state stabilized. The most surprising aspect of the doped results is the rigid bandlike evolution of the electronic structure, which indicates that the magnetism in NaOsO3 is not driven by Fermi surface nesting.

SrCU2(BO3)(2): A unique Mott Hubbard insulator

We discuss the recently discovered system SrCu2(BO3)(2), a realization of an exactly solvable model proposed two decades earlier. We propose its interpretation as a Mott Hubbard insulator. The possible superconducting phase arising from doping is explored, and its nature as well as its importance for testing the RVB theory of superconductivity are discussed.

Bose-Hubbard models in confining potentials: Inhomogeneous mean-field theory

We present an extensive study of Mott insulator (MI) and superfluid (SF) shells in Bose-Hubbard (BH) models for bosons in optical lattices with harmonic traps. For this we apply the inhomogeneous mean-field theory developed by Sheshadri et al. Phys. Rev. Lett. 75, 4075 (1995)]. Our results for the BH model with one type of spinless bosons agree quantitatively with quantum Monte Carlo simulations. Our approach is numerically less intensive than such simulations, so we are able to perform calculations on experimentally realistic, large three-dimensional systems, explore a wide range of parameter values, and make direct contact with a variety of experimental measurements. We also extend our inhomogeneous mean-field theory to study BH models with harmonic traps and (a) two species of bosons or (b) spin-1 bosons. With two species of bosons, we obtain rich phase diagrams with a variety of SF and MI phases and associated shells when we include a quadratic confining potential. For the spin-1 BH model, we show, in a representative case, that the system can display alternating shells of polar SF and MI phases, and we make interesting predictions for experiments in such systems.

Quantum Phase Diagram of One-Dimensional Spin and Hubbard Models with Transitions to Bond Order Wave Phases

Similar quantum phase diagrams and transitions are found for three classes of one-dimensional models with equally spaced sites, singlet ground states (GS), inversion symmetry at sites and a bond order wave (BOW) phase in some sectors. The models are frustrated spin-1/2 chains with variable range exchange, half-filled Hubbard models with spin-independent interactions and modified Hubbard models with site energies for describing organic charge transfer salts. In some range of parameters, the models have a first order quantum transition at which the GS expectation value of the sublattice spin < S-A(2)> of odd or even-numbered sites is discontinuous. There is an intermediate BOW phase for other model parameters that lead to two continuous quantum transitions with continuous < S-A(2)>. Exact diagonalization of finite systems and symmetry arguments provide a unified picture of familiar 1D models that have appeared separately in widely different contexts.

Quantum Phase Diagram of One-Dimensional Spin and Hubbard Models with Transitions to Bond Order Wave Phases

Similar quantum phase diagrams and transitions are found for three classes of one-dimensional models with equally spaced sites, singlet ground states (GS), inversion symmetry at sites and a bond order wave (BOW) phase in some sectors. The models are frustrated spin-1/2 chains with variable range exchange, half-filled Hubbard models with spin-independent interactions and modified Hubbard models with site energies for describing organic charge transfer salts. In some range of parameters, the models have a first order quantum transition at which the GS expectation value of the sublattice spin < S-A(2)> of odd or even-numbered sites is discontinuous. There is an intermediate BOW phase for other model parameters that lead to two continuous quantum transitions with continuous < S-A(2)>. Exact diagonalization of finite systems and symmetry arguments provide a unified picture of familiar 1D models that have appeared separately in widely different contexts.