Effect of Alkyl Chain Arrangement on Conformation and Dynamics in a Surfactant Intercalated Layered Double Hydroxide: Spectroscopic Measurements and MD Simulations

The anionic surfactant dodecyl sulfate (DDS) has been intercalated in an Mg-Al layered double hydroxide (LDH). Monolayer and bilayer arrangements of the alkyl chains of the intercalated surfactant can be engineered by tuning the Al/Mg ratio of the LDH. In both arrangements the anionic headgroup of the surfactant is tethered to the LDH sheets, and consequently translational mobility of the chains is absent. The degrees of freedom of the confined alkyl chains are restricted to changes in conformation. The effects of the arrangement of the intercalated surfactant chains on conformational order and dynamics have been,investigated by spectroscopic measurements and molecular dynamics simulations. Infrared, Raman, and C-13 NMR spectroscopies were used to investigate conformation of the alkyl chains in the monolayer and bilayer arrangements and variable contact time cross-polarization magic angle spinning (VCT CPMAS) NMR measurements to probe molecular motion. The alkyl chains in the monolayer arrangement of the intercalated DDS chains showed considerably greater conformational disorder and faster dynamics as compared to chains in the bilayer arrangement, in spite of the fact that the volume available per chain in the monolayer is smaller than that in the bilayer. Atomistic MD simulations of the two arrangements of the intercalated surfactant were carried out using an isothermal-isobaric ensemble. The simulations are able to reproduce the essential results of the experiment-greater conformational disorder and faster dynamics for the alkyl chains in the monolayer arrangement of the intercalated surfactant. The MD simulations show that these results are a consequence of the fact that the nature of conformational disorder in the two arrangements is different. In the monolayer arrangement the alkyl chains can sustain isolated gauche defects, whereas in the bilayer arrangement gauche conformers occur only as part of a kink a gauche(+) trans gauche(-) sequence.

Variability of sif and cod of stochastic structural systems

Stochastic structural systems having a stochastic distribution of material properties and stochastic external loadings in space are analysed when a crack of deterministic size is present. The material properties and external loadings are considered to constitute independent, two-dimensional, univariate, real, homogeneous stochastic fields. The stochastic fields are characterized by their means, variances, autocorrelation functions or the equivalent power spectral density functions, and scale fluctuations. The Young's modulus and Poisson's ratio are treated to be stochastic quantities. The external loading is treated to be a stochastic field in space. The energy release rate is derived using the method of virtual crack extension. The deterministic relationship is derived to represent the sensitivities of energy release rate with respect to both virtual crack extension and real system parameter fluctuations. Taylor series expansion is used and truncation is made to the first order. This leads to the determination of second-order properties of the output quantities to the first order. Using the linear perturbations about the mean values of the output quantities, the statistical information about the energy release rates, SIF and crack opening displacements are obtained. Both plane stress and plane strain cases are considered. The general expressions for the SIF in all the three fracture modes are derived and a more detailed analysis is conducted for a mode I situation. A numerical example is given.

Dynamics of Fractional Brownian Walks

We investigate the dynamics of polymers whose solution configurations are represented by fractional Brownian walks. The calculation of the two dynamical quantities considered here, the longest relaxation time tau(r) and the intrinsic viscosity [eta], is formulated in terms of Langevin equations and is carried out within the continuum approach developed in an earlier paper. Our results for tau(r) and [eta] reproduce known scaling relations and provide reasonable numerical estimates of scaling amplitudes. The possible relevance of the work to the study of globular proteins and other compact polymeric phases is discussed.

A two timescale stochastic approximation scheme for simulation-based parametric optimization

A two timescale stochastic approximation scheme which uses coupled iterations is used for simulation-based parametric optimization as an alternative to traditional "infinitesimal perturbation analysis" schemes, It avoids the aggregation of data present in many other schemes. Its convergence is analyzed, and a queueing example is presented.

Molecular-dynamics studies of systems of confined dumbbell molecules

We present the results of molecular-dynamics simulations of systems of dumbbell molecules confined by parallel molecular walls. We have carried out systematic studies of three cases: freezing, steady flows, and stick-slip friction. We find that the molecular orientational degrees of freedom cause the surface layers to deviate from a planar configuration. Nevertheless, steady flows, in a channel as narrow as 15 molecular sizes, display continuum behavior. A range of mechanisms in the dynamics of the freezing of a confined fluid is found, as a function of the wall-fluid interactions and the bond length of the dumbbell molecules. The simple order-disorder transition associated with stick-slip motion in the presence of a layer of monoatomic lubricant molecules is supplanted by more complex behavior due to rotational degrees of freedom of the diatomic molecules.

