Ab Initio Density Functional Theory Research Articles

Dissociative Adsorption of O2 on Ag3Au(111) Surface: A Density Functional Theory Study

The catalytic efficiency of oxygen reduction catalysts is notably influenced by the dissociative adsorption of O2. We conducted a systematic investigation into the dissociative adsorption of O2 on the Ag3Au(111) surface using ab initio density functional theory (DFT) calculations. Our computational findings indicate that adsorption the configuration designated t-b-t exhibits favorable energetics on the Ag3Au(111) surface. Regarding the dissociation of O2, we identified a reasonable dissociation pathway, which proceeds from the initial t-b-t state to the creation of two oxygen atoms that occupy a set of neighboring fcc sites. Furthermore, our analysis indicates that the adsorption of O2 on the Ag3Au(111) surface is less favored thermodynamically and more difficult to dissociate than that on the Ag(111) surface. This study furnishes a theoretical framework elucidating the prospective utilization of Ag-Au alloy in the capacity of oxygen reduction catalysts.

Computational Investigation of LiF Formation at Graphite-Electrolyte Interfaces.

The performance of rechargeable batteries is strongly influenced by the solid-electrolyte interphase (SEI), and a comprehensive understanding of SEI formation from the atomic level is crucial for effective battery design. The dynamics of the electrode-electrolyte interface is important and needs to be considered when evaluating the mechanism of the SEI formation. Here, we employed ab initio molecular dynamics (AIMD) and density functional theory (DFT) calculations to examine interfacial behaviors and LiF formation. Through molecular dynamics and structure sampling, we successfully constructed an electrochemical stability diagram correlating the thermodynamic free energy with the potential, which is determined by the work function of electrode surfaces. DFT calculations revealed that LiF formation at the graphite-electrolyte interfaces occurs easily via the intermediate LiHF complex. Interestingly, LiF tends to be solvated by solvents rather than directly deposited onto electrode surfaces (e.g., the Au electrode), a phenomenon we identify as a critical determinant of the porous and uneven nature of the LiF layer observed on graphite electrodes. Our finding offers new mechanistic insights into LiF formation at graphite-electrolyte interfaces.

Surface-Enriched Room-Temperature Liquid Bismuth for Catalytic CO2 Reduction.

Bismuth-based electrocatalysts are effective for carbon dioxide (CO2) reduction to formate. However, at room temperature, these materials are only available in solid state, which inevitably suffers from surface deactivation, declining current densities, and Faradaic efficiencies. Here, the formation of a liquid bismuth catalyst on the liquid gallium surface at ambient conditions is shown as its exceptional performance in the electrochemical reduction of CO2 (i.e., CO2RR). By doping a trace amount of bismuth (740ppm atomic) in gallium liquid metal, a surface enrichment of bismuth by over 400 times (30 at%) in liquid state is obtained without atomic aggregation, achieving 98% Faradic efficiency for CO2 conversion to formate over 80 h. Ab initio molecular simulations and density functional theory calculations reveal that bismuth atoms in the liquid state are the most energetically favorable sites for the CO2RR intermediates, superior to solid Bi-sites, as well as joint GaBi-sites. This study opens an avenue for fabricating high-performing liquid-state metallic catalysts that cannot be reached by elementary metals under electrocatalytic conditions.

Quinoline derivatives as corrosion inhibitors for mild steel in acid medium: Deeper insights from experimental, ab initio density functional theory modeling, and in silico ecotoxicity investigations

Two (E)-2-(((2-((7-chloroquinolin-4-yl)amino)ethyl)imino)methyl)-4-R-phenol Schiff bases derivatives CMN (R = NO2) and CMP (R = H) were successfully prepared. Then, CMN and CMP were evaluated as corrosion inhibitors for 1020 mild steel in 1.0molL–1 HCl using different methods. Gravimetric experiments have shown that CMN and CMP reach efficiencies of 80 and 91% at 1.45mmolL–1. Atomic Force Microscopy (AFM) and Scanning Electron Microscopy (SEM) measurements have indicated the formation of a protective adsorbed film on the mild steel surface. Electrochemical Impedance Spectroscopy (EIS), Linear Polarization Resistance (LPR), Potentiodynamic Polarization (PP), and Electrochemical Frequency Modulation (EFM) have also shed some light on electrochemical behavior of compounds. These data have depicted better polarization resistance and lower corrosion current density in the presence of both organic molecules, also proving that the process is controlled by charge transfer and that CMN and CMP are mixed-type corrosion inhibitors. ab initio Density Functional Theory (DFT) investigated how two quinoline-based molecules (CMN and CMP) prevent corrosion on iron surfaces. Simulations revealed CMP bonds more strongly to iron than CMN, leading to better protection. Simulations also showed that these inhibitor molecules form a protective layer on the iron, slowing down the movement of corrosive substances. Environment risk assessment through QSAR-based models have indicated that the corrosion inhibitors do not have the potential to bioaccumulate in fish. Still, the acute and chronic toxicity values at the different trophic levels investigated call for attention for environmentally friendly structural optimization.

