Nonequilibrium phase transitions in a driven sandpile model

We construct a driven sandpile slope model and study it by numerical simulations in one dimension. The model is specified by a threshold slope sigma(c), a parameter alpha, governing the local current-slope relation (beyond threshold), and j(in), the mean input current of sand. A non-equilibrium phase diagram is obtained in the alpha-j(in) plane. We find an infinity of phases, characterized by different mean slopes and separated by continuous or first-order boundaries, some of which we obtain analytically. Extensions to two dimensions are discussed.z

Velocity distribution function for a dilute granular material in shear flow

The velocity distribution function for the steady shear flow of disks (in two dimensions) and spheres (in three dimensions) in a channel is determined in the limit where the frequency of particle-wall collisions is large compared to particle-particle collisions. An asymptotic analysis is used in the small parameter epsilon, which is naL in two dimensions and na(2)L in three dimensions, where; n is the number density of particles (per unit area in two dimensions and per unit volume in three dimensions), L is the separation of the walls of the channel and a is the particle diameter. The particle-wall collisions are inelastic, and are described by simple relations which involve coefficients of restitution e(t) and e(n) in the tangential and normal directions, and both elastic and inelastic binary collisions between particles are considered. In the absence of binary collisions between particles, it is found that the particle velocities converge to two constant values (u(x), u(y)) = (+/-V, O) after repeated collisions with the wall, where u(x) and u(y) are the velocities tangential and normal to the wall, V = (1 - e(t))V-w/(1 + e(t)), and V-w and -V-w, are the tangential velocities of the walls of the channel. The effect of binary collisions is included using a self-consistent calculation, and the distribution function is determined using the condition that the net collisional flux of particles at any point in velocity space is zero at steady state. Certain approximations are made regarding the velocities of particles undergoing binary collisions :in order to obtain analytical results for the distribution function, and these approximations are justified analytically by showing that the error incurred decreases proportional to epsilon(1/2) in the limit epsilon --> 0. A numerical calculation of the mean square of the difference between the exact flux and the approximate flux confirms that the error decreases proportional to epsilon(1/2) in the limit epsilon --> 0. The moments of the velocity distribution function are evaluated, and it is found that [u(x)(2)] --> V-2, [u(y)(2)] similar to V-2 epsilon and -[u(x)u(y)] similar to V-2 epsilon log(epsilon(-1)) in the limit epsilon --> 0. It is found that the distribution function and the scaling laws for the velocity moments are similar for both two- and three-dimensional systems.

Quantum spin models with exact dimer ground states

Inspired by the exact solution of the Majumdar-Ghosh model, a family of one-dimensional, translationally invariant spin Hamiltonians is constructed. The exchange coupling in these models is antiferromagnetic, and decreases linearly with the separation between the spins. The coupling becomes identically zero beyond a certain distance. It is rigorously proved that the dimer configuration is an exact, superstable ground-state configuration of all the members of the family on a periodic chain. The ground state is twofold degenerate, and there exists an energy gap above the ground state. The Majumdar-Ghosh Hamiltonian with a twofold degenerate dimer ground state is just the first member of the family. The scheme of construction is generalized to two and three dimensions, and illustrated with the help of some concrete examples. The first member in two dimensions is the Shastry-Sutherland model. Many of these models have exponentially degenerate, exact dimer ground states.

An extended Cahn-Hilliard model for interfaces with cubic anisotropy

For studying systems with a cubic anisotropy in interfacial energy sigma, we extend the Cahn-Hilliard model by including in it a fourth-rank term, namely, gamma (ijlm) [partial derivative (2) c/(partial derivativex(i) partial derivativex(j))] [partial derivative (2) c/(partial derivativex(l) partial derivativex(m))]. This term leads to an additional linear term in the evolution equation for the composition parameter field. It also leads to an orientation-dependent effective fourth-rank coefficient gamma ([hkl]) in the governing equation for the one-dimensional composition profile across a planar interface. The main effect of a non-negative gamma ([hkl]) is to increase both sigma and interfacial width w, each of which, upon suitable scaling, is related to gamma ([hkl]) through a universal scaling function. In this model, sigma is a differentiable function of interface orientation (n) over cap, and does not exhibit cusps; therefore, the equilibrium particle shapes (Wulff shapes) do not contain planar facets. However, the anisotropy in the interfacial energy can be large enough to give rise to corners in the Wulff shapes in two dimensions. In particles of finite sizes, the corners become rounded, and their shapes tend towards the Wulff shape with increasing particle size.

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.

Phase transitions in the anchored organic bilayers of long-chain alkylammonium lead iodides (CnH2n+1NH3)(2)PbI4; n=12, 16, 18

