Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Molecular Dynamics Simulations of Melting Iron Nanoparticles with/without Defects Using a Reaxff Reactive Force Field.

ABSTRACT: Molecular dynamics simulations are performed to study thermal properties of bulk iron material and Fe nanoparticles (FNP) by using a ReaxFF reactive force field. Thermodynamic and energy properties such as radial distribution function, Lindemann index and potential energy plots are adopted to study the melting behaviors of FNPs from 300?K to 2500?K. A step-heating method is introduced to obtain equilibrium melting points. Our results show ReaxFF force field is able to detect size effect in FNP melting no matter in energy or structure evolution aspect. Extra storage energy of FNPs caused by defects (0%-10%) is firstly studied in this paper: defects will not affect the melting point of FNPs directly but increase the system energy especially when temperature reaches the melting points.

SUBMITTER: Sun J

PROVIDER: S-EPMC7042232 | biostudies-literature | 2020 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Molecular Dynamics Simulations of Melting Iron Nanoparticles with/without Defects Using a Reaxff Reactive Force Field.

Sun Junlei J Liu Pingan P Wang Mengjun M Liu Junpeng J

Scientific reports 20200225 1

Molecular dynamics simulations are performed to study thermal properties of bulk iron material and Fe nanoparticles (FNP) by using a ReaxFF reactive force field. Thermodynamic and energy properties such as radial distribution function, Lindemann index and potential energy plots are adopted to study the melting behaviors of FNPs from 300 K to 2500 K. A step-heating method is introduced to obtain equilibrium melting points. Our results show ReaxFF force field is able to detect size effect in FNP m ...[more]

PMID: 32099061

Similar Datasets

ReaxFF/AMBER-A Framework for Hybrid Reactive/Nonreactive Force Field Molecular Dynamics Simulations.

Project description:Combined quantum mechanical/molecular mechanical (QM/MM) models using semiempirical and ab initio methods have been extensively reported on over the past few decades. These methods have been shown to be capable of providing unique insights into a range of problems, but they are still limited to relatively short time scales, especially QM/MM models using ab initio methods. An intermediate approach between a QM based model and classical mechanics could help fill this time-scale gap and facilitate the study of a range of interesting problems. Reactive force fields represent the intermediate approach explored in this paper. A widely used reactive model is ReaxFF, which has largely been applied to materials science problems and is generally used as a stand-alone (i.e., the full system is modeled using ReaxFF). We report a hybrid ReaxFF/AMBER molecular dynamics (MD) tool, which introduces ReaxFF capabilities to capture bond breaking and formation within the AMBER MD software package. This tool enables us to study local reactive events in large systems at a fraction of the computational costs of QM/MM models. We describe the implementation of ReaxFF/AMBER, validate this implementation using a benzene molecule solvated in water, and compare its performance against a range of similar approaches. To illustrate the predictive capabilities of ReaxFF/AMBER, we carried out a Claisen rearrangement study in aqueous solution. In a first for ReaxFF, we were able to use AMBER's potential of mean force (PMF) capabilities to perform a PMF study on this organic reaction. The ability to capture local reaction events in large systems using combined ReaxFF/AMBER opens up a range of problems that can be tackled using this model to address both chemical and biological processes.

| S-EPMC8145783 | biostudies-literature

Development of ReaxFF Reactive Force Field for Aqueous Iron-Sulfur Clusters with Applications to Stability and Reactivity in Water.

Project description:Iron-sulfur clusters serve unique roles in biochemistry, geochemistry, and renewable energy technologies. However, a full theoretical understanding of their structures and properties is still lacking. To facilitate large-scale reactive molecular dynamics simulations of iron-sulfur clusters in aqueous environments, a ReaxFF reactive force field is developed, based on an extensive set of quantum chemical calculations. This force field compares favorably with the reference calculations on gas-phase species and significantly improves on a previous ReaxFF parametrization. We employ the new potential to study the stability and reactivity of iron-sulfur clusters in explicit water with constant-temperature reactive molecular dynamics. The aqueous species exhibit a dynamic, temperature-dependent behavior, in good agreement with previous much more costly ab initio simulations.

| S-EPMC8028049 | biostudies-literature

Development of a ReaxFF reactive force field for glycine and application to solvent effect and tautomerization.

