Ab Initio Calculations Research Articles

First principles validation of energy barriers in Ni75Al25

Precipitates in nickel-based superalloys form during heat treatment on a time scale inaccessible to direct molecular dynamics simulation, but can be studied using kinetic Monte Carlo (KMC) modelling. This requires reliable values for the barrier energies separating distinct configurations over the trajectory of the system. In this study, we validate vacancy migration barriers found with the Activation-Relaxation Technique nouveau (ARTn) method in partially ordered Ni75Al25 with a monovacancy using published potentials for the atomic interactions against first-principles methods. In a first step, we confirm that the ARTn barrier energies agree with those determined with the nudged elastic band (NEB) method. As the number of atoms used in those calculations is too great for direct ab initio calculations, we cut the cell size to 255 atoms, thus controlling finite size effects. We then use the plane-wave density functional theory code CASTEP and its inbuilt NEB method in the smaller cells. This provides us with a continuous validation chain from first principles to KMC simulations with interatomic potentials (IPs). We evaluate the barrier energies of five further IPs with NEB, demonstrating that none yields values with sufficient reliability for KMC simulations, with some of them failing completely. This is a first step towards quantifying the errors incurred in KMC simulations of precipitate formation and evolution.

Effectiveness of ab initio molecular dynamics in simulating EXAFS spectra from layered systems.

The simulation of EXAFS spectra of thin films via ab initio methods is discussed. The procedure for producing the spectra is presented as well as an application to a two-dimensional material (WSe2) where the effectiveness of this method in reproducing the spectrum and the linear dichroic response is shown. A series of further examples in which the method has been employed for the structural determination of materials are given.

High-performance transistors based on monolayer all-inorganic two-dimensional halide perovskite Cs2PbI2Cl2 and its optoelectronic application

The two-dimensional (2D) Ruddlesden−Popper phase all-inorganic halide perovskite Cs2PbI2Cl2 has attracted tremendous attention because of its high carrier mobility, strong light-absorbing ability, and excellent environmental stability, possessing enormous potential for use as a high-performance optoelectronic device. Here, we present two transistor devices based on monolayer (ML) Cs2PbI2Cl2 and reveal their performance and current transport mechanisms through ab initio quantum transport calculations. The n-type 10 nm metal-oxide-semiconductor field-effect transistor shows a high on-state current Ion of 3160 μA/μm, and excellent performance regarding subthreshold swing, intrinsic delay time (τ), and power consumption, which completely match the International Roadmap for Device and Systems (2021 version) targets for high-performance devices and low-power devices in 2028. Moreover, phototransistors based on monolayer Cs2PbI2Cl2 also display a strong response to UV light, with a light-to-dark current ratio reaching a maximum of 104. Meanwhile, a high specific detectivity of 2.45 × 1011 Jones is obtained. The results of our study suggest that ML Cs2PbI2Cl2 is a promising candidate for high-performance transistors.

Ab Initio Prediction of Vapor Pressure for Diverse Atomic Layer Deposition Precursors.

Understanding the saturated vapor pressure (Pvap) is vital for evaluating atomic layer deposition (ALD) precursors, as it directly influences the ALD temperature window and, by extension, the processability of compounds. The early estimation of vapor pressure ranges is crucial during the initial stages of novel precursor design, reducing the reliance on empirical synthesis or experimentation. However, predicting vapor pressure through computer simulations is often impeded by the scarcity of suitable empirical force fields for molecular dynamics simulations. This challenge is further compounded by the diverse chemical substances and the introduction of new elements into modern ALD processes, necessitating robust force fields that can accommodate metals, organics, and halides. In response, this study introduces a novel approach utilizing a quantum mechanically derived force field for the prediction of vapor pressure across a wide spectrum of potential ALD precursors. This approach enables the creation of system-specific force fields through parametrization based on ab initio calculations for a single molecule. We develop a comprehensive workflow to simulate both liquid and gaseous equilibrium phases, allowing the calculation of vapor pressure across a wide temperature range. Our methodology has been validated with a diverse set of ALD precursors, demonstrating its robustness in predicting Pvap at specified temperatures. The approach yields a Pearson's correlation coefficient (R2) greater than 0.9 on a logarithmic scale and a root-mean-squared deviation in self-solvation-free energies as low as 1.3 kcal mol-1. This innovative workflow, which does not require any prior experimental data, marks a significant advancement in the computer-aided design of novel ALD precursors, paving the way for accelerating developments in technology.

