Hydrodynamics of soft active matter

This review summarizes theoretical progress in the field of active matter, placing it in the context of recent experiments. This approach offers a unified framework for the mechanical and statistical properties of living matter: biofilaments and molecular motors in vitro or in vivo, collections of motile microorganisms, animal flocks, and chemical or mechanical imitations. A major goal of this review is to integrate several approaches proposed in the literature, from semimicroscopic to phenomenological. In particular, first considered are ``dry'' systems, defined as those where momentum is not conserved due to friction with a substrate or an embedding porous medium. The differences and similarities between two types of orientationally ordered states, the nematic and the polar, are clarified. Next, the active hydrodynamics of suspensions or ``wet'' systems is discussed and the relation with and difference from the dry case, as well as various large-scale instabilities of these nonequilibrium states of matter, are highlighted. Further highlighted are various large-scale instabilities of these nonequilibrium states of matter. Various semimicroscopic derivations of the continuum theory are discussed and connected, highlighting the unifying and generic nature of the continuum model. Throughout the review, the experimental relevance of these theories for describing bacterial swarms and suspensions, the cytoskeleton of living cells, and vibrated granular material is discussed. Promising extensions toward greater realism in specific contexts from cell biology to animal behavior are suggested, and remarks are given on some exotic active-matter analogs. Last, the outlook for a quantitative understanding of active matter, through the interplay of detailed theory with controlled experiments on simplified systems, with living or artificial constituents, is summarized.

Flow-enhanced pairing and other unusual effects in Fermi gases in synthetic gauge fields

Recent experiments on fermions in synthetic gauge fields result in systems with a spin-orbit coupling along one spatial axis, a detuning field, and a Zeeman field. We show theoretically that the presence of all three results in interesting and unusual phenomena in a system of interacting fermions (interactions described by a scattering length). For two fermions, bound states appear only over a certain range of the center-of-mass momenta. The deepest bound state appears at a nonzero center-of-mass momentum. For center-of-mass momenta without a bound state, the gauge field induces a resonance-like feature in the scattering continuum resulting in a large scattering phase shift. In the case of many particles, we demonstrate that the system, in a parameter range, shows flow-enhanced pairing, i.e., a Fulde-Farrell-Larkin-Ovchnnikov superfluid state made of robust pairs with a finite center-of-mass momentum. Yet another regime of parameters offers the opportunity to study strongly interacting normal states of spin-orbit-coupled fermionic systems utilizing the resonance-like feature induced by the synthetic gauge field.

Effect of Stress and Interface defects on Photo Luminescence of Si nano-crystals embedded in SiO2

Effect of stress and interface defects on photo luminescence property of a silicon nano-crystal (Si-nc) embedded in amorphous silicon dioxide (a-SiO2) are studied in this paper using a self-consistent quantum-continuum based modeling framework. Si-ncs or quantum dots show photoluminescence at room temperature. Whether its origin is due to Si-nc/a-SiO2 interface defects or quantum confinement of carriers in Si-nc is still an outstanding question. Earlier reports have shown that stresses greater than 12 GPa change the indirect energy band gap structure of bulk Si to a direct energy band gap structure. Such stresses are observed very often in nanostructures and these stresses influence the carrier confinement energy significantly. Hence, it is important to determine the effect of stress in addition to the structure of interface defects on photoluminescence property of Si-nc. In the present work, first a Si-nc embedded in a-SiO2 is constructed using molecular dynamics simulation framework considering the actual conditions they are grown so that the interface and residual stress in the structure evolves naturally during formation. We observe that the structure thus created has an interface of about 1 nm thick consisting of 41.95% of defective states mostly Sin+ (n = 0 to 3) coordination states. Further, both the Si-nc core and the embedding matrix are observed to be under a compressive strain. This residual strain field is applied in an effective mass k.p Hamiltonian formulation to determine the energy states of the carriers. The photo luminescence property computed based on the carrier confinement energy and interface energy states associated with defects will be analysed in details in the paper.

Reaction dynamics under confinement: an exact path integral treatment of a two-stage model of stochastic gene expression