Ultrafast solvation dynamics in water: Isotope effects and comparison with experimental results

A detailed theoretical study of solvation dynamics in water is presented. The motivation of the present study comes from the recent experimental observation that the dynamics of solvation of an ion in water is ultrafast and the solvation time correlation function decays with a time constant of about 55 fs. The slower decay in the long time can be described by a sum of two exponentials with time constants equal to 126 and 880 fs. The molecular theory (developed earlier) predicts a time constant equal to 52 fs for the initial Gaussian decay and time constants equal to 134 and 886 fs for the two exponential components at the long time. This nearly perfect agreement is obtained by using the most detailed dynamical information available in the literature. The present study emphasizes the importance of the intermolecular vibrational band originating from the O...O stretching mode of the O�H...O units in the initial dynamics and raises several interesting questions regarding the nature of the decay of this mode. We have also studied the effects of isotope substitution on solvation dynamics. It is predicted that a significant isotope effect may be observed in the long time. The experimental results have also been compared with the prediction of the dynamic mean spherical approximation (DMSA); the agreement is not satisfactory at the long time. It is further found that the molecular theory and the DMSA lead to virtually identical results if the translational modes of the solvent molecules are neglected in the former. DMSA has also been used to investigate the dynamics of solvation of a dipolar solute in water. It is found that the dynamics of dipolar solvation exhibit features rather different from those of ion solvation. © 1995 American Institute of Physics.

Ultrafast Solvation Dynamics Of An Ion In A Restricted Environment

The solvation dynamics of an excited coumarin dye molecule (C-480) enclosed within a restricted space have been studied using molecular hydrodynamic theory (MHT) and compared with the recent experimental findings. The solvation dynamics of the dye molecule within the cavity of a toroidal gamma-cyclodextrin molecule have been shown to be explained only in terms of the freezing of the solvent translational modes using MHT. The results of the theoretical calculation are in good agreement with the experimental results. The inertial components of the solvation time correlation function remain the same in both the restricted environment and in the free space. These results are interesting in the light of the simulation studies of Maroncelli and Fleming [J chem Phys, 89 (1988) 5044] which concludes that the participation of the different solvation shells in controlling the dynamics are much different. The earlier studies have been reviewed and the recent findings are discussed.

Vibration and Stability of Fluid Conveying Pipes with Stochastic Parameters

Flexible cantilever pipes conveying fluids with high velocity are analysed for their dynamic response and stability behaviour. The Young's modulus and mass per unit length of the pipe material have a stochastic distribution. The stochastic fields, that model the fluctuations of Young's modulus and mass density are characterized through their respective means, variances and autocorrelation functions or their equivalent power spectral density functions. The stochastic non self-adjoint partial differential equation is solved for the moments of characteristic values, by treating the point fluctuations to be stochastic perturbations. The second-order statistics of vibration frequencies and mode shapes are obtained. The critical flow velocity is-first evaluated using the averaged eigenvalue equation. Through the eigenvalue equation, the statistics of vibration frequencies are transformed to yield critical flow velocity statistics. Expressions for the bounds of eigenvalues are obtained, which in turn yield the corresponding bounds for critical flow velocities.

Stochastic Aeroelastic Analysis of Composite Helicopter Rotor

The effect of uncertainty in composite material properties on the aeroelastic response, vibratory loads, and stability of a hingeless helicopter rotor is investigated. The uncertainty impact on rotating natural frequencies of the blade is studied with Monte Carlo simulations and first-order reliability methods. The stochastic aeroelastic analyses in hover and forward flight are carried out with Monte Carlo simulations. The flap, lag, and torsion responses show considerable scatter from their baseline values, and the uncertainty impact varies with the azimuth angle. Furthermore, the blade response shows finite probability of resonance-type conditions caused by modal frequencies approaching multiples of the rotor speed. The 4/rev vibratory forces show large deviations from their baseline values. The lag mode damping shows considerable scatter due to uncertain material properties with an almost 40% probability of instability in hover.

Monte Carlo and molecular dynamics simulation of argon clusters andn-alkanes in the confined regions of zeolites

