Results for “coupling constant” — Page 66

NASA NTRS1994

Excitation of interstellar hydrogen chloride

We have computed new rate coefficients for the collisional excitation of HCl by He, in the close-coupled formalism and using an interaction potential determined recently by Willey, Choong, & DeLucia. Results have been obtained for temperatures between 10 K and 300 K. With the use of the infinite order sudden approximation, we have derived approximate expressions of general applicability which may be used to estimate how the rate constant for a transition (J to J prime) is apportioned among the various hyperfine states F prime of the final state J prime. Using these new rate coefficients, we have obtained predictions for the HCl rotational line strengths expected from a dense clump of interstellar gas, as a function of the HCl fractional abundance. Over a wide range of HCl abundances, we have found that the line luminosities are proportional to abundance(exp 2/3), a general result which can be explained using a simple analytical approximation. Our model for the excitation of HCl within a dense molecular cloud core indicates that the J = 1 goes to 0 line strengths measured by Blake, Keene, & Phillips toward the Orion Molecular Cloud (OMC-1) imply a fractional abundance n(HCl)/n(H2) approximately 2 x 10(exp -9), a value which amounts to only approximately 0.3% of the cosmic abundance of chlorine nuclei. Given a fractional abundance of 2 x 10(exp -9), the contribution of HCl emission to the total radiative cooling of a dense clump is small. For Orion, we predict a flux approximately 10(exp -19) W/sq cm for the HCl J = 3 goes to 2 line near 159.8 micrometers, suggesting that the strength of this line could be measured using the Infrared Space Observatory.

Neufild, David A.↗

NASA NTRS2023

Changing Nature of High-Impact Snowfall Events in Eastern North America

Snowstorms cause substantial disruption in the eastern United States and Canada each winter. While reductions in annual snowfall are projected over most of this region due to anthropogenic global warming, daily snowfall extremes that have the greatest impact may not decrease in the same manner. We examine changes to two extreme snowfall metrics: the 95th percentile of daily snowfall (SF95, cm) and the number of events during which 10% of the mean annual snowfall is exceeded during a single day (TC10, events yr−1). We explore changes to these metrics in two ensembles of the fifth-generation Canadian Regional Climate Model (CRCM5), including four 0.22° ( ≈ 25 km) simulations driven by different coupled general circulation models as well as the higher-resolution (0.11°, ≈ 12 km) ClimEx ensemble, driven by 50 members of a large initial-condition ensemble of one global model. We find that while mean annual snowfall is projected to decrease over our domain, SF95 is projected to remain relatively constant, suggesting that the most extreme daily snowfalls currently observed are likely to occur even in a warmer future climate. The region of the largest TC10 values exhibits a northward shift, with a larger percentage of annual snowfall occurring during a few large events along the U.S.-Canada border. These projected changes to the nature of snowfall events may have important socioeconomic consequences in this densely populated region of North America.

heavy snow↗

DOE OSTI2020

Thermodynamics of monoclinic and tetragonal hafnium dioxide (HfO 2 ) at ambient pressure

Accurate and precise thermodynamic models of oxide compounds and their phases are important for calculating the phase stability of oxide materials. We develop and use a coupled quantum mechanical and molecular dynamics approach to create thermodynamic models of hafnia (HfO2) polymorphs from 0 K to 3000 K at ambient pressure. The approach is based on the quasi-harmonic approximation below the Debye temperature and on ab-initio molecular dynamics calculations above the Debye temperature to predict constant pressure heat capacities (C-p). A Bayesian model provides interpolated values between these regimes. As a case study, we develop thermodynamic models of monoclinic and tetragonal HfO2 polymorphs. The predicted heat capacities are in excellent agreement with experiment, and the predicted temperature of the monoclinic to tetragonal phase transition (2173 K) is in good agreement with the experimental value (2078K). These results provide a comprehensive and accurate thermodynamic model of the monoclinic and tetragonal phases of hafnia on a broad range of temperatures and can serve as input to CALPHAD assessment of multi-component hafnia-based phase diagrams.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Optimization of Solid Oxide Electrolysis Cell Systems Accounting for Long-Term Performance and Health Degradation