Gene expression in living systems is inherently stochastic, and tends to produce varying numbers of proteins over repeated cycles of transcription and translation. In this paper, an expression is derived for the steady-state protein number distribution starting from a two-stage kinetic model of the gene expression process involving p proteins and r mRNAs. The derivation is based on an exact path integral evaluation of the joint distribution, P(p, r, t), of p and r at time t, which can be expressed in terms of the coupled Langevin equations for p and r that represent the two-stage model in continuum form. The steady-state distribution of p alone, P(p), is obtained from P(p, r, t) (a bivariate Gaussian) by integrating out the r degrees of freedom and taking the limit t -> infinity. P(p) is found to be proportional to the product of a Gaussian and a complementary error function. It provides a generally satisfactory fit to simulation data on the same two-stage process when the translational efficiency (a measure of intrinsic noise levels in the system) is relatively low; it is less successful as a model of the data when the translational efficiency (and noise levels) are high.

Filling-in Void and Sparse Regions in Protein Sequence Space by Protein-Like Artificial Sequences Enables Remarkable Enhancement in Remote Homology Detection Capability

Protein functional annotation relies on the identification of accurate relationships, sequence divergence being a key factor. This is especially evident when distant protein relationships are demonstrated only with three-dimensional structures. To address this challenge, we describe a computational approach to purposefully bridge gaps between related protein families through directed design of protein-like ``linker'' sequences. For this, we represented SCOP domain families, integrated with sequence homologues, as multiple profiles and performed HMM-HMM alignments between related domain families. Where convincing alignments were achieved, we applied a roulette wheel-based method to design 3,611,010 protein-like sequences corresponding to 374 SCOP folds. To analyze their ability to link proteins in homology searches, we used 3024 queries to search two databases, one containing only natural sequences and another one additionally containing designed sequences. Our results showed that augmented database searches showed up to 30% improvement in fold coverage for over 74% of the folds, with 52 folds achieving all theoretically possible connections. Although sequences could not be designed between some families, the availability of designed sequences between other families within the fold established the sequence continuum to demonstrate 373 difficult relationships. Ultimately, as a practical and realistic extension, we demonstrate that such protein-like sequences can be ``plugged-into'' routine and generic sequence database searches to empower not only remote homology detection but also fold recognition. Our richly statistically supported findings show that complementary searches in both databases will increase the effectiveness of sequence-based searches in recognizing all homologues sharing a common fold. (C) 2013 Elsevier Ltd. All rights reserved.

Local stability of a gravitating filament: a dispersion relation

Filamentary structures are ubiquitous in astrophysics and are observed at various scales. On a cosmological scale, matter is usually distributed along filaments, and filaments are also typical features of the interstellar medium. Within a cosmic filament, matter can contract and form galaxies, whereas an interstellar gas filament can clump into a series of bead-like structures that can then turn into stars. To investigate the growth of such instabilities, we derive a local dispersion relation for an idealized self-gravitating filament and study some of its properties. Our idealized picture consists of an infinite self-gravitating and rotating cylinder with pressure and density related by a polytropic equation of state. We assume no specific density distribution, treat matter as a fluid, and use hydrodynamics to derive the linearized equations that govern the local perturbations. We obtain a dispersion relation for axisymmetric perturbations and study its properties in the (kR, kz) phase space, where kR and kz are the radial and longitudinal wavenumbers, respectively. While the boundary between the stable and unstable regimes is symmetrical in kR and kz and analogous to the Jeans criterion, the most unstable mode displays an asymmetry that could constrain the shape of the structures that form within the filament. Here the results are applied to a fiducial interstellar filament, but could be extended for other astrophysical systems, such as cosmological filaments and tidal tails.

Healing time for the growth of thin films on patterned substrates

The healing times for the growth of thin films on patterned substrates are studied using simulations of two discrete models of surface growth: the Family model and the Das Sarma-Tamborenea (DT) model. The healing time, defined as the time at which the characteristics of the growing interface are ``healed'' to those obtained in growth on a flat substrate, is determined via the study of the nearest-neighbor height difference correlation function. Two different initial patterns are considered in this work: a relatively smooth tent-shaped triangular substrate and an atomically rough substrate with singlesite pillars or grooves. We find that the healing time of the Family and DT models on aL x L triangular substrate is proportional to L-z, where z is the dynamical exponent of the models. For the Family model, we also analyze theoretically, using a continuum description based on the linear Edwards-Wilkinson equation, the time evolution of the nearest-neighbor height difference correlation function in this system. The correlation functions obtained from continuum theory and simulation are found to be consistent with each other for the relatively smooth triangular substrate. For substrates with periodic and random distributions of pillars or grooves of varying size, the healing time is found to increase linearly with the height (depth) of pillars (grooves). We show explicitly that the simulation data for the Family model grown on a substrate with pillars or grooves do not agree with results of a calculation based on the continuum Edwards-Wilkinson equation. This result implies that a continuum description does not work when the initial pattern is atomically rough. The observed dependence of the healing time on the substrate size and the initial height (depth) of pillars (grooves) can be understood from the details of the diffusion rule of the atomistic model. The healing time of both models for pillars is larger than that for grooves with depth equal to the height of the pillars. The calculated healing time for both Family and DT models is found to depend on how the pillars and grooves are distributed over the substrate. (C) 2014 Elsevier B.V. All rights reserved.