Project description:Tautomerization of amino acids between the neutral form (NF) and the zwitterionic form (ZW) in water has been extensively studied, often using glycine as a model to understand this fundamental process. In spite of many advanced studies, the tautomerization reaction remains poorly understood because of the intrinsic complexities of the system, including multiple accessible reaction pathways, charge transfer, and variations of solvation structure. To establish an accurate model that can be used for molecular dynamics simulations, a ReaxFF reactive force field has been developed for glycine. A training set for the ReaxFF hydrocarbon potential was augmented with several glycine conformers and glycine-water complexes. The force field parameters were optimized to reproduce the quantum mechanically derived energies of the species in the training set. The optimized potential could accurately describe the properties of gas-phase glycine. It was applied to investigate the effect of solvation on the conformational distribution of glycine. Molecular dynamics simulations indicated significant differences in the dominant conformers in the gas phase and in water. This suggests that the tautomerization of glycine occurs through a conformational isomerization followed by the proton transfer event. The direct reaction mechanism of the NF ? ZW proton transfer reaction in water, as well as mechanisms mediated by one or two water molecules, were investigated using molecular dynamics simulations. The results suggest that the proton transfer reaction is most likely mediated by a single water molecule. The ReaxFF potential developed in this work provides an accurate description of proton transfer in glycine and thus provides a useful methodology for simulating proton transfer reactions in organic molecules in the aqueous environment.

| S-EPMC3042430 | biostudies-literature

Development and validation of a ReaxFF reactive force field for Cu cation/water interactions and copper metal/metal oxide/metal hydroxide condensed phases.

Project description:To enable large-scale reactive dynamic simulations of copper oxide/water and copper ion/water interactions we have extended the ReaxFF reactive force field framework to Cu/O/H interactions. To this end, we employed a multistage force field development strategy, where the initial training set (containing metal/metal oxide/metal hydroxide condensed phase data and [Cu(H(2)O)(n)](2+) cluster structures and energies) is augmented by single-point quantum mechanices (QM) energies from [Cu(H(2)O)(n)](2+) clusters abstracted from a ReaxFF molecular dynamics simulation. This provides a convenient strategy to both enrich the training set and to validate the final force field. To further validate the force field description we performed molecular dynamics simulations on Cu(2+)/water systems. We found good agreement between our results and earlier experimental and QM-based molecular dynamics work for the average Cu/water coordination, Jahn-Teller distortion, and inversion in [Cu(H(2)O)(6)](2+) clusters and first- and second-shell O-Cu-O angular distributions, indicating that this force field gives a satisfactory description of the Cu-cation/water interactions. We believe that this force field provides a computationally convenient method for studying the solution and surface chemistry of metal cations and metal oxides and, as such, has applications for studying protein/metal cation complexes, pH-dependent crystal growth/dissolution, and surface catalysis.

| S-EPMC4224311 | biostudies-literature

Force field bias in protein folding simulations.

Project description:Long timescale (>1 micros) molecular dynamics simulations of protein folding offer a powerful tool for understanding the atomic-scale interactions that determine a protein's folding pathway and stabilize its native state. Unfortunately, when the simulated protein fails to fold, it is often unclear whether the failure is due to a deficiency in the underlying force fields or simply a lack of sufficient simulation time. We examine one such case, the human Pin1 WW domain, using the recently developed deactivated morphing method to calculate free energy differences between misfolded and folded states. We find that the force field we used favors the misfolded states, explaining the failure of the folding simulations. Possible further applications of deactivated morphing and implications for force field development are discussed.

| S-EPMC2711430 | biostudies-literature

Parmbsc1: a refined force field for DNA simulations.

Project description:We present parmbsc1, a force field for DNA atomistic simulation, which has been parameterized from high-level quantum mechanical data and tested for nearly 100 systems (representing a total simulation time of ? 140 ?s) covering most of DNA structural space. Parmbsc1 provides high-quality results in diverse systems. Parameters and trajectories are available at http://mmb.irbbarcelona.org/ParmBSC1/.

| S-EPMC4700514 | biostudies-literature

Force field influences in beta-hairpin folding simulations.

