Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase errors”

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 613 records · Page 34

Lefschetz thimble quantum Monte Carlo for spin systems

Monte Carlo simulations are useful tools for modeling quantum systems, but in some cases they suffer from a sign problem, leading to an exponential slow down in their convergence to a value. While solving the sign problem is generically NP hard, many techniques exist for mitigating the sign problem in specific cases; in particular, the technique of deforming the Monte Carlo simulation's plane of integration onto Lefschetz thimbles (complex hypersurfaces of stationary phase) has seen significant success in the context of quantum field theories. We extend this methodology to spin systems by utilizing spin coherent state path integrals to reexpress the spin system's partition function in terms of continuous variables. Using some toy systems, we demonstrate its effectiveness at lessening the sign problem in this setting, despite the fact that the initial mapping to spin coherent states introduces its own sign problem. The standard formulation of the spin coherent path integral is known to make use of uncontrolled approximations; despite this, for large spins they are typically considered to yield accurate results, so it is somewhat surprising that our results show significant systematic errors. Furthermore, possibly of independent interest, our use of Lefschetz thimbles to overcome the intrinsic sign problem in spin coherent state path integral Monte Carlo enables a novel numerical demonstration of a breakdown in the spin coherent path integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Wave Prediction using X-Band Radar (Final Summary Report)

Model Predictive Control (MPC) is the best controls framework available to wave energy today. It allows for the maximization of energy yield from a WEC system while respecting various constraints and losses in the system. This type of constrained optimization is a critical capability for the techno-economic optimization of WEC devices and is difficult to do effectively with causal controls frameworks. MPC requires a future prediction of wave excitation forces, requiring wave prediction on a future time horizon of up to 30 seconds. The future horizon required is highly dependent on the WEC and PTO type available. To our knowledge no-one has been able to implement MPC on a WEC device at sea, due to the fact that phase-resolved wave-prediction is not a capability that has been sufficiently developed to date. The key objectives of this project are focused on developing “industry-ready” wave prediction technology building blocks that can be applied to a wide range of different wave energy conversion topologies. Our collaborative work with device developers has demonstrated that the level of improvement attainable for any particular device is heavily dependent on the device and PTO topology chosen. The range of annual average power capture improvements is on the order of 10% to over 100%. However, the important aspect is that we can get within about 10% of theoretical upper limits even when considering errors in the wave prediction and more importantly, we can optimize performance while respecting various constraints such as motion amplitude or peak structural loads. The combination of improved performance while limiting structural loads has a net effect of reducing the levelized cost of energy from these emerging technologies and allows us to design control laws that meet an economic optimum. Efforts under a previous project on controls has focused on developing a set of wave-prediction algorithms that work by leveraging a network of wave sensing measurement buoys that are equipped with real-time telemetry. The present work is focusing primarily on using X-band radar to measure the wave-field around the WEC to predict the future excitation forces on the structure and combining the measurements of these different sources to create improvements in the wave-prediction accuracy.

16 TIDAL AND WAVE POWER↗

APOGEE-2S Mg–Al anti-correlation of the metal-poor globular cluster NGC 2298

We present detailed elemental abundances and radial velocities of stars in the metal-poor globular cluster (GC) NGC 2298, based on near-infrared high-resolution ($R$ ~ 22 500) spectra of 12 members obtained during the second phase of the Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) at Las Campanas Observatory as part of the seventeenth Data Release (DR 17) of the Sloan Digital Sky Survey IV (SDSS-IV). We employed the Brussels Automatic Code for Characterizing High accuracy Spectra (BACCHUS) software to investigate abundances for a variety of species including α elements (Mg, Si, and Ca), the odd-Z element Al, and iron-peak elements (Fe and Ni) located in the innermost regions of NGC 2298. We find a mean and median metallicity [Fe/H] = -1.76 and -1.75, respectively, with a star-to-star spread of 0.14 dex, which is compatible with the internal measurement errors. Thus, we find no evidence for an intrinsic [Fe/H] abundance spread in NGC 2298. The typical $α$-element enrichment in NGC 2298 is overabundant relative to the Sun, and it follows the trend of other metal-poor GCs. We confirm the existence of an Al-enhanced population in this cluster, which is clearly anti-correlated with Mg, indicating the prevalence of the multiple-population phenomenon in NGC 2298.

