Orientational correlations in molecular liquid SnI4

The total scattering structure factor of liquid tin tetraiodide (SnI(4)) has been interpreted by means of reverse Monte Carlo (RMC) modeling. From the sets of particle coordinates provided by RMC, which are consistent with experimental results within errors, partial radial distribution functions as well as correlation functions characterizing mutual orientations of molecules as a function of distance between molecular centers can be calculated. Interestingly and very much in contrast to liquids of symmetric XCl(4) molecules, the corner-to-face (or "Apollo")-type orientation of neighboring molecules has a significant (about 20%) occurrence in liquid SnI(4). Via comparison with a reference system, obtained by hard sphere Monte Carlo simulation, we demonstrate that intermolecular two-body correlations in liquid SnI(4) are determined largely by excluded volume (steric) effects; that is, intermolecular two-body interactions play only a minor role. On the other hand, as it is manifested in the large difference between the reference and "real" systems in terms of the orientational correlations, higher order interactions are indispensable. This feature can explain the extremely rich phase behavior of SnI(4) at high pressures.     

Orientational and translational correlations of liquid methane over the nanometer-picosecond scales by molecular dynamics simulation and inelastic neutron scattering

Five models for the site-site intermolecular pair interactions of methane are compared in some detail and used to investigate both structural and dynamical properties of the dense liquid deuteromethane by means of molecular dynamics (MD) simulations. The orientational distribution probabilities of molecular pairs are carefully analyzed for each anisotropic potential model. We propose a revision of existing classification methods used to group the innumerable relative orientations of methane-methane pairs into six basic geometries. With this new approach, our results for the probability of the six basic categories as a function of the intermolecular distance are different from the ones present in the literature, where the role of the angular spread on the anisotropic interaction energy is not taken in full consideration and certain configurations with no significant change in the pair-potential are assigned to different categories. The analysis of the static orientational correlations in liquid methane and the prevalence of certain configurations in different ranges guide the subsequent discussion of the MD model-dependent results for the dynamic structure factor. Comparison with our inelastic neutron scattering results for liquid CD(4) at the nanometer and picosecond space and time scales allows us to confirm the full adequacy of the Tsuzuki, Uchimaru and Tanabe model of 1998 with respect to more recent potentials.     

Hydration of superoxide studied by molecular dynamics simulation

Molecular dynamics was used to study the hydration of superoxide (O). The Helmholtz free energy of hydration of O was estimated by the thermodynamic integration method. The diffusion of O and the water structure around O were also studied. Two water models were used in the calculations and the results were compared to experiments.     

Dimer asymmetry in superoxide dismutase studied by molecular dynamics simulation

Summary Molecular dynamics (MD) simulations of 100 ps have been carried out to study the active-site behaviour of the Cu,Zn superoxide dismutase dimer (SOD) in water. The active site of each subunit was monitored during the whole simulation by calculating the distances between functional residues and the catalytic copper. The results indicate that charge orientation is maintained at each active site but the solvent accessibility varies. Analysis of the MD simulation, carried out by using the atomic displacement covariance matrix, has shown a different intra-subunit correlation pattern for the two monomers and the presence of inter-subunit correlations. The MD simulation presented here indicates an asymmetry in the two active sites and different dynamic behaviour of the two SOD subunits.     

Backbone dynamics of proteins studied by two-dimensional heteronuclear NMR spectroscopy and molecular dynamics simulations

