Grain growth and grain boundary dynamics in colloidal polycrystals

Grain boundary dynamics and grain growth play a pivotal role in the fabrication of functional polycrystalline materials. However, not much is known about the delicate interplay between various microscopic processes that drive grain boundary motion which eventually culminates in the desired grain morphology. Colloidal systems are ideally suited to bridge the gap between the microscopic and macroscopic processes underlying grain growth, since their dynamics can be followed in real space and real time with single-particle resolution. The present review aims at highlighting contributions from colloid experiments that have led to a holistic understanding of grain growth in polycrystalline materials.

Dynamics of the North American continent

The forces that cause deformation of western North America have been debated for decades. Recent studies, primarily based on analysis of crustal stresses in the western United States, have suggested that the deformation of the region is mainly controlled by gravitational potential energy (GPE) variations and boundary loads, with basal tractions due to mantle flow playing a relatively minor role. We address these issues by modelling the deviatoric stress field over western North America from a 3-D finite element mantle circulation model with lateral viscosity variations. Our approach takes into account the contribution from both topography and shallow lithosphere structure (GPE) as well as that from deeper mantle flow in one single model, as opposed to separate lithosphere and circulation models, as has been done so far. In addition to predicting the deviatoric stresses we also jointly fit the constraints of geoid, dynamic topography and plate motion both globally and over North America, in order to ensure that the forces that arise in our models are dynamically consistent. We examine the sensitivity of the dynamic models to different lateral viscosity variations. We find that circulation models that include upper mantle slabs yield a better fit to observed plate velocities. Our results indicate that a model of GPE variations coupled with mantle convection gives the best fit to the observational constraints. We argue that although GPE variations control a large part of the deformation of the western United States, deeper mantle tractions also play a significant role. The average deviatoric stress magnitudes in the western United States range 30-40 MPa. The cratonic region exhibits higher coupling to mantle flow than the rest of the continent. We find that a relatively strong San Andreas fault gives a better fit to the observational constraints, especially that of plate velocity in western North America.

Projecting low and extensive dimensional chaos from spatio-temporal dynamics

We review the spatio-temporal dynamical features of the Ananthakrishna model for the Portevin-Le Chatelier effect, a kind of plastic instability observed under constant strain rate deformation conditions. We then establish a qualitative correspondence between the spatio-temporal structures that evolve continuously in the instability domain and the nature of the irregularity of the scalar stress signal. Rest of the study is on quantifying the dynamical information contained in the stress signals about the spatio-temporal dynamics of the model. We show that at low applied strain rates, there is a one-to-one correspondence with the randomly nucleated isolated bursts of mobile dislocation density and the stress drops. We then show that the model equations are spatio-temporally chaotic by demonstrating the number of positive Lyapunov exponents and Lyapunov dimension scale with the system size at low and high strain rates. Using a modified algorithm for calculating correlation dimension density, we show that the stress-strain signals at low applied strain rates corresponding to spatially uncorrelated dislocation bands exhibit features of low dimensional chaos. This is made quantitative by demonstrating that the model equations can be approximately reduced to space independent model equations for the average dislocation densities, which is known to be low-dimensionally chaotic. However, the scaling regime for the correlation dimension shrinks with increasing applied strain rate due to increasing propensity for propagation of the dislocation bands. The stress signals in the partially propagating to fully propagating bands turn to have features of extensive chaos.

Solvation dynamics of tryptophan in water-dimethyl sulfoxide binary mixture: In search of molecular origin of composition dependent multiple anomalies

