Accelerated Molecular Dynamics Study of Z/E Isomerization of Azobenzene: Kramers’ Theory Validation

The reaction rate in a condensed phase is retarded by solvent friction and deviates from the transition state theory (TST) in a highly viscous solvent. The Z/E isomerization kinetics of 4-dimethylamino-4′-nitroazobenzene in solution have been experimentally analyzed in a wide range of pressure in our previous studies. Herein the non-TST behaviors of the rate constant were theoretically modelled by means of Kramers’ theory and its extensions and by molecular dynamics (MD) simulations. The Z/E isomerization takes place in milliseconds or longer and so is not reachable by normal MD simulations. The accelerated MD method was employed to simulate the rare events using a biased isomerization potential barrier.     

Fluorescence spectra of organic dyes in solution: a time dependent multilevel approach

Classical all-atom molecular dynamics (MD) simulations and quantum mechanical (QM) time-dependent density functional theory (TD-DFT) calculations are employed to study the conformational and photophysical properties of the first emitter excited state of tetramethyl-rhodamine iso-thiocyanate fluorophore in aqueous solution. For this purpose, a specific and accurate force field has been parameterised from QM data to model the fluorophore's first bright excited state. During the MD simulations, the consequences of the π→π* electronic transition on the structure and microsolvation sphere of the dye has been analysed in some detail and compared to the ground state behaviour. Thereafter, fluorescence has been calculated at the TD-DFT level on configurations sampled from the simulated MD trajectories, allowing us to include time dependent solvent effects in the computed emission spectrum. The latter, when compared with the absorption spectrum, reproduces well the experimental Stokes shift, further validating the proposed multilevel computational procedure.     

Calculating the hopping times of confined fluids: two hard disks in a box

The dynamical transition between the anomalous single file diffusion of highly confined fluids and bulk normal diffusion can be described by a phenomenological model involving a particle hopping time tau(hop). We suggest a theoretical formalism that will be useful for the calculation of tau(hop) for a variety of systems and test it using a simple model consisting of two hard disks confined to a rectangular box with hard walls. In the case where the particles are moving diffusively, we find the hopping time diverges as a power law in the threshold region with an exponent of -(3/2). Under conditions where the particles move inertially, transition state theory predicts a power law behavior with an exponent of -2. Molecular dynamics simulations confirm the transition state theory result for inertial dynamics, while Brownian dynamics simulations suggest the scaling exponent is highly sensitive to the details of the algorithm.     

Selective sampling of transition paths

In this short paper, we introduce an approximate method for the quick estimate of rate constants based on a simple sampling method of reactive transition paths over high energy barriers. It makes use of the previously introduced accelerated molecular dynamics (MD) simulation method to generate initial points for trajectory shooting. The accelerated MD simulations, although with the loss of real dynamics, lead to a quick calculation of thermodynamic properties and at the same time produce an ensemble of configurations with an enhanced sampling over the phase space that is more "reactive." The forward/backward trajectory shooting as that used in the transition path sampling method is then initiated from the configurations obtained from accelerated MD simulations to generate transition paths on the original unbiased potential. This method selectively enhances sampling of successful trajectories and at the same time accelerates significantly the calculation of rate constants.     

Anomalous mixing behavior of polyisobutylene/polypropylene blends: molecular dynamics simulation study

The unusual mixing behavior of polyisobutylene (PIB) with head-to-head (hhPP) and head-to-tail polypropylene (PP) is studied using large-scale molecular dynamics (MD). The heats of mixing and Flory chi parameters were computed from MD simulations of both blends using a united atom model. The chi parameters from the simulations were estimated from the structure factors using the random phase approximation in analogy with neutron scattering (SANS) experiments. MD simulations for syndiotactic hhPP/PIB predicted a lower critical solution temperature with a chi parameter in very good agreement with SANS experiments on the atactic hhPP/PIB blend. MD simulations also predicted that the isotactic PP/PIB blend was immiscible at high molecular weight in qualitative agreement with cloud point measurements on atactic PP/PIB.     

Effect of chain conformational change on micelle structures: experimental studies and molecular dynamics simulations