Evaluating the thermal stability of chemicals and systems: A review

AbstractIn the realm of chemical processing, particularly at the industrial scale, safety is of utmost importance. A predominant factor causing accidents within the chemical industry is runaway phenomena, primarily initiated by uncontrolled exothermic reactions. This review critically examines the often‐overlooked decomposition mechanisms as a significant contributor to thermal energy release, necessitating a comprehensive revision and understanding of both experimental and theoretical strategies for assessing thermal degradation. Key to this discourse is the explication of calorimetry as the principal experimental technique, alongside ab initio quantum chemistry simulations as a robust theoretical framework for quantifying the most relevant properties. However, more than mere cognisance of these methodologies is required for a meticulous thermal stability assessment. The review emphasizes identifying and quantifying fundamental parameters through experimental and theoretical investigations. Only upon acquiring these parameters, including kinetic, thermodynamic, onset, and peak characteristics of the exothermic decomposition reactions, can one effectively mitigate risks and hazards in designing and optimizing chemical processes and apparatus. Furthermore, this review delineates qualitative and quantitative methodologies for hazard assessment, proffering strategies for estimating safe operational conditions and sizing relief devices. The paper culminates in exploring future trajectories in thermal stability assessments, focusing on emerging applications in lithium‐ion batteries, electrolyzers, electrified reactors, ionic liquids, artificial intelligence and machine learning approaches. Thus, the paper underlines the evolving landscape of thermal risk management in contemporary and future chemical industries.

Theoretical Study of Ground-State Barium-Rare Gas Van der Waals Complexes: Combining Rule Modeling and Ab Initio Calculations.

The present study aims to generate the potential energy curves (PECs) and spectroscopic constants for barium alkaline earth (AE) atoms interacting with rare gas (RG) atoms (He, Ne, Ar, Kr, and Xe). The study focuses on investigating the van der Waals bonds that characterize the interactions between alkaline-earth metals and RG atoms, with a specific emphasis on employing the Tang and Toennies (TT) potential model, known to accurately describe such interactions. The TT potential model was employed in conjunction with combining rules to calculate its parameters, which include dispersion coefficients C 2n and Born-Mayer constants A and b. Additionally, we have conducted high-level ab initio calculations at the CCSD(T) level for all Ba-RG ground states. Obtained PECs from both methods have been used to evaluate the spectroscopic properties D e, R e, ωe, B e, and ωeχe. Our findings reveal that the derived spectroscopic constants from the TT model exhibit good agreement with the results obtained from CCSD(T) calculations and with other available theoretical studies. Furthermore, to gain insights into the relative differences among AE-RG species, we calculated the κ parameter for AE-RG and AE+-RG (AE = Sr, Ca, Mg, Ba; RG = He-Xe) complexes. It is found that except for the case of Ba-RG and Ba+-RG, the κ values within the same series, AE-RG and AE+-RG, are remarkably close to each other.

Toward Realistic Models of the Electrocatalytic Oxygen Evolution Reaction.

The electrocatalytic oxygen evolution reaction (OER) supplies the protons and electrons needed to transform renewable electricity into chemicals and fuels. However, the OER is kinetically sluggish; it operates at significant rates only when the applied potential far exceeds the reversible voltage. The origin of this overpotential is hidden in a complex mechanism involving multiple electron transfers and chemical bond making/breaking steps. Our desire to improve catalytic performance has then made mechanistic studies of the OER an area of major scientific inquiry, though the complexity of the reaction has made understanding difficult. While historically, mechanistic studies have relied solely on experiment and phenomenological models, over the past twenty years ab initio simulation has been playing an increasingly important role in developing our understanding of the electrocatalytic OER and its reaction mechanisms. In this Review we cover advances in our mechanistic understanding of the OER, organized by increasing complexity in the way through which the OER is modeled. We begin with phenomenological models built using experimental data before reviewing early efforts to incorporate ab initio methods into mechanistic studies. We go on to cover how the assumptions in these early ab initio simulations─no electric field, electrolyte, or explicit kinetics─have been relaxed. Through comparison with experimental literature, we explore the veracity of these different assumptions. We summarize by discussing the most critical open challenges in developing models to understand the mechanisms of the OER.