Structural characteristics of irrational Type-II Twin interfaces

The Type-II Twin Boundary (TB) is a critical interface in functional materials whose irrational Miller-index identity has recently drawn significant research interest. This study establishes general structural characteristics of the Type-II twin interface, utilizing TBs in Shape Memory Alloys (SMAs) - TiPd, TiPt, and AuCd - as study targets. It is shown how the irrational identity of each TB is explained by the Terrace-Disconnection (T-D) structural topology. It is proposed that the terrace is the rational-plane nearest to the irrational TB in the reciprocal space, having integral Miller-indices of least magnitude. Crystallographic-registry on this terrace requires non-trivial coherence-strains. A novel kinematic-origin of the coherence-strain is proposed, coming directly from a transformation of the classical twinning deformation-gradient. This transformation revealed that the classical twinning-shear partitions into the coherence-strain and a new metric termed the “terrace-shear”. It is shown that the magnitude of shear relating the twin-structure to the matrix is the terrace-shear and not the twinning-shear, contrary to classical understanding. Furthermore, the Burgers vector of the twinning disconnection is shown to be related directly to the terrace-shear. The energy of each Type-II interface is determined from ab initio Density Functional Theory (DFT) calculations. It is shown that the energy-minimal atomic-structure on the terrace requires determination of a “lattice-offset” that is non-trivial and unknown apriori. In summary, this study expounds on T-D topological structure of Type-II twin interfaces, establishing methods to identify rational terraces, coherence strains, ab initio planar TB energies and revealing a unique partitioning of the twinning-shear exhibited by this class of interfaces.

ZrN doped with 3d transition metals: a computational investigation of main physical properties for high-temperature structural applications

In the present study, ab initio density functional theory calculations were used to assess the effect of first-row transition metals (X = Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn) on the stability of Zr0.5N0.5 nitrides. Specifically, the structural, mechanical, and electronic properties were studied to evaluate their applicability in high-temperature structural applications such as coating. The heat of formation for all X-doped Zr0.5N0.5 ternaries were found to be lower than that of the undoped Zr0.5N0.5. Specifically, Mn-doped Zr0.5N0.5 was observed to be the most thermodynamically stable structure, due to its lowest heat of formation. The density of states for both the undoped and doped Zr0.5N0.5 nitrides indicated full metallic behavior and observed that doping with 3d-transition metals reduce the density of states at the Fermi energy, thereby enhancing the electronic stability. Furthermore, mechanical stability was observed in these nitrides with increased melting temperatures expect for Zr0.5N0.5 doped Ti. Since Zr0.5N0.5 doped with X is thermodynamically, electronically, and mechanically stable, they are deemed suitable for high-temperature structural applications especially Zr0.5N0.5 doped Mn.

Synthesis and characterization of ceramic high entropy carbide thin films from the Cr-Hf-Mo-Ta-W refractory metal system

We use reactive DC magnetron sputtering to showcase synthesis strategies for multicomponent carbides with the NaCl-type fcc structure and illustrate how deposition conditions allow controlling the formation of metallic and ceramic single phases in the Cr-Hf-Mo-Ta-W system. The synthesis is performed in argon flow and different acetylene flows from 0 to 12sccm, at ambient and elevated temperatures (700 °C), respectively, hindering/promoting the adatom diffusion. Structural and microstructural investigations reveal the formation of the bcc metallic phase (a=3.188−3.209Å) in films deposited without acetylene flow, also supported by ab initio density function theory (DFT) analysis of lattice parameters as a function of the C content. Experimentally, a bcc-to-fcc phase transition is observed through the formation of an amorphous coating. Contrarily, samples deposited in higher acetylene flow show an fcc multielement carbide phase (a=4.33−4.49 Å). The crystalline films reveal columnar morphology, while the amorphous ones are very dense. We report promising mechanical properties, with hardness up to 25±1GPa. The indentation moduli reach up to 319±6GPa and show trends consistent with DFT predictions. Our study paves the path towards the preparation of Cr-Hf-Mo-Ta-W multicomponent carbides by magnetron sputtering, showing promising microstructure as well as mechanical properties.

Violation of the Cauchy–Born rule in multi-principal element alloys