This study focuses on optimizing solid oxide electrolysis cell (SOEC) systems for efficient and durable long-term hydrogen (H2) production. While the elevated operating temperatures of SOECs offer advantages in terms of efficiency, they also lead to chemical degradation, which shortens cell lifespan. To address this challenge, dynamic degradation models are coupled with a steady-state, two-dimensional, non-isothermal SOEC model and steady-state auxiliary balance of plant equipment models, within the IDAES modeling and optimization framework. A quasi-steady state approach is presented to reduce model size and computational complexity. Long-term dynamic simulations at constant H2 production rate illustrate the thermal effects of chemical degradation. Dynamic optimization is used to minimize the lifetime cost of H2 production, accounting for SOEC replacement, operating, and energy expenses. Several optimized operating profiles are compared by calculating the Levelized Cost of Hydrogen (LCOH).

Giridhar, Nishant↗

DOE OSTI2020

Magnetic Sorbent for the Removal of Selenium(IV) from Simulated Industrial Wastewaters: Determination of Column Kinetic Parameters

A novel meso- and microporous tire-derived-carbon support with magnetic iron oxide nanoparticle adsorbents that selectively adsorbs Se(IV) ions from simulated contaminated water has been developed. In this work, the physicochemical characteristics of the composite adsorbent are characterized with respect to porosity and surface area, chemical composition, and microstructure morphology. The kinetics of this composite adsorbent in a fixed-bed setting has been determined. Several column runs were conducted and analyzed by inductively coupled plasma-optical emission spectroscopy (ICP-OES) to determine the concentration gradient vs time. These results were then fit to a pseudo-second order rate law to obtain equilibrium values. Combining calculated equilibrium values with effluent concentration data, enabled the application of the Adams–Bohart model to determine reaction constants and column coefficients. Column parameters obtained from different flow rates and fittings of the Adams–Bohart model were remarkably consistent. These findings enable the application of this sorbent to fixed-bed column systems and opens up further research into mixed pollutants tests with real wastewater and scaling of selenium pollutant removal.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Optimization of Solid Oxide Electrolysis Cell Systems Accounting for Long-Term Performance and Health Degradation

This study focuses on optimizing solid oxide electrolysis cell (SOEC) systems for efficient and durable long-term hydrogen (H2) production. While the elevated operating temperatures of SOECs offer advantages in terms of efficiency, they also lead to chemical degradation, which shortens cell lifespan. To address this challenge, dynamic degradation models are coupled with a steady-state, two-dimensional, non-isothermal SOEC model and steady-state auxiliary balance of plant equipment models, within the IDAES modeling and optimization framework. A quasi-steady state approach is presented to reduce model size and computational complexity. Long-term dynamic simulations at constant H2 production rate illustrate the thermal effects of chemical degradation. Dynamic optimization is used to minimize the lifetime cost of H2 production, accounting for SOEC replacement, operating, and energy expenses. Several optimized operating profiles are compared by calculating the Levelized Cost of Hydrogen (LCOH).

Giridhar, Nishant↗

NASA NTRS1988

Cosmic-ray-modified stellar winds. II - A perturbation approach

A perturbation method is developed to describe the modification of a stellar wind flow with termination shock by the Galactic cosmic rays. The perturbation parameter epsilon used in the analysis is the ratio of the galactic cosmic-ray pressure to the thermal gas pressure at a large distance from the star. The positive galactic cosmic-ray pressure gradient serves to brake the outflowing stellar wind gas, particularly just upstream of the termination shock of the wind. A one-fluid polytropic model is used to describe the thermal stellar wind gas, with the coupling between the cosmic rays and the thermal gas being determined by an average hydrodynamical diffusion coefficient kappa describing the scattering of cosmic rays by magnetohydrodynamic waves traveling in the background flow. The analysis takes into account cosmic-ray modifications of the critical point of the wind, the thermal gas entropy constants on both sides of the shock, and the fluid velocity profile.

