Search

Coarse-grained (CG) molecular dynamics (MD) simulations have become popular for investigating systems on multiple length and time scales ranging from atomistic to mesoscales. In CGMD, several atoms are mapped onto a single CG bead and the effective interactions between CG beads are determined. Iterative coarse-graining methods, such as iterative Boltzmann inversion (IBI), are computationally expensive and can have convergence issues. In this paper, we present a direct and computationally efficient theoretical procedure for coarse-graining based on the Ornstein-Zernike (OZ) and hypernetted chain (HNC) integral equation theory. We demonstrate the OZ-HNC-based CG method by coarse-graining a bulk water system, a water-methanol mixture system, and an electrolyte system. We show that the accuracy of the CG potentials obtained from the OZ-HNC-based coarse-graining is comparable to iterative systematic coarse-graining methods. Furthermore, we show that the CG potentials from OZ-HNC can be used to reduce the number of iterations and hence the computational cost of the iterative systematic coarse-graining approaches, like IBI and relative entropy minimization. [ABSTRACT FROM AUTHOR]

In the preceding paper [I. Fukuda, J. Chem. Phys. 139, 174107 (2013)], the zero-multipole (ZM) summation method was proposed for efficiently evaluating the electrostatic Coulombic interactions of a classical point charge system. The summation takes a simple pairwise form, but prevents the electrically non-neutral multipole states that may artificially be generated by a simple cutoff truncation, which often causes large energetic noises and significant artifacts. The purpose of this paper is to judge the ability of the ZMmethod by investigating the accuracy, parameter dependencies, and stability in applications to liquid systems. To conduct this, first, the energy-functional error was divided into three terms and each term was analyzed by a theoretical error-bound estimation. This estimation gave us a clear basis of the discussions on the numerical investigations. It also gave a new viewpoint between the excess energy error and the damping effect by the damping parameter. Second, with the aid of these analyses, the ZM method was evaluated based on molecular dynamics (MD) simulations of two fundamental liquid systems, a molten sodium-chlorine ion system and a pure water molecule system. In the ion system, the energy accuracy, compared with the Ewald summation, was better for a larger value of multipole moment l currently induced until l ≲ 3 on average. This accuracy improvement with increasing l is due to the enhancement of the excess-energy accuracy. However, this improvement is wholly effective in the total accuracy if the theoretical moment l is smaller than or equal to a system intrinsic moment L. The simulation results thus indicate L ∼ 3 in this system, and we observed less accuracy in l = 4. We demonstrated the origins of parameter dependencies appearing in the crossing behavior and the oscillations of the energy error curves. With raising the moment l we observed, smaller values of the damping parameter provided more accurate results and smoother behaviors with respect to cutoff length were obtained. These features can be explained, on the basis of the theoretical error analyses, such that the excess energy accuracy is improved with increasing l and that the total accuracy improvement within l ⩽ L is facilitated by a small damping parameter. Although the accuracy was fundamentally similar to the ion system, the bulk water system exhibited distinguishable quantitative behaviors. A smaller damping parameter was effective in all the practical cutoff distance, and this fact can be interpreted by the reduction of the excess subset. A lower moment was advantageous in the energy accuracy, where l = 1 was slightly superior to l = 2 in this system. However, the method with l = 2 (viz., the zero-quadrupole sum) gave accurate results for the radial distribution function. We confirmed the stability in the numerical integration for MD simulations employing the ZM scheme. This result is supported by the sufficient smoothness of the energy function. Along with the smoothness, the pairwise feature and the allowance of the atom-based cutoff mode on the energy formula lead to the exact zero total-force, ensuring the total-momentum conservations for typical MD equations of motion. [ABSTRACT FROM AUTHOR]

This paper deals with molecular simulation of the chemical potentials in aqueous electrolyte solutions for the water solvent and its relationship to chemical potential simulation results for the electrolyte solute. We use the Gibbs-Duhem equation linking the concentration dependence of these quantities to test the thermodynamic consistency of separate calculations of each quantity. We consider aqueous NaCl solutions at ambient conditions, using the standard SPC/E force field for water and the Joung-Cheatham force field for the electrolyte. We calculate the water chemical potential using the osmotic ensemble Monte Carlo algorithm by varying the number of water molecules at a constant amount of solute. We demonstrate numerical consistency of these results in terms of the Gibbs-Duhem equation in conjunction with our previous calculations of the electrolyte chemical potential. We present the chemical potential vs molality curves for both solvent and solute in the form of appropriately chosen analytical equations fitted to the simulation data. As a byproduct, in the context of the force fields considered, we also obtain values for the Henry convention standard molar chemical potential for aqueous NaCl using molality as the concentration variable and for the chemical potential of pure SPC/E water. These values are in reasonable agreement with the experimental values. [ABSTRACT FROM AUTHOR]

In this paper, the second in a series devoted to molecular modeling of protein aggregation, a mesoscale model of proteins together with extensive discontinuous molecular dynamics simulation is used to study the phenomenon in a confined medium. The medium, as a model of a crowded cellular environment, is represented by a spherical cavity, as well as cylindrical tubes with two aspect ratios. The aggregation process leads to the formation of β sheets and eventually fibrils, whose deposition on biological tissues is believed to be a major factor contributing to many neuro-degenerative diseases, such as Alzheimer’s, Parkinson’s, and amyotrophic lateral sclerosis diseases. Several important properties of the aggregation process, including dynamic evolution of the total number of the aggregates, the mean aggregate size, and the number of peptides that contribute to the formation of the β sheets, have been computed. We show, similar to the unconfined media studied in Paper I [S. Zheng et al., J. Chem. Phys. 145, 134306 (2016)], that the computed properties follow dynamic scaling, characterized by power laws. The existence of such dynamic scaling in unconfined media was recently confirmed by experiments. The exponents that characterize the power-law dependence on time of the properties of the aggregation process in spherical cavities are shown to agree with those in unbounded fluids at the same protein density, while the exponents for aggregation in the cylindrical tubes exhibit sensitivity to the geometry of the system. The effects of the number of amino acids in the protein, as well as the size of the confined media, have also been studied. Similarities and differences between aggregation in confined and unconfined media are described, including the possibility of no fibril formation, if confinement is severe. [ABSTRACT FROM AUTHOR]

This article extends the uncertainty quantification analysis introduced in Paper I for molecular dynamics (MD) simulations of concentration driven ionic flow through a silica nanopore. Attention is now focused on characterizing, for a fixed pore diameter of D = 21 Å, the sensitivity of the system to the Lennard-Jones energy parameters, &eh;Na+ and &eh;Cl-, defining the depth of the potential well for the two ions Na+ and Cl-, respectively. A forward propagation analysis is applied to map the uncertainty in these parameters to the MD predictions of the ionic fluxes. Polynomial chaos expansions and Bayesian inference are exploited to isolate the effect of the intrinsic noise, stemming from thermal fluctuations of the atoms, and properly quantify the impact of parametric uncertainty on the target MD predictions. A Bayes factor analysis is then used to determine the most suitable regression model to represent the MD noisy data. The study shows that the response surface of the Na+ conductance can be effectively inferred despite the substantial noise level, whereas the noise partially hides the underlying trend in the Cl- conductance data over the studied range. Finally, the dependence of the conductances on the uncertain potential parameters is analyzed in terms of correlations with key bulk transport coefficients, namely, viscosity and collective diffusivities, computed using Green-Kubo time correlations. [ABSTRACT FROM AUTHOR]

