Search NASA⌕ Search

SEARCH · Search NASA

Results for “Rescaling”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Hydrostatic and Non‐Hydrostatic Baroclinic Instability in the Dynamical Core of the DOE Global Climate Model

Abstract The dynamical core of the Department of Energy global climate model is used to understand the role of non‐hydrostatic dynamics in the simulation of dry and moist baroclinic waves. To reduce computational cost, the Diabatic Acceleration and REscaling approach is adopted. Scale analysis and numerical simulations suggest that the model solution is not distorted by the change of spherical metric terms due to the change of Earth radius, neither are the inter‐scale interactions strongly altered as indicated by the analysis of spectral flux. Compared with hydrostatic simulations, the onset of baroclinic instability is delayed under non‐hydrostatic dynamics as the associated weaker vertical motions tend to increase the critical wavelength and narrow the range of waves that can be baroclinically unstable. During the development of baroclinic waves, non‐hydrostatic dynamics tends to induce vertical motions in the upper troposphere, accelerating the eastward propagation of upper‐level ridges through their impact on local vorticity tendency. These processes reduce the westward tilt of the vertical ridge axes and suppress the conversion of mean flow available potential energy to eddy kinetic energy, leading to weaker baroclinic eddies than in hydrostatic simulations. The contrasts between hydrostatic and non‐hydrostatic settings hold in supplemental experiments that vary the background flow Rossby number, the amount of water vapor content and vertical resolution. We also find that mesoscale and smaller‐scale activities can be considerably under‐represented when the vertical resolution is limited to that typically used in global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Nowcasting Earthquakes With Stochastic Simulations: Information Entropy of Earthquake Catalogs

Earthquake nowcasting has been proposed as a means of tracking the change in large earthquake potential in a seismically active area. The method was developed using observable seismic data, in which probabilities of future large earthquakes can be computed using Receiver Operating Characteristic methods. Furthermore, analysis of the Shannon information content of the earthquake catalogs has been used to show that there is information contained in the catalogs, and that it can vary in time. So an important question remains, where does the information originate? In this paper, we examine this question using stochastic simulations of earthquake catalogs. Our catalog simulations are computed using an Earthquake Rescaled Aftershock Seismicity (“ERAS”) stochastic model. This model is similar in many ways to other stochastic seismicity simulations, but has the advantage that the model has only 2 free parameters to be set, one for the aftershock (Omori-Utsu) time decay, and one for the aftershock spatial migration away from the epicenter. Generating a simulation catalog and fitting the two parameters to the observed catalog such as California takes only a few minutes of wall clock time. While clustering can arise from random, Poisson statistics, we show that significant information in the simulation catalogs arises from the “non-Poisson” power-law aftershock clustering, implying that the practice of de-clustering observed catalogs may remove information that would otherwise be useful in forecasting and nowcasting. We also show that the nowcasting method provides similar results with the ERAS model as it does with observed seismicity.

58 GEOSCIENCES↗
DOE OSTI2022

A unifying Bayesian framework for merging X-ray diffraction data

Novel X-ray methods are transforming the study of the functional dynamics of biomolecules. Key to this revolution is detection of often subtle conformational changes from diffraction data. Diffraction data contain patterns of bright spots known as reflections. To compute the electron density of a molecule, the intensity of each reflection must be estimated, and redundant observations reduced to consensus intensities. Systematic effects, however, lead to the measurement of equivalent reflections on different scales, corrupting observation of changes in electron density. Here, we present a modern Bayesian solution to this problem, which uses deep learning and variational inference to simultaneously rescale and merge reflection observations. We successfully apply this method to monochromatic and polychromatic single-crystal diffraction data, as well as serial femtosecond crystallography data. We find that this approach is applicable to the analysis of many types of diffraction experiments, while accurately and sensitively detecting subtle dynamics and anomalous scattering.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

The SRG/eROSITA All-Sky Survey: Constraints on the structure growth from cluster number counts

