A semiclassical reactive flux method for the calculation of condensed phase activated rate constants

A semi classical reactive flux algorithm for calculating thermally activated rate constants is presented which is based on a semi-classical transition state theory due to Chapman, Garrett and Miller [J. Chem. Phys. 63 (1975) 2710]. This reactive flux technique, when combined with the semiclassical TST, enables one to describe dynamical recrossings of the transition state on the same footing as tunneling effects. Most importantly, the method is readily applied to nonlinear multidimensional systems over a wide range of temperatures. It will be shown that the method works very well for a variety of existing models.     

Effective force field for liquid hydrogen fluoride from ab initio molecular dynamics simulation using the force-matching method

A recently developed force-matching method for obtaining effective force fields for condensed matter systems from ab initio molecular dynamics (MD) simulations has been applied to fit a simple nonpolarizable two-site pairwise force field for liquid hydrogen fluoride. The ab initio MD in this case was a Car-Parrinello (CP) MD simulation of 64 HF molecules at nearly ambient conditions within the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure included a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The performance of the force-match potential was examined for the gas-phase dimer and for the liquid phase at various temperatures. The model was able to reproduce correctly the bent structure and energetics of the gas-phase dimer, while the results for the structural properties, self-diffusion, vibrational spectra, density, and thermodynamic properties of liquid HF were compared to both experiment and the CP MD simulation. The force-matching model performs well in reproducing nearly all of the liquid properties as well as the aggregation behavior at different temperatures. The model is computationally cheap and compares favorably to many more computationally expensive potential energy functions for liquid HF.     

Vibrational energy relaxation dynamics of Si---H stretching modes on stepped H/Si(111) 1×1 surfaces

The vibrational energy relaxation rates of excited Si---H stretching modes on the monohydride steps of miscut H/Si(111) 1×1 surfaces are calculated using Bloch-Redfield theory combined with classical molecular dynamics (MD) simulation. The structure and vibrational frequencies of the surface are first investigated using the Car-Parrinello ab initio MD method. The calculated Si---Si---H bending frequencies and relaxed structures are then used to refine the empirical potential for the classical MD simulations. The lifetime of the excited Si---H stretching mode at the step is found to be shorter than the modes on the terrace. Both the magnitude and the trend of the calculated results agree well with the experimental measurement on the 9° monohydride stepped surface. The vibrational relaxation rate of the Si---H stretching modes on the 15° monohydride stepped surface are also calculated and predicted to have a slightly shorter lifetime than for the 9° surface.     

Ab initio molecular-dynamics simulation of aqueous proton solvation and transport revisited

The solvation and transport of the hydrated excess proton is studied using the Car-Parrinello molecular-dynamics (CPMD) simulation method. The simulations were performed using BLYP and HCTH gradient-corrected exchange-correlation energy functionals. The fictitious electronic mass was chosen to be small enough so that the underlying water structural and dynamical properties were converged with respect to this important CPMD simulation parameter. An unphysical overstructuring of liquid water in the CPMD simulations using the BLYP functional resulted in the formation of long-lived hydrogen-bonding structures involving the excess proton and a particular (special) water oxygen. The excess proton was observed to be attracted to the special oxygen through the entire length of the BLYP CPMD simulations. Consequently, the excess proton diffusion was limited by the mobility of the special oxygen in the slowly diffusing water network and, in turn, the excess proton self-diffusion coefficient was found to be significantly below the experimental value. On the other hand, the structural properties of liquid water in the HCTH CPMD simulation were seen to be in better agreement with experiment, although the water and excess proton diffusions were still well below the experimental value.     

Probing the molecular-scale lipid bilayer response to shear flow using nonequilibrium molecular dynamics

A nonequilibrium molecular dynamics simulation of the response of dimyristoylphosphatidylcholine (DMPC) bilayers to a solvent shear flow is presented. Application of shear flow to planar, stationary DMPC bilayers results in a redistribution of the membrane density profile along the bilayer normal due to the alignment of the lipids in the direction of flow and an increase in average lipid chain length. An increase in the intermolecular and intramolecular order of the lipids in response to the shear flow is also observed. This study provides groundwork for understanding the mechanism of the full response of lipid bilayers to externally imposed solvent shear flows, beginning with the response in the absence of collective lipid motions such as undulations and bilayer flow.     

Precise atomic masses of147Eu,147Gd,and151Tb derived from their decay properties

Theβ-decay energies of¹⁴⁷Eu,¹⁴⁷Gd, and¹⁵¹Tb were determined by usingγ-spectroscopical methods. The comparison of experimental with calculatedK-capture probabilities yielded theQ _EC values 1.690( _?16 ⁺²¹ )MeV and 2.203( _?13 ⁺¹⁹ )MeV for¹⁴⁷Eu and¹⁴⁷Gd, respectively. By measuring the ratio of positron decay to electron capture for two branches in¹⁴⁷Eu decay, the decay energiesQ _EC=1.702(13) MeV andQ _EC=1.709(18)MeV were derived. Also fromEC/β ⁺ ratios the valuesQ _EC=2.225(75) MeV for¹⁴⁷Gd, andQ _EC=2.566(12)MeV for¹⁵¹Tb were obtained. Earlier discrepancies in the mass adjustment of these isotopes were removed. In course of the present studiesγ-decay properties of¹⁴⁷Eu and¹⁴⁷Gd were reinvestigated.     

Unique spatial heterogeneity in ionic liquids

A multiscale coarse-graining model for ionic liquids has been extended to investigate the unique aggregation of cations in ionic liquids through computer simulation. It has been found that, with sufficiently long side chains, the tail groups of cations aggregate to form spatially heterogeneous domains, while headgroups of the cations and the anions distribute as uniformly as possible. This is understood as the result of competition between the charged electrostatic interactions between headgroups and anions and the collective short-range interactions between the neutral tail groups. This aggregation can help to explain a number of experimentally observed physical phenomena in ionic liquids.