Enhanced sampling techniques that target a set of collective variables and that use molecular dynamics as the driving engine have seen widespread application in the computational molecular sciences as a means to explore the free-energy landscapes of complex systems. The use of molecular dynamics as the fundamental driver of the sampling requires the introduction of a time step whose magnitude is limited by the fastest motions in a system. While standard multiple time-stepping methods allow larger time steps to be employed for the slower and computationally more expensive forces, the maximum achievable increase in time step is limited by resonance phenomena, which inextricably couple fast and slow motions. Recently, we introduced deterministic and stochastic resonance-free multiple time step algorithms for molecular dynamics that solve this resonance problem and allow ten- to twenty-fold gains in the large time step compared to standard multiple time step algorithms [P. Minary et al., Phys. Rev. Lett. 93, 150201 (2004); B. Leimkuhler et al., Mol. Phys. 111, 3579-3594 (2013)]. These methods are based on the imposition of isokinetic constraints that couple the physical system to Nosé-Hoover chains or Nos'e-Hoover Langevin schemes. In this paper, we show how to adapt these methods for collective variable-based enhanced sampling techniques, specifically adiabatic free-energy dynamics/temperature-accelerated molecular dynamics, unified free-energy dynamics, and by extension, metadynamics, thus allowing simulations employing these methods to employ similarly very large time steps. The combination of resonance-free multiple time step integrators with free-energy-based enhanced sampling significantly improves the efficiency of conformational exploration. [ABSTRACT FROM AUTHOR]

Numerous molecular systems, including solutions, proteins, and composite materials, can be modeled using mixed-resolution representations, of which the quantum mechanics/molecular mechanics (QM/MM) approach has become the most widely used. However, the QM/MM approach often faces a number of challenges, including the high cost of repetitive QM computations, the slow sampling even for the MM part in those cases where a system under investigation has a complex dynamics, and a difficulty in providing a simple, qualitative interpretation of numerical results in terms of the influence of the molecular environment upon the active QM region. In this paper, we address these issues by combining QM/MM modeling with the methodology of "bottom-up" coarse-graining (CG) to provide the theoretical basis for a systematic quantum-mechanical/coarse-grained molecular mechanics (QM/CG-MM) mixed resolution approach. A derivation of the method is presented based on a combination of statistical mechanics and quantum mechanics, leading to an equation for the effective Hamiltonian of the QM part, a central concept in the QM/CG-MM theory. A detailed analysis of different contributions to the effective Hamiltonian from electrostatic, induction, dispersion, and exchange interactions between the QM part and the surroundings is provided, serving as a foundation for a potential hierarchy of QM/CG-MM methods varying in their accuracy and computational cost. A relationship of the QM/CG-MM methodology to other mixed resolution approaches is also discussed. [ABSTRACT FROM AUTHOR]

Conventional molecular-dynamics (cMD) simulation has a well-known limitation in accessible time and length scales, and thus various enhanced sampling techniques have been proposed to alleviate the problem. In this paper, we explore the utility of replica exchange with solute tempering (REST) (i.e., a variant of Hamiltonian replica exchange methods) to simulate the self-assembly of a supramolecular polymer in explicit solvent and compare the performance with temperaturebased replica exchange MD (T-REMD) as well as cMD. As a test system, we consider a relatively simple all-atom model of supramolecular polymerization (namely, benzene-1,3,5-tricarboxamides in methylcyclohexane solvent). Our results show that both REST and T-REMD are able to predict highly ordered polymer structures with helical H-bonding patterns, in contrast to cMD which completely fails to obtain such a structure for the present model. At the same time, we have also experienced some technical challenge (i.e., aggregation-dispersion transition and the resulting bottleneck for replica traversal), which is illustrated numerically. Since the computational cost of REST scales more moderately than T-REMD, we expect that REST will be useful for studying the self-assembly of larger systems in solution with enhanced rearrangement of monomers. [ABSTRACT FROM AUTHOR]

As proteins contain both hydrophobic and hydrophilic amino acids, they will readily adsorb onto interfaces between water and hydrophobic fluids such as oil. This adsorption normally causes changes in the protein structure, which can result in loss of protein function and irreversible adsorption, leading to the formation of protein interfacial films. While this can be advantageous in some applications (e.g., food technology), in most cases it limits our ability to exploit protein functionality at interfaces. To understand and control protein interfacial adsorption and function, it is necessary to understand the microscopic conformation of proteins at liquid interfaces. In this paper, molecular dynamics simulations are used to investigate the adsorption and conformation of two similar proteins, lysozyme and a-lactalbumin, at a water-octane interface. While they both adsorb onto the interface, a-lactalbumin does so in a specific orientation, mediated by two amphipathic helices, while lysozyme adsorbs in a non-specific manner. Using replica exchange simulations, both proteins are found to possess a number of distinct interfacial conformations, with compact states similar to the solution conformation being most common for both proteins. Decomposing the different contributions to the protein energy at oil-water interfaces suggests that conformational change for a-lactalbumin, unlike lysozyme, is driven by favourable protein-oil interactions. Revealing these differences between the factors that govern the conformational change at interfaces in otherwise similar proteins can give insight into the control of protein interfacial adsorption, aggregation, and function. [ABSTRACT FROM AUTHOR]

Molecular Dynamics-Green's Function Reaction Dynamics (MD-GFRD) is a multiscale simulation method for particle dynamics or particle-based reaction-diffusion dynamics that is suited for systems involving low particle densities. Particles in a low-density region are just diffusing and not interacting. In this case, one can avoid the costly integration of microscopic equations of motion, such as molecular dynamics (MD), and instead turn to an event-based scheme in which the times to the next particle interaction and the new particle positions at that time can be sampled. At high (local) concentrations, however, e.g., when particles are interacting in a nontrivial way, particle positions must still be updated with small time steps of the microscopic dynamical equations. The efficiency of a multiscale simulation that uses these two schemes largely depends on the coupling between them and the decisions when to switch between the two scales. Here we present an efficient scheme for multi-scale MD-GFRDsimulations. It has been shownthatMD-GFRDschemes are more efficient than brute-force molecular dynamics simulations up to a molar concentration of 10² µM. In this paper, we showthat the choice of the propagation domains has a relevant impact on the computational performance. Domains are constructed using a local optimization of their sizes and a minimal domain size is proposed. The algorithm is shown to be more efficient than brute-force Brownian dynamics simulations up to a molar concentration of 10³ µM and is up to an order of magnitude more efficient compared with previous MD-GFRD schemes. [ABSTRACT FROM AUTHOR]