Recent advancements in methods used in wide-area surveys have demonstrated the reliability of the number density of galaxy clusters as a viable tool for precision cosmology. Beyond testing the current cosmological paradigm, cluster number counts can also be used to investigate the discrepancies currently affecting cosmological measurements. In particular, cosmological studies based on cosmic shear and other large-scale structure probes routinely find a value for the amplitude of the fluctuations in the universe S​ 8 = σ 8 (Ω m /0.3) 0.5 smaller than the one inferred from the primary cosmic microwave background. In this work, we investigate this tension by measuring structure evolution across cosmic time as probed by the number counts of massive halos with the first SRG/eROSITA All-Sky Survey cluster catalog in the western Galactic hemisphere, complemented with the overlapping Dark Energy Survey Year-3, Kilo-Degree Survey, and Hyper Suprime-Cam data for weak lensing mass calibration, by implementing two different parameterizations and a model-agnostic method. In the first model, we measured the cosmic linear growth index as γ = 1.19 ± 0.21, which is in tension with the standard value of γ = 0.55 but in good statistical agreement with other large-scale structure probes. The second model is a phenomenological scenario in which we rescale the linear matter power spectrum at low redshift to investigate a potential reduction of structure formation, and it provided similar results. Finally, in a third strategy, we considered a standard ΛCDM cosmology, but we separated the cluster catalog into five redshift bins, measuring the cosmological parameters in each and inferring the evolution of the structure formation, finding hints of a reduction. Interestingly, the S ​8 value inferred from the number counts of the cluster eRASS1 when we add a degree of freedom to the matter power spectrum recovers the value inferred by cosmic shear studies. The observed reduction in the growth rate or systematic uncertainties associated with various measurements may account for the discrepancy in the S 8 ​ values suggested between cosmic shear probes and eROSITA cluster number counts and Planck CMB measurements.

cosmological parameters↗
DOE OSTI2021

Modeling of q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the KSTAR tokamak

The EPEC code is employed to model the q95 windows for n = 2 and n = 1 resonant magnetic perturbation (RMP)-induced edge localized mode (ELM) suppression in typical KSTAR H-mode discharges. The plasma equilibria used in the study are derived by rescaling the experimental plasma equilibrium in KSTAR discharge #18594 measured at time t = 6450 ms. The n = 2 and n = 1 ELM suppression windows predicted by EPEC are comparatively narrow in q95 (i.e., Δq 95 ≃ 0.1) and are such that (assuming a relative pedestal pressure reduction of 15% is required to trigger ELM suppression), even in the middle of the windows, there is a threshold RMP coil current of about 1–3 kA/turn that must be exceeded before ELM suppression occurs. Furthermore, the n = 2 and n = 1 ELM suppression windows calculated by the EPEC code are consistent with the experimental observations and are also broadly similar to those recently calculated by the TM1 code [Hu et al., Phys. Plasmas 28, 052505 (2021)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Accelerated steady-state electrostatic particle-in-cell simulation of Langmuir probes

First-principles particle-in-cell (PIC) simulation is a powerful tool for understanding plasma behavior, but this power often comes at great computational expense. Artificially reducing the ion/electron mass ratio is a time-honored practice to reduce simulation costs. Usually, this is a severe approximation. However, for steady-state collisionless, electrostatic (Vlasov–Poisson) systems, the solution with reduced mass ratio can be scaled to the solution for the real mass ratio, with no approximation. This “scaled mass” method, which works with already-existing PIC codes, can reduce the computation time for a large class of electrostatic PIC simulations by the square root of the mass ratio. The particle distributions of the resulting steady state must be trivially rescaled to yield the true distributions, but the self-consistent electrostatic field is independent of the mass ratio. This method is equivalent to “numerical timestepping,” an approach that evolves electron and ion populations with different time steps. Numerical timestepping can be viewed as a special case of the speed-limited PIC (SLPIC) method, which is not restricted to steady-state phenomena. Although the scaled-mass approach is simplest, numerical timestepping and SLPIC more easily generalize to include other effects, such as collisions. The equivalence of these new approaches is demonstrated by applying them to simulate a cylindrical Langmuir probe in electron–argon plasma, speeding up simulation by two orders of magnitude. Methods such as SLPIC can therefore play an invaluable role in interpreting probe measurements by including geometric effects, collisions, secondary emission, and non-Maxwellian distributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Universality in diffusion-controlled nucleation and growth

Nucleation and growth are studied in a system that undergoes diffusion-controlled condensation under gradual changes in parameters, such as cooling. Here, it is demonstrated that when the Gibbs–Thompson effect becomes negligible, the system falls into a universal regime. i.e., the final droplet size distribution remains invariant under certain rescaling of system parameters. An approximate yet very accurate analytic form is obtained for the final droplet size distribution in this regime.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Energy conservation in real-time nuclear–electronic orbital Ehrenfest dynamics

