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.

Investigations on the melting and bending modulus of polymer grafted bilayers using dissipative particle dynamics

Understanding the influence of polymer grafted bilayers on the physicomechanical properties of lipid membranes is important while developing liposomal based drug delivery systems. The melting characteristics and bending moduli of polymer grafted bilayers are investigated using dissipative particle dynamics simulations as a function of the amount of grafted polymer and lipid tail length. Simulations are carried out using a modified Andersen barostat, whereby the membrane is maintained in a tensionless state. For lipids made up of four to six tail beads, the transition from the low temperature L-beta phase to the L-alpha phase is lowered only above a grafting fraction of G(f)=0.12 for polymers made up of 20 beads. Below G(f)=0.12 small changes are observed only for the HT4 bilayer. The bending modulus of the bilayers is obtained as a function of G(f) from a Fourier analysis of the height fluctuations. Using the theory developed by Marsh Biochim. Biophys. Acta 1615, 33 (2003)] for polymer grafted membranes, the contributions to the bending modulus due to changes arising from the grafted polymer and bilayer thinning are partitioned. The contributions to the changes in kappa from bilayer thinning were found to lie within 11% for the lipids with four to six tail beads, increasing to 15% for the lipids containing nine tail beads. The changes in the area stretch modulus were also assessed and were found to have a small influence on the overall contribution from membrane thinning. The increase in the area per head group of the lipids was found to be consistent with the scalings predicted by self-consistent mean field results. (C) 2010 American Institute of Physics.

Transition from dissipative to conservative dynamics in equations of hydrodynamics

We show, by using direct numerical simulations and theory, how, by increasing the order of dissipativity (alpha) in equations of hydrodynamics, there is a transition from a dissipative to a conservative system. This remarkable result, already conjectured for the asymptotic case alpha -> infinity U. Frisch et al., Phys. Rev. Lett. 101, 144501 (2008)], is now shown to be true for any large, but finite, value of alpha greater than a crossover value alpha(crossover). We thus provide a self-consistent picture of how dissipative systems, under certain conditions, start behaving like conservative systems and hence elucidate the subtle connection between equilibrium statistical mechanics and out-of-equilibrium turbulent flows.

Dynamics of Stick-Slip: Some Universal and Not So Universal Features

Stick-slip is usually observed in driven dissipative threshold systems. In these set of lectures, we discuss, some generic and system specific features of stickslip systems by considering a few examples wherein there has been some progress in understanding the associated dynamics. In most stick slip systems, both at low and high drive rates, the system slides smoothly, but within a window of drive rates, the motion becomes intermittent; the system alternately “sticks” till the stress builds up to a threshold value, and then “slips” when the stress is rapidly released. This intermittent motion can be traced to the existence of an unstable branch separating the two resistive branches in the force-drive-rate relation. While the two resistive branches are experimentally measurable, the unstable branch is usually not measurable and is only inferred.

Solvation dynamics in a Brownian dipolar lattice. Comparison between computer simulation and various molecular theories of solvation dynamics

Several recent theoretical and computer simulation studies have considered solvation dynamics in a Brownian dipolar lattice which provides a simple model solvent for which detailed calculations can be carried out. In this article a fully microscopic calculation of the solvation dynamics of an ion in a Brownian dipolar lattice is presented. The calculation is based on the non‐Markovian molecular hydrodynamic theory developed recently. The main assumption of the present calculation is that the two‐particle orientational correlation functions of the solid can be replaced by those of the liquid state. It is shown that such a calculation provides an excellent agreement with the computer simulation results. More importantly, the present calculations clearly demonstrate that the frequency‐dependent dielectric friction plays an important role in the long time decay of the solvation time correlation function. We also find that the present calculation provides somewhat better agreement than either the dynamic mean spherical approximation (DMSA) or the Fried–Mukamel theory which use the simulated frequency‐dependent dielectric function. It is found that the dissipative kernels used in the molecular hydrodynamic approach and in the Fried–Mukamel theory are vastly different, especially at short times. However, in spite of this disagreement, the two theories still lead to comparable results in good agreement with computer simulation, which suggests that even a semiquantitatively accurate dissipative kernel may be sufficient to obtain a reliable solvation time correlation function. A new wave vector and frequency‐dependent dissipative kernel (or memory function) is proposed which correctly goes over to the appropriate expressions in both the single particle and the collective limits. This form is expected to lead to better results than all the existing descriptions.

Ultrafast underdamped solvation: Agreement between computer simulation and various theories of solvation dynamics

