Hard-core bosons in a zig-zag optical superlattice

We study a system of hard-core bosons at half-filling in a one-dimensional optical superlattice. The bosons are allowed to hop to nearest-and next-nearest-neighbor sites. We obtain the ground-state phase diagram as a function of microscopic parameters using the finite-size density-matrix renormalization-group method. Depending on the sign of the next-nearest-neighbor hopping and the strength of the superlattice potential the system exhibits three different phases, namely the bond-order (BO) solid, the superlattice induced Mott insulator (SLMI), and the superfluid (SF) phase. When the signs of both hopping amplitudes are the same (the unfrustratedase), the system undergoes a transition from the SF to the SLMI at a nonzero value of the superlattice potential. On the other hand, when the two amplitudes differ in sign (the frustrated case), the SF is unstable to switching on a superlattice potential and also exists only up to a finite value of the next-nearest-neighbor hopping. This part of the phase diagram is dominated by the BO phase which breaks translation symmetry spontaneously even in the absence of the superlattice potential and can thus be characterized by a bond-order parameter. The transition from BO to SLMI appears to be first order.

Correlated Conductance Fluctuations Close to the Berezinskii-Kosterlitz-Thouless Transition in Ultrathin NbN Films

We probe the presence of long-range correlations in phase fluctuations by analyzing the higher-order spectrum of resistance fluctuations in ultrathin NbN superconducting films. The non-Gaussian component of resistance fluctuations is found to be sensitive to film thickness close to the transition, which allows us to distinguish between mean field and Berezinskii-Kosterlitz-Thouless (BKT) type superconducting transitions. The extent of non-Gaussianity was found to be bounded by the BKT and mean field transition temperatures and depends strongly on the roughness and structural inhomogeneity of the superconducting films. Our experiment outlines a novel fluctuation-based kinetic probe in detecting the nature of superconductivity in disordered low-dimensional materials.

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

We develop a strong-coupling (t << 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 the 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.

Local structural disorder and its influence on the average global structure and polar properties in Na0.5Bi0.5TiO3

Na0.5Bi0.5TiO3 (NBT) and its derivatives have prompted a great surge in interest owing to their potential as lead-free piezoelectrics. In spite of five decades since its discovery, there is still a lack of clarity on crucial issues such as the origin of significant dielectric relaxation at room temperature, structural factors influencing its depoling, and the status of the recently proposed monoclinic (Cc) structure vis-a-vis the nanosized structural heterogeneities. In this work, these issues are resolved by comparative analysis of local and global structures on poled and unpoled NBT specimens using electron, x-ray, and neutron diffraction in conjunction with first-principles calculation, dielectric, ferroelectric, and piezoelectric measurements. The reported global monoclinic (Cc) distortion is shown not to correspond to the thermodynamic equilibrium state at room temperature. The global monocliniclike appearance rather owes its origin to the presence of local structural and strain heterogeneities. Poling removes the structural inhomogeneities and establishes a long-range rhombohedral distortion. In the process the system gets irreversibly transformed from a nonergodic relaxor to a normal ferroelectric state. The thermal depoling is shown to be associated with the onset of incompatible in-phase tilted octahedral regions in the field-stabilized long range rhombohedral distortion.

Observation of Enhanced Diffusivity in Magnetically Powered Reciprocal Swimmers

We report on the development of a system of micron-sized reciprocal swimmers that can be powered with small homogeneous magnetic fields, and whose motion resembles that of a helical flagellum moving back and forth. We have measured the diffusivities of the swimmers to be higher compared to nonactuated objects of identical dimensions at long time scales, in accordance with the theoretical predictions made by Lauga Phys. Rev. Lett. 106, 178101 (2011)]. Randomness in the reciprocity of the actuation strokes was found to have a strong influence on the enhancement of the diffusivity, which has been investigated with numerical calculations.

Effects of initial height on the steady-state persistence probability of linear growth models