Project description:All-atom force fields are now routinely used for more detailed understanding of protein folding mechanisms. However, it has been pointed out that use of all-atom force fields does not guarantee more accurate representations of proteins; in fact, sometimes it even leads to biased structural distributions. Indeed, several issues remain to be solved in force field developments, such as accurate treatment of implicit solvation for efficient conformational sampling and proper treatment of backbone interactions for secondary structure propensities. In this study, we first investigate the quality of several recently improved backbone interaction schemes in AMBER for folding simulations of a beta-hairpin peptide, and further study their influences on the peptide's folding mechanism. Due to the significant number of simulations needed for a thorough analysis of tested force fields, the implicit Poisson-Boltzmann solvent was used in all simulations. The chosen implicit solvent was found to be reasonable for studies of secondary structures based on a set of simulations of both alpha-helical and beta-hairpin peptides with the TIP3P explicit solvent as benchmark. Replica exchange molecular dynamics was also utilized for further efficient conformational sampling. Among the tested AMBER force fields, ff03 and a revised ff99 force field were found to produce structural and thermodynamic data in comparably good agreement with the experiment. However, detailed folding pathways, such as the order of backbone hydrogen bond zipping and the existence of intermediate states, are different between the two force fields, leading to force field-dependent folding mechanisms.

| S-EPMC2242415 | biostudies-literature

ReaxFF Force Field Development and Application for Toluene Adsorption on MnMO<i><sub>x</sub></i> (M = Cu, Fe, Ni) Catalysts.

Project description:In numerous studies, the application of the molecular dynamics scheme based on the reactive force field (ReaxFF) method has been proven effective in modeling the catalytic behavior of metal-organic compounds. Recently, this method has been successfully applied for MxOy (M = Cu, Fe, Mn, Ni) transition-metal oxides. Yet, bimetallic metal oxides of the type MnMOx (M = Cu, Fe, Ni) were also present in the experimental system but could not be modeled since not all of the force field parameters were available at the time. To bridge this gap, the force field for modeling bimetallic metal oxides had to be developed. Here, we establish the needed force field parameter sets (namely, Cu/Mn/O, Fe/Mn/O, and Ni/Mn/O) and apply them to the problem of toluene adsorption on bimetallic oxide catalyst surfaces to verify their validity. Each training set consisted of at least 10 crystal structures containing at least Cu-Mn-O, Fe-Mn-O, or Ni-Mn-O atoms in contact obtained from the available structure databases. The parameter training has been done using the in-home-compiled version of the ReaxFF code. After training the force fields for geometry reproduction, the parameters were refined using the optimization by atom charges, comparing the ReaxFF values to those obtained for the respective structures using periodic crystal density functional theory (DFT) codes. The as-developed force fields were then applied to the process of toluene adsorption/degradation on MnMOx catalysts. Results obtained show agreement with previous experimental expectations, although some remarks are given since the initially presumed crystal structure of bimetallic oxide Mn1-xMxOy crystallites may still have an impact on theoretical predictions. The presented are, to the best of the authors' knowledge, the first applications of the ReaxFF approach to the Mn-(Cu|Fe|Ni)-O-C-H interaction.

| S-EPMC8713286 | biostudies-literature

Mixing ReaxFF parameters for transition metal oxides using force-matching method.

Project description:We present the development of the method for the refitting the ReaxFF parameters for a system consisting of the mixed transition metal oxides. We have tested several methods allowing to calculate the differences between the vectors of the forces acting on atoms obtained from the reference DFT simulation and the parameters-dependent ReaxFF. We conclude that the footrule method yields the best parameters among the investigated options. We then validate the parameters using the system consisting of the hematite supported (TiO2)n clusters. The results indicate the refitted parameters allow to obtain acceptable geometries of the clusters upon MD simulation on the ReaxFF level, and despite the short timescale lead to the stable structures.

| S-EPMC8671285 | biostudies-literature

GEM*: A Molecular Electronic Density-Based Force Field for Molecular Dynamics Simulations.

Project description:GEM*, a force field that combines Coulomb and Exchange terms calculated with Hermite Gaussians with the polarization, bonded, and modified van der Waals terms from AMOEBA is presented. GEM* is tested on an initial water model fitted at the same level as AMOEBA. The integrals required for the evaluation of the intermolecular Coulomb interactions are efficiently evaluated by means of reciprocal space methods. The GEM* water model is tested by comparing energies and forces for a series of water oligomers and MD simulations. Timings for GEM* compared to AMOEBA are presented and discussed.

| S-EPMC5207213 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data