Real-time nuclear–electronic orbital Ehrenfest (RT-NEO-Ehrenfest) dynamics methods provide a first-principles approach for describing nonadiabatic molecular processes with nuclear quantum effects. For an efficient description of proton transfer within RT-NEO-Ehrenfest dynamics, the basis function center associated with the quantum proton can be allowed to move classically. Here, this traveling proton basis (TPB) approach effectively captures proton quantum dynamics, although its energy conservation behavior is not yet fully satisfactory. Two recently proposed TPB approaches, in principle, conserve the extended energy, which includes both the system energy and the kinetic energy associated with the proton basis function center. Herein, a thermostatted TPB approach is proposed to improve the conservation of the system energy, excluding the kinetic energy associated with the proton basis function center. In this approach, the quantum proton dynamics are modulated by dynamically rescaling the proton momentum operator to maintain the system energy conservation. With the excited-state intramolecular proton transfer of o-hydroxybenzaldehyde as an example, this approach is shown to significantly improve the system energy conservation while preserving the accuracy of the quantum proton dynamics as achieved in the original TPB approach.

Ehrenfest dynamics↗
DOE OSTI2021

Laboratory earthquake forecasting: A machine learning competition

Earthquake prediction, the long-sought holy grail of earthquake science, continues to confound Earth scientists. Could we make advances by crowdsourcing, drawing from the vast knowledge and creativity of the machine learning (ML) community? We used Google’s ML competition platform, Kaggle, to engage the worldwide ML community with a competition to develop and improve data analysis approaches on a forecasting problem that uses laboratory earthquake data. The competitors were tasked with predicting the time remaining before the next earthquake of successive laboratory quake events, based on only a small portion of the laboratory seismic data. The more than 4,500 participating teams created and shared more than 400 computer programs in openly accessible notebooks. Complementing the now well-known features of seismic data that map to fault criticality in the laboratory, the winning teams employed unexpected strategies based on rescaling failure times as a fraction of the seismic cycle and comparing input distribution of training and testing data. In addition to yielding scientific insights into fault processes in the laboratory and their relation with the evolution of the statistical properties of the associated seismic data, the competition serves as a pedagogical tool for teaching ML in geophysics. The approach may provide a model for other competitions in geosciences or other domains of study to help engage the ML community on problems of significance.

58 GEOSCIENCES↗
DOE OSTI2021

Time–connectivity superposition and the gel/glass duality of weak colloidal gels

Colloidal gels result from the aggregation of Brownian particles suspended in a solvent. Gelation is induced by attractive interactions between individual particles that drive the formation of clusters, which in turn aggregate to form a space-spanning structure. We study this process in aluminosilicate colloidal gels through time-resolved structural and mechanical spectroscopy. Using the time–connectivity superposition principle a series of rapidly acquired linear viscoelastic spectra, measured throughout the gelation process by applying an exponential chirp protocol, are rescaled onto a universal master curve that spans over eight orders of magnitude in reduced frequency. This analysis reveals that the underlying relaxation time spectrum of the colloidal gel is symmetric in time with power-law tails characterized by a single exponent that is set at the gel point. The microstructural mechanical network has a dual character; at short length scales and fast times it appears glassy, whereas at longer times and larger scales it is gel-like. These results can be captured by a simple three-parameter constitutive model and demonstrate that the microstructure of a mature colloidal gel bears the residual skeleton of the original sample-spanning network that is created at the gel point. Our conclusions are confirmed by applying the same technique to another well-known colloidal gel system composed of attractive silica nanoparticles. The results illustrate the power of the time–connectivity superposition principle for this class of soft glassy materials and provide a compact description for the dichotomous viscoelastic nature of weak colloidal gels.

36 MATERIALS SCIENCE↗
DOE OSTI2025

3D pattern formation of a protein–membrane suspension