The effects of the initial height on the temporal persistence probability of steady-state height fluctuations in up-down symmetric linear models of surface growth are investigated. We study the (1 + 1)-dimensional Family model and the (1 + 1)-and (2 + 1)-dimensional larger curvature (LC) model. Both the Family and LC models have up-down symmetry, so the positive and negative persistence probabilities in the steady state, averaged over all values of the initial height h(0), are equal to each other. However, these two probabilities are not equal if one considers a fixed nonzero value of h(0). Plots of the positive persistence probability for negative initial height versus time exhibit power-law behavior if the magnitude of the initial height is larger than the interface width at saturation. By symmetry, the negative persistence probability for positive initial height also exhibits the same behavior. The persistence exponent that describes this power-law decay decreases as the magnitude of the initial height is increased. The dependence of the persistence probability on the initial height, the system size, and the discrete sampling time is found to exhibit scaling behavior.

Transverse single spin asymmetry in e+p↑→e+J/ψ+X and transverse momentum dependent evolution of the Sivers function

We extend our analysis of transverse single spin asymmetry in electroproduction of J/ψ to include the effect of the scale evolution of the transverse momentum dependent (TMD) parton distribution functions and gluon Sivers function. We estimate single spin asymmetry for JLab, HERMES, COMPASS, and eRHIC energies using the color evaporation model of charmonium production, using an analytically obtained approximate solution of TMD evolution equations discussed in the literature. We find that there is a reduction in the asymmetry compared with our predictions for the earlier case considered by us, wherein the Q2 dependence came only from DGLAP evolution of the unpolarized gluon densities and a different parametrization of the TMD Sivers function was used.

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.

Elliptical tracers in two-dimensional, homogeneous, isotropic fluid turbulence: The statistics of alignment, rotation, and nematic order

We study the statistical properties of orientation and rotation dynamics of elliptical tracer particles in two-dimensional, homogeneous, and isotropic turbulence by direct numerical simulations. We consider both the cases in which the turbulent flow is generated by forcing at large and intermediate length scales. We show that the two cases are qualitatively different. For large-scale forcing, the spatial distribution of particle orientations forms large-scale structures, which are absent for intermediate-scale forcing. The alignment with the local directions of the flow is much weaker in the latter case than in the former. For intermediate-scale forcing, the statistics of rotation rates depends weakly on the Reynolds number and on the aspect ratio of particles. In contrast with what is observed in three-dimensional turbulence, in two dimensions the mean-square rotation rate increases as the aspect ratio increases.

First-principles DFT plus GW study of oxygen vacancies in rutile TiO2

We perform first-principles calculations of the quasiparticle defect states, charge transition levels, and formation energies of oxygen vacancies in rutile titanium dioxide. The calculations are done within the recently developed combined DFT + GW formalism, including the necessary electrostatic corrections for the supercells with charged defects. We find the oxygen vacancy to be a negative U defect, where U is the defect electron addition energy. For Fermi level values below similar to 2.8 eV (relative to the valence-band maximum), we find the +2 charge state of the vacancy to be the most stable, while above 2.8 eV we find that the neutral charge state is the most stable.

Higher spin cosmology

We construct cosmological solutions of higher spin gravity in 2 + 1 dimensional de Sitter space. We show that a consistent thermodynamics can be obtained for their horizons by demanding appropriate holonomy conditions. This is equivalent to demanding the integrability of the Euclidean boundary conformal field theory partition function, and it reduces to Gibbons-Hawking thermodynamics in the spin-2 case. By using the prescription of Maldacena, we relate the thermodynamics of these solutions to those of higher spin black holes in AdS(3).

Large deviation function and fluctuation theorem for classical particle transport

We analytically evaluate the large deviation function in a simple model of classical particle transfer between two reservoirs. We illustrate how the asymptotic long-time regime is reached starting from a special propagating initial condition. We show that the steady-state fluctuation theorem holds provided that the distribution of the particle number decays faster than an exponential, implying analyticity of the generating function and a discrete spectrum for its evolution operator.

Finite-temperature dynamics of vortices in Bose-Einstein condensates