The effect of chain conformation change on the self-assembly behavior of poly(gamma-benzyl- l-glutamate)-block-poly(ethylene glycol) (PBLG-b-PEG) was studied both experimentally by transmission electron microscopy, laser light scattering, and circular dichroism and computationally using molecular dynamics (MD) simulation. It was found that, by introducing trifluoroacetic acid to the PBLG-b-PEG solution, the conformation of the PBLG chain transforms from alpha-helix to random coil, which results in a change of the micelle structures formed by PBLG-b-PEG from rod to sphere. Meanwhile, the MD simulations were performed by using Brownian dynamics on the self-assembly behavior of model AB-type diblock copolymers with various chain rigidities of the A-block. The results show that, by decreasing the fraction of rigid chain conformation of the A-block, which corresponds to the helix-coil transition in the PBLG-b-PEG sample, the aggregate structure transforms from rod to sphere. The MD simulations also provide chain packing information in the micelles. On the basis of both experimental and MD simulation results, the mechanism regarding the effect of the conformation change of the polypeptide block copolymer on its self-association behavior is suggested.     

氟化钾团簇的相变、成核和重结晶的分子动力学模拟

我们利用Born-Mayer-Huggins相互作用势函数对(KF)_N(N=108，256，500和864)团簇进行了分子动力学(MD)模拟。为了避免周期性边界条件对相变、成核和重结晶的干扰作用，对体系采用了自由边界。基于MD模拟结果，对团簇的熔化温度、熔化焓、扩散系数、成核速率、固液界面自由能、临界核大小等进行了计算和讨论。在对(KF)₈₆₄双晶团簇的热退火模拟中，观察到了固态的重结晶和晶粒的生长。经典的成核理论成功地解释了(KF)₈₆₄双晶团簇的重结晶MD模拟结果。     

Ab initio direct dynamics studies on the reaction of H atom with CH3CH2Cl

The multiple channel reaction H + CH(3)CH(2)Cl --> products has been studied by the ab initio direct dynamics method. The potential energy surface information is calculated at the MP2/6-311G(d,p) level of theory. The energies along the minimum energy path are further improved by single-point energy calculations at the PMP4(SDTQ)/6-311+G(3df,2p) level of theory. For the reaction, four reaction channels (one chlorine abstraction, one alpha-hydrogen abstraction, and two beta-hydrogen abstractions) have been identified. The rate constants for each reaction channel are calculated by using canonical variational transition state theory incorporating the small-curvature tunneling correction in the temperature range 298-5000 K. The total rate constants, which are calculated from the sum of the individual rate constants, are in good agreement with the experimental data. The calculated temperature dependence of the branching fractions indicates that for the title reaction, H-abstraction reaction is the major reaction channel in the whole temperature range 298-5000 K.     

Theoretical Investigation for the Abstraction Reaction of H with (CH_3)_3GeH

The reaction of H atom with (CH3)3GeH is considered to play important role in chemical vapor deposition (CVD) processes used in the semiconductor industry 1-2. The reaction mechanism and kinetics nature for this reaction are therefore essential input data for computer-modelling studies directed towards obtaining an understanding of the factors controlling CVD processes. However, despite its importance, the kinetics work about this reaction was very limited. Only two groups studied exper…     

Molecular dynamics simulation of the A-DNA to B-DNA transition in aqueous RbCl solution

Unrestrained molecular dynamics (MD) simulations have been carried out to characterize the stability of DNA conformations and the dynamics of A-DNA→B-DNA conformational transitions in aqueous RbCl solutions. The PARM99 force field in the AMBER8 package was used to investigate the effect of RbCl concentration on the dynamics of the A→B conformational transition in the DNA duplex d(CGCGAATTCGCG)2 . Canonical Aand B-form DNA were assumed for the initial conformation and the final conformation had a length per complete turn that matched the canonical B-DNA. The DNA structure was monitored for 3.0 ns and the distances between the C5′ atoms were obtained from the simulations. It was found that all of the double stranded DNA strands of A-DNA converged to the structure of B-form DNA within 1.0 ns during the unrestrained MD simulations. In addition, increasing the RbCl concentration in aqueous solution hindered the A→B conformational transition and the transition in aqueous RbCl solution was faster than that in aqueous NaCl solution for the same electrolyte strength. The effects of the types and concentrations of counterions on the dynamics of the A→B conformational transition can be understood in terms of the variation in water activity and the number of accumulated counterions in the major grooves of A-DNA. The rubidium ion distributions around both fixed A-DNA and B-DNA were obtained using the restrained MD simulations to help explain the effect of RbCl concentration on the dynamics of the A→B conformational transition.     