Ab initio quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulation is a useful tool to calculate thermodynamic properties such as potential of mean force for chemical reactions but intensely time consuming. In this paper, we developed a new method using the internal force correction for low-level semiempirical QM/MM molecular dynamics samplings with a predefined reaction coordinate. As a correction term, the internal force was predicted with a machine learning scheme, which provides a sophisticated force field, and added to the atomic forces on the reaction coordinate related atoms at each integration step. We applied this method to two reactions in aqueous solution and reproduced potentials of mean force at the ab initio QM/MM level. The saving in computational cost is about 2 orders of magnitude. The present work reveals great potentials for machine learning in QM/MM simulations to study complex chemical processes. [ABSTRACT FROM AUTHOR]

Fermi operator expansion (FOE) methods are powerful alternatives to diagonalization type methods for solving Kohn-Sham density functional theory (KSDFT). One example is the pole expansion and selected inversion (PEXSI) method, which approximates the Fermi operator by rational matrix functions and reduces the computational complexity to at most quadratic scaling for solving KSDFT. Unlike diagonalization type methods, the chemical potential often cannot be directly read off from the result of a single step of evaluation of the Fermi operator. Hence multiple evaluations are needed to be sequentially performed to compute the chemical potential to ensure the correct number of electrons within a given tolerance. This hinders the performance of FOE methods in practice. In this paper, we develop an efficient and robust strategy to determine the chemical potential in the context of the PEXSI method. The main idea of the new method is not to find the exact chemical potential at each self-consistent-field (SCF) iteration but to dynamically and rigorously update the upper and lower bounds for the true chemical potential, so that the chemical potential reaches its convergence along the SCF iteration. Instead of evaluating the Fermi operator for multiple times sequentially, our method uses a two-level strategy that evaluates the Fermi operators in parallel. In the regime of full parallelization, the wall clock time of each SCF iteration is always close to the time for one single evaluation of the Fermi operator, even when the initial guess is far away from the converged solution. We demonstrate the effectiveness of the new method using examples with metallic and insulating characters, as well as results from ab initio molecular dynamics. [ABSTRACT FROM AUTHOR]

The manipulation of spin states at the molecular scale is of fundamental importance for the development of molecular spintronic devices. One of the feasible approaches for the modification of a molecular spin state is through the adsorption of certain specific atoms or molecules including H, NO, CO, NH3, and O2. In this paper, we demonstrate that the local spin state of an individual iron phthalocyanine (FePc) molecule adsorbed on an Au(111) surface exhibits controllable switching by hydrogen adsorption, as evidenced by using first-principles calculations based on density functional theory. Our theoretical calculations indicate that different numbers of hydrogen adsorbed at the pyridinic N sites of the FePc molecule largely modify the structural and electronic properties of the FePc/Au(111) composite by forming extra N-H bonds. In particular, the adsorption of one or up to three hydrogen atoms induces a redistribution of charge (spin) density within the FePc molecule, and hence a switching to a lowspin state (S = 1/2) from an intermediate spin state (S = 1) is achieved, while the adsorption of four hydrogen atoms distorts the molecular conformation by increasing Fe-N bond lengths in FePc and thus breaks the ligand field exerted on the Fe 3d orbitals via stronger hybridization with the substrate, leading to an opposite switching to a high-spin state (S = 2). These findings obtained from the theoretical simulations could be useful for experimental manipulation or design of single-molecule spintronic devices. [ABSTRACT FROM AUTHOR]

Sets of atoms collectively behaving as rigid bodies are often used in molecular dynamics to model entire molecules or parts thereof. This is a coarse-graining strategy that eliminates degrees of freedom and supposedly admits larger time steps without abandoning the atomistic character of a model. In this paper, we rely on a particular factorization of the rotation matrix to simplify the mechanical formulation of systems containing rigid bodies. We then propose a new derivation for the exact solution of torque-free rotations, which are employed as part of a symplectic numerical integration scheme for rigid-body dynamics. We also review methods for calculating pressure in systems of rigid bodies with pairwise-additive potentials and periodic boundary conditions. Finally, simulations of liquid phases, with special focus on water, are employed to analyze the numerical aspects of the proposed methodology. Our results show that energy drift is avoided for time step sizes up to 5 fs, but only if a proper smoothing is applied to the interatomic potentials. Despite this, the effects of discretization errors are relevant, even for smaller time steps. These errors induce, for instance, a systematic failure of the expected equipartition of kinetic energy between translational and rotational degrees of freedom. [ABSTRACT FROM AUTHOR]

In this paper we perform molecular dynamics simulations of a spherical polyelectrolyte brush and counterions in a salt-free medium. The dielectric discontinuity on the grafted nanoparticle surface is taken into account by the method of image charges. Properties of the polyelectrolyte brush are obtained for different parameters, including valency of the counterions, radius of the nanoparticle, and the brush total charge. The monovalent counterions density profiles are obtained and compared with a simple mean-field theoretical approach. The theory allows us to obtain osmotic properties of the system. [ABSTRACT FROM AUTHOR]

The shear thinning of a lubricant significantly affects lubrication film generation at high shear rates. The critical shear rate, defined at the onset of shear thinning, marks the transition of lubricant behaviors. It is challenging to capture the entire shear-thinning curve by means of molecular dynamics (MD) simulations owing to the low signal-to-noise ratio or long calculation time at comparatively low shear rates (104-106 s-1), which is likely coincident with the shear rates of interest for lubrication applications. This paper proposes an approach that correlates the shear-thinning phenomenon with the change in the molecular conformation characterized by the radius of gyration of the molecule. Such a correlation should be feasible to capture the major mechanism of shear thinning for smallto moderate-sized non-spherical molecules, which is shear-induced molecular alignment. The idea is demonstrated by analyzing the critical shear rate for squalane (C30H62) and 1-decene trimer (C30H62); it is then implemented to study the behaviors of different molecular weight poly-α-olefin (PAO) structures. Time-temperature-pressure superpositioning (TTPS) is demonstrated and it helps further extend the ranges of the temperature and pressure for shear-thinning behavior analyses. The research leads to a relationship between molecular weight and critical shear rate for PAO structures, and the results are compared with those from the Einstein-Debye equation. [ABSTRACT FROM AUTHOR]

