Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time Constant”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Estimating time in quantum chaotic systems and black holes

We characterize new universal features of the dynamics of chaotic quantum many-body systems, by considering a hypothetical task of "time estimation". Most macroscopic observables in a chaotic system equilibrate to nearly constant late-time values. Intuitively, it should become increasingly difficult to estimate the precise value of time by making measurements on the state. We use a quantity called the Fisher information from quantum metrology to quantify the minimum uncertainty in estimating time. Due to unitarity, the uncertainty in the time estimate does not grow with time if we have access to optimal measurements on the full system. Restricting the measurements to act on a small subsystem or to have low computational complexity leads to results expected from equilibration, where the time uncertainty becomes large at late times. With optimal measurements on a subsystem larger than half of the system, we regain the ability to estimate the time very precisely, even at late times. Hawking's calculation for the reduced density matrix of the black hole radiation in semiclassical gravity contradicts our general predictions for unitary quantum chaotic systems. Hawking's state always has a large uncertainty for attempts to estimate the time using the radiation, whereas our general results imply that the uncertainty should become small after the Page time. This gives a new version of the black hole information loss paradox in terms of the time estimation task. By restricting to simple measurements on the radiation, the time uncertainty becomes large. This indicates from a new perspective that the observations of computationally bounded agents are consistent with the semiclassical effective description of gravity.

Black holes↗
DOE OSTI2025

Huge ensembles – Part 1: Design of ensemble weather forecasts using spherical Fourier neural operators

Abstract. Simulating low-likelihood high-impact extreme weather events in a warming world is a significant and challenging task for current ensemble forecasting systems. While these systems presently use up to 100 members, larger ensembles could enrich the sampling of internal variability. They may capture the long tails associated with climate hazards better than traditional ensemble sizes. Due to computational constraints, it is infeasible to generate huge ensembles (comprised of 1000–10 000 members) with traditional, physics-based numerical models. In this two-part paper, we replace traditional numerical simulations with machine learning (ML) to generate hindcasts of huge ensembles. In Part 1, we construct an ensemble weather forecasting system based on spherical Fourier neural operators (SFNOs), and we discuss important design decisions for constructing such an ensemble. The ensemble represents model uncertainty through perturbed-parameter techniques, and it represents initial condition uncertainty through bred vectors, which sample the fastest-growing modes of the forecast. Using the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS) as a baseline, we develop an evaluation pipeline composed of mean, spectral, and extreme diagnostics. With large-scale, distributed SFNOs with 1.1 billion learned parameters, we achieve calibrated probabilistic forecasts. As the trajectories of the individual members diverge, the ML ensemble mean spectra degrade with lead time, consistent with physical expectations. However, the individual ensemble members' spectra stay constant with lead time. Therefore, these members simulate realistic weather states during the rollout, and the ML ensemble passes a crucial spectral test in the literature. The IFS and ML ensembles have similar extreme forecast indices, and we show that the ML extreme weather forecasts are reliable and discriminating. These diagnostics ensure that the ensemble can reliably simulate the time evolution of the atmosphere, including low-likelihood high-impact extremes. In Part 2, we generate a huge ensemble initialized each day in summer 2023, and we characterize the simulations of extremes.

Mahesh, Ankur↗
DOE OSTI2024

Understanding the Origin of Negative Temperature Dependence and Activity of N-Coordinated Cobalt Sites During Ethylene Dimerization

