Search NASA⌕ Search

SEARCH · Search NASA

Results for “cluster expansion”

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

Unifying principles for catalytic hydrotreating processes (Final Technical Report)

This project builds on the hypothesis that the hydrotreating processes for the removal of oxygen and sulfur are fundamentally similar at the atomic-scale and existing knowledge from the treatment of petroleum derived feedstock can be leveraged for the design of novel catalysts for the upgrade of bio-oil. We tested this hypothesis by comparing computed potential energy diagrams for hydrodesulfurization (HDS) of thiophene over MoS 2 with hydrodeoxygenation (HDO) of furan over MoO 3 and concluded that certain aspects, such as catalyst promotion with transition metals, are valid strategies for both reactions. On the other hand, we also noticed significant differences in the mechanism for hydrogen (H 2 ) activation, which requires sites with metallic character. While MoS 2 is known to have metallic edge states that can catalyze H 2 dissociation, this elementary step is prohibitively slow on defect-free oxides. Only in the presence of vacancies or by creating metal/oxide interfaces can efficient H 2 activation sites during HDO be formed. The need for bifunctional catalyst when it comes to efficient and selective HDO or dehydrogenation reactions was further corroborated in joint experimental and theoretical studies of the Guerbet reaction for the coupling of biomass derived oxygenates over PdCu alloys, nitrate reduction over In-promoted Pd nanoparticles, and ethylene dehydroaromatization over Ga-exchanged ZSM-5 zeolites. All of these catalytic systems have in common that catalytic sites with distinct functional requirements are needed to create a working catalyst. Detailed computational studies were carried out for HDO of m -cresol and phenol on Ru-modified TiO 2 surfaces, which allowed us to attribute catalytic activity to the metal/oxide interface. A surprising finding was that heterolytic cleavage of the H-H bond across the Ru/TiO 2 interface was critically important, despite lower barriers for homolytic H 2 activation on Ru metal. The explanation lies in the high barriers for hydrogen spillover from Ru to TiO 2 , which becomes unnecessary in the heterolytic activation pathway. Moreover, we also reported that proton transfer steps between metal and oxide sites are mediated by weakly adsorbed surface water. During attempts to develop and validate a kinetic Monte Carlo (kMC) model for HDO reactions at the Ru/TiO 2 interface, it became clear that lateral interactions are paramount to describe realistic surface chemistry and without these interactions, the reduction and hydroxylation behavior from our simulations was inconsistent with reported experiments. To assess the importance of lateral interactions in popular computational catalyst design strategies relying on the identification of reactivity descriptors, which can be used along with Brønsted–Evans–Polanyi (BEP) and scaling relations as input to a microkinetic model (MKM) to make predictions for activity or selectivity trends, we compared predicted trends with those obtained from descriptor-based kMC models. We critically evaluated the benefits of kMC over MKM in terms of trend predictions and computational cost when using only a small set of input parameters. After confirming that in the absence of lateral interactions the kMC and MKM approaches yield identical trends and mechanistic information, we observed substantial differences between the two kinetic models when lateral interactions were introduced. The mean-field implementation applies coverage corrections directly to the descriptors, causing an artificial overprediction of the activity of strongly binding metals. In contrast, the cluster expansion in kMC implementation can differentiate among the highly active metals but it is very sensitive to the set of included interaction parameters. Considering that computational screening relies on a minimal set of descriptors, for which MKM makes reasonable trend predictions at a ca. three orders of magnitude lower computational cost than kMC, we concluded that the MKM approach does provide an overall better entry point for computational catalyst design. Overall, this project has led to 11 peer-reviewed publications, and their scientific is impact is well illustrated by their combined 702 citations.

08 HYDROGEN↗

Modeling Neutron Star Matter in the Age of Multimessenger Astrophysics

