Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weak form”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Molecular dynamics simulations of uranyl and plutonyl cations in a task-specific ionic liquid

Ionic liquids (ILs) are a unique class of solvents with potential applications in advanced separation technologies relevant to the nuclear industry. ILs are salts with low melting points and a wide range of tunable physical properties, such as viscosity, hydrophobiciy, conductivity, and liquidus range. ILs have negligible vapor pressure, are often non-flammable, and can have high thermal stability and a wide electrochemical window, making them attractive for use in separations processes relevant to the nuclear industry. Metal salts generally have a low solubility in ILs; however, by incorporating new functional groups into the IL cation or anion that promote complexation with the metal, the solubility can be greatly increased. One such task-specific ionic liquid (TSIL) is 1-carboxy-N, N, N-trimethylglycine bis(trifluoromethylsulfonyl)imide ([Hbet][Tf 2 N]). Water, which is detrimental for electrochemical separations, is a common impurity in ILs and can coordinate with actinyl cations, particularly in ILs containing only weakly coordinating components. Understanding the behavior of actinides in TSIL/water mixtures on a molecular level is vital for designing improved separations processes. Classical molecular dynamics simulations of uranyl(VI) and plutonyl(VI) in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]) with deprotonated Hbet (betaine) and water have been performed to understand the coordination and dynamics of the actinyl cations. We find that betaine is a much stronger ligand than water and prefers to coordinate the metal in a bidentate manner. Potential of mean force simulations yield a relative free energy for betaine coordination of approximately -120 to -90 kJ/mol in mixtures with water. As the amount of betaine coordinated to the actinide increases, the diffusion coefficient of the actinyl cation decreases. Moreover, the betaine ligand is able to bridge between two metal centers, resulting in dimeric complexes with actinide–actinide distances of ~5 Å. Potential of mean force simulations show that these structures are stable, with relative free energies of up to -40 kJ/mol. The crystal structure for [(UO 2 ) 2 (bet) 6 (H 2 O) 2 ][Tf 2 N] 4 shows that the betaine bridges between two uranium atoms to form dimeric complexes similar to those found in our simulations.

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

Barium stars as tracers of s -process nucleosynthesis in AGB stars

Barium (Ba) stars help to verify asymptotic giant branch (AGB) star nucleosynthesis models since they experienced pollution from an AGB binary companion and thus their spectra carry the signatures of the slow neutron capture process (s process). For a large number (180) of Ba stars, we searched for AGB stellar models that match the observed abundance patterns. We aim to uncover any systematic deviations of the sample abundances from the predictions of the nucleosynthesis models. We employed three machine learning algorithms as classifiers: a Random Forest method, developed for this work, and the two classifiers used in our previous study. Compared to that work, we also expanded our observational sample with 11 Ba stars available in the supersolar metallicity range. We studied the statistical behaviour of the different s-process elements in the observational sample to investigate if the AGB models systematically under- or overpredict the abundances observed in the Ba stars and show the results in the form of violin plots of the residuals between spectroscopic abundances and model predictions. We inspected the correlations between the observed [Fe/H], the s-process elemental abundances, and the residuals. We employed the [Zr/Fe] and [Nb/Fe] abundances as a thermometer to constrain the operational temperature that rules the production of these elements in the sample stars, assuming a steady-state s process. We also investigated the mass distribution of the identified polluter AGB stars and the behaviour of the δ parameter, which describes the fraction of accreted AGB material relative to the Ba star envelope. We find a significant trend in the residuals that implies an underproduction of the elements just after the first s-process peak (Nb, Mo, and Ru) in the models relative to the observations. This may originate from a neutron-capture process (e.g. the intermediate neutron-capture process, i process) not yet included in the AGB models of metallicity from solar to roughly 1/5 solar, corresponding to the range of the Ba stars. Correlations are found between the residuals of these peculiar elements, suggesting a common origin for the deviations from the models. In addition, there is a weak metallicity dependence of the residuals of these elements. The s-process temperatures derived with the [Zr/Fe] – [Nb/Fe] thermometer have an unrealistic value for the majority of our stars. The most likely explanation is that at least a fraction of these elements are not produced in a steady-state s process, and instead may be due to processes not included in the AGB models. The mass distribution of the identified models confirms that our sample of Ba stars was polluted by low-mass AGB stars (< 4 M ⊙ ). Most of the matching AGB models require low accreted mass, but a few systems with high accreted mass are needed to explain the observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