A theoretical analysis of the three currently popular microscopic theories of solvation dynamics, namely, the dynamic mean spherical approximation (DMSA), the molecular hydrodynamic theory (MHT), and the memory function theory (MFT) is carried out. It is shown that in the underdamped limit of momentum relaxation, all three theories lead to nearly identical results when the translational motions of both the solute ion and the solvent molecules are neglected. In this limit, the theoretical prediction is in almost perfect agreement with the computer simulation results of solvation dynamics in the model Stockmayer liquid. However, the situation changes significantly in the presence of the translational motion of the solvent molecules. In this case, DMSA breaks down but the other two theories correctly predict the acceleration of solvation in agreement with the simulation results. We find that the translational motion of a light solute ion can play an important role in its own solvation. None of the existing theories describe this aspect. A generalization of the extended hydrodynamic theory is presented which, for the first time, includes the contribution of solute motion towards its own solvation dynamics. The extended theory gives excellent agreement with the simulations where solute motion is allowed. It is further shown that in the absence of translation, the memory function theory of Fried and Mukamel can be recovered from the hydrodynamic equations if the wave vector dependent dissipative kernel in the hydrodynamic description is replaced by its long wavelength value. We suggest a convenient memory kernel which is superior to the limiting forms used in earlier descriptions. We also present an alternate, quite general, statistical mechanical expression for the time dependent solvation energy of an ion. This expression has remarkable similarity with that for the translational dielectric friction on a moving ion.

Stabilized SPH-based simulations of impact dynamics using acceleration-corrected artificial viscosity

Artificial viscosity in SPH-based computations of impact dynamics is a numerical artifice that helps stabilize spurious oscillations near the shock fronts and requires certain user-defined parameters. Improper choice of these parameters may lead to spurious entropy generation within the discretized system and make it over-dissipative. This is of particular concern in impact mechanics problems wherein the transient structural response may depend sensitively on the transfer of momentum and kinetic energy due to impact. In order to address this difficulty, an acceleration correction algorithm was proposed in Shaw and Reid (''Heuristic acceleration correction algorithm for use in SPH computations in impact mechanics'', Comput. Methods Appl. Mech. Engrg., 198, 3962-3974) and further rationalized in Shaw et al. (An Optimally Corrected Form of Acceleration Correction Algorithm within SPH-based Simulations of Solid Mechanics, submitted to Comput. Methods Appl. Mech. Engrg). It was shown that the acceleration correction algorithm removes spurious high frequency oscillations in the computed response whilst retaining the stabilizing characteristics of the artificial viscosity in the presence of shocks and layers with sharp gradients. In this paper, we aim at gathering further insights into the acceleration correction algorithm by further exploring its application to problems related to impact dynamics. The numerical evidence in this work thus establishes that, together with the acceleration correction algorithm, SPH can be used as an accurate and efficient tool in dynamic, inelastic structural mechanics. (C) 2011 Elsevier Ltd. All rights reserved.

Self-organized criticality in dynamics without branching

We demonstrate the phenomenon of self-organized criticality (SOC) in a simple random walk model described by a random walk of a myopic ant, i.e., a walker who can see only nearest neighbors. The ant acts on the underlying lattice aiming at uniform digging, i.e., reduction of the height profile of the surface but is unaffected by the underlying lattice. In one, two, and three dimensions we have explored this model and have obtained power laws in the time intervals between consecutive events of "digging." Being a simple random walk, the power laws in space translate to power laws in time. We also study the finite size scaling of asymptotic scale invariant process as well as dynamic scaling in this system. This model differs qualitatively from the cascade models of SOC.

Shock propagation in gas dynamics: explicit form of higher order compatibility conditions

With the use of tensor analysis and the method of singular surfaces, an infinite system of equations can be derived to study the propagation of curved shocks of arbitrary strength in gas dynamics. The first three of these have been explicitly given here. This system is further reduced to one involving scalars only. The choice of dependent variables in the infinite system is quite important, it leads to coefficients free from singularities for all values of the shock strength.

Experiments on a plume with off-source heating: Implications for cloud fluid dynamics

