Polymer reversal rate calculated via locally scaled diffusion map

A recent study on the dynamics of polymer reversal inside a nanopore by Huang and Makarov [J. Chem. Phys. 128, 114903 (2008)] demonstrated that the reaction rate cannot be reproduced by projecting the dynamics onto a single empirical reaction coordinate, a result suggesting the dynamics of this system cannot be correctly described by using a single collective coordinate. To further investigate this possibility we have applied our recently developed multiscale framework, locally scaled diffusion map (LSDMap), to obtain collective reaction coordinates for this system. Using a single diffusion coordinate, we obtain a reversal rate via Kramers expression that is in good agreement with the exact rate obtained from the simulations. Our mathematically rigorous approach accounts for the local heterogeneity of molecular configuration space in constructing a diffusion map, from which collective coordinates emerge. We believe this approach can be applied in general to characterize complex macromolecular dynamics by providing an accurate definition of the collective coordinates associated with processes at different time scales.     

An application of coupled reference interaction site model/molecular dynamics to the conformational analysis of the alanine dipeptide

We present an application of our recently proposed coupled reference interaction site model (RISM) molecular dynamics (MD) solvation free energy methodology [Freedman and Truong, Chem. Phys. Lett. 381, 362 (2003); J. Chem. Phys. 121, 2187 (2004)] to study the conformational stability of alanine dipeptide in aqueous solution. In this methodology, radial distribution functions obtained from a single MD simulation are substituted into a RISM expression for solvation free energy. Consequently, iterative solution of the RISM equation is not needed. The relative solvation free energies of seven different conformations of the alanine dipeptide in aqueous solution are calculated. Results from the coupled RISM/MD methodology are in good agreement with those from earlier simulations using the accurate free energy perturbation approach, showing that the alphaR conformation is most stabilized by solution. This study establishes a framework for applying this coupled RISM/MD method to larger biological systems.     

The role of the excited electronic states in the C+ + H2O reaction

The electronic excited states of the [COH2]+ system have been studied in order to establish their role in the dynamics of the C+ + H2O-->[COH]+ +H reaction, which is a prototypical ion-molecule reaction. The most relevant minima and saddle points of the lowest excited state have been determined and energy profiles for the lowest excited doublet and quartet electronic states have been computed along the fragmentation and isomerization coordinates. Also, nonadiabatic coupling strengths between the ground and the first excited state have been computed where they can be large. Our analysis suggests that the first excited state could play an important role in the generation of the formyl isomer, which has been detected in crossed beam experiments [D. M. Sonnenfroh et al., J. Chem. Phys. 83, 3985 (1985)], but could not be explained in quasiclassical trajectory computations [Y. Ishikawa et al., Chem. Phys. Lett. 370, 490 (2003); J. R. Flores, J. Chem. Phys. 125, 164309 (2006)].     

Modeling real dynamics in the coarse-grained representation of condensed phase systems

This work presents a systematic multiscale methodology to provide a more faithful representation of real dynamics in coarse-grained molecular simulation models. The theoretical formalism is based on the recently developed multiscale coarse-graining (MS-CG) method [S. Izvekov and G. A. Voth, J. Phys. Chem. B. 109, 2469 (2005); J. Chem. Phys. 123, 134105 (2005)] and relies on the generalized Langevin equation approach and its simpler Langevin equation limit. The friction coefficients are determined in multiscale fashion from the underlying all-atom molecular dynamics simulations using force-velocity and velocity-velocity correlation functions for the coarse-grained sites. The diffusion properties in the resulting CG Brownian dynamics simulations are shown to be quite accurate. The time dependence of the velocity autocorrelation function is also well-reproduced relative to the all-atom model if sufficient resolution of the CG sites is implemented.     

Reliable molecular simulations of solute-solvent systems with a minimum number of solvent shells

