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

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

Laser induced reentrant freezing in two-dimensional attractive colloidal systems

The effects of an externally applied one-dimensional periodic potential on the freezing/melting behaviour of two-dimensional systems of colloidal particles with a short-range attractive interaction are studied using Monte Carlo simulations. In such systems, incommensuration results when the periodicity of the external potential does not match the length-scale at which the minimum of the attractive potential occurs. To study the effects of this incommensuration, we consider two different models for the system. Our simulations for both these models show the phenomenon of reentrant freezing as the strength of the periodic potential is varied. Our simulations also show that different exotic phases can form when the strength of the periodic potential is high, depending on the length-scale at which the minimum of the attractive pair-potential occurs.Comment: 24 pages (including figures) in preprint forma

Correlation between Local Structure Distortions and Martensitic Transformation in Ni-Mn-In alloys

The local structural distortions arising as a consequence of increasing Mn content in Ni_2Mn_1+xIn_1-x (x=0, 0.3, 0.4, 0.5 and 0.6) and its effect on martensitic transformation have been studied using Extended X-ray Absorption Fine Structure (EXAFS) spectroscopy. Using the room temperature EXAFS at the Ni and Mn K-edges in the above compositions, the changes associated with respect to the local structure of these absorbing atoms are compared. It is seen that in the alloys exhibiting martensitic transformation ($x \ge 0.4$) there is a significant difference between the Ni-In and Ni-Mn bond lengths even in the austenitic phase indicating atomic volume to be the main factor in inducing martensitic transformation in Ni-Mn-In Heusler alloys.Comment: 8 pages, 2 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