Ko, C. M.↗

NASA NTRS2016

Characterization of the High-Albedo NEA 3691 Bede

Characterization of NEAs provides important inputs to models for atmospheric entry, risk assessment and mitigation. Diameter is a key parameter because diameter translates to kinetic energy in atmospheric entry. Diameters can be derived from the absolute magnitude, H(PA=0deg), and from thermal modeling of observed IR fluxes. For both methods, the albedo (pv) is important - high pv surfaces have cooler temperatures, larger diameters for a given Hmag, and shallower phase curves (larger slope parameter G). Thermal model parameters are coupled, however, so that a higher thermal inertia also results in a cooler surface temperature. Multiple parameters contribute to constraining the diameter. Observations made at multiple observing geometries can contribute to understanding the relationships between and potentially breaking some of the degeneracies between parameters. We present data and analyses on NEA 3691 Bede with the aim of best constraining the diameter and pv from a combination of thermal modeling and light curve analyses. We employ our UKIRT+Michelle mid-IR photometric observations of 3691 Bede's thermal emission at 2 phase angles (27&43 deg 2015-03-19 & 04-13), in addition to WISE data (33deg 2010-05-27, Mainzer+2011). Observing geometries differ by solar phase angles and by moderate changes in heliocentric distance (e.g., further distances produce somewhat cooler surface temperatures). With the NEATM model and for a constant IR beaming parameter (eta=constant), there is a family of solutions for (diameter, pv, G, eta) where G is the slope parameter from the H-G Relation. NEATM models employing Pravec+2012's choice of G=0.43, produce D=1.8 km and pv≈0.4, given that G=0.43 is assumed from studies of main belt asteroids (Warner+2009). We present an analysis of the light curve of 3691 Bede to constrain G from observations. We also investigate fitting thermophysical models (TPM, Rozitis+11) to constrain the coupled parameters of thermal inertia (Gamma) and surface roughness, which in turn affect diameter and pv. Surface composition can be related to pv. This study focuses on understanding and characterizing the dependency of parameters with the aim of constraining diameter, pv and thermal inertia for 3691 Bede.

Wooden, Diane H.↗

NASA NTRS2016

Characterization of the High-Albedo NEA 3691 Bede

Characterization of NEAs provides important inputs to models for atmospheric entry, risk assessment and mitigation. Diameter is a key parameter because diameter translates to kinetic energy in atmospheric entry. Diameters can be derived from the absolute magnitude, H(PA=0deg), and from thermal modeling of observed IR fluxes. For both methods, the albedo (pv) is important - high pv surfaces have cooler temperatures, larger diameters for a given Hmag, and shallower phase curves (larger slope parameter G). Thermal model parameters are coupled, however, so that a higher thermal inertia also results in a cooler surface temperature. Multiple parameters contribute to constraining the diameter. Observations made at multiple observing geometries can contribute to understanding the relationships between and potentially breaking some of the degeneracies between parameters. We present data and analyses on NEA 3691 Bede with the aim of best constraining the diameter and pv from a combination of thermal modeling and light curve analyses. We employ our UKIRT+Michelle mid-IR photometric observations of 3691 Bede's thermal emission at 2 phase angles (27&43 deg 2015-03-19 & 04-13), in addition to WISE data (33deg 2010-05-27, Mainzer+2011). Observing geometries differ by solar phase angles and by moderate changes in heliocentric distance (e.g., further distances produce somewhat cooler surface temperatures). With the NEATM model and for a constant IR beaming parameter (eta=constant), there is a family of solutions for (diameter, pv, G, eta) where G is the slope parameter from the H-G Relation. NEATM models employing Pravec+2012's choice of G=0.43, produce D=1.8 km and pv≈0.4, given that G=0.43 is assumed from studies of main belt asteroids (Warner+2009). We present an analysis of the light curve of 3691 Bede to constrain G from observations. We also investigate fitting thermophysical models (TPM, Rozitis+11) to constrain the coupled parameters of thermal inertia (Gamma) and surface roughness, which in turn affect diameter and pv. Surface composition can be related to pv. This study focuses on understanding and characterizing the dependency of parameters with the aim of constraining diameter, pv and thermal inertia for 3691 Bede.