Abstract The interpretation of available and forthcoming data obtained from multimessenger astrophysical observations—potentially providing unprecedented access to neutron star properties—will require the development of novel, accurate theoretical models of dense matter. Of great importance, in this context, will be the capability for devising a description of thermal effects applicable to the study of quantities other than the equation of state, such as the transport coefficients and the neutrino mean free path in the nuclear medium. A formalism based on correlated basis states and the cluster expansion technique has been previously employed to derive a well-behaved effective interaction—suitable for use in standard perturbation theory—from a state-of-the-art nuclear Hamiltonian, including phenomenological two- and three-nucleon potentials. Here, we provide a comprehensive and self-contained account of the extension of this approach to the treatment of finite-temperature effects, and report the results of numerical calculations of a number of properties of nuclear matter with arbitrary neutron excess and temperature up to 50 MeV.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistical Mechanics on Lattices (smol) v0.0.1

smol is a Python implementation of the Cluster Expansion Method and extensions of this methodology; which is used to fit applied lattice models from first principle calculations for subsequent use in Monte Carlo simulations for computing thermodynamic and statistical properties of atomic configuration. smol has a several notable advantages over other similar software packages. The first being its Python implementation which makes it easy to use even for users with little software and coding skills. Despite being implemented in Python, essential Monte Carlo routines are implemented in Cython such that performance is not compromised and is competitive even to similar software implemented in C/C++. Furthermore, the package has specific functionality for handling complex ionic materials (such as cathodes and electrolytes) that is not readily available in other packages. Lastly, the software has a flexible and modular design with the intention of making it fast and efficient to develop and extend the methodology.

Barroso-Luque, Luis↗

Vacancy-mediated transport and segregation tendencies of solutes in FCC Ni under diffusional creep

Nickel-based alloys are known for their excellent corrosion resistance and high-temperature strength. However, at high temperatures, polycrystals experience grain boundary-stress-driven vacancy diffusion leading to diffusion-driven creep. Understanding the fundamentals of defect migration is essential to predict diffusion creep behavior. In this talk, we explore the Nabarro-Herring diffusion creep and the associated solute segregation that is likely to result from the different diffusivities of the solute and solvent elements. Density functional theory calculations are performed to determine the formation, migration and binding energies of the solute and vacancy defects in stress-free and stressed states. The dilute-limit transport coefficients are then calculated using the self-consistent mean field theory as implemented in the Kinetic Cluster Expansion (KineCluE) code. Predictions of diffusion creep strain rates and solute segregation are then presented for various solutes, stress and temperatures. This work opens research avenues for exploration of diffusion creep via characterization of elemental segregation to grain boundaries.

36 MATERIALS SCIENCE↗

Understanding grain boundary segregation under irradiation: A mesoscale modeling study

Grain boundaries (GBs) are an important component in the design of irradiation resistant materials as they play a crucial role as point defect sinks for non-equilibrium vacancies and self-interstitials. However, in the presence of sustained point defect fluxes, the coupled transport of atomic components with inherently different diffusivities often leads to radiation-induced segregation (RIS), with undesired effects on properties such as intergranular corrosion and stress-assisted cracking. In this talk, we investigate the fundamental mechanisms leading to GB segregation in BCC FeCr-based alloys using multiscale modeling and experiments. A phase-field model is employed to account for both thermodynamic and kinetic segregation mechanisms. The full Onsager transport matrix is derived using kinetic Monte Carlo, molecular dynamics, and kinetic cluster expansion (KineCluE) calculations. The thermodynamics of GB segregation is described via an atomic density-based modification to the CALPHAD free energy. The modeling predictions are then compared against segregation data from atom probe characterization for various temperatures, compositions, and neutron irradiation conditions. We will conclude by discussing the utility of the multiscale modeling approach to shed light on the fundamental mechanisms leading to solute-GB interaction.

36 MATERIALS SCIENCE↗

Spectral function for 4 He using the Chebyshev expansion in coupled-cluster theory

Here, we compute spectral function for 4 He by combining coupled-cluster theory with an expansion of integral transforms into Chebyshev polynomials. Our method allows us to estimate the uncertainty of spectral reconstruction. The properties of the Chebyshev polynomials make the procedure numerically stable and considerably lower in memory usage than the typically employed Lanczos algorithm. We benchmark our predictions with other calculations in the literature and with electron-scattering data in the quasi-elastic peak. The spectral function formalism allows one to extend ab initio lepton-nucleus cross sections into the relativistic regime. This makes it a promising tool for modeling this process at higher-energy transfers. The results we present open the door for studies of heavier nuclei, important for the neutrino oscillation programs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Spin Decoherence Dynamics of Er 3+ in CeO 2 Films