We compute four-dimensional diabatic potential energy surfaces and transition dipole moment surfaces of O2-O2, relevant for the theoretical description of collision-induced absorption in the forbidden X3∑-g → a1Δg and X3∑-g →b1∑+g bands at 7883 cm-1 and 13 122 cm-1, respectively. We compute potentials at the multi-reference configuration interaction (MRCI) level and dipole surfaces at the MRCI and complete active space self-consistent field (CASSCF) levels of theory. Potentials and dipole surfaces are transformed to a diabatic basis using a recent multiple-property-based diabatization algorithm. We discuss the angular expansion of these surfaces, derive the symmetry constraints on the expansion coefficients, and present working equations for determining the expansion coefficients by numerical integration over the angles. We also present an interpolation scheme with exponential extrapolation to both short and large separations, which is used for representing the O2-O2 distance dependence of the angular expansion coefficients. For the triplet ground state of the complex, the potential energy surface is in reasonable agreement with previous calculations, whereas global excited state potentials are reported here for the first time. The transition dipole moment surfaces are strongly dependent on the level of theory at which they are calculated, as is also shown here by benchmark calculations at high symmetry geometries. Therefore, ab initio calculations of the collision-induced absorption spectra cannot become quantitatively predictive unless more accurate transition dipole surfaces can be computed. This is left as an open question for method development in electronic structure theory. The calculated potential energy and transition dipole moment surfaces are employed in quantum dynamical calculations of collision-induced absorption spectra reported in Paper II [T. Karman et al., J. Chem. Phys. 147, 084307 (2017)]. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations are useful in obtaining thermodynamic and kinetic properties of bio-molecules, but they are limited by the time scale barrier. That is, we may not obtain properties' efficiently because we need to run microseconds or longer simulations using femtosecond time steps. To overcome this time scale barrier, we can use the weighted ensemble (WE) method, a powerful enhanced sampling method that efficiently samples thermodynamic and kinetic properties. However, the WE method requires an appropriate partitioning of phase space into discrete macrostates, which can be problematic when we have a high-dimensional collective space or when little is known a priori about the molecular system. Hence, we developed a new WE-based method, called the "Concurrent Adaptive Sampling (CAS) algorithm," to tackle these issues. The CAS algorithm is not constrained to use only one or two collective variables, unlike most reaction coordinate-dependent methods. Instead, it can use a large number of collective variables and adaptive macrostates to enhance the sampling in the high-dimensional space. This is especially useful for systems in which we do not know what the right reaction coordinates are, in which case we can use many collective variables to sample conformations and pathways. In addition, a clustering technique based on the committor function is used to accelerate sampling the slowest process in the molecular system. In this paper, we introduce the new method and show results from two-dimensional models and bio-molecules, specifically penta-alanine and a triazine trimer. [ABSTRACT FROM AUTHOR]

Modeling fluid cycloalkanes with molecular dynamics simulations has proven to be a very challenging task partly because of lacking a reliable force field based on quantum chemistry calculations. In this paper, we construct an ab initio force field for fluid cyclopropane using the second-order Møller-Plesset perturbation theory. We consider 15 conformers of the cyclopropane dimer for the orientation sampling. Single-point energies at important geometries are calibrated by the coupled cluster with single, double, and perturbative triple excitation method. Dunning's correlation consistent basis sets (up to aug-cc-pVTZ) are used in extrapolating the interaction energies at the complete basis set limit. The force field parameters in a 9-site Lennard-Jones model are regressed by the calculated interaction energies without using empirical data. With this ab initio force field, we perform molecular dynamics simulations of fluid cyclopropane and calculate both the structural and dynamical properties. We compare the simulation results with those using an empirical force field and obtain a quantitative agreement for the detailed atom-wise radial distribution functions. The experimentally observed gross radial distribution function (extracted from the neutron scattering measurements) is well reproduced in our simulation. Moreover, the calculated self-diffusion coefficients and shear viscosities are in good agreement with the experimental data over a wide range of thermodynamic conditions. To the best of our knowledge, this is the first ab initio force field which is capable of competing with empirical force fields for simulating fluid cyclopropane. [ABSTRACT FROM AUTHOR]

In this paper, broadband dielectric spectroscopy was applied to investigate molecular dynamics of 1,2,3,4,6-penta-O-(trimethylsilyl)-D-glucopyranose (S-GLU) at ambient and elevated pressures. Our studies showed that apart from the structural relaxation, one well resolved asymmetric secondary process (initially labeled as β) is observed in the spectra measured at p = 0.1 MPa. Analysis with the use of the coupling model and criterion proposed by Ngai and Capaccioli indicated that the β-process in S-GLU is probably a Johari-Goldstein relaxation of intermolecular origin. Further high pressure experiments demonstrated that there are in fact two secondary processes contributing to the β-relaxation. Therefore, one can postulate that the coupling model is a necessary, but not sufficient criterion to identify the true nature of the given secondary relaxation process. The role of pressure experiments in better understanding of the molecular origin of local mobility seems to be much more important. Interestingly, our research also revealed that the structural relaxation in S-GLU is very sensitive to compression. It was reflected in an extremely high pressure coefficient of the glass transition temperature (dTg/dp = 412 K/GPa). According to the literature data, such a high value of dTg/dp has not been obtained so far for any H-bonded, van der Waals, or polymeric glassformers. [ABSTRACT FROM AUTHOR]

The calculation of caloric properties such as heat capacity, Joule-Thomson coefficients, and the speed of sound by classical force-field-based molecular simulation methodology has received scant attention in the literature, particularly for systems composed of complex molecules whose force fields (FFs) are characterized by a combination of intramolecular and intermolecular terms. The calculation of a thermodynamic property for a system whose molecules are described by such a FF involves the calculation of the residual property prior to its addition to the corresponding ideal-gas property, the latter of which is separately calculated, either using thermochemical compilations or nowadays accurate quantum mechanical calculations. Although the simulation of a volumetric residual property proceeds by simply replacing the intermolecular FF in the rigid molecule case by the total (intramolecular plus intermolecular) FF, this is not the case for a caloric property. We describe the correct methodology required to perform such calculations and illustrate it in this paper for the case of the internal energy and the enthalpy and their corresponding molar heat capacities. We provide numerical results for cP, one of the most important caloric properties. We also consider approximations to the correct calculation procedure previously used in the literature and illustrate their consequences for the examples of the relatively simple molecule 2-propanol, CH3CH(OH)CH3, and for the more complex molecule monoethanolamine, HO(CH2)2NH2, an important fluid used in carbon capture. [ABSTRACT FROM AUTHOR]

The multiscale coarse-graining (MS-CG) method is a method for constructing a coarse-grained (CG) model of a system using data obtained from molecular dynamics simulations of the corresponding atomically detailed model. The formal statistical mechanical derivation of the method shows that the potential energy function extracted from an MS-CG calculation is a variational approximation for the true potential of mean force of the CG sites, one that becomes exact in the limit that a complete basis set is used in the variational calculation if enough data are obtained from the atomistic simulations. Most applications of the MS-CG method have employed a representation for the nonbonded part of the CG potential that is a sum of all possible pair interactions. This approach, despite being quite successful for some CG models, is inadequate for some others. Here we propose a systematic method for including three body terms as well as two body terms in the nonbonded part of the CG potential energy. The current method is more general than a previous version presented in a recent paper of this series [L. Larini, L. Lu, and G. A. Voth, J. Chem. Phys. 132, 164107 (2010)], in the sense that it does not make any restrictive choices for the functional form of the three body potential. We use hierarchical multiresolution functions that are similar to wavelets to develop very flexible basis function expansions with both two and three body basis functions. The variational problem is solved by a numerical technique that is capable of automatically selecting an appropriate subset of basis functions from a large initial set. We apply the method to two very different coarse-grained models: a solvent free model of a two component solution made of identical Lennard-Jones particles and a one site model of SPC/E water where a site is placed at the center of mass of each water molecule. These calculations show that the inclusion of three body terms in the nonbonded CG potential can lead to significant improvement in the accuracy of CG potentials and hence of CG simulations. [ABSTRACT FROM AUTHOR]