We study the dynamics of a single vortex and a pair of vortices in quasi two-dimensional Bose-Einstein condensates at finite temperatures. To this end, we use the stochastic Gross-Pitaevskii equation, which is the Langevin equation for the Bose-Einstein condensate. For a pair of vortices, we study the dynamics of both the vortex-vortex and vortex-antivortex pairs, which are generated by rotating the trap and moving the Gaussian obstacle potential, respectively. Due to thermal fluctuations, the constituent vortices are not symmetrically generated with respect to each other at finite temperatures. This initial asymmetry coupled with the presence of random thermal fluctuations in the system can lead to different decay rates for the component vortices of the pair, especially in the case of two corotating vortices.

Structure-function hierarchies and von Karman-Howarth relations for turbulence in magnetohydrodynamical equations

We generalize the method of A. M. Polyakov, Phys. Rev. E 52, 6183 (1995)] for obtaining structure-function relations in turbulence in the stochastically forced Burgers equation, to develop structure-function hierarchies for turbulence in three models for magnetohydrodynamics (MHD). These are the Burgers analogs of MHD in one dimension Eur. Phys. J.B 9, 725 (1999)], and in three dimensions (3DMHD and 3D Hall MHD). Our study provides a convenient and unified scheme for the development of structure-function hierarchies for turbulence in a variety of coupled hydrodynamical equations. For turbulence in the three sets of MHD equations mentioned above, we obtain exact relations for third-order structure functions and their derivatives; these expressions are the analogs of the von Karman-Howarth relations for fluid turbulence. We compare our work with earlier studies of such relations in 3DMHD and 3D Hall MHD.

Free turbulent shear layer in a point vortex gas as a problem in nonequilibrium statistical mechanics

This paper attempts to unravel any relations that may exist between turbulent shear flows and statistical mechanics through a detailed numerical investigation in the simplest case where both can be well defined. The flow considered for the purpose is the two-dimensional (2D) temporal free shear layer with a velocity difference Delta U across it, statistically homogeneous in the streamwise direction (x) and evolving from a plane vortex sheet in the direction normal to it (y) in a periodic-in-x domain L x +/-infinity. Extensive computer simulations of the flow are carried out through appropriate initial-value problems for a ``vortex gas'' comprising N point vortices of the same strength (gamma = L Delta U/N) and sign. Such a vortex gas is known to provide weak solutions of the Euler equation. More than ten different initial-condition classes are investigated using simulations involving up to 32 000 vortices, with ensemble averages evaluated over up to 10(3) realizations and integration over 10(4)L/Delta U. The temporal evolution of such a system is found to exhibit three distinct regimes. In Regime I the evolution is strongly influenced by the initial condition, sometimes lasting a significant fraction of L/Delta U. Regime III is a long-time domain-dependent evolution towards a statistically stationary state, via ``violent'' and ``slow'' relaxations P.-H. Chavanis, Physica A 391, 3657 (2012)], over flow time scales of order 10(2) and 10(4)L/Delta U, respectively (for N = 400). The final state involves a single structure that stochastically samples the domain, possibly constituting a ``relative equilibrium.'' The vortex distribution within the structure follows a nonisotropic truncated form of the Lundgren-Pointin (L-P) equilibrium distribution (with negatively high temperatures; L-P parameter lambda close to -1). The central finding is that, in the intermediate Regime II, the spreading rate of the layer is universal over the wide range of cases considered here. The value (in terms of momentum thickness) is 0.0166 +/- 0.0002 times Delta U. Regime II, extensively studied in the turbulent shear flow literature as a self-similar ``equilibrium'' state, is, however, a part of the rapid nonequilibrium evolution of the vortex-gas system, which we term ``explosive'' as it lasts less than one L/Delta U. Regime II also exhibits significant values of N-independent two-vortex correlations, indicating that current kinetic theories that neglect correlations or consider them as O(1/N) cannot describe this regime. The evolution of the layer thickness in present simulations in Regimes I and II agree with the experimental observations of spatially evolving (3D Navier-Stokes) shear layers. Further, the vorticity-stream-function relations in Regime III are close to those computed in 2D Navier-Stokes temporal shear layers J. Sommeria, C. Staquet, and R. Robert, J. Fluid Mech. 233, 661 (1991)]. These findings suggest the dominance of what may be called the Kelvin-Biot-Savart mechanism in determining the growth of the free shear layer through large-scale momentum and vorticity dispersal.