Developing telecom-compatible spin-photon interfaces is essential towards scalable quantum networks. Erbium ions (Er 3+ ) exhibit a unique combination of a telecom (1.5 mu m) optical transition and an effective spin-1=2 ground state, but identifying a host that enables heterogeneous device integration while preserving long optical and spin coherence remains an open challenge. In this work, we explore the potential of Er 3+ :CeO 2 films on silicon and study the Er 3+ spin coherence, offering low nuclear spin density and the potential for on-chip integration. We demonstrate a 38.8 mu s spin coherence, which can be extended to 176.4 mu s with dynamical decoupling. Pairing experiments with cluster correlation expansion calculations, we identify spectral diffusion induced by bath Er 3+ spin flip-flops as the dominant decoherence mechanism and provide pathways to millisecond-scale coherence.

36 MATERIALS SCIENCE↗

High-level coupled-cluster energetics by merging moment expansions with selected configuration interaction

Inspired by our earlier semi-stochastic work aimed at converging high-level coupled-cluster (CC) energetics, we propose a novel form of the CC(P; Q) theory in which the stochastic Quantum Monte Carlo propagations, used to identify dominant higher-than-doubly excited determinants, are replaced by the selected configuration interaction (CI) approach using the perturbative selection made iteratively (CIPSI) algorithm. The advantages of the resulting CIPSI-driven CC(P; Q) methodology are illustrated by a few molecular examples, including the dissociation of F 2 and the automerization of cyclobutadiene, where we recover the electronic energies corresponding to the CC calculations with a full treatment of singles, doubles, and triples based on the information extracted from compact CI wave functions originating from relatively inexpensive Hamiltonian diagonalizations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Whole-Genome Comparisons of Ergot Fungi Reveals the Divergence and Evolution of Species within the Genus Claviceps Are the Result of Varying Mechanisms Driving Genome Evolution and Host Range Expansion

The genus Claviceps has been known for centuries as an economically important fungal genus for pharmacology and agricultural research. Only recently have researchers begun to unravel the evolutionary history of the genus, with origins in South America and classification of four distinct sections through ecological, morphological, and metabolic features (Claviceps sects. Citrinae, Paspalorum, Pusillae, and Claviceps). The first three sections are additionally characterized by narrow host range, whereas section Claviceps is considered evolutionarily more successful and adaptable as it has the largest host range and biogeographical distribution. However, the reasons for this success and adaptability remain unclear. Our study elucidates factors influencing adaptability by sequencing and annotating 50 Claviceps genomes, representing 21 species, for a comprehensive comparison of genome architecture and plasticity in relation to host range potential. Our results show the trajectory from specialized genomes (sects. Citrinae and Paspalorum) toward adaptive genomes (sects. Pusillae and Claviceps) through colocalization of transposable elements around predicted effectors and a putative loss of repeat-induced point mutation resulting in unconstrained tandem gene duplication coinciding with increased host range potential and speciation. Alterations of genomic architecture and plasticity can substantially influence and shape the evolutionary trajectory of fungal pathogens and their adaptability. Furthermore, our study provides a large increase in available genomic resources to propel future studies of Claviceps in pharmacology and agricultural research, as well as, research into deeper understanding of the evolution of adaptable plant pathogens.

59 BASIC BIOLOGICAL SCIENCES↗

Effect of Fluorination on Lithium Transport and Short-Range Order in Disordered-Rocksalt-Type Lithium-Ion Battery Cathodes

