A new parametrization of the widely used Stillinger-Weber potential is proposed for silicon, allowing for an improved modelling of defects and plasticity-related properties. The performance of the new potential is compared to the original version, as well as to another parametrization (Vink et al 2001 J. Non-Cryst. Solids, 282 248), in the case of several situations: point defects and dislocation core stability, threshold displacement energies, bulk shear, generalized stacking fault energy surfaces, fracture, melting temperature, amorphous structure, and crystalline phase stability. A significant improvement is obtained in the case of dislocation cores, bulk behaviour under high shear stress, the amorphous structure, and computation of threshold displacement energies, while most of the features of the original version (elastic constants, point defects) are retained. However, despite a slight improvement, a complex process like fracture remains difficult to model.
Download full-text PDF |
Source |
---|---|
http://dx.doi.org/10.1088/0953-8984/25/5/055801 | DOI Listing |
J Phys Chem Lett
August 2023
Pacific Northwest National Laboratory, Richland, Washington 99352, United States.
Reactive force fields (RFFs) are an expedient approach to sample chemical reaction paths in complex systems, relative to density functional theory. However, there is continued need to improve efficiencies, specifically in systems that have slow transverse degrees of freedom, as in highly viscous and superconcentrated solutions. Here, we present an RFF that is differentiated from current models (e.
View Article and Find Full Text PDFJ Chem Phys
June 2022
Energy Technologies Area, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA.
In this paper, we consider the problem of quantifying parametric uncertainty in classical empirical interatomic potentials (IPs) using both Bayesian (Markov Chain Monte Carlo) and frequentist (profile likelihood) methods. We interface these tools with the Open Knowledgebase of Interatomic Models and study three models based on the Lennard-Jones, Morse, and Stillinger-Weber potentials. We confirm that IPs are typically sloppy, i.
View Article and Find Full Text PDFJ Chem Theory Comput
May 2020
Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany.
Current machine learning (ML) models aimed at learning force fields are plagued by their high computational cost at every integration time step. We describe a number of practical and computationally efficient strategies to parametrize traditional force fields for molecular liquids from ML: the particle decomposition ansatz to two- and three-body force fields, the use of kernel-based ML models that incorporate physical symmetries, the incorporation of switching functions close to the cutoff, and the use of covariant meshing to boost the training set size. Results are presented for model molecular liquids: pairwise Lennard-Jones, three-body Stillinger-Weber, and bottom-up coarse-graining of water.
View Article and Find Full Text PDFPhys Chem Chem Phys
January 2018
CNR.IPCF, via Moruzzi 1, 56124 Pisa, Italy.
Melting and sintering of silicon nanoparticles are investigated by means of classical molecular dynamics simulations to disclose the dependence of modelling on the system type, the simulation procedure and interaction potential. The capability of our parametrization of a reactive force field ReaxFF to describe such processes is assessed through a comparison with formally simpler Stillinger-Weber and Tersoff potentials, which are frequently used for simulating silicon-based materials. A substantial dependence of both the predicted melting point and its variation as a function of the nanoparticle size on the simulation model is also highlighted.
View Article and Find Full Text PDFJ Chem Phys
June 2016
Chemistry and Chemical Biology, University of California, Merced, 5200 North Lake Road, Merced, California 95343, USA.
Three empirical interatomic force fields are parametrized using structural, elastic, and phonon dispersion data for bulk CdSe and their predictions are then compared for the structures and phonons of CdSe quantum dots having average diameters of ˜2.8 and ˜5.2 nm (˜410 and ˜2630 atoms, respectively).
View Article and Find Full Text PDFEnter search terms and have AI summaries delivered each week - change queries or unsubscribe any time!