Engineering Papers⌕ Search

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.

At least 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Coordination-Induced N–H Bond Weakening in a Molybdenum Pyrrolidine Complex: Isotopic Labeling Provides Insight into the Pathway for H 2 Evolution

The synthesis and characterization of a cationic molybdenum pyrrolidine complex is described that exhibits significant coordination-induced N–H bond weakening. Here, the N–H bond dissociation free energy (BDFE) of the coordinated pyrrolidine in [( Ph Tpy)(PPh 2 Me) 2 Mo(NH(pyrr))][BArF 24 ] ([1-NH(pyrr)] + ; PhTpy = 4'-Ph-2,2',6',2''-terpyridine, NH(pyrr) = pyrrolidine, ArF 24 = [C 6 H 3 -3,5-(CF 3 ) 2 ] 4 ) was determined to be between 41–51 kcal mol -1 by thermochemical analysis and supported by a density functional theory (DFT) computed value of 48 kcal mol -1 . The complex [1-NH(pyrr)] + underwent proton-coupled electron transfer (PCET) to 2,4,6-tri-tert-butylphenoxyl radical, as well as spontaneous H 2 evolution upon gentle heating to furnish the corresponding molybdenum pyrrolidide complex, [( Ph Tpy)(PPh 2 Me) 2 Mo(N(pyrr))][BArF 24 ] ([1-N(pyrr)] + ). Thermolysis of the deuterated isotopolog, [1-ND(pyrr)] + still produced H 2 with concomitant incorporation of the isotopic label into the pyrrolidide ligand in the product [(1-N(pyrr-d n )] + (n = 0–2), consistent with an H 2 evolution pathway involving intramolecular H–H bond formation followed by an intermolecular product-forming PCET step. These observations provide context for understanding H 2 evolution in the nonclassical ammine complex [( Ph Tpy)(PPh 2 Me) 2 Mo(NH 3 )][BArF 24 ] ([1-NH 3 ] + ) and are supported by DFT-computed reaction thermochemistry. Overall, these studies offer rare insight into the H 2 formation pathway in nonclassical amine complexes with N–H BDFEs below the thermodynamic threshold for H 2 evolution and inform the development of well-defined, thermody-namically potent PCET reagents.

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

Synergistic Ru Co atomic pair with enhanced activity toward levulinic acid hydrogenation

Development of efficient metal-based catalysts is of great importance for levulinic acid (LA) hydrogenation to γ-valerolactone (GVL). The widely employed Ru-based catalysts are advantageous for H 2 dissociation, however, the steric hindrance for large Ru particles hampers their coordination to C=O moiety in LA, and thereby decreasing the activity. Herein, we report a Ru 1 Co 1 -N-C double single-atom catalyst (DSAC) with synergistic Ru and Co atomic pairs for LA hydrogenation into GVL. The Ru and Co doped zeolitic imidazole frameworks (RuCo-doped ZIF-8) precursor was rationally designed ((Ru+Co)/(Zn+Ru+Co) = 2 at.%), where the Zn node spatially isolates Ru and Co species, expanding the adjacent Ru-Co distance and facilitating the formation of the Ru-Co atomic pair upon pyrolysis, with each atom coordinated with three nitrogen atoms (N 3 -Ru 1 Co 1 -N 3 ). The Ru 1 Co 1 -N-C catalyst exhibits outstanding catalytic activity, with a turnover frequency (TOF) of 1980 h –1 , surpassing previously reported Ru-based catalysts. Experimental investigation and density functional theory (DFT) calculations reveal that the electron-rich Ru induced by less electronegative Co facilitates H 2 dissociation, while atomic Ru in dual-atomic pairs promotes C=O activation, Ru and Co atomic pairs synergistically enhancing LA conversion to GVL. In conclusion, this research will shed light on the precise control of active sites at atomic scale, and also provides a new concept for designing high-performance Ru-based catalysts towards LA hydrogenation to GVL.

Double single-atom catalysts↗
DOE OSTI2026

Task Sharing of Proton Incorporation in Vertically Aligned Nanocomposite Triple Conductors: Growth, Structure, and Surface Exchange Kinetics

As protonic ceramic electrolysis cells emerge for efficient H 2 production, there is a need to develop air electrode materials enabling fast, durable steam splitting and proton incorporation. Single-phase triple conductors may fail to satisfy the myriad performance/stability requirements, and their critical charge-carriers (holes, oxygen vacancies, and protons) are in competition, limiting their concentrations. Instead, we propose task-sharing, vertically aligned nanocomposites (VANs), comprising a proton conductor (BaZr 0.9 Y 0.1 O 3-δ ) and a redox-active mixed ionic electronic conductor (Ce 0.9 Pr 0.1 O 2-δ ), that may enable rapid proton surface exchange at the solid–gas interface and transport along the solid–solid heterointerfaces. We grew VANs by pulsed laser deposition and investigated the interplay between their processing conditions, structure, and proton and oxygen surface exchange kinetics. We varied the substrate temperature, laser repetition rate, laser fluence, and processing oxygen pressure. The crystallinity and phases were characterized by grazing-incidence X-ray diffraction, and the strain and structural order as a function of depth were evaluated by angle-dependent synchrotron X-ray pair distribution function analysis. To evaluate the potential for interdiffusion, the formation energies of substitutional defects were simulated with density functional theory. Corresponding structural analysis and elemental mapping were performed by scanning/transmission electron microscopy, energy-dispersive X-ray spectroscopy, and electron energy-loss spectroscopy, indicating distinct nanoscale compositional regions with a hierarchical structure embedded in individual VANs columns and minimal interdiffusion across a bilayer film. Proton and oxygen surface exchange coefficients (k H , k O ) and polarization resistances were evaluated by electrical and optical relaxations and impedance spectroscopy of VAN-incorporated protonic ceramic electrochemical cells, respectively, at 400–500 °C, demonstrating values comparable to some of the best-known triple and mixed conductors.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Competing polar phases in 2D ferroelectric transition metal thio- and selenophosphates

The aim of this Perspective is to provide an in-depth discussion of a certain class of 2D ferroelectrics that feature a van der Waals gap where more than one polar phase can exist energetically close to the ground state. Polar phases of interest can include different ferroelectric, antiferroelectric, and paraelectric phases, which can be transformed into each other through external stimuli, offering unprecedented control over polar material properties. With this level of control, important size limitations and integration problems can be overcome, and material platforms with tunable polar properties for multi-functional devices and sensors can be developed for the next generation of information and energy technologies. We highlight the co-existence of polar properties that have been found on local scales using advanced local characterization techniques as well as methods to control polar properties in transition metal thio-/selenophosphates. In conclusion, we give an overview of future challenges and opportunities in this emerging field of energetically close polar phases.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Density functional study of atomic arrangements in CrMnFeCoNi high-entropy alloy and their impact on vacancy formation energy and segregation

Here, using the density functional theory-coupled Monte Carlo approach, we explored the chemical short-range order (SRO) and element segregation in equimolar CrMnFeCoNi alloy. We found that state-of-the-art approximation of random element distribution is only applicable at > 1100 K close to the melting temperature, while the Cr-Cr repulsion driving the system stabilization and accompanying the formation of cubic Cr sublattice, and mild Ni-Ni attraction are the most prominent pair interactions at lower temperatures. Chemical potential and vacancy formation energy calculations indicate that Cr is most sensitive to the local chemical environment, making Cr atoms most stabilized when the preferred SRO is introduced. While the vacancy formation is predicted equally probable among five constituting elements in the random solid solution, Cr and Ni atoms show the lowest vacancy formation energies in the structure with SRO. Furthermore, distinct element segregation was predicted in the vicinity of planar defects, including symmetric tilt grain boundary and stacking fault, which we correlated to the site- and chemistry-dependent atomic volume and bond lengths. It suggests that the local mechanical strain and bond energy induce the SRO development and element segregation: Namely, the system takes advantage of segregation of Ni atoms having large atomic volume or Cr-Cr pairs having elongated bond lengths to fill in the excess volume at defects that relaxes the mechanical strain field and optimizes bond energy distribution. The correlation between the SRO and properties of CrMnFeCoNi alloy needs further investigations, which is expected to greatly help understand and control the properties of high-entropy alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A Transferable Force Field for Simulating Adsorption in Metal–Organic Frameworks with Open Metal Sites Based on the 12–6–4 Lennard-Jones Potential

Metal−organic frameworks (MOFs) that contain coordinatively unsaturated open metal sites (OMSs) provide strong host− guest interactions, making them promising sorbents for low-concentration gas adsorption applications such as direct air capture and atmospheric water harvesting. However, accurately modeling host−guest interactions involving OMSs remains challenging for classical force fields (FFs) based on the 12−6 Lennard−Jones (LJ) potential, as the polarization effect of the guest molecule induced by the positively charged OMS is not considered. Here, we introduce an FF based on the 12−6−4 LJ potential, which incorporates charge−induced dipole interactions and is parametrized against a diverse set of host−guest potential energy surfaces (PESs) obtained from density functional theory (DFT). The resulting FF, trained on a generic trimetallic cluster, performs well in both host−guest binding energetics and gas adsorption isotherms across different OMS-containing MOFs, including MOF-74 series and Cu-BTC. These results highlight the excellent transferability of our approach and its potential to enhance the accuracy and robustness of high-throughput MOF discovery workflows, particularly for gas adsorption and separation in large and diverse MOF databases.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Measurement of the sound velocity and Grüneisen parameter of polystyrene at inertial confinement fusion conditions

The principal Hugoniot, sound velocity, and Grüneisen parameter of polystyrene were measured at conditions relevant to shocks in inertial confinement fusion implosions, from 100 to 1000 GPa. The sound velocity is in good agreement with quantum molecular dynamics calculations and all tabular equation of state models at pressures below 200 GPa. Above 200 GPa, the experimental results agree with two of the examined tables, but do not agree with the most recent table developed for design of ICF experiments. The Grüneisen parameter increases with density below ~3.1 g/cm 3 and approaches the asymptotic value for an ideal gas after complete dissociation. This behavior is in good agreement with quantum molecular dynamics results and previous work but is not represented by any of the tabular models. Finally, the discrepancy between tabular models and experimental measurement of the sound velocity and Grüneisen parameter is sufficient to impact simulations of ICF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ab initio quantum dynamics of plasmonic charge carriers

Surface plasmons respond strongly to electric fields and generate energetic (hot) charge carriers that can be used in optoelectronic devices. However, utilization of hot carriers has to outcompete fast carrier relaxation in metals. Often, the reported efficiencies of nanoscale devices based on plasmon excitations are low and the mechanisms of device operation remain inconsistent within the field. Further developments hinge on fundamental understanding of the nature and kinetics of plasmon decay processes. In this work, we provide an overview of time domain ab initio modeling of hot carrier dynamics, demonstrating that atomistic details of the surface structure of plasmonic materials, and their chemical interaction with semiconducting and molecular charge acceptors and substrates, play key roles in governing the dynamics pathways. Considered case studies represent various dynamics channels and illustrate solar energy and optoelectronic applications.

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

Functional Renormalization Group Study of Superconductivity in Rhombohedral Trilayer Graphene

We employ a functional renormalization group approach to ascertain the pairing mechanism and symmetry of the superconducting phase observed in rhombohedral trilayer graphene. Superconductivity in this system occurs in a regime of carrier density and displacement field with a weakly distorted annular Fermi sea. We find that repulsive Coulomb interactions can induce electron pairing on the Fermi surface by taking advantage of momentum-space structure associated with the finite width of the Fermi sea annulus. Furthermore, the degeneracy between spin-singlet and spin-triplet pairing is lifted by valley-exchange interactions that strengthen under the RG flow and develop nontrivial momentum-space structure. We find that the leading pairing instability is d-wave-like and spin-singlet, and that the theoretical phase diagram versus carrier density and displacement field agrees qualitatively with experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Direct aromatization of CO 2 via combined CO 2 hydrogenation and zeolite-based acid catalysis

Aromatics, including benzene, toluene, and xylenes (BTX), are essential chemical building blocks and are widely used as solvents, fuel additives, and polymers. With the recent development in CO 2 capture technologies and the progress made in producing H 2 using renewable energy, direct hydrogenation of CO 2 to aromatics via heterogeneous catalysis has emerged as a promising pathway to accomplish the production of aromatics with simultaneous utilization of waste CO 2 . In this review, we focus on recent advances in the nascent field of direct CO 2 aromatization, whereby tandem catalysts composed of CO 2 hydrogenation and aromatization functionalities are designed and deployed. We review two categories of tandem catalysts: catalysts integrating Fe-based/H-ZSM-5 components following RWGS (reverse water-gas shift of CO 2 to CO)-FT (Fischer-Tropsch synthesis of lower olefins)-aromatization pathways, and catalysts combining metal oxide/H-ZSM-5 domains following CO 2 to methanolaromatization pathways. The key parameters that determine the catalytic performance, such as the composition and structure of the Fe-based or metal oxide-based CO 2 conversion catalysts, the properties of H-ZSM-5, and the synergy between the two components, are analyzed to provide insights for the design of efficient tandem catalysts for CO 2 aromatization. In parallel, thermodynamic analyses, mechanistic studies, and density functional theory (DFT) computations for the relevant reaction routes and pathways are discussed to offer improved understanding of CO 2 activation, reaction intermediates, and product formation. In closing, the challenges and prospects for these tandem reactions are addressed to provide suggested paths forward for future research.

10 SYNTHETIC FUELS↗
DOE OSTI2024

Characterization of Deformational Isomerization Potential and Interconversion Dynamics with Ultrafast X-ray Solution Scattering

Dimeric complexes composed of d 8 square planar metal centers and rigid bridging ligands provide model systems to understand the interplay between attractive dispersion forces and steric strain in order to assist the development of reliable methods to model metal dimer complexes more broadly. [Ir 2 (dimen) 4 ] 2+ (dimen = para-diisocyanomenthane) presents a unique case study for such phenomena, as distortions of the optimal structure of a ligand with limited conformational flexibility counteract the attractive dispersive forces from the metal and ligand to yield a complex with two ground state deformational isomers. Here, in this work, we use ultrafast X-ray solution scattering (XSS) and optical transient absorption spectroscopy (OTAS) to reveal the nature of the equilibrium distribution and the exchange rate between the deformational isomers. The two ground state isomers have spectrally distinct electronic excitations that enable the selective excitation of one isomer or the other using a femtosecond duration pulse of visible light. We then track the dynamics of the nonequilibrium depletion of the electronic ground state population—often termed the ground state hole—with ultrafast XSS and OTAS, revealing a restoration of the ground state equilibrium in 2.3 ps. This combined experimental and theoretical study provides a critical test of various density functional approximations in the description of bridged d 8 –d 8 metal complexes. The results show that density functional theory calculations can reproduce the primary experimental observations if dispersion interactions are added, and a hybrid functional, which includes exact exchange, is used.

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

Chemical Design of Spin Frustration to Realize Topological Spin Glasses

Patterning spins to generate collective behavior is at the core of condensed matter physics. Physicists develop techniques, including the fabrication of magnetic nanostructures and precision layering of materials specifically to engender frustrated lattices. As chemists, we can access such exotic materials through targeted chemical synthesis and create new lattice types by chemical design. Here, in this study, we introduce a new approach to induce magnetic frustration on a modified honeycomb lattice through a competition of alternating antiferromagnetic (AFM) and ferromagnetic (FM) nearest-neighbor interactions. By subtly modulating these two types of interactions through facile synthetic modifications, we created two systems: (1) a topological spin glass and (2) a frustrated spin-canted magnet with low-temperature exchange bias. To design this unconventional magnetic lattice, we used a metal–organic framework (MOF) platform, Ni 3 (pymca) 3 X 3 (NipymcaX where pymca = pyrimidine-2-carboxylato and X = Cl, Br). We isolated two MOFs, NipymcaCl and NipymcaBr, featuring canted Ni 2+ -based moments. Despite this similarity, differences in the single-ion anisotropies of the Ni 2+ spins result in distinct magnetic properties for each material. NipymcaCl is a topological spin glass, while NipymcaBr is a rare frustrated magnet with low-temperature exchange bias. Density functional theory calculations and Monte Carlo simulations on the NipymcaX lattice support the presence of magnetic frustration as a result of alternating AFM and FM interactions. Our calculations enabled us to determine the ground-state spin configuration and the distribution of spin–spin correlations relative to paradigmatic kagomé and triangular lattices. This modified honeycomb lattice is similar to the electronic Kekulé-O phase in graphene and provides a highly tunable platform to realize unconventional spin physics.

Exchange bias↗
DOE OSTI2024

Segregation of Al and its effect on coercivity in Nd-Fe-B

Doping a small amount of Al can effectively enhance coercivity in Nd-Fe-B magnets. Here we investigated the partitioning behavior of Al and its effect on coercivity in Nd-Fe-B using first principles DFT (density functional theory) calculation and micromagnetic simulation. The calculated substitution energies of Fe by Al are negative at the crystallographic sites of 4c and 8j 2 while they are positive values at the other sites in Nd 2 Fe 14 B (2:14:1), implying a small solubility of Al in 2:14:1. Further, Al prefers to segregate at grain boundary (GB) and stabilize the Nd-Fe-Al phase with a Nd 6 Fe 13 Si-type tetragonal structure (6:13:1). The formation of the antiferromagnetic or weak ferrimagnetic 6:13:1-like phase depletes Fe and reduces the amount of ferromagnetic Nd–Fe type grain boundary phase (GBP), which weakens the inter-grain magnetic interaction. Micromagnetic simulations indicate that the 6:13:1-like GBP increases the pinning field of magnetic domain wall at GB and suppresses the nucleation of reversal magnetic domain on the grain surface of 2:14:1 during demagnetization process. The formation of Al-rich shell on 2:14:1 grain surface can further moderately increase the domain pinning field at GB and the nucleation field of 2:14:1 grain. Developing novel processing method to tailor Al segregation and promote formation of 6:13:1-like phase at GB can be a promising approach to improve coercivity in Nd-Fe-B magnet.

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

Interplay between electron localization, magnetic order, and Jahn-Teller distortion dictates LiMnO2 phase stability

The development of manganese (Mn)-rich cathodes for Li-ion batteries promises to alleviate potential supply chain bottlenecks in battery manufacturing. Fundamental challenges in Mn-rich cathodes arise from phenomena such as structural changes due to cooperative Jahn-Teller (JT) distortions of in octahedral environments, Mn migration, and phase transformations to spinel-like order, all of which affect the electrochemical performance. These physically complex phenomena motivate an re-examination of the Li-Mn-O rock-salt space, with a focus on the thermodynamics of the prototypical, polymorphs. It is found that the generalized gradient approximation (GGA-PBEsol) and meta-GGA ( ) density functionals with empirically fitted on-site Hubbard corrections yield spurious stable phases for , such as predicting a phase with -like order ( ) to be the ground state instead of the orthorhombic (Pmmn) phase, which is the experimentally known ground state. Accounting for antiferromagnetic order in each structure is shown to have a substantial effect on the total energies and resulting phase stability. By using hybrid-GGA (HSE06) and GGA with self-consistent Hubbard parameters (on-site and inter-site ) calculated from linear response theory, the experimentally observed phase stability trends are recovered. The calculated on-site between Mn- states in the experimentally observed orthorhombic, layered, and spinel phases are significantly smaller than in and disordered layered structures, by within GGA. The smaller values of are shown to be correlated with a collinear ordering of JT distortions, in which all orbitals are oriented in the same direction. This cooperative JT effect can lead to greater electron delocalization from Mn along the states due to increased Mn-O covalency, which contributes to the greater electronic stability compared to the phases with noncollinear JT arrangements. The structures with collinear ordering of JT distortions also generate greater vibrational entropy, which helps stabilize these phases at high temperature. These phases are shown to be strongly insulating with large calculated band gaps , which are computed using HSE06 and .

Kam, Ronald L↗
DOE OSTI2022

Machine learning for design principles for single atom catalysts towards electrochemical reactions

Machine learning (ML) integrated density functional theory (DFT) calculations have recently been used to accelerate the design and discovery of heterogeneous catalysts such as single atom catalysts (SACs) through the establishment of deep structure–activity relationships. Here, this review provides recent progress in the ML-aided rational design of heterogeneous catalysts with the focus on SACs in terms of structure–activity relationships, feature importance analysis, high-throughput screening, stability, and metal–support interactions for electrochemistry. Support vector machine (SVM), random forest regression (RFR), and deep neural networks (DNN) along with atomic properties are mainly used for the design of SACs. The ML results have shown that the number of electrons in the d orbital, oxide formation enthalpy, ionization energy, Bader charge, d-band center, and enthalpy of vaporization are mainly the most important parameters for the defining of the structure–activity relationships for electrochemistry. However, the black-box nature of ML techniques occasionally makes a physical interpretation of descriptors, such as the Bader charge, d-band center, and enthalpy of vaporization, non-trivial. At the current stage, ML application is limited by the lack of a large and high-quality database. Future prospects for the development of a large database and a generalized ML algorithm for SAC design are discussed to give insights for further studies in this field.

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

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Chemistry↗
DOE OSTI2023

Probing lithium mobility at a solid electrolyte surface

Abstract Solid-state electrolytes overcome many challenges of present-day lithium ion batteries, such as safety hazards and dendrite formation 1,2 . However, detailed understanding of the involved lithium dynamics is missing due to a lack of in operando measurements with chemical and interfacial specificity. Here we investigate a prototypical solid-state electrolyte using linear and nonlinear extreme-ultraviolet spectroscopies. Leveraging the surface sensitivity of extreme-ultraviolet-second-harmonic-generation spectroscopy, we obtained a direct spectral signature of surface lithium ions, showing a distinct blueshift relative to bulk absorption spectra. First-principles simulations attributed the shift to transitions from the lithium 1 s state to hybridized Li- s /Ti- d orbitals at the surface. Our calculations further suggest a reduction in lithium interfacial mobility due to suppressed low-frequency rattling modes, which is the fundamental origin of the large interfacial resistance in this material. Our findings pave the way for new optimization strategies to develop these electrochemical devices via interfacial engineering of lithium ions.

25 ENERGY STORAGE↗
DOE OSTI2021

COMPASO : A new halo finder for competitive assignment to spherical overdensities

We describe a new method (compaso) for identifying groups of particles in cosmological N-body simulations. compaso builds upon existing spherical overdensity (SO) algorithms by taking into consideration the tidal radius around a smaller halo before competitively assigning halo membership to the particles. In this way, the compaso finder allows for more effective deblending of haloes in close proximity as well as the formation of new haloes on the outskirts of larger ones. This halo-finding algorithm is used in the abacussummit suite of N-body simulations, designed to meet the cosmological simulation requirements of the Dark Energy Spectroscopic Instrument (DESI) survey. compaso is developed as a highly efficient on-the-fly group finder, which is crucial for enabling good load-balancing between the GPU and CPU and the creation of high-resolution merger trees. In this paper, we describe the halo-finding procedure and its particular implementation in abacus, accompanying it with a qualitative analysis of the finder. We test the robustness of the compaso catalogues before and after applying the cleaning method described in an accompanying paper and demonstrate its effectiveness by comparing it with other validation techniques. We then visualize the haloes and their density profiles, finding that they are well fit by the NFW formalism. Finally, we compare other properties such as radius–mass relationships and two-point correlation functions with that of another widely used halo finder, rockstar.

79 ASTRONOMY AND ASTROPHYSICS↗