Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Non-Electricity Based Renewable Fuels: Theory and Computation for Solar Thermochemical Hydrogen

Dominated by photovoltaics and wind, current renewable energy sources generate mostly electricity, but 80% of the global final energy consumption occurs in form of fuels. Therefore, direct solar fuel generation would be a major breakthrough for the energy transition. Solar thermochemical hydrogen (STCH) is one of the very few potential routes towards scalable renewable fuels, but currently suffers from lack of an oxide working material that could optimally perform energy conversion within the thermodynamic boundary conditions. Theory and computation can contribute in two distinct ways, through materials search and discovery, but also by providing detailed mechanistic models for specific systems so to advance our understanding of possible design strategies. To enable high-throughput materials screening, we developed a defect graph neural network (dGNN) machine learning approach,[1] which accelerates the prediction of defect formation energies by replacing the tedious density functional theory (DFT) supercell calculations for all possible defect sites. This approach enables high-throughput database screening of oxides, which was integrated with thermodynamic modeling to extract the reduction entropies as additional selection criterion for STCH. Once potential candidate materials are identified, detailed models can guide materials design by predicting performance characteristics. One challenge is to quantitatively predict thermochemical equilibria at high concentrations when the redox active defects start to interact with each other, thereby impeding the formation of additional defects. Introducing a model for the free energy of defect interaction, parametrized on the basis of DFT data, we simulated the complete STCH redox cycle for (Sr,Ce)MnO3 alloys, achieving near-quantitative agreement with experimental data.[2] The analysis of these simulations reveals how defect interactions diminish the reduction entropy and H2 yield, suggesting to include these interactions in design considerations. Finally, we revisit the popular van't Hoff method for analyzing reduction enthalpies and entropies. This method is not ideal, as it involves a temperature-dependent convolution of gas-phase and solid-state entropies, causing uncertainties in the same order of magnitude as the physical quantities of interest. To avoid this problem, we suggest a simple alternative approach which can be applied to experimental and simulated data alike.

first-principles calculations↗
DOE OSTI2025

Single-Walled Carbon Nanotubes with Confined Chalcogens as the Catalysts and Electrodes for Oxygen Reduction Reaction in Fuel Cells

The goal of this project is to synthesize and characterize a new non-metal electro-catalyst for oxygen reduction reaction (ORR) for fuel cell applications. The intended catalyst is a composite material composed of sulfur chains encapsulated in narrow diameter single-walled carbon nanotubes (S@SWNTs). S@SWNTs were successfully synthesized through sulfur vapor infusion method, and validated with Raman spectroscopy. However, the electrochemical analysis on the ORR catalytic activity of S@SWNTs indicated that it had low ORR catalytic activity. Our theoretical study based on density functional theory (DFT) revealed that the poor oxygen adsorption (low binding energy) on the surface of S@SWNTs was the bottleneck of the entire catalytic reaction. The focus of the project was subsequently pivoted to the development of a new non-noble metal ORR catalyst that could provide durability in acidic electrolyte. The hypothesis was to encapsulate small iron (Fe) clusters in SWNTs (Fe@SWNTs) can provide ORR electro-catalytic activity and long durability in acidic environment. DFT-based computational studies were carried out to explore the feasibility of the Fe@SWNTs catalyst. The theoretical study indicated that Fe@SWNTs indeed could catalyze the ORR with lower theoretical overpotential than platinum (Pt). However, its weaker bonding energy with oxygen was the bottleneck of the overall reaction. On the other hand, Fe clusters (composed of 7 Fe atoms) encapsulated in nitrogen-doped SWNTs (Fe 7 @N 4 WSNTs) showed proper oxygen adsorption by the Fe cluster and low theoretical overpotential of ORR. Comparing to the Fe single atom catalyst on N-doped SWNTs (Fe-N 4 SWNTs), which is one the best non-noble metal ORR catalysts reported in the literature, the Fe 7 @N 4 WSNTs showed lower overpotential and better resistance to acidic environment. Fe encapsulated N-doped SWNTs were synthesized with ferrocene as the Fe precursors through vapor infusion method, and experimental validate is underway. This study theoretically demonstrated the feasibility of a new type of non-noble metal electro-catalyst for ORR that could have high catalytic activity and long durability.