In this work, the mean field (MF) method, a continuum-based model designed for treating complex molecular systems, such as liquids and solutions, recently presented by Brancato et al. [J. Chem. Phys. 122, 154109 (2005)], has been further developed and improved especially in the treatment of the electrostatics. The revised model has been used to investigate the size effects on several physical properties of various solute-solvent systems by increasing the number of explicitly included solvent molecules from few tens up to thousands. Results on simple ions, such as sodium and chloride ions, and on a small peptide, such as alanine dipeptide analog (AcAlaNHMe), have shown that solvation structures and dynamics, as well as solvent-induced changes in the solute conformation, can be correctly reproduced by the MF model, providing that only two or three solvent layers are treated explicitly.     

Excess electron kinetics by dispersive conductivity: A proposition

A proposition for the simultaneous experimental determination of mobility and scavenger rate coefficient of localized excess electrons is exposed. Based on an earlier theory [Schiller, R., Balog, J., Nagy, G., 2005. Continuous-time random-walk theory of interfering diffusion and chemical reaction with an application to electrochemical impedance spectra of oxidized Zr–1%Nb. J. Chem. Phys. 123, 094704-1-7], it is expected that the steady state conductivity of an insulating liquid, irradiated or illuminated in the presence of some scavenger, depends on frequency. By analyzing the frequency dependence one can obtain the values of the above material properties.     

A molecular dynamics study of water nucleation using the TIP4P/2005 model

Extensive molecular dynamics simulations were conducted using the TIP4P/2005 water model of Abascal and Vega [J. Chem. Phys. 123, 234505 (2005)] to investigate its condensation from supersaturated vapor to liquid at 330 K. The mean first passage time method [J. Wedekind, R. Strey, and D. Reguera, J. Chem. Phys. 126, 134103 (2007); L. S. Bartell and D. T. Wu, 125, 194503 (2006)] was used to analyze the influence of finite size effects, thermostats, and charged species on the nucleation dynamics. We find that the Nose?-Hoover thermostat and the one proposed by Bussi et al. [J. Chem. Phys. 126, 014101 (2007)] give essentially the same averages. We identify the maximum thermostat coupling time to guarantee proper thermostating for these simulations. The presence of charged species has a dramatic impact on the dynamics, inducing a marked change towards a pure growth regime, which highlights the importance of ions in the formation of liquid droplets in the atmosphere. It was found a small but noticeable sign preference at intermediate cluster sizes (between 5 and 30 water molecules) corresponding mostly to the formation of the second solvation shell around the ion. The TIP4P/2005 water model predicts that anions induce faster formation of water clusters than cations of the same magnitude of charge.     

Rovibrational energy levels of hydrogen peroxide, studied by MULTIMODE with a reaction path Hamiltonian

Recently, Carter and Handy [J. Chem. Phys. 113 (2000) 987] have introduced the theory of the reaction path Hamiltonian (RPH) [J. Chem. Phys. 72 (1980) 99] into the variational scheme MULTIMODE, for the calculation of the J=0 vibrational levels of polyatomic molecules, which have a single large-amplitude motion. In this theory the reaction path coordinate s becomes the fourth dimension of the moment-of-inertia tensor, and must be treated separately from the remaining 3N-7 normal coordinates in the MULTIMODE program. In the modified program, complete integration is performed over s, and the M-mode MULTIMODE coupling approximation for the evaluation of the matrix elements applies only to the 3N-7 normal coordinates. In this paper the new algorithm is extended to the calculation of rotational-vibration energy levels (i.e. J>0) with the RPH, following from our analogous implementation for rigid molecules [Theoret. Chem. Acc. 100 (1998) 191]. The full theory is given, and all extra terms have been included to give the exact kinetic energy operator. In order to validate the new code, we report studies on hydrogen peroxide (H2O2), where the reaction path is equivalent to torsional motion. H2O2 has previously been studied variationally using a valence coordinate Hamiltonian; complete agreement for calculated rovibrational levels is obtained between the previous results and those from the new code, using the identical potential surface. MULTIMODE is now able to calculate rovibrational levels for polyatomic molecules which have one large-amplitude motion.     