The on-demand production of short-chain linear alpha olefins (LAOs; C4-C8) via C2H4 dimerization and oligomerization is industrially attractive, prompting extensive research on designing active, selective, and stable catalysts for industrial use. Cobalt supported on ammoniated carbon (Co(NH3)x/C) catalysts have shown remarkable activity and selectivity in this process. However, critical aspects such as the active phase, active site structure, the role of the catalyst support, cobalt loading effects, and the inverse correlation of the reaction rate with temperature remain inadequately understood. This study systematically explores these factors using a combination of steady-state differential catalytic tests, in situ molecular characterization including diffuse reflectance UV-Vis (DR-UV-Vis), Infrared, and Raman spectroscopies, and ex situ X-ray diffraction (XRD) and high annular aberration-corrected dark field transmission electron microscopy (HAADF-STEM). Various supports (SiO2, Al2O3, NH4-ZSM-5, g-C3N4, and C) and cobalt loadings (1.0-3.0 Co nm-2) were studied to determine the optimal catalyst composition and identify the active phase and sites. Carbon-supported catalysts uniquely produce C4-8 LAOs during C2H4 dimerization, with site-time-yield remaining constant (~10-3 s-1) for 1.0-4.0 Co nm-2 at prolonged reaction times (24-48?h time-on-stream). At higher loadings of 6.0 Co nm-2, the formation of crystalline CoO and Co3O4 phases reduces catalytic activity and LAO selectivity. Our findings show that active catalysts lack crystalline cobalt oxides and instead feature dispersed Co2+ sites, tetra-coordinated to a mix of N/NH3 and O/H2O ligands, which catalyze C2H4 dimerization via the Cossee-Arlman mechanism, exhibiting 1st order dependence on C2H4 concentration. The observed inverse rate-temperature correlation is attributed to compensation effects (i.e., presence of Cremer-Constable relationship) linked to changes in adsorption enthalpic and entropic factors.

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

Effects of Critical Compression Ratio on Rating Gasoline Knock Propensity

It is common practice in the automotive industry to explore the knock limits of fuels on an engine by a comparison of the knock limited spark advance (KLSA) at threshold knock intensity. However, the knock propensity of gasolines can be rated by changing one of three metrics on a variable compression ratio Cooperative Fuels Research (CFR) octane rating engine while holding the other two variables constant: knock intensity, spark timing, and critical compression ratio. The operational differences between the standard research octane number (RON) rating and modern engine operation have been explored in three parts. The first part focused on the effects of lambda and knock characterization. The second part studied the effects of spark timing. This third part explores the knock ratings of several gasolines by comparing the critical compression ratios at constant combustion phasing and knock intensity. The threshold knock intensity was based on the standard octane rating D1 pickup or by maximum amplitude of pressure oscillations (MAPO) measured by a piezoelectric cylinder pressure transducer. Several Fuels for Advanced Combustion Engines (FACE) gasolines, primary reference fuels (PRFs), and toluene standardization fuels (TSFs) were tested on a CFR octane rating engine with advanced data acquisition equipment and a piezoelectric cylinder pressure transducer. These tests deviated from the ASTM D2699 standard octane rating procedure. For each test fuel, the CFR engine was operated at stoichiometry at a constant combustion phasing (CA50) and the compression ratio was modified until a threshold knock intensity was realized. It was found that the chemical composition of the fuels affected the relationship of critical compression ratios between the D1 knockmeter and piezoelectric pressure transducer knock intensity thresholds, as well as the measured combustion maximum pressure rise rate and spark timing setting for constant CA50. For highly aromatic fuels tested at a constant MAPO knock intensity threshold, it was found that the maximum pressure rise rate was two to three times higher than that of highly paraffinic fuels with similar RON and the spark advance was several crank angle degrees less for constant combustion phasing.

Kolodziej, Christopher P↗
DOE OSTI2026

Constrained variational optimization of counting-time allocation in sequential scattering measurements: Application to Bonse–Hart USANS