Ab initio Calculation and Kinetic Study for the Abstraction Reaction of H with Sihcl3

ＴｒｉｃｈｌｏｒｏｓｉｌａｎｅｉｓａｎｉｍｐｏｒｔａｎｔｍａｔｅｒｉａｌｉｎｐｌａｓｍａＣｈｅｍｉｃａｌＶａｐｏｒＤｅｐｏｓｉｔｉｏｎ (ＣＶＤ)ａｎｄｉｎｓｅｍｉｃｏｎｄｕｃｔｏｒｄｅｖｉｃｅｐｒｏｃｅｓｓ .1 4 Ｔｈｅｒｅａｃｔｉｏｎｏｆｔｒｉｃｈｌｏｒｏｓｉｌａｎｅｗｉｔｈａｔｏｍｉｃｈｙｄｒｏｇｅｎ ,ｔｈｅｓｉｍｐｌｅｓｔｆｒｅｅ ｒａｄｉｃａｌｓｐｅｃｉｅｓ,ｈａｓｄｒａｗｎｃｏｎｓｉｄｅｒａｂｌｅａｔｔｅｎｔｉｏｎ :ｋｉｎｅｔｉｃｐａｒａｍｅｔｅｒｓｆｏｒＨ ａｔｏｍｒｅａｃｔｉｏｎａｒｅｄ…     

H原子和SiH2Cl2抽提反应的理论研究

对H+SiH2Cl2反应进行了详细的理论研究,理论证明了抽提氢的通道是唯一可行的反应通道。并在从头算给出的电子结构信息基础上,用变分过渡态理论(CVT)加小曲率隧道效应校正(SCT)等方法对该反应进行了直接的动力学研究,得到该反应的理论速率常数,并详细讨论了各动力学参数沿反应坐标的变化。在较宽的温度范围内,反应速率常数表现出非Arrhenius行为,用三参数公式似合了速-温关系式,为k(T)=(1.32×10^-22)T^3.67exp(-26/T)。理论计算的速率常数与实验数值符合得很好。     

Mechanistic study and kinetic properties of the CF3CHO + Cl reaction

Theoretical investigations have been carried out on the mechanism and kinetics for the reaction of CF 3 CHO + Cl using duallevel direct dynamics method. The potential energy surface information was obtained at the MCQCISD/3//MP2/cc-pVDZ level and the kinetic calculations were done using variational transition state theory with interpolated single-point energy (VTST-ISPE) approach. The calculated results show that the reaction proceeds primarily via the H-abstraction channel, while the Cl-addition channel is unfavorable due to the higher barriers. The improved canonical variational transition-state theory (ICVT) with the small-curvature tunneling correction (SCT) was used to calculate the rate constants. The theoretical rate constants at room temperature are in general agreement with the experimental values. A three-parameter rate constant expression was fitted over a wide temperature range of 200-2000 K.     

Application of Fokker-Planck-Kramers equation treatment for short-time dynamics of diffusion-controlled reaction in supercritical Lennard-Jones fluids over a wide density range

The validity of a Fokker-Planck-Kramers equation (FPKE) treatment of the rate of diffusion-controlled reaction at short times [K. Ibuki and M. Ueno, J. Chem. Phys. 119, 7054 (2003)] is tested in a supercritical Lennard-Jones fluid over a wide density range by comparing it with the Langevin dynamics and molecular dynamics simulations and other theories. The density n range studied is 0.323n(c)< or =n< or =2.58n(c) and the temperature 1.52T(c), where n(c) and T(c) are the critical density and temperature, respectively. For the rate of bimolecular reactions, the transition between the collision-limited and diffusion-limited regimes is expected to take place in this density range. The simulations show that the rate constant decays with time extensively at high densities, and that the magnitude of decay decreases gradually with decreasing density. The decay profiles of the rate constants obtained by the simulations are reproduced reasonably well by the FPKE treatment in the whole density range studied if a continuous velocity distribution is used in solving the FPKE approximately. If a discontinuous velocity distribution is used instead of the continuous one, the FPKE treatment leads to a rate constant much larger than the simulation results at medium and low densities. The rate constants calculated from the Smoluchowski-Collins-Kimball (SCK) theory based on the diffusion equation are somewhat smaller than the simulation results in medium and low densities when the intrinsic rate constant is chosen to adjust the steady state rate constant in the low density limit to that derived by the kinetic collision theory. The discrepancy is relatively small, so that the SCK theory provides a useful guideline for a qualitative discussion of the density effect on the rate constant.     

