Ultrafast Solvation Dynamics of an Ion in the gamma-Cyclodextrin Cavity: The Role of Restricted Environment

A detailed study of the solvation dynamics of a charged coumarin dye molecule in gamma-cyclodextrin/water has been carried out by using two different theoretical approaches. The first approach is based on a multishell continuum model (MSCM). This model predicts the time scales of the dynamics rather well, provided an accurate description of the frequency-dependent dielectric function is supplied. The reason for this rather surprising agreement is 2-fold. First, there is a cancellation of errors, second, the two-zone model mimics the heterogeneous microenvironment surrounding the ion rather well. The second approach is based on the molecular hydrodynamics theory (MI-IT). In this molecular approach, the solvation dynamics has been studied by restricting the translational motion of the solvent molecules enclosed within the cavity. The results from the molecular theory are also in good agreement with the experimental results. Our study indicates that, in the present case, the restricted environment affects only the long time decay of the solvation time correlation function. The short time dynamics is still governed by the librational (and/or vibrational) modes present in bulk water.

Structure and dynamics of two-dimensional sheared granular flows

The structure and dynamics of the two-dimensional linear shear flow of inelastic disks at high area fractions are analyzed. The event-driven simulation technique is used in the hard-particle limit, where the particles interact through instantaneous collisions. The structure (relative arrangement of particles) is analyzed using the bond-orientational order parameter. It is found that the shear flow reduces the order in the system, and the order parameter in a shear flow is lower than that in a collection of elastic hard disks at equilibrium. The distribution of relative velocities between colliding particles is analyzed. The relative velocity distribution undergoes a transition from a Gaussian distribution for nearly elastic particles, to an exponential distribution at low coefficients of restitution. However, the single-particle distribution function is close to a Gaussian in the dense limit, indicating that correlations between colliding particles have a strong influence on the relative velocity distribution. This results in a much lower dissipation rate than that predicted using the molecular chaos assumption, where the velocities of colliding particles are considered to be uncorrelated.

A conditionally linearized Monte Carlo filter in non-linear structural dynamics

State and parameter estimations of non-linear dynamical systems, based on incomplete and noisy measurements, are considered using Monte Carlo simulations. Given the measurements. the proposed method obtains the marginalized posterior distribution of an appropriately chosen (ideally small) subset of the state vector using a particle filter. Samples (particles) of the marginalized states are then used to construct a family of conditionally linearized system of equations and thus obtain the posterior distribution of the states using a bank of Kalman filters. Discrete process equations for the marginalized states are derived through truncated Ito-Taylor expansions. Increased analyticity and reduced dispersion of weights computed over a smaller sample space of marginalized states are the key features of the filter that help achieve smaller sample variance of the estimates. Numerical illustrations are provided for state/parameter estimations of a Duffing oscillator and a 3-DOF non-linear oscillator. Performance of the filter in parameter estimation is also assessed using measurements obtained through experiments on simple models in the laboratory. Despite an added computational cost, the results verify that the proposed filter generally produces estimates with lower sample variance over the standard sequential importance sampling (SIS) filter.

Dynamics of dense sheared granular flows. Part 1. Structure and diffusion