Sequential scattering measurements are often performed under a fixed experimental-time budget, even though the expected count rate varies strongly across the measured coordinate. When the dwell time at each measurement position can be controlled independently, this variation creates a general resource-allocation problem: how should the available time be distributed to minimize the uncertainty of the reconstructed profile? We formulate this problem as a constrained variational optimization for measurements governed by Poisson counting statistics. When each measurement is treated independently, minimizing the averaged squared relative uncertainty yields an inverse-square-root intensity allocation. The formulation is then generalized to include correlations between neighboring measurements and an instrumental resolution operator, leading to an allocation criterion that equalizes the marginal reduction in posterior uncertainty per unit measurement time. Bonse–Hart ultra-small-angle neutron scattering (USANS), in which reciprocal space is sampled sequentially through analyzer-angle stepping, provides an experimentally grounded application. Computational benchmarking shows that the optimized allocation outperforms uniform-time and constant-relative-error strategies, while application to an experimentally measured graphite USANS profile from the Spallation Neutron Source, using Poisson resampling under alternative schedules, demonstrates how counting time should be redistributed toward weak-intensity regions under an identical total duration. The resulting framework applies to sequential scattering and related scanning measurements whenever local dwell times are adjustable and directly determine the measurement uncertainties, and when the relevant correlation and instrumental-response models are available.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗
DOE OSTI2026

Rotational coherence dominates early-time dynamics and produces long-time revivals in the S 2 state of azulene

Here, the ultrafast dynamics of azulene have been debated for decades, with reported picosecond decay constants variously attributed to intramolecular vibrational redistribution (IVR), internal conversion, or rotational dephasing. Using polarization- and femtosecond time-resolved resonance-enhanced multiphoton ionization spectroscopy with a nanosecond delay window, we disentangle this long-standing inconsistency and show that the early 2–5 ps decay component arises entirely from the rotational dephasing of an excited-state wavepacket. Identical time constants extracted from the decay of the parallel signal and the rise of the perpendicular signal across multiple vibronic origins provide an unambiguous rotational anisotropy signature, eliminating the need for IVR-based interpretations. Extending the measurement window to 1.3 ns reveals well-structured J-type and C-type rotational coherence revivals in S 2 azulene on top of the well-documented fluorescence decay, demonstrating that both the short- and long-time dynamics contain information about the coherent rotational dynamics. These results show that azulene, and by extension polycyclic aromatic hydrocarbons, can sustain structured rotational coherence deep into the nanosecond regime, positioning PAHs as model systems for quantum-coherent wavepacket dynamics and providing a framework for probing coherence, decoherence, and rotational control in electronically rich molecular systems.