79 ASTRONOMY AND ASTROPHYSICS↗

High-speed Focal Plane Wave Front Sensing with an Optical Chopper

Focal plane wave front sensing and control is a critical approach to reducing noncommon path errors between a conventional astronomical adaptive optics (AO) wave front sensor (WFS) detector and a science camera. However, in addition to mitigating noncommon path errors, recent focal plane wave front sensing techniques have been developed to operate at speeds fast enough to enable "multi-WFS" AO, where residual atmospheric errors are further corrected by a focal plane WFS. Although a number of such techniques have been recently developed for coronagraphic imaging, here we present one designed for noncoronagraphic imaging. Utilizing conventional AO system components, this concept additionally requires (1) a detector imaging the focal plane of the WFS light source and (2) a pupil plane optical chopper device that is the noncommon path to the first WFS and is synchronized to the focal plane imager readout. These minimal hardware requirements enable the temporal amplitude modulation to resolve the sine ambiguity of even wave-front modes for low, mid, and high wave front spatial frequencies. Similar capabilities have been demonstrated with classical phase diversity by defocusing the detector, but such techniques are incompatible with simultaneous science observations. This optical chopping technique, however, enables science imaging at up to 50% duty cycle. In conclusion, we present both simulations and laboratory validation of this concept on SEAL, the Santa Cruz Extreme AO Laboratory testbed.

79 ASTRONOMY AND ASTROPHYSICS↗

Discovering Blind Geothermal Systems in the Great Basin Region: An Integrated Geologic and Geophysical Approach for Establishing Geothermal Play Fairways: All Phases

Most geothermal resources in the Great Basin region of the western USA are blind, and thus the discovery of new commercial-grade systems requires synthesis of favorable characteristics for geothermal activity. The geothermal play fairway concept involves integration of multiple parameters indicative of geothermal activity to identify promising areas for new development. This project integrated multiple datasets to apply the play fairway concept and assess geothermal potential in a large region of the Great Basin in Nevada. It is therefore referred to as the Nevada play fairway project. This project was a strong collaborative effort between several organizations, led by the Nevada Bureau of Mines and Geology at the University of Nevada, Reno, but with key support from the U.S. Geological Survey, ATLAS Geosciences, Inc,, Hi-Q Geophysical, Inc., Lawrence Berkeley National Laboratory, Utah Geological Survey, and Innovative Geothermal Ltd. In Budget Period 1 of this project, available data for nine geologic, geochemical, and geophysical parameters were initially synthesized to produce a new detailed geothermal potential map of 96,000 km2 from west-central to eastern Nevada. These parameters were grouped into subsets and individually weighted to delineate rankings for local permeability, intermediate permeability, regional permeability, and thermal potential, which collectively defined geothermal play fairways (i.e., most likely locations for significant geothermal fluid flow). This initial work was aimed at reducing the risks in regional exploration and therefore facilitating discovery of new commercial-grade systems in blind settings, as well as in areas with surface expressions of geothermal activity. Budget Period 2 of the project involved detailed analysis of some of the most promising areas identified in Phase 1. Twenty-four highly prospective areas, including both known undeveloped systems and previously undiscovered potential blind systems, were identified for further analysis. After reconnaissance of these areas, five of the most promising sites were selected for detailed studies. Multiple techniques were employed in the detailed studies, including geologic mapping, shallow temperature surveys, gravity surveys, Lidar, geochemical studies, seismic reflection analysis, and 3D modeling. The goal of the detailed studies was to identify specific areas with the highest likelihood for high permeability and thermal fluids, such that drill sites could be targeted. Three main sets of predictive maps were generated for each detailed study area: 1) play fairway maps, 2) play fairway error maps, and 3) direct evidence maps. Local- and intermediate-scale permeability models were revised to reflect results of the detailed geologic, geophysical, and geochemical analyses. Budget Period 3 of the project involved more detailed geophysical analyses and temperature-gradient (TG) drilling in southeastern Gabbs Valley and northern Granite Springs Valley, deemed the two most promising sites, with the goal of providing preliminary validation of the play fairway methodology. In southeastern Gabbs Valley, the collocation of a favorable structural setting (displacement transfer zone and fault intersections), Quaternary faults, intersecting and terminating gravity gradients, magnetic low, shallow (2 m) temperature anomaly, low resistivity anomaly, and promising geothermometry from nearby water wells provided evidence for a blind system. Drilling of six TG holes defines an apparent geothermal system at this locality with temperatures as high as 124°C at 152 m. This system is blind, with no surface hot springs, fumaroles, or paleo-geothermal deposits. For northern Granite Springs Valley, a favorable structural setting (termination of a major Quaternary normal fault), terminating gravity gradient, magnetic gradient, newly discovered sinter deposits, nearby warm water wells, previously drilled TG holes in the vicinity, and promising geothermometry suggest a hidden system. Drilling of six new TG holes yields temperatures of ~96°C at ~250 m, suggesting the presence of a geothermal system. Major lessons learned in the course of this project include: 1) initially identified sites commonly include multiple favorable structural settings at a finer scale; 2) promising sites in Cenozoic basins cannot be recognized without detailed geophysical surveys; and 3) play fairway analysis should be refined as the exploration program vectors into the most promising sites and finer-scale data are acquired. In addition to producing copious amounts of data, this project resulted in 16 published papers, 10 abstracts, more than 40 presentations across the U.S. and abroad (including several keynote addresses), 2 Masters theses, and 7 media reports.