We present a short overview of pressure calculation in molecular dynamics or Monte Carlo simulations. The emphasis is given to polarizable models in order to resolve the controversy caused by the paper of M. J. Louwerse and E. J. Baerends [Chem. Phys. Lett. 421, 138 (2006)] about pressure calculation in systems with periodic boundaries. We systematically derive expressions for the pressure and show that despite the lack of explicit pairwise additivity, the pressure formula for polarizable models is identical with that of nonpolarizable ones. However, a strict condition for using this formula is that the induced dipole should be in perfect mechanical equilibrium prior to pressure calculation. The perfect convergence of induced dipoles ensures conservation of energy as well. We demonstrate using more cumbersome but exact methods that the derived expressions for the polarizable model of water provide correct numerical results. We also show that the inaccuracy caused by imperfect convergence of the induced dipoles correlates with the inaccuracy of the calculated pressure. [ABSTRACT FROM AUTHOR]

In the quest towards coarse-grained potentials and new water models, we present an extension of the force matching technique to parameterize an all-atom force field for rigid water. The methodology presented here allows to improve the matching procedure by first optimizing the weighting exponents present in the objective function. A new gauge for unambiguously evaluating the quality of the fit has been introduced; it is based on the root mean square difference of the distributions of target properties between reference data and fitted potentials. Four rigid water models have been parameterized; the matching procedure has been used to assess the role of the ghost atom in TIP4P-like models and of electrostatic damping. In the former case, burying the negative charge inside the molecule allows to fit better the torques. In the latter, since short-range interactions are damped, a better fit of the forces is obtained. Overall, the best performing model is the one with a ghost atom and with electrostatic damping. The approach shown in this paper is of general validity and could be applied to any matching algorithm and to any level of coarse graining, also for non-rigid molecules. [ABSTRACT FROM AUTHOR]

Photo-induced electron detachment spectroscopy of anionic boron clusters, B4- and B5-, is theoretically investigated by performing electronic structure calculations and nuclear dynamics simulations. While the electronic potential energy surfaces (X1Ag, ã3B2u, b3B1u, Ã1B2u, c3B2g, and B1B2g of neutral B4 and X2B2, Ã2A1, B2B2, C2A1, D2B1, and E2A1 of neutral B5) and their coupling surfaces are constructed in this paper, the details of the nuclear dynamics on these electronic states are presented in Paper II. Electronic structure calculations are carried out at the complete active space self-consistent field - multi-reference configuration interaction level of theory employing the correlation consistent polarized valance triple zeta basis set. Using the calculated electronic structure data suitable vibronic Hamiltonians are constructed utilizing a diabatic electronic basis and displacement coordinates of the normal vibrational modes. The theoretical results are discussed in relation to those recorded in recent experiments. [ABSTRACT FROM AUTHOR]

We propose a rigorous method for removing rigid-body motions from a given molecular dynamics trajectory of a flexible macromolecule. The method becomes exact in the limit of an infinitesimally small sampling step for the input trajectory. In a recent paper [G. Kneller, J. Chem. Phys. 128, 194101 (2008)], one of us showed that virtual internal atomic displacements for small time increments can be derived from Gauss' principle of least constraint, which leads to a rotational superposition problem for the atomic coordinates in two consecutive time frames of the input trajectory. Here, we demonstrate that the accumulation of these displacements in a molecular-fixed frame, which evolves in time according to the virtual rigid-body motions, leads to the desired trajectory for internal motions. The atomic coordinates in the input and output trajectory are related by a roto-translation, which guarantees that the internal energy of the molecule is left invariant. We present a convenient implementation of our method, in which the accumulation of the internal displacements is performed implicitly. Two numerical examples illustrate the difference to the classical approach for removing macromolecular rigid-body motions, which consists of aligning its configurations in the input trajectory with a fixed reference structure. [ABSTRACT FROM AUTHOR]

In this paper, the nanostructural organization and subpicosecond intermolecular dynamics in the mixtures of CS2 and the room temperature ionic liquid (IL) 1-pentyl-3-methylimidazolium bis{(trifluoromethane)sulfonyl}amide ([C5mim][NTf2]) were studied as a function of concentration using molecular dynamics (MD) simulations and optical heterodyne-detected Raman-induced Kerr effect spectroscopy. At low CS2 concentrations (<10 mol.% CS2/IL), the MD simulations indicate that the CS2 molecules are localized in the nonpolar domains. In contrast, at higher concentrations (≥10 mol.% CS2/IL), the MD simulations show aggregation of the CS2 molecules. The optical Kerr effect (OKE) spectra of the mixtures are interpreted in terms of an additivity model with the components arising from the subpicosecond dynamics of CS2 and the IL. Comparison of the CS2-component with the OKE spectra of CS2 in alkane solvents is consistent with CS2 mainly being localized in the nonpolar domains, even at high CS2 concentrations, and the local CS2 concentration being higher than the bulk CS2 concentration. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations were conducted to investigate the dynamic properties of melts of nonconcatenated ring polymers and compared to melts of linear polymers. The longest rings were composed of N = 1600 monomers per chain which corresponds to roughly 57 entanglement lengths for comparable linear polymers. The ring melts were found to diffuse faster than their linear counterparts, with both architectures approximately obeying a D ∼ N-2.4 scaling law for large N. The mean-square displacement of the center-of-mass of the rings follows a sub-diffusive behavior for times and distances beyond the ring extension g2>, neither compatible with the Rouse nor the reptation model. The rings relax stress much faster than linear polymers, and the zero-shear viscosity was found to vary as η0 ∼ N1.4 ± 0.2 which is much weaker than the N3.4 behavior of linear chains, not matching any commonly known model for polymer dynamics when compared to the observed mean-square displacements. These findings are discussed in view of the conformational properties of the rings presented in the preceding paper [J. D. Halverson, W. Lee, G. S. Grest, A. Y. Grosberg, and K. Kremer, J. Chem. Phys. 134, 204904 (2011)]. [ABSTRACT FROM AUTHOR]

On the mesoscale, the molecular motion in a microporous material can be represented as a sequence of hops between different pore locations and from one pore to the other. On the same scale, the memory effects in the motion of a tagged particle are embedded in the displacement autocorrelation function (DACF), the discrete counterpart of the velocity autocorrelation function (VACF). In this paper, a mesoscopic hopping model, based on a lattice-gas automata dynamics, is presented for the coarse-grained modeling of the DACF in a microporous material under conditions of thermodynamic equilibrium. In our model, that we will refer to as central cell model, the motion of one tagged particle is mimicked through probabilistic hops from one location to the other in a small lattice of cells where all the other particles are indistinguishable; the cells closest to the one containing the tagged particle are simulated explicitly in the canonical ensemble, whereas the border cells are treated as mean-field cells in the grand-canonical ensemble. In the present paper, numerical simulation of the central cell model are shown to provide the same results as a traditional lattice-gas simulation. Along with this a mean-field theory of self-diffusion which incorporates time correlations is discussed. [ABSTRACT FROM AUTHOR]