SEARCHING FOR MESONIC DARK MATTER WITH THE HEAVY PHOTON SEARCH EXPERIMENT

Several highly-sensitive astrophysical experiments over the past couple of decades have demonstrated that the current abundance of visible Standard Model matter cannot explain galactic rotation curves, the expansion history of the Universe, or the apparent warping of light in empty space. Instead, one finds strong agreement with this body of experimental results upon positing the existence of an invisible particulate field, dark matter. Namely, a cold, weakly interacting dark matter component can explain all these phenomena. A number of accelerator-based experiments have been developed to search for the weak couplings/interactions of these particles, many of them concentrating on particle models with masses of tens to thousands of GeV. A relatively new, well-motivated model is a dark sector coupled to the Standard Model via a dark photon. The current abundance of dark matter can be obtained if one assumes that dark matter is coupled to light by a MeV to GeV particle with a U(1) symmetry. The parameter space of these models remains largely unexplored because they are difficult to probe experimentally. In this thesis, I analyze data from the Heavy Photon Search (HPS) detector, whose two detector halves closely surround the electron beam, providing acceptance to far-forward boosted interactions. This forward acceptance to highly boosted particles yields unprecedented sensitivity to MeV-scale invariant masses. I exhaustively optimize the offline reconstruction of the HPS detector. Each reconstruction object, from Silicon Vertex Tracker hits to tracks, is studied to maximize acceptance of dark matter events. I then use the 2021 run data to search for one model of dark-photon-mediated matter, the Strongly Interacting Massive Particle (SIMP). SIMP models provide self-interacting dark matter candidates that can form bound states resembling dark mesons. HPS can detect SIMPs through the decay of a dark vector boson (either a dark ¿ or ¿) into e+e- pairs. I obtain exclusion contours for SIMPs using both an optimized cuts-based selection and a machine-learning-based selection, advancing our knowledge of the nature of dark matter.

O'Dwyer, Rory [Stanford Univ., CA (United States).↗
DOE OSTI2025

Computational Investigation of the Chemical Bond between An(III) Ions and Soft-Donor Ligands

The chemical bonding of actinide ions with arene and borohydride ligands is explored via quantum chemical methods to understand how the transuranium elements interact with softdonor ligands. Specifically, the [An(C 6 Me 6 )(BH 4 ) 3 ] complexes (An = U, Np, and Pu) and their reduced congeners are studied. Density functional theory (DFT) shows that the metal–ligand interactions in the neutral complexes are governed by electrostatic interactions. Both DFT and complete active space (CASSCF) results show that as one moves from U to Pu, the 5f-orbitals are stabilized leading to a poorer energy match with the ligand orbitals. This contributes to progressively weaker metal-arene and metal-borohydride interactions across the series due to a decrease in energy-driven covalency. A reduction in orbital contributions to bonding is obtained for the transuranium-arene interactions as well. Upon reduction, the arene is reduced, forming a δ-bond. This causes the An–arene distances to contract by 0.1–0.2 Å compared to the neutral complexes. The ground state is assigned as the intermediate-spin state where the arene radical is antiferromagnetically coupled to the metal-centered f-electrons in Np and Pu. On the other hand, the ferromagnetically and antiferromagnetically coupled states are close in energy in the uranium complex, but do not mix when spin– orbit coupling is included using a state-interaction approach (SO-CASPT2). The population of the CASSCF δ*-antibonding natural orbital increases from U to Pu consistent with the increased An–arene distances, weaker interactions, and decreasing covalency across the series. Although the An–B distance increases by ca. 0.06 Å upon reduction, both the neutral and reduced species involve an An(III)–borohydride bond and as such are qualitatively similar. The Np complexes can be assigned to have slightly weaker bonding than the uranium analogs but are overall “uranium-like”. The Pu complexes are predicted to have less covalent contributions to bonding in both the Pu–arene and Pu–borohydride interactions; however, the Pu–arene interaction is predicted to be particularly weak.

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

