Isospin- and momentum-dependent interaction in isospin-dependent quantum molecular dynamics

An isospin-dependent quantum molecular dynamical model (IQMD) is developed, with the isospin degree of freedom in the momentum-dependent interaction(MDI) included in IQMD, to obtain an isospin- and momentum-dependent interaction (IMDI) in IQMD. We investigate the effect of IMDI on the isospin fractionation ratio and its dynamical mechanism in the intermediate energy heavy ion collisions. It is found that the IMDI induces the significant reductions in the isospin fractionation ratio for all of beam energies, impact parameters, neutron-proton ratios and mass number of colliding systems. However, the strong dependence of isospin fractionation ratio on the symmetrical potential is preserved, with the isospin degree of freedom included in the MDI, i.e. the isospin fractionation ratio is still a good probe for extracting the information about the equation of state of isospin asymmetrical nuclear matter.

Molecular dynamics and control following excitation with an ultra-short intense laser pulse

Vibronic excitations of the tri-atomic molecule OClO (A(2)A(2)(nu(1), nu(2), nu(3)) <-- (XB1)-B-2 (0, 0, 0)) with weak and strong ultra-short laser fields are studied within full quantum wavepacket dynamics in hyperspherical coordinates. Different dynamics is observed following excitation with laser pulses of different intensities. With a strong laser pulse, many vibrational states are excited and a spatially more localised wavepacket arises. The numerical results show that the population of different vibrational states of the wavepacket on the excited potential energy surface is altered by the intensity of the laser pulse. The numerical results also suggest a related effect on the phase of the wavepacket. These interesting phenomena can be understood by an analysis of the corresponding results for two model diatomic molecules. The possible physical mechanisms of control of chemical processes using strong laser fields are discussed. (C) 2004 Elsevier B.V. All rights reserved.

Study on structure and molecular dynamics of starch/poly(sodium acrylate)grafted superabsorbent by C-13 solid state NMR

The domain-structure of samples containing a series of starch/poly(sodium acrylate)-grafted superabsorbents, pure starch, pure poly(sodium acrylate), and blend of starch/poly(sodium acrylate) has been studied by high-resolution solid-state C-13 NMR spectroscopy at room temperature. The result shows that the crystallinity of starch decreases greatly in the grafted and blended samples.

Molecular dynamics study on the phase diagrams of linear and branched chain molecules

The particle transfer molecular dynamics is used to study the phase equilibria of linear and branched chain molecules. The scaling of the critical temperature versus chain length is obtained and the critical densities are found to decrease with increasing chain length, which are in agreement with the results of experiment and theory. The phase diagrams of the linear and the branched chain molecules nearly overlap with each other. Moreover, the radial distribution functions of linear and branched chain molecules in gas phase are very similar, but in the liquid phase, they are different for different kinds of chains.

Calculating the equation of state parameters and predicting the spinodal curve of isotactic polypropylene/poly(ethylene-co-octene) blend by molecular dynamics simulations combined with Sanchez-Lacombe lattice fluid theory

Molecular dynamics simulations are adopted to calculate the equation of state characteristic parameters P*, rho*, and T* of isotactic polypropylene (iPP) and poly(ethylene-co-octene) (PEOC), which can be further used in the Sanchez-Lacombe lattice fluid theory (SLLFT) to describe the respective physical properties. The calculated T* is a function of the temperature, which was also found in the literature. To solve this problem, we propose a Boltzmann fitting of the data and obtain T* at the high-temperature limit. With these characteristic parameters, the pressure-volume-temperature (PVT) data of iPP and PEOC are predicted by the SLLFT equation of state. To justify the correctness of our results, we also obtain the PVT data for iPP and PEOC by experiments. Good agreement is found between the two sets of data. By integrating the Euler-Lagrange equation and the Cahn-Hilliard relation, we predict the density profiles and the surface tensions for iPP and PEOC, respectively. Furthermore, a recursive method is proposed to obtain the characteristic interaction energy parameter between iPP and PEOC. This method, which does not require fitting to the experimental phase equilibrium data, suggests an alternative way to predict the phase diagrams that are not easily obtained in experiments.

Temperature dependence of the distribution of the first passage time: Results from discontinuous molecular dynamics simulations of an all-atom model of the second beta-hairpin fragment of protein G