15 GEOTHERMAL ENERGY↗

Modeling Subsurface Explosions Recorded on a Distributed Fiber Optic Sensor

Fiber optic distributed acoustic sensors (DAS) are becoming a widely used tool for seismic sensing. Here we examine recordings of two subsurface chemical explosions, DAG-1 and DAG-3, each of which was about one metric ton (TNT equivalent), that were recorded from a helical fiber installed in two boreholes 80 m away from the source location. Several clear phases including the initial P wave, a weak S wave, and a surface reflected P wave are observed on the helical DAS data. We estimate a velocity model using arrival times measured from the fiber. The DAS waveform data were compared with colocated accelerometers at specific depths in both frequency and time domains. Furthermore, the spectra of the DAS data matched spectra estimated from the accelerometer records. Comparisons of observed waveform shape between the accelerometer records and the fiber measurements (strain-rate) show reasonable agreement except for the data near the event depth. The DAS data and the accelerometer agreed in relative amplitudes but we had difficulties in matching absolute amplitudes, possibly due to errors in metadata. Synthetic strain-rate waveforms were calculated using a 2D wavenumber algorithm and matched the waveform shape and relative amplitudes. In general, DAS is effective at recording strong ground motions at high spatial density. Comparison of the synthetic seismograms with observed data indicate that the waveforms are not consistent with a pure isotropic explosion source and that the observed S waves originate from very near the source region.

58 GEOSCIENCES↗

Oak Ridge National Laboratory Current Magnitude and Frequency Response of SIPROTEC Relays

This project examined the current magnitude and frequency response of 7SJ61 and 7SJ62 SIPROTEC relays. The measured sample frequencies for the 7SJ61 and 7SJS62 relays were 800 and 960 Hz, respectively. The relays were tested at frequencies of up to 1 and 20.16 kHz. Two experiments were run in a lab environment. The first experiment was run at frequencies of up to 1 kHz every 10 Hz. The second one was run at frequencies of up to 20.16 kHz every 200 and 240 Hz (common factors of sample rate frequency) for the 7SJ61 and 7SJ62 SIPROTEC relays, respectively. Both experiments were set at a current magnitude of 0.1 and 0.5 Arms. A testbed with a Relay Test Equipment (RTE) was set to run the experiment at frequencies of up to 1 kHz, while a testbed with a RTE and High Frequency Power Source/ Amplifier (HFPS/A) was set to run at frequencies of up to 20.16 kHz. Currents of 0.1 and 0.5 Arms at different frequencies were connected at Phase A of the relays, and a current of 1 Arms (60 Hz) was connected at Phase B of relays. The 1 Arms current allowed the relays to trip and saved their record events by setting both relays at a definite-time overcurrent protection function. The record events were collected after the 7SJ61 and 7SJ62 relays tripped, and the current magnitude and frequency were plotted. The record events from the relays were collected using the DIGSI® 4 software. The SYNCHROWAVE® Event and MATLAB® software were then used to plot the current magnitude and frequency, respectively, from the relay record events. The current magnitude and frequency were plotted at time and frequency domain, respectively.The 7SJ61 and 7SJ62 relays measured source frequencies of up to 400 and 480 Hz, respectively. The percentage error of frequency for the relays increased after 400 and 480 Hz, respectively. The 7SJ62 relay with 0.1 and 0.5 Arms measured the alias frequency up to 1440 Hz. The 7SJ61 relay with 0.1 and 0.5 Arms measured the alias frequency up to 1000 and 1200 Hz, respectively. At frequencies higher than the alias frequency, the relays measured a DC amplitude current instead of an AC amplitude current. In the DC amplitude current zone, non-desired frequencies were measured by both relays. The current magnitude intensity (dB) for the 7SJ61 and 7SJ62 relays from 1 to 20.16 kHz were plotted. The current magnitude intensity for the 7SJ61 relay decreased up to 2600 and 1200 Hz for 0.5 and 0.1 Arms tests, respectively. However, the current magnitude intensity for the 7SJ62 relay decreased up to 4800 and 2640 Hz for 0.5 and 0.1 Arms tests, respectively. The current magnitude and frequency response study for the 7SJ61 and 7SJ62 relays was performed satisfactorily, observing the source, alias, and non-desired measured frequency zones for both relays.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A first principles framework to predict the transient performance of latent heat thermal energy storage