Shear flows of inelastic spheres in three dimensions in the Volume fraction range 0.4-0.64 are analysed using event-driven simulations.Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to the line Joining the centres). Here, we have considered both e(t) = +1 and e(t) = e(n) (rough particles) and e(t) =-1 (smooth particles), and the normal coefficient of restitution e(n) was varied in the range 0.6-0.98. Care was taken to avoid inelastic collapse and ensure there are no particle overlaps during the simulation. First, we studied the ordering in the system by examining the icosahedral order parameter Q(6) in three dimensions and the planar order parameter q(6) in the plane perpendicular to the gradient direction. It was found that for shear flows of sufficiently large size, the system Continues to be in the random state, with Q(6) and q(6) close to 0, even for volume fractions between phi = 0.5 and phi = 0.6; in contrast, for a system of elastic particles in the absence of shear, the system orders (crystallizes) at phi = 0.49. This indicates that the shear flow prevents ordering in a system of sufficiently large size. In a shear flow of inelastic particles, the strain rate and the temperature are related through the energy balance equation, and all time scales can be non-dimensionalized by the inverse of the strain rate. Therefore, the dynamics of the system are determined only by the volume fraction and the coefficients of restitution. The variation of the collision frequency with volume fraction and coefficient of estitution was examined. It was found, by plotting the inverse of the collision frequency as a function of volume fraction, that the collision frequency at constant strain rate diverges at a volume fraction phi(ad) (volume fraction for arrested dynamics) which is lower than the random close-packing Volume fraction 0.64 in the absence of shear. The volume fraction phi(ad) decreases as the coefficient of restitution is decreased from e(n) = 1; phi(ad) has a minimum of about 0.585 for coefficient of restitution e(n) in the range 0.6-0.8 for rough particles and is slightly larger for smooth particles. It is found that the dissipation rate and all components of the stress diverge proportional to the collision frequency in the close-packing limit. The qualitative behaviour of the increase in the stress and dissipation rate are well Captured by results derived from kinetic theory, but the quantitative agreement is lacking even if the collision frequency obtained from simulations is used to calculate the pair correlation function used In the theory.

Dynamics of dense sheared granular flows. Part II. The relative velocity distributions

The distribution of relative velocities between colliding particles in shear flows of inelastic spheres is analysed in the Volume fraction range 0.4-0.64. Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to line joining the centres). The distribution or pre-collisional normal relative velocities (along the line Joining the centres of the particles) is Found to be an exponential distribution for particles with low normal coefficient of restitution in the range 0.6-0.7. This is in contrast to the Gaussian distribution for the normal relative velocity in all elastic fluid in the absence of shear. A composite distribution function, which consists of an exponential and a Gaussian component, is proposed to span the range of inelasticities considered here. In the case of roughd particles, the relative velocity tangential to the surfaces at contact is also evaluated, and it is found to be close to a Gaussian distribution even for highly inelastic particles.Empirical relations are formulated for the relative velocity distribution. These are used to calculate the collisional contributions to the pressure, shear stress and the energy dissipation rate in a shear flow. The results of the calculation were round to be in quantitative agreement with simulation results, even for low coefficients of restitution for which the predictions obtained using the Enskog approximation are in error by an order of magnitude. The results are also applied to the flow down an inclined plane, to predict the angle of repose and the variation of the volume fraction with angle of inclination. These results are also found to be in quantitative agreement with previous simulations.

Dissipative dynamics of a harmonic oscillator: A nonperturbative approach

Starting from a microscopic theory, we derive a master equation for a harmonic oscillator coupled to a bath of noninteracting oscillators. We follow a nonperturbative approach, proposed earlier by us for the free Brownian particle. The diffusion constants are calculated analytically and the positivity of the master equation is shown to hold above a critical temperature. We compare the long time behavior of the average kinetic and potential energies with known thermodynamic results. In the limit of vanishing oscillator frequency of the system, we recover the results of the free Brownian particle.

A lower bound to the survival probability and an approximate first passage time distribution for Markovian and non-Markovian dynamics in phase space

We derive a very general expression of the survival probability and the first passage time distribution for a particle executing Brownian motion in full phase space with an absorbing boundary condition at a point in the position space, which is valid irrespective of the statistical nature of the dynamics. The expression, together with the Jensen's inequality, naturally leads to a lower bound to the actual survival probability and an approximate first passage time distribution. These are expressed in terms of the position-position, velocity-velocity, and position-velocity variances. Knowledge of these variances enables one to compute a lower bound to the survival probability and consequently the first passage distribution function. As examples, we compute these for a Gaussian Markovian process and, in the case of non-Markovian process, with an exponentially decaying friction kernel and also with a power law friction kernel. Our analysis shows that the survival probability decays exponentially at the long time irrespective of the nature of the dynamics with an exponent equal to the transition state rate constant.

Diffusion of flexible, charged, nanoscopic molecules in solution: Size and pH dependence for PAMAM dendrimer

