Supersolid and solitonic phases in the one-dimensional extended Bose-Hubbard model

We report our findings on the quantum phase transitions in cold bosonic atoms in a one-dimensional optical lattice using the finite-size density-matrix renormalization-group method in the framework of the extended Bose-Hubbard model. We consider wide ranges of values for the filling factors and the nearest-neighbor interactions. At commensurate fillings, we obtain two different types of charge-density wave phases and a Mott insulator phase. However, departure from commensurate fillings yields the exotic supersolid phase where both the crystalline and the superfluid orders coexist. In addition, we obtain the signatures for the solitary waves and the superfluid phase.

Bond-order wave phase, spin solitons, and thermodynamics of a frustrated linear spin-1/2 Heisenberg antiferromagnet

The linear spin-1/2 Heisenberg antiferromagnet with exchanges J(1) and J(2) between first and second neighbors has a bond-order wave (BOW) phase that starts at the fluid-dimer transition at J(2)/J(1)=0.2411 and is particularly simple at J(2)/J(1)=1/2. The BOW phase has a doubly degenerate singlet ground state, broken inversion symmetry, and a finite-energy gap E-m to the lowest-triplet state. The interval 0.4 < J(2)/J(1) < 1.0 has large E-m and small finite-size corrections. Exact solutions are presented up to N = 28 spins with either periodic or open boundary conditions and for thermodynamics up to N = 18. The elementary excitations of the BOW phase with large E-m are topological spin-1/2 solitons that separate BOWs with opposite phase in a regular array of spins. The molar spin susceptibility chi(M)(T) is exponentially small for T << E-m and increases nearly linearly with T to a broad maximum. J(1) and J(2) spin chains approximate the magnetic properties of the BOW phase of Hubbard-type models and provide a starting point for modeling alkali-tetracyanoquinodimethane salts.

Dust Deposition on Single Cylinders

The inertial impaction of Lycopodium spores on single wires lying transverse to the direction of flow has been studied. The equations of particle motion in a potential flow field have been modified for the case when Stokes' law is inapplicable. Solutions to the above equations have been obtained by digital computation. Rec, the Reynolds number based on cylinder diameter, varied from 4 to 240; particle trajectories in a flow field at Rec = 10 have been determined for inertia parameter K = 1, 2, 4, 6, and 10. Ten trajectories were developed for the above cases by the numerical stepwise method. Experiments were performed by depositing Lycopodium spores on adhesive-coated wires of various diameters and at different velocities. The weight of dust deposited was determined with a microbalance. The experimental conditions were:. Wire diameters: 345, 457, 1500 μ. Particle diameter: 35 μ. Air velocities: 20-250 cm/sec. Inertia parameter: 1-60. The particle was considered as a point mass in the theoretical analysis. But in the experiments the ratio of particle size to wire size was not negligible (rp/rc = 0·1) and hence the effect of finite size of particle on collection efficiency due to the direct interception effect has been estimated. The effect of particle size distribution on collection efficiency has also been estimated. The experimental efficiencies obtained compare well with the calculated efficiencies at Rec = 10 when direct interception is taken into account.

Dynamical cluster approximation: Nonlocal dynamics of correlated electron systems

We recently introduced the dynamical cluster approximation (DCA), a technique that includes short-ranged dynamical correlations in addition to the local dynamics of the dynamical mean-field approximation while preserving causality. The technique is based on an iterative self-consistency scheme on a finite-size periodic cluster. The dynamical mean-field approximation (exact result) is obtained by taking the cluster to a single site (the thermodynamic limit). Here, we provide details of our method, explicitly show that it is causal, systematic, Phi derivable, and that it becomes conserving as the cluster size increases. We demonstrate the DCA by applying it to a quantum Monte Carlo and exact enumeration study of the two-dimensional Falicov-Kimball model. The resulting spectral functions preserve causality, and the spectra and the charge-density-wave transition temperature converge quickly and systematically to the thermodynamic limit as the cluster size increases.

Interfacial shear rheology of highly confined glassy polymers

We present results on interfacial shear rheology measurements on Langmuir monolayers of two different polymers, poly(vinyl acetate) and poly(methyl methacrylate) as a function of surface concentration and temperature. While for the high glass transition poly(methyl methacrylate) polymer we find a systematic transition from a viscous dominated regime to an elastic dominated regime as surface concentration is increased, monolayers of the low glass transition polymer, poly(vinyl acetate), remain viscous even at very high surface concentrations. We further interpret the results in terms of the soft glassy rheology model of Sollich et al. P. Sollich, F. C. Lequeux, P. Hebraud and M. E. Cates, Phys. Rev. Lett., 1997, 78, 2020-2023] and provide evidence of possible reduction in glass transition temperatures in both poly(methyl methacrylate) and poly(vinyl acetate) monolayers due to finite size effects.