Zhan, Jie [University of Georgia, Athens, GA (Unit↗
DOE OSTI2024

Reconstruction of the BNB and NuMI Neutrino Bunch Structure with ICARUS

ICARUS serves as the Far Detector of the Short Baseline Neutrino (SBN) program at Fermilab, sitting on-axis on the Booster Neutrino Beam (BNB) and 6$^\circ$ off-axis from the Neutrinos at the Main Injector (NuMI) beam. Neutrinos from both beams inherit the timing sub-structure of their parent proton spills, which is in turn derived from either the Booster's or the Main Injector's synchrotron acceleration. Since neutrino propagation introduces only a constant offset, their timing structure is preserved as they travel. Identifying this structure in data represents a powerful tool for selecting neutrino events and searching for physics beyond the Standard Model (BSM). This poster presents the preliminary reconstruction of the BNB and NuMI neutrino bunch structure with ICARUS data, exploiting only the precise timing of ICARUS optical readout system to both locate and assign a time to each interaction.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Reconstruction of the BNB and NuMI Neutrino Bunch Structure with ICARUS

ICARUS serves as the Far Detector of the Short Baseline Neutrino (SBN) program at Fermilab, sitting on-axis on the Booster Neutrino Beam (BNB) and 6$^\circ$ off-axis from the Neutrinos at the Main Injector (NuMI) beam. Neutrinos from both beams inherit the timing sub-structure of their parent proton spills, which is in turn derived from either the Booster's or the Main Injector's synchrotron acceleration. Since neutrino propagation introduces only a constant offset, their timing structure is preserved as they travel. Identifying this structure in data represents a powerful tool for selecting neutrino events and searching for physics beyond the Standard Model (BSM). This poster presents the preliminary reconstruction of the BNB and NuMI neutrino bunch structure with ICARUS data, exploiting only the precise timing of ICARUS optical readout system to both locate and assign a time to each interaction.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Rapid curing dynamics of PEG-thiol-ene resins allow facile 3D bioprinting and in-air cell-laden microgel fabrication

Thiol-norbornene photoclick hydrogels are highly efficient in tissue engineering applications due to their fast gelation, cytocompatibility, and tunability. In this work, we utilized the advantageous features of polyethylene glycol (PEG)-thiol-ene resins to enable fabrication of complex and heterogeneous tissue scaffolds using 3D bioprinting and in-air drop encapsulation techniques. We demonstrated that photoclickable PEG-thiol-ene resins could be tuned by varying the ratio of PEG-dithiol to PEG norbornene to generate a wide range of mechanical stiffness (0.5–12 kPa) and swelling ratios. Importantly, all formulations maintained a constant, rapid gelation time (<0.5 s). We used this resin in biological projection microstereolithography (BioPµSL) to print complex structures with geometric fidelity and demonstrated biocompatibility by printing cell-laden microgrids. Moreover, the rapid gelling kinetics of this resin permitted high-throughput fabrication of tunable, cell-laden microgels in air using a biological in-air drop encapsulation apparatus (BioIDEA). We demonstrated that these microgels could support cell viability and be assembled into a gradient structure. This PEG-thiol-ene resin, along with BioPµSL and BioIDEA technology, will allow rapid fabrication of complex and heterogeneous tissues that mimic native tissues with cellular and mechanical gradients. The engineered tissue scaffolds with a controlled microscale porosity could be utilized in applications including gradient tissue engineering, biosensing, and in vitro tissue models.

36 MATERIALS SCIENCE↗
DOE OSTI2024

ML Classifier Fusion for Three Data Streams with Quality Inversely Proportional to Time Resolution

We consider a monitoring scenario of phenomenon using three different streams of measurements whose quality is proportional to their constant inter-arrival times. Each measurement of a stream needs to be binary-classified to reflect the state of interest of the phenomenon. A set of classifiers is separately trained and fused for each stream at its time resolution using measurements collected under known states. We present a machine learning method to fuse the outputs of these fusers to provide a final classification at the finest time resolution. We show that this fused-fusers method provides decisions with likely superior classification probability compared to the best individual classifiers and fused-classifiers. We derive generalization equations that guarantee a superior classification probability of fused-fusers with a confidence probability specified by the classifiers’ generalization equations. We apply these results to study a practical problem of classifying Pu/Np target dissolution events at a radiochemical processing facility using gamma spectral measurements of effluent flows.

Rao, Nageswara↗
DOE OSTI2024

Slow electron-phonon relaxation controls the dynamics of the superconducting resistive transition

Here, we investigate the temporal and spatial scales of resistance fluctuations (𝑅 fluctuations) at the superconducting resistive transition accessed through voltage fluctuation measurements in thin epitaxial TiN films. This material is characterized by slow electron-phonon relaxation, which puts it far beyond the applicability range of the textbook scenario of superconducting fluctuations. The measured Lorentzian spectrum of the 𝑅 fluctuations identifies their correlation time, which is nearly constant across the transition region and has no relation to the conventional Ginzburg-Landau timescale. Instead, the correlation time coincides with the energy relaxation time determined by a combination of the electron-phonon relaxation and the relaxation via diffusion into reservoirs. Our data are quantitatively consistent with the model of spontaneous temperature fluctuations and highlight the lack of understanding of the resistive transition in materials with slow electron-phonon relaxation.

critical phenomena↗
DOE OSTI2024

On the Stability of Power Transmission Systems Under Persistent Inverter Attacks: A Bi-Linear Matrix Approach

We investigate the stability and robustness properties of a power transmission system under persistent deceiving attacks on inverter-interfaced energy resources. The attacks can corrupt the damping coefficients in the inverters' controllers and measurements of the frequency at the points of coupling. Leveraging tools from hybrid dynamical systems theory, we characterize a broad family of persistent (and not necessarily periodic) attacks acting on the inverters, under which the stability properties of the transmission system can be shown to not be compromised. To address potentially conservative conditions identified through conventional bounding techniques, sufficient conditions on the average activation time of the attacks are identified via Lyapunov theory, as well as the formulation and solution of a class of bilinear matrix inequalities (BMI). The results are obtained for constant and slowly time-varying loads via input-to-state stability (ISS) tools. Numerical simulations on the IEEE 39-bus test system are also presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Ligand field exciton annihilation in bulk CrCl 3

The layered van der Waals material CrCl 3 exhibits very strongly bound ligand field excitons that control optoelectronic applications and are connected with magnetic ordering by virtue of their d-orbital origin. Time-resolved photoluminescence of these exciton populations at room temperature shows that their relaxation is dominated by exciton–exciton annihilation and that the spontaneous decay lifetime is very long. Furthermore, these observations allow the rough quantification of the exciton annihilation rate constant and contextualization in light of a recent theory of universal scaling behavior of the annihilation process.

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

AI-Batt (Autonomous Identification of Battery Life Models) [SWR 21-36]

Autonomous Identification of Battery Life Models (AI-Batt) AI-Batt is a MATLAB code base for developing lifetime models for batteries from accelerated aging data. The code base provides many functions for processing, visualizing, and modeling battery aging data, making the data processing, exploration, and modeling workflow substantially faster. These tools are tailored for working with battery aging data sets, which usually consist of many separate time-series for each cell, with many test conditions and possible replicates at each condition, which makes it difficult to simply process or visualize the data set. Complex modeling tasks, such as cross-validation, sensitivity analysis, and uncertainty quantification have been implemented to enable thorough statistical investigation of model predictions. Additionally, several machine-learning algorithms are implemented to autonomously identify suitable models via symbolic regression. Data processing functions automatically cast data from the struct data type, which is commonly used to store experimental data, but is not an acceptable input for most algorithms, to the table data type, which can be easily used as input to any optimization algorithm. Also, the data can be separated into time-invariant and time-variant data tables, which is helpful for exploring the data set as well as developing separate models for time-variant and time-invariant aging mechanisms. For example, in aging tests with constant temperature, temperature is a time-invariant experimental condition. Visualization tools enable plotting of data, model fits, and model simulations possible with single-line function calls, empowering data exploration of complex data sets with both time-varying and time-invariant trends. Plots can be automatically generated for the whole data set, or separated by data group (groups of test replicates) or individual data series. Data points or data series can be automatically colored by the value of a variable with a variety of color maps, and model predictions can also be colored by the value of a fit statistic. Comparisons between data sets and the predictions/simulations of different models on the same data set can be easily plotted as well. Distributions of parameter values from bootstrap resampling can be plotted to visualize the reliability of parameter estimation, or determine any correlations between parameters. Modeling tools handle the complex task of creating and parsing symbolic equations for modeling battery lifetime. Equations are parsed to grab relevant data variables, parameter values, or specified sub-models for input into optimization, evaluation, or simulation functions. Models can be optimized locally (one set of parameters for each data series), bi-level (some parameters shared across the data set), or globally (single set of parameters for all data). Functions implementing symbolic regression algorithms help users to discover effective model equations, even in poorly sampled, high-dimensional data.

Smith, Kandler [National Renewable Energy Lab. (NR↗
DOE OSTI2024

Effect of 60 Co Irradiation on Boehmite Dissolution in Caustic Solutions

Here, in this work, we examine how radiation impacts the dissolution behavior of boehmite by subjecting dry nanoparticles of different sizes to 60 Co γ radiation and subsequently analyzing their dissolution behavior in caustic solutions as a function of temperature. The measured kinetics show that irradiation with an amount 228.24 Mrad significantly slows the dissolution rate, particularly for smaller sizes at lower temperatures. Specifically, the temperature-dependent dissolution rates of irradiated 20 nm boehmite versus pristine material in 3 M NaOH solutions were several times lower (e.g., rate constant of 0.026 vs 0.075 h –1 at 60 °C), with an apparent activation energy 40 kJ mol –1 higher. Although various imaging techniques and X-ray diffraction measurements consistently revealed no obvious differences between pristine and irradiated samples, after irradiation significant binding energy shifts were detected in the X-ray photoelectron Spectroscopy peaks of Al 2p and O 1s, and a change in their relative intensities indicated a lower O/Al ratio. This suggests that γ-irradiation may stabilize boehmite particle surfaces by driving their chemistry and structure toward more stable aluminum oxide forms. This finding may help explain slower dissolution rates of boehmite in nuclear waste and may be useful for the development of more robust predictive models and effective strategies for waste processing.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Control of Excitonic Energy Transfer in RGB Quantum Dot:Polymer Composites for Tunable White Emission

Tint-controlled white light is crucial for both illumination systems and display applications. Here, in this study, we demonstrate solution-processed quantum-dot light-emitting diodes (QD-LEDs) featuring a red–green–blue (RGB) QD–poly(methyl methacrylate) (PMMA) composite emissive layer (EML) for tunable white electroluminescence (EL). In this composite EML, PMMA functions as a dispersion matrix that modulates the interdot spacing (d), thereby controlling Förster resonance energy transfer (FRET) between QDs. By adjustment of the PMMA content, the balance of R, G, and B emissions is controlled, enabling systematic and continuous tuning of the EL color from greenish to reddish white at a fixed RGB ratio and constant driving bias. Time-resolved photoluminescence measurements confirm that the variation in the exciton lifetime with the PMMA fraction is the primary factor for tuning the EL color. Notably, nearly pure white EL with CIE coordinates close to (0.33, 0.33) is achieved using a diluted PMMA matrix without significant degradation of the electrical properties. Our results demonstrate d as an independent design parameter for decoupling color tuning from RGB composition and electrical operation, providing a versatile design framework for high-quality white- or tint-controlled QD-LEDs toward advanced solid-state lighting and display technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Pressure induced modification of electronic and magnetic properties of MnCrNbAl and MnCrTaAl

Spin-gapless semiconductor (SGS) is a new class of material that has been studied recently for potential applications in spintronics. This material behaves as an insulator for one spin channel, and as a gapless semiconductor for the opposite spin. In this work, we present results of a computational study of two quaternary Heusler alloys, MnCrNbAl and MnCrTaAl that have been recently reported to exhibit spin-gapless semiconducting electronic structure. In particular, using density functional calculations we analyze the effect of external pressure on electronic and magnetic properties of these compounds. It is shown that while these two alloys exhibit nearly SGS behavior at optimal lattice constants and at negative pressure (expansion), they are half-metals at equilibrium, and magnetic semiconductors at larger lattice constant. At the same time, reduction of the unit cell volume has a detrimental effect on electronic properties of these materials, by modifying the exchange splitting of their electronic structure and ultimately destroying their half-metallic/semiconducting behavior. Thus, our results indicate that both MnCrNbAl and MnCrTaAl may be attractive for practical device applications in spin-based electronics, but a potential compression of the unit cell volume (e.g. in thin-film applications) should be avoided.

Materials Science↗
DOE OSTI2024

TDCOSMO - XVI. Measurement of the Hubble constant from the lensed quasar WGD 2038–4008

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant (H0). The TDCOSMO Collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038−4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain DΔ t = 1.68−0.38+0.40 Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat ΛCDM cosmology, we constrain the Hubble constant to be H0 = 65−14+23 km s−1 Mpc−1. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.Key words: gravitational lensing: strong / cosmological parameters / distance scale⋆ Corresponding author; kcwong19@gmail.com.⋆⋆ NHFP Einstein fellow.

79 ASTRONOMY AND ASTROPHYSICS↗