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.

Light scattering in random media with wavelength-scale structures : astronomical and industrial applications

Light scattering, or scattering and absorption of electromagnetic waves, is an important tool in all remote-sensing observations. In astronomy, the light scattered or absorbed by a distant object can be the only source of information. In Solar-system studies, the light-scattering methods are employed when interpreting observations of atmosphereless bodies such as asteroids, atmospheres of planets, and cometary or interplanetary dust. Our Earth is constantly monitored from artificial satellites at different wavelengths. With remote sensing of Earth the light-scattering methods are not the only source of information: there is always the possibility to make in situ measurements. The satellite-based remote sensing is, however, superior in the sense of speed and coverage if only the scattered signal can be reliably interpreted. The optical properties of many industrial products play a key role in their quality. Especially for products such as paint and paper, the ability to obscure the background and to reflect light is of utmost importance. High-grade papers are evaluated based on their brightness, opacity, color, and gloss. In product development, there is a need for computer-based simulation methods that could predict the optical properties and, therefore, could be used in optimizing the quality while reducing the material costs. With paper, for instance, pilot experiments with an actual paper machine can be very time- and resource-consuming. The light-scattering methods presented in this thesis solve rigorously the interaction of light and material with wavelength-scale structures. These methods are computationally demanding, thus the speed and accuracy of the methods play a key role. Different implementations of the discrete-dipole approximation are compared in the thesis and the results provide practical guidelines in choosing a suitable code. In addition, a novel method is presented for the numerical computations of orientation-averaged light-scattering properties of a particle, and the method is compared against existing techniques. Simulation of light scattering for various targets and the possible problems arising from the finite size of the model target are discussed in the thesis. Scattering by single particles and small clusters is considered, as well as scattering in particulate media, and scattering in continuous media with porosity or surface roughness. Various techniques for modeling the scattering media are presented and the results are applied to optimizing the structure of paper. However, the same methods can be applied in light-scattering studies of Solar-system regoliths or cometary dust, or in any remote-sensing problem involving light scattering in random media with wavelength-scale structures.

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.

Trans-Scale Mechanics: Looking For The Missing Links Between Continuum And Micro/Nanoscopic Reality

Problems involving coupled multiple space and time scales offer a real challenge for conventional frameworks of either particle or continuum mechanics. In this paper, four cases studies (shear band formation in bulk metallic glasses, spallation resulting from stress wave, interaction between a probe tip and sample, the simulation of nanoindentation with molecular statistical thermodynamics) are provided to illustrate the three levels of trans-scale problems (problems due to various physical mechanisms at macro-level, problems due to micro-structural evolution at macro/micro-level, problems due to the coupling of atoms/molecules and a finite size body at micro/nano-level) and their formulations. Accordingly, non-equilibrium statistical mechanics, coupled trans-scale equations and simultaneous solutions, and trans-scale algorithms based on atomic/molecular interaction are suggested as the three possible modes of trans-scale mechanics.

Local stresses and strains in axially cracked cylindrical shells

A perturbation solution is obtained for the local stress-strain fields in an axially cracked cylindrical shell. The tenth-order differential equations are used that take into account the transverse shear deformation. The perturbation of a curvature parameter, λ, is employed, where . The stress intensity factors for finite size cylindrical shells subjected to bending and internal pressure are evaluated. Sufficient accuracy can be obtained without using fine mesh sizes in regions near the crack tip. Also analyzed are the influence of cylinder diameter and shearing stiffness on bulging.