1
|
Selloni A. Aqueous Titania Interfaces. Annu Rev Phys Chem 2024; 75:47-65. [PMID: 38271659 DOI: 10.1146/annurev-physchem-090722-015957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2024]
Abstract
Water-metal oxide interfaces are central to many phenomena and applications, ranging from material corrosion and dissolution to photoelectrochemistry and bioengineering. In particular, the discovery of photocatalytic water splitting on TiO2 has motivated intensive studies of water-TiO2 interfaces for decades. So far, a broad understanding of the interaction of water vapor with several TiO2 surfaces has been obtained. However, much less is known about liquid water-TiO2 interfaces, which are more relevant to many practical applications. Probing these complex systems at the molecular level is experimentally challenging and is sometimes possible only through computational studies. This review summarizes recent advances in the atomistic understanding, mostly through computational simulations, of the structure and dynamics of interfacial water on TiO2 surfaces. The main focus is on the nature, molecular or dissociated, of water in direct contact with low-index defect-free crystalline surfaces. The hydroxyls resulting from water dissociation are essential in the photooxidation of water and critically affect the surface chemistry of TiO2.
Collapse
Affiliation(s)
- Annabella Selloni
- Department of Chemistry, Princeton University, Princeton, New Jersey, USA;
| |
Collapse
|
2
|
Liu M, Wang J, Hu J, Liu P, Niu H, Yan X, Li J, Yan H, Yang B, Sun Y, Chen C, Kresse G, Zuo L, Chen XQ. Layer-by-layer phase transformation in Ti 3O 5 revealed by machine-learning molecular dynamics simulations. Nat Commun 2024; 15:3079. [PMID: 38594273 PMCID: PMC11004112 DOI: 10.1038/s41467-024-47422-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 03/28/2024] [Indexed: 04/11/2024] Open
Abstract
Reconstructive phase transitions involving breaking and reconstruction of primary chemical bonds are ubiquitous and important for many technological applications. In contrast to displacive phase transitions, the dynamics of reconstructive phase transitions are usually slow due to the large energy barrier. Nevertheless, the reconstructive phase transformation from β- to λ-Ti3O5 exhibits an ultrafast and reversible behavior. Despite extensive studies, the underlying microscopic mechanism remains unclear. Here, we discover a kinetically favorable in-plane nucleated layer-by-layer transformation mechanism through metadynamics and large-scale molecular dynamics simulations. This is enabled by developing an efficient machine learning potential with near first-principles accuracy through an on-the-fly active learning method and an advanced sampling technique. Our results reveal that the β-λ phase transformation initiates with the formation of two-dimensional nuclei in the ab-plane and then proceeds layer-by-layer through a multistep barrier-lowering kinetic process via intermediate metastable phases. Our work not only provides important insight into the ultrafast and reversible nature of the β-λ transition, but also presents useful strategies and methods for tackling other complex structural phase transitions.
Collapse
Affiliation(s)
- Mingfeng Liu
- Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang, 110016, China
- School of Materials Science and Engineering, University of Science and Technology of China, Shenyang, 110016, China
| | - Jiantao Wang
- Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang, 110016, China
- School of Materials Science and Engineering, University of Science and Technology of China, Shenyang, 110016, China
| | - Junwei Hu
- State Key Laboratory of Solidification Processing, International Center for Materials Discovery, School of Materials Science and Engineering, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Peitao Liu
- Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang, 110016, China.
| | - Haiyang Niu
- State Key Laboratory of Solidification Processing, International Center for Materials Discovery, School of Materials Science and Engineering, Northwestern Polytechnical University, Xi'an, 710072, China.
| | - Xuexi Yan
- Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang, 110016, China
| | - Jiangxu Li
- Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang, 110016, China
| | - Haile Yan
- Key Laboratory for Anisotropy and Texture of Materials (Ministry of Education), School of Materials Science and Engineering, Northeastern University, Shenyang, 110819, China
| | - Bo Yang
- Key Laboratory for Anisotropy and Texture of Materials (Ministry of Education), School of Materials Science and Engineering, Northeastern University, Shenyang, 110819, China
| | - Yan Sun
- Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang, 110016, China
| | - Chunlin Chen
- Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang, 110016, China
| | - Georg Kresse
- University of Vienna, Faculty of Physics and Center for Computational Materials Science, Kolingasse 14-16, A-1090, Vienna, Austria
| | - Liang Zuo
- Key Laboratory for Anisotropy and Texture of Materials (Ministry of Education), School of Materials Science and Engineering, Northeastern University, Shenyang, 110819, China
| | - Xing-Qiu Chen
- Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang, 110016, China
| |
Collapse
|
3
|
Sheldon C, Paier J, Usvyat D, Sauer J. Hybrid RPA:DFT Approach for Adsorption on Transition Metal Surfaces: Methane and Ethane on Platinum (111). J Chem Theory Comput 2024; 20:2219-2227. [PMID: 38330551 PMCID: PMC10938501 DOI: 10.1021/acs.jctc.3c01308] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2023] [Revised: 01/17/2024] [Accepted: 01/18/2024] [Indexed: 02/10/2024]
Abstract
The hybrid QM:QM approach is extended to adsorption on transition metal surfaces. The random phase approximation (RPA) as the high-level method is applied to cluster models and, using the subtractive scheme, embedded in periodic models which are treated with density functional theory (DFT) that is the low-level method. The PBE functional, both without dispersion and augmented with the many-body dispersion (MBD), is employed. Adsorption of methane and ethane on the Pt(111) surface is studied. For methane in a 2 × 2 surface cell, the hybrid RPA:PBE and RPA:PBE+MBD results, -14.3 and -16.0 kJ mol-1, respectively, are in close agreement with the periodic RPA value of -13.8 kJ mol-1 at significantly reduced computational cost (factor of ∼50). For methane and ethane, the RPA:PBE results (-14.3 and -17.8 kJ mol-1, respectively) indicate underbinding relative to energies derived from experimental desorption barriers for relevant loadings (-15.6 ± 1.6 and -27.2 ± 2.9 kJ mol-1, respectively), whereas the hybrid RPA:PBE+MBD results (-16.0 and -24.9 kJ mol-1, respectively) agree with the experiment well within experimental uncertainty limits (deviation of -0.4 ± 1.5 and +2.3 ± 2.9 kJ mol-1, respectively). Finding a cluster that adequately and robustly represents the adsorbate at the bulk surface is important for the success of the RPA-based QM:QM scheme for metals.
Collapse
Affiliation(s)
- Christopher Sheldon
- Institut
für Chemie, Humboldt-Universität
zu Berlin, Unter den Linden 6, Berlin 10099, Germany
- Fritz-Haber-Institut
der Max-Planck-Gesellschaft, Faradayweg 4, Berlin 14195, Germany
| | - Joachim Paier
- Institut
für Chemie, Humboldt-Universität
zu Berlin, Unter den Linden 6, Berlin 10099, Germany
- Lehrstuhl
für Theoretische Chemie, Friedrich-Alexander-Universität
Erlangen-Nürnberg, Egerlandstrasse 3, Erlangen 91058, Germany
| | - Denis Usvyat
- Institut
für Chemie, Humboldt-Universität
zu Berlin, Unter den Linden 6, Berlin 10099, Germany
| | - Joachim Sauer
- Institut
für Chemie, Humboldt-Universität
zu Berlin, Unter den Linden 6, Berlin 10099, Germany
| |
Collapse
|
4
|
Butler PV, Hafizi R, Day GM. Machine-Learned Potentials by Active Learning from Organic Crystal Structure Prediction Landscapes. J Phys Chem A 2024; 128:945-957. [PMID: 38277275 PMCID: PMC10860135 DOI: 10.1021/acs.jpca.3c07129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2023] [Revised: 01/04/2024] [Accepted: 01/11/2024] [Indexed: 01/28/2024]
Abstract
A primary challenge in organic molecular crystal structure prediction (CSP) is accurately ranking the energies of potential structures. While high-level solid-state density functional theory (DFT) methods allow for mostly reliable discrimination of the low-energy structures, their high computational cost is problematic because of the need to evaluate tens to hundreds of thousands of trial crystal structures to fully explore typical crystal energy landscapes. Consequently, lower-cost but less accurate empirical force fields are often used, sometimes as the first stage of a hierarchical scheme involving multiple stages of increasingly accurate energy calculations. Machine-learned interatomic potentials (MLIPs), trained to reproduce the results of ab initio methods with computational costs close to those of force fields, can improve the efficiency of the CSP by reducing or eliminating the need for costly DFT calculations. Here, we investigate active learning methods for training MLIPs with CSP datasets. The combination of active learning with the well-developed sampling methods from CSP yields potentials in a highly automated workflow that are relevant over a wide range of the crystal packing space. To demonstrate these potentials, we illustrate efficiently reranking large, diverse crystal structure landscapes to near-DFT accuracy from force field-based CSP, improving the reliability of the final energy ranking. Furthermore, we demonstrate how these potentials can be extended to more accurately model structures far from lattice energy minima through additional on-the-fly training within Monte Carlo simulations.
Collapse
Affiliation(s)
| | - Roohollah Hafizi
- School of Chemistry, University
of Southampton, Southampton SO17 1BJ, U.K.
| | - Graeme M. Day
- School of Chemistry, University
of Southampton, Southampton SO17 1BJ, U.K.
| |
Collapse
|
5
|
Cheng H, Jiao P, Wang J, Qing M, Deng Y, Liu JM, Bellaiche L, Wu D, Yang Y. Tunable and parabolic piezoelectricity in hafnia under epitaxial strain. Nat Commun 2024; 15:394. [PMID: 38195734 PMCID: PMC10776838 DOI: 10.1038/s41467-023-44207-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2023] [Accepted: 12/04/2023] [Indexed: 01/11/2024] Open
Abstract
Piezoelectrics are a class of functional materials that have been extensively used for application in modern electro-mechanical and mechatronics technologies. The sign of longitudinal piezoelectric coefficients is typically positive but recently a few ferroelectrics, such as ferroelectric polymer poly(vinylidene fluoride) and van der Waals ferroelectric CuInP2S6, were experimentally found to have negative piezoelectricity. Here, using first-principles calculation and measurements, we show that the sign of the longitudinal linear piezoelectric coefficient of HfO2 can be tuned from positive to negative via epitaxial strain. Nonlinear and even parabolic piezoelectric behaviors are further found at tensile epitaxial strain. This parabolic piezoelectric behavior implies that the polarization decreases when increasing the magnitude of either compressive or tensile longitudinal strain, or, equivalently, that the strain increases when increasing the magnitude of electric field being either parallel or antiparallel to the direction of polarization. The unusual piezoelectric effects are from the chemical coordination of the active oxygen atoms. These striking piezoelectric features of positive and negative sign, as well as linear and parabolical behaviors, expand the current knowledge in piezoelectricity and broaden the potential of piezoelectric applications towards electro-mechanical and communications technology.
Collapse
Affiliation(s)
- Hao Cheng
- Laboratory of Solid State Microstructures, Nanjing University, Nanjing, 210093, China
- Jiangsu Key Laboratory of Artificial Functional Materials, Department of Materials Science and Engineering, Nanjing University, Nanjing, 210093, China
| | - Peijie Jiao
- Laboratory of Solid State Microstructures, Nanjing University, Nanjing, 210093, China
- Jiangsu Key Laboratory of Artificial Functional Materials, Department of Materials Science and Engineering, Nanjing University, Nanjing, 210093, China
| | - Jian Wang
- Laboratory of Solid State Microstructures, Nanjing University, Nanjing, 210093, China
- Jiangsu Key Laboratory of Artificial Functional Materials, Department of Materials Science and Engineering, Nanjing University, Nanjing, 210093, China
| | - Mingkai Qing
- Laboratory of Solid State Microstructures, Nanjing University, Nanjing, 210093, China
- Jiangsu Key Laboratory of Artificial Functional Materials, Department of Materials Science and Engineering, Nanjing University, Nanjing, 210093, China
| | - Yu Deng
- Laboratory of Solid State Microstructures, Nanjing University, Nanjing, 210093, China
- Jiangsu Key Laboratory of Artificial Functional Materials, Department of Materials Science and Engineering, Nanjing University, Nanjing, 210093, China
| | - Jun-Ming Liu
- Laboratory of Solid State Microstructures, Nanjing University, Nanjing, 210093, China
| | - Laurent Bellaiche
- Physics Department, Institute for Nanoscience and Engineering, University of Arkansas, Fayetteville, AR, 72701, USA.
| | - Di Wu
- Laboratory of Solid State Microstructures, Nanjing University, Nanjing, 210093, China.
- Jiangsu Key Laboratory of Artificial Functional Materials, Department of Materials Science and Engineering, Nanjing University, Nanjing, 210093, China.
| | - Yurong Yang
- Laboratory of Solid State Microstructures, Nanjing University, Nanjing, 210093, China.
- Jiangsu Key Laboratory of Artificial Functional Materials, Department of Materials Science and Engineering, Nanjing University, Nanjing, 210093, China.
| |
Collapse
|
6
|
Stoppelman JP, Wilkinson AP, McDaniel JG. Equation of state predictions for ScF3 and CaZrF6 with neural network-driven molecular dynamics. J Chem Phys 2023; 159:084707. [PMID: 37638627 DOI: 10.1063/5.0157615] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Accepted: 08/09/2023] [Indexed: 08/29/2023] Open
Abstract
In silico property prediction based on density functional theory (DFT) is increasingly performed for crystalline materials. Whether quantitative agreement with experiment can be achieved with current methods is often an unresolved question, and may require detailed examination of physical effects such as electron correlation, reciprocal space sampling, phonon anharmonicity, and nuclear quantum effects (NQE), among others. In this work, we attempt first-principles equation of state prediction for the crystalline materials ScF3 and CaZrF6, which are known to exhibit negative thermal expansion (NTE) over a broad temperature range. We develop neural network (NN) potentials for both ScF3 and CaZrF6 trained to extensive DFT data, and conduct direct molecular dynamics prediction of the equation(s) of state over a broad temperature/pressure range. The NN potentials serve as surrogates of the DFT Hamiltonian with enhanced computational efficiency allowing for simulations with larger supercells and inclusion of NQE utilizing path integral approaches. The conclusion of the study is mixed: while some equation of state behavior is predicted in semiquantitative agreement with experiment, the pressure-induced softening phenomenon observed for ScF3 is not captured in our simulations. We show that NQE have a moderate effect on NTE at low temperature but does not significantly contribute to equation of state predictions at increasing temperature. Overall, while the NN potentials are valuable for property prediction of these NTE (and related) materials, we infer that a higher level of electron correlation, beyond the generalized gradient approximation density functional employed here, is necessary for achieving quantitative agreement with experiment.
Collapse
Affiliation(s)
- John P Stoppelman
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Angus P Wilkinson
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
- School of Materials Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0245, USA
| | - Jesse G McDaniel
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
7
|
Stenczel TK, El-Machachi Z, Liepuoniute G, Morrow JD, Bartók AP, Probert MIJ, Csányi G, Deringer VL. Machine-learned acceleration for molecular dynamics in CASTEP. J Chem Phys 2023; 159:044803. [PMID: 37497818 DOI: 10.1063/5.0155621] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2023] [Accepted: 07/03/2023] [Indexed: 07/28/2023] Open
Abstract
Machine learning (ML) methods are of rapidly growing interest for materials modeling, and yet, the use of ML interatomic potentials for new systems is often more demanding than that of established density-functional theory (DFT) packages. Here, we describe computational methodology to combine the CASTEP first-principles simulation software with the on-the-fly fitting and evaluation of ML interatomic potential models. Our approach is based on regular checking against DFT reference data, which provides a direct measure of the accuracy of the evolving ML model. We discuss the general framework and the specific solutions implemented, and we present an example application to high-temperature molecular-dynamics simulations of carbon nanostructures. The code is freely available for academic research.
Collapse
Affiliation(s)
- Tamás K Stenczel
- Engineering Laboratory, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
| | - Zakariya El-Machachi
- Department of Chemistry, Inorganic Chemistry Laboratory, University of Oxford, Oxford OX1 3QR, United Kingdom
| | - Guoda Liepuoniute
- Engineering Laboratory, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
| | - Joe D Morrow
- Department of Chemistry, Inorganic Chemistry Laboratory, University of Oxford, Oxford OX1 3QR, United Kingdom
| | - Albert P Bartók
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
- Warwick Centre for Predictive Modelling, School of Engineering, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Matt I J Probert
- School of Physics, Engineering and Technology, University of York, York YO10 5DD, United Kingdom
| | - Gábor Csányi
- Engineering Laboratory, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
| | - Volker L Deringer
- Department of Chemistry, Inorganic Chemistry Laboratory, University of Oxford, Oxford OX1 3QR, United Kingdom
| |
Collapse
|