magnitude↗

DOE OSTI2020

An error-controlled adaptive time-stepping method for particle advancement in coupled CFD-DEM simulations

Coupled Computational-Fluid-Dynamics (CFD) and Discrete-Element-Method (DEM) models provide an accurate description of multiphase physical systems where a solid granular particle phase exists in an underlying gaseous continuous medium. The time integration of the granular phase in these simulations is typically handled using an explicit scheme with a constant time-step among all particles that is invariant in time to resolve inter-particle collisions. A locally third-order accurate adaptive time integration technique for particles that employs an embedded locally second-order scheme for error determination is presented in this work. The particle time-step size is dynamically adapted based on solution error, thus leading to significant savings in computational time. The efficacy of our scheme is quantified using four test cases of varying complexity (binary collision, homogeneous cooling system, fluidized bed and hopper discharge). The adaptive time-stepping method exhibits improved performance (~ 2–3 times in most of the cases studied) compared to three commonly used non-adaptive time-step methods (first-order Euler-explicit, second-order Adams-Bashforth and third-order Runge-Kutta schemes), while maintaining the same level of accuracy and parallel scalability.

42 ENGINEERING↗

DOE OSTI2025

Characterization and thermometry of dissipatively stabilized steady states

In this work we study the properties of dissipatively stabilized steady states of noisy quantum algorithms, exploring the extent to which they can be well approximated as thermal distributions, and proposing methods to extract the effective temperature T. We study an algorithm called the relaxational quantum eigensolver (RQE), which is one of a family of algorithms that attempt to find ground states and balance error in noisy quantum devices. In RQE, we weakly couple a second register of auxiliary ‘shadow’ qubits to the primary system in Trotterized evolution, thus engineering an approximate zero-temperature bath by periodically resetting the auxiliary qubits during the algorithm’s runtime. Balancing the infinite temperature bath of random gate error, RQE returns states with an average energy equal to a constant fraction of the ground state. We probe the steady states of this algorithm for a range of base error rates, using several methods for estimating both T and deviations from thermal behavior. In particular, we both confirm that the steady states of these systems are often well-approximated by thermal distributions, and show that the same resources used for cooling can be adopted for thermometry, yielding a fairly reliable measure of the temperature. These methods could be readily implemented in near-term quantum hardware, and for stabilizing and probing Hamiltonians where simulating approximate thermal states is hard for classical computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Interplay between magnetism and band topology in the kagome magnets R Mn 6 Sn 6

Kagome-lattice magnets R Mn 6 Sn 6 recently emerged as a new platform to exploit the interplay between magnetism and topological electronic states. Some of the most exciting features of this family are the dramatic dependence of the easy magnetization direction on the rare-earth specie, despite other magnetic and electronic properties being essentially unchanged, and the kagome geometry of the Mn planes that in principle can generate flat bands and Dirac points; gapping of the Dirac points by spin-orbit coupling has been suggested recently to be responsible for the observed anomalous Hall response in the member TbMn 6 Sn 6 . In this paper, we address both issues with ab initio calculations. We have discovered the significant role played by higher-order crystal-field parameters and rare-earth magnetic anisotropy constants in these systems. Further, we demonstrate that the microscopic origin of rare-earth magnetic anisotropy can also be quantified and understood at various levels: ab initio, phenomenological, and analytical. In particular, using a simple and physically transparent analytical model based on perturbation theory, we are able to explain, with full quantitative agreement, the evolution of rare-earth magnetic anisotropy across the series. We analyze in detail the topological properties of Mn-dominated bands and demonstrate how they emerge from the multiorbital planar kagome model. We further show that, despite this fact, most of the topological features at the Brillouin zone corner K are strongly 3D and therefore cannot explain the observed quasi-2D anomalous Hall effect, while the most pronounced quasi-2D dispersion are too far removed from the Fermi level. By employing self-consistent calculations with ab initio many-body approaches, we demonstrate that the exchange-correlation effects beyond the density functional theory for itinerant Mn- d electrons do not significantly alter the obtained electronic and magnetic structure. Therefore, we conclude that, contrary to previous claims, the most pronounced 2D kagome-derived topological band features bear little relevance to transport in RMn 6 Sn 6 , albeit they may possibly be brought to focus by electron or hole doping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Cross-Cap Defects and Fault-Tolerant Logical Gates in the Surface Code and the Honeycomb Floquet Code