Baryonic effects in the Effective Field Theory of Large-Scale Structure and an analytic recipe for lensing in CMB-S4

Upcoming Large-Scale Structure surveys will likely become the next leading sources of cosmological information, making it crucial to have a precise understanding of the influence of baryons on cosmological observables. The Effective Field Theory of Large-Scale Structure (EFTofLSS) provides a consistent way to predict the clustering of dark matter and baryons on large scales, where their leading corrections in perturbation theory are given by a simple and calculable functional form even after the onset of baryonic processes. In this paper, we extend the two-fluid-like system up to two-loop order in perturbation theory. Along the way, we show that a new linear counterterm proportional to the relative velocity of the fluids could generically be present, but we show that its effects are expected to be small in our universe. Regardless, we show how to consistently perform perturbation theory in the presence of this new term. We find that the EFTofLSS at two-loop order can accurately account for the details of baryonic processes on large scales. Additionally, we compare our results to a hydrodynamical N-body simulation at many redshifts and find that the counterterms associated with the leading corrections to dark matter and baryons start to differ between redshifts z ≈ 3 and z ≈ 2, signaling the onset of star-formation physics. We then use these fits to compute the lensing power spectrum, show that the understanding of baryonic processes will be important for analyzing CMB-S4 data, and show that the two-loop EFTofLSS accurately captures these effects for ℓ ≲ 2000. Our results are also potentially of interest for current and future weak lensing surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Polyhedral Distortions and Unusual Magnetic Order in Spinel FeMn 2 O 4

Spinel compounds AB 2 X 4 consist of both tetrahedral (AX 4 ) and octahedral (BX 6) environments with the former forming a diamond lattice and the latter a geometrically frustrated pyrochlore lattice. Exploring the fascinating physical properties and their correlations with structural features is critical in understanding these materials. FeMn 2 O 4 has been reported to exhibit one structural transition and two successive magnetic transitions. In this work, we report the polyhedral distortions and their correlations to the structural and two magnetic transitions in FeMn 2 O 4 by employing the high-resolution neutron powder diffraction. The cation distribution is found to be ($Mn_{0.9}^{2+}Fe_{0.1}^{3+}$)$_{A}$($Mn^{3+}Fe_{0.9}^{3+}Mn_{0.1}^{2+}$)$_{B}$O 4 . While large trigonal distortion is found even in the high-temperature cubic phase, the first-order cubic-tetragonal structural transition associated with the elongation of both tetrahedra and octahedra with shared oxygen atoms along the c axis occurs at T S ≈ 750 K, driven by the Jahn–Teller effect of the orbital active B-site Mn 3+ cation. Strong magnetoelastic coupling is unveiled at T N1 ≈ 400 K as manifested by the appearance of Néel-type collinear ferrimagnetic order, an anomaly in both tetrahedral and octahedral distortions, as well as an anomalous decrease of the lattice constants c and a weak anomaly of a. Upon cooling to T N2 ≈ 65 K, it evolves to a noncollinear ferrimagnetic order accompanied by the different moments at the split magnetic sites B1 and B2. Only one-half of the B-site Mn 3+ /Fe 3+ spins, i.e., the B2-site spins in the pyrochlore lattice, are canted, which is a unique magnetic order among spinels. The canting angle between A-site and B2-site moments is ~25°, but the B1-site moment stays antiparallel to the A-site moment even at 10 K. This noncollinear order is accompanied by a modification of the O–B–O bond angles in the octahedra without significant change in lattice constants or tetrahedral/octahedral distortion parameters, indicating a distinct magnetoelastic coupling. We demonstrate distinct roles of the A-site and B-site magnetic cations in the structural and magnetic properties of FeMn 2 O 4 . Our study indicates that FeMn 2 O 4 is a wonderful platform to unveil interesting magnetic order and to investigate their correlations with polyhedral distortions and lattice.

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