08 HYDROGEN↗
DOE OSTI2025

Intrinsic Defect-Induced Local Semiconducting-to-Metallic Regions Within Monolayer 1T-TiS2 Displayed by First-Principles Calculations and Scanning Tunneling Microscopy

Using density functional theory (DFT) and scanning tunneling microscopy (STM), the intrinsic point defects, formation energy, and electronic structure of 1T-TiS2 were investigated. Defect systems include single-atom vacancies, interstitial and adatom additions, and direct atomic substitution. Using a collective approach for analyzing realistic systems for point defect investigation, we provide a more straightforward comparison to the experimental measurements, reproducing more realistic environmental conditions related to thin film growth. STM images are compared to computationally simulated electron density images to identify specific geometries that result from favorable point defects. DFT suggests that titanium interstitials are the most energetically favorable intrinsic defect, and sulfur vacancies are more likely to form than titanium vacancies within this realistic analysis, which is in agreement with STM data. A pristine, stoichiometric monolayer system is calculated to have a direct band gap of 0.422 eV, which varies based on local point defects. Local semiconducting-to-metallic electronic transitions are predicted to occur based on the presence of Ti interstitials.

Keeney, P. J.↗
DOE OSTI2025

Electronic density of states as the descriptor of elastic bond strength, ductility, and local lattice distortion in BCC refractory alloys

Although electronic density of states (DOS) is fundamental to materials properties, its general relationship to mechanical properties of alloys is not well established. In this paper, using density functional theory (DFT) calculations, we show that the electronic occupancy at the Fermi level, N(E f ), obtained from DOS is a key descriptor of alloy strength and ductility. Our comprehensive analysis of numerous body centered cubic (BCC) refractory high entropy alloys (RHEAs) shows an overwhelming correlation that low N(E f ) indicates strong bonds that have high stiffness resulting in high elastic constants. High bond stiffness indicates presence of covalent nature of bonds that are directional in nature resulting in resistance to deformation leading to high bulk (B) and shear (G) moduli. Consequently, N(E f ) provides a direct correlation to the tendency of alloy ductility evidenced in the Pugh ratio (G/B). As stiffer bonds result in lower local lattice distortion (LLD), N(E f ) are LLD are also found to be corelated which opens up a correlation to solid solution strengthening and yield strength. Thus, this work unveils fundamental correlations between N(E f ) and (1) elastic bond strength, (2) ductility, and (3) LLD. These correlations open opportunities for the design of high strength high ductile RHEAs.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Polyphase-discrete Fourier transform spectrum analysis for the Search for Extraterrestrial Intelligence sky survey

The sensitivity of a matched filter-detection system to a finite-duration continuous wave (CW) tone is compared with the sensitivities of a windowed discrete Fourier transform (DFT) system and an ideal bandpass filter-bank system. These comparisons are made in the context of the NASA Search for Extraterrestrial Intelligence (SETI) microwave observing project (MOP) sky survey. A review of the theory of polyphase-DFT filter banks and its relationship to the well-known windowed-DFT process is presented. The polyphase-DFT system approximates the ideal bandpass filter bank by using as few as eight filter taps per polyphase branch. An improvement in sensitivity of approx. 3 dB over a windowed-DFT system can be obtained by using the polyphase-DFT approach. Sidelobe rejection of the polyphase-DFT system is vastly superior to the windowed-DFT system, thereby improving its performance in the presence of radio frequency interference (RFI).

Zimmerman, G. A.↗
DOE OSTI2025

Investigation of the Effect of Framework Flexibility on CO 2 Adsorption in SIFSIX-3-Cu Using a Machine-Learned Force Field