Towards an understanding of many-particle effects in hydrophobic association in methane solutions

This paper applies the multiscale coarse-graining method [S. Izvekov and G. A. Voth, J. Phys. Chem. B 109, 2469 (2005); J. Chem. Phys. 123, 134105 (2005)] to analyze many-body effects in concentrated methane solutions. Pairwise decompositions of N-particle solute-solute potentials of mean force (PMFs), and the respective solvent cavity potentials, enthalpic, entropic, and heat capacity of hydrophobic association, are calculated directly from unconstrained molecular-dynamics simulations of methane solutions at different molar fractions, with the highest being 0.055. The many-body effects in hydrophobic hydration are further studied using N-methane PMFs, which are explicitly dependent on solvent coordinates.     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.     

The use of atomic intrinsic polarizabilities in the evaluation of the dispersion energy

The recent approach presented by Becke and Johnson [J. Chem. Phys. 122, 154104 (2005); 123, 024101 (2005); 123, 154101 (2005); 124, 174104 (2006); 124, 014104 (2006)] for the evaluation of dispersion interactions based on the properties of the exchange-hole dipole moment is combined with a Hirshfeld-type partitioning for the molecular polarizabilities into atomic contributions, recently presented by some of the present authors [A. Krishtal et al., J. Chem. Phys. 125, 034312 (2006)]. The results on a series of nine dimers, involving neon, methane, ethene, acetylene, benzene, and CO(2), taken at their equilibrium geometry, indicate that when the C(6), C(8), and C(10) terms are taken into account, the resulting dispersion energies can be obtained deviating 3% or 8% from high level literature data [E. R. Johnson and A. D. Becke, J. Chem. Phys. 124, 174104 (2006)], without the use of a damping function, the only outlier being the parallel face-to-face benzene dimer.     

Electrostatics of proteins in dielectric solvent continua. II. First applications in molecular dynamics simulations

In the preceding paper by Stork and Tavan, [J. Chem. Phys. 126, 165105 (2007)], the authors have reformulated an electrostatic theory which treats proteins surrounded by dielectric solvent continua and approximately solves the associated Poisson equation [B. Egwolf and P. Tavan, J. Chem. Phys. 118, 2039 (2003)]. The resulting solution comprises analytical expressions for the electrostatic reaction field (RF) and potential, which are generated within the protein by the polarization of the surrounding continuum. Here the field and potential are represented in terms of Gaussian RF dipole densities localized at the protein atoms. Quite like in a polarizable force field, also the RF dipole at a given protein atom is induced by the partial charges and RF dipoles at the other atoms. Based on the reformulated theory, the authors have suggested expressions for the RF forces, which obey Newton's third law. Previous continuum approaches, which were also built on solutions of the Poisson equation, used to violate the reactio principle required by this law, and thus were inapplicable to molecular dynamics (MD) simulations. In this paper, the authors suggest a set of techniques by which one can surmount the few remaining hurdles still hampering the application of the theory to MD simulations of soluble proteins and peptides. These techniques comprise the treatment of the RF dipoles within an extended Lagrangian approach and the optimization of the atomic RF polarizabilities. Using the well-studied conformational dynamics of alanine dipeptide as the simplest example, the authors demonstrate the remarkable accuracy and efficiency of the resulting RF-MD approach.     

Semiclassical calculation of nonadiabatic thermal rate constants: application to condensed phase reactions