Wrinkling of Atomic Planes in Ultrathin Au Nanowires

A detailed understanding of structure and stability of nanowires is critical for applications. Atomic resolution imaging of ultrathin single crystalline Au nanowires using aberration-corrected microscopy reveals an intriguing relaxation whereby the atoms in the close-packed atomic planes normal to the growth direction are displaced in the axial direction leading to wrinkling of the (111) atomic plane normal to the wire axis. First-principles calculations of the structure of such nanowires confirm this wrinkling phenomenon, whereby the close-packed planes relax to form saddle-like surfaces. Molecular dynamics studies of wires with varying diameters and different bounding surfaces point to the key role of surface stress on the relaxation process. Using continuum mechanics arguments, we show that the wrinkling arises due to anisotropy in the surface stresses and in the elastic response, along with the divergence of surface-induced bulk stress near the edges of a faceted structure. The observations provide new understanding on the equilibrium structure of nanoscale systems and could have important implications for applications in sensing and actuation.

Linewidths and line shapes in the vicinity of graphene

It is well known that graphene, by virtue of its pi-cloud delocalization, has a continuum of electronic energy states and thus behaves nearly like a metal. Instances involving quenching of electronic energy excitation in fluorophores placed in the proximity of graphene sheets are well documented. In this paper, we perform theoretical investigations on the broadening of vibrational and electronic transitions in the vicinity of graphene. We find that for CO vibrations in the vicinity of undoped graphene, the broadening at a distance of 5 angstrom is similar to 0.008 cm(-1)((kappa) over tilde = 2, (kappa) over tilde being the effective dielectric constant). In comparison, for electronic transitions, the linewidth is much larger, being of the order of several cm(-1). Also, if the transition dipole were parallel to the graphene sheet, the linewidth would be reduced to half the value for the case where it is perpendicular, an observation which should be easy to check experimentally for electronic transitions. This should be observable for the f - f transitions (which are rather narrow) of Lanthanide complexes placed within a distance of a few nanometers from a graphene sheet. Further the linewidth would have a (distance)(-4) dependence as one varies the distance from graphene. (C) 2014 AIP Publishing LLC.

VIBRATIONAL CHARACTERISTICS OF ZIGZAG, ARMCHAIR AND CHIRAL CANTILEVER SINGLE-WALLED CARBON NANOTUBES

Finite element analysis has been performed to study vibrational characteristics of cantilever single walled carbon nanotubes. Finite element models are generated by specifying the C-C bond rigidities, which are estimated by equating energies from molecular mechanics and continuum mechanics. Bending, torsion, and axial modes are identified based on effective mass for armchair, zigzag and chiral cantilever single walled carbon nanotubes, whose Young's modulus is evaluated from the bending frequency. Empirical relations are provided for frequencies of bending, torsion, and axial modes.

A novel method for studying the buckling of nanotubes considering geometrical imperfections

Buckling of nanotubes has been studied using many methods such as molecular dynamics (MD), molecular mechanics, and continuum-based shell theories. In MD, motion of the individual atoms is tracked under applied temperature and pressure, ensuring a reliable estimate of the material response. The response thus simulated varies for individual nanotubes and is only as accurate as the force field used to model the atomic interactions. On the other hand, there exists a rich literature on the understanding of continuum mechanics-based shell theories. Based on the observations on the behavior of nanotubes, there have been a number of shell theory-based approaches to study the buckling of nanotubes. Although some of these methods yield a reasonable estimate of the buckling stress, investigation and comparison of buckled mode shapes obtained from continuum analysis and MD are sparse. Previous studies show that the direct application of shell theories to study nanotube buckling often leads to erroneous results. The present study reveals that a major source of this error can be attributed to the departure of the shape of the nanotube from a perfect cylindrical shell. Analogous to the shell buckling in the macro-scale, in this work, the nanotube is modeled as a thin-shell with initial imperfection. Then, a nonlinear buckling analysis is carried out using the Riks method. It is observed that this proposed approach yields significantly improved estimate of the buckling stress and mode shapes. It is also shown that the present method can account for the variation of buckling stress as a function of the temperature considered. Hence, this can prove to be a robust method for a continuum analysis of nanosystems taking in the effect of variation of temperature as well.