Metal–organic frameworks (MOFs) offer promise as selective CO 2 sorbents, but successful MOF sorbent materials need high CO 2 binding affinity and selectivity for CO 2 over water. This work focuses on the use of machine-learned force fields (MLFFs) to model CO 2 adsorption in flexible MOFs, with a focus on SIFSIX-3-Cu, an anion-pillared MOF known for its high CO 2 affinity. A preliminary high-throughput screening of over 900 anion-pillared MOFs was performed using rigid UFF+DDEC6 force fields to predict zero-loading heats of adsorption for CO 2 and H 2 O. SIFSIX-3-Cu was selected for further computational study due to its predicted CO 2 heat of adsorption and experimental relevance. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies and forces, with an iterative sampling scheme combining molecular dynamics, geometry optimization, random geometric insertion, and NVT Monte Carlo-based configuration generation to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included, contrasting with previous models that approximated the MOF as rigid. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations with the MLFF produced CO 2 adsorption isotherms in good agreement with experimental data at direct air capture (DAC) pressures (e.g., 40 Pa), in contrast to previous overestimations of CO 2 sorption by models with rigid structures. Bond and angle histogram analysis showed that MOF flexibility increased the variance of fluorine–fluorine diagonal distances at adsorption sites, resulting in a lower predicted sorption for flexible, asymmetric SIFSIX-3-Cu pore geometries compared to the rigid, symmetric DFT-optimized SIFSIX-3-Cu pore geometry. A detailed description of flexibility afforded by the MLFF resulted in an accurately predicted CO 2 uptake (0.88 mmol/g) at low pressure (40 Pa) compared to the experimentally measured value (1.24 mmol/g). In conclusion, these results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low-pressure applications.

adsorption↗
DOE OSTI2026

How Does Metal Spin State Affect Electronic Communication in Mixed-Valence Dimers? Insights from Ultrafast Near-Infrared and Soft X-ray Transient Absorption Spectroscopy

Recent advancements in photocatalysis, photovoltaics, and quantum information science take advantage of electron spin, and determining how spin multiplicity affects electron transfer is key to understanding these phenomena. Here, in this study, we examine how metal spin state affects electronic communication in an organometallic mixed-valence dimer, ferrocenyl cobaltocenium ([Fe II Cp 2 Co III Cp 2 ] + ). This complex can be photoexcited from its low-spin singlet Fe II ground state to form intermediate-spin triplet Fe II and high-spin quintet Fe II excited states. Using femtosecond optical transient absorption (OTA) spectroscopy with visible (vis), near-infrared (NIR), and short-wave IR (SWIR) probes, supported by time-dependent density functional theory (TD-DFT) calculations, we measure Fe II Co III →Fe III Co II intervalence charge transfer (IVCT) bands in each of the Fe II spin states. Mulliken–Hush analysis of the excited-state IVCT bands was used to compute the electronic coupling between the metal centers in all three spin states, which increased as quintet < triplet < singlet. Meanwhile, the peak energy of the bands, and thus the ΔG of the IVCT transition, increased as triplet < quintet < singlet. Then, to directly probe the electronic structure at each metal center, we employed picosecond soft X-ray transient absorption (XTA) spectroscopy at the Fe and Co L 3 edges. Our results show that the low-spin and high-spin states of [Fe II Cp 2 Co III Cp 2 ] + are valence-localized, while the intermediate-spin state is partially delocalized. The differences in charge delocalization are attributed to differences in orbital occupation and geometry that affect the free energy and superexchange coupling.