Fluorine substitution is a critical enabler for improving the cycle life and energy density of disordered rocksalt (DRX) Li-ion battery cathode materials which offer prospects for high energy density cathodes, without the reliance on limited mineral resources. Due to the strong Li–F interaction, fluorine also is expected to modify the short-range cation order in these materials which is critical for Li-ion transport. In this work, density functional theory and Monte Carlo simulations are combined to investigate the impact of Li–F short-range ordering on the formation of Li percolation and diffusion in DRX materials. Additionally, the modeling reveals that F substitution is always beneficial at sufficiently high concentrations and can, surprisingly, even facilitate percolation in compounds without Li excess, giving them the ability to incorporate more transition metal redox capacity and thereby higher energy density. It is found that for F levels below 15%, its effect can be beneficial or disadvantageous depending on the intrinsic short-range order in the unfluorinated oxide, while for high fluorination levels the effects are always beneficial. Using extensive simulations, a map is also presented showing the trade-off between transition-metal capacity, Li-transport, and synthetic accessibility, and two of the more extreme predictions are experimentally confirmed.

25 ENERGY STORAGE↗

Design, Processing, and Properties of WTaCrV-Hf Multi-principal Element Alloys

Refractory multi-principal element alloys are candidates for high-temperature structural components due, in part, to their high strength and high melting points. Single-phase materials are initially preferred for isotropic material properties as a function of time and temperature in service conditions. This work outlines a computational rank-ordering and experimental validation methodology for single-phase body-centered-cubic phase stability in WTaCrV-Hf alloys using order–disorder transition temperature. Eight compositions were fabricated by arc-melting and heat-treated at 1400 °C for 24 hrs. X-ray diffraction, energy-dispersive x-ray spectroscopy, and Vickers hardness testing showed alloys with order–disorder transition temperatures below 600 °C formed a single-phase body-centered-cubic structure during solidification and remained single-phase after heat-treatment. The sample possessing the lowest order–disorder transition temperature exhibited slip traces suggestive of room-temperature plastic deformation under Vickers indentation, with both heat-treated single-phase samples exhibiting hardnesses over 800 HV with little cracking compared to tungsten. These results establish order–disorder transition temperature as a viable predictive parameter for multi-principal element alloy phase stability. The methodology outlined in this work provides a framework for future design, fabrication, and characterization of high-temperature structural multi-principal element alloys.

CALPHAD↗

Microstructure modeling of nuclear structural materials: Recent progress and future directions

Modeling and simulation of microstructures are essential to understand the complex responses and behaviors of nuclear materials in extreme environments. The needs to assess the extended life operation as well as the growing interest in accelerating nuclear materials development and qualification have stimulated the use of high-fidelity multiscale models aided by empirical and ab initio data. This paper reviews the role of various models across different length and time scales in investigating irradiation effects on microstructure evolution and degradation, in particular the embrittlement caused by radiation induced or enhanced formation of nanoscale chemical heterogeneities. The strength and limitations of these models, including classical rate theories, cluster dynamics, phase-field methods, and atomistic models informed by ab initio energies, are discussed with seminal examples. Challenges regarding the lack of thermo-kinetic data and theoretical treatments considering chemical complexities and magnetic excitations, as well as the stabilizing effect by excess point defects in nuclear structural materials are presented, along with potential solutions based on ab initio informed surrogate energy models and statistical sampling by Monte Carlo simulations. Further, the review then highlights the opportunities to leverage the advantages of different methods by establishing hybrid models by shared variables or coupled codes and applications. Finally, the review concludes with forward-looking remarks on how the use of physics-based models can aid the improvement of machine-learning models of property degradation and vice versa.

36 MATERIALS SCIENCE↗

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗

Spin Decoherence in VOPc@graphene Nanoribbon Complexes

Carbon nanoribbons or nanographene qubit arrays can facilitate quantum-to-quantum transduction between light, charge, and spin, making them an excellent testbed for fundamental science in quantum coherent systems and for the construction of higher-level qubit circuits. Here, in this work, we study spin decoherence due to coupling with a surrounding nuclear spin bath of an electronic molecular spin of a vanadyl phthalocyanine (VOPc) molecule integrated onto an armchair-edged graphene nanoribbon (GNR). Density functional theory (DFT) is used to obtain ground-state atomic configurations. Decay of spin coherence in Hahn echo experiments is then simulated using the cluster correlation expansion method with a spin Hamiltonian involving hyperfine and electric field gradient tensors calculated from DFT. We find that the decoherence time T 2 is anisotropic with respect to magnetic field orientation and determined only by the hydrogen nuclear spins on both VOPc and GNR. Large electron spin echo envelope modulation (ESEEM) due to nitrogen and vanadium nuclear spins is present at specific field ranges and can be completely suppressed by tuning the magnetic field. The relation between these field ranges and the hyperfine interactions is analyzed. The effects of interactions with the nuclear quadrupole moments are also studied, validating the applicability and limitations of the spin Hamiltonian when they are disregarded.