Trotter's theorem forms the theoretical basis of most modern molecular dynamics. In essence this theorem states that a time displacement operator (a Lie operator) constructed by exponentiating a sum of noncommuting operators can be approximated by a product of single operators provided the time interval is 'very small.' In theory 'very small' implies infinitesimally small (at which point the approximate product becomes exact), while in practical analysis a finite time interval is divided into several small subintervals or steps. It follows, therefore, that the larger the number of steps the better the approximation to the exact time displacement operator. The question therefore arises: How many steps are sufficient? For bounded operators, standard theorems are available to provide the answer. In this paper we show that a very simple combinatorial formula can be derived which allows the computation of the global differences (as a function of the number of steps) between the Taylor coefficients of the exact time displacement operator and an approximate one constructed by using a finite number of steps. The formula holds for both bounded and nonbounded operators and shows, quantitatively, what is qualitatively expected-that the error decreases with increasing number of steps. Furthermore, the formula applies irrespective of the complexity of the system, boundary conditions, or the thermodynamic ensemble employed for averaging the initial conditions. The analysis yields explicit expressions for the Taylor coefficients which are then used to compute the errors. In the case of the algorithmically based practical numerical simulations in which fixed, albeit small, steps are repeatedly applied, the rise in the number of steps does not reduce the size of the steps but increases the total time of interest. The combinatorial formula shows that, here, the errors diverge. Furthermore, this work can be used to supplement other efforts such as the use of shadow Hamiltonians where the truncation of the series expansion of the latter will produce errors in the higher order propagator moments. [ABSTRACT FROM AUTHOR]

In a recent paper, we have developed an efficient implementation of the ring polymer molecular dynamics (RPMD) method for calculating bimolecular chemical reaction rates in the gas phase, and illustrated it with applications to some benchmark atom-diatom reactions. In this paper, we show that the same methodology can readily be used to treat more complex polyatomic reactions in their full dimensionality, such as the hydrogen abstraction reaction from methane, H + CH4→ H2 + CH3. The present calculations were carried out using a modified and recalibrated version of the Jordan-Gilbert potential energy surface. The thermal rate coefficients obtained between 200 and 2000 K are presented and compared with previous results for the same potential energy surface. Throughout the temperature range that is available for comparison, the RPMD approximation gives better agreement with accurate quantum mechanical (multiconfigurational time-dependent Hartree) calculations than do either the centroid density version of quantum transition state theory (QTST) or the quantum instanton (QI) model. The RPMD rate coefficients are within a factor of 2 of the exact quantum mechanical rate coefficients at temperatures in the deep tunneling regime. These results indicate that our previous assessment of the accuracy of the RPMD approximation for atom-diatom reactions remains valid for more complex polyatomic reactions. They also suggest that the sensitivity of the QTST and QI rate coefficients to the choice of the transition state dividing surface becomes more of an issue as the dimensionality of the reaction increases. [ABSTRACT FROM AUTHOR]

In the present paper, we report a molecular dynamics simulation of two-dimensional Lennard-Jones system with a simple square start configuration. Mean square displacement was computed showing interesting dependence on high pressure conditions in short time scale, corresponding to an abrupt restructurization. This paper is the first to report the qualitative and quantitative details of this phenomenon. [ABSTRACT FROM AUTHOR]

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 [ABSTRACT FROM AUTHOR]

A bistable, dipolar stator-rotor molecular system-candidate for molecular electronics is investigated. We demonstrate that it is possible to control the intramolecular torsional states and dynamics in this system by applying an appropriate additional electric field (instead of biasing one), achieving fine tuning and modulation of the relevant properties. The electric field effects on the quantities responsible for torsional dynamics (potential energy surface, potential barrier height, quantum and classical transition probabilities, correlation time, HOMO-LUMO gap) are studied from first principles. Our results indicate that it is possible to artificially stabilize the metastable conformational state of the studied molecule. The importance of this is evident, as the current-voltage characteristics of the metastable state are clearly distinguishable from the current-voltage characteristics of the two stable states. We report for the first time exact calculations related to the possibilities to control the thermally induced stochastic switching, and reduce the noise in a practical application. Thus, we believe that the molecule studied in this paper could operate as a field-switchable molecular device under real conditions. [ABSTRACT FROM AUTHOR]

Attempts to simulate crystal growth of ice from liquid water and to provide a consistent microscopic description of this process have been challenging tasks. In this paper we have adapted our previously developed molecular dynamics simulation methodology to enable the investigation of steady-state directional crystal growth/melting of ice. Specifically, we examine ice/water systems of the (001), (110), and (111) faces of ice Ic and the (0001), [formula], and [formula] faces of ice Ih, where the TIP4P, TIP4P-Ew, and SPC/E water models have been utilized. The influence of different growth/melting conditions (temperature gradients and growth velocities) is investigated. Profile functions of properties of interest across the interface are obtained from nonequilibrium steady-state simulations and provide consistent descriptions of ice/water interfaces. The widths of the various crystallographic faces are found to increase in the apparent order Ic111, Ih0001 < [formula] < [formula] < Ic001 < Ic110. The observed growth rates were in agreement with experimental values and the possible dependence on the various faces is explored. The melting temperatures obtained with the present methodology for the different models are in good agreement with estimates from other work. [ABSTRACT FROM AUTHOR]

Over the last number of years several simulation methods have been introduced to study rare events such as nucleation. In this paper we examine the crystal nucleation rate of hard spheres using three such numerical techniques: molecular dynamics, forward flux sampling, and a Bennett-Chandler-type theory where the nucleation barrier is determined using umbrella sampling simulations. The resulting nucleation rates are compared with the experimental rates of Harland and van Megen [Phys. Rev. E 55, 3054 (1997)], Sinn et al. [Prog. Colloid Polym. Sci. 118, 266 (2001)], Schätzel and Ackerson [Phys. Rev. E 48, 3766 (1993)], and the predicted rates for monodisperse and 5% polydisperse hard spheres of Auer and Frenkel [Nature 409, 1020 (2001)]. When the rates are examined in units of the long-time diffusion coefficient, we find agreement between all the theoretically predicted nucleation rates, however, the experimental results display a markedly different behavior for low supersaturation. Additionally, we examined the precritical nuclei arising in the molecular dynamics, forward flux sampling, and umbrella sampling simulations. The structure of the nuclei appears independent of the simulation method, and in all cases, the nuclei contains on average significantly more face-centered-cubic ordered particles than hexagonal-close-packed ordered particles. [ABSTRACT FROM AUTHOR]

In this paper, we examine how and when quantum evolution can be approximated in terms of (generalized) classical dynamics in calculations of correlation functions, with a focus on the symmetrized time correlation function introduced by Schofield. To that end, this function is expressed as a path integral in complex time and written in terms of sum and difference path variables. Taylor series expansion of the path integral's exponent to first and second order in the difference variables leads to two original developments. The first order expansion is used to obtain a simple, path integral based, derivation of the so-called Schofield's quantum correction factor. The second order result is employed to show how quantum mechanical delocalization manifests itself in the approximation of the correlation function and hinders, even in the semiclassical limit, the interpretation of the propagators in terms of sets of guiding classical trajectories dressed with appropriate weights. [ABSTRACT FROM AUTHOR]