Role of multicritical points in addressing some key problems concerning condensed matter science

We illustrate the potential of using higher order critical points in the deeper understanding of several interesting problems of condensed matter science, e.g. critical adsorption, finite size effects, morphology of critical fluctuations, reversible aggregation of colloids, dynamics of the ordering process, etc.

The collective dynamics of self-propelled particles

We propose a method for the dynamic simulation of a collection of self-propelled particles in a viscous Newtonian fluid. We restrict attention to particles whose size and velocity are small enough that the fluid motion is in the creeping flow regime. We propose a simple model for a self-propelled particle, and extended the Stokesian Dynamics method to conduct dynamic simulations of a collection of such particles. In our description, each particle is treated as a sphere with an orientation vector p, whose locomotion is driven by the action of a force dipole Sp of constant magnitude S0 at a point slightly displaced from its centre. To simplify the calculation, we place the dipole at the centre of the particle, and introduce a virtual propulsion force Fp to effect propulsion. The magnitude F0 of this force is proportional to S0. The directions of Sp and Fp are determined by p. In isolation, a self-propelled particle moves at a constant velocity u0 p, with the speed u0 determined by S0. When it coexists with many such particles, its hydrodynamic interaction with the other particles alters its velocity and, more importantly, its orientation. As a result, the motion of the particle is chaotic. Our simulations are not restricted to low particle concentration, as we implement the full hydrodynamic interactions between the particles, but we restrict the motion of particles to two dimensions to reduce computation. We have studied the statistical properties of a suspension of self-propelled particles for a range of the particle concentration, quantified by the area fraction φa. We find several interesting features in the microstructure and statistics. We find that particles tend to swim in clusters wherein they are in close proximity. Consequently, incorporating the finite size of the particles and the near-field hydrodynamic interactions is of the essence. There is a continuous process of breakage and formation of the clusters. We find that the distributions of particle velocity at low and high φa are qualitatively different; it is close to the normal distribution at high φa, in agreement with experimental measurements. The motion of the particles is diffusive at long time, and the self-diffusivity decreases with increasing φa. The pair correlation function shows a large anisotropic build-up near contact, which decays rapidly with separation. There is also an anisotropic orientation correlation near contact, which decays more slowly with separation. Movies are available with the online version of the paper.

Study of Low Temperature Magnetic Properties of a Single Chain Magnet with Alternate Isotropic and Non-collinear Anisotropic Units

Here we study thermodynamic properties of an important class of single-chain magnets (SCMs), where alternate units are isotropic and anisotropic with anisotropy axes being non-collinear. This class of SCMs shows slow relaxation at low temperatures which results from the interplay of two different relaxation mechanisms, namely dynamical and thermal. Here anisotropy is assumed to be large and negative, as a result, anisotropic units behave like canted spins at low temperatures; but even then simple Ising-type model does not capture the essential physics of the system due to quantum mechanical nature of the isotropic units. We here show how statistical behavior of this class of SCMs can be studied using a transfer matrix (TM) method. We also, for the first time, discuss in detail how weak inter-chain interactions can be treated by a TM method. The finite size effect is also discussed which becomes important for low temperature dynamics. At the end of this paper, we apply this technique to study a real helical chain magnet.

Cavitation in materials with distributed weak zones: Implications on the origin of brittle fracture in metallic glasses

Several experimental studies have shown that fracture surfaces in brittle metallic glasses (MGs) generally exhibit nanoscale corrugations which may be attributed to the nucleation and coalescence of nanovoids during crack propagation. Recent atomistic simulations suggest that this phenomenon is due to large spatial fluctuations in material properties in a brittle MG, which leads to void nucleation in regions of low atomic density and then catastrophic fracture through void coalescence. To explain this behavior, we propose a model of a heterogeneous solid containing a distribution of weak zones to represent a brittle MG. Plane strain continuum finite element analysis of cavitation in such an elastic-plastic solid is performed with the weak zones idealized as periodically distributed regions having lower yield strength than the background material. It is found that the presence of weak zones can significantly reduce the critical hydrostatic stress for the onset of cavitation which is controlled uniquely by the local yield properties of these zones. Also, the presence of weak zones diminishes the sensitivity of the cavitation stress to the volume fraction of a preexisting void. These results provide plausible explanations for the observations reported in recent atomistic simulations of brittle MGs. An analytical solution for a composite, incompressible elastic-plastic solid with a weak inner core is used to investigate the effect of volume fraction and yield strength of the core on the nature of cavitation bifurcation. It is shown that snap-cavitation may occur, giving rise to sudden formation of voids with finite size, which does not happen in a homogeneous plastic solid. (c) 2012 Elsevier Ltd. All rights reserved.

