-
Oxygen Vacancy Formation Energy in Metal Oxides: High Throughput Computational Studies and Machine Learning Predictions
Authors:
Bianca Baldassarri,
Jiangang He,
Abhijith Gopakumar,
Sean Griesemer,
Adolfo J. A. Salgado-Casanova,
Tzu-Chen Liu,
Steven B. Torrisi,
Chris Wolverton
Abstract:
The oxygen vacancy formation energy ($ΔE_{vf}$) governs defect dynamics and is a useful metric to perform materials selection for a variety of applications. However, density functional theory (DFT) calculations of $ΔE_{vf}$ come at a greater computational cost than the typical bulk calculations available in materials databases due to the involvement of multiple vacancy-containing supercells. As a…
▽ More
The oxygen vacancy formation energy ($ΔE_{vf}$) governs defect dynamics and is a useful metric to perform materials selection for a variety of applications. However, density functional theory (DFT) calculations of $ΔE_{vf}$ come at a greater computational cost than the typical bulk calculations available in materials databases due to the involvement of multiple vacancy-containing supercells. As a result, available repositories of direct calculations of $ΔE_{vf}$ remain relatively scarce, and the development of machine learning models capable of delivering accurate predictions is of interest. In the present, work we address both such points. We first report the results of new high-throughput DFT calculations of oxygen vacancy formation energies of the different unique oxygen sites in over 1000 different oxide materials, which together form the largest dataset of directly computed oxygen vacancy formation energies to date, to our knowledge. We then utilize the resulting dataset of $\sim$2500 $ΔE_{vf}$ values to train random forest models with different sets of features, examining both novel features introduced in this work and ones previously employed in the literature. We demonstrate the benefits of including features that contain information specific to the vacancy site and account for both cation identity and oxidation state, and achieve a mean absolute error upon prediction of $\sim$0.3 eV/O, which is comparable to the accuracy observed upon comparison of DFT computations of oxygen vacancy formation energy and experimental results. Finally, we demonstrate the predictive power of the developed models in the search for new compounds for solar-thermochemical water-splitting applications, finding over 250 new AA$^{\prime}$BB$^{\prime}$O$_6$ double perovskite candidates.
△ Less
Submitted 3 September, 2023;
originally announced September 2023.
-
Complexity of Many-Body Interactions in Transition Metals via Machine-Learned Force Fields from the TM23 Data Set
Authors:
Cameron J. Owen,
Steven B. Torrisi,
Yu Xie,
Simon Batzner,
Kyle Bystrom,
Jennifer Coulter,
Albert Musaelian,
Lixin Sun,
Boris Kozinsky
Abstract:
This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases of d-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE)…
▽ More
This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases of d-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupied d states near the Fermi level, we determine that the large, sharp d density of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.
△ Less
Submitted 26 September, 2023; v1 submitted 25 February, 2023;
originally announced February 2023.
-
High-Throughput $\textit{Ab Initio}$ Design of Atomic Interfaces using InterMatch
Authors:
Eli Gerber,
Steven B. Torrisi,
Sara Shabani,
Eric Seewald,
Jordan Pack,
Jennifer E. Hoffman,
Cory R. Dean,
Abhay N. Pasupathy,
Eun-Ah Kim
Abstract:
Forming a hetero-interface is a materials-design strategy that can access an astronomically large phase space. However, the immense phase space necessitates a high-throughput approach for optimal interface design. Here we introduce a high-throughput computational framework, InterMatch, for efficiently predicting charge transfer, strain, and superlattice structure of an interface by leveraging the…
▽ More
Forming a hetero-interface is a materials-design strategy that can access an astronomically large phase space. However, the immense phase space necessitates a high-throughput approach for optimal interface design. Here we introduce a high-throughput computational framework, InterMatch, for efficiently predicting charge transfer, strain, and superlattice structure of an interface by leveraging the databases of individual bulk materials. Specifically, the algorithm reads in the lattice vectors, density of states, and the stiffness tensors for each material in their isolated form from the Materials Project. From these bulk properties, InterMatch estimates the interfacial properties. We benchmark InterMatch predictions for the charge transfer against experimental measurements and supercell density-functional theory calculations. We then use InterMatch to predict promising interface candidates for doping transition metal dichalcogenide MoSe$_2$. Finally, we explain experimental observation of factor of 10 variation in the supercell periodicity within a few microns in graphene/$α$-RuCl$_3$ by exploring low energy superlattice structures as a function of twist angle using InterMatch. We anticipate our open-source InterMatch algorithm accelerating and guiding ever-growing interfacial design efforts. Moreover, the interface database resulting from the InterMatch searches presented in this paper can be readily accessed through https://contribs.materialsproject.org/projects/intermatch/ .
△ Less
Submitted 29 November, 2022;
originally announced November 2022.
-
Effects of Structural Distortions on the Electronic Structure of T-type Transition Metal Dichalcogenides
Authors:
Daniel T. Larson,
Wei Chen,
Steven B. Torrisi,
Jennifer Coulter,
Shiang Fang,
Efthimios Kaxiras
Abstract:
Single-layer transition metal dichalcogenides (TMDCs) can adopt two distinct structures corresponding to different coordination of the metal atoms. TMDCs adopting the T-type structure exhibit a rich and diverse set of phenomena, including charge density waves (CDW) in a $\sqrt{13}\times\sqrt{13}$ supercell pattern in TaS$_2$ and TaSe$_2$, and a possible excitonic insulating phase in TiSe$_2$. Thes…
▽ More
Single-layer transition metal dichalcogenides (TMDCs) can adopt two distinct structures corresponding to different coordination of the metal atoms. TMDCs adopting the T-type structure exhibit a rich and diverse set of phenomena, including charge density waves (CDW) in a $\sqrt{13}\times\sqrt{13}$ supercell pattern in TaS$_2$ and TaSe$_2$, and a possible excitonic insulating phase in TiSe$_2$. These properties make the T-TMDCs desirable components of layered heterostructure devices. In order to predict the emergent properties of combinations of different layered materials, one needs simple and accurate models for the constituent layers which can take into account potential effects of lattice mismatch, relaxation, strain, and structural distortion. Previous studies have developed ab initio tight-binding Hamiltonians for H-type TMDCs [arXiv:1709.07510]. Here we extend this work to include T-type TMDCs. We demonstrate the capabilities of our model using three example systems: a 1-dimensional sinusoidal ripple, the 2$\times$2 CDW in TiSe$_2$, and the $\sqrt{13}\times\sqrt{13}$ CDW in TaS$_2$. Using the technique of band unfolding we compare the electronic structure of the distorted crystals to the pristine band structure and find excellent agreement with direct DFT calculations, provided the magnitude of the distortions remains in the linear regime.
△ Less
Submitted 23 March, 2020;
originally announced March 2020.
-
Electronic structure calculations of twisted multi-layer graphene superlattices
Authors:
Georgios A. Tritsaris,
Stephen Carr,
Ziyan Zhu,
Yiqi Xie,
Steven B. Torrisi,
Jing Tang,
Marios Mattheakis,
Daniel Larson,
Efthimios Kaxiras
Abstract:
Quantum confinement endows two-dimensional (2D) layered materials with exceptional physics and novel properties compared to their bulk counterparts. Although certain two- and few-layer configurations of graphene have been realized and studied, a systematic investigation of the properties of arbitrarily layered graphene assemblies is still lacking. We introduce theoretical concepts and methods for…
▽ More
Quantum confinement endows two-dimensional (2D) layered materials with exceptional physics and novel properties compared to their bulk counterparts. Although certain two- and few-layer configurations of graphene have been realized and studied, a systematic investigation of the properties of arbitrarily layered graphene assemblies is still lacking. We introduce theoretical concepts and methods for the processing of materials information, and as a case study, apply them to investigate the electronic structure of multi-layer graphene-based assemblies in a high-throughput fashion. We provide a critical discussion of patterns and trends in tight binding band structures and we identify specific layered assemblies using low-dispersion electronic bands as indicators of potentially interesting physics like strongly correlated behavior. A combination of data-driven models for visualization and prediction is used to intelligently explore the materials space. This work more generally aims to increase confidence in the combined use of physics-based and data-driven modeling for the systematic refinement of knowledge about 2D layered materials, with implications for the development of novel quantum devices.
△ Less
Submitted 30 January, 2020;
originally announced January 2020.
-
Two-Dimensional Forms of Robust CO$_2$ Reduction Photocatalysts
Authors:
Steven B. Torrisi,
Arunima K. Singh,
Joseph H. Montoya,
Tathagata Biswas,
Kristin A. Persson
Abstract:
Novel photoelectrocatalysts that use sunlight to power the CO$_2$ reduction reaction will be crucial for carbon-neutral power and energy-efficient industrial processes. Scalable photoelectrocatalysts must satisfy a stringent set of criteria, such as stability under operating conditions, product selectivity, and efficient light absorption. Two-dimensional materials can offer high specific surface a…
▽ More
Novel photoelectrocatalysts that use sunlight to power the CO$_2$ reduction reaction will be crucial for carbon-neutral power and energy-efficient industrial processes. Scalable photoelectrocatalysts must satisfy a stringent set of criteria, such as stability under operating conditions, product selectivity, and efficient light absorption. Two-dimensional materials can offer high specific surface area, tunability, and potential for heterostructuring, providing a fresh landscape of candidate catalysts. From a set of promising bulk CO$_2$ reduction photoelectrocatalysts, we screen for candidate mono-layers of these materials, then study their catalytic feasibility and suitability. For stable monolayer candidates, we verify the presence of visible-light band gaps, check that band edges can support CO$_2$ reduction, determine exciton binding energies, and compute surface reactivity. We find for SiAs, ZnTe, and ZnSe monolayers, visible light absorption is possible, and reaction selectivity biases towards CO production. We thus identify SiAs, ZnTe, and ZnSe monolayers as targets for further investigation, expanding the chemical space for CO$_2$ photoreduction.
△ Less
Submitted 19 December, 2019;
originally announced December 2019.
-
On-the-Fly Active Learning of Interpretable Bayesian Force Fields for Atomistic Rare Events
Authors:
Jonathan Vandermause,
Steven B. Torrisi,
Simon Batzner,
Yu Xie,
Lixin Sun,
Alexie M. Kolpak,
Boris Kozinsky
Abstract:
Machine learned force fields typically require manual construction of training sets consisting of thousands of first principles calculations, which can result in low training efficiency and unpredictable errors when applied to structures not represented in the training set of the model. This severely limits the practical application of these models in systems with dynamics governed by important ra…
▽ More
Machine learned force fields typically require manual construction of training sets consisting of thousands of first principles calculations, which can result in low training efficiency and unpredictable errors when applied to structures not represented in the training set of the model. This severely limits the practical application of these models in systems with dynamics governed by important rare events, such as chemical reactions and diffusion. We present an adaptive Bayesian inference method for automating the training of interpretable, low-dimensional, and multi-element interatomic force fields using structures drawn on the fly from molecular dynamics simulations. Within an active learning framework, the internal uncertainty of a Gaussian process regression model is used to decide whether to accept the model prediction or to perform a first principles calculation to augment the training set of the model. The method is applied to a range of single- and multi-element systems and shown to achieve a favorable balance of accuracy and computational efficiency, while requiring a minimal amount of ab initio training data. We provide a fully open-source implementation of our method, as well as a procedure to map trained models to computationally efficient tabulated force fields.
△ Less
Submitted 19 November, 2019; v1 submitted 3 April, 2019;
originally announced April 2019.
-
Relaxation and Domain Formation in Incommensurate 2D Heterostructures
Authors:
Stephen Carr,
Daniel Massatt,
Steven B. Torrisi,
Paul Cazeaux,
Mitchell Luskin,
Efthimios Kaxiras
Abstract:
We introduce configuration space as a natural representation for calculating the mechanical relaxation patterns of incommensurate two-dimensional (2D) bilayers, bypassing supercell approximations to encompass aperiodic relaxation patterns. The approach can be applied to a wide variety of 2D materials through the use of a continuum model in combination with a generalized stacking fault energy for i…
▽ More
We introduce configuration space as a natural representation for calculating the mechanical relaxation patterns of incommensurate two-dimensional (2D) bilayers, bypassing supercell approximations to encompass aperiodic relaxation patterns. The approach can be applied to a wide variety of 2D materials through the use of a continuum model in combination with a generalized stacking fault energy for interlayer interactions. We present computational results for small-angle twisted bilayer graphene and molybdenum disulfide (MoS$_2$), a representative material of the transition metal dichalcogenide (TMDC) family of 2D semiconductors. We calculate accurate relaxations for MoS$_2$ even at small twist-angle values, enabled by the fact that our approach does not rely on empirical atomistic potentials for interlayer coupling. The results demonstrate the efficiency of the configuration space method by computing relaxations with minimal computational cost for twist angles down to $0.05^\circ$, which is smaller than what can be explored by any available real space techniques. We also outline a general explanation of domain formation in 2D bilayers with nearly-aligned lattices, taking advantage of the relationship between real space and configuration space.
△ Less
Submitted 17 May, 2018;
originally announced May 2018.
-
Perpendicular laser cooling with a rotating wall potential in a Penning trap
Authors:
Steven B. Torrisi,
Joseph W. Britton,
Justin G. Bohnet,
John J. Bollinger
Abstract:
We investigate the impact of a rotating wall potential on perpendicular laser cooling in a Penning ion trap. By including energy exchange with the rotating wall, we extend previous Doppler laser cooling theory and show that low perpendicular temperatures are more readily achieved with a rotating wall than without. Detailed numerical studies determine optimal operating parameters for producing low…
▽ More
We investigate the impact of a rotating wall potential on perpendicular laser cooling in a Penning ion trap. By including energy exchange with the rotating wall, we extend previous Doppler laser cooling theory and show that low perpendicular temperatures are more readily achieved with a rotating wall than without. Detailed numerical studies determine optimal operating parameters for producing low temperature, stable 2-dimensional crystals, important for quantum information processing experiments employing Penning traps.
△ Less
Submitted 29 May, 2016; v1 submitted 16 February, 2016;
originally announced February 2016.