We consider the Z 2 toric code, surface code, and Floquet code defined on a nonorientable surface, which can be considered as families of codes extending Shor’s nine-qubit code. We investigate the fault-tolerant logical gates of the Z 2 toric code in this setup, which corresponds to e ↔ m exchanging symmetry of the underlying Z 2 gauge theory. We find that nonorientable geometry provides a new way for the emergent symmetry to act on the code space, and discover the new realization of the fault-tolerant Hadamard gate of the two-dimensional surface code with a single cross cap connecting the vertices nonlocally along a slit, dubbed a nonorientable surface code. This Hadamard gate can be realized by a constant-depth local unitary circuit modulo nonlocality caused by a cross cap. Via folding, the nonorientable surface code can be turned into a bilayer local quantum code, where the folded cross cap is equivalent to a bilayer twist terminated on a gapped boundary and the logical Hadamard only contains local gates with intralayer couplings when being away from the cross cap, as opposed to the interlayer couplings on each site needed in the case of the folded surface code. We further obtain the complete logical Clifford gate set for a stack of nonorientable surface codes and similarly for codes defined on Klein-bottle geometries. We then construct the honeycomb Floquet code in the presence of a single cross cap, and find that the period of the sequential Pauli measurements acts as a H Z logical gate on the single logical qubit, where the cross cap enriches the dynamics compared with the orientable case. We find that the dynamics of the honeycomb Floquet code is precisely described by a condensation operator of the Z 2 gauge theory, and illustrate the exotic dynamics of our code in terms of a condensation operator supported at a nonorientable surface. Published by the American Physical Society 2024

Physics↗

DOE OSTI2025

Covariance-Free Bifidelity Control Variates Importance Sampling for Rare Event Reliability Analysis

Multifidelity modeling has been steadily gaining attention as a tool to address the problem of exorbitant model evaluation costs that makes the estimation of failure probabilities a significant computational challenge for complex real-world problems, particularly when failure is a rare event. To implement multifidelity modeling, estimators that efficiently combine information from multiple models/sources are necessary. In past works, the variance reduction techniques of control variates (CV) and importance sampling (IS) have been leveraged for this task. In this paper, we present the CVIS framework—a creative take on a coupled CV and IS estimator for bifidelity reliability analysis. The framework addresses some of the practical challenges of the CV method by using an estimator for the control variate mean and sidestepping the need to estimate the covariance between the original estimator and the control variate through a clever choice for the tuning constant. Furthermore, the task of selecting an efficient IS distribution is also considered, with a view towards maximally leveraging the bifidelity structure and maintaining expressivity. Additionally, a diagnostic is provided that indicates both the efficiency of the algorithm as well as the relative predictive quality of the models utilized. Finally, the behavior and performance of the framework is explored through analytical and numerical examples.

Markov chain Monte Carlo↗

DOE OSTI2021

Collective Phonon–Polaritonic Modes in Silicon Carbide Subarrays

