-
Accelerating ab initio melting property calculations with machine learning: Application to the high entropy alloy TaVCrW
Authors:
Li-Fang Zhu,
Fritz Koermann,
Qing Chen,
Malin Selleby,
Joerg Neugebauer,
and Blazej Grabowski
Abstract:
Melting properties are critical for designing novel materials, especially for discovering high-performance, high-melting refractory materials. Experimental measurements of these properties are extremely challenging due to their high melting temperatures. Complementary theoretical predictions are, therefore, indispensable. The conventional free energy approach using density functional theory (DFT)…
▽ More
Melting properties are critical for designing novel materials, especially for discovering high-performance, high-melting refractory materials. Experimental measurements of these properties are extremely challenging due to their high melting temperatures. Complementary theoretical predictions are, therefore, indispensable. The conventional free energy approach using density functional theory (DFT) has been a gold standard for such purposes because of its high accuracy. However,it generally involves expensive thermodynamic integration using ab initio molecular dynamic simulations. The high computational cost makes high-throughput calculations infeasible. Here, we propose a highly efficient DFT-based method aided by a specially designed machine learning potential. As the machine learning potential can closely reproduce the ab initio phase space, even for multi-component alloys, the costly thermodynamic integration can be fully substituted with more efficient free energy perturbation calculations. The method achieves overall savings of computational resources by 80% compared to current alternatives. We apply the method to the high-entropy alloy TaVCrW and calculate its melting properties, including melting temperature, entropy and enthalpy of fusion, and volume change at the melting point. Additionally, the heat capacities of solid and liquid TaVCrW are calculated. The results agree reasonably with the calphad extrapolated values.
△ Less
Submitted 16 August, 2024;
originally announced August 2024.
-
Dissecting physics of carbon ordering in bcc iron
Authors:
Sam Waseda,
Tilmann Hickel,
Julien Morthomas,
Patrice Chantrenne,
Michel Perez,
Jörg Neugebauer
Abstract:
Zener ordering is a phenomenon that octahedral interstitial atoms such as carbon occupy the same sublattice inside bcc matrix such as iron. The original formulation relies on a mean field theory, which is still most in use today. We employ multiple methods, such as Molecular Dynamics, Metropolis Monte Carlo, Mean Field Theory with chemical interactions and finite temperature effects to show that t…
▽ More
Zener ordering is a phenomenon that octahedral interstitial atoms such as carbon occupy the same sublattice inside bcc matrix such as iron. The original formulation relies on a mean field theory, which is still most in use today. We employ multiple methods, such as Molecular Dynamics, Metropolis Monte Carlo, Mean Field Theory with chemical interactions and finite temperature effects to show that the Zener ordering for iron carbon systems is governed by local chemical interactions and finite temperature effects and less of mean field nature as described originally by Zener.
△ Less
Submitted 30 July, 2024;
originally announced July 2024.
-
Development of an Atomic Cluster Expansion potential for iron and its oxides
Authors:
Baptiste Bienvenu,
Mira Todorova,
Jörg Neugebauer,
Dierk Raabe,
Matous Mrovec,
Yury Lysogorskiy,
Ralf Drautz
Abstract:
The combined structural and electronic complexity of iron oxides poses many challenges to atomistic modeling. To leverage limitations in terms of the accessible length and time scales, one requires a physically justified interatomic potential which is accurate to correctly account for the complexity of iron-oxygen systems. Such a potential is not yet available in the literature. In this work, we p…
▽ More
The combined structural and electronic complexity of iron oxides poses many challenges to atomistic modeling. To leverage limitations in terms of the accessible length and time scales, one requires a physically justified interatomic potential which is accurate to correctly account for the complexity of iron-oxygen systems. Such a potential is not yet available in the literature. In this work, we propose a machine-learning potential based on the Atomic Cluster Expansion for modeling the iron-oxygen system, which explicitly accounts for magnetism. We test the potential on a wide range of properties of iron and its oxides, and demonstrate its ability to describe the thermodynamics of systems spanning the whole range of oxygen content and including magnetic degrees of freedom.
△ Less
Submitted 18 July, 2024;
originally announced July 2024.
-
On the origin of univalent Mg$^+$ ions in solution and their role in anomalous anodic hydrogen evolution
Authors:
Florian Deißenbeck,
Sudarsan Surendralal,
Mira Todorova,
Stefan Wippermann,
Jörg Neugebauer
Abstract:
Aqueous metal corrosion is a major economic concern in modern society. A phenomenon that has puzzled generations of scientists in this field is the so-called anomalous hydrogen evolution: the violent dissolution of magnesium under electron-rich (anodic) conditions, accompanied by strong hydrogen evolution, and a key mechanism hampering Mg technology. Experimental studies have indicated the presenc…
▽ More
Aqueous metal corrosion is a major economic concern in modern society. A phenomenon that has puzzled generations of scientists in this field is the so-called anomalous hydrogen evolution: the violent dissolution of magnesium under electron-rich (anodic) conditions, accompanied by strong hydrogen evolution, and a key mechanism hampering Mg technology. Experimental studies have indicated the presence of univalent Mg$^+$ in solution, but these findings have been largely ignored because they defy our common chemical understanding and evaded direct experimental observation. Using recent advances in the \emph{ab initio} description of solid-liquid electrochemical interfaces under controlled potential conditions, we described the full reaction path of Mg atom dissolution from a kinked Mg surface under anodic conditions. Our study reveals the formation of a solvated [Mg$^{2+}$(OH)$^-$]$^+$ ion complex, challenging the conventional assumption of Mg$^{2+}$ ion. This insight provides an intuitive explanation for the postulated presence of (coulombically) univalent Mg$^+$ ions and the absence of protective oxide/hydroxide layers normally formed under anodic/oxidizing conditions. The discovery of this unexpected and unconventional reaction mechanism is crucial for identifying new strategies for corrosion prevention and can be transferred to other metals.
△ Less
Submitted 18 July, 2024;
originally announced July 2024.
-
From electrons to phase diagrams with classical and machine learning potentials: automated workflows for materials science with pyiron
Authors:
Sarath Menon,
Yury Lysogorskiy,
Alexander L. M. Knoll,
Niklas Leimeroth,
Marvin Poul,
Minaam Qamar,
Jan Janssen,
Matous Mrovec,
Jochen Rohrer,
Karsten Albe,
Jörg Behler,
Ralf Drautz,
Jörg Neugebauer
Abstract:
We present a comprehensive and user-friendly framework built upon the pyiron integrated development environment (IDE), enabling researchers to perform the entire Machine Learning Potential (MLP) development cycle consisting of (i) creating systematic DFT databases, (ii) fitting the Density Functional Theory (DFT) data to empirical potentials or MLPs, and (iii) validating the potentials in a largel…
▽ More
We present a comprehensive and user-friendly framework built upon the pyiron integrated development environment (IDE), enabling researchers to perform the entire Machine Learning Potential (MLP) development cycle consisting of (i) creating systematic DFT databases, (ii) fitting the Density Functional Theory (DFT) data to empirical potentials or MLPs, and (iii) validating the potentials in a largely automatic approach. The power and performance of this framework are demonstrated for three conceptually very different classes of interatomic potentials: an empirical potential (embedded atom method - EAM), neural networks (high-dimensional neural network potentials - HDNNP) and expansions in basis sets (atomic cluster expansion - ACE). As an advanced example for validation and application, we show the computation of a binary composition-temperature phase diagram for Al-Li, a technologically important lightweight alloy system with applications in the aerospace industry.
△ Less
Submitted 8 March, 2024;
originally announced March 2024.
-
Roadmap on Data-Centric Materials Science
Authors:
Stefan Bauer,
Peter Benner,
Tristan Bereau,
Volker Blum,
Mario Boley,
Christian Carbogno,
C. Richard A. Catlow,
Gerhard Dehm,
Sebastian Eibl,
Ralph Ernstorfer,
Ádám Fekete,
Lucas Foppa,
Peter Fratzl,
Christoph Freysoldt,
Baptiste Gault,
Luca M. Ghiringhelli,
Sajal K. Giri,
Anton Gladyshev,
Pawan Goyal,
Jason Hattrick-Simpers,
Lara Kabalan,
Petr Karpov,
Mohammad S. Khorrami,
Christoph Koch,
Sebastian Kokott
, et al. (36 additional authors not shown)
Abstract:
Science is and always has been based on data, but the terms "data-centric" and the "4th paradigm of" materials research indicate a radical change in how information is retrieved, handled and research is performed. It signifies a transformative shift towards managing vast data collections, digital repositories, and innovative data analytics methods. The integration of Artificial Intelligence (AI) a…
▽ More
Science is and always has been based on data, but the terms "data-centric" and the "4th paradigm of" materials research indicate a radical change in how information is retrieved, handled and research is performed. It signifies a transformative shift towards managing vast data collections, digital repositories, and innovative data analytics methods. The integration of Artificial Intelligence (AI) and its subset Machine Learning (ML), has become pivotal in addressing all these challenges. This Roadmap on Data-Centric Materials Science explores fundamental concepts and methodologies, illustrating diverse applications in electronic-structure theory, soft matter theory, microstructure research, and experimental techniques like photoemission, atom probe tomography, and electron microscopy. While the roadmap delves into specific areas within the broad interdisciplinary field of materials science, the provided examples elucidate key concepts applicable to a wider range of topics. The discussed instances offer insights into addressing the multifaceted challenges encountered in contemporary materials research.
△ Less
Submitted 1 May, 2024; v1 submitted 1 February, 2024;
originally announced February 2024.
-
Understanding atom probe's analytical performance for iron oxides using correlation histograms and ab initio calculations
Authors:
Se-Ho Kim,
Shalini Bhatt,
Daniel K. Schreiber,
Jörg Neugebauer,
Christoph Freysoldt,
Baptiste Gault,
Shyam Katnagallu
Abstract:
Field evaporation from ionic or covalently bonded materials often leads to the emission of molecular ions. The metastability of these molecular ions, particularly under the influence of the intense electrostatic field (1010 Vm-1), makes them prone to dissociation with or without an exchange of energy amongst them. These processes can affect the analytical performance of atom probe tomography (APT)…
▽ More
Field evaporation from ionic or covalently bonded materials often leads to the emission of molecular ions. The metastability of these molecular ions, particularly under the influence of the intense electrostatic field (1010 Vm-1), makes them prone to dissociation with or without an exchange of energy amongst them. These processes can affect the analytical performance of atom probe tomography (APT). For instance, neutral species formed through dissociation may not be detected at all or with a time of flight no longer related to their mass, causing their loss from the analysis. Here, we evaluated the changes in the measured composition of FeO, Fe2O3 and Fe3O4 across a wide range of analysis conditions. Possible dissociation reactions are predicted by density-functional theory (DFT) calculations considering the spin states of the molecules. The energetically favoured reactions are traced on to the multi-hit ion correlation histograms, to confirm their existence within experiments, using an automated Python-based routine. The detected reactions are carefully analysed to reflect upon the influence of these neutrals from dissociation reactions on the performance of APT for analysing iron oxides.
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
Triplet Excitation-Energy Transfer Couplings from Subsystem Time-Dependent Density-Functional Theory
Authors:
Sabine Käfer,
Niklas Niemeyer,
Johannes Tölle,
Johannes Neugebauer
Abstract:
We present an implementation of Triplet Excitation-Energy Transfer (TEET) couplings based on subsystem-based Time-Dependent Density-Functional Theory (sTDDFT). TEET couplings are systematically investigated by comparing "exact" and approximate variants of sTDDFT. We demonstrate that, while sTDDFT utilizing explicit approximate Non-Additive Kinetic Energy (NAKE) density functionals is well-suited f…
▽ More
We present an implementation of Triplet Excitation-Energy Transfer (TEET) couplings based on subsystem-based Time-Dependent Density-Functional Theory (sTDDFT). TEET couplings are systematically investigated by comparing "exact" and approximate variants of sTDDFT. We demonstrate that, while sTDDFT utilizing explicit approximate Non-Additive Kinetic Energy (NAKE) density functionals is well-suited for describing Singlet Excitation-Energy Transfer (SEET) processes, it is inadequate for characterizing TEET. However, we show that Projection-based Embedding (PbE)-based sTDDFT addresses the challenges faced by NAKE-sTDDFT and emerges as a promising method for accurately describing electronic couplings in TEET processes.
We also introduce the mixed PbE-/NAKE-embedding procedure to investigate TEET effects in solvated pairs of chromophores. This approach offers a good balance between accuracy and efficiency, enabling comprehensive studies of TEET processes in complex environments.
△ Less
Submitted 12 December, 2023;
originally announced December 2023.
-
Accelerating Analytic-Continuation GW Calculations with a Laplace Transformation and Natural Auxiliary Functions
Authors:
Johannes Tölle,
Niklas Niemeyer,
Johannes Neugebauer
Abstract:
We present a simple and accurate GW implementation based on a combination of a Laplace transformation (LT) and other acceleration techniques used in post-SCF quantum chemistry, namely, natural auxiliary functions and the frozen-core approximation. The LT-GW approach combines three major benefits: (a) a small prefactor for the computational scaling, (b) easy integration into existing molecular GW i…
▽ More
We present a simple and accurate GW implementation based on a combination of a Laplace transformation (LT) and other acceleration techniques used in post-SCF quantum chemistry, namely, natural auxiliary functions and the frozen-core approximation. The LT-GW approach combines three major benefits: (a) a small prefactor for the computational scaling, (b) easy integration into existing molecular GW implementations, and (c) significant performance improvements for a wide range of possible applications. Illustrating these advantages for systems consisting of up to 352 atoms and 7412 basis functions, we further demonstrate the benefits of this approach combined with an efficient implementation of the Bethe-Salpeter equation.
△ Less
Submitted 15 November, 2023; v1 submitted 10 July, 2023;
originally announced July 2023.
-
A machine learning framework for quantifying chemical segregation and microstructural features in atom probe tomography data
Authors:
Alaukik Saxena,
Nikita Polin,
Navyanth Kusampudi,
Shyam Katnagallu,
Leopoldo Molina-Luna,
Oliver Gutfleisch,
Benjamin Berkels,
Baptiste Gault,
Jörg Neugebauer,
Christoph Freysoldt
Abstract:
Atom probe tomography (APT) is ideally suited to characterize and understand the interplay of chemical segregation and microstructure in modern multicomponent materials. Yet, the quantitative analysis typically relies on human expertise to define regions of interest. We introduce a computationally efficient, multistage machine learning strategy to identify chemically distinct domains in a semi aut…
▽ More
Atom probe tomography (APT) is ideally suited to characterize and understand the interplay of chemical segregation and microstructure in modern multicomponent materials. Yet, the quantitative analysis typically relies on human expertise to define regions of interest. We introduce a computationally efficient, multistage machine learning strategy to identify chemically distinct domains in a semi automated way, and subsequently quantify their geometric and compositional characteristics. In our algorithmic pipeline, we first coarse grain the APT data into voxels, collect the composition statistics, and decompose it via clustering in composition space. The composition classification then enables the real space segmentation via a density based clustering algorithm, thus revealing the microstructure at voxel resolution. Our approach is demonstrated for a Sm(Co,Fe)ZrCu alloy. The alloy exhibits two precipitate phases with a plate-like, but intertwined morphology. The primary segmentation is further refined to disentangle these geometrically complex precipitates into individual plate like parts by an unsupervised approach based on principle component analysis, or a U-Net based semantic segmentation trained on the former. Following the chemical and geometric analysis, detailed chemical distribution and segregation effects relative to the predominant plate-like geometry can be readily mapped without resorting to the initial voxelization.
△ Less
Submitted 18 April, 2023;
originally announced April 2023.
-
Machine learning-enabled tomographic imaging of chemical short-range atomic ordering
Authors:
Yue Li,
Timoteo Colnaghi,
Yilun Gong,
Huaide Zhang,
Yuan Yu,
Ye Wei,
Bin Gan,
Min Song,
Andreas Marek,
Markus Rampp,
Siyuan Zhang,
Zongrui Pei,
Matthias Wuttig,
Jörg Neugebauer,
Zhangwei Wang,
Baptiste Gault
Abstract:
In solids, chemical short-range order (CSRO) refers to the self-organisation of atoms of certain species occupying specific crystal sites. CSRO is increasingly being envisaged as a lever to tailor the mechanical and functional properties of materials. Yet quantitative relationships between properties and the morphology, number density, and atomic configurations of CSRO domains remain elusive. Here…
▽ More
In solids, chemical short-range order (CSRO) refers to the self-organisation of atoms of certain species occupying specific crystal sites. CSRO is increasingly being envisaged as a lever to tailor the mechanical and functional properties of materials. Yet quantitative relationships between properties and the morphology, number density, and atomic configurations of CSRO domains remain elusive. Herein, we showcase how machine learning-enhanced atom probe tomography (APT) can mine the near-atomically resolved APT data and jointly exploit the technique's high elemental sensitivity to provide a 3D quantitative analysis of CSRO in a CoCrNi medium-entropy alloy. We reveal multiple CSRO configurations, with their formation supported by state-of-the-art Monte-Carlo simulations. Quantitative analysis of these CSROs allows us to establish relationships between processing parameters and physical properties. The unambiguous characterization of CSRO will help refine strategies for designing advanced materials by manipulating atomic-scale architectures.
△ Less
Submitted 23 March, 2023;
originally announced March 2023.
-
Solute Induced Defect Phase Transformations in Mg Grain Boundaries
Authors:
P. Mathews,
S. Zhang,
C. Scheu,
R. Janisch,
J. Neugebauer,
T. Hickel
Abstract:
The study of defect phases is important for designing nanostructured metals and alloys. Grain boundaries (GBs) form one class of defects that directly influence materials properties, such as deformability and strength. At the same time, alloying can introduce GB phase transformations and can therewith alter mechanical performance. In this work, the defect phases of a $Σ$7 $(12\bar{3}0)$ [0001] 21.…
▽ More
The study of defect phases is important for designing nanostructured metals and alloys. Grain boundaries (GBs) form one class of defects that directly influence materials properties, such as deformability and strength. At the same time, alloying can introduce GB phase transformations and can therewith alter mechanical performance. In this work, the defect phases of a $Σ$7 $(12\bar{3}0)$ [0001] 21.78$^\circ$ symmetric tilt GB in hcp Mg are investigated. Ab-initio simulations as a function of stress and temperature (using quasi-harmonic approximation) are performed, and different types of phase transformations are revealed. To this end, the influence of the chemical degree of freedom on the defect phases is studied for the example of Ga addition, using an efficient screening approach that combines empirical potentials and accurate ab-initio calculations. By exploiting the concept of defect phase diagrams, a phase transformation from the T to the A structural type and as well as a systematic transition of the segregation site preference is revealed. The results qualitatively agree well with experimental observations from scanning transition electron microscopy. The underlying physical mechanisms have an impact on grain-boundary engineering in metallic alloys.
△ Less
Submitted 16 March, 2023;
originally announced March 2023.
-
Phase stability and defect studies of Mg-based Laves phases using defect phase diagrams
Authors:
A. Tehranchi,
M. Lipinska-Chwalek,
J. Mayer,
J. Neugebauer,
T. Hickel
Abstract:
Laves phases often form as secondary phases in metallic alloys and have a significant effect on their structural properties. Thus, phase stability studies for these chemically and structurally complex phases in addition to mechanical behavior studies are of great interest. In this work, we use the concept of metastable bulk phase and defect phase diagrams to augment the understanding of the bulk p…
▽ More
Laves phases often form as secondary phases in metallic alloys and have a significant effect on their structural properties. Thus, phase stability studies for these chemically and structurally complex phases in addition to mechanical behavior studies are of great interest. In this work, we use the concept of metastable bulk phase and defect phase diagrams to augment the understanding of the bulk phase and defect phase stability in Laves phases in Mg-based alloys. In this way, we resolve the discrepancy between bulk phase diagrams and experimental observations regarding the formation of Mg-rich C14 and Al-rich C15 Laves phases in MgAlCa alloys at moderate temperatures. Moreover, the effect of the thermodynamic state of alloys on the competition between solute-rich hcp-like planar defects and stoichiometric basal stacking faults is clarified, which determines the brittleness of these alloys. \end{abstract}
△ Less
Submitted 16 March, 2023;
originally announced March 2023.
-
Constructing phase diagrams for defects by correlated atomic-scale characterization
Authors:
Xuyang Zhou,
Prince Mathews,
Benjamin Berkels,
Saba Ahmad,
Amel Shamseldeen Ali Alhassan,
Philipp Keuter,
Jochen M. Schneider,
Dierk Raabe,
Jörg Neugebauer,
Gerhard Dehm,
Tilmann Hickel,
Christina Scheu,
Siyuan Zhang
Abstract:
Phase transformations and crystallographic defects are two essential tools to drive innovations in materials. Bulk materials design via tuning chemical compositions has been systematized using phase diagrams. We show here that the same thermodynamic concept can be applied to understand the chemistry at defects. We present a combined experimental and modelling approach to scope and build phase diag…
▽ More
Phase transformations and crystallographic defects are two essential tools to drive innovations in materials. Bulk materials design via tuning chemical compositions has been systematized using phase diagrams. We show here that the same thermodynamic concept can be applied to understand the chemistry at defects. We present a combined experimental and modelling approach to scope and build phase diagrams for defects. The discovery was enabled by triggering phase transformations of individual defects through local alloying, and sequentially imaging the structural and chemical changes using atomic-resolution scanning transmission electron microscopy. By observing atomic-scale phase transformations of a Mg grain boundary through Ga alloying, we exemplified the method to construct a grain boundary phase diagram using ab initio simulations and thermodynamic principles. The methodology enables a systematic development of defect phase diagrams to propel a new paradigm for materials design utilizing chemical complexity and phase transformations at defects.
△ Less
Submitted 31 August, 2023; v1 submitted 16 March, 2023;
originally announced March 2023.
-
Metastable defect phase diagrams as a tool to describe chemically driven defect formation: Application to planar defects
Authors:
A. Tehranchi,
S. Zhang,
A. Zendegani,
C. Scheu,
T. Hickel,
J. Neugebauer
Abstract:
Thermodynamic bulk phase diagrams have become the roadmap used by researchers to identify alloy compositions and process conditions that result in novel materials with tailored microstructures. Recent experimental studies show that changes in the alloy composition can drive not only transitions in the bulk phases present in a material, but also in the concentration and type of defects they contain…
▽ More
Thermodynamic bulk phase diagrams have become the roadmap used by researchers to identify alloy compositions and process conditions that result in novel materials with tailored microstructures. Recent experimental studies show that changes in the alloy composition can drive not only transitions in the bulk phases present in a material, but also in the concentration and type of defects they contain. Defect phase diagrams in combination with density functional theory provide a natural route to study these chemically driven defects. Our results show, however, that direct application of thermodynamic approaches can fail to reproduce the experimentally observed defect formation. Therefore, we extend the concept to metastable defect phase diagrams to account for kinetic limitations that prevent the system from reaching equilibrium. We successfully applied this concept to explain the formation of large concentrations of planar defects in supersaturated Fe-Nb solid solutions and to identify in a joint study with experiments conditions in Mg-Al-Ca alloys for defect phase occurrence. The concept offers new avenues for designing materials with tailored defect structures.
△ Less
Submitted 13 March, 2023;
originally announced March 2023.
-
Quantification of electronic and magnetoelastic mechanisms of first-order magnetic phase transitions from first principles: application to caloric effects in La(Fe$_x$Si$_{1-x}$)$_{13}$
Authors:
Eduardo Mendive Tapia,
Christopher E. Patrick,
Tilmann Hickel,
Jörg Neugebauer,
Julie B. Staunton
Abstract:
La(Fe$_x$Si$_{1-x}$)$_{13}$ and derived quaternary compounds are well-known for their giant, tunable, magneto- and barocaloric responses around a first-order paramagnetic-ferromagnetic transition near room temperature with low hysteresis. Remarkably, such a transition shows a large spontaneous volume change together with itinerant electron metamagnetic features. While magnetovolume effects are wel…
▽ More
La(Fe$_x$Si$_{1-x}$)$_{13}$ and derived quaternary compounds are well-known for their giant, tunable, magneto- and barocaloric responses around a first-order paramagnetic-ferromagnetic transition near room temperature with low hysteresis. Remarkably, such a transition shows a large spontaneous volume change together with itinerant electron metamagnetic features. While magnetovolume effects are well-established mechanisms driving first-order transitions, purely electronic sources have a long, subtle history and remain poorly understood. Here we apply a disordered local moment picture to quantify electronic and magnetoelastic effects at finite temperature in La(Fe$_x$Si$_{1-x}$)$_{13}$ from first-principles. We obtain results in very good agreement with experiment and demonstrate that the magnetoelastic coupling, rather than purely electronic mechanisms, drives the first-order character and causes at the same time a huge electronic entropy contribution to the caloric response.
△ Less
Submitted 24 October, 2023; v1 submitted 13 February, 2023;
originally announced February 2023.
-
Accurate computation of chemical contrast in field ion microscopy
Authors:
Shalini Bhatt,
Shyam Katnagallu,
Jörg Neugebauer,
Christoph Freysoldt
Abstract:
We present a computational approach to simulate local contrast observed in Field Ion Microscopy (FIM). It is based on density-functional theory utilizing the Tersoff-Hamann approach as done in Scanning Tunneling Microscopy (STM). A key requirement is the highly accurate computation of the surface states' wave-function tails. To refine the Kohn-Sham states from standard iterative global solvers we…
▽ More
We present a computational approach to simulate local contrast observed in Field Ion Microscopy (FIM). It is based on density-functional theory utilizing the Tersoff-Hamann approach as done in Scanning Tunneling Microscopy (STM). A key requirement is the highly accurate computation of the surface states' wave-function tails. To refine the Kohn-Sham states from standard iterative global solvers we introduce and discuss the EXtrapolation of Tails via Reverse integration Algorithm (EXTRA). The decaying tails are obtained by reverse integration (from outside in) using a Numerov-like algorithm. The starting conditions are then iteratively adapted to match the values of plane-wave Kohn-Sham wave functions close to the surface. We demonstrate the performance of the proposed algorithm by analysing and showing the chemical contrast for Ta at Ni surface.
△ Less
Submitted 10 May, 2023; v1 submitted 21 December, 2022;
originally announced December 2022.
-
Computationally accelerated experimental materials characterization -- drawing inspiration from high-throughput simulation workflows
Authors:
Markus Stricker,
Lars Banko,
Nik Sarazin,
Niklas Siemer,
Jörg Neugebauer,
Alfred Ludwig
Abstract:
Computational materials science is increasingly benefitting from data management, automation, and algorithm-based decision-making in controlling simulations. Experimental materials science is also undergoing a change and increasingly more `machine learning' is incorporated in materials discovery campaigns. The obvious benefits include automation, reproducibility, data provenance, and reusability o…
▽ More
Computational materials science is increasingly benefitting from data management, automation, and algorithm-based decision-making in controlling simulations. Experimental materials science is also undergoing a change and increasingly more `machine learning' is incorporated in materials discovery campaigns. The obvious benefits include automation, reproducibility, data provenance, and reusability of managed data, however, is not widely available. We demonstrate an implementation of a Gaussian Process Regression directly controlling an experimental measurement device in pyiron, a framework designed for high-throughput simulations, as a first step to increasingly combine experimental and simulated data in one framework. With data from both in the same framework, a heretofore untapped and much-needed potential for the acceleration of materials characterization and materials discovery campaigns becomes available.
△ Less
Submitted 9 December, 2022;
originally announced December 2022.
-
Tailoring negative pressure by crystal defects: Crack induced hydride formation in Al alloys
Authors:
A. Tehranchi,
P. Chakraborty,
M. L. Freixes,
E. McEniry,
B. Gault,
T. Hickel,
J. Neugebauer
Abstract:
Climate change motivates the search for non-carbon-emitting energy generation and storage solutions. Metal hydrides show promising characteristics for this purpose. They can be further stabilized by tailoring the negative pressure of microstructural and structural defects. Using systematic ab initio and atomistic simulations, we demonstrate that an enhancement in the formation of hydrides at the n…
▽ More
Climate change motivates the search for non-carbon-emitting energy generation and storage solutions. Metal hydrides show promising characteristics for this purpose. They can be further stabilized by tailoring the negative pressure of microstructural and structural defects. Using systematic ab initio and atomistic simulations, we demonstrate that an enhancement in the formation of hydrides at the negatively pressurized crack tip region is feasible by increasing the mechanical tensile load on the specimen. The theoretical predictions have been used to reassess and interpret atom probe tomography experiments for a high-strength 7XXX-aluminium alloy that show a substantial enhancement of hydrogen concentration at structural defects near a stress-corrosion crack tip. These results contain important implications for enhancing the capability of metals as H-storage materials.
△ Less
Submitted 20 September, 2022;
originally announced September 2022.
-
Effect of Sn on generalized stacking fault energy surfaces in zirconium and its hydrides
Authors:
P. Chakraborty,
I. Mouton,
B. Gault,
A. Tehranchi,
J. Neugebauer,
T. Hickel
Abstract:
Hydrogen embrittlement in Zr alloy fuel cladding is a primary safety concern for water based nuclear reactors. Here we investigated the stabilisation of planar defects within the forming hydrides by Sn, the primary alloying element of Zircaloy-4 used in the cladding. In order to explain formation of hydrides and planar defects observed in our experiments, we performed atomic-scale ab initio calcul…
▽ More
Hydrogen embrittlement in Zr alloy fuel cladding is a primary safety concern for water based nuclear reactors. Here we investigated the stabilisation of planar defects within the forming hydrides by Sn, the primary alloying element of Zircaloy-4 used in the cladding. In order to explain formation of hydrides and planar defects observed in our experiments, we performed atomic-scale ab initio calculations focusing on the solute interactions with generalized stacking faults in hcp $α$-Zr and fcc zirconium hydrides. Our calculations showed that an increase in Sn concentration leads to a stabilisation of stacking faults in both $α$-Zr and hydride phases. However, the solution enthalpy of Sn is lower in the $α$-Zr as compared to the other hydride phases indicative of two competing processes of Sn depletion/enrichment at the Zr hydride/matrix interface. This is corroborated by experimental findings, where Sn is repelled by hydrides and is mostly found trapped at interfaces and planar defects indicative of stacking faults inside the hydride phases. Our systematic investigation enables us to understand the presence and distribution of solutes in the hydride phases, which provides a deeper insight into the microstructural evolution of such alloy's properties during its service lifetime.
△ Less
Submitted 19 September, 2022;
originally announced September 2022.
-
Ab initio vacancy formation energies and kinetics at metal surfaces under high electric field
Authors:
Shyam Katnagallu,
Christoph Freysoldt,
Baptiste Gault,
Jörg Neugebauer
Abstract:
Recording field ion microscope images under field evaporating conditions and subsequently reconstructing the underlying atomic configuration, called three-dimensional field ion microscopy (3D-FIM) is one of the few techniques capable of resolving crystalline defects at an atomic scale. However, the quantification of the observed vacancies and their origins are still a matter of debate. It was sugg…
▽ More
Recording field ion microscope images under field evaporating conditions and subsequently reconstructing the underlying atomic configuration, called three-dimensional field ion microscopy (3D-FIM) is one of the few techniques capable of resolving crystalline defects at an atomic scale. However, the quantification of the observed vacancies and their origins are still a matter of debate. It was suggested that high electric fields (1-5 V/Å) used in 3D-FIM could introduce artefact vacancies. To investigate such effects, we used density functional theory (DFT) simulations. Stepped Ni and Pt surfaces with kinks were modelled in the repeated slab approach with a (971) surface orientation. An electrostatic field of up to 4 V/Å was introduced on one side of the slab using the generalized dipole correction. Contrary to what was proposed, we show that the formation of vacancies on the electrified metal surface is more difficult compared to a field-free case. We also find that the electric field can introduce kinetic barriers to a potential vacancy-annihilation mechanism. We rationalize these findings by comparing to insights from field evaporation models.
△ Less
Submitted 14 September, 2022;
originally announced September 2022.
-
Construction and analysis of surface phase diagrams to describe segregation and dissolution behavior of Al and Ca in Mg alloys
Authors:
Jing Yang,
K. B. Sravan Kumar,
Mira Todorova,
Jörg Neugebauer
Abstract:
Segregation and dissolution behavior of Mg alloyed with Ca and Al are studied by performing density functional theory calculations considering an extensive set of surface structures and compositions. Combining ab initio surface science approaches with cluster expansion for ordered surface structures we construct surface phase diagrams for these alloys. We utilize these diagrams to study segregatio…
▽ More
Segregation and dissolution behavior of Mg alloyed with Ca and Al are studied by performing density functional theory calculations considering an extensive set of surface structures and compositions. Combining ab initio surface science approaches with cluster expansion for ordered surface structures we construct surface phase diagrams for these alloys. We utilize these diagrams to study segregation phenomena and chemical trends for surfaces in contact with a dry environment or with an aqueous electrolyte. We show that the presence of water dramatically impacts the stability and chemical composition of the considered metallic surfaces. We furthermore find that the two alloying elements behave qualitatively different: whereas Ca strongly segregates to the surface and becomes dissolved upon exposure of the surface to water, Al shows an anti-segregation behavior, i.e., it remains in Mg bulk. These findings provide an explanation for the experimentally observed increase/decrease in corrosion rates when alloying Mg with Al/Ca.
△ Less
Submitted 10 July, 2023; v1 submitted 20 July, 2022;
originally announced July 2022.
-
Systematic Atomic Structure Datasets for Machine Learning Potentials: Application to Defects in Magnesium
Authors:
Marvin Poul,
Liam Huber,
Erik Bitzek,
Jörg Neugebauer
Abstract:
We present a physically motivated strategy for the construction of training sets for transferable machine learning interatomic potentials. It is based on a systematic exploration of all possible space groups in random crystal structures, together with deformations of cell shape, size, and atomic positions. The resulting potentials turn out to be unbiased and generically applicable to studies of bu…
▽ More
We present a physically motivated strategy for the construction of training sets for transferable machine learning interatomic potentials. It is based on a systematic exploration of all possible space groups in random crystal structures, together with deformations of cell shape, size, and atomic positions. The resulting potentials turn out to be unbiased and generically applicable to studies of bulk defects without including any defect structures in the training set or employing any additional Active Learning. Using this approach we construct transferable potentials for pure Magnesium that reproduce the properties of hexagonal closed packed (hcp) and body centered cubic (bcc) polymorphs very well. In the process we investigate how different types of training structures impact the properties and the predictive power of the resulting potential.
△ Less
Submitted 9 February, 2023; v1 submitted 8 July, 2022;
originally announced July 2022.
-
Lyapunov-Type Inequalities for a Fractional Boundary Value Problem with a Fractional Boundary Condition
Authors:
Sougata Dhar,
Jeffrey T. Neugebauer
Abstract:
In this paper, we consider a linear fractional differential equation with fractional boundary conditions. First, by obtaining Green's function, we derive the Lyapunov-type inequalities for such boundary value problems. Furthermore, we use the contraction mapping theorem to study the existence of a unique solution to a nonlinear problem.
In this paper, we consider a linear fractional differential equation with fractional boundary conditions. First, by obtaining Green's function, we derive the Lyapunov-type inequalities for such boundary value problems. Furthermore, we use the contraction mapping theorem to study the existence of a unique solution to a nonlinear problem.
△ Less
Submitted 4 May, 2022;
originally announced May 2022.
-
Existence of positive solutions of a Hammerstein integral equation using the layered compression-expansion fixed point theorem
Authors:
Sougata Dhar,
Jeffrey W. Lyons,
Jeffrey T. Neugebauer
Abstract:
In this paper, we show the existence of a positive solution of a Hammerstein integral equation under certain conditions on the kernel. We apply the recent Layered Compression-Expansion Fixed Point Theorem. Finally, we provide corollaries to help in application of the main results and present an example.
In this paper, we show the existence of a positive solution of a Hammerstein integral equation under certain conditions on the kernel. We apply the recent Layered Compression-Expansion Fixed Point Theorem. Finally, we provide corollaries to help in application of the main results and present an example.
△ Less
Submitted 4 May, 2022;
originally announced May 2022.
-
Controlled doping of electrocatalysts through engineering impurities
Authors:
Se-Ho Kim,
Su-Hyun Yoo,
Sangyong Shin,
Ayman A. El-Zoka,
Olga Kasian,
Joohyun Lim,
Jiwon Jeong,
Christina Scheu,
Joerg Neugebauer,
Hyunjoo Lee,
Mira Todorova,
Baptiste Gault
Abstract:
Fuel cells recombine water from H2 and O2 thereby powering e.g. cars or houses with no direct carbon emission. In anion-exchange membrane fuel cells (AEMFCs), to reach high power densities, operating at high pH is an alternative to using large volumes of noble metals catalysts at the cathode, where the oxygen-reduction reaction occurs. However, the sluggish kinetics of the hydrogen-oxidation react…
▽ More
Fuel cells recombine water from H2 and O2 thereby powering e.g. cars or houses with no direct carbon emission. In anion-exchange membrane fuel cells (AEMFCs), to reach high power densities, operating at high pH is an alternative to using large volumes of noble metals catalysts at the cathode, where the oxygen-reduction reaction occurs. However, the sluggish kinetics of the hydrogen-oxidation reaction (HOR) hinders upscaling despite promising catalysts. Here, we observe an unexpected ingress of B into Pd nano-catalysts synthesised by wet-chemistry, gain control over this B-doping, and report on its influence on the HOR activity in alkaline conditions. We rationalize our findings using ab-initio calculations of both H- and OH-adsorption on B-doped Pd. Using this "impurity engineering" approach, we thus design Pt-free catalysts as required in electrochemical energy conversion devices, e.g. next generations of AEMFCs, that satisfy the economic and environmental constraints, i.e. reasonable operating costs and long-term stability, to enable the hydrogen economy.
△ Less
Submitted 1 March, 2022;
originally announced March 2022.
-
Machine learning-enabled high-entropy alloy discovery
Authors:
Ziyuan Rao,
PoYen Tung,
Ruiwen Xie,
Ye Wei,
Hongbin Zhang,
Alberto Ferrari,
T. P. C. Klaver,
Fritz Körmann,
Prithiv Thoudden Sukumar,
Alisson Kwiatkowski da Silva,
Yao Chen,
Zhiming Li,
Dirk Ponge,
Jörg Neugebauer,
Oliver Gutfleisch,
Stefan Bauer,
Dierk Raabe
Abstract:
High-entropy alloys are solid solutions of multiple principal elements, capable of reaching composition and feature regimes inaccessible for dilute materials. Discovering those with valuable properties, however, relies on serendipity, as thermodynamic alloy design rules alone often fail in high-dimensional composition spaces. Here, we propose an active-learning strategy to accelerate the design of…
▽ More
High-entropy alloys are solid solutions of multiple principal elements, capable of reaching composition and feature regimes inaccessible for dilute materials. Discovering those with valuable properties, however, relies on serendipity, as thermodynamic alloy design rules alone often fail in high-dimensional composition spaces. Here, we propose an active-learning strategy to accelerate the design of novel high-entropy Invar alloys in a practically infinite compositional space, based on very sparse data. Our approach works as a closed-loop, integrating machine learning with density-functional theory, thermodynamic calculations, and experiments. After processing and characterizing 17 new alloys (out of millions of possible compositions), we identified 2 high-entropy Invar alloys with extremely low thermal expansion coefficients around 2*10-6 K-1 at 300 K. Our study thus opens a new pathway for the fast and automated discovery of high-entropy alloys with optimal thermal, magnetic and electrical properties.
△ Less
Submitted 28 February, 2022;
originally announced February 2022.
-
Ab initio calculation of the magnetic Gibbs free energy of materials using magnetically constrained supercells
Authors:
Eduardo Mendive-Tapia,
Jörg Neugebauer,
Tilmann Hickel
Abstract:
We present a first-principles approach for the computation of the magnetic Gibbs free energy of materials using magnetically constrained supercell calculations. Our approach is based on an adiabatic approximation of slowly varying local moment orientations, the so-called finite-temperature disordered local moment picture. It describes magnetic phase transitions and how electronic and/or magnetostr…
▽ More
We present a first-principles approach for the computation of the magnetic Gibbs free energy of materials using magnetically constrained supercell calculations. Our approach is based on an adiabatic approximation of slowly varying local moment orientations, the so-called finite-temperature disordered local moment picture. It describes magnetic phase transitions and how electronic and/or magnetostructural mechanisms generate a discontinuous (first-order) character. We demonstrate that the statistical mechanics of the local moment orientations can be described by an affordable number of supercell calculations containing noncollinear magnetic configurations. The applicability of our approach is illustrated by firstly studying the ferromagnetic state in bcc Fe. We then investigate the temperature-dependent properties of a triangular antiferromagnetic state stabilizing in two antiperovskite systems Mn$_3$AN (A = Ga, Ni). Our calculations provide the negative volume expansion of these materials as well as the ab initio origin of the discontinuous character of the phase transitions, electronic and/or magnetostructural, in good agreement with experiment.
△ Less
Submitted 23 February, 2022;
originally announced February 2022.
-
Photoinduced structural dynamics of multiferroic TbMnO$_3$
Authors:
Elsa Abreu,
Matteo Savoini,
Larissa Boie,
Paul Beaud,
Vincent Esposito,
Martin Kubli,
Martin J. Neugebauer,
Michael Porer,
Urs Staub,
Bulat Burganov,
Chris Dornes,
Angel Rodriguez-Fernandez,
Lucas Huber,
Gabriel Lantz,
José R. L. Mardegan,
Sergii Parchenko,
Jochen Rittmann,
Cris Svetina,
Gerhard Ingold,
Steven L. Johnson
Abstract:
We use time-resolved hard x-ray diffraction to investigate the structural dynamics of the multiferroic insulator TbMnO$_3$ in the low temperature antiferromagnetic and ferroelectrically ordered phase. The lattice response following photoexcitation at 1.55 eV is detected by measuring the (0 2 4) and (1 3 -5) Bragg reflections. A 0.02% tensile strain, normal to the surface, is seen to arise within 2…
▽ More
We use time-resolved hard x-ray diffraction to investigate the structural dynamics of the multiferroic insulator TbMnO$_3$ in the low temperature antiferromagnetic and ferroelectrically ordered phase. The lattice response following photoexcitation at 1.55 eV is detected by measuring the (0 2 4) and (1 3 -5) Bragg reflections. A 0.02% tensile strain, normal to the surface, is seen to arise within 20 - 30 ps. The magnitude of this transient strain is over an order of magnitude lower than that predicted from laser-induced heating, which we attribute to a bottleneck in the energy transfer between the electronic and lattice subsystems. The timescale for the transient expansion is consistent with that of previously reported demagnetization dynamics. We discuss a possible relationship between structural and demagnetization dynamics in TbMnO$_3$, in which photoinduced atomic motion modulates the exchange interaction, leading to a destruction of the magnetic order in the system.
△ Less
Submitted 24 January, 2022;
originally announced January 2022.
-
Automated optimization of convergence parameters in plane wave density functional theory calculations via a tensor decomposition-based uncertainty quantification
Authors:
Jan Janssen,
Edgar Makarov,
Tilmann Hickel,
Alexander V. Shapeev,
Jörg Neugebauer
Abstract:
First principles approaches have revolutionized our ability in using computers to predict, explore and design materials. A major advantage commonly associated with these approaches is that they are fully parameter free. However, numerically solving the underlying equations requires to choose a set of convergence parameters. With the advent of high-throughput calculations it becomes exceedingly imp…
▽ More
First principles approaches have revolutionized our ability in using computers to predict, explore and design materials. A major advantage commonly associated with these approaches is that they are fully parameter free. However, numerically solving the underlying equations requires to choose a set of convergence parameters. With the advent of high-throughput calculations it becomes exceedingly important to achieve a truly parameter free approach. Utilizing uncertainty quantification (UQ) and tensor decomposition we derive a numerically highly efficient representation of the statistical and systematic error in the multidimensional space of the convergence parameters. Based on this formalism we implement a fully automated approach that requires as input the target accuracy rather than convergence parameters. The performance and robustness of the approach are shown by applying it to a large set of elements crystallizing in a cubic fcc lattice.
△ Less
Submitted 7 December, 2021;
originally announced December 2021.
-
Origins of the hydrogen signal in atom probe tomography
Authors:
Su-Hyun Yoo,
Se-Ho Kim,
Eric Woods,
Baptiste Gault,
Mira Todorova,
Jörg Neugebauer
Abstract:
Atom Probe Tomography (APT) analysis is being actively used to provide near-atomic-scale information on the composition of complex materials in three-dimensions. In recent years, there has been a surge of interest in the technique to investigate the distribution of hydrogen in metals. However, the presence of hydrogen in the analysis of almost all specimens from nearly all material systems has cau…
▽ More
Atom Probe Tomography (APT) analysis is being actively used to provide near-atomic-scale information on the composition of complex materials in three-dimensions. In recent years, there has been a surge of interest in the technique to investigate the distribution of hydrogen in metals. However, the presence of hydrogen in the analysis of almost all specimens from nearly all material systems has caused numerous debates as to its origins and impact on the quantitativeness of the measurement. It is often perceived that most H arises from residual gas ionization, therefore affecting primarily materials with a relatively low evaporation field. In this work, we perform systematic investigations to identify the origin of H residuals in APT experiments by combining density-functional theory (DFT) calculations and APT measurements on an alkali and a noble metal, namely Na and Pt, respectively. We report that no H residual is found in Na metal samples, but in Pt, which has a higher evaporation field, a relatively high signal of H is detected. These results contradict the hypothesis of the H signal being due to direct ionization of residual H$_2$ without much interaction with the specimen's surface. Based on DFT, we demonstrate that alkali metals are thermodynamically less likely to be subject to H contamination under APT-operating conditions compared to transition or noble metals. These insights indicate that the detected H-signal is not only from ionization of residual gaseous H$_2$ alone, but is strongly influenced by material-specific physical properties. The origin of H residuals is elucidated by considering different conditions encountered during APT experiments, specifically, specimen-preparation, transportation, and APT-operating conditions by taking thermodynamic and kinetic aspects into account.
△ Less
Submitted 30 September, 2021;
originally announced September 2021.
-
Multi-State Formulation of the Frozen-Density Embedding Quasi-Diabatization Approach
Authors:
Patrick Eschenbach,
Denis G. Artiukhin,
Johannes Neugebauer
Abstract:
We present a multi-state implementation of the recently developed FDE-diab methodology [J. Chem. Phys., 148 (2018), 214104] in the Serenity program. The new framework extends the original approach such that any number of charge-localized quasi-diabatic states can be coupled, giving an access to calculations of ground and excited state spin-density distributions as well as to excitation energies. W…
▽ More
We present a multi-state implementation of the recently developed FDE-diab methodology [J. Chem. Phys., 148 (2018), 214104] in the Serenity program. The new framework extends the original approach such that any number of charge-localized quasi-diabatic states can be coupled, giving an access to calculations of ground and excited state spin-density distributions as well as to excitation energies. We show that it is possible to obtain results similar to those from correlated wave function approaches such as the complete active space self-consistent field method at much lower computational effort. Additionally, we present a series of approximate computational schemes, which further decrease the overall computational cost and systematically converge to the full FDE-diab solution. The proposed methodology enables computational studies on spin-density distributions and related properties for large molecular systems of biochemical interest.
△ Less
Submitted 4 September, 2021;
originally announced September 2021.
-
Solvation Free Energies in Subsystem Density Functional Theory
Authors:
Moritz Bensberg,
Paul L. Türtscher,
Jan P. Unsleber,
Markus Reiher,
Johannes Neugebauer
Abstract:
For many chemical processes the accurate description of solvent effects are vitally important. Here, we describe a hybrid ansatz for the explicit quantum mechanical description of solute-solvent and solvent-solvent interactions based on subsystem density functional theory and continuum solvation schemes. Since explicit solvent molecules may compromise scalability of the model and transferability o…
▽ More
For many chemical processes the accurate description of solvent effects are vitally important. Here, we describe a hybrid ansatz for the explicit quantum mechanical description of solute-solvent and solvent-solvent interactions based on subsystem density functional theory and continuum solvation schemes. Since explicit solvent molecules may compromise scalability of the model and transferability of the predicted solvent effect, we aim to retain both, for different solutes as well as for different solvents. The key for the transferability is the consistent subsystem decomposition of solute and solvent. The key for the scalability is the performance of subsystem DFT for increasing numbers of subsystems. We investigate molecular dynamics and stationary point sampling of solvent configurations and compare the resulting (Gibbs) free energies to experiment and theoretical methods. We can show that with our hybrid model reaction barriers and reaction energies are accurately reproduced compared to experimental data.
△ Less
Submitted 25 August, 2021;
originally announced August 2021.
-
Understanding alkali contamination in colloidal nanomaterials to unlock grain boundary impurity engineering
Authors:
Se-Ho Kim,
Su-Hyun Yoo,
Poulami Chakraborty,
Jiwon Jeong,
Joohyun Lim,
Ayman A. El-Zoka,
Xuyang Zhou,
Leigh T. Stephenson,
Tilmann Hickel,
Jörg Neugebauer,
Christina Scheu,
Mira Todorova,
Baptiste Gault
Abstract:
Metal nano-aerogels combine a large surface area, a high structural stability, and a high catalytic activity towards a variety of chemical reactions. The performance of such nanostructures is underpinned by the atomic-level distribution of their constituents. Yet monitoring their sub-nanoscale structure and composition to guide property optimization remains extremely challenging. Here, we synthesi…
▽ More
Metal nano-aerogels combine a large surface area, a high structural stability, and a high catalytic activity towards a variety of chemical reactions. The performance of such nanostructures is underpinned by the atomic-level distribution of their constituents. Yet monitoring their sub-nanoscale structure and composition to guide property optimization remains extremely challenging. Here, we synthesized Pd nano-aerogels from a K2PdCl4 precursor and two different NaBH4 reductant concentrations in distilled water. Atom probe tomography reveals that the aerogel is poly-crystalline and that impurities (Na, K) are integrated from the solution into grain boundaries. Ab initio calculations indicate that these impurities preferentially bound to the Pd-metal surface and are ultimately found in grain boundaries forming as the particles coalesce during synthesis, with Na atoms thermodynamically equilibrating with the surrounding solution and K atoms remaining between growing grains. If controlled, impurity integration, i.e. grain boundary decoration, may offer opportunities for designing new nano-aerogels.
△ Less
Submitted 4 November, 2021; v1 submitted 11 June, 2021;
originally announced June 2021.
-
Atomistic deformation behavior of single and twin crystalline Cu nanopillars with preexisting dislocations
Authors:
Won-Seok Ko,
Alexander Stukowski,
Raheleh Hadian,
Ali Nematollahi,
Jong Bae Jeon,
Won Seok Choi,
Gerhard Dehm,
Jörg Neugebauer,
Christoph Kirchlechner,
Blazej Grabowski
Abstract:
Molecular dynamics simulations are performed to investigate the role of a coherent Σ3 (111) twin boundary on the plastic deformation behavior of Cu nanopillars. Our work reveals that the mechanical response of pillars with and without the twin boundary is decisively driven by the characteristics of initial dislocation sources. In the condition of comparably large pillar size and abundant initial m…
▽ More
Molecular dynamics simulations are performed to investigate the role of a coherent Σ3 (111) twin boundary on the plastic deformation behavior of Cu nanopillars. Our work reveals that the mechanical response of pillars with and without the twin boundary is decisively driven by the characteristics of initial dislocation sources. In the condition of comparably large pillar size and abundant initial mobile dislocations, overall yield and flow stresses are controlled by the longest, available mobile dislocation. An inverse correlation of the yield and flow stresses with the length of the longest dislocation is established, and its extrapolation agrees well with experimental yield stress data. The experimentally reported subtle differences in yield and flow stresses between pillars with and without the twin boundary are thus likely related to the maximum lengths of the mobile dislocations. In the condition of comparably small pillar size, for which a reduction of mobile dislocations during heat treatment and mechanical loading occurs, the mechanical response of pillars with and without the twin boundary can be clearly distinguished. Dislocation starvation during deformation is more clearly present in pillars without the twin boundary than in pillars with the twin boundary because the twin boundary acts as a pinning surface for the dislocation network.
△ Less
Submitted 8 April, 2021; v1 submitted 8 April, 2021;
originally announced April 2021.
-
Revealing atomic-scale vacancy-solute interaction in nickel
Authors:
Felipe F. Morgado,
Shyam Katnagallu,
Christoph Freysoldt,
Benjamin Klaes,
François Vurpillot,
Jörg Neugebauer,
Dierk Raabe,
Steffen Neumeier,
Baptiste Gault,
Leigh T. Stephenson
Abstract:
Imaging individual vacancies in solids and revealing their interactions with solute atoms remains one of the frontiers in microscopy and microanalysis. Here we study a creep-deformed binary Ni-2 at.% Ta alloy. Atom probe tomography reveals a random distribution of Ta. Field ion microscopy, with contrast interpretation supported by density-functional theory and time-of-flight mass spectrometry, evi…
▽ More
Imaging individual vacancies in solids and revealing their interactions with solute atoms remains one of the frontiers in microscopy and microanalysis. Here we study a creep-deformed binary Ni-2 at.% Ta alloy. Atom probe tomography reveals a random distribution of Ta. Field ion microscopy, with contrast interpretation supported by density-functional theory and time-of-flight mass spectrometry, evidences a positive correlation of tantalum with vacancies. Our results support solute-vacancy binding, which explains improvement in creep resistance of Ta-containing Ni-based superalloys and helps guide future material design strategies.
△ Less
Submitted 16 June, 2021; v1 submitted 2 March, 2021;
originally announced March 2021.
-
Ab initio investigations of point and complex defect structures in B2-FeAl
Authors:
Halil İbrahim Sözen,
Tilmann Hickel,
Jörg Neugebauer
Abstract:
In this work we have studied the defect structure and corresponding defect concentration investigations through the theoretical, experimental and computational works on B2-type Fe-Al alloys. We have used ab initio framework in order to investigate the defect structure. To have a proper explanation for high defect concentration in B2-FeAl, we did not confine with point defect, but extend the work o…
▽ More
In this work we have studied the defect structure and corresponding defect concentration investigations through the theoretical, experimental and computational works on B2-type Fe-Al alloys. We have used ab initio framework in order to investigate the defect structure. To have a proper explanation for high defect concentration in B2-FeAl, we did not confine with point defect, but extend the work on defect complexes. The possible defect formation energies were calculated with the dependence of chemical potential and carefully investigated against supercell size and the effect of magnetism. The calculations revealed that the double Fe antisite at Fe rich condition, the single Fe vacancy at intermediate region (i.e in the stoichiometry) and the double Al antisite is the dominant defect close to Al rich condition, where mainly Al rich region was unstable. From the obtained defect formation energies, defect concentrations were calculated at different temperatures with respect to Al concentration for B2-FeAl. It has been found that increasing Al content and temperature gradually leads to increase in the vacancy content. It has also seen that the dominant defect for all temperature ranges was the single Fe vacancy at the exact stoichiometry and the highest single Fe vacancy content detected with 1.6 % at 1450 K.
△ Less
Submitted 15 January, 2021;
originally announced January 2021.
-
Comparison of coherent phonon generation by electronic and ionic Raman scattering in LaAlO$_3$
Authors:
Martin J. Neugebauer,
Dominik M. Juraschek,
Matteo Savoini,
Pascal Engeler,
Larissa Boie,
Elsa Abreu,
Nicola A. Spaldin,
Steven L. Johnson
Abstract:
In ionic Raman scattering, infrared-active phonons mediate a scattering process that results in the creation or destruction of a Raman-active phonon. This mechanism relies on nonlinear interactions between phonons and has in recent years been associated with a variety of emergent lattice-driven phenomena in complex transition-metal oxides, but the underlying mechanism is often obscured by the pres…
▽ More
In ionic Raman scattering, infrared-active phonons mediate a scattering process that results in the creation or destruction of a Raman-active phonon. This mechanism relies on nonlinear interactions between phonons and has in recent years been associated with a variety of emergent lattice-driven phenomena in complex transition-metal oxides, but the underlying mechanism is often obscured by the presence of multiple coupled order parameters in play. Here, we use time-resolved spectroscopy to compare coherent phonons generated by ionic Raman scattering with those created by more conventional electronic Raman scattering on the nonmagnetic and non-strongly-correlated wide band-gap insulator LaAlO$_3$. We find that the oscillatory amplitude of the low-frequency Raman-active $E_g$ mode exhibits a sharp peak when we tune our pump frequency into resonance with the high-frequency infrared-active $E_u$ mode, consistent with first-principles calculations. Our results suggest that ionic Raman scattering can strongly dominate electronic Raman scattering in wide band-gap insulating materials. We also see evidence of competing scattering channels at fluences above 28~mJ/cm$^2$ that alter the measured amplitude of the coherent phonon response.
△ Less
Submitted 24 October, 2020;
originally announced October 2020.
-
Ab-initio investigation of lattice distortions in response to van der Waals interactions in FeSe
Authors:
Felix Lochner,
Ilya M. Eremin,
Tilmann Hickel,
Jörg Neugebauer
Abstract:
The electronic structure in unconventional superconductors holds a key to understand the momentum-dependent pairing interactions and the resulting superconducting gap function. In superconducting Fe-based chalcogenides, there have been controversial results regarding the importance of the $k_z$ dependence of the electronic dispersion, the gap structure and the pairing mechanisms of iron-based supe…
▽ More
The electronic structure in unconventional superconductors holds a key to understand the momentum-dependent pairing interactions and the resulting superconducting gap function. In superconducting Fe-based chalcogenides, there have been controversial results regarding the importance of the $k_z$ dependence of the electronic dispersion, the gap structure and the pairing mechanisms of iron-based superconductivity. Here, we present a detailed investigation of the van der Waals interaction in FeSe and its interplay with magnetic disorder and real space structural properties. Using density functional theory we show that they need to be taken into account upon investigation of the 3-dimensional effects, including non-trivial topology, of FeSe$_{1-x}$Te$_x$ and FeSe$_{1-x}$S$_x$ systems. In addition, the impact of paramagnetic (PM) disorder is considered within the spin-space average approach. Our calculations show that the PM relaxed structure supports the picture of different competing ordered magnetic states in the nematic regime, yielding magnetic frustration.
△ Less
Submitted 8 July, 2020;
originally announced July 2020.
-
Correlation between electronic and structural orders in 1T-TiSe2
Authors:
Hiroki Ueda,
Michael Porer,
José R. L. Mardegan,
Sergii Parchenko,
Namrata Gurung,
Federica Fabrizi,
Mahesh Ramakrishnan,
Larissa Boie,
Martin Josef Neugebauer,
Bulat Burganov,
Max Burian,
Steven Lee Johnson,
Kai Rossnagel,
Urs Staub
Abstract:
The correlation between electronic and crystal structures of 1T-TiSe2 in the charge density wave (CDW) state is studied by x-ray diffraction. Three families of reflections are used to probe atomic displacements and the orbital asymmetry in Se. Two distinct onset temperatures are found, TCDW and a lower T* indicative for an onset of Se out-of-plane atomic displacements. T* coincides with a DC resis…
▽ More
The correlation between electronic and crystal structures of 1T-TiSe2 in the charge density wave (CDW) state is studied by x-ray diffraction. Three families of reflections are used to probe atomic displacements and the orbital asymmetry in Se. Two distinct onset temperatures are found, TCDW and a lower T* indicative for an onset of Se out-of-plane atomic displacements. T* coincides with a DC resistivity maximum and the onset of the proposed gyrotropic electronic structure. However, no indication for chirality is found. The relation between the atomic displacements and the transport properties is discussed in terms of Ti 3d and Se 4p states that only weakly couple to the CDW order.
△ Less
Submitted 16 June, 2020;
originally announced June 2020.
-
Dielectric properties of nano-confined water: a canonical thermopotentiostat approach
Authors:
F. Deißenbeck,
C. Freysoldt,
M. Todorova,
J. Neugebauer,
S. Wippermann
Abstract:
We introduce a novel approach to sample the canonical ensemble at constant temperature and applied electric potential. Our approach can be straightforwardly implemented into any density-functional theory code. Using thermopotentiostat molecular dynamics simulations allows us to compute the dielectric constant of nano-confined water without any assumptions for the dielectric volume. Compared to the…
▽ More
We introduce a novel approach to sample the canonical ensemble at constant temperature and applied electric potential. Our approach can be straightforwardly implemented into any density-functional theory code. Using thermopotentiostat molecular dynamics simulations allows us to compute the dielectric constant of nano-confined water without any assumptions for the dielectric volume. Compared to the commonly used approach of calculating dielectric properties from polarization fluctuations, our thermopotentiostat technique reduces the required computational time by two orders of magnitude.
△ Less
Submitted 14 August, 2020; v1 submitted 18 March, 2020;
originally announced March 2020.
-
Ab initio Description of Bond-Breaking in Large Electric Fields
Authors:
Michael Ashton,
Arpit Mishra,
Joerg Neugebauer,
Christoph Freysoldt
Abstract:
Strong ($10^{10}$ V/m) electric fields capable of inducing atomic bond-breaking represent a powerful tool for surface chemistry. However, their exact effects are difficult to predict due to a lack of suitable tools to probe their associated atomic-scale mechanisms. Here we introduce a generalized dipole correction for charged repeated-slab models that controls the electric field on both sides of t…
▽ More
Strong ($10^{10}$ V/m) electric fields capable of inducing atomic bond-breaking represent a powerful tool for surface chemistry. However, their exact effects are difficult to predict due to a lack of suitable tools to probe their associated atomic-scale mechanisms. Here we introduce a generalized dipole correction for charged repeated-slab models that controls the electric field on both sides of the slab, thereby enabling direct theoretical treatment of field-induced bond-breaking events. As a prototype application, we consider field evaporation from a kinked W surface. We reveal two qualitatively different desorption mechanisms that can be selected by the magnitude of the applied field.
△ Less
Submitted 21 February, 2020; v1 submitted 7 February, 2020;
originally announced February 2020.
-
Segregation--Assisted Spinodal and Transient Spinodal Phase Separation at Grain Boundaries
Authors:
Reza Darvishi Kamachali,
Alisson Kwiatkowski da Silva,
Eunan McEniry,
Dirk Ponge,
Baptiste Gault,
Joerg Neugebauer,
Dierk Raabe
Abstract:
Segregation to grain boundaries affects their cohesion, corrosion and embrittlement and plays a critical role in heterogeneous nucleation. In order to quantitatively study segregation and phase separation at grain boundaries, we derive a density-based phase-field model. In this model, we describe the grain boundary free energy based on available bulk thermodynamic data while an atomic grain bounda…
▽ More
Segregation to grain boundaries affects their cohesion, corrosion and embrittlement and plays a critical role in heterogeneous nucleation. In order to quantitatively study segregation and phase separation at grain boundaries, we derive a density-based phase-field model. In this model, we describe the grain boundary free energy based on available bulk thermodynamic data while an atomic grain boundary density is obtained using atomistic simulations. To benchmark the performance of the model, we study Mn grain boundary segregation in the Fe--Mn system. 3D simulation results are compared against atom probe tomography measurements. We show that a continuous increase in the alloy composition results in a discontinuous jump in the Mn grain boundary segregation. This jump corresponds to an interfacial spinodal phase separation. For alloy compositions above the interfacial spinodal, we found a transient spinodal phase separation phenomenon which opens opportunities for knowledge-based microstructure design through the chemical manipulation of grain boundaries. The proposed density-based model provides a powerful tool to study thermodynamics and kinetics of segregation and phase separation at grain boundaries.
△ Less
Submitted 14 November, 2019; v1 submitted 20 May, 2019;
originally announced May 2019.
-
Imaging individual solute atoms at crystalline imperfections in metals
Authors:
Shyam Katnagallu,
Leigh T. Stephenson,
Isabelle Mouton,
Christoph Freysoldt,
Aparna P. A. Subramanyam,
Jan Jenke,
Alvin N. Ladines,
Steffen Neumeier,
Thomas Hammerschmidt,
Ralf Drautz,
Jörg Neugebauer,
François Vurpillot,
Dierk Raabe,
Baptiste Gault
Abstract:
Directly imaging all atoms constituting a material and, maybe more importantly, crystalline defects that dictate materials' properties, remains a formidable challenge. Here, we propose a new approach to chemistry-sensitive field-ion microscopy (FIM) combining contrast interpretation from density-functional theory (DFT) and elemental identification enabled by time-of-flight mass-spectrometry and da…
▽ More
Directly imaging all atoms constituting a material and, maybe more importantly, crystalline defects that dictate materials' properties, remains a formidable challenge. Here, we propose a new approach to chemistry-sensitive field-ion microscopy (FIM) combining contrast interpretation from density-functional theory (DFT) and elemental identification enabled by time-of-flight mass-spectrometry and data mining. Analytical-FIM has true atomic resolution and we demonstrate how the technique can reveal the presence of individual solute atoms at specific positions in the microstructure. The performance of this new technique is showcased in revealing individual Re atoms at crystalline defects formed in Ni during creep deformation. The atomistic details offered by A-FIM allowed us to directly compare our results with simulations, and to tackle a long-standing question of how Re extends lifetime of Ni-based superalloys in service at high-temperature.
△ Less
Submitted 31 March, 2019; v1 submitted 8 March, 2019;
originally announced March 2019.
-
Ab initio vibrational free energies including anharmonicity for multicomponent alloys
Authors:
Blazej Grabowski,
Yuji Ikeda,
Fritz Körmann,
Christoph Freysoldt,
Andrew Ian Duff,
Alexander Shapeev,
Jörg Neugebauer
Abstract:
A density-functional-theory based approach to efficiently compute numerically exact vibrational free energies - including anharmonicity - for chemically complex multicomponent alloys is developed. It is based on a combination of thermodynamic integration and a machine-learning potential. We demonstrate the performance of the approach by computing the anharmonic free energy of the prototypical five…
▽ More
A density-functional-theory based approach to efficiently compute numerically exact vibrational free energies - including anharmonicity - for chemically complex multicomponent alloys is developed. It is based on a combination of thermodynamic integration and a machine-learning potential. We demonstrate the performance of the approach by computing the anharmonic free energy of the prototypical five-component VNbMoTaW refractory high entropy alloy.
△ Less
Submitted 28 February, 2019;
originally announced February 2019.
-
Optical control of vibrational coherence triggered by an ultrafast phase transition
Authors:
M. J. Neugebauer,
T. Huber,
M. Savoini,
E. Abreu,
V. Esposito,
M. Kubli,
L. Rettig,
E. Bothschafter,
S. Grübel,
T. Kubacka,
J. Rittmann,
G. Ingold,
P. Beaud,
D. Dominko,
J. Demsar,
S. L. Johnson
Abstract:
Femtosecond time-resolved x-ray diffraction is employed to study the dynamics of the periodic lattice distortion (PLD) associated with the charge-density-wave (CDW) in K0.3MoO3. Using a multi-pulse scheme we show the ability to extend the lifetime of coherent oscillations of the PLD about the undistorted structure through re-excitation of the electronic states. This suggests that it is possible to…
▽ More
Femtosecond time-resolved x-ray diffraction is employed to study the dynamics of the periodic lattice distortion (PLD) associated with the charge-density-wave (CDW) in K0.3MoO3. Using a multi-pulse scheme we show the ability to extend the lifetime of coherent oscillations of the PLD about the undistorted structure through re-excitation of the electronic states. This suggests that it is possible to enter a regime where the symmetry of the potential energy landscape corresponds to the high symmetry phase but the scattering pathways that lead to the damping of coherent dynamics are still controllable by altering the electronic state population. The demonstrated control over the coherence time offers new routes for manipulation of coherent lattice states.
△ Less
Submitted 1 February, 2019;
originally announced February 2019.
-
Ultrafast transient increase of oxygen octahedral rotations in a perovskite
Authors:
M. Porer,
M. Fechner,
M. Kubli,
M. J. Neugebauer,
S. Parchenko,
V. Esposito,
A. Narayan,
N. A. Spaldin,
R. Huber,
M. Radovic,
E. M. Bothschafter,
J. M. Glownia,
T. Sato,
S. Song,
S. L. Johnson,
U. Staub
Abstract:
The ability to control the structure of a crystalline solid on ultrafast timescales bears enormous potential for information storage and manipulation or generating new functional states of matter [1]. In many materials where the ultrafast control of crystalline structures has been explored, optical excitation pushes materials towards their less ordered high temperature phase [2{9] as electronicall…
▽ More
The ability to control the structure of a crystalline solid on ultrafast timescales bears enormous potential for information storage and manipulation or generating new functional states of matter [1]. In many materials where the ultrafast control of crystalline structures has been explored, optical excitation pushes materials towards their less ordered high temperature phase [2{9] as electronically driven ordered phases melt and possible concomitant structural modifications relax. Nonetheless, for a few select materials it has been shown that photoexcitation can slightly enhance the amplitude of an electronic ordering phenomenon (i.e. its electronic order parameter) [9{13]. Here we show via femtosecond hard X-ray diffraction that photodoping of the perovskite EuTiO3 transiently increases the order parameter associated with a purely structural [14] phase transition represented by the antiferrodistortive rotation of the oxygen octahedra. This can be understood from an ultrafast charge-transfer induced reduction of the Goldschmidt tolerance factor [15], which is a fundamental control parameter for the properties of perovskites
△ Less
Submitted 24 January, 2019;
originally announced January 2019.
-
Impact of lattice relaxations on phase transitions in a high-entropy alloy studied by machine-learning potentials
Authors:
Tatiana Kostiuchenko,
Fritz Körmann,
Jörg Neugebauer,
Alexander Shapeev
Abstract:
Recently, high-entropy alloys (HEAs) have attracted wide attention due to their extraordinary materials properties. A main challenge in identifying new HEAs is the lack of efficient approaches for exploring their huge compositional space. Ab initio calculations have emerged as a powerful approach that complements experiment. However, for multicomponent alloys existing approaches suffer from the ch…
▽ More
Recently, high-entropy alloys (HEAs) have attracted wide attention due to their extraordinary materials properties. A main challenge in identifying new HEAs is the lack of efficient approaches for exploring their huge compositional space. Ab initio calculations have emerged as a powerful approach that complements experiment. However, for multicomponent alloys existing approaches suffer from the chemical complexity involved. In this work we propose a method for studying HEAs computationally. Our approach is based on the application of machine-learning potentials based on ab initio data in combination with Monte Carlo simulations. The high efficiency and performance of the approach are demonstrated on the prototype bcc NbMoTaW HEA. The approach is employed to study phase stability, phase transitions, and chemical short-range order. The importance of including local relaxation effects is revealed: they significantly stabilize single-phase formation of bcc NbMoTaW down to room temperature. Finally, a so-far unknown mechanism that drives chemical order due to atomic relaxation at ambient temperatures is discovered.
△ Less
Submitted 6 May, 2019; v1 submitted 25 October, 2018;
originally announced October 2018.
-
{110} planar faults in strained bcc metals: Origins and implications of a commonly observed artifact of classical potentials
Authors:
Johannes J. Möller,
Matous Mrovec,
Ivan Bleskov,
Jörg Neugebauer,
Thomas Hammerschmidt,
Ralf Drautz,
Christian Elsässer,
Tilmann Hickel,
Erik Bitzek
Abstract:
Large-scale atomistic simulations with classical potentials can provide valuable insights into microscopic deformation mechanisms and defect-defect interactions in materials. Unfortunately, these assets often come with the uncertainty of whether the observed mechanisms are based on realistic physical phenomena or whether they are artifacts of the employed material models. One such example is the o…
▽ More
Large-scale atomistic simulations with classical potentials can provide valuable insights into microscopic deformation mechanisms and defect-defect interactions in materials. Unfortunately, these assets often come with the uncertainty of whether the observed mechanisms are based on realistic physical phenomena or whether they are artifacts of the employed material models. One such example is the often reported occurrence of stable planar faults (PFs) in body-centered cubic (bcc) metals subjected to high strains, e.g., at crack tips or in strained nano-objects. In this paper, we study the strain dependence of the generalized stacking fault energy (GSFE) of {110} planes in various bcc metals with material models of increasing sophistication, i.e., (modified) embedded atom method, angular-dependent, Tersoff, and bond-order potentials as well as density functional theory. We show that under applied tensile strains the GSFE curves of many classical potentials exhibit a local minimum which gives rise to the formation of stable PFs. These PFs do not appear when more sophisticated material models are used and have thus to be regarded as artifacts of the potentials. We demonstrate that the local GSFE minimum is not formed for reasons of symmetry and we recommend including the determination of the strain-dependent (110) GSFE as a benchmark for newly developed potentials.
△ Less
Submitted 3 October, 2018;
originally announced October 2018.
-
Ultrafast relaxation dynamics of the antiferrodistortive phase in Ca doped SrTiO3
Authors:
M. Porer,
M. Fechner,
E. Bothschafter,
L. Rettig,
M. Savoini,
V. Esposito,
J. Rittmann,
M. Kubli,
M. J. Neugebauer,
E. Abreu,
T. Kubacka,
T. Huber,
G. Lantz,
S. Parchenko,
S. Grübel,
A. Paarmann,
J. Noack,
P. Beaud,
G. Ingold,
U. Aschauer,
S. L. Johnson,
U. Staub
Abstract:
The ultrafast dynamics of the octahedral rotation in Ca:SrTiO3 is studied by time resolved x-ray diffraction after photo excitation over the band gap. By monitoring the diffraction intensity of a superlattice reflection that is directly related to the structural order parameter of the soft-mode driven antiferrodistortive phase in Ca:SrTiO3, we observe a ultrafast relaxation on a 0.2 ps timescale o…
▽ More
The ultrafast dynamics of the octahedral rotation in Ca:SrTiO3 is studied by time resolved x-ray diffraction after photo excitation over the band gap. By monitoring the diffraction intensity of a superlattice reflection that is directly related to the structural order parameter of the soft-mode driven antiferrodistortive phase in Ca:SrTiO3, we observe a ultrafast relaxation on a 0.2 ps timescale of the rotation of the oxygen octahedron, which is found to be independent of the initial temperaure despite large changes in the corresponding soft-mode frequency. A further, much smaller reduction on a slower picosecond timescale is attributed to thermal effects. Time-dependent density-functional-theory calculations show that the fast response can be ascribed to an ultrafast displacive modification of the soft-mode potential towards the normal state, induced by holes created in the oxygen 2p states.
△ Less
Submitted 1 May, 2018;
originally announced May 2018.