Experimental and simulation studies have uncovered at least two anomalous concentration regimes in water-dimethyl sulfoxide (DMSO) binary mixture whose precise origin has remained a subject of debate. In order to facilitate time domain experimental investigation of the dynamics of such binary mixtures, we explore strength or extent of influence of these anomalies in dipolar solvation dynamics by carrying out long molecular dynamics simulations over a wide range of DMSO concentration. The solvation time correlation function so calculated indeed displays strong composition dependent anomalies, reflected in pronounced non-exponential kinetics and non-monotonous composition dependence of the average solvation time constant. In particular, we find remarkable slow-down in the solvation dynamics around 10%-20% and 35%-50% mole percentage. We investigate microscopic origin of these two anomalies. The population distribution analyses of different structural morphology elucidate that these two slowing down are reflections of intriguing structural transformations in water-DMSO mixture. The structural transformations themselves can be explained in terms of a change in the relative coordination number of DMSO and water molecules, from 1DMSO:2H(2)O to 1H(2)O:1DMSO and 1H(2)O:2DMSO complex formation. Thus, while the emergence of first slow down (at 15% DMSO mole percentage) is due to the percolation among DMSO molecules supported by the water molecules (whose percolating network remains largely unaffected), the 2nd anomaly (centered on 40%-50%) is due to the formation of the network structure where the unit of 1DMSO:1H(2)O and 2DMSO:1H(2)O dominates to give rise to rich dynamical features. Through an analysis of partial solvation dynamics an interesting negative cross-correlation between water and DMSO is observed that makes an important contribution to relaxation at intermediate to longer times.

Relation Between the Diffusivity, Viscosity, and Ionic Radius of LiCl in Water, Methanol, and Ethylene Glycol: A Molecular Dynamics Simulation

A molecular dynamics (MD) investigation of LiCl in water, methanol, and ethylene glycol (EG) at 298 K is reported. Several; structural and dynamical properties of the ions as well as the solvent such as self-diffusivity, radial distribution functions, void and neck distributions, velocity autocorrelation functions, and mean residence times of solvent in the first solvation shell have been computed. The results show that the reciprocal relationship between the self-diffusivity of the ions and the viscosity is valid in almost all solvents with the exception of water. From an analysis of radial distribution functions and coordination numbers the nature of hydrogen bonding within the solvent and its influence on the void and neck distribution becomes evident. It is seen that the solvent solvent interaction is important in EG while solute solvent interactions dominate in water and methanol. From Voronoi tessellation, it is seen that the voids and necks within methanol are larger as compared to those within water or EG. On the basis of the void and neck distributions obtained from MD simulations and literature experimental data of limiting ion conductivity for various ions of different sizes we show that there is a relation between the void and neck radius on e one hand and dependence of conductivity on the ionic radius on the other. It is shown that the presence of large diameter voids and necks in methanol is responsible for maximum in limiting ion conductivity (lambda(0)) of TMA(+), while in water in EG, the maximum is seen for Rb+. In the case of monovalent anions, maximum in lambda(0) as a function ionic radius is seen for Br- in water EG but for the larger ClO4- ion in methanol. The relation between the void and neck distribution and the variation in lambda(0) with ionic radius arises via the Levitation effect which is discussed. These studies show the importance of the solvent structure and the associated void structure.

Dynamics of prolate spheroidal mass distributions with varying eccentricity

In this paper we calculate the potential for a prolate spheroidal distribution as in a dark matter halo with a radially varying eccentricity. This is obtained by summing up the shell-by-shell contributions of isodensity surfaces, which are taken to be concentric and with a common polar axis and with an axis ratio that varies with radius. Interestingly, the constancy of potential inside a shell is shown to be a good approximation even when the isodensity contours are dissimilar spheroids, as long as the radial variation in eccentricity is small as seen in realistic systems. We consider three cases where the isodensity contours are more prolate at large radii, or are less prolate or have a constant eccentricity. Other relevant physical quantities like the rotation velocity, the net orbital and vertical frequency due to the halo and an exponential disc of finite thickness embedded in it are obtained. We apply this to the kinematical origin of Galactic warp, and show that a prolate-shaped halo is not conducive to making long-lived warps - contrary to what has been proposed in the literature. The results for a prolate mass distribution with a variable axis ratio obtained are general, and can be applied to other astrophysical systems, such as prolate bars, for a more realistic treatment.

Dynamics of a flexible splitter plate in the wake of a circular cylinder

