Development of a ReaxFF Force Field for Cu/S/C/H and Reactive MD Simulations of Methyl Thiolate Decomposition on Cu (100)

We present a hierarchical Bayesian framework for the selection of force fields in molecular dynamics (MD) simulations. The framework associates the variability of the optimal parameters of the MD potentials under different environmental conditions with the corresponding variability in experimental data. The high computational cost associated with the hierarchical Bayesian framework is reduced by orders of magnitude through a parallelized Transitional Markov Chain Monte Carlo method combined with the Laplace Asymptotic Approximation. The suitability of the hierarchical approach is demonstrated by performing MD simulations with prescribed parameters to obtain data for transport coefficients under different conditions, which are then used to infer and evaluate the parameters of the MD model. We demonstrate the selection of MD models based on experimental data and verify that the hierarchical model can accurately quantify the uncertainty across experiments; improve the posterior probability density function estimation of the parameters, thus, improve predictions on future experiments; identify the most plausible force field to describe the underlying structure of a given dataset. The framework and associated software are applicable to a wide range of nanoscale simulations associated with experimental data with a hierarchical structure.

Download Full-text

Using C-DFT to Develop an e-ReaxFF Force Field for Acetophenone Radical Anion

The Journal of Chemical Physics ◽

10.1063/5.0064705 ◽

2021 ◽

Author(s):

Katheryn Anne Penrod ◽

Maximiliano Aldo Burgess ◽

Dooman Akbarian ◽

Ismaila Dabo ◽

William Henry Hunter Woodward ◽

...

Keyword(s):

Force Field ◽

Radical Anion ◽

Reaxff Force Field

Download Full-text

Simulating Absorption Spectra of Flavonoids in Aqueous Solution: A Polarizable QM/MM Study

Molecules ◽

10.3390/molecules25245853 ◽

2020 ◽

Vol 25 (24) ◽

pp. 5853

Author(s):

Sulejman Skoko ◽

Matteo Ambrosetti ◽

Tommaso Giovannini ◽

Chiara Cappelli

Keyword(s):

Molecular Dynamics ◽

Aqueous Solution ◽

Hydrogen Bonding ◽

Force Field ◽

Absorption Spectra ◽

Computational Study ◽

Md Simulations ◽

Quantum Mechanical ◽

Hydrogen Bonding Interactions

We present a detailed computational study of the UV/Vis spectra of four relevant flavonoids in aqueous solution, namely luteolin, kaempferol, quercetin, and myricetin. The absorption spectra are simulated by exploiting a fully polarizable quantum mechanical (QM)/molecular mechanics (MM) model, based on the fluctuating charge (FQ) force field. Such a model is coupled with configurational sampling obtained by performing classical molecular dynamics (MD) simulations. The calculated QM/FQ spectra are compared with the experiments. We show that an accurate reproduction of the UV/Vis spectra of the selected flavonoids can be obtained by appropriately taking into account the role of configurational sampling, polarization, and hydrogen bonding interactions.

Download Full-text

A force field for MD simulations on rhenium organometallic compounds developed from enthalpy of sublimation and X-ray diffraction measurements

The Journal of Chemical Thermodynamics ◽

10.1016/j.jct.2019.01.016 ◽

2019 ◽

Vol 133 ◽

pp. 60-69 ◽

Cited By ~ 1

Author(s):

Carlos E.S. Bernardes ◽

Mariana T. Donato ◽

M. Fátima M. Piedade ◽

Hermínio P. Diogo ◽

José N. Canongia Lopes ◽

...

Keyword(s):

Force Field ◽

Organometallic Compounds ◽

Md Simulations ◽

Enthalpy Of Sublimation ◽

X Ray Diffraction ◽

X Ray

Download Full-text

Force field development and simulations of senior dialkyl sulfoxides

Physical Chemistry Chemical Physics ◽

10.1039/c5cp08006a ◽

2016 ◽

Vol 18 (15) ◽

pp. 10507-10515 ◽

Cited By ~ 8

Author(s):

Vitaly V. Chaban

Keyword(s):

Molecular Dynamics ◽

Force Field ◽

Ab Initio ◽

Md Simulations ◽

Field Development ◽

Force Field Development ◽

Structure And Dynamics ◽

Ethyl Methyl ◽

Methyl Sulfoxide ◽

Diethyl Sulfoxide

Thermodynamics, structure, and dynamics of diethyl sulfoxide (DESO) and ethyl methyl sulfoxide (EMSO) were investigated using ab initio calculations and non-polarizable potential based molecular dynamics (MD) simulations.

Download Full-text