The nonadiabatic transition state theory proposed recently by Zhao et al. [J. Chem. Phys. 121, 8854 (2004)] is extended to calculate rate constants of complex systems by using the Monte Carlo and umbrella sampling methods. Surface hopping molecular dynamics technique is incorporated to take into account the dynamic recrossing effect. A nontrivial benchmark model of the nonadiabatic reaction in the condensed phase is used for the numerical test. It is found that our semiclassical results agree well with those produced by the rigorous quantum mechanical method. Comparing with available analytical approaches, we find that the simple statistical theory proposed by Straub and Berne [J. Chem. Phys. 87, 6111 (1987)] is applicable for a wide friction region although their formula is obtained using Landau-Zener [Phys. Z. Sowjetunion 2, 46 (1932); Proc. R. Soc. London, Ser. A 137, 696 (1932)] nonadiabatic transition probability along a one-dimensional diffusive coordinate. We also investigate how the nuclear tunneling events affect the dependence of the rate constant on the friction.     

Projection of two-dimensional diffusion in a curved midline and narrow varying width channel onto the longitudinal dimension

This study focuses on the derivation of a general effective diffusion coefficient to describe the two-dimensional (2D) diffusion in a narrow and smoothly asymmetric channel of varying width, in the simple diffusional motion of noninteracting pointlike particles under no external field. We present a generalization to the case of an asymmetric channel using the projection method introduced earlier by Kalinay and Percus [J. Chem. Phys. 122, 204701 (2005); and Phys. Rev. E 74, 041203 (2006)] to project the 2D diffusion equation into an effective one-dimensional generalized Fick-Jacobs equation. The expression for the diffusion coefficient given in Eq. (23) is our main result. This expression is a more general effective diffusion coefficient for narrow channels in 2D, which contains the well-known previous results as special cases, namely, those obtained by Bradley [Phys. Rev. E 80, 061142 (2009)], and more recently by Berezhkovskii and Szabo [J. Chem. Phys. 135, 074108 (2011)]. Finally, we study some specific 2D asymmetric channel configurations to test and show the broader applicability of Eq. (23).     

Use of a nondirect-product basis for treating singularities in triatomic rotational-vibrational calculations

A technique has been developed which in principle allows the determination of the full rotational-vibrational eigenspectrum of triatomic molecules by treating the important singularities present in the triatomic rotational-vibrational kinetic energy operator given in Jacobi coordinates and the R(1) embedding. The singular term related to the diatom-type coordinate, R(1), deemed to be unimportant for spectroscopic applications, is given no special attention. The work extends a previous [J. Chem. Phys., 2005, 122, 024101] vibration-only approach and employs a generalized finite basis representation (GFBR) resulting in a nonsymmetric Hamiltonian matrix [J. Chem. Phys., 2006, 124, 014110]. The basis set to be used is obtained by taking the direct product of a 1-D DVR basis, related to R(1), with a 5-D nondirect-product basis, the latter formed by coupling Bessel-DVR functions depending on the distance-type coordinate causing the singularity, associated Legendre polynomials depending on the Jacobi angle, and rotational functions depending on the three Euler angles. The robust implicitly restarted Arnoldi method within the ARPACK package is used for the determination of a number of eigenvalues of the nonsymmetric Hamiltonian matrix. The suitability of the proposed approach is shown by the determination of the rotational-vibrational energy levels of the ground electronic state of H(3)(+) somewhat above its barrier to linearity. Convergence of the eigenenergies is checked by an alternative approach, employing a Hamiltonian expressed in Radau coordinates, a standard direct-product basis, and no treatment of the singularities.     

Lattice-Boltzmann simulations of the dynamics of polymer solutions in periodic and confined geometries