Ab initio study of surfaces of lead and tin based metal halide perovskite structures

Abstract Hybrid perovskite materials, known for their potential in cost-effective optoelectronic applications, face a knowledge gap in crucial areas, particularly the atomic-level properties of the surface. This study addresses this challenge by refining ab initio methods for characterizing surface structures of cubic methylammonium lead bromide and methylammonium tin bromide (MAMeBr3 with Me = Sn, Pb), avoiding superficial restrictions in atomic movement during geometry optimization. The resulting structures confirmed nearly random MA+ molecule alignment, comparable to real-world experimental conditions. Calculating surface energies for these structures with crystal orientations {100} and {110}, each with different terminations, provides valuable insights into structural properties. Using a carefully chosen thermodynamic reference state, mimicking experimental conditions enables a thermodynamic discussion and facilitates the modulation of the MeBr2 component’s chemical potential. This modulation, in turn, allows for the prediction of crystal morphologies, as illustrated by Wulff’s construction. This approach establishes a crucial link between theoretical predictions and experimental conditions, shedding light on the complexities of hybrid perovskite materials.

Thermodynamic study of tin tetraiodide (SnI4) sublimation by effusion techniques

This study presents the first investigation of the sublimation behavior of tin tetraiodide, SnI4, using effusion-based techniques, within a low temperature range (313–340) K. The temperature range covered in the experiments was lower than in previously reported studies based on static methods. Knudsen Effusion Mass Loss (KEML) measurements were performed in the range of (317.1–339.6) K using effusion cells with different orifice sizes. The vapor pressures were measured in the range (0.13–1.13) Pa and were found to be independent of the orifice size. The standard molar enthalpy and Gibbs energy of sublimation at 298.15 K obtained by the Clarke and Glew fit of experimental data are (88.1 ± 0.9) kJ⋅mol−1 and (38.96 ± 0.08) kJ⋅mol−1, respectively. Knudsen Effusion Mass Spectrometry (KEMS) experiments were also performed in the range (313.3–331.7) K, resulting in a sublimation enthalpy value in good agreement with the KEML values and not negligibly higher vapor pressure values. KEMS vapor pressure data were also analyzed by the third-law method. A comparison of our experimental results with the literature data available for both sublimation and evaporation properties of SnI4 is reported. Additionally, ancillary DFT and ab initio calculations were performed to estimate the molecular properties of SnI4(g) and the extent of the gas-phase dissociation to SnI2 and I2.

Unraveling the Geometry-Driven C═C Epoxidation and C-H Hydroxylation Reactivity of Tetra-Coordinated Nonheme Iron(IV)-Oxo Complexes.