Hamiltonians↗

Capturing Surface Coverage Effects in Heterogeneous Catalysis

Adsorbate–adsorbate lateral interactions at relevant surface coverages have a significant effect on chemical kinetics, thereby influencing the activity of a heterogeneous catalyst. Coverage-dependent kinetic and thermodynamic parameters therefore must be included in studies of such complex systems to properly predict the turnover frequencies and kinetic trends. Thus, it becomes extremely important to accurately capture the strength of lateral interactions between neighboring species under realistic reaction conditions. In this Perspective, we discuss the various existing computational and experimental methods for determining adspecies coverage and configurational effects. The choice of the tools and methods employed in such studies depends on factors such as time, length scales, computational cost, the presence of solvents, and reaction conditions. The applications of each method and the respective challenges are also discussed here. As a result, we discuss the recent developments and future of the state-of-the-art for inclusion of surface coverage and configuration into a holistic picture for accurate predictions of catalytic behavior.

09 BIOMASS FUELS↗

Decoherence of $V^{–}_{B}$ spin defects in monoisotopic hexagonal boron nitride

Spin defects in hexagonal boron nitride (hBN) are promising quantum systems for the design of flexible two-dimensional quantum sensing platforms. Here we rely on hBN crystals isotopically enriched with either 10 B or 11 B to investigate the isotope-dependent properties of a spin defect featuring a broadband photoluminescence signal in the near infrared. By analyzing the hyperfine structure of the spin defect while changing the boron isotope, we first confirm that it corresponds to the negatively charged boron-vacancy center ($V^{–}_{B}$). We then show that its spin coherence properties are slightly improved in 10B-enriched samples. This is supported by numerical simulations employing cluster correlation expansion methods, which reveal the importance of the hyperfine Fermi contact term for calculating the coherence time of point defects in hBN. Using cross-relaxation spectroscopy, we finally identify dark electron spin impurities as an additional source of decoherence. This work provides new insights into the properties of $V^{–}_{B}$ spin defects, which are valuable for the future development of hBN-based quantum sensing foils.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin–spin interactions in defects in solids from mixed all-electron and pseudopotential first-principles calculations

Abstract Understanding the quantum dynamics of spin defects and their coherence properties requires an accurate modeling of spin-spin interaction in solids and molecules, for example by using spin Hamiltonians with parameters obtained from first principles calculations. We present a real-space approach based on density functional theory for the calculation of spin-Hamiltonian parameters, where only selected atoms are treated at the all-electron level, while the rest of the system is described with the pseudopotential approximation. Our approach permits calculations for systems containing more than 1000 atoms, as demonstrated for defects in diamond and silicon carbide. We show that only a small number of atoms surrounding the defect needs to be treated at the all-electron level, in order to obtain an overall all-electron accuracy for hyperfine and zero-field splitting tensors. We also present results for coherence times, computed with the cluster correlation expansion method, highlighting the importance of accurate spin-Hamiltonian parameters for quantitative predictions of spin dynamics.

36 MATERIALS SCIENCE↗

Generalized scaling of spin qubit coherence in over 12,000 host materials

Significance Atomic defects in solid-state materials are promising candidates as quantum bits, or qubits. New materials are actively being investigated as hosts for new defect qubits; however, there are no unifying guidelines that can quantitatively predict qubit performance in a new material. One of the most critical property of qubits is their quantum coherence. While cluster correlation expansion (CCE) techniques are useful to simulate the coherence of electron spins in defects, they are computationally expensive to investigate broad classes of stable materials. Using CCE simulations, we reveal a general scaling relation between the electron spin coherence time and the properties of qubit host materials that enables rapid and quantitative exploration of new materials hosting spin defects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