Ab Initio Density Functional Theory Research Articles

Perfluorinated Amines: Accelerating Lithium Electrodeposition by Tailoring Interfacial Structure and Modulated Solvation for High-Performance Batteries.

Modulating interfacial electrochemistry represents a prevalent approach for mitigating lithium dendrite growth and enhancing battery performance. Nevertheless, while most additives exhibit inhibitory characteristics, the accelerating effects on interfacial electrochemistry have garnered limited attention. In this work, perfluoromorpholine (PFM) with facilitated kinetics is utilized to preferentially adsorb on the lithium metal interface. The PFM molecules disrupt the solvation structure of Li+ and enhance the migration of Li+. Combined with the benzotrifluoride, a synergistic acceleration-inhibition system is formed. The ab initio molecular dynamics (AIMD) and density functional theory (DFT) calculation of the loose outer solvation clusters and the key adsorption-deposition step supports the fast diffusion and stable interface electrochemistry with an accelerated filling mode with C─F and C─H groups. The approach induces the uniform lithium deposition. Excellent cycling performance is achieved in Li||Li symmetric cells, and even after 200 cycles in Li||NCM811 full cells, 80% of the capacity is retained. This work elucidates the accelerated electrochemical processes at the interface and expands the design strategies of acceleration fluorinated additives for lithium metal batteries.

Comparison of lattice thermal conductivity using ab-initio DFT, machine learning interatomic potentials, and temperature dependent effective potential: a case study of hexagonal BN and BP bilayer

Discovering high thermal conductivity materials is essential for various practical applications, particularly in electronic cooling. The significance of two-dimensional (2D) materials lies in their unique properties that emerge due to their reduced dimensionality, making them highly promising for a wide range of applications. Hexagonal boron nitride (BN), both monolayer and bilayer forms, has garnered attention for its fascinating properties. In this work, we focus on bilayer boron phosphide (BP), which is isostructural to its BN analogue. The lattice thermal conductivity of both bilayer BN and BP have been calculated using ab-initio density functional theory, machine learning with the moment tensor potential method, and the temperature-dependent effective-potential method (TDEP). The TDEP approach gives more accurate results for both BN and BP materials. The lattice thermal conductivity of bilayer BP is lower than that of bilayer BN at room temperature, attributed to increased phonon anharmonicity. This study highlights the importance of understanding phonon scattering mechanisms in determining the thermal conductivity of 2D materials, contributing to the broader understanding and potential applications of these materials in future technologies.

Quantum Simulations and Experimental Insights into Glyphosate Adsorption Using Graphene-Based Nanomaterials.

The increasing global demand for food and agrarian development brings to light a dual issue concerning the use of substances that are crucial for increasing productivity yet can be harmful to human health and the environment when misused. Herein, we combine insights from high-level quantum simulations and experimental findings to elucidate the fundamental physicochemical mechanisms behind developing graphene-based nanomaterials for the adsorption of emerging contaminants, with a specific focus on pesticide glyphosate (GLY). We conducted a comprehensive theoretical and experimental investigation of graphene-based supports as promising candidates for detecting, sensing, capturing, and removing GLY applications. By combining ab initio molecular dynamics and density functional theory calculations, we explored several chemical environments encountered by GLY during its interaction with graphene-based substrates, including pristine and punctual defect regions. Our results unveiled distinct interaction behaviors: physisorption in pristine and doped graphene regions, chemisorption leading to molecular dissociation in vacancy-type defect regions, and complex transformations involving the capture of N and O atoms from impurity-adsorbed graphene, resulting in the formation of new GLY-derived compounds. The theoretical findings were substantiated by FTIR and Raman spectroscopy, which proposed a mechanism explaining GLY adsorption in graphene-based nanomaterials. The comprehensive evaluation of adsorption energies and associated properties provides valuable insights into the intricate nature of these interactions, shedding light on potential applications and guiding future experimental investigations of graphene-based nanofilters for water decontamination.

High pressure melt line of nickel using a generalized embedded atomic method potential

As the second most abundant metal in the Earth's core, nickel plays an important role in determining the structure and temperature of the Earth's core. Yet, the melt line of Ni at pressures corresponding to the Earth's core has not been explored in the literature. Many previous experimental and simulation efforts have reported the melting point of Ni at pressures below 100 GPa, but there exist large discrepancies, most of which have persisted due to various experimental and simulation bottlenecks in handling extreme pressure and temperature conditions. We adopted the generalized embedded atom method, which overcomes the limitations of existing interatomic potentials, to probe phase stability and phase boundaries of Ni at pressures between 50 and 500 GPa. The potential was validated by comparing the cold curves, phonon dispersion curves, and enthalpies of fusion with ab initio density functional theory calculations. Our analysis shows that face centered cubic (FCC) is stable, and the hexagonal close packed (HCP) and body centered cubic (BCC) phases are metastable close to the melt line. Melting temperatures at different pressures were obtained from two-phase co-existence simulations and take the following functional form: Tm=1969.23+19.15P−0.012P2. In contrast to iron, differences between the melting points of the stable and metastable phases of Ni are less than 250 K at 300 GPa, and the difference in melting points of the metastable BCC and HCP phases changes sign at 500 GPa, which implies that the phase transition mechanisms during solidification can be very complex.

Impact of Rh, Ru, and Pd Leads and Contact Topologies on Performance of WSe2 FETs: A First Comparative Ab Initio Study.

