scholarly journals Employing Hybrid Lennard-Jones and Axilrod-Teller Potentials to Parametrize Force Fields for the Simulation of Materials’ Properties

Materials ◽  
2021 ◽  
Vol 14 (21) ◽  
pp. 6352
Author(s):  
Danilo de Camargo Branco ◽  
Gary J. Cheng

The development of novel materials has challenges besides their synthesis. Materials such as novel MXenes are difficult to probe experimentally due to their reduced size and low stability under ambient conditions. Quantum mechanics and molecular dynamics simulations have been valuable options for material properties determination. However, computational materials scientists may still have difficulty finding specific force field models for their simulations. Force fields are usually hard to parametrize, and their parameters’ determination is computationally expensive. We show the Lennard-Jones (2-body interactions) combined with the Axilrod-Teller (3-body interactions) parametrization process’ applicability for metals and new classes of materials (MXenes). Because this parametrization process is simple and computationally inexpensive, it allows users to predict materials’ behaviors under close-to-ambient conditions in molecular dynamics, independent of pre-existing potential files. Using the process described in this work, we have made the Ti2C parameters set available for the first time in a peer-reviewed work.

Author(s):  
Toshihiro Kaneko ◽  
Kenji Yasuoka ◽  
Ayori Mitsutake ◽  
Xiao Cheng Zeng

Multicanonical molecular dynamics simulations are applied, for the first time, to study the liquid-solid and solid-solid transitions in Lennard-Jones (LJ) clusters. The transition temperatures are estimated based on the peak position in the heat capacity versus temperature curve. For LJ31, LJ58 and LJ98, our results on the solid-solid transition temperature are in good agreement with previous ones. For LJ309, the predicted liquid-solid transition temperature is also in agreement with previous result.


2020 ◽  
Vol 9 (1) ◽  
pp. 700-715 ◽  
Author(s):  
Wei Jian ◽  
David Hui ◽  
Denvid Lau

AbstractRecent advances in biomedicine largely rely on the development in nanoengineering. As the access to unique properties in biomaterials is not readily available from traditional techniques, the nanoengineering becomes an effective approach for research and development, by which the performance as well as the functionalities of biomaterials has been greatly improved and enriched. This review focuses on the main materials used in biomedicine, including metallic materials, polymers, and nanocomposites, as well as the major applications of nanoengineering in developing biomedical treatments and techniques. Research that provides an in-depth understanding of material properties and efficient enhancement of material performance using molecular dynamics simulations from the nanoengineering perspective are discussed. The advanced techniques which facilitate nanoengineering in biomedical applications are also presented to inspire further improvement in the future. Furthermore, the potential challenges of nanoengineering in biomedicine are evaluated by summarizing concerned issues and possible solutions.


1993 ◽  
Vol 317 ◽  
Author(s):  
N.A. Marks ◽  
P. Guan ◽  
D.R. Mckenzie ◽  
B.A. PailThorpe

ABSTRACTMolecular dynamics simulations of nickel and carbon have been used to study the phenomena due to ion impact. The nickel and carbon interactions were described using the Lennard-Jones and Stillinger-Weber potentials respectively. The phenomena occurring after the impact of 100 e V to 1 keV ions were studied in the nickel simulations, which were both two and three-dimensional. Supersonic focussed collision sequences (or focusons) were observed, and associated with these focusons were unexpected sonic bow waves, which were a major energy loss mechanism for the focuson. A number of 2D carbon films were grown and the stress in the films as a function of incident ion energy was Measured. With increasing energy the stress changed from tensile to compressive and reached a maximum around 50 eV, in agreement with experiment.


Author(s):  
Francesca Ferrari ◽  
Maicol Bissaro ◽  
Simone Fabbian ◽  
Jessica de Almeida Roger ◽  
Stefano Mammi ◽  
...  

