1
|
Zongo K, Sun H, Ouellet-Plamondon C, Béland LK. A unified moment tensor potential for silicon, oxygen, and silica. NPJ COMPUTATIONAL MATERIALS 2024; 10:218. [PMID: 39282246 PMCID: PMC11399103 DOI: 10.1038/s41524-024-01390-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/19/2023] [Accepted: 08/17/2024] [Indexed: 09/18/2024]
Abstract
Si and its oxides have been extensively explored in theoretical research due to their technological importance. Simultaneously describing interatomic interactions within both Si and SiO2 without the use of ab initio methods is considered challenging, given the charge transfers involved. Herein, this challenge is overcome by developing a unified machine learning interatomic potentials describing the Si/SiO2/O system, based on the moment tensor potential (MTP) framework. This MTP is trained using a comprehensive database generated using density functional theory simulations, encompassing diverse crystal structures, point defects, extended defects, and disordered structure. Extensive testing of the MTP is performed, indicating it can describe static and dynamic features of very diverse Si, O, and SiO2 atomic structures with a degree of fidelity approaching that of DFT.
Collapse
Affiliation(s)
- Karim Zongo
- Département de génie de la construction, École de technologie supérieure, Université du Québec, Montréal, QC Canada
| | - Hao Sun
- Department of Mechanical and Materials Engineering, Queen's university, Kingston, ON Canada
| | - Claudiane Ouellet-Plamondon
- Département de génie de la construction, École de technologie supérieure, Université du Québec, Montréal, QC Canada
| | - Laurent Karim Béland
- Department of Mechanical and Materials Engineering, Queen's university, Kingston, ON Canada
| |
Collapse
|
2
|
Luo H, O'Rourke JG, Deng J. Radiogenic heating sustains long-lived volcanism and magnetic dynamos in super-Earths. SCIENCE ADVANCES 2024; 10:eado7603. [PMID: 39270025 PMCID: PMC11397497 DOI: 10.1126/sciadv.ado7603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Accepted: 08/09/2024] [Indexed: 09/15/2024]
Abstract
Radiogenic heat production is fundamental to the energy budget of planets. Roughly half of the heat that Earth loses through its surface today comes from the three long-lived, heat-producing elements (potassium, thorium, and uranium). These three elements have long been believed to be highly lithophile and thus concentrate in the mantle of rocky planets. However, our study shows that they all become siderophile under the pressure and temperature conditions relevant to the core formation of large rocky planets dubbed super-Earths. Mantle convection in super-Earths is then primarily driven by heating from the core rather than by a mix of internal heating and cooling from above as in Earth. Partitioning these sources of radiogenic heat into the core remarkably increases the core-mantle boundary (CMB) temperature and the total heat flow across the CMB in super-Earths. Consequently, super-Earths are likely to host long-lived volcanism and strong magnetic dynamos.
Collapse
Affiliation(s)
- Haiyang Luo
- Department of Geosciences, Princeton University, Princeton, NJ, USA
| | - Joseph G O'Rourke
- School of Earth and Space Exploration, Arizona State University, Tempe, AZ, USA
| | - Jie Deng
- Department of Geosciences, Princeton University, Princeton, NJ, USA
| |
Collapse
|
3
|
Berger F, Rybicki M, Sauer J. Molecular Dynamics with Chemical Accuracy─Alkane Adsorption in Acidic Zeolites. ACS Catal 2023. [DOI: 10.1021/acscatal.2c05493] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Fabian Berger
- Institut für Chemie, Humboldt-Universität zu Berlin, D-10099Berlin, Germany
| | - Marcin Rybicki
- Institut für Chemie, Humboldt-Universität zu Berlin, D-10099Berlin, Germany
| | - Joachim Sauer
- Institut für Chemie, Humboldt-Universität zu Berlin, D-10099Berlin, Germany
- Department of Physical and Macromolecular Chemistry & Charles University Center of Advanced Materials, Charles University, Hlavova 8, 128 43Prague 2, Czech Republic
| |
Collapse
|
4
|
De Wispelaere K, Plessow PN, Studt F. Toward Computing Accurate Free Energies in Heterogeneous Catalysis: a Case Study for Adsorbed Isobutene in H-ZSM-5. ACS PHYSICAL CHEMISTRY AU 2022; 2:399-406. [PMID: 36855690 PMCID: PMC9955322 DOI: 10.1021/acsphyschemau.2c00020] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Herein, we propose a novel computational protocol that enables calculating free energies with improved accuracy by combining the best available techniques for enthalpy and entropy calculation. While the entropy is described by enhanced sampling molecular dynamics techniques, the energy is calculated using ab initio methods. We apply the method to assess the stability of isobutene adsorption intermediates in the zeolite H-SSZ-13, a prototypical problem that is computationally extremely challenging in terms of calculating enthalpy and entropy. We find that at typical operating conditions for zeolite catalysis (400 °C), the physisorbed π-complex, and not the tertiary carbenium ion as often reported, is the most stable intermediate. This method paves the way for sampling-based techniques to calculate the accurate free energies in a broad range of chemistry-related disciplines, thus presenting a big step forward toward predictive modeling.
Collapse
Affiliation(s)
- Kristof De Wispelaere
- Center
for Molecular Modeling, Ghent University, Technologiepark 46, B-9052 Ghent, Belgium,
| | - Philipp N. Plessow
- Institute
of Catalysis Research and Technology, Karlsruhe
Institute of Technology, Hermann-von-Helmholtz-Platz 1, D-76344 Eggenstein-Leopoldshafen, Germany,
| | - Felix Studt
- Institute
of Catalysis Research and Technology, Karlsruhe
Institute of Technology, Hermann-von-Helmholtz-Platz 1, D-76344 Eggenstein-Leopoldshafen, Germany,Institute
for Chemical Technology and Polymer Chemistry, Karlsruhe Institute of Technology, D-76131 Karlsruhe, Germany,
| |
Collapse
|
5
|
Löffelsender S, Schwerdtfeger P, Grimme S, Mewes JM. It's Complicated: On Relativistic Effects and Periodic Trends in the Melting and Boiling Points of the Group 11 Coinage Metals. J Am Chem Soc 2021; 144:485-494. [PMID: 34965098 DOI: 10.1021/jacs.1c10881] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
While the color of metallic gold is a prominent and well-investigated example for the impact of relativistic effects, much less is known regarding the influence on its melting and boiling point (MP/BP). To remedy this situation, this work takes on the challenging task of exploring the phase transitions of the Group 11 coinage metals Cu, Ag, and Au through nonrelativistic (NR) and scalar/spin-orbit relativistic (SR/SOR) Gibbs energy calculations with λ-scaled density-functional theory (λDFT). At the SOR level, the calculations provide BPs in excellent agreement with experimental values (1%), while MPs exhibit more significant deviations (2-10%). Comparing SOR calculations to those conducted in the NR limit reveals some remarkably large and, at the same time, some surprisingly small relativistic shifts. Most notably, the BP of Au increases by about 800 K due to relativity, which is in line with the strong relativistic increase of the cohesive energy, whereas the MP of Au is very similar at the SOR and NR levels, defying the typically robust correlation between MP and cohesive energy. Eventually, an inspection of thermodynamic quantities traces the trend-breaking behavior of Au back to phase-specific effects in liquid Au, which render NR Au more similar to SOR Ag, in line with a half-a-century-old hypothesis of Pyykkö.
Collapse
Affiliation(s)
- Sarah Löffelsender
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University, Auckland Campus, 0632 Auckland, New Zealand
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Jan-Michael Mewes
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| |
Collapse
|
6
|
Accurate acid dissociation constant (pK a) calculation for the sulfachloropyridazine and similar molecules. J Mol Model 2021; 27:233. [PMID: 34324066 DOI: 10.1007/s00894-021-04851-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2021] [Accepted: 07/05/2021] [Indexed: 10/20/2022]
Abstract
Accurate calculation of the acid dissociation constant (pKa) has fundamental importance for the description of molecular systems with pharmacological activities. The search for a more appropriate procedure for its determination is always welcome and has aroused increasing interest from the scientific community. In this sense, this work presents a computational study involving the combination of ten DFT functionals (M062X, M06L, B3LYP, BLYP, PBEPBE, BP86, LC-BLYP, SPBE, CAM-B3LYP, LC-PBEPBE) and HF method, eight basis set functions (6-311G, 6-311 + G, 6-311G(d,p), 6-311 + G(d,p), 6-311+ +G(d,p), 6-311(2d,2p), 6-311+ +G(2d,2p), and aug-cc-pVDZ), and three solvation models (SMD, PCM, and CPCM) for an accurate sulfachloropyridazine (SCR) pKa determination. It was found that the smallest deviation (0.02 unit of pKa) between the current study and experimental result was achieved with the BLYP/6-311 + G(d,p)/PCM combination. Therefore, this combination was extended to calculate the pKa of six SCR similar molecules selected through the eletroshape similarity method. For all these molecules, the difference between the obtained results and experimental data ranged between 0.14 and 0.69 units of pKa. This feature suggests that the obtained combination can determine pKa with experimental precision for complexes that are formed by sulfonamide functional group (SO2NHR). Graphical Abstract A computational study involving the combination of different levels of theory, basis sets and solvation models for an accurate sulfanamide pKa determination.
Collapse
|
7
|
Hu L, Huang B, Liu F. Atomistic Mechanism Underlying the Si(111)-(7×7) Surface Reconstruction Revealed by Artificial Neural-Network Potential. PHYSICAL REVIEW LETTERS 2021; 126:176101. [PMID: 33988452 DOI: 10.1103/physrevlett.126.176101] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/07/2020] [Revised: 10/09/2020] [Accepted: 03/31/2021] [Indexed: 06/12/2023]
Abstract
The 7×7 reconstruction of the Si(111) surface represents arguably the most fascinating surface reconstruction so far observed in nature. Yet, the atomistic mechanism underpinning its formation remains unclear after it was discovered sixty years ago. Experimentally, it is observed post priori so that analysis of its formation mechanism can only be carried out in analogy with archaeology. Theoretically, density-functional theory (DFT) correctly predicts the Si(111)-(7×7) ground state but is impractical to simulate its formation process; while empirical potentials failed to produce it as the ground state. Developing an artificial neural-network potential of DFT quality, we carried out accurate large-scale simulations to unravel the formation of the Si(111)-(7×7) surface. We reveal a possible step-mediated atom-pop rate-limiting process that triggers massive nonconserved atomic rearrangements, most remarkably, a critical process of collective vacancy diffusion that mediates a sequence of selective dimer, corner-hole, stacking-fault, and dimer-line pattern formation, to fulfill the 7×7 reconstruction. Our findings may not only solve the long-standing mystery of this famous surface reconstruction but they also illustrate the power of machine learning in studying complex structures.
Collapse
Affiliation(s)
- Lin Hu
- Beijing Computational Science Research Center, Beijing 100193, China
| | - Bing Huang
- Beijing Computational Science Research Center, Beijing 100193, China
| | - Feng Liu
- Department of Materials Science and Engineering, University of Utah, Salt Lake City, Utah 84112, USA
| |
Collapse
|
8
|
Adams DJ, Wang L, Steinle-Neumann G, Passerone D, Churakov SV. Anharmonic effects on the dynamics of solid aluminium from ab initiosimulations. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2021; 33:175501. [PMID: 33176283 DOI: 10.1088/1361-648x/abc972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/07/2020] [Accepted: 11/11/2020] [Indexed: 06/11/2023]
Abstract
Two approaches to simulations of phonon properties of solids beyond the harmonic approximation, the self-consistentab initiolattice dynamics (SCAILD) and decoupled anharmonic mode approximation (DAMA) are critically benchmarked against each other and molecular dynamics simulations using a density-functional-theory description of electronic states, and compared to experimental data for fcc aluminium. The temperature-dependence of phonon dispersion and the phonon density-of-states, heat capacity, and the mean atomic displacement for fcc aluminium are examined with these approaches at ambient pressure. A comparison of results obtained with the harmonic approximation to the ones predicted by SCAILD and DAMA reveal a negligible anharmonic contribution to phonon frequencies, a small, but significant influence on heat capacity, and a strong effect on atomic mean-square displacement. The phase space accessed with SCAILD and DAMA is reduced relative to molecular and harmonic lattice dynamics simulations. In particular the DAMA results are in good agreement with displacement amplitudes determined by the Debye-Waller factor in x-ray diffraction experiments.
Collapse
Affiliation(s)
| | - Lin Wang
- Bayerisches Geoinstitut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | | | - Daniele Passerone
- Empa-Swiss Federal Laboratories for Materials Science and Technology, CH-8600 Dübendorf, Switzerland
| | - Sergey V Churakov
- University of Bern, CH-3012 Bern, Switzerland
- Laboratory for Waste Management, Paul Scherrer Institute, CH-5232 Villigen-PSI, Switzerland
| |
Collapse
|
9
|
Mewes J, Schwerdtfeger P. Ausschließlich relativistisch: Periodische Trends in den Schmelz‐ und Siedepunkten der Gruppe 12. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202100486] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Jan‐Michael Mewes
- Mulliken Center for Theoretical Chemistry University of Bonn Beringstraße 4 53115 Bonn Deutschland
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics The New Zealand Institute for Advanced Study Massey University Auckland 0632 Auckland Neuseeland
| |
Collapse
|
10
|
Mewes JM, Schwerdtfeger P. Exclusively Relativistic: Periodic Trends in the Melting and Boiling Points of Group 12. Angew Chem Int Ed Engl 2021; 60:7703-7709. [PMID: 33576164 PMCID: PMC8048430 DOI: 10.1002/anie.202100486] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2021] [Indexed: 02/01/2023]
Abstract
First-principles simulations can advance our understanding of phase transitions but are often too costly for the heavier elements, which require a relativistic treatment. Addressing this challenge, we recently composed an indirect approach: A precise incremental calculation of absolute Gibbs energies for the solid and liquid with a relativistic Hamiltonian that enables an accurate determination of melting and boiling points (MPs and BPs). Here, we apply this approach to the Group 12 elements Zn, Cd, Hg, and Cn, whose MPs and BPs we calculate with a mean absolute deviation of only 5 % and 1 %, respectively, while we confirm the previously predicted liquid aggregate state of Cn. At a non-relativistic level of theory, we obtain surprisingly similar MPs and BPs of 650±30 K and 1250±20 K, suggesting that periodic trends in this group are exclusively relativistic in nature. Ultimately, we discuss these results and their implication for Groups 11 and 14.
Collapse
Affiliation(s)
- Jan-Michael Mewes
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstrasse 4, 53115, Bonn, Germany
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, 0632, Auckland, New Zealand
| |
Collapse
|
11
|
Jinnouchi R, Karsai F, Verdi C, Kresse G. First-principles hydration free energies of oxygenated species at water-platinum interfaces. J Chem Phys 2021; 154:094107. [PMID: 33685177 DOI: 10.1063/5.0036097] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The hydration free energy of atoms and molecules adsorbed at liquid-solid interfaces strongly influences the stability and reactivity of solid surfaces. However, its evaluation is challenging in both experiments and theories. In this work, a machine learning aided molecular dynamics method is proposed and applied to oxygen atoms and hydroxyl groups adsorbed on Pt(111) and Pt(100) surfaces in water. The proposed method adopts thermodynamic integration with respect to a coupling parameter specifying a path from well-defined non-interacting species to the fully interacting ones. The atomistic interactions are described by a machine-learned inter-atomic potential trained on first-principles data. The free energy calculated by the machine-learned potential is further corrected by using thermodynamic perturbation theory to provide the first-principles free energy. The calculated hydration free energies indicate that only the hydroxyl group adsorbed on the Pt(111) surface attains a hydration stabilization. The observed trend is attributed to differences in the adsorption site and surface morphology.
Collapse
Affiliation(s)
| | - Ferenc Karsai
- VASP Software GmbH, Sensengasse 8/16, 1090 Vienna, Austria
| | - Carla Verdi
- Computational Materials Physics, Faculty of Physics, University of Vienna, Sensengasse 8/12, 1090 Vienna, Austria
| | - Georg Kresse
- VASP Software GmbH, Sensengasse 8/16, 1090 Vienna, Austria
| |
Collapse
|
12
|
Smits OR, Mewes J, Jerabek P, Schwerdtfeger P. Oganesson: A Noble Gas Element That Is Neither Noble Nor a Gas. Angew Chem Int Ed Engl 2020; 59:23636-23640. [PMID: 32959952 PMCID: PMC7814676 DOI: 10.1002/anie.202011976] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2020] [Indexed: 11/07/2022]
Abstract
Oganesson (Og) is the last entry into the Periodic Table completing the seventh period of elements and group 18 of the noble gases. Only five atoms of Og have been successfully produced in nuclear collision experiments, with an estimate half-life for294 118 Og of 0. 69 + 0 . 64 - 0 . 22 ms.[1] With such a short lifetime, chemical and physical properties inevitably have to come from accurate relativistic quantum theory. Here, we employ two complementary computational approaches, namely parallel tempering Monte-Carlo (PTMC) simulations and first-principles thermodynamic integration (TI), both calibrated against a highly accurate coupled-cluster reference to pin-down the melting and boiling points of this super-heavy element. In excellent agreement, these approaches show Og to be a solid at ambient conditions with a melting point of ≈325 K. In contrast, calculations in the nonrelativistic limit reveal a melting point for Og of 220 K, suggesting a gaseous state as expected for a typical noble gas element. Accordingly, relativistic effects shift the solid-to-liquid phase transition by about 100 K.
Collapse
Affiliation(s)
- Odile R. Smits
- The New Zealand Institute for Advanced Study and the Institute for Natural and Mathematical ScienceMassey University (Albany)0632AucklandNew Zealand
| | - Jan‐Michael Mewes
- Mulliken Center for Theoretical ChemistryUniversity of BonnBeringstr. 453115BonnGermany
| | - Paul Jerabek
- Nanotechnology DepartmentHelmholtz-Zentrum GeesthachtMax-Planck-Straße 121502GeesthachtGermany
| | - Peter Schwerdtfeger
- The New Zealand Institute for Advanced Study and the Institute for Natural and Mathematical ScienceMassey University (Albany)0632AucklandNew Zealand
| |
Collapse
|
13
|
Smits OR, Mewes J, Jerabek P, Schwerdtfeger P. Oganesson: Ein Edelgas, das weder edel noch ein Gas ist. Angew Chem Int Ed Engl 2020. [DOI: 10.1002/ange.202011976] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Odile R. Smits
- The New Zealand Institute for Advanced Study and the Institute for Natural and Mathematical Science Massey University (Albany) 0632 Auckland Neuseeland
| | - Jan‐Michael Mewes
- Mulliken Center for Theoretical Chemistry University of Bonn Beringstr. 4 53115 Bonn Deutschland
| | - Paul Jerabek
- Nanotechnology Department Helmholtz-Zentrum Geesthacht Max-Planck-Straße 1 21502 Geesthacht Deutschland
| | - Peter Schwerdtfeger
- The New Zealand Institute for Advanced Study and the Institute for Natural and Mathematical Science Massey University (Albany) 0632 Auckland Neuseeland
| |
Collapse
|
14
|
Mewes JM, Smits OR. Accurate elemental boiling points from first principles. Phys Chem Chem Phys 2020; 22:24041-24050. [PMID: 33078780 DOI: 10.1039/d0cp02884c] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
The normal boiling point (NBP) is a fundamental property of liquids and marks the intersection of the Gibbs energies of the liquid and the gas-phase at ambient pressure. This work provides the first comprehensive demonstration of the calculation of boiling points of atomic liquids through first-principles molecular-dynamics simulations. To this end, thermodynamic integration (TDI) and perturbation theory (TPT) are combined with a density-functional theory (DFT) Hamiltonian, which provides absolute Gibbs energies, internal energies, and entropies of atomic liquids with an accuracy of a few meV/atom. Linear extrapolation to the intersection with the Gibbs energy of a non-interacting gas-phase eventually pins-down the NBPs. While these direct results can already be quite accurate, they are susceptible to a systematic over or underbinding of the employed density functional. It is shown how this dependency can be strongly reduced and the robustness of the method increased through a simple linear correction termed λ-scaling. Eventually, by carefully tuning of the technical parameters of the approach, the walltime per element is reduced from weeks to about a day (10-20k core-hours), enabling extensive testing for B, Al, Na, K, Ca, Sr, Ba, Mn, Cu, Xe, and Hg. This comprehensive benchmark demonstrates the excellent performance and robustness of the approach with a mean absolute deviation (MAD) of less than 2% from experimental NBPs and very similar accuracy for liquid entropies (MAD 2.3 J (mol K)-1, 2% relative). In some cases, the uncertainties in the predictions are several times smaller than the variation between literature values, allowing us to clear out long-standing ambiguities in the NBPs of B and Ba.
Collapse
Affiliation(s)
- Jan-Michael Mewes
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany. and Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, 0632 Auckland, New Zealand
| | - Odile R Smits
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, 0632 Auckland, New Zealand
| |
Collapse
|
15
|
Miao L, Wang LW. Liquid to crystal Si growth simulation using machine learning force field. J Chem Phys 2020; 153:074501. [DOI: 10.1063/5.0011163] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Ling Miao
- School of Optical and Electronic Information, Huazhong University of Science and Technology, Wuhan 430074, People’s Republic of China
| | - Lin-Wang Wang
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
16
|
Abstract
Elemental gallium possesses several intriguing properties, such as a low melting point, a density anomaly and an electronic structure in which covalent and metallic features coexist. In order to simulate this complex system, we construct an ab initio quality interaction potential by training a neural network on a set of density functional theory calculations performed on configurations generated in multithermal–multibaric simulations. Here we show that the relative equilibrium between liquid gallium, α-Ga, β-Ga, and Ga-II is well described. The resulting phase diagram is in agreement with the experimental findings. The local structure of liquid gallium and its nucleation into α-Ga and β-Ga are studied. We find that the formation of metastable β-Ga is kinetically favored over the thermodinamically stable α-Ga. Finally, we provide insight into the experimental observations of extreme undercooling of liquid Ga. Exploring nucleation processes of gallium by molecular simulation is extremely challenging due to its structural complexity. Here the authors demonstrate a neural network potential trained on a multithermal–multibaric DFT data for the study of the phase diagram of gallium in a wide temperature and pressure range.
Collapse
|
17
|
Adamyan VM, Bondarev VN, Zavalniuk VV. Graphene thermal break-down induced by anharmonic bending mode. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2019; 31:465401. [PMID: 31341096 DOI: 10.1088/1361-648x/ab3477] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
A novel approach is proposed allowing us to prove, self-consistently, that free-standing graphene reaching a certain temperature loses its mechanical stability resulting in abrupt breakdown, which can be interpreted as melting. Our study is based on the idea of the crucial role of the anomalously soft bending 'sound' mode in the jump transition of graphene from the state with relatively small bending fluctuations to a state with fluctuations close in amplitude to the graphene lattice constant. The acme of the developed theory is in establishing a quantitative relationship connecting the graphene elastic moduli of second, third (negative!), and fourth orders at the melting temperature T m that permits us to calculate T m. The results obtained lay a theoretical foundation for an analog of Lindemann criterion for graphene.
Collapse
Affiliation(s)
- V M Adamyan
- Department of Theoretical Physics, Odessa I.I. Mechnikov National University, 2 Dvoryanska St., Odessa 65026, Ukraine
| | | | | |
Collapse
|
18
|
Deringer VL, Caro MA, Csányi G. Machine Learning Interatomic Potentials as Emerging Tools for Materials Science. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2019; 31:e1902765. [PMID: 31486179 DOI: 10.1002/adma.201902765] [Citation(s) in RCA: 209] [Impact Index Per Article: 41.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/30/2019] [Revised: 07/26/2019] [Indexed: 05/22/2023]
Abstract
Atomic-scale modeling and understanding of materials have made remarkable progress, but they are still fundamentally limited by the large computational cost of explicit electronic-structure methods such as density-functional theory. This Progress Report shows how machine learning (ML) is currently enabling a new degree of realism in materials modeling: by "learning" electronic-structure data, ML-based interatomic potentials give access to atomistic simulations that reach similar accuracy levels but are orders of magnitude faster. A brief introduction to the new tools is given, and then, applications to some select problems in materials science are highlighted: phase-change materials for memory devices; nanoparticle catalysts; and carbon-based electrodes for chemical sensing, supercapacitors, and batteries. It is hoped that the present work will inspire the development and wider use of ML-based interatomic potentials in diverse areas of materials research.
Collapse
Affiliation(s)
- Volker L Deringer
- Department of Engineering, University of Cambridge, Cambridge, CB2 1PZ, UK
- Department of Chemistry, University of Cambridge, Cambridge, CB2 1EW, UK
| | - Miguel A Caro
- Department of Electrical Engineering and Automation and Department of Applied Physics, Aalto University, Espoo, 02150, Finland
| | - Gábor Csányi
- Department of Engineering, University of Cambridge, Cambridge, CB2 1PZ, UK
| |
Collapse
|
19
|
Mewes J, Smits OR, Kresse G, Schwerdtfeger P. Copernicium: A Relativistic Noble Liquid. Angew Chem Int Ed Engl 2019. [DOI: 10.1002/ange.201906966] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Jan‐Michael Mewes
- Centre for Theoretical Chemistry and PhysicsThe New Zealand Institute for Advanced StudyMassey University Auckland 0632 Auckland New Zealand
- Mulliken Center for Theoretical ChemistryUniversity of Bonn Beringstr. 4 53115 Bonn Germany
| | - Odile R. Smits
- Centre for Theoretical Chemistry and PhysicsThe New Zealand Institute for Advanced StudyMassey University Auckland 0632 Auckland New Zealand
| | - Georg Kresse
- University of ViennaFaculty of Physics and Center for Computational Materials Sciences Sensengasse 8/12 1090 Wien Austria
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and PhysicsThe New Zealand Institute for Advanced StudyMassey University Auckland 0632 Auckland New Zealand
| |
Collapse
|
20
|
Mewes JM, Smits OR, Kresse G, Schwerdtfeger P. Copernicium: A Relativistic Noble Liquid. Angew Chem Int Ed Engl 2019; 58:17964-17968. [PMID: 31596013 PMCID: PMC6916354 DOI: 10.1002/anie.201906966] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2019] [Revised: 10/01/2019] [Indexed: 11/05/2022]
Abstract
The chemical nature and aggregate state of superheavy copernicium (Cn) have been subject of speculation for many years. While strong relativistic effects render Cn chemically inert, which led Pitzer to suggest a noble-gas-like behavior in 1975, Eichler and co-workers in 2008 reported substantial interactions with a gold surface in atom-at-a-time experiments, suggesting a metallic character and a solid aggregate state. Herein, we explore the physicochemical properties of Cn by means of first-principles free-energy calculations, which confirm Pitzer's original hypothesis: With predicted melting and boiling points of 283±11 K and 340±10 K, Cn is indeed a volatile liquid and exhibits a density very similar to that of mercury. However, in stark contrast to mercury and the lighter Group 12 metals, we find bulk Cn to be bound by dispersion and to exhibit a large band gap of 6.4 eV, which is consistent with a noble-gas-like character. This non-group-conforming behavior is eventually traced back to strong scalar-relativistic effects, and in the non-relativistic limit, Cn appears as a common Group 12 metal.
Collapse
Affiliation(s)
- Jan-Michael Mewes
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, 0632, Auckland, New Zealand.,Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115, Bonn, Germany
| | - Odile R Smits
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, 0632, Auckland, New Zealand
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Sciences, Sensengasse 8/12, 1090, Wien, Austria
| | - Peter Schwerdtfeger
- Centre for Theoretical Chemistry and Physics, The New Zealand Institute for Advanced Study, Massey University Auckland, 0632, Auckland, New Zealand
| |
Collapse
|
21
|
Chehaibou B, Badawi M, Bučko T, Bazhirov T, Rocca D. Computing RPA Adsorption Enthalpies by Machine Learning Thermodynamic Perturbation Theory. J Chem Theory Comput 2019; 15:6333-6342. [DOI: 10.1021/acs.jctc.9b00782] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Bilal Chehaibou
- Université de Lorraine, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
- CNRS, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
| | - Michael Badawi
- Université de Lorraine, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
- CNRS, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
| | - Tomáš Bučko
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Mlynská Dolina, Ilkovičova 6, SK-84215 Bratislava, Slovakia
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Dúbravská cesta 9, SK-84236 Bratislava, Slovakia
| | - Timur Bazhirov
- Exabyte Inc., San Francisco, California 94103, United States
| | - Dario Rocca
- Université de Lorraine, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
- CNRS, LPCT, UMR 7019, 54506 Vandoeuvre-lès-Nancy, France
| |
Collapse
|
22
|
Bruneval F. Assessment of the Linearized GW Density Matrix for Molecules. J Chem Theory Comput 2019; 15:4069-4078. [DOI: 10.1021/acs.jctc.9b00333] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Fabien Bruneval
- DEN, Service de Recherches de Métallurgie Physique, CEA, Université Paris-Saclay, F-91191 Gif-sur-Yvette, France
| |
Collapse
|
23
|
Abstract
We propose modifications to the functional form of the Strongly Constrained and Appropriately Normed (SCAN) density functional to eliminate numerical instabilities. This is necessary to allow reliable, automatic generation of pseudopotentials (including projector augmented-wave potentials). The regularized SCAN is designed to match the original form very closely, and we show that its performance remains comparable.
Collapse
Affiliation(s)
- Albert P Bartók
- Rutherford Appleton Laboratory, Scientific Computing Department Science and Technology Facilities Council, Didcot OX11 0QX, United Kingdom
| | - Jonathan R Yates
- Department of Materials University of Oxford, Oxford OX1 3PH, United Kingdom
| |
Collapse
|