Quantum phases and phase transitions of frustrated hard-core bosons on a triangular ladder

Kinetically frustrated bosons at half filling in the presence of a competing nearest-neighbor repulsion support a wide supersolid regime on the two-dimensional triangular lattice. We study this model on a two-leg ladder using the finite-size density-matrix renormalization-group method, obtaining a phase diagram which contains three phases: a uniform superfluid (SF), an insulating charge density wave (CDW) crystal, and a bond ordered insulator (BO). We show that the transitions from SF to CDW and SF to BO are continuous in nature, with critical exponents varying continuously along the phase boundaries, while the transition from CDW to BO is found to be first order. The phase diagram is also found to contain an exactly solvable Majumdar Ghosh point, and reentrant SF to CDW phase transitions.

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.

Supersolid in a one-dimensional model of hard-core bosons

We study a system of hard-core boson on a one-dimensional lattice with frustrated next-nearest-neighbor hopping and nearest-neighbor interaction. At half filling, for equal magnitude of nearest- and next-nearest-neighbor hopping, the ground state of this system exhibits a first-order phase transition from a bond-ordered solid to a charge-density-wave solid as a function of the nearest- neighbor interaction. Moving away from half filling we investigate the system at incommensurate densities, where we find a supersolid phase which has concurrent off-diagonal long-range order and density-wave order which is unusual in a system of hard-core bosons in one dimension. Using the finite-size density-matrix renormalization group method, we obtain the complete phase diagram for this model.

Renyi entropies of free bosons on the torus and holography

We analytically evaluate the Renyi entropies for the two dimensional free boson CFT. The CFT is considered to be compactified on a circle and at finite temperature. The Renyi entropies S-n are evaluated for a single interval using the two point function of bosonic twist fields on a torus. For the case of the compact boson, the sum over the classical saddle points results in the Riemann-Siegel theta function associated with the A(n-1) lattice. We then study the Renyi entropies in the decompactification regime. We show that in the limit when the size of the interval becomes the size of the spatial circle, the entanglement entropy reduces to the thermal entropy of free bosons on a circle. We then set up a systematic high temperature expansion of the Renyi entropies and evaluate the finite size corrections for free bosons. Finally we compare these finite size corrections both for the free boson CFT and the free fermion CFT with the one-loop corrections obtained from bulk three dimensional handlebody spacetimes which have higher genus Riemann surfaces as its boundary. One-loop corrections in these geometries are entirely determined by quantum numbers of the excitations present in the bulk. This implies that the leading finite size corrections contributions from one-loop determinants of the Chern-Simons gauge field and the Dirac field in the dual geometry should reproduce that of the free boson and the free fermion CFT respectively. By evaluating these corrections both in the bulk and in the CFT explicitly we show that this expectation is indeed true.

Evolution of sex ratios in social hymenoptera: consequences of finite brood size

Evolutionarily stable sex ratios are determined for social hymenoptera under local mate competition (LMC) and when the brood size is finite. LMC is modelled by the parameter d. Of the reproductive progeny from a single foundress nest, a fraction d disperses (outbreeding), while (1-d) mate amongst themselves (sibmating). When the brood size is finite, d is taken to be the probability of an offspring dispersing, and similarly, r, the proportion of male offspring, the probability of a haploid egg being laid. Under the joint influence of these two stochastic processes, there is a nonzero probability that some females remain unmated in the nest. As a result, the optimal proportion of males (corresponding to the evolutionarily stable strategy, ESS) is higher than that obtained when the brood size is infinite. When the queen controls the sex ration, the ESS becomes more female biased under increased inbreeding (lower d), However, the ESS under worker control shows an unexpected pattern, including an increase in the proportion of males with increased inbreeding. This effect is traced to the complex interaction between inbreeding and local mate competition.

Size-Dependent Melting of Finite-Length Nanowires

We report on the size-dependent melting of nanowires with finite length based on the thermodynamic as well as liquid drop model. It has been inferred that the length dependency cannot be ignored, unlike the case of infinite length nanowires. To validate the length dependency, we have analyzed a few experimental results reported in the literature.