We report here on a series of laboratory experiments on plumes, undertaken with the object of simulating the effect of the heat release that occurs in clouds on condensation of water vapor. The experimental technique used for this purpose relies on ohmic heating generated in an electrically conducting plume fluid subjected to a suitable alternating voltage across specified axial stations in the plume flow [Bhat et al., 1989]. The present series of experiments achieves a value of the Richardson number that is toward the lower end of the range that characteristics cumulus clouds. It is found that the buoyancy enhancement due to heating disrupts the eddy structures in the flow and reduces the dilution owing to entrainment of ambient fluid that would otherwise have occurred in the central region of the plume. Heating also reduces the spread rate of the plume, but as it accelerates the flow as well, the overall specific mass flux in the plume does not show a very significant change at the heat input employed in the experiment. However, there is some indication that the entrainment rate (proportional to the streamwise derivative of the mass flux) is slightly higher immediately after heat injection and slightly lower farther downstream. The measurements support a previous proposal for a cloud scenario [Bhat and Narasimha, 1996] and demonstrate how fresh insights into certain aspects of the fluid dynamics of clouds may be derived from the experimental techniques employed here.

Cell Dynamics of Model Proteins

We design rapidly folding sequences by assigning the strongest couplings to the contacts present in a target native state in a two dimensional model of heteropolymers. The pathways to folding and their dependence on the temperature are illustrated via a mapping of the dynamics into motion within the space of the maximally compact cells.

Dynamics of ribonuclease A and ribonuclease S: Computational and experimental studies

RNase S is a complex consisting of two proteolytic fragments of RNase A: the S peptide (residues 1-20) and S protein (residues 21-124). RNase S and RNase A have very similar X-ray structures and enzymatic activities. previous experiments have shown increased rates of hydrogen exchange and greater sensitivity to tryptic cleavage for RNase S relative to RNase A. It has therefore been asserted that the RNase S complex is considerably more dynamically flexible than RNase A. In the present study we examine the differences in the dynamics of RNaseS and RNase A computationally, by MD simulations, and experimentally, using trypsin cleavage as a probe of dynamics. The fluctuations around the average solution structure during the simulation were analyzed by measuring the RMS deviation in coordinates. No significant differences between RNase S and RNase A dynamics were observed in the simulations. We were able to account for the apparent discrepancy between simulation and experiment by a simple model, According to this model, the experimentally observed differences in dynamics can be quantitatively explained by the small amounts of free S peptide and S protein that are present in equilibrium with the RNase S complex. Thus, folded RNase A and the RNase S complex have identical dynamic behavior, despite the presence of a break in polypeptide chain between residues 20 and 21 in the latter molecule. This is in contrast to what has been widely believed for over 30 years about this important fragment complementation system.

Self-Absorption and Focusing Effects in Carrier Relaxation Dynamics under Depletion Conditions

The effect of a one-dimensional field (1) on the self-absorption characteristics and (2) when we have a finite numerical aperture for the objective lens that focuses the laser beam on the solid are considered here. Self-absorption, in particular its manifestation as an inner filter for the emitted signal, has been observed in luminescence experiments. Models for this effect exist and have been analyzed, but only in the absence of space charge. Using our previous results on minority carrier relaxation in the presence of a field, we obtain expressions incorporating inner filter effects. Focusing of a light beam on the sample, by an objective lens, results in a three-dimensional source and consequently a three-dimensional continuity equation to be solved for the minority carrier concentration. Assuming a one-dimensional electric field and employing Fourier-Bessel transforms, we recast the problem of carrier relaxation and solve the same via an identity that relates it to solutions obtained in the absence of focusing effects. The inner filter effect as well as focusing introduces new time scales in the problem of carrier relaxation. The interplay between the electric field and the parameters which characterize these effects and the consequent modulation of the intensity and time scales of carrier decay signals are analyzed and discussed.

Vortex dynamics at rf frequencies in Bi2Sr2CaCu2O8 single crystals

We have studied the low magnetic field high temperature region of the H-T phase diagram of Bi2Sr2CaCu2O8 single crystals using the technique of non-resonant rf response at a frequency of 20 MHz. With H(rf)parallel to a, H parallel to c, the isothermal magnetic field scans below T-c show that the frequency f(H) of the tank circuit decreases continuously with increase in H before saturating at H similar to H-D(T). Such a decrease in f(H) reflects increasing rf penetration into the weakly screened region between CuO bilayers. The saturation of f(H) at its lowest value for H similar to H-D(T) indicates complete rf penetration land hence the disappearance of field dependence) due to the vanishing of the screening rf currents I-rf(c) in those regions or equivalently when the phase coherence between adjacent superconducting layers vanishes. Therefore H,(T) represents the decoupling of the adjacent superconducting bilayers, and hence also a 3D to 2D decoupling transition of the vortex structure. Simultaneous monitoring of the field dependent rf power dissipation P(H) shows a maximum in dP/dH at H-D(T). The observed H-D(T) line in many crystals is in excellent agreement with the (l/t-1) behavior proposed for decoupling.