Thermal energy storage (TES) is increasingly recognized as an essential component of efficient Combined Heat and Power (CHP), Concentrated Solar Power (CSP), Heating Ventilation and Air Conditioning (HVAC), and refrigeration as it reduces peak demand while helping to manage intermittent availability of energy (e.g., from solar or wind). Latent Heat Thermal Energy Storage (LHTES) is a viable option because of its high energy storage density. Parametric analysis of LHTES in terms of dimensionless numbers is highly desired as a tool to model LHTES systems. One approach is to develop a model equation so as to minimize the error between the model and data obtained from experiments or simulations. While this approach can produce an accurate correlation applicable within the range of data used for its creation, it does not provide physical understanding of the rate-limiting process controlling the transient behavior of the device. In this paper we present an alternative approach whereby the potential rate-limiting processes are identified from first principles and then the key process is determined as a function of time as a LHTES device is charged. For example, in a simple geometry, the melt-fraction can be expected to vary linearly in time if the heat transfer rate is limited by natural convection of the phase changing material and we show it scales with the PCM Grashof number as $Gr^1_p$ and PCM Prandtl number as $Pr_p^{(1/3)}$. On the other hand, if surface area of solid PCM limits the heat transfer rate, the melt fraction increases asymptotically to reach full melting. The existence of these linear and asymptotic regions and the $Gr^1_pP r^{1/3}_p$ shape of the melt fraction curve is verified using our database of 64 simulations. Of practical importance in designing LHTES devices is the melt fraction at which the heat transfer rate ceases to be limited by convection, after which the heat storage rate deteriorates. For our geometry, this is found to be about 90%. This test case of our methodology shows the value of our approach, that predicting heat storage rate based on the rate-limiting physical phenomenon as a function of time is an effective approach to modeling LHTES devices.

25 ENERGY STORAGE↗

Worldwide benchmarking of cost-effective radiometers for direct and diffuse irradiance

Solar energy projects can benefit from direct normal irradiance (DNI) and diffuse horizontal irradiance (DHI) measurements during all project phases. Several commercial measurement systems for DNI and DHI are available. Sun trackers with pyranometers and pyrheliometers can provide highly accurate measurements but are often impractical in solar energy applications. For less expensive and more robust sensors, it is often unclear which accuracy can be expected under a project site's specific atmospheric conditions. We address this challenge through our dedicated experimental comparison of relevant sensor systems (rotating shadowband irradiometer [short RSI], Delta-T SPN1, EKO MS-90, PyranoCam, Sunto CaptPro, Kipp & Zonen CSD3) at up to six sites worldwide. The RSI systems (rRMSD 3 to 8.6%, DNI; 4.8 to 7.6%, DHI) and PyranoCam (rRMSD 2.6 to 5.2%, DNI; 4.4 to 5.8%, DHI) exhibit similar error metrics and are the most accurate systems in the test. Delta-T SPN1 and EKO MS-90 (rRMSD 6.8 to 15%, DNI; 10.6 to 20.1%, DHI) but especially Kipp & Zonen CSD3 and Sunto CaptPro show significant deviations (rRMSD 17.7 to 20%, DNI; 33 to 58%, DHI). We evaluate the influence of relevant atmospheric parameters on the sensors' accuracies by a rather unique measurement setup. MS-90's DNI errors depend on DNI itself, with overestimations for low reference DNI. The deviations of SPN1's DHI and DNI measurements increase sharply in situations with high circumsolar irradiance. Also CaptPro and CSD3's increased measurement errors are related to circumsolar irradiance. For RSI and PyranoCam, only moderate influences on the measurements are identified, indicating a general applicability of these instruments.