Thin-thick surface phase coexistence and boundary tension of the square-well fluid on a weak attractive surface

Prewetting transition is studied for the square-well fluid of attractive-well diameter lambda(ff)sigma(ff)=1.5 in the presence of a homogeneous surface modeled by the square-well potential of attractive well from 0.8sigma(ff) to 1.8sigma(ff). We investigate surface phase coexistence of thin-thick film transition using grand-canonical transition matrix Monte Carlo (GC-TMMC) and histogram reweighting techniques. Molecular dynamics (MD) and GC-TMMC are utilized to predict the properties of the fluid for various surface fluid affinities. Occurrences of prewetting transition with the variation of surface affinity are observed for a domain of reduced temperature from T(*)=0.62 to 0.75. We have used MD and GC-TMMC+finite size scaling (FSS) simulations to calculate the boundary tension as a function of temperature as well as surface affinity. Boundary tensions via MD and GC-TMMC+FSS methods are in good agreement. The boundary tension increases with the decrease of wall-fluid affinity. Prewetting critical properties are calculated using rectilinear diameter approach and scaling analysis. We found that critical temperature and density increase with the decrease of wall-fluid affinity.     

Ab initio and kinetic calculations for the reactions of Cl with SiH(n)Cl(4-n) (n=1,2,3,4)

The direct hydrogen abstraction reactions of Cl atom with SiH(n)Cl(4-n) (n=1,2,3,4) have been studied systematically using ab initio molecular orbital theory. Geometries have been optimized at the MP2 level with 6-311+G(d) basis set, QCISD(T)/6-311+G(d,p) has been used in the final single point energy calculation. The kinetic calculations of these reactions have been explored using the canonical variational transition (CVT) state theory method with small-curvature tunneling (SCT) effect correction over the temperature range of 200-2000 K. The CVT/SCT rate constants exhibit typical non-Arrhenius behavior and three-parameter rate-temperature formulas have been fitted for the reactions of Cl with SiH4, SiH3Cl, SiH2Cl2, and SiHCl3, respectively (in unit of cm(3) molecule(-1) s(-1)). The calculated CVT/SCT rate constants are in agreement with the available experimental values.     

Seven dimensional quantum dynamics study of the H2+NH2-->H+NH3 reaction

Initial state-selected time-dependent wave packet dynamics calculations have been performed for the H2+NH2-->H+NH3 reaction using a seven dimensional model on an analytical potential energy surface based on the one developed by Corchado and Espinosa-Garcia [J. Chem. Phys. 106, 4013 (1997)]. The model assumes that the two spectator NH bonds are fixed at their equilibrium values and nonreactive NH2 group keeps C2v symmetry and the rotation-vibration coupling in NH2 is neglected. The total reaction probabilities are calculated when the two reactants are initially at their ground states, when the NH2 bending mode is excited, and when H2 is on its first vibrational excited state, with total angular momentum J=0. The converged cross sections for the reaction are also reported for these initial states. Thermal rate constants and equilibrium constants are calculated for the temperature range of 200-2000 K and compared with transition state theory results and the available experimental data. The study shows that (a) the reaction is dominated by ground-state reactivity and the main contribution to the thermal rate constants is thought to come from this state, (b) the excitation energy of H2 was used to enhance reactivity while the excitation of the NH2 bending mode hampers the reaction, (c) the calculated thermal rate constants are very close to the experimental data and transition state theory results at high and middle temperature, while they are ten times higher than that of transition state theory at low temperature (T=200 K), and (d) the equilibrium constants results indicate that the approximations applied may have different roles in the forward and reverse reactions.