The electronic structure and reactivity of tetra-coordinated nonheme iron(IV)-oxo complexes have remained unexplored for years. The recent synthesis of a closed-shell iron(IV)-oxo complex [(quinisox)FeIV(O)]+ (1) has set up a platform to understand how such complexes compare with the celebrated open-shell iron-oxo chemistry. Herein, using density functional theory and ab initio calculations, we present an in-depth electronic structure investigation of the C═C epoxidation [oxygen atom transfer (OAT)] and C-H hydroxylation [hydrogen atom transfer (HAT)] reactivity of 1. Using a solvent-coordinated geometry of 1 (1') and other potential tetra-coordinated iron(IV)-oxo complexes bearing rigid ligands (2 and 3), we established the geometric origin of spin-state energetics and reactivity of 1. Complex 1 featuring a strong Fe-O bond exhibits OAT and HAT reactivity in its quintet state. The lowest quintet OAT pathway has a lower barrier by ∼4 kcal/mol than the quintet HAT pathway, corroborating the experimentally observed gas-phase OAT reactivity preference. A conventional HAT reactivity preference for 2 and a comparable OAT and HAT reactivity for 3 are observed. This further supports the geometry-driven reactivity preference for 1. Noncovalent interaction analyses reveal a pronounced π-π interaction between the substrate and ligand in the OAT transition state, rationalizing the origin of the observed reactivity preference for 1.

Pressure-Induced Re-Entrant Superconductivity in Transition Metal Dichalcogenide TiSe2.

Transition metal dichalcogenide TiSe2 exhibits a superconducting dome within a low pressure range of 2-4GPa, which peaks with the maximal transition temperature Tc of ≈1.8 K. Here it is reported that applying high pressure induces a new superconducting state in TiSe2, which starts at ≈16GPa with a substantially higher Tc that reaches 5.6 K at ≈21.5GPa with no sign of decline. Combining high-throughput first-principles structure search, X-ray diffraction, and Raman spectroscopy measurements up to 30GPa, It is found that TiSe2 undergoes a first-order structural transition from the 1T phase under ambient pressure to a new 4O phase under high pressure. Comparative ab initio calculations reveal that while the conventional phonon-mediated pairing mechanism may account for the superconductivity observed in 1T-TiSe2 under low pressure, the electron-phonon coupling of 4O-TiSe2 is too weak to induce a superconducting state whose transition temperature is as high as 5.6 K under high pressure. The new superconducting state found in pressurized TiSe2 requires further study on its underlying mechanism.

Effect of Dynamical Motion in ab Initio Calculations of Solid-State Nuclear Magnetic and Nuclear Quadrupole Resonance Spectra.

Solid-state nuclear magnetic resonance (SSNMR) and nuclear quadrupole resonance (NQR) spectra provide detailed information about the electronic and atomic structure of solids. Modern ab initio methods such as density functional theory (DFT) can be used to calculate NMR and NQR spectra from first-principles, providing a meaningful avenue to connect theory and experiment. Prediction of SSNMR and NQR spectra from DFT relies on accurate calculation of the electric field gradient (EFG) tensor associated with the potential of electrons at the nuclear centers. While static calculations of EFGs are commonly seen in the literature, the effects of dynamical motion of atoms in molecules and solids have been less explored. In this study, we develop a method to calculate EFGs of solids while taking into account the dynamics of atoms through DFT-based molecular dynamics simulations. The method we develop is general, in the sense that it can be applied to any material at any desired temperature and pressure. Here, we focus on application of the method to NaNO2 and study in detail the EFGs of 14N, 17O, and 23Na. We find in the cases of 14N and 17O that the dynamical motion of the atoms can be used to calculate mean EFGs that are in closer agreement with experiments than those of static calculations. For 23Na, we find a complex behavior of the EFGs when atomic motion is incorporated that is not at all captured in static calculations. In particular, we find a distribution of EFGs that is influenced strongly by the local (changing) bond environment, with a pattern that reflects the coordination structure of 23Na. We expect the methodology developed here to provide a path forward for understanding materials in which static EFG calculations do not align with experiments.

Charge transfer effects in (HfNbTiVZr)C—Shown by ab initio calculations and X‐ray photoelectron spectroscopy

AbstractConsidering charge transfer effects and the variability of the bonding between elements with different electronegativity opens up a deeper understanding of the electronic structure and as a result many of the properties in high entropy‐related materials. This study investigates the importance of the diverse bonding and chemical environments when discussing multicomponent carbide materials. A combination of ab initio calculations and X‐ray photoelectron spectroscopy (XPS) was used to investigate the electronic structure of multicomponent thin films based on the (HfNbTiVZr)C system. The charge transfer was quantified theoretically using relaxed and nonrelaxed multicomponent as well as binary carbide reference structures, employing a fixed sphere model. High‐resolution XPS spectra from (HfNbTiVZr)C magnetron‐sputtered thin films displayed core‐level binding energy shifts and broadening effects as a result of the complex chemical environment. Charge transfer effects and a changed electronic structure in the multicomponent material, compared with the reference binary carbides, are observed both experimentally and in the density functional theory (DFT) simulations. The observed effects loosely follow electronegativity considerations, leading to a deviation from an ideal solid solution structure assuming nondistinguishable chemically equivalent environments.

The x-ray absorption spectrum of the tert-butyl radical: An experimental and computational investigation.

We report the x-ray absorption spectrum (XAS) of the tert-butyl radical, C4H9. The radical was generated pyrolytically from azo-tert-butane, and the XAS of the pure radical was obtained by subtraction of spectra recorded at different temperatures. The bands in the XAS were assigned by abinitio calculations that are in very good agreement with the experimental data. The lowest energy signal in the XAS is assigned to the C1s electron transition from the central carbon atom to the singly occupied molecular orbital (SOMO), while higher transitions correspond to C1s excitations from terminal carbon atoms. Furthermore, we investigated the fragmentation of the radical following resonant C1s excitation by electron-ion-coincidence spectroscopy. Several fragmentation channels were identified. The C1s excitation of the terminal carbons is associated with a stronger fragmentation tendency compared to the lowest C1s excitation of the central carbon into the SOMO. For this core excited state, we still observe an intact parent ion, C4H9+, and a comparatively higher tendency to dissociate into CH3+ + C3H6+.

Electrochemical Stability of MXenes in Water Based on Constant Potential AIMD Simulations.

MXene has been recently explored as promising electrocatalytic materials to accelerate the electrocatalytic process for hydrogen evolution, but their dynamic stability under electrochemical conditions remains elusive. Here we performed first-principle ab initio molecular dynamics calculations to reveal the electrochemical stability of Ti2CTx MXene in different aqueous environments. The results revealed the high vulnerability of the pure and vacancy-defected Ti2CO2 MXene towards water attack, leading to surface oxidation of MXene under neutral electrochemical condition that formed adsorbed oxygen species to Ti and dissociated proton in solution. The surface oxidation of Ti2CO2 could be prevented in the acid condition or in the neutral condition under the negative potential. Differently, the fully F- or OH-functionalized Ti2CF2 and Ti2C(OH)2 as well as the mixed functionalized Ti2C(O0.5OH0.5)2 and Ti2CO1.12F0.88 are highly stable under various electrochemical conditions, which can effectively prevent close contact between water and surface Ti atoms via electronic repulsion or steric hindrance. These findings provide atomic level understanding of the aqueous stability of MXene and provide useful strategies to prevent degradation and achieve highly stable MXenes.

Isomeric effects on the electron impact scattering for selected five membered heterocyclic ring molecules

We present electron scattering cross-sectional data for selected isomers of the five-membered ring molecules, specifically oxazole and isoxazole, thiazole and isothiazole, and imidazole and pyrazole. The ab-initio R-matrix method, incorporating static exchange polarization approximations, is employed for this calculation from energy range 0.1 eV to 20 eV. Three shape resonances are identified and characterized in each system, consisting of two π* and one σ* shape resonances. Notably, thiazole and isothiazole exhibit additional σ* resonance which is which is at a lower energy than other σ* resonance. The calculated resonance positions align well with available experimental as well as theoretical data. We also performed electronic structure calculations to aid in characterization of resonances. Comparison of the scattering cross-sectional data with their respective isomers reveal marginal differences in the magnitude of the elastic cross section particularly below 1 eV. This discrepancy may be attributed to variations in the long-range electronic dipole contribution to the electron-molecule interaction process. Further, to understand dependency of dipole moment on elastic cross section, we performed a fitting procedure for the elastic cross sections at 0.1 eV with square of the dipole moment of all the five membered ring molecules studied here. The fitting formula so obtained was used to estimate the cross section for other three five membered ring molecules: Furan, Thiophene and Pyrrole. Additionally, total ionization cross sections are computed using the binary-encounter-Bethe (BEB) model, demonstrating nearly identical results when compared with their isomers. Our findings offer insights into the influence of structural changes on the electron scattering data, providing guidance for other theoretical investigations.

Enhancement of superconductivity and structural phase transitions under pressure in FeTe0.89S0.11

The ternary iron chalcogenide, FeTe0.89S0.11 is a prominent member of the family of Fe-based superconductors with an ambient pressure Tc of around 8 K and a simple structure formed by layers of edge-sharing distorted tetrahedra separated by a van der Waals gap. Upon compression to 6 GPa, the Tc increases monotonically to 11 K. The superconductivity disappeared at about 6 GPa which was correlated to a softening of a Raman mode. High-resolution synchrotron X-ray diffraction shows that the FeTe0.89S0.11 started to transform from the monoclinic phase described by space group Cm to a disordered triclinic phase described by space group P1 at around 10 GPa (theoretical) - 14 GPa (experimental). Due to the phase transition the studied material transforms from a layered compound to a polymer, which is accompanied by a volume collapse of 4.65 %, indicating a discontinuous first-order phase transition. The disappearance of Raman modes, along with pressure independent resistivity, and ab initio calculations results support the structural and electronics phase transition. According to our results Fe(Te,S) may be useful for future high-field power applications being formed of relatively inexpensive and nontoxic elements.

Comparative analysis of KXH3(X= Mg, be) hydride cubic perovskites for hydrogen storage properties: A computational approach

The burgeoning potential of hydrogen as a future clean energy source has spurred intense research into hydrogen storage solutions. Novel hydride perovskite materials aside metal hydride are at the forefront of this effort, engaging researchers worldwide. This study utilizes Ab-initio calculations based on density functional theory (DFT) implemented in VASP to compare the electrical, elastic, thermodynamic, and optical properties of cubic KXH3 (X = Mg, Be) compounds. Generalized gradient approximation (GGA) with the Perdew-Burke-Ernzrhof (PBE) scheme is employed for optimization. Our calculations agree well with previously reported results. However, the tolerance factor suggests KBeH3 may not be stable in the cubic form. The electronic structure reveals KMgH3 to be semiconductive, while KBeH3 exhibits metallic behavior. Both compounds exhibit mechanical stability within the study. Optical properties are also investigated, showing suitability for hydrogen storage in the low energy range. Additionally, the mechanical parameters satisfy the Born stability criterion, indicating good mechanical stability for both compounds. Finally, the gravimetric ratios for hydrogen storage are calculated as 5.86 wt% and 4.52 wt% for KMgH3 and KBeH3, respectively. These high values suggest both compounds hold significant promise for diverse renewable energy applications and long-term hydrogen fuel storage.

Ab initio framework for deciphering trade-off relationships in multi-component alloys

While first-principles methods have been successfully applied to characterize individual properties of multi-principal element alloys (MPEA), their use in searching for optimal trade-offs between competing properties is hampered by high computational demands. In this work, we present a framework to explore Pareto-optimal compositions by integrating advanced ab initio-based techniques into a Bayesian multi-objective optimization workflow, complemented by a simple analytical model providing straightforward analysis of trends. We benchmark the framework by applying it to solid solution strengthening and ductility of refractory MPEAs, with the parameters of the strengthening and ductility models being efficiently computed using a combination of the coherent-potential approximation method, accounting for finite-temperature effects, and actively-learned moment-tensor potentials parameterized with ab initio data. Properties obtained from ab initio calculations are subsequently used to extend predictions of all relevant material properties to a large class of refractory alloys with the help of the analytical model validated by the data and relying on a few element-specific parameters and universal functions that describe bonding between elements. Our findings offer crucial insights into the traditional strength-vs-ductility dilemma of refractory MPEAs. The proposed framework is versatile and can be extended to other materials and properties of interest, enabling a predictive and tractable high-throughput screening of Pareto-optimal MPEAs over the entire composition space.

Enhancing the Electrochemical Activity of 2D Materials Edges through Oriented Electric Fields.

The edges of 2D materials have emerged as promising electrochemical catalyst systems, yet their performance still lags behind that of noble metals. Here, we demonstrate the potential of oriented electric fields (OEFs) to enhance the electrochemical activity of 2D materials edges. By atomically engineering the edge of a fluorographene/graphene/MoS2 heterojunction nanoribbon, strong and localized OEFs were realized as confirmed by simulations and spatially resolved spectroscopy. The observed fringing OEF results in an enhancement of the heterogeneous charge transfer rate between the edge and the electrolyte by 2 orders of magnitude according to impedance spectroscopy. Ab initio calculations indicate a field-induced decrease in the reactant adsorption energy as the origin of this improvement. We apply the OEF-enhanced edge reactivity to hydrogen evolution reactions (HER) and observe a significantly enhanced electrochemical performance, as evidenced by a 30% decrease in Tafel slope and a 3-fold enhanced turnover frequency. Our findings demonstrate the potential of OEFs for tailoring the catalytic properties of 2D material edges toward future complex reactions.