This is the first in a series of two papers investigating the effect of electron-phonon coupling in two-dimensional Fourier transformed electronic spectroscopy. We present a series of one- and two-dimensional nonlinear spectroscopic techniques for studying a dye molecule in solution. Ultrafast laser pulse excitation of an electronic transition coupled to vibrational modes induces a propagating vibrational wave packet that manifests itself in oscillating signal intensities and line shapes. For the two-dimensional electronic spectra we can attribute the observed modulations to periodic enhancement and decrement of the relative amplitudes of rephasing and nonrephasing contributions to the total response. Different metrics of the two-dimensional signals are shown to relate to the frequency-frequency correlation function which provides the connection between experimentally accessible observations and the underlying microscopic molecular dynamics. A detailed theory of the time-dependent two-dimensional spectral line shapes is presented in the accompanying paper [T. Mančal et al., J. Chem. Phys. 132, 184515 (2010)]. [ABSTRACT FROM AUTHOR]

Capillary waves have been observed in systems ranging from the surfaces of ordinary fluids to interfaces in biological membranes and have been one of the most studied areas in the physics of fluids. Recent advances in fluorescence microscopy and imaging enabled quantitative measurements of thermally driven capillary waves in lipid monolayers and bilayers, which resulted in accurate measurements of the line tension in monolayer domains. Even though there has been a considerable amount of work on the statics and dynamics of capillary waves in three dimensions, to the best of our knowledge, there is no detailed theoretical analysis for two-dimensional droplet morphologies. In this paper, we derive the dynamic correlation function for two-dimensional fluid droplets using linear response theory and verify our results using a novel particle-based simulation technique for binary mixtures. [ABSTRACT FROM AUTHOR]

Gas-phase atomic anions lack bound electronic excited states, yet in solution many of these anions exhibit intense absorption bands due to the presence of excited states, referred to as charge-transfer-to-solvent (CTTS) states that are bound only by the presence of the solvent. CTTS spectra thus serve as delicate probes of solute-solvent interactions, but the fact that they are created by the interactions of a solute with many solvent molecules makes them a challenge to describe theoretically. In this paper, we use mixed quantum/classical molecular dynamics with the two-electron Fourier-grid (2EFG) electronic structure method presented in the previous paper [W. J. Glover, R. E. Larsen, and B. J. Schwartz, J. Chem. Phys. 132, 144101 (2010)] to simulate the CTTS states of a sodium anion in liquid tetrahydrofuran, Na-/THF. Since our 2EFG method is based on configuration interaction with single and double excitations in a grid basis, it allows for an exact treatment of the two valence electrons of the sodium anion. To simulate Na-/THF, we first develop a new electron–THF pseudopotential, and we verify the accuracy of this potential by reproducing the experimental absorption spectrum of an excess electron in liquid THF with near quantitative accuracy. We also are able to reproduce the CTTS spectrum of Na-/THF and find that the CTTS states of Na- exhibit a Rydberg-like progression due to the pre-existing long-range solvent polarization around the anion. We also find that the CTTS states are highly mixed with the disjoint electronic states supported by naturally occurring solvent cavities that exist in liquid THF. This mixing explains why the solvated electrons that are ejected following CTTS excitation appear with their equilibrium absorption spectrum. The mixing of the CTTS and solvent-cavity states also explains why the recombination of the electron and its geminate Na0 partner occurs on slower time scales when photoexciting in the blue compared to in the red of the CTTS band: blue excitation accesses CTTS states whose charge densities lies further from the position of the anion, whereas red excitation accesses CTTS states that lie primarily within the anion’s first solvation shell. Finally, we see that the radial character of the CTTS states near the Na+ core matches that of Na0, explaining why the spectrum of this species appears instantly after photoexciting Na-. [ABSTRACT FROM AUTHOR]

In this paper we investigate the possibility to detect the hypothesized liquid-liquid critical point of water in supercooled aqueous solutions of salts. Molecular dynamics computer simulations are conducted on bulk TIP4P water and on an aqueous solution of sodium chloride in TIP4P water, with concentration c=0.67 mol/kg. The liquid-liquid critical point is found both in the bulk and in the solution. Its position in the thermodynamic plane shifts to higher temperature and lower pressure for the solution. Comparison with available experimental data allowed us to produce the phase diagrams of both bulk water and the aqueous solution as measurable in experiments. Given the position of the liquid-liquid critical point in the solution as obtained from our simulations, the experimental determination of the hypothesized liquid-liquid critical point of water in aqueous solutions of salts appears possible. [ABSTRACT FROM AUTHOR]

In this paper, based on the integrated tempering sampling we introduce a selective integrated tempering sampling (SITS) method for the efficient conformation sampling and thermodynamics calculations for a subsystem in a large one, such as biomolecules solvated in aqueous solutions. By introducing a potential surface scaled with temperature, the sampling over the configuration space of interest (e.g., the solvated biomolecule) is selectively enhanced but the rest of the system (e.g., the solvent) stays largely unperturbed. The applications of this method to biomolecular systems allow highly efficient sampling over both energy and configuration spaces of interest. Comparing to the popular and powerful replica exchange molecular dynamics (REMD), the method presented in this paper is significantly more efficient in yielding relevant thermodynamics quantities (such as the potential of mean force for biomolecular conformational changes in aqueous solutions). It is more important that SITS but not REMD yielded results that are consistent with the traditional umbrella sampling free energy calculations when explicit solvent model is used since SITS avoids the sampling of the irrelevant phase space (such as the boiling water at high temperatures). [ABSTRACT FROM AUTHOR]

The molecular dynamics expression of heat flux, originally derived by Irving and Kirkwood [J. Chem. Phys. 18, 817 (1950)] for pairwise potentials, is generalized in this paper for systems with many-body potentials. The original formula consists of a kinetic part and a potential part, and the latter term is found in the present study to be expressible as a summation of contributions from all the many-body potentials defined in the system. The energy transfer among a set of sites for which a many-body potential is defined is discussed and evaluated by the rate of increase in the kinetic energy of each site due to the potential, and its accumulation over all the potentials in the system is shown to make up the potential part of the generalized expression. A molecular dynamics simulation for liquid n-octane was performed to demonstrate the applicability of the new expression obtained in this study to measure the heat flux and to elucidate the contributions of inter- and intramolecular potentials to heat conduction. [ABSTRACT FROM AUTHOR]