To investigate the backbone dynamics of proteins ¹⁵N longitudinal and transverse relaxation experiments combined with {¹H, ¹⁵N{ NOE measurements together with molecular dynamics simulations were carried out using ribonuclease T₁ and the complex of ribonuclease T₁ with 2′GMP as a model protein. The intensity decay of individual amide cross peaks in a series of (¹H, ¹⁵N)HSQC spectra with appropriate relaxation periods was fitted to a single exponential by using a simplex algorithm in order to obtain ¹⁵N T₁ and T₂ relaxation times. The relaxation times were analyzed in terms of the “model-free” approach introduced by Lipari and Szabo. In addition, a nanosecond molecular dynamics (MD ) simulation of ribonuclease T₁ and its 2′GMP complex in water was carried out. The angular reorientations of the backbone amide groups were classified with several coordinate frames following a transformation of NH vector trajectories. In this study, NH librations and backbone dihedral angle fluctuations were distinguished. The NH bond librations were found to be similar for all amides as characterized by correlation times of librational motions in a subpicosecond scale. The angular amplitudes of these motions were found to be about 10°–12° for out-of-plane displacements and 3°–5° for the in-plane displacement. The contributions from the much slower backbone dihedral angle fluctuations strongly depend on the secondary structure. The dependence of the amplitude of local motion on the residue location in the backbone is in good agreement with the results of NMR relaxation measurements and the X-ray data. The protein dynamics is characterized by a highly restricted local motion of those parts of the backbone with defined secondary structure as well as by a high flexibility in loop regions. Comparison of the MD and NMR data of the free liganded enzyme ribonuclease T₁ clearly indicates a restriction of the mobility within certain regions of the backbone upon inhibitor binding. © 1996 John Wiley & Sons, Inc.     

Stability of two-dimensional crystalline aggregates of a protein studied by molecular dynamics

Two-dimensional protein (ferritin) aggregates with a square lattice symmetry, which were formed within a thin liquid layer on a mercury surface, were studied by molecular dynamics (MD) simulation. For the simulation, the ferritin molecule was modeled by an assembly of 49 spheres, and the intermolecular interactions were given by simple formulae. During the simulation, molecules were confined within a layer, which corresponds to the thin liquid layer. An annealing MD simulation was done starting from a random molecular configuration within the layer, and aggregates with the square lattice symmetry were also obtained. To study the stability of aggregates, dissociation processes of the aggregates were analyzed using MD simulations at room temperature. Interactions between the nearest-neighbor molecules were regarded as bonds. Mean bond energies and correlation coefficients between the bond energies were calculated from the MD trajectories. A decay profile according to the dissociation was obtained, yielding a dissociation rate constant. Buried bonds were stronger than peripheral bonds. The larger the aggregate size, the stronger the bond for each of the buried and peripheral bonds. A simple theoretical account, which is applicable to a general bonded network, was introduced to analyze the dynamics of the aggregates. © 1994 by John Wiley & Sons, Inc.     

Intermolecular polarizability dynamics of aqueous formamide liquid mixtures studied by molecular dynamics simulations

A molecular dynamics simulation study is presented for the relaxation of the polarizability anisotropy in liquid mixtures of formamide and water, using a dipolar induction scheme that involves the intrinsic polarizability and first hyperpolarizability tensors of the molecules, and the dipole-quadrupole polarizability of water species. The long time diffusive decay of the collective polarizability anisotropy correlations exhibits a substantial slowing down as the formamide mole fraction increases in the mixture. The diffusive times for the polarizability relaxation obtained from the authors' simulations are in good agreement with optical Kerr effect experimental data, and they are found to correlate nearly linearly with the estimated mean lifetimes of the hydrogen bonds within the mixture, suggesting that the relaxation of the hydrogen bond network is responsible to some extent for the collective relaxation of the polarizability anisotropy of the mixture. The short time behavior of the polarizability anisotropy relaxation was investigated by computing the nuclear response function, R(t), which is very rapidly dominated by the formamide contribution as it is added to water, due to the much larger polarizability anisotropy of formamide molecules compared to that of water. Several contributions to the Raman spectrum were also analyzed as a function of composition, and the dynamical origin of the different bands was determined.     

Clustering of glycine molecules in aqueous solution studied by molecular dynamics simulation

The nature of glycine--glycine interactions in aqueous solution has been studied using molecular dynamics simulations at four different concentrations and, in each case, four different temperatures. Although evidence is found for formation of small, transient hydrogen-bonded clusters of glycine molecules, the main type of interaction between glycine molecules is found to be single NH...OC hydrogen bonds. Double-hydrogen-bonded "dimers", which have often been cited as a significant species present in aqueous solutions of glycine, are only observed infrequently. When double-hydrogen-bonded dimers are formed, they dissociate quickly (typically within less than ca. 4 ps), although the broken hydrogen bonds have a higher than average probability of reforming. Several aspects of the clustering of glycine molecules are investigated as a function of both temperature and concentration, including the size distribution of glycine clusters, the radii of gyration of the clusters, and aspects of the lifetimes of glycine-glycine hydrogen bonding by means of hydrogen-bond correlation functions. Diffusion coefficients for the glycine clusters and water molecules are also investigated and provide results in realistic agreement with experimental results.     

Orientational ordering of ferroelectric thiobenzoate liquid crystals

The orientational ordering of a series of ferroelectric thiobenzoate liquid crystals was studied by natural abundance ¹³C NMR spectroscopy. The technique used was a combination of variable angle spinning (VAS) and separated local field spectroscopy (SLF). With rapid sample spinning about an axis forming an angle of c. 45° with respect to the magnetic field, the smectic director aligns parallel to the spinning axis, leading to narrow peaks in the ¹³C NMR spectrum. The two-dimensional NMR spectroscopic method SLF allows the observation of decoupled ¹³C signals in the ω2 dimension and first-order Csbnd;H splitting patterns in the ω1 dimension, from which the dipolar C-H coupling constants for individual bonds can be obtained. The order parameters for different molecular segments of eight different compounds, all containing two phenyl rings linked by a thioester group, were obtained this way. A considerable influence of length, branching and chirality of the aliphatic chain on the order parameters was observed.     

Orientational dynamics of water trapped between two nanoscopic hydrophobic solutes: A molecular dynamics simulation study

We investigate thoroughly the effect of confinement and solute topology on the orientational dynamics of water molecule in the interplate region between two nanoscopic hydrophobic paraffinlike plates. Results are obtained from molecular dynamics simulations of aqueous solutions of paraffinlike plates in the isothermal-isobaric ensemble. An analysis of survival time auto correlation function shows that the residence time of the water molecule in the confined region between two model nanoscopic hydrophobic plates depends on solute surface topology (intermolecular distance within the paraffinlike plate). As expected, the extent of confinement also changes the residence time of water molecules considerably. Orientational dynamics was analyzed along three different directions, viz., dipole moment, HH, and perpendicular to molecular plane vectors. It has been demonstrated that the rotational dynamics of the confined water does not follow the Debye rotational diffusion model, and surface topology of the solute plate and the extent of confinement have considerable effect on the rotational dynamics of the confined water molecules.     

DNA strand break: structural and electrostatic properties studied by molecular dynamics simulation

Due to their lethal consequences and a relatively high probability of introduction of repair errors and mutations, single and double strand breaks are among the most important and dangerous DNA lesions. However, the mechanisms of their recognition and repair processes are only poorly known at present. This work defines and analyzes a DNA with single strand break as a template study for future complex analyses of biologically serious double strand break damage and its enzymatic repair mechanisms. Besides a non-damaged DNA serving as a reference system with no surprising results, system with open valences of the atoms at the strand break ends as well as a system with filled valences were simulated. In both cases during the first few nanoseconds the broken ends of strand breaks are significantly exposed to the outside of the molecule. However, with increasing time, the system with single strand break with open valences is partially disrupted. On the contrary, the system with filled valences shows stable conformation with newly created hydrogen bond between the two strand break endings. Moreover, these endings are steadily situated in the inner part of the molecule, thus making the recognition and docking process of a repair enzyme more complicated in the case of filled valences.     

First principles molecular dynamics simulation of liquid rubidium

We present a first principles molecular dynamics simulation of liquid rubidium. The atomic forces are obtained from a quantum mechanical calculation of its electronic structure within the local density approximation of the density functional formalism and using the pseudopotential plane-wave method. We compare our results with the structure and dynamics predicted by classical molecular dynamics simulations.     

Car-parrinello molecular dynamics simulation of liquid formic acid

First-principles molecular dynamics has been used to investigate the structural, vibrational, and energetic properties of formic acid, formic acid-formate anion dimers, and liquid formic acid in a periodically repeated box with 32 formic acid molecules. We found that in liquid formic acid the hydrogen-bonded clusters mainly consist of linear branching chains. From our simulation, we got good agreement with the available structural and dynamical data. We also studied the proton transfer in the cis-formic acid-formate anion dimer, and we showed that this proton transfer does not have any potential barrier. The hydrogen bonding statistics as well as the mean lifetime of the hydrogen bonds are analyzed.     

Computer simulations of molecular motion in liquid crystals by the method of Brownian dynamics

Computer simulations of the molecular motion of polymer chains in the presence of a strong nematic field were carried out by the method of Brownian dynamics. Two models were studied: the first model (linear liquid crystal) is a freely jointed chain with rigid bonds, the second model (comb-like liquid crystal) is a chain with fixed bond angles and rigid side groups. The influence of ordering on chain conformations, orientational and translational mobility and spectra of relaxation times was investigated.     

Orientational optical nonlinearity of nematic liquid crystals induced by high-molecular-mass azo-containing compounds

Processes of light-induced reorientation of nematic liquid-crystalline molecules induced by the addition of low concentrations (0.1–2.0 wt %) of comb-shaped polymers and carbosilane dendrimers containing azobenzene fragments are studied. When the molecular structure of the above compounds becomes more complicated, the induced orientational nonlinearity increases. The introduction of 2G and 3G dendrimers into a nematic has for the first time made it possible to visualize and study a purely optical first-order Freedericksz transition in the field of a linearly polarized wave.     

Orientational dynamics and energy landscape features of thermotropic liquid crystals: An analogy with supercooled liquids

Recent optical kerr effect (OKE) studies have revealed that orientational relaxation of rodlike nematogens near the isotropic-nematic (I-N) phase boundary and also in the nematic phase exhibit temporal power law decay at intermediate times. Such behaviour has drawn an intriguing analogy with supercooled liquids. Here, we have investigated the single-particle and collective orientational dynamics of a family of model system of thermotropic liquid crystals using extensive computer simulations. Several remarkable features of glassy dynamics are on display including non-exponential relaxation, dynamical heterogeneity, and non-Arrhenius temperature dependence of the orientational relaxation time. Over a temperature range near the I-N phase boundary, the system behaves like a fragile glass-forming liquid. Using proper scaling, we construct the usual relaxation time versus inverse temperature plot and explicitly demonstrate that one can successfully define a density dependent fragility of liquid crystals. The fragility of liquid crystals shows a temperature and density dependence which is remarkably similar to the fragility of glass forming supercooled liquids. Energy landscape analysis of inherent structures shows that the breakdown of the Arrhenius temperature dependence of relaxation rate occurs at a temperature that marks the onset of the growth of the depth of the potential energy minima explored by the system.     

The plastic and liquid phases of CCl3Br studied by molecular dynamics simulations

We present a molecular dynamics study of the liquid and plastic crystalline phases of CCl(3)Br. We investigated the short-range orientational order using a recently developed classification method and we found that both phases behave in a very similar way. The only differences occur at very short molecular separations, which are shown to be very rare. The rotational dynamics was explored using time correlation functions of the molecular bonds. We found that the relaxation dynamics corresponds to an isotropic diffusive mode for the liquid phase but departs from this behavior as the temperature is decreased and the system transitions into the plastic phase.     

Twist viscosities and flow alignment of biaxial nematic liquid crystal phases of a soft ellipsoid-string fluid studied by molecular dynamics simulation

We have calculated the twist viscosity and the alignment angle between the director and the stream lines in shear flow of a liquid crystal model system, which forms biaxial nematic liquid crystals, as functions of the density, from the Green-Kubo relations by equilibrium molecular dynamics simulation and by a nonequilibrium molecular dynamics algorithm, where a torque conjugate to the director angular velocity is applied to rotate the director. The model system consists of a soft ellipsoid-string fluid where the ellipsoids interact according a repulsive version of the Gay-Berne potential. Four different length-to-width-to-breadth ratios have been studied. On compression, this system forms discotic or calamitic uniaxial nematic phases depending on the dimensions of the molecules, and on further compression a biaxial nematic phase is formed. In the uniaxial nematic phase there is one twist viscosity and one alignment angle. In the biaxial nematic phase there are three twist viscosities and three alignment angles corresponding to the rotation around the various directors and the different alignments of the directors relative to the stream lines, respectively. It is found that the smallest twist viscosity arises by rotation around the director formed by the long axes, the second smallest one arises by rotation around the director formed by the normals of the broadsides, and the largest one by rotation around the remaining director. The first twist viscosity is rather independent of the density whereas the last two ones increase strongly with density. One finds that there is one stable director alignment relative to the streamlines, namely where the director formed by the long axes is almost parallel to the stream lines and where the director formed by the normals of the broadsides is almost parallel to the shear plane. The relative magnitudes of the components of the twist viscosities span a fairly wide interval so this model should be useful for parameterisation experimental data.     

A molecular dynamics simulation study of nanoscale liquid threads

The properties of liquid threads in angular direction are studied. On the basis of Gibbs theory, a thermodynamic model is constructed and a formula of surface of tension is derived. For the determination of surface tension, different methods are presented. We investigate seven different systems (the numbers of molecules N are 1600, 2240, 2880, 3360, 4000, 4800, and 5280) by molecular dynamics simulations. We analyze the surface tension obtained by various routes, then whether the classical Laplace equation is applicable in nanoscale can be determined. The results obtained by molecular dynamics simulations support that surface tension is dependent on the dividing surface curvature.