The single perovskite slab alkylammonium lead iodides (CnH2n+1NH3)(2)PbI4, n = 12, 16, 18, display two phase transitions, just above room temperature, associated with changes in the alkylammonium chains. We have followed these two phase transitions using scanning calorimetry, X-ray powder diffraction, and IR and Raman spectroscopies. We find the first phase transition to be associated with symmetry changes arising from a dynamic rotational disordering of the ammonium headgroup of the chain whereas the second transition, the melting of the chains in two dimensions, is characterized by an increased conformational disorder of the methylene units of the alkyl chains. We examine these phase transitions in light of the interesting optical properties of these materials, as well as the relevance of these systems as models for phase transitions in lipid bilayers.

{C}^0$ penalty methods for the fully nonlinear Monge-Ampère equation

In this paper, we develop and analyze C(0) penalty methods for the fully nonlinear Monge-Ampere equation det(D(2)u) = f in two dimensions. The key idea in designing our methods is to build discretizations such that the resulting discrete linearizations are symmetric, stable, and consistent with the continuous linearization. We are then able to show the well-posedness of the penalty method as well as quasi-optimal error estimates using the Banach fixed-point theorem as our main tool. Numerical experiments are presented which support the theoretical results.

Optimality of orthogonal sequences with multi-dimensional signaling

The sum capacity on a symbol-synchronous CDMA system having processing gain N and supporting K power constrained users is achieved by employing any set of N orthogonal sequences if a few users are allowed to signal along multiple dimensions. Analogously, the minimum received power (energy-per-chip) on the symbolsynchronous CDMA system supporting K users that demand specified data rates is attained by employing any set of N orthogonal sequences. At most (N - 1) users need to be split and if there are no oversized users, these split users need to signal only in two dimensions each. These results show that sum capacity or minimum sum power can be achieved with minimal downlink signaling.

Interaction anisotropy and shear instability of aspirin polymorphs established by nanoindentation

Nanoindentation is applied to the two polymorphs of aspirin to examine and differentiate their interaction anisotropy and shear instability. Aspirin provides an excellent test system for the technique because: (i) polymorphs I and II exhibit structural similarity in two dimensions, thereby facilitating clear examination of the differences in mechanical response in relation to well-defined differences between the two crystal structures; (ii) single crystals of the metastable polymorph II have only recently become accessible; (iii) shear instability has been proposed for II. Different elastic moduli and hardness values determined for the two polymorphs are correlated with their crystal structures, and the interpretation is supported by measured thermal expansion coefficients. The stress-induced transformation of the metastable polymorph II to the stable polymorph I can be brought about rapidly by mechanical milling, and proceeds via a slip mechanism. This work establishes that nanoindentation provides ``signature'' responses for the two aspirin polymorphs, despite their very similar crystal structures. It also demonstrates the value of the technique to quantify stability relationships and phase transformations in molecular crystals, enabling a deeper understanding of polymorphism in the context of crystal engineering.

Selective double quantum resolved correlation experiment for the complete separation of entire proton NMR spectra of enantiomers

The present study reports a two dimensional NMR experiment which separates single quantum spectra of enantiomers from that of a racemic mixture. This is a blend of selective double quantum refocusing, for resolving couplings and chemical shift interactions along two dimensions followed by correlation of the selectively excited protons to the entire coupled spin network. The concept is solely based on the presence of distinct intra methyl dipolar couplings of different enantiomers when dissolved in chiral orienting media. The analysis of single enantiomer spectrum obtained from respective F-2 cross sections yield all the spectral information. (C) 2011 Elsevier Inc. All rights reserved.

Mechanical properties of ZnS nanowires and thin films: Microscopic origin of the dependence on size and growth direction

Mechanical properties of ZnS nanowires and thin films are studied as a function of size and growth direction using all-atom molecular dynamics simulations. Using the stress-strain relationship we extract Young's moduli of nanowires and thin films at room temperature. Our results show that Young's modulus of 0001] nanowires has strong size dependence. On the other hand, 01 (1) over bar0] nanowires do not exhibit a strong size dependence of Young's modulus in the size range we have investigated. We provide a microscopic understanding of this behavior on the basis of bond stretching and contraction due to the rearrangement of atoms in the surface layers. The ultimate tensile strengths of the nanowires do not show much size dependence. To investigate the mechanical behavior of ZnS in two dimensions, we calculate Young's modulus of thin films under tensile strain along the 0001] direction. Young's modulus of thin films converges to the bulk value more rapidly than that of the 0001] nanowire.

Nonequilibrium arrhythmic states and transitions in a mathematical model for diffuse fibrosis in human cardiac tissue