Rigid splitter plates in the wake of bluff bodies are known to suppress the primary vortex shedding. In the present work, we experimentally study the problem of a flexible splitter plate in the wake of a circular cylinder. In this case, the splitter plate is free to continuously deform along its length due to the fluid forces acting on it; the flexural rigidity (EI) of the plate being an important parameter. Direct visualizations of the splitter plate motions, for very low values of flexural rigidity (EI), indicate periodic traveling wave type deformations of the splitter plate with maximum tip amplitudes of the order of I cylinder diameter. As the Reynolds number based on cylinder diameter is varied, two regimes of periodic splitter plate motions are found that are referred to as mode I and mode II, with a regime of aperiodic motions between them. The frequency of plate motions in both periodic modes is found to be close to the plane cylinder Strouhal number of about 0.2, while the average frequencies in the non-periodic regime are substantially lower. The measured normalized phase speed of the traveling wave for both periodic modes is also close to the convection speed of vortices in the plane cylinder wake. As the flexural rigidity of the plate (EI) is increased, the response of the plate was found to shift to the right when plotted with flow speed or Re. To better capture the effect of varying EI, we define and use a non-dimensional bending stiffness, K*, similar to the ones used in the flag flutter problem, K*=EI/(0.5 rho(UL3)-L-2), where U is the free-stream velocity and L is the splitter plate length. Amplitude data for different EI cases when plotted against this parameter appear to collapse on to a single curve for a given splitter plate length. Measurements of the splitter plate motions for varying splitter plate lengths indicate that plates that are substantially larger than the formation length of the plane cylinder wake have similar responses, while shorter plates show significant differences.

Spiral-Wave Dynamics in a Mathematical Model of Human Ventricular Tissue with Myocytes and Fibroblasts

Cardiac fibroblasts, when coupled functionally with myocytes, can modulate the electrophysiological properties of cardiac tissue. We present systematic numerical studies of such modulation of electrophysiological properties in mathematical models for (a) single myocyte-fibroblast (MF) units and (b) two-dimensional (2D) arrays of such units; our models build on earlier ones and allow for zero-, one-, and two-sided MF couplings. Our studies of MF units elucidate the dependence of the action-potential (AP) morphology on parameters such as E-f, the fibroblast resting-membrane potential, the fibroblast conductance G(f), and the MF gap-junctional coupling G(gap). Furthermore, we find that our MF composite can show autorhythmic and oscillatory behaviors in addition to an excitable response. Our 2D studies use (a) both homogeneous and inhomogeneous distributions of fibroblasts, (b) various ranges for parameters such as G(gap), G(f), and E-f, and (c) intercellular couplings that can be zero-sided, one-sided, and two-sided connections of fibroblasts with myocytes. We show, in particular, that the plane-wave conduction velocity CV decreases as a function of G(gap), for zero-sided and one-sided couplings; however, for two-sided coupling, CV decreases initially and then increases as a function of G(gap), and, eventually, we observe that conduction failure occurs for low values of G(gap). In our homogeneous studies, we find that the rotation speed and stability of a spiral wave can be controlled either by controlling G(gap) or E-f. Our studies with fibroblast inhomogeneities show that a spiral wave can get anchored to a local fibroblast inhomogeneity. We also study the efficacy of a low-amplitude control scheme, which has been suggested for the control of spiral-wave turbulence in mathematical models for cardiac tissue, in our MF model both with and without heterogeneities.

Dynamics of edge dislocations in a sheared lamellar mesophase