Multi-principal element alloys are a novel class of materials that are formed by combining multiple elements in high concentrations and show exceptional properties compared to conventional alloys. These alloys have high configurational entropy due to inherent atomic disorder. The Cauchy–Born rule is a popular homogenization method for linking atomistic to continuum length scales. In this Letter, we use ab initio density functional theory calculations to report that the Cauchy–Born rule, which holds in conventional alloys in the absence of defects, is not valid in multi-principal element alloys. The violation of the Cauchy–Born rule in these alloys is attributed to the presence of atomic disorder due to which the deformations are inhomogeneous. Our results also show the deviations in deformation are related to the magnitude of the stretch and shear by power laws.

Direct Evidence of Reversible Changes in Electrolyte and its Interplay with LiO2 Intermediate in Li-O2 Batteries.

Lithium-oxygen batteries show promising energy storage potential with high theoretical energy density; however, further investigation of chemical reactions is required. In this study, experimental Raman and theoretical analyzes are performed for a Li-O2 battery with LiClO4/dimethyl sulfoxide (DMSO) electrolyte and carbon cathode to understand the role of intermediate species in the reactional mechanism of the cell using a high donor number solvent. Operando Raman results reveal reversible changes in the DMSO bands, in addition to the formation and decomposition of Li2O2. On discharge, a decrease in DMSO polarizability is observed and bands of DMSO-Li+-anion interactions are evidenced and supported by abinitio density functional theory (DFT) calculations. Molecular dynamics (MD) force field simulations and operando Raman show that DMSO interacts with LiO2(sol), highlighting the stability of the electrolyte compared to the interaction with reactive . On charging, the presence of Li+ indicates the formation of a lithium-deficient phase, followed by the release of Li+ and oxygen. Therefore, this study contributes to understanding the discharge/charge chemistry of a Li-O2 cell, employing a common carbon cathode and DMSO electrolyte. The combination of a simple characterization technique in operando mode and theoretical studies provides essential information on the mechanism of Li-O2system.

Transport coefficients of warm dense matter from Kohn-Sham density functional theory

Transport coefficients of warm dense matter from Kohn-Sham density functional theory

Hydrogen diffusion in BCC-Fe: DFT study of tensorial stress effects and interactions with point defects

Hydrogen embrittlement is a multifaceted phenomenon that can significantly compromise the toughness of susceptible metals. We study the critical process of hydrogen diffusion within the body-centered cubic lattice of iron (BCC-Fe or α−Fe) using ab-initio density functional theory (DFT). Beyond standard investigations on the effect of hydrostatic stress, we extend our study to incorporate the influence of uniaxial and pure shear stress states. We find substantial alterations in diffusion barriers. Moreover, we study scenarios involving paired point (zero-dimensional) defects, including one vacancy and an additional interstitial hydrogen atom. Calculations comprehend various configurations and transitions between them, providing insight into the interplay between the chemical environment and mechanical fields. The barrier values here determined offer essential data for complementary techniques, including Monte-Carlo models, to accurately describe hydrogen diffusion. This research contributes to a deeper understanding of hydrogen embrittlement, offering insight that can inform the design of safer materials used in structural engineering.

Towards atomistic modelling of solid Pb-O formation and dissolution in liquid lead coolant: Interatomic potential development

Microscopic description of solid Pb-O formation and dissolution in liquid lead coolant is important for the modelling of fast neutron reactors. For this purpose, in this work we develop interatomic potential models for lead melt with dissolved oxygen. Potentials fitting is based on ab initio density functional theory (DFT) calculations and quantum molecular dynamics. Two different potential models are trained on the DFT data using the force-matching method: a computationally cheap embedded atom method (EAM) potential and a computationally expensive machine-learned moment tensor potential (MTP) that allows reproducing DFT data very accurately. The potentials are further validated to reproduce the properties of lead melt and the coordination of oxygen atoms in molten lead. The diffusion coefficient for oxygen in molten lead is calculated and compared with the available experimental data. We analyse how these interatomic potential models describe the crystal structure of the solid lead oxide to consider PbO nucleation and dissolution in molten lead. We extensively discuss the arising problems with the description of long-range electrostatic interactions. The possibility of matching the experimental data on the oxygen solubility limit by tuning the EAM potential is shown.

Dzyaloshinskii-Moriya interactions, Néel skyrmions and V4 magnetic clusters in multiferroic lacunar spinel GaV4S8