The adsorption, desorption, and diffusion kinetics of N2 on thick (up to ∼9 μm) porous films of amorphous solid water (ASW) films were studied using molecular beam techniques and temperature programmed desorption. Porous ASW films were grown on Pt(111) at low temperature (<30 K) from a collimated H2O beam at glancing incident angles. In thin films (<1 μm), the desorption kinetics are well described by a model that assumes rapid and uniform N2 distribution throughout the film. In thicker films (>1 μm), N2 adsorption at 27 K results in a nonuniform distribution, where most of N2 is trapped in the outer region of the film. Redistribution of N2 can be induced by thermal annealing. The apparent activation energy for this process is ∼7 kJ/mol, which is approximately half of the desorption activation energy at the corresponding coverage. Preadsorption of Kr preferentially adsorbs onto the highest energy binding sites, thereby preventing N2 from trapping in the outer region of the film which facilitates N2 transport deeper into the porous film. Despite the onset of limited diffusion, the adsorption kinetics are efficient, precursor mediated, and independent of film thickness. An adsorption mechanism is proposed, in which a high-coverage N2 front propagates into a pore by the rapid transport of physisorbed second layer N2 species on top of the first surface bound layer. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations of NaCl fluid are used to understand the behavior of ionic fluid to screen the field generated by charges on the ionic crystal surfaces in absence of any external electric field. The NaCl fluid in the strongly coupled regime (corresponding to the melt) in contact with the charged octopolar (111) NaCl surface shows that the spatial correlations decay in an oscillatory manner, with a screening length λQ given by the envelope of the damped oscillations. By contrast to the Debye-Huckel theory, in the strongly coupled regime, λQ increases with increasing coupling strength (also seen in bulk ionic simulations). The NaCl fluid confined between neutral (100) NaCl surfaces also shows weak oscillatory charge decay near the surface. Similar oscillatory exponential decay was seen when the NaCl fluid was confined between two analytically smooth neutral walls. The origin of these oscillations was due to the difference in ion sizes. NaCl fluid confined between neutral octopolar (110) and dipolar (110) surface show stronger density oscillations than (100) surface but comparatively very weak charge oscillations. This paper shows that the strength of the charges on the crystal surfaces is enough to induce a characteristic spatial distribution of charges in the contacting fluid and the extent of distribution depends on the type of surface. [ABSTRACT FROM AUTHOR]

Markovian state models (MSMs) are a convenient and efficient means to compactly describe the kinetics of a molecular system as well as a formalism for using many short simulations to predict long time scale behavior. Building a MSM consists of grouping the conformations into states and estimating the transition probabilities between these states. In a previous paper, we described an efficient method for calculating the uncertainty due to finite sampling in the mean first passage time between two states. In this paper, we extend the uncertainty analysis to derive similar closed-form solutions for the distributions of the eigenvalues and eigenvectors of the transition matrix, quantities that have numerous applications when using the model. We demonstrate the accuracy of the distributions on a six-state model of the terminally blocked alanine peptide. We also show how to significantly reduce the total number of simulations necessary to build a model with a given precision using these uncertainty estimates for the blocked alanine system and for a 2454-state MSM for the dynamics of the villin headpiece. [ABSTRACT FROM AUTHOR]

A numerical solution for the quantum dynamics of the spin boson problem is obtained using the semiclassical initial value series representation approach to the quantum dynamics. The zeroth order term of the series is computed using the new forward-backward representation for correlation functions presented in the preceding adjacent paper. This leads to a rapid convergence of the Monte Carlo sampling, as compared to previous attempts. The zeroth order results are already quite accurate. The first order term of the series is small, demonstrating the rapid convergence of the semiclassical initial value representation series. This is the first time that the first order term in the semiclassical initial value representation series has been converged for systems with the order of 50 degrees of freedom. [ABSTRACT FROM AUTHOR]

In this paper, we present a path integral hybrid Monte Carlo (PIHMC) method for rotating molecules in quantum fluids. This is an extension of our PIHMC for correlated Bose fluids [S. Miura and J. Tanaka, J. Chem. Phys. 120, 2160 (2004)] to handle the molecular rotation quantum mechanically. A novel technique referred to be an effective potential of quantum rotation is introduced to incorporate the rotational degree of freedom in the path integral molecular dynamics or hybrid Monte Carlo algorithm. For a permutation move to satisfy Bose statistics, we devise a multilevel Metropolis method combined with a configurational-bias technique for efficiently sampling the permutation and the associated atomic coordinates. Then, we have applied the PIHMC to a helium-4 cluster doped with a carbonyl sulfide molecule. The effects of the quantum rotation on the solvation structure and energetics were examined. Translational and rotational fluctuations of the dopant in the superfluid cluster were also analyzed. [ABSTRACT FROM AUTHOR]

Multiblock polymers in aqueous solution, where one or several blocks are hydrophobic, exhibit a rich variety of phases and states of aggregation. In this paper, we investigate a pentablock system ABCBA, where the B block is always hydrophilic and the A and C blocks have varying degrees of hydrophobicity depending on external conditions. We report coarse-grained molecular-dynamics simulations where the solvent is included explicitly and monomers interact via a 6-9 Lennard Jones potential function. The hydrophobic interaction is modeled by tuning the parameter controlling the strength of the interaction between the hydrophobic monomers and the solvent. We investigate the structure and morphology of the micelles for two concrete situations representing changes in temperature and the pH level. The simulated system is directly relevant to a recently synthesized pentablock system consisting of a triblock Pluronic® with an added pH-sensitive end group [B. C. Anderson et al., Macromolecules 36, 1670 (2003)]. [ABSTRACT FROM AUTHOR]

We present an adaptive time stepping scheme based on the extrapolative method of Barth and Schlick [LN, J. Chem. Phys. 109, 1633 (1998)] to numerically integrate the Langevin equation with a molecular-dynamics potential. This approach allows us to use (on average) a time step for the strong nonbonded force integration corresponding to half the period of the fastest bond oscillation, without compromising the slow degrees of freedom in the problem. We show with simple examples how the dynamic step size stabilizes integration operators, and discuss some of the limitations of such stability. The method introduced uses a slightly more accurate inner integrator than LN to accommodate the larger steps. The adaptive time step approach reproduces temporal features of the bovine pancreatic trypsin inhibitor (BPTI) test system (similar to the one used in the original introduction of LN) compared to short-time integrators, but with energies that are shifted with respect to both LN, and traditional stochastic versions of Verlet. Although the introduction of longer steps has the effect of systematically heating the bonded components of the potential, the temporal fluctuations of the slow degrees of freedom are reproduced accurately. The purpose of this paper is to display a mechanism by which the resonance traditionally associated with using time steps corresponding to half the period of oscillations in molecular dynamics can be avoided. This has theoretical utility in terms of designing numerical integration schemes—the key point is that by factoring a propagator so that time steps are not constant one can recover stability with an overall (average) time step at a resonance frequency. There are, of course, limitations to this approach associated with the complicated, nonlinear nature of the molecular-dynamics (MD) potential (i.e., it is not as straightforward as the linear test problem we use to motivate the method). While the basic notion remains in the full Newtonian problem, it is easier to see the effects when damping is considered to be physical—that is, we do not view our method as a perturbation of Newtonian dynamics, we associate the damping with the environment, for example, a water bath (with γ≈90 ps-1) [Zagrovic and Pande, J. Comp. Chem. 24, 1432 (2003)]. All stochastic approaches to MD are stabilized by large physical damping, but here, we are really using it only to show that the resonance frequency can be obtained. Another simplifying assumption used in this paper is “heavy” hydrogen (we take the hydrogen mass to be 10 amu)—the view here is that we are interested primarily in the slowest degrees of freedom, and this approach has effects similar to bond freezing and united atom treatments of hydrogen. So from the point of view of biomolecular applications, the method described here is best suited to studies in which water is not explicit (so that damping in the problem can really be viewed as environmental interaction), and the interest is in slow dynamics where the effects of hydrogen are neglectable. There are a number of parameters in the LN method and the one derived here, and we cannot in a short paper address all adjustments, so our primary goal as a first pass is to show that stability can be recovered for a set of numerically forced (and hence artificial) bond oscillations, and compare stability to fixed-step methods. [ABSTRACT FROM AUTHOR]