In order to understand self-diffusion (D) of a charged, flexible, and porous nanoscopic molecule in water, we carry out very long, fully atomistic molecular dynamics simulation of PAMAM dendrimer up to eight generations in explicit salt water under varying pH. We find that while the radius of gyration (R-g) varies as N-1/3, the self-diffusion constant (D) scales, surprisingly, as N-alpha, with alpha=0.39 at high pH and 0.5 at neutral pH, indicating a dramatic breakdown of Stokes-Einstein relation for diffusion of charged nanoscopic molecules. The variation in D as a function of radius of gyration demonstrates the importance of treating water and ions explicitly in the diffusion process of a flexible nanoscopic molecule. In agreement with recent experiments, the self-diffusion constant increases with pH, revealing the importance of dielectric friction in the diffusion process. The shape of a dendrimer is found to fluctuate on a nanosecond time scale. We argue that this flexibility (and also the porosity) of the dendrimer may play an important role in determining the mean square displacement of the dendrimer and the breakdown of the Stokes-Einstein relation between diffusion constant and the radius.

Dynamics after a sweep through a quantum critical point

The coherent quantum evolution of a one-dimensional many-particle system after slowly sweeping the Hamiltonian through a critical point is studied using a generalized quantum Ising model containing both integrable and nonintegrable regimes. It is known from previous work that universal power laws of the sweep rate appear in such quantities as the mean number of excitations created by the sweep. Several other phenomena are found that are not reflected by such averages: there are two different scaling behaviors of the entanglement entropy and a relaxation that is power law in time rather than exponential. The final state of evolution after the quench is not characterized by any effective temperature, and the Loschmidt echo converges algebraically for long times, with cusplike singularities in the integrable case that are dynamically broadened by nonintegrable perturbations.

Melting of an Anchored Bilayer: Molecular Dynamics Simulations of the Structural Transition in (CnH2n+1NH3)(2)PbI4 (n=12, 14, 16, 18)

The thermally driven Structural phase transition in the organic-inorganic hybrid perovskite (CnH2n+1NH3)(2)PbI4 has been investigated using molecular dynamics (MD) simulations. This system consists of positively charged alkyl-amine chains anchored to a rigid negatively charged PbI4 sheet with the chains organized as bilayers with a herringbone arrangement. Atomistic simulations were performed using ail isothermal-isobaric ensemble over a wide temperature range from 65 to 665 K for different alkyl chain lengths, n = 12, 14, 16, and 18. The simulations are able to reproduce the essential Features of the experimental observations of this system, including the existence of a transition, the linear variation of the transition temperature with alkyl chain length, and the expansion of the bilayer thickness at the transition. By use of the distance fluctuation Criteria, it is Shown that the transition is associated With a Melting of the alkyl chains of the anchored bilayer. Ail analysis of the conformation of the alkyl chains shows increased disorder in the form of gauche defects above due melting transition. Simulations also show that the melting transition is characterized by the complete disappearance of all-trans alkyl chains in the anchored bilayer, in agreement with experimental observations. A conformationally disordered chain has a larger effective cross-sectional area, and above due transition a uniformly tilted arrangement of the anchored chains call no longer be Sustained. At the melt the angular distribution of the orientation of the chains are 110 longer uniform; the chains are splayed allowing for increased space for individual chains of the anchored bilayer. This is reflected in a sharp rise in the ratio of the mean head-to-head to tail-to-tail distance of the chains of the bilayer at the transition resulting in in expansion of the bilayer thickness. The present MD simulations provide a simple explanation as to how changes in conformation of individual alkyl-chains gives rise to the observed increase in the interlayer lattice spacing of (CnH2n+1NH3)(2)PbI4 at the melting transition.

A class of one-dimensional steady-state flows in radiation-gas-dynamics

The study of steady-state flows in radiation-gas-dynamics, when radiation pressure is negligible in comparison with gas pressure, can be reduced to the study of a single first-order ordinary differential equation in particle velocity and radiation pressure. The class of steady flows, determined by the fact that the velocities in two uniform states are real, i.e. the Rankine-Hugoniot points are real, has been discussed in detail in a previous paper by one of us, when the Mach number M of the flow in one of the uniform states (at x=+∞) is greater than one and the flow direction is in the negative direction of the x-axis. In this paper we have discussed the case when M is less than or equal to one and the flow direction is still in the negative direction of the x-axis. We have drawn the various phase planes and the integral curves in each phase plane give various steady flows. We have also discussed the appearance of discontinuities in these flows.