Using ab initio density functional theory with static mean-field correlations, we calculate the Heisenberg and Dzyaloshinskii-Moriya interactions (DMI) for an atomistic spin Hamiltonian for the lacunar spinel, GaV4S8. The parameters describing these interactions are used in atomistic spin dynamics and micromagnetic simulations. The magnetic properties of the lacunar spinel GaV4S8, a material well-known from experiment to host magnetic skyrmions of Néel character, are simulated with these ab initio calculated parameters. The Dzyaloshinskii-Moriya contribution to the micromagnetic energy is a sum of two Lifshitz invariants, supporting the formation of Néel skyrmions and its symmetry agrees with what is usually expected for C3ν-symmetric systems. There are several conclusions one may draw from this work. One concerns the quantum nature of the magnetism, where we show that the precise magnetic state of the V4 cluster is crucial for understanding quantitatively the magnetic phase diagram. In particular, we demonstrate that a distributed-moment state of each V4 cluster explains well a variety of properties of GaV4S8, such as the band gap, observed Curie temperature and especially the stability of Néel skyrmions in the experimentally relevant temperature and magnetic-field range. In addition, we find that electronic correlations change visibly the calculated value of the DMI.

Ab-initio DFT＋[formula omitted] study on the electronic correlation in multiferroic XY-antiferromagnetic Ba2CoGe2O7

Ab-initio density functional theory (DFT)＋U has been used to study the spin-resolved electronic band structure and the partial density of states of the multiferroic XY-antiferromagnetic Ba2CoGe2O7. Using the constrained random-phase approximation (cRPA) method, the value of Hubbard U and Hund’s coupling JH has been determined and plugged into the further dynamical mean field theory (DMFT) calculations. The strong metal-ligand hybridization between Co-3d and O-2p orbitals in the CoO4 tetrahedron has been demonstrated via DFT and DFT＋DMFT calculations. The calculated indirect band gap in the antiferromagnetic state without and with U (=5 eV) correction is found to be 0.67 and 2.71 eV, respectively. The highest calculated ordered magnetic moment of 2.7 μB/Co2+ is close to the reported experimental value. The presented DFT calculations based on the different orbital models brings the underlying strong electronic correlation in the Ba2CoGe2O7.

Physical origin of enhanced electrical conduction in aluminum-graphene composites

The electronic and transport properties of aluminum-graphene composite materials were investigated using the ab initio plane wave density functional theory. The interfacial structure is reported for several configurations. In some cases, the face-centered aluminum (111) surface relaxes in a nearly ideal registry with graphene, resulting in a remarkably continuous interface structure. The Kubo–Greenwood formula and space-projected conductivity were employed to study electronic conduction in aluminum single- and double-layer graphene-aluminum composite models. The electronic density of states at the Fermi level is enhanced by the graphene for certain aluminum–graphene interfaces, thus improving electronic conductivity. In double-layer graphene composites, conductivity varies non-monotonically with temperature, showing an increase between 300 and 400 K at short aluminum-graphene distances, unlike the consistent decrease in single-layer composites.

Study of the Crystal Structure and Hydrogen Bonding during Cold Crystallization of Poly(trimethylene 2,5-furandicarboxylate).

Here, we present a detailed description of the in situ isothermal crystallization of poly(trimethylene 2,5-furandicarboxylate)(PTF) as revealed by real-time Fourier transform infrared spectroscopy (FTIR) and grazing incidence wide-angle X-ray scattering (GIWAXS). From FTIR experiments, the evolution of hydrogen bonding with crystallization time can be monitored in real time, while from GIWAXS, crystal formation can be followed. Density functional theory (DFT) calculations have been used to simulate FTIR spectra for different theoretical structures, enabling a precise band assignment. In addition, based on DFT ab initio calculations, the influence of hydrogen bonding on the evolution with crystallization time can be understood. Moreover, from DFT calculations and comparison with both FTIR and GIWAXS experiments, a crystalline structure of poly(trimethylene 2,5-furandicarboxylate) is proposed. Our results demonstrate that hydrogen bonding is present in both the crystalline and the amorphous phases and its rearrangement can be considered as a significant driving force for crystallization of poly(alkylene 2,5-furanoate)s.

Characterizing ABS–copper chemistry-dependent adhesion: From the atomic to macro level

The electroplating of copper layers on acrylonitrile–butadiene–styrene (ABS) polymer surfaces is a common need in industrial applications. The development of new eco-friendly techniques to promote ABS–copper adhesion is crucial due to the harsh and expensive chemicals used in the existing processes. The initial path involves establishing a better understanding on chemistry-dependent ABS–copper adhesion mechanisms. By formulating the material with more functional groups, the adhesion between ABS and the coated copper layer was characterized through a multiscale approach. At the atomic scale, we performed ab initio density functional theory simulations to obtain the variations in the adsorption energy. At the nano and micro scales, we adopted atomic force microscopy to measure the adhesion forces on the bulk polymer surfaces. At the macro scale, we compared these adhesion results with both 90-degree peeling tests and finite element model analyses. This study presents the first comprehensive approach for bridging information gaps across the scales (from nano to macro) with computational and experimental data, facilitating the design of formulated ABS additives.