Many essential cellular processes, including cell division and the establishment of cell polarity during embryogenesis, are regulated by pattern-forming proteins. These proteins often need to bind to a substrate, such as the cell membrane, onto which they interact and form two-dimensional (2D) patterns. It is unclear how the membrane’s continuity and dimensionality impact pattern formation. Here, we address this gap using the MinDE system, a prototypical example of pattern-forming membrane proteins. We show that when the lipid substrate is fragmented into submicrometer-sized diffusive liposomes, adenosine triphosphate-driven protein–protein interactions generate three-dimensional (3D) spatially extended patterns, despite the complete loss of membrane continuity. Remarkably, these 3D patterns emerge at scales four orders of magnitude larger than the individual liposomes. By systematically varying protein concentration, liposome size, and density, we observed and characterized a variety of 3D dynamical patterns not seen on continuous 2D membranes, including traveling waves, dynamical spirals, and a coexistence phase. Simulations and linear stability analysis of a coarse-grained model revealed that the physical properties of the dispersed membrane effectively rescale both the protein–membrane binding rates and diffusion, two key parameters governing pattern formation and wavelength selection. These findings highlight the robustness of Min’s pattern-forming ability, suggesting that protein–membrane suspensions could serve as an adaptable template for studying out-of-equilibrium self-organization in 3D, beyond in vivo contexts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the stepped pressure equilibrium code (SPEC) is efficient when the domain is a hollow torus, but ill-conditioning of the discretised linear equations occurs in the (solid) toroid due to the artificially singular coordinate parameterisation near the axis. Here, we propose an improved choice for the reference axis to prevent coordinates surfaces from overlapping. Then, we examine the parity and asymptotics of the magnetic vector potential near the axis and suggest the use of recombined and rescaled Zernike radial basis functions. The maximum relative error in the magnetic field of the Wendelstein 7-X geometry is shown to reach 10 –9 at high resolution in a series of convergence tests and benchmarks against the boundary integral equation solver for Taylor states. The new method is also reported to significantly improve the accuracy of multi-volume SPEC calculations. A comparison between free-boundary SPEC and the analytical Dommaschk potential is presented with higher-than-usual Fourier resolution. It is illustrated that we are able to resolve low amplitude current sheets when an interface is placed where there is no flux surface in the analytic solution. This was previously concealed because of insufficient numerical resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A sound horizon-free measurement of H 0 in DESI 2024

The physical size of the sound horizon at recombination is a powerful source of information for early-time measurements of the Hubble constant H 0 , and many proposed solutions to the Hubble tension therefore involve modifications to this scale. In light of this, there has been growing interest in measuring H 0 independently of the sound horizon. We present the first such measurement to use data from the Dark Energy Spectroscopic Instrument (DESI), jointly analyzing the full-shape galaxy power spectra of DESI luminous red galaxies, emission line galaxies, quasars, and the bright galaxy sample, in a total of six redshift bins. Information from the sound horizon scale is removed from our constraints via a rescaling procedure at the power spectrum level, with our sound horizon-marginalized measurement being driven instead primarily by the matter-radiation equality scale. This measurement is then combined with additional sound horizon-free information from Planck+ACT CMB lensing, uncalibrated type Ia supernovae, and the DESI Lyman-α forest. We agnostically combine with the DESY5, Pantheon+, and Union3 supernova datasets, with our tightest respective constraints being H 0 = 66.7 +1.7 -1.9 , 67.9 +1.9 -2.1 , and 67.8 +2.0 -2.2 km s -1 Mpc -1 . This corresponds to a sub-3% sound horizon-free constraint of the Hubble constant, and is the most precise measurement of its kind to date. Even without including information from the sound horizon, our measurement is still in 2.2-3.0σ tension with SH0ES. Additionally, the consistency between our result and other measurements that do rely on the sound horizon scale provides no evidence for new early-Universe physics (e.g. early dark energy). Future DESI data releases will allow unprecedented measurements of H 0 and place strong constraints on models that use beyond-ΛCDM physics to ameliorate the Hubble tension.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

H 0 without the sound horizon (or supernovae): A 2% measurement in DESI DR1

The sound horizon scale rs is a key source of information for measurements of H 0 from early-time data, and is therefore a common target of new physics proposed to solve the Hubble tension. We present a sub-2% measurement of the Hubble constant that is independent of this scale, using data from the first data release of the Dark Energy Spectroscopic Instrument (DESI DR1). Building on previous work, we remove dependency on the sound horizon size using a heuristic rescaling procedure at the power spectrum level. A key innovation is the inclusion of uncalibrated (agnostic to r s ) post-reconstruction BAO measurements from DESI DR1, as well as using the CMB acoustic scale θ* as a high-redshift anchor. Uncalibrated type-Ia supernovae are often included as an independent source of Ω m information; here we demonstrate the robustness of our results by additionally considering two supernova-independent alternative datasets. We find somewhat higher values of H 0 relative to our previous work: 69.2 +1.3 -1.4 , 70.3 +1.4 -1.2 , and 69.6 +1.3 -1.8 km s -1 Mpc -1 respectively when including measurements from i) Planck/ACT CMB lensing × unWISE galaxies, ii) the DES Year 3 6×2pt analysis, and iii) Planck/ACT CMB lensing + the DES Year 5 supernova analysis. These remarkably consistent constraints achieve better than 2% precision; they are among the most stringent sound horizon-independent measurements from LSS to date, and provide a powerful avenue for probing the origin of the Hubble tension.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: N -body mock challenge for the quasar sample