Cavitation in brittle metallic glasses - Effects of stress state and distributed weak zones

Experimental studies and atomistic simulations have shown that brittle metallic glasses fail by a cavitation mechanism whose origin has been traced to the presence of intrinsic atomic density fluctuations which give rise to weak zones with reduced yield strength. It has been shown recently through continuum analysis that the presence of these zones can lower the cavitation stress considerably under equibiaxial loading. The objective of the present work is to study the effect of the applied stress state on the cavitation behavior of such a heterogeneous plastic solid with distributed weak zones. To this end, 2D plane strain finite element simulations are performed by subjecting a unit cell containing a weak zone to different (biaxiality) stress ratios. The volume fraction and yield strength of the weak zone are varied over a wide range. The results show that unlike in a homogeneous plastic solid, the cavitation stress of the heterogeneous aggregate does not reduce appreciably as the stress ratio decreases from unity when the yield strength of the weak zone is low. It is found that a non-dimensional parameter characterizing the stress state prevailing in the weak zone and its yield properties uniquely control the cavitation stress. The nature of cavitation bifurcation may change from unstable bifurcation to the left at sufficiently low stress ratio to one involving snap cavitation at high stress ratio. (C) 2014 Elsevier Ltd. All rights reserved.

Monopoles, Dirac operator, and index theory for fuzzy SU(3) / (U(1) x U(1))

The intersection of the ten-dimensional fuzzy conifold Y-F(10) with S-F(5) x S-F(5) is the compact eight-dimensional fuzzy space X-F(8). We show that X-F(8) is (the analogue of) a principal U(1) x U(1) bundle over fuzzy SU(3) / U(1) x U(1)) ( M-F(6)). We construct M-F(6) using the Gell-Mann matrices by adapting Schwinger's construction. The space M-F(6) is of relevance in higher dimensional quantum Hall effect and matrix models of D-branes. Further we show that the sections of the monopole bundle can be expressed in the basis of SU(3) eigenvectors. We construct the Dirac operator on M-F(6) from the Ginsparg-Wilson algebra on this space. Finally, we show that the index of the Dirac operator correctly reproduces the known results in the continuum.

Optimal Linear Glauber Model

Contrary to the actual nonlinear Glauber model, the linear Glauber model (LGM) is exactly solvable, although the detailed balance condition is not generally satisfied. This motivates us to address the issue of writing the transition rate () in a best possible linear form such that the mean squared error in satisfying the detailed balance condition is least. The advantage of this work is that, by studying the LGM analytically, we will be able to anticipate how the kinetic properties of an arbitrary Ising system depend on the temperature and the coupling constants. The analytical expressions for the optimal values of the parameters involved in the linear are obtained using a simple Moore-Penrose pseudoinverse matrix. This approach is quite general, in principle applicable to any system and can reproduce the exact results for one dimensional Ising system. In the continuum limit, we get a linear time-dependent Ginzburg-Landau equation from the Glauber's microscopic model of non-conservative dynamics. We analyze the critical and dynamic properties of the model, and show that most of the important results obtained in different studies can be reproduced by our new mathematical approach. We will also show in this paper that the effect of magnetic field can easily be studied within our approach; in particular, we show that the inverse of relaxation time changes quadratically with (weak) magnetic field and that the fluctuation-dissipation theorem is valid for our model.

Soft-spring wall based non-periodic boundary conditions for non-equilibrium molecular dynamics of dense fluids

Non-equilibrium molecular dynamics (MD) simulations require imposition of non-periodic boundary conditions (NPBCs) that seamlessly account for the effect of the truncated bulk region on the simulated MD region. Standard implementation of specular boundary conditions in such simulations results in spurious density and force fluctuations near the domain boundary and is therefore inappropriate for coupled atomistic-continuum calculations. In this work, we present a novel NPBC model that relies on boundary atoms attached to a simple cubic lattice with soft springs to account for interactions from particles which would have been present in an untruncated full domain treatment. We show that the proposed model suppresses the unphysical fluctuations in the density to less than 1% of the mean while simultaneously eliminating spurious oscillations in both mean and boundary forces. The model allows for an effective coupling of atomistic and continuum solvers as demonstrated through multiscale simulation of boundary driven singular flow in a cavity. The geometric flexibility of the model enables straightforward extension to nonplanar complex domains without any adverse effects on dynamic properties such as the diffusion coefficient. (c) 2015 AIP Publishing LLC.