Quasi-universal 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 quasi-universal 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 quasi-universal behavior of the current and spectral function for different values of the hopping.Comment: (5 pages, 5 figures, ReVTeX

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\times 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 non-equilibrium setups without mass transport.Comment: 14 pages, 9 figure

Strong-coupling expansion for ultracold bosons in an optical lattice at finite temperatures in the presence of superfluidity

We develop a strong-coupling ($t \ll U$) expansion technique for calculating the density profile for bosonic atoms trapped in an optical lattice with an overall harmonic trap at finite temperature and finite on site interaction in the presence of superfluid regions. Our results match well with quantum Monte Carlo simulations at finite temperature. We also show that the superfluid order parameter never vanishes in the trap due to proximity effect. Our calculations for the scaled density in the vacuum to superfluid transition agree well with the experimental data for appropriate temperatures. We present calculations for the entropy per particle as a function of temperature which can be used to calibrate the temperature in experiments. We also discuss issues connected with the demonstration of universal quantum critical scaling in the experiments.Comment: 11 pages, 9 figure

Two-channel Kondo physics in odd impurity chains

We study odd-membered chains of spin-(1/2) impurities, with each end connected to its own metallic lead. For antiferromagnetic exchange coupling, universal two-channel Kondo (2CK) physics is shown to arise at low energies. Two overscreening mechanisms are found to occur depending on coupling strength, with distinct signatures in physical properties. For strong inter-impurity coupling, a residual chain spin-(1/2) moment experiences a renormalized effective coupling to the leads; while in the weak-coupling regime, Kondo coupling is mediated via incipient single-channel Kondo singlet formation. We also investigate models where the leads are tunnel-coupled to the impurity chain, permitting variable dot filling under applied gate voltages. Effective low-energy models for each regime of filling are derived, and for even-fillings where the chain ground state is a spin singlet, an orbital 2CK effect is found to be operative. Provided mirror symmetry is preserved, 2CK physics is shown to be wholly robust to variable dot filling; in particular the single-particle spectrum at the Fermi level, and hence the low-temperature zero-bias conductance, is always pinned to half-unitarity. We derive a Friedel-Luttinger sum rule and from it show that, in contrast to a Fermi liquid, the Luttinger integral is non-zero and determined solely by the `excess' dot charge as controlled by gate voltage. The relevance of the work to real quantum dot devices, where inter-lead charge-transfer processes fatal to 2CK physics are present, is also discussed. Physical arguments and numerical renormalization group techniques are used to obtain a detailed understanding of these problems.Comment: 21 pages, 19 figure

Supernova Ejecta in the Youngest Galactic Supernova Remnant G1.9+0.3

G1.9+0.3 is the youngest known Galactic supernova remnant (SNR), with an estimated supernova (SN) explosion date of about 1900, and most likely located near the Galactic Center. Only the outermost ejecta layers with free-expansion velocities larger than about 18,000 km/s have been shocked so far in this dynamically young, likely Type Ia SNR. A long (980 ks) Chandra observation in 2011 allowed spatially-resolved spectroscopy of heavy-element ejecta. We denoised Chandra data with the spatio-spectral method of Krishnamurthy et al., and used a wavelet-based technique to spatially localize thermal emission produced by intermediate-mass elements (IMEs: Si and S) and iron. The spatial distribution of both IMEs and Fe is extremely asymmetric, with the strongest ejecta emission in the northern rim. Fe Kalpha emission is particularly prominent there, and fits with thermal models indicate strongly oversolar Fe abundances. In a localized, outlying region in the northern rim, IMEs are less abundant than Fe, indicating that undiluted Fe-group elements (including 56Ni) with velocities larger than 18,000 km/s were ejected by this SN. But in the inner west rim, we find Si- and S-rich ejecta without any traces of Fe, so high-velocity products of O-burning were also ejected. G1.9+0.3 appears similar to energetic Type Ia SNe such as SN 2010jn where iron-group elements at such high free-expansion velocities have been recently detected. The pronounced asymmetry in the ejecta distribution and abundance inhomogeneities are best explained by a strongly asymmetric SN explosion, similar to those produced in some recent 3D delayed-detonation Type Ia models.Comment: 6 pages, 3 figures, submitted to ApJ Letter

Nonuniform Expansion of the Youngest Galactic Supernova Remnant G1.9+0.3

We report measurements of X-ray expansion of the youngest Galactic supernova remnant, G1.9+0.3, using Chandra observations in 2007, 2009, and 2011. The measured rates strongly deviate from uniform expansion, decreasing radially by about 60% along the X-ray bright SE-NW axis from 0.84% +/- 0.06% per yr to 0.52% +/- 0.03% per yr. This corresponds to undecelerated ages of 120-190 yr, confirming the young age of G1.9+0.3, and implying a significant deceleration of the blast wave. The synchrotron-dominated X-ray emission brightens at a rate of 1.9% +/- 0.4% per yr. We identify bright outer and inner rims with the blast wave and reverse shock, respectively. Sharp density gradients in either ejecta or ambient medium are required to produce the sudden deceleration of the reverse shock or the blast wave implied by the large spread in expansion ages. The blast wave could have been decelerated recently by an encounter with a modest density discontinuity in the ambient medium, such as found at a wind termination shock, requiring strong mass loss in the progenitor. Alternatively, the reverse shock might have encountered an order-of-magnitude density discontinuity within the ejecta, such as found in pulsating delayed-detonation Type Ia models. We demonstrate that the blast wave is much more decelerated than the reverse shock in these models for remnants at ages similar to G1.9+0.3. Similar effects may also be produced by dense shells possibly associated with high-velocity features in Type Ia spectra. Accounting for the asymmetry of G1.9+0.3 will require more realistic 3D Type Ia models.Comment: 6 pages, 4 figures, accepted for publication in ApJ Letters, minor revision

Equilibrium glassy phase in a polydisperse hard sphere system

The phase diagram of a polydisperse hard sphere system is examined by numerical minimization of a discretized form of the Ramakrishnan-Yussouff free energy functional. Crystalline and glassy local minima of the free energy are located and the phase diagram in the density-polydispersity plane is mapped out by comparing the free energies of different local minima. The crystalline phase disappears and the glass becomes the equilibrium phase beyond a "terminal" value of the polydispersity. A crystal to glass transition is also observed as the density is increased at high polydispersity. The phase diagram obtained in our study is qualitatively similar to that of hard spheres in a quenched random potential.Comment: 4 pages, 4 figure