<p>In this manuscript, for the first time, we presented a fragment library and we validated its performance by comparison with a well-established technique for fragment screening as solution NMR. We were able to screen 400 different fragments producing a total of 1200 independent fragment-protein recognition pathways. As far as we know, this represents the largest screening based on Molecular dynamics ever reported. Our simulations successfully detected the true binders in the library in a prospective study, showing a notable agreement with a state-of-art screening we performed by NMR on the same dataset.</p>


2021 ◽  
Author(s):  
Xiangyun Lei ◽  
Andrew Medford

Abstract Molecular dynamics simulations are an invaluable tool in numerous scientific fields. However, the ubiquitous classical force fields cannot describe reactive systems, and quantum molecular dynamics are too computationally demanding to treat large systems or long timescales. Reactive force fields based on physics or machine learning can be used to bridge the gap in time and length scales, but these force fields require substantial effort to construct and are highly specific to a given chemical composition and application. A significant limitation of machine learning models is the use of element-specific features, leading to models that scale poorly with the number of elements. This work introduces the Gaussian multipole (GMP) featurization scheme that utilizes physically-relevant multipole expansions of the electron density around atoms to yield feature vectors that interpolate between element types and have a fixed dimension regardless of the number of elements present. We combine GMP with neural networks to directly compare it to the widely used Behler-Parinello symmetry functions for the MD17 dataset, revealing that it exhibits improved accuracy and computational efficiency. Further, we demonstrate that GMP-based models can achieve chemical accuracy for the QM9 dataset, and their accuracy remains reasonable even when extrapolating to new elements. Finally, we test GMP-based models for the Open Catalysis Project (OCP) dataset, revealing comparable performance to graph convolutional deep learning models. The results indicate that this featurization scheme fills a critical gap in the construction of efficient and transferable machine-learned force fields.


2021 ◽  
Author(s):  
Martin P. Lautenschlaeger ◽  
Hans Hasse

It was shown recently that using the two-gradient method, thermal, caloric, and transport properties of fluids under quasi-equilibrium conditions can be determined simultaneously from nonequilibrium molecular dynamics simulations. It is shown here that the influence of shear stresses on these properties can also be studied using the same method. The studied fluid is described by the Lennard-Jones truncated and shifted potential with the cut-off radius r*c = 2.5σ. For a given temperature T and density ρ, the influence of the shear rate on the following fluid properties is determined: pressure p, internal energy u, enthalpy h, isobaric heat capacity cp, thermal expansion coefficient αp, shear viscosity η, and self-diffusion coefficient D. Data for 27 state points in the range of T ∈ [0.7, 8.0] and ρ ∈ [0.3, 1.0] are reported for five different shear rates (γ ̇ ∈ [0.1,1.0]). Correlations for all properties are provided and compared with literature data. An influence of the shear stress on the fluid properties was found only for states with low temperature and high density. The shear-rate dependence is caused by changes in the local structure of the fluid which were also investigated in the present work. A criterion for identifying the regions in which a given shear stress has an influence on the fluid properties was developed. It is based on information on the local structure of the fluid. For the self-diffusivity, shear-induced anisotropic effects were observed and are discussed.


2019 ◽  
Author(s):  
Thiago José Pinheiro dos Santos ◽  
Charlles Abreu ◽  
Bruno Horta ◽  
Frederico W. Tavares

Mass transport coefficients play an important role in process design and in compositional grading of oil reservoirs. As experimental measurements of these properties can be costly and hazardous, Molecular Dynamics simulations emerge as an alternative approach. In this work, we used Molecular Dynamics to calculate the self-diffusion coefficients of methane/n-hexane mixtures at different conditions, in both liquid and supercritical phases. We evaluated how the finite box size and the choice of the force field affect the calculated properties at high pressures. Results show a strong dependency between self-diffusion and the simulation box size. The Yeh-Hummer analytical correction [J. Phys. Chem. B, 108, 15873 (2004)] can attenuate this effect, but sometimes makes the results depart from experimental data due to issues concerning the force fields. We have also found that different all-atom and united-atom models can produce biased results due to caging effects and to different dihedral configurations of the n-alkane.


Sign in / Sign up

Export Citation Format

Share Document