SEARCH · Engineering Papers
Results for “density functional theory development”
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.
A Universal Machine Learning Model for Elemental Grain Boundary Energies
The grain boundary (GB) energy has a profound influence on the grain growth and properties of polycrystalline metals. Here, we show that the energy of a GB, normalized by the bulk cohesive energy, can be described purely by four geometric features. By machine learning on a large computed database of 361 small Σ (Σ<10) GBs of more than 50 metals, we develop a model that can predict the grain boundary energies to within a mean absolute error of 0.13 J m –2 . More importantly, this universal GB energy model can be extrapolated to the energies of high Σ GBs without loss in accuracy. These results highlight the importance of capturing fundamental scaling physics and domain knowledge in the design of interpretable, extrapolatable machine learning models for materials science.
Role of Fluorine in Chemomechanics of Cation-Disordered Rocksalt Cathodes
The role of nanoscale chemomechanical behavior in the macroscopic performance of functional materials is well recognized. For lithium-ion battery cathodes, tremendous effort has been devoted to the development of new chemistry and structure, yet fundamental understanding of the correlation between redox processes and mechanical properties of the novel materials lags behind. In the present study, we prepare large discrete single grains of Li-excess cation-disordered rocksalts (DRX) and investigate their chemomechanical behavior at the particle level, using nanoresolution X-ray and electron-based spectro-imaging and chemical mapping techniques. While irregular cracking upon lithium extraction leads to the eventual breakdown of the baseline DRX oxide (Li 1.2 Ti 0.4 Mn 0.4 O 2 ) particles at a high delithiation state, the fluorinated-DRX (Li 1.3 Ti 0.3 Mn 0.4 O 1.7 F 0.3 ) clearly displays aligned cracking along the <001> direction. The resulting periodicity in the cracking pattern enables the particles to retain their integrity and, consequently, improved electrochemical stability. Density functional theory (DFT) calculations showed that fluorination leads to increased concentration of Li+ on the (001) planes and preferential Li movements along the <001>-family directions, revealing the underlying mechanism for directional cracking. Overall, our study demonstrates the unique role of fluorine in modulating nanoscale chemomechanics, which in turn influences the evolution of charge and strain heterogeneity at the particle level. These insights provide important design guidelines in further improving DRX cathode materials.
Energy density functional and sensitivity of energies of giant resonances to bulk nuclear matter properties
The development of a modern and more realistic nuclear energy density functional (EDF) for accurate predictions of properties of nuclei is the subject of enhanced activity, since it is very important for the study of properties of nuclear matter (NM), giant resonances and, in particular, of properties of rare nuclei with unusual neutron-to-proton ratios. Here, we provide a short review of the current status of the nuclear EDF and the theoretical results obtained for properties of nuclei and nuclear matter. We will first describe a method for determining the parameters of the EDF, associated with the Skyrme type effective interaction, by carrying out a Hartree-Fock based fit to extensive set of data of ground state properties and constraints. Next we will describe the fully self-consistent Hartree-Fock (HF) based random-phase-approximation (RPA) theory for calculating the strength functions S(E) and centroid energies E CEN of giant resonances and the folding model (FM) distorted wave Вorn approximation (DWBA) to calculate the excitation cross section of giant resonances by α scattering. Then we will provide results for: (i) the Skyrme parameters of the KDE0v1 EDF; (ii) consequences of violation of self-consistently in HF-based RPA; (iii) FM-DWBA calculation of excitation cross section; (iv) values of the E CEN of isoscalar and isovector giant resonances of multipolarities L=0–3 for a wide range of spherical nuclei, using 33 EDFs associated with standard form of the Skyrme type interactions, commonly employed in the literature; and (v) the sensitivities E CEN of the giant resonances to bulk properties of NM. We also determine constraints on NM properties, such as the incompressibility coefficient and effective mass, by comparing with experimental data on E CEN of giant resonances.
Buffered Coordination Modulation as a Means of Controlling Crystal Morphology and Molecular Diffusion in an Anisotropic Metal-Organic Framework
We know significant advances have been made in the synthesis of chemically selective environments within metal-organic frameworks, yet materials development and industrial implementation have been hindered by the inability to predictively control crystallite size and shape. One common strategy to control crystal growth is the inclusion of coordination modulators, which are molecular species designed to compete with the linker for metal coordination during synthesis. However, these modulators can simultaneously alter the pH of the reaction solution, an effect that can also significantly influence crystal morphology. Herein, noncoordinating buffers are used to independently control reaction pH during metal-organic framework synthesis, enabling direct interrogation of the role of the coordinating species on crystal growth. We demonstrate the efficacy of this strategy in the synthesis of low-dispersity single-crystals of the framework Co 2 (dobdc) (dobdc 4- = 2,5-dioxido-1,4-benzenedicarboxylate) in a pH 7-buffered solution using cobalt(II) acetate as the metal source. Density functional theory calculations reveal that acetate competitively binds to Co during crystallization, and by using a series of cobalt(II) salts with carboxylate anions of varying coordination strength, it is possible to control crystal growth along the c -direction. Finally, we use zero length column chromatography to show that crystal morphology has a direct impact on guest diffusional path length for the industrially important hydrocarbon m -xylene. Together, these results provide molecular-level insight into the use of modulators in governing crystallite morphology and a powerful strategy for the control of molecular diffusion rates within metal-organic frameworks.
Synergistic Doping and Stabilization of Magnetically Tunable LnTi3(Sb,Sn)4 (Ln: Ce–Gd) Kagome Metals
Here, we present our synthesis and characterization of the LnTi3(Sb,Sn)4 (Ln: Ce, Pr, Nd, Sm, Gd) family of cleavable kagome metals. While these materials are isostructural to the LnTi3Bi4 family, they only form as (Sb,Sn) solid-solutions with no corresponding LnTi3Sb4 or LnTi3Sn4 phases. We use a combination of first-principles density functional theory (DFT) and crystal orbital Hamilton population (COHP) calculations to show that (Sb,Sn) alloying has a stabilizing effect on the structure by adjusting the Fermi level, filling bonding states, depopulating antibonding states, and adjusting the density of states (DOS) toward local minima, an effect we call “synergistic doping.” Through a detailed characterization of the SmTi3(Sb,Sn)4 series, we further demonstrate that the tunable Fermi level has a profound effect on the physical properties. We observe multiple magnetic phases that stem from a competition between antiferromagnetic (AFM) and ferromagnetic-like (FM) ground states. Furthermore, the (Sb,Sn) ratio allows us to tune from direct AFM-FM competition toward a complex admixture with properties reminiscent of the recently discovered TbTi3Bi4. Ultimately, our work demonstrates how the concept of synergistic doping provides a means to stabilize new structures while developing systems with intrinsically tunable chemical, magnetic, and electronic properties.
Highly dense atomic Fe–Ni dual metal sites for efficient CO 2 to CO electrolyzers at industrial current densities
Carbon-supported, atomically dispersed, nitrogen-coordinated metal sites (e.g., Fe and Ni) are arguably the most promising catalysts for the electrochemical reduction of CO 2 to CO due to their unique catalytic properties and the use of earth-abundant elements. However, conventional single metal sites are constrained by their structural simplicity, causing either too weak or too strong absorption/desorption of multiple critical intermediates (e.g., *COOH and *CO). Current catalysts also suffer from ultra-low loadings (<1.0 wt%) of atomic metal active sites in catalysts, leading to inadequate performance for CO 2 -to-CO conversion. Here, we develop dual Ni/Fe metal site catalysts with significantly increased atomically dispersed metal loadings (up to 4.8 wt%). A gas-phase chemical vapor deposition (CVD) approach to introducing single Ni sites was integrated with Fe 2 O 3 /ZIF-8 precursors, followed by an optimal thermal activation. The optimized CVD-Ni/Fe–N–C catalyst exhibited remarkable electrocatalytic performance for the CO 2 reduction to CO in a continuous membrane-electrode-assembly electrolyzer, achieving a maximum CO faradaic efficiency (FE CO ) of 96% at a current density of 700 mA cm −2 in a near-neutral electrolyte. Furthermore, a desirable but challenging acidic flow-cell electrolyzer was designed using this dual metal site catalyst to improve CO 2 utilization, accomplishing a FE CO of up to 95% at a CO partial current density close to 600 mA cm −2 . Density functional theory (DFT) calculations suggest a synergetic effect between Fe–Ni pairs facilitating *COOH intermediate formation and *CO desorption simultaneously during CO 2 to CO conversion. This is key to breaking the linear scaling relationship of conventional single-metal site catalysts during the CO 2 reduction reaction.
Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts
Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.
Synergic Effects of Surface Chemistry and Applied Potentials on the Kinetics of the Electrocatalytic Oxygen Evolution Reaction in IrO 2
Improving the efficacy of the oxygen evolution reaction (OER) through water oxidation is critical for advancing photoelectrochemical water splitting. Among the catalysts, IrO 2 exhibits a high OER catalytic activity and stability under acidic conditions. The OER mechanism in this system has been a topic of intense research; however, many mechanistic understandings are lacking. In this work, we carried out first-principles calculations with an implicit solvation model at a constant potential to study the OER behavior on the IrO 2 (110) surface. We find that the surface hydrogen coverage has a significant effect on OER kinetics and transition states. We then develop a microkinetics model that accounts for the continuous evolution of both OER activation energy and hydrogen coverage as a function of an applied potential. We show that this inclusion leads to significant improvement in the simulated Tafel plot compared to available experiments. In conclusion, our results point to a complex interplay between surface chemistry and the applied potential on OER kinetics.
Structure and Oxygen Evolution Activity of β-NiOOH: Where Are the Protons?
Ni oxides and oxyhydroxides (NiO x ) have been studied for a long time as cathode materials for alkaline batteries and electrocatalysts for the oxygen evolution reaction (OER). Yet, understanding of the connection between their atomic and electronic structures and electrochemical performance or stability is still incomplete. In this work, we use first-principles density functional theory (DFT) calculations to revisit the structure, electronic properties, and OER activity of β-NiOOH, the catalytically active phase of NiO x . Following extensive DFTbased screening, we identify a hitherto overlooked structure characterized by a uniform distribution of H atoms on the NiO 2 layers. All the Ni 3+ cations in this structure exhibit an identical t g 6 e g 1 electronic configuration with an occupied 3d z2 orbital. Comparison of the calculated bond lengths with extended X-ray absorption fine structure (EXAFS) data unequivocally supports this structure relative to all other low-energy configurations. Based on this structure, we uncover and detail defect-dominated OER mechanisms on the basal β-NiOOH (001) surface, with overpotentials as low as 0.39 V. The present results should provide a valuable contribution to ongoing efforts for understanding and developing enhanced transition-metal hydroxide catalysts for the OER.
Spin-informed universal graph neural networks for simulating magnetic ordering
The screening and discovery of magnetic materials are hindered by the computational cost of first-principles density-functional theory (DFT) calculations required to find the ground state magnetic ordering. Although universal machine-learning interatomic potentials (uMLIPs), also known as atomistic foundation models, offer high-fidelity models of many atomistic systems with significant speedup, they currently lack the inputs required for predicting magnetic ordering. In this work, we present a data-efficient, spin-informed graph neural network framework that incorporates spin degrees of freedom as inputs and preserves physical symmetries, extending the functionality of uMLIPs to simulate magnetic orderings. This framework speeds up DFT calculations through better initial guesses for magnetic moments, determines the ground-state ordering of bulk materials and even generalizes to magnetic ordering in surfaces. Furthermore, we implement a closed-loop anomaly detection approach that effectively addresses the classic "chicken-and-egg" problem of creating a high-quality dataset while developing a uMLIP, unearthing anomalies in large benchmark datasets and boosting model accuracy.
Effective Field Theory for jet substructure in heavy ion collisions
I develop an Effective Field Theory (EFT) framework to compute jet substructure observables for heavy ion collision experiments. As an example, I consider dijet events that accompany the formation of a weakly coupled long lived Quark Gluon Plasma (QGP) medium in a heavy ion collision and look at an observable insensitive to jet selection bias: the simultaneous measurement of jet mass along with the transverse momentum imbalance between the jets that are groomed to remove soft radiation. Treating the jet as an open quantum system, I write down a factorization formula within the SCET (Soft Collinear Effective Theory) framework in the forward scattering regime. The physics of the medium is encoded in a universal soft field correlator while the jet-medium interaction is captured by a medium induced jet function. The factorization formula leads to a Lindblad type equation for the evolution of the reduced density matrix of the jet in the Markovian approximation. The solution for this equation allows a resummation of large logarithms that arise due to the final state measurements imposed while simultaneously summing over multiple incoherent interactions of the jet with the medium.
Excited State Dynamics of a Conformationally Fluxional Copper Coordination Complex
The conversion of solar energy into chemical fuel represents a capstone goal of the 21st century and has the potential to supply terawatts of power in a globally distributed manner. However, the disparate time scales of photodriven charge separation (~fs) and steps in chemical reactions (~μs) represent an inherent bottleneck in solar-to-fuels technology. Here to address this discrepancy, we are developing earth-abundant coordination complexes that undergo light-induced conformational rearrangements such that charge separation (CS) is hastened, while charge recombination (CR) is slowed. To these ends, we report the preparation and characterization of a new series of conformationally fluxional copper coordination complexes that contain a twisted intramolecular charge transfer (TICT) fluorophore as part of their ligand scaffold. Structural and spectroscopic characterization of the Cu(I) and Cu(II) complexes formed with these ligands in their ground states establish oxidation state-dependent conformational dynamicity, while time-resolved emission and transient absorption spectroscopies define the photopysical parameters of photo-induced excited states. Building on initial reports with a related set of molecules, the improved ligand design presented here greatly simplifies the observed photophysics, effectively shutting down unwanted ligand-centered excited states previously observed. Finally, time-dependent density functional theory (TDDFT) analyses reveal an unusual metal-to-TICT electronic transition only reported once before, and though the formation of a CS state is not observed directly through experiments, TDDFT geometry optimizations in the excited states support the formation of transient Cu(II) CS species, lending credence to the potential success of our approach. These studies establish a clear model for the excited state dynamics at play in proof-of-concept systems and clarify key design parameters for future optimizations toward achieving long-lived CS via photoinduced conformational gating.
Photoluminescence spectra of point defects in semiconductors: Validation of first-principles calculations
Optically and magnetically active point defects in semiconductors are interesting platforms for the development of solid state quantum technologies. Their optical properties are usually probed by measuring photoluminescence spectra, which provide information on excitation energies and on the interaction of electrons with lattice vibrations. We present a combined computational and experimental study of photoluminescence spectra of defects in diamond and SiC, aimed at assessing the validity of theoretical and numerical approximations used in first-principles calculations, including the use of the Franck-Condon principle and the displaced harmonic oscillator approximation. We focus on prototypical examples of solid state qubits, the divacancy centers in SiC and the nitrogen-vacancy in diamond, and we report computed photoluminescence spectra as a function of temperature that are in very good agreement with the measured ones. As expected we find that the use of hybrid functionals leads to more accurate results than semilocal functionals. Interestingly our calculations show that constrained density functional theory (CDFT) and time-dependent hybrid DFT perform equally well in describing the excited state potential energy surface of triplet states; our findings indicate that CDFT, a relatively cheap computational approach, is sufficiently accurate for the calculations of photoluminescence spectra of the defects studied here. Finally, we find that only by correcting for finite-size effects and extrapolating to the dilute limit can one obtain a good agreement between theory and experiment. Our results provide a detailed validation protocol of first-principles calculations of photoluminescence spectra, necessary both for the interpretation of experiments and for robust predictions of the electronic properties of point defects in semiconductors.
The Fe 2 O 3 (0001) Surface Under Electroreduction Conditions: A DFT Study of L-Cysteine Adsorption
Local catalyst surface structure and environment can play a significant role towards catalytic activity and selectivity. Surface functionalization using organic additives, such as amino acid chains or peptides, can alter surface properties. Density Functional Theory calculations are used to evaluate the potential dependent surface stability of different terminations of the Fe 2 O 3 (0001) surface. Adsorption of L-Cysteine in different redox states and through different binding modes (carboxylic: O-Fe, amine: N-Fe and thiol: S-Fe) is evaluated. At moderate electrochemical reducing conditions, Fe 2 O 3 (0001) exposes a partially reduced termination with both surface H atoms and undercoordinated Fe atoms in the outermost layer. L-Cysteine adsorption occurs most preferentially through carboxylic acid, O-Fe, binding and does not significantly alter the relative surface stability of different surface terminations. A partially reduced surface with L-Cysteine functionalization will be stable under electroreduction conditions. As a result, stable functionalization of an oxide material through amino acid chains or peptide adsorption may provide an additional design lever to develop improved catalytic systems.
Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U
Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.
2025 ASMS Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS
Title (20 words): Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS Introduction (120 words): One of the challenges facing wide-spread adoption of nuclear power is the development of efficient separation processes for used nuclear fuel. The molecules in separation processes are subjected to an extreme environment due to the high radiation fields from the used fuel and highly acidic media used for fuel dissolution, which results in significant molecular degradation, leading to reduced process efficiency. These degradation products must be identified and studied so mitigation strategies can be developed to maintain process efficiency. However, complex systems can have many degradation products, complicating identification. Untargeted analysis tools could be used to understand radiation chemistry in complex systems. However, this would necessitate improved understanding of the gas-phase fragmentation mechanisms of fuel cycle molecules like tetraoctyldiglycolamide (TODGA). Methods (120 words): The gas-phase fragmentation of protonated TODGA was investigated using collision-induced dissociation (CID), resonance ejection, and infrared ion spectroscopy (IRIS). CID and resonance ejection experiments were conducted using a Bruker Daltonics (Bremen, Gemany) SolariX XR fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. IRIS spectra of protonated TODGA and its two CID fragmentation products were measured using a modified Bruker amaZon Speed ETD 3D quadrupole ion trap mass spectrometer coupled to the Free Electron Lasers for Infrared eXperiments (FELIX) free electron laser. Measured spectra were compared with density functional theory (DFT) calculations using the Gaussian 16, Revision C.02 software package with the ?B97X-D functional and def2-TZVPP basis sets. Candidate structures were generated using the CREST 3.0 conformational sampling software tool. Preliminary Data (300 words): Collision-induced dissociation of protonated TODGA ([C36H73N2O3]+, m/z=581.562) results two fragment ions, one at m/z=340.285 assigned as [C20H38NO3]+ and the other at m/z=312.290, assigned as [C19H38NO2]+. Based on the assigned formula and the structure of protonated TODGA, the fragment at m/z=340.285 is likely formed from elimination of neutral dioctylamine. Comparison of the IRIS spectrum of m/z=340.285 with DFT predictions suggests it contains a ring structure, and is assigned as N-octyl-N-(6-oxo-1,4-dioxan-2-ylidene)octan-1-aminium. Based on this structure and the structure of protonated TODGA, we hypothesize this fragment formed from elimination of neutral dioctylamine followed by a ring closure mechanism. Comparison of the IRIS spectrum of the fragment at m/z=312.290 with DFT predictions also indicated the presence of a ring structure, assigned as N-(1,3-dioxolan-4-ylidene)-N-octyloctan-1-aminium. This product could be formed from elimination of carbon monoxide from the ring of m/z=340.285 as a sequential fragmentation or formed directly from protonated TODGA via elimination of neutral N,N-dioctylformamide followed by a ring closure. Resonance ejection experiments where m/z=340 was continuously ejected from the IRC cell showed no decrease in intensity of m/z=312.290 across several collision energies, suggesting that the later, direct formation mechanism, dominates. The location of the ionizing proton in protonated TODGA is important for modeling the fragmentation mechanisms. DFT calculations suggested that the position of bands involving the coupled vibrations of the amide C—N and C=O bonds in TODGA are the most sensitive to proton location. Evaluation of the IRIS spectrum of protonated TODGA suggests that the ionizing proton is located between the two amid oxygens. This protonation location was calculated to lie approximately 30 kJ/mol lower in energy than the next lowest energy location, with the proton located solely on one of the amide oxygens. Novel aspect (20 words): Infrared ion spectroscopy combined with resonance ejection experiments and density functional theory to probe the collision-induced dissociation mechanism of tetraoctyldiglycolamide.
Linear magnetoresistance with a universal energy scale in the strong-coupling superconductor Mo8Ga41 without quantum criticality
The recent discovery of a nonsaturating linear magnetoresistance in several correlated electron systems near a quantum critical point has revealed an interesting interplay between the linear magnetoresistance and the zero-field linear-in-temperature resistivity. These studies suggest a possible role of quantum criticality on the observed linear magnetoresistance. Here we report our discovery of a nonsaturating, linear magnetoresistance in Mo 8 Ga 41 , a nearly isotropic strong electron-phonon coupling superconductor with a linear-in-temperature resistivity from the transition temperature to ~ 55 K. The growth of the resistivity in field is comparable to that in temperature, provided that both quantities are measured in the energy unit. Overall, our data sets are remarkably similar to magnetoresistance data of the optimally doped La 2 - x Sr x CuO 4 , despite the clearly different crystal and electronic structures, and the apparent absence of quantum critical physics in Mo 8 Ga 41 . A new empirical scaling formula is developed, which is able to capture the key features of the low-temperature magnetoresistance data of Mo 8 Ga 41 , as well as the data of La 2 - x Sr x CuO 4 .