Visualizing the Low-Energy Electronic Structure of Prototypical Hybrid Halide Perovskite through Clear Band Measurements.

Organic-inorganic hybrid perovskites (OIHPs) are a promising class of materials that rival conventional semiconductors in various optoelectronic applications. However, unraveling the precise nature of their low-energy electronic structures continues to pose a significant challenge, primarily due to the absence of clear band measurements. Here, we investigate the low-energy electronic structure of CH3NH3PbI3 (MAPI3) using angle-resolved photoelectron spectroscopy combined with ab initio density functional theory. We successfully visualize the electronic structure of MAPI3 near the bulk valence band maximum by using a laboratory photon source (He Iα, 21.2 eV) at low temperature and explore its fundamental properties. The observed valence band exhibits a highly isotropic and parabolic band characterized by small effective masses of 0.20-0.21 me, without notable spectral signatures associated with a large polaron or the Rashba effect, subjects that are intensely debated in the literature. Concurrently, our spin-resolved measurements directly disprove the giant Rashba scenario previously suggested in a similar perovskite compound by establishing an upper limit for the Rashba parameter (αR) of 0.28 eV Å. Our results unveil the unusually complex nature of the low-energy electronic structure of OIHPs, thereby advancing our fundamental understanding of this important class of materials.

Influence of temperature on bandgap shifts, optical properties and photovoltaic parameters of GaAs/AlAs and GaAs/AlSb p–n heterojunctions: insights from ab-initio DFT + NEGF studies

The III–V group semiconductors are highly promising absorbers for heterojunctions based solar cell devices due to their high conversion efficiency. In this work, we explore the solar cell properties and the role of electron–phonon coupling (EPC) on the solar cell parameters of GaAs/AlSb and GaAs/AlAs p–n heterojunctions using non-equilibrium Green function method (NEGF) in combination of ab-initio density functional theory (DFT). In addition, the band offsets at the heterointerfaces, optical absorption and bandgap shifts (BGSs) due to temperature are estimated using DFT + NEGF approach. The interface band gaps in heterostructures are found to be lower than bulk band gaps leading to a shift in optical absorption coefficient towards lower energy side that results in stronger photocurrent. The temperature dependent electronic BGS is significantly influenced by the phonon density and phonon energy via EPC. The phonon influenced BGS is found to change the optical absorption, photocurrent density and open-circuit voltage. In case of GaAs/AlSb junction, the interface phonons are found to have significantly higher energies as compared to the bulk phonons and thereby may have important implications for photovoltaic (PV) properties. Overall, the present study reveals the influence of EPC on the optical absorption and PV properties of GaAs/AlSb and GaAs/AlSb p–n heterojunctions. Furthermore, the study shows that the DFT + NEGF method can be successfully used to obtain the reasonable quantitative estimates of temperature dependent BGSs, optical absorption and PV properties of p–n heterojunctions.

Multiscale Modeling of Magnetoelectric Nanoparticles for the Analysis of Spatially Selective Neural Stimulation.

The growing field of nanoscale neural stimulators offers a potential alternative to larger scale electrodes for brain stimulation. Nanoelectrodes made of magnetoelectric nanoparticles (MENPs) can provide an alternative to invasive electrodes for brain stimulation via magnetic-to-electric signal transduction. However, the magnetoelectric effect is a complex phenomenon and challenging to probe experimentally. Consequently, quantifying the stimulation voltage provided by MENPs is difficult, hindering precise regulation and control of neural stimulation and limiting their practical implementation as wireless nanoelectrodes. The work herein develops an approach to determine the stimulation voltage for MENPs in a finite element analysis (FEA) model. This model is informed by atomistic material properties from ab initio Density Functional Theory (DFT) calculations and supplemented by experimentally obtainable nanoscale parameters. This process overcomes the need for experimentally inaccessible characteristics for magnetoelectricity, and offers insights into the effect of the more manageable variables, such as the driving magnetic field. The model's voltage is compared to in vivo experimental data to assess its validity. With this, a predictable and controllable stimulation is simulated by MENPs, computationally substantiating their spatial selectivity. This work proposes a generalizable and accessible method for evaluating the stimulation capability of magnetoelectric nanostructures, facilitating their realization as wireless neural stimulators in the future.