14 SOLAR ENERGY↗

Critical role of magnetic moments in heavy-fermion materials: Revisiting mysteries of SmB 6

Heavy-fermion family exhibits fascinating and often puzzling properties due to the presence of open-shell f ions and the complexity of the associated charge, orbital, and spin degrees of freedom. SmB 6 is a prototypical heavy-fermion compound that is electrically insulating but yet it displays quantum oscillations, which are a telltale signature of the metallic state. Adding to the enigma is the possibility that SmB 6 is a topological Kondo insulator. Here, by treating the spin degree of freedom on an equal footing with other degrees of freedom using the parameter-free strongly-constrained and appropriately-normed (SCAN) density functional, we explore the ground-state electronic structure of SmB6. A number of competing magnetic phases lying very closely in energy are found, indicating the key role of spin fluctuations in the material. The computed band structure, crystal-field splittings in the f-electron complex, the heavy effective electron mass at the Fermi energy, and the large specific heat are all in good agreement with the corresponding experimental results. In particular, our predicted FS explains the experimentally observed bulk quantum oscillations as well as the low electrical conductivity of SmB 6 . The topological Kondo state of SmB6 is shown to be robust regardless of its magnetic configuration. The excellent performance of SCAN in heavy-fermion systems is explained in terms of its ability to treat self-interaction errors and symmetry breaking within the framework of the density functional theory. Our study provides a new approach for modeling heavy-fermion materials

36 MATERIALS SCIENCE↗

Correlation between Spin and Orbital Dynamics during Laser-Induced Femtosecond Demagnetization

Spin and orbital angular momenta are two intrinsic properties of an electron and are responsible for the physics of a solid. How the spin and orbital evolve with respect to each other on several hundred femtoseconds is largely unknown, but it is at the center of laser-induced ultrafast demagnetization. In this paper, we introduce a concept of the spin–orbital correlation diagram, where spin angular momentum is plotted against orbital angular momentum, much like the position-velocity phase diagram in classical mechanics. We use four sets of highly accurate time-resolved X-ray magnetic circular dichroism data to construct four correlation diagrams for iron and cobalt. To our surprise, a pattern emerges. The trace on the correlation diagram for iron is an arc, and at the end of demagnetization, it has a pronounced cusp. The correlation diagram for cobalt is different and appears more linear but with kinks. We carry out first-principles calculations with two different methods: time-dependent density functional theory (TDDFT) and time-dependent Liouville density functional theory. These two methods agree that the experimental findings for both Fe and Co are not due to experimental errors. It is the spin–orbit coupling that correlates the spin dynamics to the orbital dynamics. Microscopically, Fe and Co have different orbital occupations, which leads to distinctive correlation diagrams. We believe that this correlation diagram presents a useful tool to better understand spin and orbital dynamics on an ultrafast time scale. A brief discussion on the magnetic anisotropy energy is also provided.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Graphical Gaussian Process Regression Model for Aqueous Solvation Free Energy Prediction of Organic Molecules in Redox Flow Battery

The solvation free energy of organic molecules is a critical parameter in determining emergent properties such as solubility, liquid-phase equilibrium constants, and pKa and redox potentials in an organic redox flow battery. In this work, we present a machine learning (ML) model that can learn and predict the aqueous solvation free energy of an organic molecule using Gaussian process regression method based on a new molecular graph kernel. To investigate the performance of the ML model on electrostatic interaction, the nonpolar interaction contribution of solvent and the conformational entropy of solute in solvation free energy, three data sets with implicit or explicit water solvent models, and contribution of conformational entropy of solute are tested. We demonstrate that our ML model can predict the solvation free energy of molecules at chemical accuracy with a mean absolute error of less than 1 kcal/mol for subsets of the QM9 dataset and the Freesolv database. To solve the general data scarcity problem for a graph-based ML model, we propose a dimension reduction algorithm based on the distance between molecular graphs, which can be used to examine the diversity of the molecular data set. It provides a promising way to build a minimum training set to improve prediction for certain test sets where the space of molecular structures is predetermined.