Localized surface phonon polaritons (LSPhPs) can be implemented to engineer light–matter interactions through nanoscale patterning for a range of midinfrared application spaces. However, the polar material systems studied to date have mainly focused on simple designs featuring a single element in the periodic unit cell. Increasing the complexity of the unit cell can serve to modify the resonant near-fields and intra- and inter-unit-cell coupling as well as to dictate spectral tuning in the far-field. In this work, we exploit more complicated unit-cell structures to realize LSPhP modes with additional degrees of design freedom, which are largely unexplored. Collectively excited LSPhP modes with distinctly symmetric and antisymmetric near-fields are supported in these subarray designs, which are based on nanopillars that are scaled by the number of subarray elements to ensure a constant unit-cell size. Moreover, we observe an anomalous mode-matching of the collective symmetric mode in our fabricated subarrays that is robust to changing numbers of pillars within the subarrays as well as to defects intentionally introduced in the form of missing pillars. This work therefore illustrates the hierarchical design of tailored LSPhP resonances and modal near-field profiles simultaneously for a variety of IR applications such as surface-enhanced spectroscopies and biochemical sensing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1980

A dual charge-coupled device /CCD/, astronomical spectrometer and direct imaging camera. I - Optical and detector systems

The MASCOT (MIT Astronomical Spectrometer/Camera for Optical Telescopes), an instrument capable of simultaneously performing both direct imaging and spectrometry of faint objects, is examined. An optical layout is given of the instrument which uses two CCD's mounted on the same temperature regulated detector block. Two sources of noise on the signal are discussed: (1) the CCD readout noise, which results in a constant uncertainty in the number of electrons collected from each pixel; and (2) the photon counting noise. The sensitivity of the device is limited by the sky brightness, the overall quantum efficiency, the resolution, and the readout noise of the CCD. Therefore, total system efficiency is calculated at about 15%.

Meyer, S. S.↗

DOE OSTI2022

Pulsating one-dimensional detonation in ammonia-hydrogen–air mixtures

In this study, one-dimensional detonations in ammonia/hydrogen-air mixtures are numerically investigated by solving the fully compressible Navier-Stokes equations with detailed chemistry. Pulsating instabilities with single-mode are observed during the detonation wave propagation, accompanied by periodic coupling and decoupling of the lead shock wave and the reactive front. The ratio between driver pressure and initial pressure determines the overdrive degree and thus the oscillatory mode of detonation for a premixture with certain composition. The effects of hydrogen dilution and mixture equivalence ratio on pulsating detonations are also examined under a constant driver pressure. The growing hydrogen fraction in fuel blends significantly increases the oscillation frequency. In addition, the pulsating detonation frequency rises with increasing equivalence ratio under fuel-lean conditions, peaks under stoichiometric conditions, and falls under fuel-rich conditions as the equivalence ratio increases further. Evolutions of reactants, main intermediate radicals, and products are analysed in both fuel-lean and fuel-rich conditions. A chemical explosive mode analysis further confirms the highly-autoignitive nature of the mixture in the induction zone between reaction front and shock front where thermal diffusion plays a negligible role.

08 HYDROGEN↗

DOE OSTI2023

Impact of fusion reactor neutronics modeling for transmutation and thermal feedback *

Abstract Fusion neutronics calculations provide important metrics pertinent to fusion device operations, such as tritium breeding ratios (TBRs) and data on heat deposition, material activation, and damage. Because of the high computational burden required to generate a high-fidelity Monte Carlo simulation of a 3D fusion device, various assumptions are made to reduce computational time by simplifying the reactor model or the calculation iteration. This paper explores the impact of fusion neutronics metrics such as the TBR and decay heat of structural materials based on assumptions of material composition in the fusion reactor and temperature modeling of materials. Results show that for compact tokamaks with high power and long operational cycles, the transmutation of structural materials is significant enough to cause a substantial change in the flux spectrum and decrease the TBR by 1.68% after 2 years of full power operation. Additionally, assuming a constant temperature and material density can impact the TBR calculations up to 3%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