More than 22 000 folding kinetic simulations were performed to study the temperature dependence of the distribution of first passage time (FPT) for the folding of an all-atom Go-like model of the second beta-hairpin fragment of protein G. We find that the mean FPT (MFPT) for folding has a U (or V)-shaped dependence on the temperature with a minimum at a characteristic optimal folding temperature T-opt*. The optimal folding temperature T-opt* is located between the thermodynamic folding transition temperature and the solidification temperature based on the Lindemann criterion for the solid. Both the T-opt* and the MFPT decrease when the energy bias gap against nonnative contacts increases. The high-order moments are nearly constant when the temperature is higher than T-opt* and start to diverge when the temperature is lower than T-opt*. The distribution of FPT is close to a log-normal-like distribution at T* greater than or equal to T-opt*. At even lower temperatures, the distribution starts to develop long power-law-like tails, indicating the non-self-averaging intermittent behavior of the folding dynamics. It is demonstrated that the distribution of FPT can also be calculated reliably from the derivative of the fraction not folded (or fraction folded), a measurable quantity by routine ensemble-averaged experimental techniques at dilute protein concentrations.

Molecular dynamics simulation of sub-transition for polyethersulfone

Molecular dynamics (MD) simulations of a polyethersulfone (PES) chain are carried out in the amorphous state by using the Dreiding 2.21 force field at four temperatures. Two types of molecular motion, i.e, rotations of phenylene rings and torsions of large segments containing two oxygen atoms, two sulfur atoms, and five phenylene rings on the backbone, are simulated. The modeling results show that the successive phenylene rings should be in-phase cooperative rotations, whereas the successive large segments should be out-of-phase cooperative torsions. By calculating the diffusion coefficient for the phenylene ring rotations, it is found that this rotation contributes to the beta -transition of PES.

Phase transition temperatures of side chain liquid crystal polymer from molecular dynamics simulation

In this paper, phase transition temperatures of side chain liquid crystal polymer were predicted by molecular dynamics simulation. We analyzed the change of energy and the degree of similarity(S) with the temperature varying. The simulated phase transition temperatures agree with the experimental values in a proper deviation.

Molecular dynamics simulation of polystyrene-block-poly(methyl methacrylate)

Molecular dynamics is applied to the system of polystyrene-block-poly(methyl methacrylate). The simulation shows that for the block copolymer system, a layered structure, which reflects microphase separation, is obtained and this structure is stable. In order to elucidate that the formation of the layered structure is reasonable, some static properties such as the radial distribution function and the dipole moment are analyzed in some detail.

Docking and molecular dynamics study on the inhibitory activity of coumarins on aldose reductase

In order to explore the inhibitory mechanism of coumarins toward aldose reductase (ALR2), AutoDock and Gromacs software were used for docking and molecular dynamics studies on 14 coumarins (CM) and ALR2 protease. The docking results indicate that residues TYR48, HIS110, and TRP111 construct the active pocket of ALR2 and, besides van der Waals and hydrophobic interaction, CM mainly interact with ALR2 by forming hydrogen bonds to cause inhibitory behavior. Except for CM1, all the other coumarins take the lactone part as acceptor to build up the hydrogen bond network with active-pocket residues. Unlike CM3, which has two comparable binding modes with ALR2, most coumarins only have one dominant orientation in their binding sites. The molecular dynamics calculation, based on the docking results, implies that the orientations of CM in the active pocket show different stabilities. Orientation of CM1 and CM3a take an unstable binding mode with ALR2; their conformations and RMSDs relative to ALR2 change a lot with the dynamic process. While the remaining CM are always hydrogen-bonded with residues TYR48 and HIS110 through the carbonyl O atom of the lactone group during the whole process, they retain the original binding mode and gradually reach dynamic equilibrium.

Docking and molecular dynamics studies toward the binding of new natural phenolic marine inhibitors and aldose reductase

Phenolic marine natural product is a kind of new potential aldose reductase inhibitors (ARIs). In order to investigate the binding mode and inhibition mechanism, molecular docking and dynamics studies were performed to explore the interactions of six phenolic inhibitors with human aldose reductase (hALR2). Considering physiological environment, all the neutral and other two ionized states of each phenolic inhibitor were adopted in the simulation. The calculations indicate that all the inhibitors are able to form stable hydrogen bonds with the hALR2 active pocket which is mainly constructed by residues TYR48, HIS110 and TRP111, and they impose the inhibition effect by occupying the active space. In all inhibitors, only La and its two ionized derivatives La_ion1 and La_ion2, in which neither of the ortho-hydrogens of 3-hydroxyl is substituted by Br, bind with hALR2 active residues using the terminal 3-hydroxyl. While, all the other inhibitors, at least one of whose ortho-sites of 3- and 6-hydroxyls are substituted by Br substituent which take much electron-withdrawing effect and steric hindrance, bind with hALR2 through the lactone group. This means that the Br substituent can effectively regulate the binding modes of phenolic inhibitors. Although the lactone bound inhibitors have relatively high RMSD values, our dynamics study shows that both binding modes are of high stability. For each inhibitor molecule, the ionization does not change its original binding mode, but it does gradually increase the binding free energy, which reveals that besides hydrogen bonds, the electrostatic effect is also important to the inhibitor–hALR2 interaction.