ABSTRACT The growth rate and expansion history of the Universe can be measured from large galaxy redshift surveys using the Alcock–Paczynski effect. We validate the Redshift Space Distortion models used in the final analysis of the Sloan Digital Sky Survey (SDSS) extended Baryon Oscillation Spectroscopic Survey (eBOSS) Data Release 16 quasar clustering sample, in configuration and Fourier space, using a series of halo occupation distribution mock catalogues generated using the OuterRim N-body simulation. We test three models on a series of non-blind mocks, in the OuterRim cosmology, and blind mocks, which have been rescaled to new cosmologies, and investigate the effects of redshift smearing and catastrophic redshifts. We find that for the non-blind mocks, the models are able to recover fσ8 to within 3 per cent and α∥ and α⊥ to within 1 per cent. The scatter in the measurements is larger for the blind mocks, due to the assumption of an incorrect fiducial cosmology. From this mock challenge, we find that all three models perform well, with similar systematic errors on fσ8, α∥, and α⊥ at the level of $\sigma _{f\sigma _8}=0.013$, $\sigma _{\alpha _\parallel }=0.012$, and $\sigma _{\alpha _\bot }=0.008$. The systematic error on the combined consensus is $\sigma _{f\sigma _8}=0.011$, $\sigma _{\alpha _\parallel }=0.008$, and $\sigma _{\alpha _\bot }=0.005$, which is used in the final DR16 analysis. For baryon acoustic oscillation fits in configuration and Fourier space, we take conservative systematic errors of $\sigma _{\alpha _\parallel }=0.010$ and $\sigma _{\alpha _\bot }=0.007$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

OzDES Reverberation Mapping Programme: Mg ii lags and R − L relation

ABSTRACT The correlation between the broad line region radius and continuum luminosity (R–L relation) of active galactic nuclei (AGNs) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At z ∼ 1–2, where AGN activity peaks, the R–L relation is constrained by the reverberation mapping (RM) lags of the Mg ii line. We present 25 Mg ii lags from the Australian Dark Energy Survey RM project based on 6 yr of monitoring. We define quantitative criteria to select good lag measurements and verify their reliability with simulations based on both the damped random walk stochastic model and the rescaled, resampled versions of the observed light curves of local, well-measured AGN. Our sample significantly increases the number of Mg ii lags and extends the R–L relation to higher redshifts and luminosities. The relative iron line strength $\mathcal {R}_{\rm Fe}$ has little impact on the R–L relation. The best-fitting Mg iiR–L relation has a slope α = 0.39 ± 0.08 with an intrinsic scatter $\sigma _{\rm rl} = 0.15^{+0.03}_{-0.02}$ . The slope is consistent with previous measurements and shallower than the H β R–L relation. The intrinsic scatter of the new R–L relation is substantially smaller than previous studies and comparable to the intrinsic scatter of the H β R–L relation. Our new R–L relation will enable more precise single-epoch mass estimates and SMBH demographic studies at cosmic noon.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The impact of baryonic potentials on the gravothermal evolution of self-interacting dark matter haloes

The presence of a central baryonic potential can have a significant impact on the gravothermal evolution of self-interacting dark matter (SIDM) haloes. We extend a semi-analytical fluid model to incorporate the influence of a static baryonic potential and calibrate it using controlled N-body simulations. We construct benchmark scenarios with varying baryon concentrations and different SIDM models, including constant and velocity-dependent self-interacting cross-sections. The presence of the baryonic potential induces changes in SIDM halo properties, including central density, core size, and velocity dispersion, and it accelerates the halo’s evolution in both expansion and collapse phases. Furthermore, we observe a quasi-universality in the gravothermal evolution of SIDM haloes with the baryonic potential, resembling a previously known feature in the absence of the baryons. By appropriately rescaling the physical quantities that characterize the SIDM haloes, the evolution of all our benchmark cases exhibits remarkable similarity. In conclusion, our findings offer a framework for testing SIDM predictions using observations of galactic systems where baryons play a significant dynamical role.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Lattice study of correlators of chromoelectric fields for heavy quarkonium dynamics in the quark-gluon plasma

We perform a lattice calculation of the correlators of two chromoelectric fields in the adjoint representation connected by adjoint Wilson lines at nonzero temperature. These correlators arise in the study of quarkonium dynamics and of adjoint heavy quark diffusion in deconfined matter. We work in SU(3) gauge theory using either gradient flow or multilevel algorithms for noise reduction, and discuss the renormalization of the correlators on the lattice. We find that a Casimir factor rescaling relates the adjoint correlators corresponding to the diffusion of an adjoint heavy quark and the octet-octet quarkonium transitions to the chromoelectric correlator in the fundamental representation describing the diffusion of a heavy quark.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