Burke, John H. [Univ. of Illinois at Urbana-Champa↗
DOE OSTI2025

Decoding the Desorption Mechanism of 2LiH:1Mg(NH2)2 Using Metal Borohydrides

The complex metal hydride 2LiH:1Mg­(NH2)2 has emerged as a promising material for stationary hydrogen storage applications, such as seasonal storage or energy backup systems, due to its high volumetric and gravimetric capacities and robust reversibility. However, its widespread adoption is hindered by sluggish reaction rates, performance degradation upon cycling, and improper end-use cases. To address these problems and better understand the desorption pathway, we used metal borohydrides (MBH4; M = Li–Cs) as chemical probes. A thorough analysis of the bulk behavior of all six materials, including hydrogen cycling experiments, X-ray absorption spectroscopy, FTIR, pXRD, solid-state NMR, and ab initio DFT simulations, shows that the borohydride additives decrease the activation energy of hydrogen release by about 20 kJ/mol for MBH4@2:1 materials versus pristine. Furthermore, more surface-sensitive studies show that the amide-to-imide desorption pathway in these materials, while essentially complete in the bulk, is incomplete in the near-surface region, suggesting an “inverse core–shell” desorption mechanism for amide dehydrogenation to imide. The kinetic enhancements produced by MBH4 additives (M = K, Rb, and Cs) are attributed to the destabilization of the amide N–H bond and interaction with the LiH/Mg­(NH2)2 interface to promote H–H bond formation. An inverse core-shell mechanism is also operative in the hydrogen desorption for the 2LiH:1LiNH2 system, suggesting this may be a general feature of amides. Given the fast dehydrogenation rate and large gravimetric capacity, these materials satisfy these requirements for telecom backups and seasonal microgrid storage applications.

Absorption↗
DOE OSTI2025

Multifunctional Polar 2D Lead Iodide Perovskites Exhibiting Persistent Spin Texture

Ferroelectric Rashba semiconductors are a rare class of multifunctional materials promising for spin–orbitronics due to the possibility of electrically switchable spin textures. Hybrid organic–inorganic metal halide perovskites offer exceptional tunability that can be harnessed to target noncentrosymmetry; however, the complex interactions between organic and inorganic components are not rationally understood. Here, we use an asymmetric spacer cation with a strong molecular dipole moment (2-fluorobenzylammonium, 2FBZ) and increase the intrinsic quantum well thickness to n > 1 to synthesize three new polar 2D lead iodide perovskites (2FBZ) 2 (A) n-1 Pb n I 3n+1 with A = methylammonium (MA) or formamidinium (FA) and n = 2, 3. Single-crystal structure analysis reveals the C 2v symmetry of these crystal structures and substantial structural distortions, which enable Rashba and Dresselhaus band splitting and persistent spin texture, confirmed by DFT calculations. Measurements of the resulting symmetry-dependent properties, such as second harmonic generation, switchable photovoltaic effect, ferroelectric polarization, and low temperature photoluminescence show switchable ferroelectric semiconducting properties. Furthermore, these results introduce a general chemical design strategy toward polar symmetry by exploiting the complex interplay between asymmetric spacer cations and A-site cations in quasi-2D (n > 1) halide perovskites to realize new classes of multifunctional materials for future spin–orbitronic applications.

Cations↗
DOE OSTI2024

Fermi surface topology and magnetotransport properties of superconducting Pd 3 Bi 2 Se 2

Pd 3 Bi 2 Se 2 is a rare realization of a superconducting metal with a non-zero topological invariant. Here, in this study, we report the growth of high-quality single crystals of layered Pd 3 Bi 2 Se 2 with a superconducting transition at T c ≈ 0.80 K and upper critical fields of ~10 mT and ~5 mT for the in plane and out-of-plane directions, respectively. Our density functional theory (DFT) calculations reveal three pairs of doubly degenerate bands crossing the Fermi level all displaying clear three dimensional dispersion consistent with the overall low electronic anisotropy (<2). The multiband electronic nature of Pd 3 Bi 2 Se 2 is evident in magneto-transport measurements, yielding a sign changing Hall resistivity at low temperatures. The magnetoresistance is non-saturating and follows Kohler’s scaling rule. We interpret the magneto-transport data in terms of open orbits that are revealed in the DFT calculated Fermi surface. de Haas-van Alphen (dHvA) oscillation measurements using torque magnetometry on single crystals yield four frequencies for out-of-plane fields: F α = (150 ± 26) T, F β = (293 ± 10) T, F γ = (375 ± 20) T and F η = (1017 ± 12) T, with the low frequency dominating the spectrum. Through the measurement of angular dependent dHvA oscillations and DFT calculations we identify the F α frequency with an approximately ellipsoidal electron pocket centered on the L 2 point of the Brillouin zone. Lifshitz-Kosevich analysis of the dHvA oscillations reveals a small cyclotron effective mass m* = (0.11 ± 0.02)m 0 and a nontrivial Berry phase for the dominant orbit. The presence of nontrivial topology in a bulk superconductor positions Pd 3 Bi 2 Se 2 as a potential candidate for exploring topological superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Revealing the electronic structure of van der Waals antiferromagnetic NiPS 3 through synchrotron-based 𝜇-ARPES and alkali metal dosing

Antiferromagnetic NiPS 3 has recently emerged as a quantum material of considerable interest, thanks to the discovery of multiple new couplings involving electrons, spins, orbitals, phonons, and magnons. However, controversies and open questions persist concerning the fundamental origins of these couplings. A critical piece of information required to advance the understanding is the precise electronic band structure of NiPS 3 . Angle-resolved photoemission spectroscopy (ARPES), combined with alkali metal dosing (AMD), can enable us to directly observe the subtle electronic states that appear around the Fermi surface, offering valuable insights into the intriguing quantum properties and interplays of the examined material. Here, in this study, we present a comprehensive characterization and analysis of the band structure of van der Waals layered antiferromagnet NiPS 3 , leveraging state-of-the-art μ-ARPES measurements supported by density functional theory (DFT) calculations. Theoretical DFT results identify the orbital contributions to the observed bands, providing a precise understanding of the experimental ARPES data. Crucially, AMD enables the observation of conduction band and defect-related states above the valence band maximum in NiPS 3 . Furthermore, temperature dependent ARPES results across the Néel transition temperature of NiPS 3 reveal that the paramagnetic and antiferromagnetic phases have nearly identical band structures, underlining the highly localized character of Ni d states. These findings substantially deepen our understanding of the electronic properties of NiPS 3 and lay a vital foundation for exploring the intriguing quantum phenomena it exhibits.

Cao, Yifeng [Boston Univ., MA (United States); Law↗
NASA NTRS1991

The behavior of quantization spectra as a function of signal-to-noise ratio

An expression for the spectrum of quantization error in a discrete-time system whose input is a sinusoid plus white Gaussian noise is derived. This quantization spectrum consists of two components: a white-noise floor and spurious harmonics. The dithering effect of the input Gaussian noise in both components of the spectrum is considered. Quantitative results in a discrete Fourier transform (DFT) example show the behavior of spurious harmonics as a function of the signal-to-noise ratio (SNR). These results have strong implications for digital reception and signal analysis systems. At low SNRs, spurious harmonics decay exponentially on a log-log scale, and the resulting spectrum is white. As the SNR increases, the spurious harmonics figure prominently in the output spectrum. A useful expression is given that roughly bounds the magnitude of a spurious harmonic as a function of the SNR.

Flanagan, M. J.↗
DOE OSTI2025

Hydrogen and water interactions with CrMnFeCoNi alloy from density functional theory calculations

High entropy alloys (HEAs) are a promising class of materials with remarkable mechanical and catalytic properties. Among these, the quinary CrMnFeCoNi alloy (also called “Cantor alloy”) has attracted considerable attention given its thermodynamic stability and remarkable mechanical properties under different temperatures. Given that various degradation mechanisms involve multiple contaminants, such as hydrogen and water in hydrogen embrittlement and surface poisoning, respectively, understanding their interactions with the Cantor alloy is critical for its practical applications as structural, nuclear, or hydrogen storage material. In this work, we perform first-principles calculations based on Density Functional Theory (DFT) to investigate such interactions when considering various microstructures, including bulk materials and those containing certain defects, such as grain boundaries, stacking faults, and vacancies. We also employ Global Sensitivity Analysis to identify the importance of different factors in the stability of the impurities. We find that the accuracy of the H formation energy is significantly affected by spin polarization and chemical short-range order. The study also identifies a strong tendency for hydrogen interstitials to segregate to Σ5(210)/[001] symmetric tilt grain boundary, even when H concentrations are high, suggesting that a certain type of grain boundaries acts as H sinks within the alloy. Further, this result is reinforced by the low formation energy of vacancy-hydrogen complexes, which can contain multiple hydrogen atoms. Finally, the surface reactivity analysis reveals that the adsorption energy of oxygen and hydroxyl groups is highly sensitive to the specific metal atom involved in the binding, with a clear preference for chromium atoms, which could have implications for the alloy’s oxidation and corrosion behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigation of the iodate sorption mechanism by CoAl LDH through experiments and ab initio molecular dynamics simulations

Radioiodine released during the nuclear-fuel cycle constitutes a persistent radiological hazard. In this study, the IO 3 – uptake mechanism of CoAl LDH was resolved by combining pH-controlled sorption experiments, synchrotron XAFS, and DFT-based AIMD simulations. At pH close to 6, approximately 90% of IO 3 – was removed, and the equilibrium distribution coefficient reached about 1.7 × 10 4 mL g –1 . EXAFS analysis indicated an average iodine–oxygen bond length of 1.81 Å and a coordination number near 3, with the fit R-factor equal to 0.002. The simulations faithfully reproduced the experimental spectrum and revealed transient proton hopping events that generated metastable I–O–H species inside the interlayer, thereby confirming nitrate-to-iodate exchange as the controlling capture pathway. Atomic density profiles and radial distribution functions further showed that IO 3 – adopt an end-on orientation perpendicular to the hydroxide sheets, while water molecules mediate proton migration without disturbing the host lattice. In conclusion, the integrated experimental–computational evidence demonstrates that CoAl LDH can rapidly and selectively sequester IO 3 – under near-neutral conditions, offering atomic scale guidance for the rational engineering of layered sorbents for advanced radioactive-waste treatment.

Kang, Jaehyuk [Jeju National Univ. (Korea, Republi↗
DOE OSTI2025

Interplay between Cation “Coloring” and Stereochemically Active Lone Pairs in AgBiS 2 Thin Films

Solution-processed AgBiS 2 thin films were fabricated using novel thiol− amine precursor inks to investigate the stereochemical activity of Bi 3+ 6s 2 lone pairs and their impact on the structure. A dual-space analysis combining Bragg diffraction and hard X-ray photoelectron spectroscopy (HAXPES) revealed a rock salt-like average structure with local distortions linked to cation coloring. Density functional theory (DFT) and crystal orbital Hamilton population (COHP) analyses confirmed that local Bi-rich and Ag-rich nanodomains amplify stereochemical activity, whereas more mixed and cation-order nanodomains are less stereochemically active. This local, nanoscopic mixing of segregated and ordered domains would indeed explain an average Fm3̅m structure that is rock salt-like and that does not manifest the full anharmonicity and noncentrosymmetry evidenced in canonical structures with stereochemical expression. These findings provide insights into the local structural and electronic complexities governing the optoelectronic properties of AgBiS 2 thin films.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

A Solid State Zwitterionic Plastic Crystal with High Static Dielectric Constant

The dielectric data in Figure 3, Figure 4, Figure S6 of the published paper was extracted from 2EOIMTSA-BDS-DATA .txt file. This file can be directly opened using a text file editor. It can also be imported to Excel/ Origin for further plotting and analysis. The G' and G'' in Figure 3 of the publihsed paper was plotted from data in file 2EOImTSA-temperature-sweep.xlsx. This file can be directly opend using Excel. The details of DFT simulations mentioned in Figure 2, Figure 7, and Figure S9 of the published paper are included in the DFT.zip file.

Huang, Zitan [Pennsylvania State University]↗
DOE OSTI2025

From Micro-environments to Macroscopic Effects: How the Alkaline Hydrogen Evolution Reaction Drives Cu Cathodic Corrosion

Cathodic corrosion of copper (Cu) has posed a significant challenge for over a century, impeding various technological progresses such as electrochemical conversion of CO 2 (eCO 2 RR) into fuels and other value-added carbon products. Here, in this study, employing a combined Density Functional Theory (DFT) and kinetic Monte Carlo (kMC) simulation approach, we delve into the atomistic level mechanism driving this phenomenon in Cu. Our hypothesis posits the pivotal role of alkaline hydrogen evolution reaction (HER) in facilitating cathodic corrosion in Cu. We rigorously develop a pH-dependent hydroxide (OH) adsorption mechanism and calculate the equilibrium OH coverage (𝜃 OH ) at varying pH levels, the thermodynamic stability of subsurface oxygen (O sub ), as well as the Cu-vacancy mediated diffusion of subsurface oxygens (O sub ). Through comprehensive analysis, we establish correlation among various microenvironments, including oxygen diffusion in subsurface layers, pH-dependent OH adsorption, and Cu dissolution into the electrolyte as (Cu-OH) complexes. Furthermore, our investigation explores the correlation between surface coordination environment of active sites and cathodic corrosion of Cu. Finally, by integrating DFT-derived thermodynamic data into a kMC model, we successfully predict the formation of experimentally observed corrosion pits on Cu-surfaces. This combined approach not only advances our fundamental understanding of Cu cathodic corrosion but also offers insights crucial for developing effective corrosion mitigation strategies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