Ab Initio Density Functional Theory Research Articles

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.

A Blyholder mechanism in the chemisorption of N2O on Ni(111) – studied with Auger-photoelectron coincidence spectroscopy

In heterogeneous catalysis the surface-adsorbate bond strength is critical for the function of the system. Here we study a series consisting of multilayer, bilayer and monolayer N2O on Ni(111) and employ Auger-photoelectron coincidence spectroscopy (APECS) to study the interaction between the molecule and the substrate directly. We observe intensity in the nitrogen Auger spectra that arise from the interaction between molecule and surface (not observed in free molecules) whereas the oxygen spectra are thickness-independent. Since the two nitrogen atoms of N2O are chemically inequivalent we can assign the intensity present in the bilayer and monolayer cases to orbitals centered on the terminal nitrogen which is closest to the Ni(111) surface. Using ab initio, molecular dynamics and solid-state density functional theory calculations we infer a Blyholder model of the surface bond as consisting of donation from the terminal nitrogen lone-pair valence orbital with back-donation from the metal into the unoccupied orbitals on that nitrogen. This coincidence technique can readily be used to study substrate–adsorbate interactions directly with chemical and orbital specificity — this opens up prospects to study fundamental steps of molecular adsorption and heterogeneous catalysis with unprecedented detail.

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

We present a comprehensive study of transport coefficients including DC electrical conductivity and related optical properties, electrical contribution to the thermal conductivity, and the shear viscosity via ab initio molecular dynamics and density functional theory calculations on the “priority 1” cases from the “Second Charged-Particle Transport Coefficient Workshop” [Stanek et al., Phys. Plasmas (to be published 2024)]. The purpose of this work is to carefully document the entire workflow used to generate our reported transport coefficients, up to and including our definitions of finite size and statistical convergence, extrapolation techniques, and choice of thermodynamic ensembles. In pursuit of accurate optical properties, we also present a novel, simple, and highly accurate algorithm for evaluating the Kramers–Kronig relations. These heuristics are often not discussed in the literature, and it is hoped that this work will facilitate the reproducibility of our data.

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.

MoS2 Field-Effect Transistor Performance Enhancement by Contact Doping and Defect Passivation via Fluorine Ions and Its Cyclic Field-Assisted Activation.

MoS2-based field-effect transistors (FETs) and, in general, transition metal dichalcogenide channels are fundamentally limited by high contact resistance (RC) and intrinsic defects, which results in low drive current and lower carrier mobilities, respectively. This work addresses these issues using a technique based on CF4 plasma treatment in the contacts and further cyclic field-assisted drift and activation of the fluorine ions (F-), which get introduced into the contact region during the CF4 plasma treatment. The F- ions are activated using cyclic pulses applied across the source-drain (S/D) contacts, which leads to their migration to the contact edges via the channel. Further, using ab initio molecular dynamics and density functional theory simulations, these F- ions are found to bond at sulfur (S) vacancies, resulting in their passivation and n-type doping in the channel and near the S/D contacts. An increase in doping results in the narrowing of the Schottky barrier width and a reduction in RC by ∼90%. Additionally, the passivation of S vacancies in the channel enhances the mobility of the FET by ∼150%. The CF4 plasma treatment in contacts and further cyclic field-assisted activation of F- ions resulted in an ON-current (ION) improvement by ∼90% and ∼480% for exfoliated and CVD-grown MoS2, respectively. Moreover, this improvement in ION has been achieved without any deterioration in the ION/IOFF, which was found to be >7-8 orders.

Structure and molecular interactions in 1‑ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ion pair conformers: abinitio DFT based study

Density functional theory (DFT) calculations with Grimme’s B97D functional including dispersion are employed to provide both structural and electronic insight into the multifold interactions occurring in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic liquid [EMI][TFSI], which is currently being targeted for applications in next-generation Li-ion battery electrolytes. The geometric structure, interaction energy and natural bond orbital analysis of the ion-pair conformers of [EMI][TFSI] show differences from the conventional H-bonds. A total of six stable [EMI][TFSI] ion pair conformers were observed. Our results show that a doubly ionic H-bond exists between [EMI]+ and [TFSI]− ion pair conformers. The greater the number of multiple H-bond type interactions, the greater the absolute value of the interaction energy. Moreover, results from the NBO analysis show that, for [EMI][TFSI] ion pair conformers, charge transfer occurs mainly from the lone pairs of oxygen and nitrogen atom to the σ-type anti-bonding orbital of the C–H and π-type anti-bonding orbitals of N–C bonds. This was evident from the values of the stabilization energy E(2) associated with each electron delocalization from the donor to acceptor orbitals. There exist large numbers of multiple concomitant hydrogen bonds, but the values of the stabilization energy E(2) are generally small (E(2)n→\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ o$$\\end{document}σ* < 3 kcal/mol) for the individual H-bond interactions. The shorter the C–H···O and/or C–H···N bond is, the larger charge transfer, and the larger the stabilization energy E(2) associated with electron delocalization from donor to acceptor. Generally, the [EMI][TFSI] ion pair conformers tend to form multiple but bent H bonds, reducing the strength of the individual H bonds from a potential (linear) maximum.

Novel Polytype of III-VI Metal Chalcogenides Nano Crystals Realized in Epitaxially Grown InTe.

III-VI metal chalcogenides have garnered considerable research attention as a novel group of layered van der Waals materials because of their exceptional physical properties and potential technological applications. Here, the epitaxial growth and stacking sequences of InTe is reported, an essential and intriguing material from III-VI metal chalcogenides. Aberration-corrected scanning transmission electron microscopy (STEM) is utilized to directly reveal the interlayer stacking modes and atomic structure, leading to a discussion of a new polytype. Furthermore, correlations between the stacking sequences and interlayer distances are substantiated by atomic-resolution STEM analysis, which offers evidence for strong interlayer coupling of the new polytype. It is proposed that layer-by-layer deposition is responsible for the formation of the unconventional stacking order, which is supported by ab initio density functional theory calculations. The results thus establish molecular beam epitaxy as a viable approach for synthesizing novel polytypes. The experimental validation of the InTe polytype here expands the family of materials in the III-VI metal chalcogenides while suggesting the possibility of new stacking sequences for known materials in this system.

Theoretical study of the complex antiferromagnetic Heusler alloys Pd[formula omitted]MnAl, Pd[formula omitted]MnIn, and Pt[formula omitted]MnAl

The complex antiferromagnetic (AFM) phases displayed by the Heuser alloys Pd2MnAl, Pd2MnIn, and Pt2MnAl, are theoretically studied. The crystalline structure, electronic, and magnetic properties are calculated within the ab initio density functional theory. The super unit cell considered to describe the various magnetic arrangements contains 128 atoms. We analyzed the ferromagnetic (FM) and the three different antiferromagnetic (AFM) phases and found that the Pd2MnAl crystalline structure is orthorhombic with dimensions 6.055, 6.461, 6.026 Å. The AFM state with the lowest energy is the AFM3 with a Mn magnetic moment of 3.813 μB. The Pd2MnIn crystallizes in an orthorhombic structure with lattice parameters 6.249, 6.763, 6.227 Å. The ground state antiferromagnetic arrangement is also of type 3 with a Mn magnetic moment of 3.932 μB. The Pt alloy crystallizes in an orthorhombic structure with lattice dimensions of 6.029, 6.906, and 5.683 Å. In this case, the antiferromagnetic ground state is of type 1 with a Mn magnetic moment of 3.716 μB. The Pd and Pt d states are deep in energy and do not interact much with the Mn d states. This produces, in Pd and Pt, no magnetic polarization.