Predictive Modelling » Publications
https://warwick.ac.uk/fac/sci/wcpm/research/publications/
The latest from Predictive Modelling » Publicationsen-GB(C) 2022 University of WarwickWed, 21 Jul 2021 07:36:22 GMThttp://blogs.law.harvard.edu/tech/rssSiteBuilder2, University of Warwick, http://go.warwick.ac.uk/sitebuilderabartokpartayanandashahckumarasinghecortnerdcrevillenganandhsadeghijrkermodeklawohnldesousamaldegundemicardimmousavinstocksnzabarasonatpbrommerpgraziosipgrigorevsabolfathisullivanvvargiamidiswleeUntaggedModelling defects in Ni–Al with EAM and DFT calculations
https://warwick.ac.uk/fac/sci/wcpm/research/publications/?newsItem=094d434554536416015461780a31533b
<p>F. Bianchini, J.R. Kermode, and A. De Vita, <a href="http://dx.doi.org/10.1088/0965-0393/24/4/045012">Modelling defects in Ni–Al with EAM and DFT calculations</a>, Modell. Simul. Mater. Sci. Eng. <b>24</b>, 045012 (2016).</p>
<p>We present detailed comparisons between the results of embedded atom model (EAM) and density functional theory (DFT) calculations on defected Ni alloy systems. We find that the EAM interatomic potentials reproduce low-temperature structural properties in both the γ and ${{\gamma}^{\prime}}$ phases, and yield accurate atomic forces in bulk-like configurations even at temperatures as high as  ~1200 K. However, they fail to describe more complex chemical bonding, in configurations including defects such as vacancies or dislocations, for which we observe significant deviations between the EAM and DFT forces, suggesting that derived properties such as (free) energy barriers to vacancy migration and dislocation glide may also be inaccurate. Testing against full DFT calculations further reveals that these deviations have a local character, and are typically severe only up to the first or second neighbours of the defect. This suggests that a QM/MM approach can be used to accurately reproduce QM observables, fully exploiting the EAM potential efficiency in the MM zone. This approach could be easily extended to ternary systems for which developing a reliable and fully transferable EAM parameterisation would be extremely challenging e.g. Ni alloy model systems with a W or Re-containing QM zone.</p>jrkermodeFri, 29 Apr 2016 10:02:17 GMT094d434554536416015461780a31533bA universal preconditioner for simulating condensed phase materials
https://warwick.ac.uk/fac/sci/wcpm/research/publications/?newsItem=094d434554536587015453aff2640214
<p>D. Packwood, J. Kermode, L. Mones, N. Bernstein, J. Woolley, N. Gould, C. Ortner, and G. Csányi, <a href="http://dx.doi.org/10.1063/1.4947024"></a><a style="font-family: Calibri, 'Bitstream Vera Sans', 'Trebuchet MS', Trebuchet, sans-serif; font-size: 14px; font-weight: normal; line-height: 21px; word-spacing: normal;" href="http://dx.doi.org/10.1063/1.4947024">A universal preconditioner for simulating condensed phase materials</a>, J. Chem. Phys. <b>144</b>, 164109 (2016)</p>
<p>We introduce a universal sparse preconditioner that accelerates geometry optimisation and saddle point search tasks that are common in the atomic scale simulation of materials. Our preconditioner is based on the neighbourhood structure and we demonstrate the gain in computational efficiency in a wide range of materials that include metals, insulators, and molecular solids. The simple structure of the preconditioner means that the gains can be realised in practice not only when using expensive electronic structuremodels but also for fast empirical potentials. Even for relatively small systems of a few hundred atoms, we observe speedups of a factor of two or more, and the gain grows with system size. An open source Python implementation within the Atomic Simulation Environment is available, offering interfaces to a wide range of atomistic codes.</p>jrkermodeTue, 26 Apr 2016 17:48:40 GMT094d434554536587015453aff2640214Development of an exchange–correlation functional with uncertainty quantification capabilities for density functional theory
http://www.sciencedirect.com/science/article/pii/S0021999116000425
<p>Manuel Aldegunde, James R. Kermode, and Nicholas Zabaras, <em>J. Comput. Phys.</em> <strong>311</strong>, 173-195, doi:<a href="http://dx.doi.org/10.1016/j.jcp.2016.01.034">10.1016/j.jcp.2016.01.034</a></p>
<p>This paper presents the development of a new exchange–correlation functional from the point of view of machine learning. Using atomization energies of solids and small molecules, we train a linear model for the exchange enhancement factor using a Bayesian approach which allows for the quantification of uncertainties in the predictions. A relevance vector machine is used to automatically select the most relevant terms of the model. We then test this model on atomization energies and also on bulk properties. The average model provides a mean absolute error of only 0.116 eV for the test points of the G2/97 set but a larger 0.314 eV for the test solids. In terms of bulk properties, the prediction for transition metals and monovalent semiconductors has a very low test error. However, as expected, predictions for types of materials not represented in the training set such as ionic solids show much larger errors.</p>nzabarasmaldegundejrkermodeTue, 16 Feb 2016 09:32:19 GMT094d434552cfc2200152e96c5f6a05beLow Speed Crack Propagation via Kink Formation and Advance on the Silicon (110) Cleavage Plane
http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.115.135501
<p>James R. Kermode, Anna Gleizer, Guy Kovel, Lars Pastewka, Gábor Csányi, Dov Sherman, and Alessandro De Vita, <em>Phys. Rev. Lett.</em> <strong>115</strong>, 135501, doi:<a href="http://dx.doi.org/10.1103/PhysRevLett.115.135501">10.1103/PhysRevLett.115.135501</a><br />
</p>
<p>We present density functional theory based atomistic calculations predicting that slow fracturing of silicon is possible at any chosen crack propagation speed under suitable temperature and load conditions. We also present experiments demonstrating fracture propagation on the Si(110) cleavage plane in the <span class="aps-inline-formula"><span class="MathJax_Preview"></span><span class="MathJax"><nobr style="box-sizing: inherit; transition: none; border: 0px; padding: 0px; margin: 0px; max-width: none; max-height: none; min-width: 0px; min-height: 0px; vertical-align: 0px;"><span class="math"><span class="mrow"><span class="mrow"><span class="mo">∼</span><span class="mn">100</span><span class="mtext"> </span><span class="mtext"> </span><span class="mi">m</span><span class="mo">/</span><span class="mi">s</span></span></span></span></nobr></span></span> speed range, consistent with our predictions. These results suggest that many other brittle crystals could be broken arbitrarily slowly in controlled experiments.</p>jrkermodeThu, 24 Sep 2015 07:34:44 GMT094d43f54fcc38e2014ffe469d080a38Classical interaction potentials for diverse materials from ab initio data: a review of potfit
http://dx.doi.org/10.1088/0965-0393/23/7/074002
<p class="mb-0"><span class="nowrap">Peter Brommer</span>, <span class="nowrap">Alexander Kiselev</span>, <span class="nowrap">Daniel Schopf</span>, <span class="nowrap">Philipp Beck</span>, <span class="nowrap">Johannes Roth</span> and <span class="nowrap">Hans-Rainer Trebin, </span><em>Modelling Simul. Mater. Sci. Eng.</em> <b>23 </b>074002 (2015). doi:<a href="http://dx.doi.org/10.1088/0965-0393/23/7/074002" target="_blank">10.1088/0965-0393/23/7/074002</a></p>
<p class="mb-0">Force matching is an established technique to generate effective potentials for molecular dynamics simulations from first-principles data. This method has been implemented in the open source code potfit . Here, we present a review of the method and describe the main features of the code. Particular emphasis is placed on the features added since the initial release: interactions represented by analytical functions, differential evolution as optimization method, and a greatly extended set of interaction models. Beyond the initially present pair and embedded-atom method potentials, potfit can now also optimize angular dependent potentials, charge and dipolar interactions, and electron-temperature-dependent potentials. We demonstrate the functionality of these interaction models using three example systems: phonons in type I clathrates, fracture of α-alumina, and laser-irradiated silicon</p>
<p class="mb-0"></p>pbrommerFri, 18 Sep 2015 12:42:34 GMT094d43f54fcc3779014fe07a46c114d1Effect of oxygen and nitrogen functionalization on the physical and electronic structure of graphene
http://dx.doi.org/10.1007/s12274-015-0768-0
<p>Alexander J. Marsden, Peter Brommer, James J. Mudd, M. Adam Dyson, Robert Cook, María Asensio, Jose Avila, Ana Levy, Jeremy Sloan, David Quigley, Gavin R. Bell, and Neil R. Wilson, <em>Nano Research</em> (2015). doi: <a href="http:/dx.doi.org/10.1007/s12274-015-0768-0" target="_blank">10.1007/s12274-015-0768-0</a></p>
<p>Covalent functionalization of graphene offers opportunities for tailoring its properties and is an unavoidable consequence of some graphene synthesis techniques. However, the changes induced by the functionalization are not well understood. By using atomic sources to control the extent of the oxygen and nitrogen functionalization, we studied the evolution in the structure and properties at the atomic scale. Atomic oxygen reversibly introduces epoxide groups whilst, under similar conditions, atomic nitrogen irreversibly creates diverse functionalities including substitutional, pyridinic, and pyrrolic nitrogen. Atomic oxygen leaves the Fermi energy at the Dirac point (i.e., undoped), whilst atomic nitrogen results in a net n-doping; however, the experimental results are consistent with the dominant electronic effect for both being a transition from delocalized to localized states, and hence the loss of the signature electronic structure of graphene.</p>pbrommerMon, 20 Jul 2015 10:03:03 GMT094d43454e772a27014eaaeaac7b3da2Diffusion of point defects in crystalline silicon using the kinetic activation-relaxation technique method
http://link.aps.org/doi/10.1103/PhysRevB.91.224106
<p class="pub-info">Mickaël Trochet, Laurent Karim Béland, Jean-François Joly, Peter Brommer, and Normand Mousseau, <em>Phys. Rev. B</em> <b>91</b>, 224106 (2015). doi:<a href="http://dx.doi.org/10.1103/PhysRevB.91.224106">10.1103/PhysRevB.91.224106</a></p>
<p class="pub-info">We study point-defect diffusion in crystalline silicon using the kinetic activation-relaxation technique (k-ART), an off-lattice kinetic Monte Carlo method with on-the-fly catalog building capabilities based on the activation-relaxation technique (ART nouveau), coupled to the standard Stillinger-Weber potential. We focus more particularly on the evolution of crystalline cells with one to four vacancies and one to four interstitials in order to provide a detailed picture of both the atomistic diffusion mechanisms and overall kinetics. We show formation energies, activation barriers for the ground state of all eight systems, and migration barriers for those systems that diffuse. Additionally, we characterize diffusion paths and special configurations such as dumbbell complex, di-interstitial (IV-pair+2I) superdiffuser, tetrahedral vacancy complex, and more. This study points to an unsuspected dynamical richness even for this apparently simple system that can only be uncovered by exhaustive and systematic approaches such as the kinetic activation-relaxation technique.</p>pbrommerThu, 18 Jun 2015 10:51:57 GMT094d43454dfc7ab8014e064bf2bd6834A framework for machine-learning-augmented multiscale atomistic simulations on parallel supercomputers
http://dx.doi.org/10.1002/qua.24952
<p>Caccin, M., Li, Z., Kermode, J. R. & De Vita, A. A framework for machine-learning-augmented multiscale atomistic simulations on parallel supercomputers. <i>Int. J. Quantum Chem.</i> (2015). doi:10.1002/qua.24952</p>
<p>Recent advances in quantum mechanical (QM)-based molecular dynamics (MD) simulations have used machine-learning (ML) to predict, rather than recalculate, QM-accurate forces in atomic configurations sufficiently similar to previously encountered ones. Here, we discuss how ML approaches can be deployed within large-scale QM/MM materials simulations on massively parallel supercomputers, making QM zones of ≳1000 atoms routinely attainable. We argue that the ML approach allows computational effort to be concentrated on the most chemically active subregions of the QM zone, significantly improving the overall efficiency of the simulation. We thus propose a novel method to partition large QM regions into multiple subregions, which can be computed in parallel to achieve optimal scaling. Then we review a recently proposed QM/ML MD scheme (Z. Li, J.R. Kermode, A. De Vita Phys. Rev. Lett., 2015, 114, 096405), discussing how this could be efficiently combined with QM-zone partitioning.</p>jrkermodeTue, 09 Jun 2015 09:07:28 GMT094d43454daffde4014dd7930efc2976Molecular Dynamics with On-the-Fly Machine Learning of Quantum-Mechanical Forces
http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.114.096405
<p>Zhenwei Li, James R. Kermode and Alessandro De Vita, Phys. Rev. Lett. <strong>114</strong>, 096405 (2015)<br />
doi:10.1103/PhysRevLett.114.096405</p>
<p>We present a molecular dynamics scheme which combines first-principles and machine-learning (ML) techniques in a single information-efficient approach. Forces on atoms are either predicted by Bayesian inference or, if necessary, computed by on-the-fly quantum-mechanical (QM) calculations and added to a growing ML database, whose completeness is, thus, never required. As a result, the scheme is accurate and general, while progressively fewer QM calls are needed when a new chemical process is encountered for the second and subsequent times, as demonstrated by tests on crystalline and molten silicon.</p>jrkermodeFri, 06 Mar 2015 16:08:21 GMT094d43454bef8b8b014befd83ea90a9aFollowing atomistic kinetics on experimental timescales with the kinetic Activation–Relaxation Technique
http://dx.doi.org/10.1016/j.commatsci.2014.11.047
<p>N. Mousseau, L.K. Béland, P. Brommer, F. El-Mellouhi, J.-F. Joly, G.K. N'Tsouaglo, O. Restrepo, M. Trochet, <em>Comp. Mater. Sci. </em><strong>100 B</strong>, 111–123 (2015), doi:<a href="http://dx.doi.org/10.1016/j.commatsci.2014.11.047">10.1016/j.commatsci.2014.11.047</a></p>
<p>The properties of materials, even at the atomic level, evolve on macroscopic time scales. Following this evolution through simulation has been a challenge for many years. For lattice-based activated diffusion, kinetic Monte Carlo has turned out to be an almost perfect solution. Various accelerated molecular dynamical schemes, for their part, have allowed the study on long time scale of relatively simple systems. There is still a need, however, for methods able to handle complex materials such as alloys and disordered systems. Here, we review the kinetic Activation–Relaxation Technique (k-ART), one of a handful of off-lattice kinetic Monte Carlo methods, with on-the-fly cataloging, that have been proposed in the last few years.</p>pbrommerMon, 02 Mar 2015 15:34:12 GMT094d43454bda134f014bdb1f8cb70e33