skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Benchmarking Adaptive Steered Molecular Dynamics (ASMD) on CHARMM Force Fields
Abstract The potentials of mean force (PMFs) along the end‐to‐end distance of two different helical peptides have been obtained and benchmarked using the adaptive steered molecular dynamics (ASMD) method. The results depend strongly on the choice of force field driving the underlying all‐atom molecular dynamics, and are reported with respect to the three most popular CHARMM force field versions: c22, c27 and c36. Two small peptides,and 1PEF, serve as the particular case studies. The comparisons between the versions of the CHARMM force fields provides both a qualitative and quantitative look at their performance in forced unfolding simulations in which peptides undergo large changes in structural conformations. We find that ASMD with the underlying c36 force field provides the most robust results for the selected benchmark peptides.  more » « less
Award ID(s):
2102455
PAR ID:
10370711
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
ChemPhysChem
Volume:
23
Issue:
17
ISSN:
1439-4235
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We study the stability and nonlinear local dynamics of spectrally stable periodic wave trains of the Korteweg‐de Vries/Kuramoto‐Sivashinsky equation when subjected to classes of periodic perturbations. It is known that for each, such a‐periodic wave train is asymptotically stable to‐periodic, i.e., subharmonic, perturbations, in the sense that initially nearby data will converge asymptotically to a small Galilean boost of the underlying wave, with exponential rates of decay. However, both the allowable size of initial perturbations and the exponential rates of decay depend onand, in fact, tend to zero as, leading to a lack of uniformity in such subharmonic stability results. Our goal here is to build upon a recent methodology introduced by the authors in the reaction–diffusion setting and achieve a subharmonic stability result, which is uniform in. This work is motivated by the dynamics of such wave trains when subjected to perturbations that are localized (i.e., integrable on the line). 
    more » « less
  2. Abstract Strain engineering is an effective strategy in modulating activity of electrocatalysts, but the effect of strain on electrochemical stability of catalysts remains poorly understood. In this work, we combineab initiothermodynamics and molecular dynamics simulations to examine the role of compressive and tensile strain in the interplay between activity and stability of metal oxides considering RuOand IrOas exemplary catalysts. We reveal that although compressive strain leads to improved activity via the adsorbate‐evolving mechanism of the oxygen evolution reaction, even small strains should substantially destabilize these catalysts promoting dissolution of transition metals. In contrast, our results show that the metal oxides requiring tensile strain to promote their catalytic activity may also benefit from enhanced stability. Importantly, we also find that the detrimental effect of strain on electrochemical stability of atomically flat surfaces could be even greater than that of surface defects. 
    more » « less
  3. Abstract We consider the mapping properties of the integral operator arising in nonlocal slender body theory (SBT) for the model geometry of a straight, periodic filament. It is well known that the classical singular SBT integral operator suffers from high wavenumber instabilities, making it unsuitable for approximating theslender body inverse problem, where the fiber velocity is prescribed and the integral operator must be inverted to find the force density along the fiber. Regularizations of the integral operator must therefore be used instead. Here, we consider two regularization methods: spectral truncation and the‐regularization of Tornberg and Shelley (2004). We compare the mapping properties of these approximations to the underlying partial differential equation (PDE) solution, which for the inverse problem is simply the Stokes Dirichlet problem with data constrained to be constant on cross sections. For the straight‐but‐periodic fiber with constant radius, we explicitly calculate the spectrum of the operator mapping fiber velocity to force for both the PDE and the approximations. We prove that the spectrum of the original SBT operator agrees closely with the PDE operator at low wavenumbers but differs at high frequencies, allowing us to define a truncated approximation with a wavenumber cutoff. For both the truncated and‐regularized approximations, we obtain rigorous‐based convergence to the PDE solution as: A fiber velocity withregularity givesconvergence, while a fiber velocity with at leastregularity yieldsconvergence. Moreover, we determine the dependence of the‐regularized error estimate on the regularization parameter. 
    more » « less
  4. Abstract The mechanisms underlying observed global patterns of partitioning precipitation () to evapotranspiration () and runoff () are controversially debated. We test the hypothesis that asynchrony between climatic water supply and demand is sufficient to explain spatio‐temporal variability of water availability. We developed a simple analytical model forthat is determined by four dimensionless characteristics of intra‐annual water supply and demand asynchrony. The analytical model, populated with gridded climate data, accurately predicted global runoff patterns within 2%–4% of independent estimates from global climate models, with spatial patterns closely correlated to observations (). The supply‐demand asynchrony hypothesis provides a physically based explanation for variability of water availability using easily measurable characteristics of climate. The model revealed widespread responsiveness of water budgets to changes in climate asynchrony in almost every global region. Furthermore, the analytical model using global averages independently reproduced the Budyko curve () providing theoretical foundation for this widely used empirical relationship. 
    more » « less
  5. Abstract The apparent end of the internally generated Martian magnetic field at 3.6–4.1 Ga is a key event in Martian history and has been linked to insufficient core cooling. We investigate the thermal and magnetic evolution of the Martian core and mantle using parameterized models and considered three improvements on previous studies. First, our models account for thermal stratification in the core. Second, the models are constrained by estimates for the present‐day areotherm. Third, we consider core thermal conductivity,, values in the range 5–40 Was suggested by recent experiments on iron alloys at Mars core conditions. The majority of our models indicate that the core of Mars is fully conductive at present with core temperatures greater than 1940 K. All of our models are consistent with the range ofW. Models with an activation volume of 6 (0)require a mantle reference viscosity of Pa s. 
    more » « less