Growth, self-assembly and dynamics of nano-scale films at fluid interfaces

Ultrathin films at fluid interfaces are important not only from a fundamental point of view as 2D complex fluids but have also become increasingly relevant in the development of novel functional materials. There has been an explosion in the synthesis work in this area over the last decade, giving rise to many exotic nanostructures at fluid interfaces. However, the factors controlling particle nucleation, growth and self-assembly at interfaces are poorly understood on a quantitative level. We will outline some of the recent attempts in this direction. Some of the selected investigations examining the macroscopic mechanical properties of molecular and particulate films at fluid interfaces will be reviewed. We conclude with a discussion of the electronic properties of these films that have potential technological and biological applications.

Effect of coordinated ions on structure and flexibility of parallel G-quadruplexes: A molecular dynamics study

ingle tract guanine residues can associate to form stable parallel quadruplex structures in the presence of certain cations. Nanosecond scale molecular dynamics simulations have been performed on fully solvated fibre model of parallel d(G(7)) quadruplex structures with Na+ or K+ ions coordinated in the cavity formed by the O6 atoms of the guanine bases. The AMBER 4.1 force field and Particle Mesh Ewald technique for electrostatic interactions have been used in all simulations. There quadruplex structures are stable during the simulation, with the middle four base tetrads showing root mean square deviation values between 0.5 to 0.8 Angstrom from the initial structure as well the high resolution crystal structure. Even in the absence of any coordinated ion in the initial structure, the G-quadruplex structure remains intact throughout the simulation. During the 1.1 ns MD simulation, one Nai counter ion from the solvent as well as several water molecules enter the central cavity to occupy the empty coordination sites within the parallel quadruplex and help stabilize the structure. Hydrogen bonding pattern depends on the nature of the coordinated ion, with the G-tetrad undergoing local structural variation to accommodate cations of different sizes. in the absence of any coordinated ion. due to strong mutual repulsion, O6 atoms within G-tetrad are forced farther apart from each other, which leads to a considerably different hydrogen bonding scheme within the G-tetrads and very favourable interaction energy between the guanine bases constituting a G-tetrad. However, a coordinated ion between G-tetrads provides extra stacking energy for the G-tetrads and makes the quadruplex structure more rigid. Na+ ions, within the quadruplex cavity, are more mobile than coordinated K+ ions. A number of hydrogen bonded water molecules are observed within the grooves of all quadruplex structures.

Monte carlo simulation for molecular gas dynamics

The dynamics of low-density flows is governed by the Boltzmann equation of the kinetic theory of gases. This is a nonlinear integro-differential equation and, in general, numerical methods must be used to obtain its solution. The present paper, after a brief review of Direct Simulation Monte Carlo (DSMC) methods due to Bird, and Belotserkovskii and Yanitskii, studies the details of theDSMC method of Deshpande for mono as well as multicomponent gases. The present method is a statistical particle-in-cell method and is based upon the Kac-Prigogine master equation which reduces to the Boltzmann equation under the hypothesis of molecular chaos. The proposed Markoff model simulating the collisions uses a Poisson distribution for the number of collisions allowed in cells into which the physical space is divided. The model is then extended to a binary mixture of gases and it is shown that it is necessary to perform the collisions in a certain sequence to obtain unbiased simulation.

A molecular dynamics study of cage-to-cage migration in sodium Y zeolite: role of surface-mediated diffusion

The details of cage-to-cage migration have been obtained from an analysis of the molecular dynamics trajectory of a probe adsorbate. It is observed that particles utilize the region within a radius of 2 angstrom from the window center but with diffusion taking place predominantly at 1.6 angstrom from the window center and a potential energy of nearly -12 kJ/mol. A barrier of about 0.5 kJ/mol is observed for surface-mediated diffusion. Surprisingly, for diffusion without surface mediation for a particle going from one cage center to another, there is an attractive well near the window instead of a barrier. At low adsorbate concentrations and room temperature, the predominant mode for cage-to-cage migration is surface-mediated diffusion. The analysis suggests that particles slide along the surface of the inner walls of the alpha-cages during migration from one cage to another.