The dynamics and interactions of edge dislocations in a nearly aligned sheared lamellar mesophase is analysed to provide insights into the relationship between disorder and rheology. First, the mesoscale permeation and momentum equations for the displacement field in the presence of external forces are derived from the model H equations for the concentration and momentum field. The secondary flow generated due to the mean shear around an isolated defect is calculated, and the excess viscosity due to the presence of the defect is determined from the excess energy dissipation due to the secondary flow. The excess viscosity for an isolated defect is found to increase with system size in the cross-stream direction as L-3/2 for an isolated defect, though this divergence is cut-off due to interactions in a defect suspension. As the defects are sheared past each other due to the mean flow, the Peach-Koehler force due to elastic interaction between pairs of defects is found to cause no net displacement relative to each other as they approach from large separation to the distance of closest approach. The equivalent force due to viscous interactions is found to increase the separation for defects of opposite sign, and decrease the separation for defects of same sign. During defect interactions, we find that there is no buckling instability due to dilation of layers for systems of realistic size. However, there is another mechanism, which is the velocity difference generated across a slightly deformed bilayer due to the mean shear, which could result in the creation of new defects. (C) 2013 AIP Publishing LLC.

Layerwise decomposition of water dynamics in reverse micelles: A simulation study of two-dimensional infrared spectrum

We present computer simulation study of two-dimensional infrared spectroscopy (2D-IR) of water confined in reverse micelles (RMs) of various sizes. The present study is motivated by the need to understand the altered dynamics of confined water by performing layerwise decomposition of water, with an aim to quantify the relative contributions of different layers water molecules to the calculated 2D-IR spectrum. The 0-1 transition spectra clearly show substantial elongation, due to in-homogeneous broadening and incomplete spectral diffusion, along the diagonal in the surface water layer of different sized RMs. Fitting of the frequency fluctuation correlation functions reveal that the motion of the surface water molecules is sub-diffusive and indicate the constrained nature of their dynamics. This is further supported by two peak nature of the angular analogue of van Hove correlation function. With increasing system size, the water molecules become more diffusive in nature and spectral diffusion almost completes in the central layer of the larger size RMs. Comparisons between experiments and simulations establish the correspondence between the spectral decomposition available in experiments with the spatial decomposition available in simulations. Simulations also allow a quantitative exploration of the relative role of water, sodium ions, and sulfonate head groups in vibrational dephasing. Interestingly, the negative cross correlation between force on oxygen and hydrogen of O-H bond in bulk water significantly decreases in the surface layer of each RM. This negative cross correlation gradually increases in the central water pool with increasing RMs size and this is found to be partly responsible for the faster relaxation rate of water in the central pool. (C) 2013 AIP Publishing LLC.

Switching and release dynamics of an electrostatically actuated MEMS switch under the influence of squeeze-film damping

This paper reports analytical modeling, simulation and experimental validation for switching and release times of an electrostatically actuated micromachined switch. Presented work is an extension of our earlier work [1] that analytically argued, and numerically and experimentally demonstrated, why pull-in time is larger that pull-up time when the actuation voltage is less than twice of the pull-in voltage. In this paper, switching dynamics is investigated under the influence of squeeze-film damping. Tests were performed on SOI (silicon-on-insulator) based parallel beams structures.Typical voltage requirement for actuation is in the range of 10-30 V. All the experiments were performed in normal atmospheric pressure. Measurement results confirm that the quality factor Q has appreciable effect on the release time compared to the switching time. The quality factor Q is extracted from the response measurement and compared with the ANSYS simulation result. In addition, the dynamic pull-in effect has also been studied and reported in this paper. A contribution of this work includes the effect of various phenomena such as squeeze-film damping, dynamic pull-in, and frequency pull-in effects on the switching dynamics of a MEMS switch.

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.

Capturability of augmented proportional navigation (APN) guidance with nonlinear engagement dynamics

Proportional Navigation (PN) and its variants are widely used guidance philosophies. However, in the presence of target maneuver, PN guidance law is effective only for a restrictive set of initial geometries. To account for target maneuvers, the concept of Augmented Proportional Navigation (APN) guidance law was introduced and analyzed in a linearized interceptor-target engagement framework presented in literature. However, there is no work in the literature, that addresses the capturability performance of the APN guidance law in a nonlinear engagement framework. This paper presents such an analysis and obtains the conditions for capturability. It also shows that a shorter time of interception is obtained when APN is formulated in the nonlinear framework as proposed in this paper. Simulation results are given to support the theoretical findings.