We present a comprehensive numerical study of spiral-and scroll-wave dynamics in a state-of-the-art mathematical model for human ventricular tissue with fiber rotation, transmural heterogeneity, myocytes, and fibroblasts. Our mathematical model introduces fibroblasts randomly, to mimic diffuse fibrosis, in the ten Tusscher-Noble-Noble-Panfilov (TNNP) model for human ventricular tissue; the passive fibroblasts in our model do not exhibit an action potential in the absence of coupling with myocytes; and we allow for a coupling between nearby myocytes and fibroblasts. Our study of a single myocyte-fibroblast (MF) composite, with a single myocyte coupled to N-f fibroblasts via a gap-junctional conductance G(gap), reveals five qualitatively different responses for this composite. Our investigations of two-dimensional domains with a random distribution of fibroblasts in a myocyte background reveal that, as the percentage P-f of fibroblasts increases, the conduction velocity of a plane wave decreases until there is conduction failure. If we consider spiral-wave dynamics in such a medium we find, in two dimensions, a variety of nonequilibrium states, temporally periodic, quasiperiodic, chaotic, and quiescent, and an intricate sequence of transitions between them; we also study the analogous sequence of transitions for three-dimensional scroll waves in a three-dimensional version of our mathematical model that includes both fiber rotation and transmural heterogeneity. We thus elucidate random-fibrosis-induced nonequilibrium transitions, which lead to conduction block for spiral waves in two dimensions and scroll waves in three dimensions. We explore possible experimental implications of our mathematical and numerical studies for plane-, spiral-, and scroll-wave dynamics in cardiac tissue with fibrosis.

Versatility of azide in serendipitous assembly of Copper(II) Magnetic Polyclusters

Engineering at the molecular level is one of the most exciting new developments for the generation of functional materials. However, the concept of designing polynuclear extended structures from bottom up is still not mature. Although progress has been made with secondary building units (SBUs) in metal organic frameworks (MOFs), the control seems to be just an illusion when it comes to bridging ligands such as the azide ion. When we say that the azido ligand is versatile in its bridging capabilities, what we mean is that it would be difficult to predict or control its bridging properties. However, this kind of serendipity is not always bad news. For example, scientists have shown that the azido ligand can mediate magnetic exchanges between paramagnetic metals in a predictable fashion (usually depending upon the bonding geometries). Therefore, it is a well-respected ligand in polynuclear assemblies. Serendipitous assemblies offer new magnetic structures that we may not otherwise even think about synthesizing. The azido ligand forms a variety of complexes with copper(II) using different blocking amines or pyridine based ligands. Its structural nature changes upon changing the substitution on amine, as well as the amount of blocking ligand. In principle, if we take any of these complexes and provide more coordination sites to the bridging azido ligands by removing a fraction of the blocking ligands, we can get new complexes with intricate structural networks and therefore different magnetic properties with the same components as used for the parent complex. In this Account, we mainly discuss the development of a number of new topological and magnetic exchange systems synthesized using this concept. Not all of these new complexes can be grouped according to their basic building structures or even by the ratio of the metal to blocking ligand. Therefore, we divided the discussion by the nuclearity of the basic building structures. Some of the complexes with the same nuclearities have very similar or even almost identical basic structures. However, the way these building units are joined together (by the azido bridges) to form the overall extended structures differ almost in every case. The complexes having the Cu-6 core are particularly interesting from a structural point of view. Although they have almost identical basic structures, some of them are extended in three dimensions, but two of them are extended in two dimensions by two different bridging networks. In the complexes having linear Cu-4 basic units, we find that using similar ligands does not always give the same bridging networks even within the basic building structures. These complexes have also enriched the field of molecular magnetism. One of the complexes with a Cu-3 building unit has provided us with the opportunity to study the competing behavior of two different kinds of magnetic exchange mechanism (ferromagnetic and antiferromagnetic) acting simultaneously between two metal ions. Through density functional theory calculations, we showed how they work independently and their additive nature to produce the overall effect. The exciting methodology for the generation of copper(II) polyclusters presented in this Account will provide the opportunity to explore analogous serendipitous assembly of diverse structures with interesting magnetic behavior using other transition metal ions having more than one unpaired electrons.

Topological Excitations in Semiconductor Heterostructures

Topological defects play an important role in the melting phenomena in two-dimensions. In this work, we report experimental observation of topological defect induced melting in two-dimensional electron systems (2DES) in the presence of strong Coulomb interaction and disorder. The phenomenon is characterised by measurement of conductivity which goes to zero in a Berezinskii-Kosterlitz-Thouless like transition. Further evidence is provided via low-frequency conductivity noise measurements.

Higher spin entanglement entropy from CFT

We consider free fermion and free boson CFTs in two dimensions, deformed by a chemical potential mu for the spin-three current. For the CFT on the infinite spatial line, we calculate the finite temperature entanglement entropy of a single interval perturbatively to second order in mu in each of the theories. We find that the result in each case is given by the same non-trivial function of temperature and interval length. Remarkably, we further obtain the same formula using a recent Wilson line proposal for the holographic entanglement entropy, in holomorphically factorized form, associated to the spin-three black hole in SL(3, R) x SL(3, R) Chern-Simons theory. Our result suggests that the order mu(2) correction to the entanglement entropy may be universal for W-algebra CFTs with spin-three chemical potential, and constitutes a check of the holographic entanglement entropy proposal for higher spin theories of gravity in AdS(3).