The Water-Fe-Pressure dependent single-crystal elastic properties of wadsleyite: Implications for the seismic anisotropy in the upper Mantle Transition Zone

Recent seismic studies suggest an anisotropic Mantle Transition Zone (MTZ) in areas adjacent to subducted slabs. Wadsleyite is the main anisotropy contributor in the upper MTZ, therefore the interpretation of these seismic observations requires the knowledge of single-crystal elastic moduli (C ij s) and the deformation-induced lattice preferred orientation (LPO) of wadsleyite. Wadsleyite can host up to 3 wt% water in its crystal structure as point defects in the form of hydroxyl groups, however, the combined effect of water content, Fe content, and pressure on the C ij s of wadsleyite remains unclear. In this study, we measured the high-pressure single-crystal elasticity of a synthetic hydrous Fe-bearing wadsleyite (0.14(5) wt% water, Fe#=9.4, Fe 3+ /ΣFe=0.3) up to 18.2 (2) GPa. In combination with previous experimental data, we separated the effects of pressure, water, and Fe contents on the Cijs and intrinsic elastic anisotropy of wadsleyite. Overall, our results suggest that the intrinsic elastic anisotropy of wadsleyite decreases with pressure, water, and Fe contents. At 15 GPa, increasing the water content by 0.1 wt% or Fe# by 1 decreases the V P and V S anisotropy of wadsleyite by ~1.1-1.3%, and ~0.8-1.3% in average, respectively. Combined with the LPO determined in previous deformation experiments, we modeled the seismic anisotropy in the upper MTZ generated by a sub-vertical mantle flow near cold subducted slabs and a sub-horizontal mantle flow in the ambient mantle. In both scenarios, the LPO of wadsleyite leads to V SV (vertically polarized shear wave velocity) > V SH (horizontally polarized shear wave velocity). Our results suggest that wadsleyite may account for a weak anisotropic MTZ (<1%) on the global scale. Considering the fact that water decreases the elastic anisotropy but promotes LPO of wadsleyite, seismic anisotropy may not be a good water sensor.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The State of the Molecular Gas in Post-starburst Galaxies

The molecular gas in galaxies traces both the fuel for star formation and the processes that can enhance or suppress star formation. Observations of the molecular gas state can thus point to when and why galaxies stop forming stars. In this study, we present Atacama Large Millimeter/submillimeter Array observations of the molecular gas in galaxies evolving through the post-starburst phase. These galaxies have low current star formation rates (SFRs), regardless of the SFR tracer used, with recent starbursts ending within the last 600 Myr. We present CO (3–2) observations for three post-starburst galaxies, and dense gas HCN/HCO+/HNC (1–0) observations for six (four new) post-starburst galaxies. The post-starbursts have low excitation traced by the CO spectral-line energy distribution up to CO (3–2), more similar to early-type than starburst galaxies. The low excitation indicates that lower density rather than high temperatures may suppress star formation during the post-starburst phase. One galaxy displays a blueshifted outflow traced by CO (3–2). MaNGA observations show that the ionized gas velocity is disturbed relative to the stellar velocity field, with a blueshifted component aligned with the molecular gas outflow, suggestive of a multiphase outflow. Low ratios of HCO + /CO, indicating low fractions of dense molecular gas relative to the total molecular gas, are seen throughout post-starburst phase, except for the youngest post-starburst galaxy considered here. These observations indicate that the impact of any feedback or quenching processes may be limited to low excitation and weak outflows in the cold molecular gas during the post-starburst phase.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Neutron Skin Measurement of208Pb Using Parity-Violating Electron Scattering

Nuclear stability requires large nuclei to have more neutrons than protons. While the central core of the nucleus is composed of symmetric nuclear matter, the extra neutrons in heavy nuclei are pushed out to the surface of the nucleus, forming a pure outer neutron skin layer. While the proton distribution inside a nucleus has been accurately measured using electromagnetic probes, previous measurements of neutron distributions in complex nuclei have suffered a lack of systematic precision. PREX-2 and CREX use Parity-Violating electron Scattering (PVeS), as an electroweak interaction probe, to measure the neutron distributions in208Pb and48Ca, respectively. Making use of elastic scattering of longitudinally polarized electrons from unpolarized (isotopically pure) targets, PREX-2 and CREX measure the asymmetry in scattering rates between opposite beam polarization states. Since both experiments share the same exact measurement concept and apparatus, much of this dissertation applies to both PREX-2 and CREX. However, the presented asymmetry analysis and neutron skin extraction are only for PREX-2. PREX-2 measures the parity-violating asymmetry, APV, for 953 MeV electrons scattered elastically from208Pb at?5°in the lab. The final result is APV= 550.00±16.09(stat.)±8.16(syst.)ppb. From this measurement, we report the neutron skin thickness in208Pb nuclei to be, Rskin= 0.278±0.078(exp.)±0.012(model) fm. This measurement has broad implications throughout the physics community, increasing our knowledge in neutron star structure, the equation of state of nuclear matter, nuclear baryon density, nuclear electric dipole polarizability, and more. Key Words: Neutron skin, electroweak interaction, PREX-2, CREX, parity-violating asymmetry, neutron star, quartz detector, weak charge radius, symmetry energy, Hall A, Jefferson Lab

Adhikari, Devi↗
DOE OSTI2024

Evaluation of Global Fire Simulations in CMIP6 Earth System Models

Fire is the primary form of terrestrial ecosystem disturbance on a global scale and an important Earth system process. Most Earth system models (ESMs) have incorporated fire modeling, with 19 of them submitting model outputs of fire-related variables to the Coupled Model Intercomparison Project Phase 6 (CMIP6). This study provides the first comprehensive evaluation of CMIP6 historical fire simulations by comparing them with multiple satellite-based products and charcoal-based historical reconstructions. Our results show that most CMIP6 models simulate the present-day global burned area and fire carbon emissions within the range of satellite-based products. They also capture the major features of observed spatial patterns and seasonal cycles, the relationship of fires with precipitation and population density, and the influence of the El Niño–Southern Oscillation (ENSO) on the interannual variability of tropical fires. Regional fire carbon emissions simulated by the CMIP6 models from 1850 to 2010 generally align with the charcoal-based reconstructions, although there are regional mismatches, such as in southern South America and eastern temperate North America prior to the 1910s and in temperate North America, eastern boreal North America, Europe, and boreal Asia since the 1980s. The CMIP6 simulations have addressed three critical issues identified in CMIP5: (1) the simulated global burned area being less than half of that of the observations, (2) the failure to reproduce the high burned area fraction observed in Africa, and (3) the weak fire seasonal variability. Furthermore, the CMIP6 models exhibit improved accuracy in capturing the observed relationship between fires and both climatic and socioeconomic drivers and better align with the historical long-term trends indicated by charcoal-based reconstructions in most regions worldwide. However, the CMIP6 models still fail to reproduce the decline in global burned area and fire carbon emissions observed over the past 2 decades, mainly attributed to an underestimation of anthropogenic fire suppression, and the spring peak in fires in the Northern Hemisphere midlatitudes, mainly due to an underestimation of crop fires. In addition, the model underestimates the fire sensitivity to wet–dry conditions, indicating the need to improve fuel wet-ness estimation. Based on these findings, we present specific guidance for fire scheme development and suggest a postprocessing methodology for using CMIP6 multi-model outputs to generate reliable fire projection products.

Wildfire, Earth system models↗
DOE OSTI2021

SDSS-IV MaNGA: the physical origin of off-galaxy H α blobs in the local Universe

ABSTRACT H α blobs are off-galaxy emission-line regions with weak or no optical counterparts. They are mostly visible in H α line, appearing as concentrated blobs. Such unusual objects have been rarely observed and studied, and their physical origin is still unclear. We have identified 13 H α blobs in the public data of MaNGA survey, by visually inspecting both the optical images and the spatially resolved maps of H α line for ∼4600 galaxy systems. Among the 13 H α blobs, 2 were reported in previously MaNGA-based studies and 11 are newly discovered. This sample, though still small in size, is by far the largest sample with both deep imaging and integral field spectroscopy. Therefore, for the first time we are able to perform statistical studies to investigate the physical origin of H α blobs. We examine the physical properties of these H α blobs and their associated galaxies, including their morphology, environments, gas-phase metallicities, kinematics of ionized gas, and ionizing sources. We find that the H α blobs in our sample can be broadly divided into two groups. One is associated with interacting/merging galaxy systems, of which the ionization is dominated by shocks or diffuse ionized gas. It is likely that these H α blobs used to be part of their nearby galaxies, but were stripped away at some point due to tidal interactions. The other group is found in gas-rich systems, appearing as low-metallicity star-forming regions that are visually detached from the main galaxy. These H α blobs could be associated with faint discs, spiral arms, or dwarf galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Science Validation for Dark Energy Research with Optical Imaging Surveys

The Universe has been expanding at an accelerating rate over the past several billion years, as though an unknown form of dark energy permeates all of space. The ultimate scientific goal of the proposed research is to distinguish between different physical mechanisms that could account for this observed accelerated expansion, for example, the zero-point energy of the vacuum, a dynamical form of energy that varies in time and/or space, or a modification to our theory of gravity. Wide- field optical imaging surveys of the night sky can test these competing models by measuring both the cosmic expansion history and the growth of large-scale structure. To this end, the Dark Energy Survey (DES) has cataloged several hundred million galaxies and thousands of supernovae. The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), will enlarge the census to billions of galaxies and hundreds of thousands of supernovae. A critical question for these dark energy experiments is whether systematic uncertainties can continue to be controlled at a level to keep pace with the statistical precision offered by such enormous datasets. The immediate research objectives of this project were (1) to prepare and validate input datasets that are the foundation of cosmological analyses with DES, and (2) to prepare for value-added characterization of Rubin Observatory commissioning data to inform early operations and accelerate the realization of dark energy science from LSST data products. For DES, we assembled and curated cosmology-ready data releases that include value-added components such as enhanced photometric and astrometric calibrations, alternative source extraction algorithms, maps of the survey coverage and survey conditions, object classifications, object quality selections, galaxy shapes, and photometric redshifts. We used the galaxy clustering technique to validate the photometric redshift distributions of various galaxy samples to be used as lenses in combined studies of galaxy clustering and weak gravitational lensing. The galaxy clustering redshift analysis was enhanced by use of a larger sample of reference galaxies from the eBOSS spectroscopic survey that extends to higher redshifts. For LSST, we prepared for science validation studies of commissioning data aimed at dark energy science capability that extend beyond the normative system-level tests to be done by the Rubin Observatory Construction Project. We identified a set of proposed survey strategies and candidate target fields that could be observed during the commissioning period to enhance science validation activities related to studies of dark energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Noise and error analysis and optimization in particle-based kinetic plasma simulations