A numerical method to simulate the dynamics of polymer solutions in confined geometries has been implemented and tested. The method combines a fluctuating lattice-Boltzmann model of the solvent [Ladd, Phys. Rev. Lett. 70, 1339 (1993)] with a point-particle model of the polymer chains. A friction term couples the monomers to the fluid [Ahlrichs and Dunweg, J. Chem. Phys. 111, 8225 (1999)], providing both the hydrodynamic interactions between the monomers and the correlated random forces. The coupled equations for particles and fluid are solved on an inertial time scale, which proves to be surprisingly simple and efficient, avoiding the costly linear algebra associated with Brownian dynamics. Complex confined geometries can be represented by a straightforward mapping of the boundary surfaces onto a regular three-dimensional grid. The hydrodynamic interactions between monomers are shown to compare well with solutions of the Stokes equations down to distances of the order of the grid spacing. Numerical results are presented for the radius of gyration, end-to-end distance, and diffusion coefficient of an isolated polymer chain, ranging from 16 to 1024 monomers in length. The simulations are in excellent agreement with renormalization group calculations for an excluded volume chain. We show that hydrodynamic interactions in large polymers can be systematically coarse-grained to substantially reduce the computational cost of the simulation. Finally, we examine the effects of confinement and flow on the polymer distribution and diffusion constant in a narrow channel. Our results support the qualitative conclusions of recent Brownian dynamics simulations of confined polymers [Jendrejack et al., J. Chem. Phys. 119, 1165 (2003) and Jendrejack et al., J. Chem. Phys. 120, 2513 (2004)].     

State-to-state reaction probabilities for the H+O2(v,j)-->O+OH(v',j') reaction on three potential energy surfaces

We report state-to-state and total reaction probabilities for J=0 and total reaction probabilities for J=2 and 4 for the title reaction, both for ground-state and initially rovibrationally excited reactants. The results for three different potential energy surfaces are compared and contrasted. The potential energy surfaces employed are the DMBE IV surface by Pastrana et al. [J. Phys. Chem. 94, 8073 (1990)], the surface by Troe and Ushakov (TU) [J. Chem. Phys. 115, 3621 (2001)], and the new XXZLG ab initio surface by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. Our results show that the total reaction probabilities from both the TU and XXZLG surfaces are much smaller in magnitude for collision energies above 1.2 eV compared to the DMBE IV surface. The three surfaces also show different behavior with regards to the effect of initial state excitation. The reactivity is increased on the XXZLG and the TU surfaces and decreased on the DMBE IV surface. Vibrational and rotational product state distributions for the XXZLG and the DMBE IV surface show different behaviors for both types of distributions. Our results show that for energies above 1.25 eV the dynamics on the DMBE IV surface are not statistical. However, there is also evidence that the dynamics on the XXZLG surface are not purely statistical for energies above the onset of the first excited product vibrational state v'=1. The magnitude of the total reaction probability is decreased for J>0 for the DMBE IV and the XXZLG surfaces for ground-state reactants. However, for initially rovibrationally excited reactants, the total reaction probability does not decrease as expected for both surfaces. As a result the total cross section averaged over all Boltzmann accessible rotational states may well be larger than the cross section reported in the literature for j=1.     

Breakdown of the pseudopotential approximation for magnetizabilities and electric multipole moments. II. The importance of gauge invariance for large-core semi-local pseudopotentials

In a previous paper we pointed out that core contributions to the static magnetizability are non-negligible, and can therefore lead to erroneous results within the pseudopotential approximation [P. Schwerdtfeger, B. Assadollahzadeh, U. Rohrmann, R. Scha?fer, and J. R. Cheeseman, J. Chem. Phys. 134, 204102 (2011)]. In a recent paper van Wu?llen showed that additional terms arising from the gauge-invariant condition for the semi-local part of a pseudopotential operator can lead to non-negligible contributions to the magnetizability tensor, which are sensitive to the gauge origin and basis set chosen [C. van Wu?llen, J. Chem. Phys. 136, 114110 (2012)]. These terms were neglected in previous calculations as they were assumed to be small. In this paper we analyze the importance of the gauge-dependent semi-local pseudopotential correction term in detail for AuF and clusters of Sn showing that it leads indeed to very large corrections to the paramagnetic term for large-core pseudopotentials. Without this correction the results become very sensitive to the basis set applied. This now resolves some of the unusual large paramagnetic contributions reported before for both AuF and Sn(2).