25 ENERGY STORAGE↗

Finite and infinite matrix product states for Gutzwiller projected mean-field wave functions

Matrix product states (MPS) and “dressed” ground states of quadratic mean fields (e.g., Gutzwiller projected Slater determinants) are both important classes of variational wave functions. This latter class has played important roles in understanding superconductivity and quantum spin liquids. We present a method to obtain both the finite and infinite MPS (iMPS) representation of the ground state of an arbitrary fermionic quadratic mean-field Hamiltonian (which in the simplest case is a Slater determinant and in the most general case is a Pfaffian). We also show how to represent products of such states (e.g., determinants times Pfaffians). From this representation one can project to single occupancy and evaluate the entanglement spectra after Gutzwiller projection. We then obtain the MPS and iMPS representation of Gutzwiller projected mean-field states that arise from the variational slave-fermion approach to the S = 1 bilinear-biquadratic quantum spin chain. To accomplish this, we develop an approach to orthogonalize degenerate iMPS to find all the states in the degenerate ground-state manifold. We find the energies of the MPS and iMPS states match the variational energies closely, indicating the method is accurate and there is minimal loss due to truncation error. We then present an exploration of the entanglement spectra of projected slave-fermion states, exploring their qualitative features and finding good qualitative agreement with the respective exact ground-state spectra found from density matrix renormalization group.

1-dimensional spin chains↗

Multiresolution GPC-Structured Control of a Single-Loop Cold-Flow Chemical Looping Testbed

Chemical looping is a near-zero emission process for generating power from coal. It is based on a multi-phase gas-solid flow and has extremely challenging nonlinear, multi-scale dynamics with jumps, producing large dynamic model uncertainty, which renders traditional robust control techniques, such as linear parameter varying H ∞ design, largely inapplicable. This process complexity is addressed in the present work through the temporal and the spatiotemporal multiresolution modeling along with the corresponding model-based control laws. Namely, the nonlinear autoregressive with exogenous input model structure, nonlinear in the wavelet basis, but linear in parameters, is used to identify the dominant temporal chemical looping process dynamics. The control inputs and the wavelet model parameters are calculated by optimizing a quadratic cost function using a gradient descent method. The respective identification and tracking error convergence of the proposed self-tuning identification and control schemes, the latter using the unconstrained generalized predictive control structure, is separately ascertained through the Lyapunov stability theorem. The rate constraint on the control signal in the temporal control law is then imposed and the control topology is augmented by an additional control loop with self-tuning deadbeat controller which uses the spatiotemporal wavelet riser dynamics representation. The novelty of this work is three-fold: (1) developing the self-tuning controller design methodology that consists in embedding the real-time tunable temporal highly nonlinear, but linearly parametrizable, multiresolution system representations into the classical rate-constrained generalized predictive quadratic optimal control structure, (2) augmenting the temporal multiresolution loop by a more complex spatiotemporal multiresolution self-tuning deadbeat control loop, and (3) demonstrating the effectiveness of the proposed methodology in producing fast recursive real-time algorithms for controlling highly uncertain nonlinear multiscale processes. The latter is shown through the data from the implemented temporal and augmented spatiotemporal solutions of a difficult chemical looping cold flow tracking control problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling Analysis of Ball-Milling Process for Battery-Electrode Synthesis

The mechanical alloying process is a promising method for synthesizing electrode materials for batteries owing to its benefits such as the ability to produce nanostructured, high-performing electrode alloys, no adverse effects on the solid electrolyte for solid-state batteries, stable production of thick electrodes, simple processing steps, and low processing costs. It is gaining intensive attention in the battery industry as one of the best methods to replace the conventional wet-slurry-solvent method, and its application is rapidly increasing these days. However, the operation is currently conducted purely based on trial-and-error methods without fully utilizing the features of its functions. Here, this may be attributed to a lack of understanding of the effect of operating parameters on the alloying process and final products. Surprisingly, there is a scarcity of the literature conducting fundamental research to comprehend the underlying physics of the entire mechanical alloying process, resulting in a significant knowledge gap. To address this knowledge gap, extensive research was conducted. The existing literature on mechanical alloying was reviewed to comprehend the current state of understanding and to discuss the direction for future research. Mathematical expressions were developed to create physics-based models capable of capturing the entire mechanical alloying process, including milling kinetics and defect-enhanced phase evolution. These methods were then applied to investigate the impact of operating parameters such as milling frequency, initial mole ratio of the alloyed materials, density of grinding balls, and energy required for the powders to become amorphous (i.e., the amorphization energy threshold). This research aimed not only to comprehend the direct effects of these operating parameters but also to unveil the physics underlying the ball-milling process. The results of our study can serve as crucial information for the battery industry in designing or operating the ball-milling process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unifying thermochemistry concepts in computational heterogeneous catalysis