2D field-effect transistors (FETs) fabricated with transition metal dichalcogenide (TMD) materials are a potential replacement for the silicon-based CMOS. However, the lack of advancement in p-type contact is also a key factor hindering TMD-based CMOS applications. The less investigated path towards improving electrical characteristics based on contact geometries with low contact resistance (RC) has also been established. Moreover, finding contact metals to reduce the RC is indeed one of the significant challenges in achieving the above goal. Our research provides the first comparative analysis of the three contact configurations for a WSe2 monolayer with different noble metals (Rh, Ru, and Pd) by employing ab initio density functional theory (DFT) and non-equilibrium Green's function (NEGF) methods. From the perspective of the contact topologies, the RC and minimum subthreshold slope (SSMIN) of all the conventional edge contacts are outperformed by the novel non-van der Waals (vdW) sandwich contacts. These non-vdW sandwich contacts reveal that their RC values are below 50 Ω∙μm, attributed to the narrow Schottky barrier widths (SBWs) and low Schottky barrier heights (SBHs). Not only are the RC values dramatically reduced by such novel contacts, but the SSMIN values are lower than 68 mV/dec. The new proposal offers the lowest RC and SSMIN, irrespective of the contact metals. Further considering the metal leads, the WSe2/Rh FETs based on the non-vdW sandwich contacts show a meager RC value of 33 Ω∙μm and an exceptional SSMIN of 63 mV/dec. The two calculated results present the smallest-ever values reported in our study, indicating that the non-vdW sandwich contacts with Rh leads can attain the best-case scenario. In contrast, the symmetric convex edge contacts with Pd leads cause the worst-case degradation, yielding an RC value of 213 Ω∙μm and an SSMIN value of 95 mV/dec. While all the WSe2/Ru FETs exhibit medium performances, the minimal shift in the transfer curves is interestingly advantageous to the circuit operation. Conclusively, the low-RC performances and the desirable SSMIN values are a combination of the contact geometries and metal leads. This innovation, achieved through noble metal leads in conjunction with the novel contact configurations, paves the way for a TMD-based CMOS with ultra-low RC and rapid switching speeds.

Chemical Environment and Structural Variations in High Entropy Oxide Thin Film Probed with Electron Microscopy.

We employ analytical transmission electron microscopy (TEM) to correlate the structural and chemical environment variations within a stacked epitaxial thin film of the high entropy oxide (HEO) Mg0.2Co0.2Ni0.2Cu0.2Zn0.2O (J14), with two layers grown at different substrate temperatures (500 and 200 °C) using pulsed laser deposition (PLD). Electron diffraction and atomically resolved STEM imaging reveal the difference in out-of-plane lattice parameters in the stacked thin film, which is further quantified on a larger scale using four-dimensional STEM (4D-STEM). In the layer deposited at a lower temperature, electron energy loss spectroscopy (EELS) mapping indicates drastic changes in the oxidation states and bonding environment for Co ions, and energy-dispersive X-ray spectroscopy (EDX) mapping detects more significant cation deficiency. Ab initio density functional theory (DFT) calculations validate that vacancies on the cation sublattice of J14 result in significant electronic and structural changes. The experimental and computational analyses indicate that low temperatures during film growth result in cation deficiency, an altered chemical environment, and reduced lattice parameters while maintaining a single phase. Our results demonstrate that the complex correlation of configurational entropy, kinetics, and thermodynamics can be utilized for accessing a range of metastable configurations in HEO materials without altering cation proportions, enabling further engineering of functional properties of HEO materials.

Novel vacancy-ordered RbKGeCl6 and RbKGeBr6 double perovskites for optoelectronic and thermoelectric applications: an ab-initio DFT study

The present study examines the key characteristics of new vacancy-ordered halide double perovskites, RbKGeCl6 and RbKGeBr6, encompassing the elastic, structural, mechanical, optoelectronic, and thermoelectric properties. The Density Functional Theory (DFT) was employed to perform the calculation of the properties, facilitating the evaluation of their potential applications in optoelectronic and thermoelectric devices. The DFT calculation was conducted using the Quantum Espresso package alongside the thermo_pw tool and the BoltzTraP codes. The results revealed that the two proposed compounds possess both chemical and mechanical stability with optimized lattice constants recorded at 10.14 Å and 10.72 Å for RbKGeCl6 and RbKGeBr6, respectively. The evaluation of the elastic properties of the materials suggested reasonably high mechanical moduli of the materials. Based on the calculated electronic properties, the materials are classified as direct gap semiconductors, with energy gap values of 2.11 eV for RbKGeCl6 and 0.80 eV for RbKGeBr6 using the GGA-PBE functional. Furthermore, the use of the SCAN approximation yields more reliable energy gap of 2.51 eV and 1.08 eV for the respective compounds. The materials exhibited a high absorption coefficient and a significantly low reflectivity within the visible-ultraviolet energy spectrum. These findings strongly suggest the promising properties of the materials under study for optoelectronic applications. Furthermore, the calculated thermoelectric properties of the materials, particularly the figure of merit, revealed the materials’ potential use as thermoelectric materials. The calculated figure of merit values of RbKGeCl6 and RbKGeBr6 were found to range from 0.73 to 0.75, respectively, between 300 K and 800 K. Despite being lower, these values are comparable to those of some well-established materials including SiGe alloys (0.95), Bi2Te3 (≈0.90), and PbTe (≈0.80).

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.

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.0 mol L–1 HCl using different methods. Gravimetric experiments have shown that CMN and CMP reach efficiencies of 80 and 91% at 1.45 mmol L–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.

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

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.