Geometry and energy of argon clusters confined in zeolite NaCaA are compared with those of free clusters. Results indicate the possible existence of magic numbers among the confined clusters. Spectra obtained from instantaneous normal mode analysis of free and confined clusters give a larger percentage of imaginary frequencies for the latter indicating that the confined cluster atoms populate the saddle points of the potential energy surface significantly. The variation of the percentage of imaginary frequencies with temperature during melting is akin to the variation of other properties. It is shown that confined clusters might exhibit inverse surface melting, unlike medium-to-large-sized free clusters that exhibit surface melting. Configurational-bias Monte Carte (CBMC) simulations of n-alkanes in zeolites Y and A are reported. CBMC method gives reliable estimates of the properties relating to the conformation of molecules. Changes in the conformational properties of n-butane and other longer n-alkanes such as n-hexane and n-heptane when they are confined in different zeolites are presented. The changes in the conformational properties of n-butane and n-hexane with temperature and concentration is discussed. In general, in zeolite Y as well as A, there is significant enhancement of the gauche population as compared to the pure unconfined fluid.

Multiscale stochastic approximation for parametric optimization of hidden Markov models

A two-time scale stochastic approximation algorithm is proposed for simulation-based parametric optimization of hidden Markov models, as an alternative to the traditional approaches to ''infinitesimal perturbation analysis.'' Its convergence is analyzed, and a queueing example is presented.

A molecular dynamics study and molecular level explanation of pressure dependence of ionic conductivity of potassium chloride in water

Experimental ionic conductivity of different alkali ions in water shows markedly different dependences on pressure. Existing theories such as that of Hubbard-Onsager are unable to explain these dependences on pressure of the ionic conductivity for all ions. We report molecular dynamics investigation of potassium chloride solution at low dilution in water at several pressures between 1 bar and 2 kbar. Two different potential models have been employed. One of the models successfully reproduces the experimentally observed trend in ionic conductivity of K+ ions in water over the 0.001-2 kbar range. We also propose a theoretical explanation, albeit at a qualitative level, to account for the dependence of ionic conductivity on pressure in terms of the previously studied Levitation Effect. It also provides a microscopic picture in terms of the pore network in liquid water.

Solvation dynamics in monohydroxy alcohols: Agreement between theory and different experiments

Recently three different experimental studies on ultrafast solvation dynamics in monohydroxy straight-chain alcohols (C-1-C-4) have been carried out, with an aim to quantify the time constant (and the amplitude) of the ultrafast component. The results reported are, however, rather different from different experiments. In order to understand the reason for these differences, we have carried out a detailed theoretical study to investigate the time dependent progress of solvation of both an ionic and a dipolar solute probe in these alcohols. For methanol, the agreement between the theoretical predictions and the experimental results [Bingemann and Ernsting J. Chem. Phys. 1995, 102, 2691 and Horng et al. J: Phys, Chern, 1995, 99, 17311] is excellent. For ethanol, propanol, and butanol, we find no ultrafast component of the time constant of 70 fs or so. For these three liquids, the theoretical results are in almost complete agreement with the experimental results of Horng et al. For ethanol and propanol, the theoretical prediction for ionic solvation is not significantly different from that of dipolar solvation. Thus, the theory suggests that the experiments of Bingemann and Ernsting and those of Horng et al. studied essentially the polar solvation dynamics. The theoretical studies also suggest that the experimental investigations of Joo et al. which report a much faster and larger ultrafast component in the same series of solvents (J. Chem. Phys. 1996, 104, 6089) might have been more sensitive to the nonpolar part of solvation dynamics than the polar part. In addition, a discussion on the validity of the present theoretical approach is presented. In this theory the ultrafast component arises from almost frictionless inertial motion of the individual solvent molecules in the force field of its neighbors.

Solvation dynamics of a charge bubble in water

A microscopic theory is used to calculate the solvation-time correlation function, (S(t)), of a light, non-stationary charge bubble in water. The calculated correlation function is found to be similar to the energy-time correlation function of a solvated electron. The ionic mobility of a charge bubble of the size of the hydrated electron is also calculated. It is found that the mobility of the charge plays a very important role in its own solvation.

Strong shock approximation in the new theory of shock dynamics

This paper investigates the propagation of a strong shock into an inhomogeneous medium using the new theory of shock dynamics. The equations are simple to solve and involve no trial-and-error method commonly used in this case. The results compare favourably with earlier results obtained in the case of self-similar flows, which arise as a special case of this theory.