Thermophysical properties of adsorbates and gas-phase species define the free energy landscape of heterogeneously catalyzed processes and are pivotal for an atomistic understanding of the catalyst performance. These thermophysical properties, such as the free energy or the enthalpy, are typically derived from density functional theory (DFT) calculations. Enthalpies are species-interdependent properties that are only meaningful when referenced to other species. The widespread use of DFT has led to a proliferation of new energetic data in the literature and databases. However, there is a lack of consistency in how DFT data is referenced and how the associated enthalpies or free energies are stored and reported, leading to challenges in reproducing or utilizing the results of prior work. Additionally, DFT suffers from exchange–correlation errors that often require corrections to align the data with other global thermochemical networks, which are not always clearly documented or explained. In this review, we introduce a set of consistent terminology and definitions, review existing approaches, and unify the techniques using the framework of linear algebra. This set of terminology and tools facilitates the correction and alignment of energies between different data formats and sources, promoting the sharing and reuse of ab initio data. Standardization of thermochemistry concepts in computational heterogeneous catalysis reduces computational cost and enhances fundamental understanding of catalytic processes, which will accelerate the computational design of optimally performing catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimating List-Mode Data Sensitivities to Nuclear Data with MCNP6

Nuclear data are a vital component of predictive simulations used in applications like experiment design, stockpile stewardship, nuclear nonproliferation/safeguards, health physics, and criticality safety. A singular simulation requires the coalescence of different areas of nuclear data such as cross sections, angular distributions, and energy distributions of emitted neutrons for different materials and energy ranges. Improving nuclear data and thus reducing the uncertainty in simulated parameters could enable smaller, better-informed safety factors and ultimately reduce operational and procedural costs. There is a constant effort to garner a better understanding of the physical quantities represented by nuclear data through experiments. Integral experiment benchmarks use simulated and measured results to validate current nuclear data values. In the past, benchmarks primarily focused on the effective multiplication factor (k eff ); however, this limited scope has caused compensating errors and areas of nuclear data that lack validation. Compensating errors are inaccuracies in nuclear data that are obfuscated by cancellation when observing integrated values such as k eff . Diverse integral benchmark experiments that look for quantities of interest other than k eff and include multiple responses minimize the possibility of compensating errors and provides validation to areas of nuclear data previously lacking experimental validation. Benchmark experiments can be optimized during the design process to be highly dependent on specific areas of nuclear data. The dependence of a response in an experiment to a specific area/type of nuclear data is defined as sensitivity. A larger sensitivity means that nuclear data uncertainties will play a larger role in the response(s) resulting in larger bias. Currently, the sensitivity capabilities of the Monte Carlo N-Particle (MCNP ®1 ) transport code are limited to responses of k eff and tallied values (e.g., flux, surface current). As a part of the EUCLID project, this work explores estimating list-mode nuclear data sensitivities that can be used to design experiments aimed to constrain and reduce compensating errors in nuclear data by focusing on responses other than k eff . Tallied values are ideal quantities that are estimated with detectors during experiments. List-mode data (a list of neutron collection times) are the direct output of detector systems in subcritical neutron noise experiments. Expanding MCNP sensitivity capabilities to include the sensitivity of responses estimated from list-mode data, such as the prompt neutron decay constant (α) and multiplicity estimates (S and D), enables more direct comparison of simulated and measured experimental quantities. Additionally, deterministic tools such as SENSMG are capable of obtaining sensitivities to a wide variety of responses; however, these tools cannot handle complex geometries due to the assumptions made in discretizing the phase-space variables of the Boltzman transport equation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