In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix C(x,y) is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫C(x,y)dy = 0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, p(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. Next, we repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. We find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. Furthermore, we show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Unified Universal Control and Coordination of Inverter-Based Resources, and Validation for a PV + Battery Hybrid Plant

As renewable energy deployment grows, hybrid power plants (HPPs) combining photovoltaic (PV) and battery systems must evolve to offer both energy and grid stability services. These systems typically include a mix of grid-following (GFL) and grid-forming (GFM) inverters, presenting unique coordination and control challenges. This Department of Energy–funded project developed and validated a Unified Universal Control and Coordination (UUCC) framework for such PV + battery hybrid plants, enabling seamless and stable operation, including ultrafast black start, autonomous synchronization, and robust frequency and voltage regulation, under different grid conditions. The project significantly advanced the understanding of inverter-based resource (IBR) control by developing and validating three complementary system-level approaches for hybrid GFL/GFM operation: 1. A combined Virtual Resistance (VR)-based GFL and Virtual Oscillator Control (VOC)-based GFM method, where each inverter type is governed by a specialized control strategy. Together, these achieve stable, fast-response coordination, eliminating inrush current and enabling smooth black start and grid synchronization across a wide range of grid strengths. 2. A Deadbeat-based UUCC strategy, which uses discrete-time, switching-cycle-level control for both GFL and GFM inverters. This approach replaces traditional PI/PLL control with a control parameter-free, high-bandwidth framework that supports stable LVRT and instantaneous synchronization under all conditions. 3. A benchmark comparison with Siemens’ commercial GFM microgrid controller, which provided a fast baseline platform. The commercial approach decoupled v & f control was implemented on a commercial microgrid controller.The baseline commercial benchmark helped highlight superior transient response and black start performance offered by the deadbeat and VOC approaches. These technical contributions offer substantial improvements over conventional inverter control schemes, which often rely on slow phase-locked loop (PLL)-based synchronization, require careful control parameters tuning, and prone to unstable in weak grids with GFL inverters and in stiff grid with GFM inverters therefore challenging for hybrid GFL+GFM under all grid conditions. The deadbeat-based UUCC framework enables simpler, faster, and more robust operation of hybrid IBR systems using wide-bandgap (WBG) devices such as SiC power semiconductors. The rapid expansion of hybrid distributed energy resources (DERs), including residential and commercial PV-BESS installations such as Tesla Powerwall, PV with vehicle-to-grid (V2G) capability, and other integrated configurations, presents complex operational challenges for medium-voltage radial distribution feeders. These networks are subject to frequent disturbances such as faults, switching operations, rapid reclosing sequences, and feeder reconfigurations, all of which introduce dynamic stress on IBRs. In addition, planned feeder segmentation and deliberate islanding for resilience will require DERs that can autonomously perform blackstart, establish voltage and frequency references, and resynchronize with the main grid. The advanced deadbeat-based UUCC control and blackstart functionalities developed in this project directly address these requirements, enabling decentralized and autonomous operation of inverter-dominated DERs in distribution systems under a wide range of fault and reconfiguration scenarios. From a public benefit perspective, these innovations enable more reliable and cost-effective integration of renewable energy into distribution networks. The ability to autonomously black start and stabilize grids under varying grid conditions support accelerates recovery from outages and support decentralized resilient energy systems. By reducing system complexity and improving performance, this project lays critical groundwork for future inverter-dominated power grids that are clean, reliable, and accessible to all.

14 SOLAR ENERGY↗
DOE OSTI2020

Spontaneously supersaturated nucleation strategy for high reproducible and efficient perovskite solar cells

High-quality perovskite films are predominantly prepared by anti-solvent assisted spin-coating method. How- ever, small device size and narrow operation window limit their applications in practical and scalable production. Here, a spontaneously supersaturated nucleation strategy is developed to fabricate high-quality perovskite films without any antisolvents. The rapid precipitation behavior of the precursor solutions could be triggered by the volatilization of methylamine ethanol at the initial film formation stage, forming uniform and compact intermediate phase films. Mirror-like large area (10 cm × 10 cm) methylammonium lead triiodide (MAPbI 3 ) perovskite films with micrometer-sized grains can be achieved by phase transition via annealing. The average power conversion efficiency (PCE) of 19.4% is achieved for 0.1 cm 2 perovskite solar cells (PSCs). The 1 cm 2 device shows a PCE of 17.57% and retains 96.8% after aging for 39 days in the air with a relative humidity of 5–30%. The perovskite module based on the perovskite films prepared by our strategy exhibits a PCE of 13.13%. Here, the perovskite film fabrication method developed in this work shows high reproducibility of the fabrication process due to the spontaneous nucleation process. Meanwhile, it also exhibits significant potential in the application of scalable solution processing deposition techniques toward the commercialization of PSCs.

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

Time-Energy Uncertainty Relation for Noisy Quantum Metrology

Detection of very weak forces and precise measurement of time are two of the many applications of quantum metrology to science and technology. To sense an unknown physical parameter, one prepares an initial state of a probe system, allows the probe to evolve as governed by a Hamiltonian 𝐻 for some time 𝑡, and then measures the probe. If 𝐻 is known, we can estimate 𝑡 by this method; if 𝑡 is known, we can estimate classical parameters on which 𝐻 depends. The accuracy of a quantum sensor can be limited by either intrinsic quantum noise or by noise arising from the interactions of the probe with its environment. In this work, we introduce and study a fundamental trade-off, which relates the amount by which noise reduces the accuracy of a quantum clock to the amount of information about the energy of the clock that leaks to the environment. Specifically, we consider an idealized scenario in which a party Alice prepares an initial pure state of the clock, allows the clock to evolve for a time that is not precisely known, and then transmits the clock through a noisy channel to a party Bob. Meanwhile, the environment (Eve) receives any information about the clock that is lost during transmission. We prove that Bob’s loss of quantum Fisher information about the elapsed time is equal to Eve’s gain of quantum Fisher information about a complementary energy parameter. We also prove a similar, but more general, trade-off that applies when Bob and Eve wish to estimate the values of parameters associated with two noncommuting observables. We derive the necessary and sufficient conditions for the accuracy of the clock to be unaffected by the noise, which form a subset of the Knill-Laflamme error-correction conditions. A state and its local time-evolution direction, if they satisfy these conditions, are said to form a metrological code. We provide a scheme to construct metrological codes in the stabilizer formalism. We show that there are metrological codes that cannot be written as a quantum error-correcting code with similar distance in which the Hamiltonian acts as a logical operator, potentially offering new schemes for constructing states that do not lose any sensitivity upon application of a noisy channel. We discuss applications of the trade-off relation to sensing using a quantum many-body probe subject to erasure or amplitude-damping noise.

metrology↗
DOE OSTI2020

The second-generation Shifted Boundary Method and its numerical analysis

Recently, the Shifted Boundary Method (SBM) was proposed within the class of unfitted (or immersed, or embedded) finite element methods. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we present enhanced variational SBM formulations for the Poisson and Stokes problems with improved flexibility and robustness. These simplified variational forms allow to relax some of the assumptions required by the mathematical proofs of stability and convergence of earlier implementations. First, we show that these new SBM implementations can be proved asymptotically stable and convergent even without the rather restrictive assumption that the inner product between the normals to the true and surrogate boundaries is positive. Second, we show that it is not necessary to introduce a stabilization term involving the tangential derivatives of the solution at Dirichlet boundaries, therefore avoiding the calibration of an additional stabilization parameter. Finally, we prove enhanced L 2 -estimates without the cumbersome assumption – of earlier proofs – that the surrogate domain is convex. Instead we rely on a conventional assumption that the boundary of the true domain is smooth, which can also be replaced by requiring convexity of the true domain. The aforementioned improvements open the way to a more general and efficient implementation of the Shifted Boundary Method, particularly in complex three-dimensional geometries. We complement these theoretical developments with numerical experiments in two and three dimensions.

42 ENGINEERING↗