Search NASA⌕ Search

SEARCH · Search NASA

Results for “spectral fitting”

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

Multimodal, microspectroscopic speciation of legacy phosphorus in two US mid-Atlantic agricultural soils

To understand phosphorus (P) mobility in agricultural soils and its potential environmental risk, it is essential to directly measure solid phase P speciation. Often, bulk P K-edge X-ray absorption near edge structure (XANES) spectroscopy followed by linear combination fitting (LCF) is utilized to determine the solid P phases in soil. However, this method may limit results to only a few major phases. Additionally, XANES spectra for different P species may have very similar features, leading to an over- or underestimate of their contribution to LCF. Here, an improved P speciation by pairing multimodal microbeam-X-ray fluorescence (µ-XRF) mapping coupled with µ-XANES (microbeam-X-ray absorption near edge structure) analysis to directly speciate major and minor P phases on the micron scale is provided. We combined maps of both tender (P, sulfur, aluminum, and silicon) and hard energy (calcium, iron [Fe], and manganese) elements to evaluate the elemental co-locations with P. To better account for uncertainty assigning XANES peaks to individual compounds, a more quantitative fingerprinting by “spectral feature analysis” was completed. With this analysis, an R-factor is reported for the fit. These results were compared to traditional LCF. Pre-edge fitting results revealed the presence of a two-component pre-edge feature for phosphate adsorbed to ferrihydrite. Additionally, phytate co-precipitated with ferrihydrite (Phytate-Fe-Cop) had a pre-edge feature, indicating direct association with Fe. Lastly, a unique P species associated with manganese oxide was identified in the soil via multimodal mapping and µ-XANES. These results allow for better prediction of P dissolution and mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Robust Automatic EXAFS First-Shell Fits

Extended X-ray absorption fine structure (EXAFS) is a widely used technique for atomic structure determination. Fourier transformation connects EXAFS in k space and R space. However, determining the appropriate k-range for the transformation can be challenging, but critical for the first-shell fit. In this study, we present an automatic method to determine the k-range using the Larch package and a Python program. The first step is to estimate spectral noise across a series of k-ranges with a fixed minimum value and identify the optimal maximum value in the k-range (k max ). The k max is determined by an empirical noise threshold that marks the point where the noise level in the Fourier transformed spectrum changes dramatically. Using the obtained k max value, the first shell is modeled to determine the minimum k value (k min ) by optimizing the background function through alignment of the spectrum with theory. The optimal k min corresponds to the point of the minimum R-factor, which quantifies the difference between the experimental and fitted spectrum. Our method was tested on various typical datasets and yielded suitable k-ranges for Fourier transformation and accurate first-shell fits. This approach helps avoid unreliable, irreproducible data analysis, especially for noisy data from diluted samples, and enables robust automatic first-shell EXAFS fitting.

EXAFS analysis↗
DOE OSTI2024

The putative center in NGC 1052

Context. Many active galaxies harbor powerful relativistic jets, however, the detailed mechanisms of their formation and acceleration remain poorly understood.Aims. To investigate the area of jet acceleration and collimation with the highest available angular resolution, we study the innermost region of the bipolar jet in the nearby low-ionization nuclear emission-line region (LINER) galaxy NGC 1052.Methods. We combined observations of NGC 1052 taken with VLBA, GMVA, and EHT over one week in the spring of 2017. Our study is focused on the size and continuum spectrum of the innermost region containing the central engine and the footpoints of both jets. We employed a synchrotron-self absorption model to fit the continuum radio spectrum and we combined the size measurements from close to the central engine out to ∼1 pc to study the jet collimation.Results. For the first time, NGC 1052 was detected with the EHT, providing a size of the central region in-between both jet bases of 43 μas perpendicular to the jet axes, corresponding to just around 250 RS (Schwarzschild radii). This size estimate supports previous studies of the jets expansion profile which suggest two breaks of the profile at around 3 × 103 RS and 1 × 104 RS distances to the core. Furthermore, we estimated the magnetic field to be 1.25 Gauss at a distance of 22 μas from the central engine by fitting a synchrotron-self absorption spectrum to the innermost emission feature, which shows a spectral turn-over at ∼130 GHz. Assuming a purely poloidal magnetic field, this implies an upper limit on the magnetic field strength at the event horizon of 2.6 × 104 Gauss, which is consistent with previous measurements.Conclusions. The complex, low-brightness, double-sided jet structure in NGC 1052 makes it a challenge to detect the source at millimeter (mm) wavelengths. However, our first EHT observations have demonstrated that detection is possible up to at least 230 GHz. This study offers a glimpse through the dense surrounding torus and into the innermost central region, where the jets are formed. This has enabled us to finally resolve this region and provide improved constraints on its expansion and magnetic field strength.Key words: methods: observational / techniques: high angular resolution / techniques: interferometric / galaxies: active / galaxies: jets / galaxies: Seyfert⋆⋆ NASA Hubble Fellowship Program, Einstein Fellow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

XRF-MAPS

x-ray fluorescence (XRF) imaging typically involves the creation and analysis of 3D data sets, where at each scan position a full energy dispersive x-ray spectrum is recorded. This allows one to later process the data in a variety of different approaches, e.g., by spectral region of interest (ROI) summation with or without background subtraction, principal component analysis, or fitting. XRF-Maps is a C++ open source software package that implements these functions to provide a tool set for the analysis of XRF data sets.

GLOWACKI, ARTHURT↗
DOE OSTI2026

Improved Spectral Classification of Local K/M Dwarfs in SDSS Surveys. I. Chemodynamic Validation and Metallicity Calibration

An examination of 109,276 spectra of low-mass stars in the Sloan Digital Sky Survey (SDSS) data archive, collected pre-2010, provides a broad collection of K/M (sub)dwarfs tracing the local (d ≲ 250 pc) population of the thin-disk, thick-disk, and halo, based on 3D kinematics. These populations have distinct metallicities and kinematics, which should be reflected in lower-mass members. One complication is measuring metallicities of M dwarfs from low-resolution spectroscopy or photometry alone remains a challenge, even with the availability of Gaia data. To better characterize physical parameters of low-mass stars observed by SDSS, we define a set of 536 improved, empirical K/M dwarf classification templates to more accurately determine spectral subtypes (K5.5–M8.5) and morphology classes (MCs) (0.5–12.5). We select the best-fit template to every archival SDSS spectrum in the range 5000–8000 Å based on minimum χ 2 value. We then use secondary metallicity estimators to calibrate the most likely [Fe/H] values corresponding to the assigned MC. We confirm that the most metal-rich K/M dwarfs have thin-disk kinematics, and we observe the known effect of Galactic radial metallicity migration, which we confirm is strongly imprinted in the chemodynamics of the local disk M dwarfs. Confirmation of these chemodynamic trends validates the precision of our improved templates for estimating [Fe/H] for M (sub)dwarfs, notably at the low-metallicity end (–3 < [Fe/H] < –1). Substructure in the chemodynamics of the most metal-poor K/M subdwarfs begins to distinguish between local Gaia-Enceladus and local in situ halo objects. This methodology opens the door to use the ubiquitous, local K/M (sub)dwarfs to trace local Galactic chemodynamics with the highest possible resolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Spectro-Microscopic Analysis of Soot Particle Composition and Source Attribution

Ambient soot particles significantly impact Earth’s radiative balance, human health, and atmospheric visibility. Their microstructural properties depend on formation and aging mechanisms, which vary by emission source and atmospheric processes. Hence, accurately identifying sources of soot enhances our understanding of their physicochemical properties and atmospheric implications. This study used a multi-modal approach to characterize and attribute sources of submicron soot particles collected in Israel during new particle formation events, biomass burning episodes, and background atmospheric conditions. Synchrotron-based X-ray microscopy was used to map soot (elemental carbon), organic carbon, and inorganic species. Implemented atomic force microscopy showed highly diverse phase states, with soot consistently exhibiting a solid-like phase. Automated µ-Raman analysis was subsequently performed on ~690 particles, identifying three soot classes based on spectral features corresponding to the "Defect" (D) and "Graphite" (G) bands of soot. We applied two-peak and five-peak fitting approaches to deconvolute the “Defect” peaks (D1, D2, D3, and D4) and G band from average Raman, revealing varying degrees of graphitic order. The degree of graphitic order was determined from metrics such as the D3 peak area, often observed when soot was internally mixed with organic material. Raman spectral features, along with temporal variations in particle classes contributions, suggest that Particle Type 1 corresponds to traffic related soot and Particle Type 2 to less graphitic soot from biomass burning, while Particle Type 3 is associated with more heterogeneous particulate representative of soot-OC mixtures emitted during new particle formation and biomass burning episodes.

Rivera-Adorno, Felipe (ORCID:0000000273557999)↗
DOE OSTI2022

Revisiting stellar properties of star-forming galaxies with stellar and nebular spectral modelling

Spectral synthesis is a powerful tool for interpreting the physical properties of galaxies by decomposing their spectral energy distributions (SEDs) into the main luminosity contributors (e.g. stellar populations of distinct age and metallicity or ionised gas). However, the impact nebular emission has on the inferred properties of star-forming (SF) galaxies has been largely overlooked over the years, with unknown ramifications to the current understanding of galaxy evolution. The objective of this work is to estimate the relations between stellar properties (e.g. total mass, mean age, and mean metallicity) of SF galaxies by simultaneously fitting the stellar and nebular continua and comparing them to the results derived through the more common purely stellar spectral synthesis approach. The main galaxy sample from SDSS DR7 was analysed with two distinct population synthesis codes: FADO, which estimates self-consistently both the stellar and nebular contributions to the SED, and the original version of STARLIGHT, as representative of purely stellar population synthesis codes. Differences between codes regarding average mass, mean age and mean metallicity values can go as high as ~0.06 dex for the overall population of galaxies and ~0.12 dex for SF galaxies (galaxies with EW(Hα) > 3 Å), with the most prominent difference between both codes in the two populations being in the light-weighted mean stellar age. FADO presents a broader range of mean stellar ages and metallicities for SF galaxies than STARLIGHT, with the latter code preferring metallicity solutions around the solar value (Z ⊙ = 0.02). A closer look into the average light- and mass-weighted star formation histories of intensively SF galaxies (EW(Hα) > 75 Å) reveals that the light contributions of simple stellar populations (SSPs) younger than ≤10 7 (10 9 ) years in STARLIGHT are higher by ~5.41% (9.11%) compared to FADO. Moreover, FADO presents higher light contributions from SSPs with metallicity ≤ Z ⊙ /200 (Z ⊙ /50) of around 8.05% (13.51%) when compared with STARLIGHT. This suggests that STARLIGHT is underestimating the average light-weighted age of intensively SF galaxies by up to ~0.17 dex and overestimating the light-weighted metallicity by up to ~0.13 dex compared to FADO (or vice versa). The comparison between the average stellar properties of passive, SF and intensively SF galaxy samples also reveals that differences between codes increase with increasing EW(Hα) and decreasing total stellar mass. Moreover, comparing SF results from FADO in a purely stellar mode with the previous results qualitatively suggests that differences between codes are primarily due to mathematical and statistical differences and secondarily due to the impact of the nebular continuum modelling approach (or lack thereof). However, it is challenging to adequately quantify the relative role of each factor since they are likely interconnected. This work finds indirect evidence that a purely stellar population synthesis approach negatively impacts the inferred stellar properties (e.g. mean age and mean metallicity) of galaxies with relatively high star formation rates (e.g. dwarf spirals, ‘green peas’, and starburst galaxies). In turn, this can bias interpretations of fundamental relations such as the mass-age or mass-metallicity, which are factors worth bearing in mind in light of future high-resolution spectroscopic surveys at higher redshifts (e.g. MOONS and 4MOST-4HS).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

time-resolved spectroscopy fit (trspecfit) v0.01

Analyze 2D time- and energy-resolved data, such as from a pump-probe spectroscopy experiment. User can select and input different peak shapes/ functions and background types to first fit a ground state/ unperturbed spectrum. This would be similar to how standard spectroscopy data is fit. Subsequently, to describe the time domain, users can choose functions that describe the temporal dynamics of one or more spectral features, such as a peak amplitude, peak position, etc. These time dynamics functions can be added or convoluted (e.g. describing an instrument response function) with each other. Functionality to integrate implicit variables leading to distributions of certain parameters/ functions is in development. Alternatively, 2D data can be analyzed one time step at a time to get an idea of the time dynamics of the system before deploying the global 2D fit described above. Typically people write custom software for this purpose. During my PhD I've seen five internal LBL and external researchers write one-off code in different languages to analyze time- and energy-resolved spectra. While this was specifically was for a laser pump - X-ray probe spectroscopy experiment, I'm trying to write a general package for the time-resolved spectroscopy community.

Mahl, Johannes [Lawrence Berkeley National Laborat↗
DOE OSTI2025

Initial measurement of reactor antineutrino oscillation at SNO+

The SNO+ collaboration reports its first spectral analysis of long-baseline reactor antineutrino oscillation using 114 tonne-years of data. Fitting the neutrino oscillation probability to the observed energy spectrum yields constraints on the neutrino mass-squared difference Δ$m^2_{21}$. In the ranges allowed by previous measurements, the best-fit Δ$m^2_{21}$ is ($8.85^{+1.10}_{-1.33}$) × 10 -5 eV 2 . This measurement is continuing in the next phases of SNO+ and is expected to surpass the present global precision on Δ$m^2_{21}$ with about three years of data. In the following, we first describe the configuration of the SNO+ detector when it was partially filled with scintillator, and a characterization of the detector response using intrinsic radioactivity. Next, we detail the event selection and expectations for reactor IBDs and (α, n) reactions. Then, we present the results of an energy spectrum analysis using 114 tonne-years of data. We conclude with prospects of future results from the SNO+ detector, which has been operating fully-filled with 780 tonnes of scintillator

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

DIGS: deep inference of galaxy spectra with neural posterior estimation

Abstract With the advent of billion-galaxy surveys with complex data, the need of the hour is to efficiently model galaxy spectral energy distributions (SEDs) with robust uncertainty quantification. The combination of simulation-based inference (SBI) and amortized neural posterior estimation (NPE) has been successfully used to analyse simulated and real galaxy photometry both precisely and efficiently. In this work, we utilise this combination and build on existing literature to analyse simulated noisy galaxy spectra. Here, we demonstrate a proof-of-concept study of spectra that is (a) an efficient analysis of galaxy SEDs and inference of galaxy parameters with physically interpretable uncertainties; and (b) amortized calculations of posterior distributions of said galaxy parameters at the modest cost of a few galaxy fits with Markov chain Monte Carlo (MCMC) methods. We utilise the SED generator and inference framework Prospector to generate simulated spectra, and train a dataset of 2 × 10 6 spectra (corresponding to a five-parameter SED model) with NPE. We show that SBI—with its combination of fast and amortized posterior estimations—is capable of inferring accurate galaxy stellar masses and metallicities. Our uncertainty constraints are comparable to or moderately weaker than traditional inverse-modelling with Bayesian MCMC methods (e.g. 0.17 and 0.26 dex in stellar mass and metallicity for a given galaxy, respectively). We also find that our inference framework conducts rapid SED inference (0.9–1.2 × 10 5 galaxy spectra via SBI/NPE at the cost of 1 MCMC-based fit). With this work, we set the stage for further work that focuses of SED fitting of galaxy spectra with SBI, in the era of JWST galaxy survey programs and the wide-field Roman Space Telescope spectroscopic surveys.

spectroscopy↗
DOE OSTI2022

Direct cosmological inference from three-dimensional correlations of the Lyman α forest

ABSTRACT When performing cosmological inference, standard analyses of the Lyman α (Ly α) three-dimensional correlation functions only consider the information carried by the distinct peak produced by baryon acoustic oscillations (BAO). In this work, we address whether this compression is sufficient to capture all the relevant cosmological information carried by these functions. We do this by performing a direct fit to the full shape, including all physical scales without compression, of synthetic Ly α autocorrelation functions and cross-correlations with quasars at effective redshift zeff = 2.3, assuming a DESI-like survey, and providing a comparison to the classic method applied to the same data set. Our approach leads to a $3.5{{\ \rm per\ cent}}$ constraint on the matter density ΩM, which is about three to four times better than what BAO alone can probe. The growth term fσ8(zeff) is constrained to the $10{{\ \rm per\ cent}}$ level, and the spectral index ns to $\sim 3-4{{\ \rm per\ cent}}$. We demonstrate that the extra information resulting from our ‘direct fit’ approach, except for the ns constraint, can be traced back to the Alcock–Paczynski effect and redshift space distortion information.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The gMOSS: the galaxy survey and galaxy populations of the large homogeneous field

ABSTRACT We present the gMOSS (Galaxies of Medium-band One-meter Schmidt telescope Survey) catalogue of ∼19 000 galaxies in 20 filters (4 broad-band SDSS and 16 medium-band filters). We observed 2.386 deg2 on the central part of the HS47.5-22 field with the 1-m Schmidt telescope of the Byurakan Astrophysical Observatory. The gMOSS is a complete flux-limited sample of galaxies with a threshold magnitude of r SDSS ≤ 22.5 AB. From photometric measurements with 16 medium-band filters and u SDSS, we get spectral energy distributions for each object in the field, which are used for further analysis. Galaxy classification and photometric redshift estimation based on spectral template matching with zebra software. The obtained redshift accuracy is σNMAD < 0.0043. Using the SED-fitting cigale code, we obtained the main properties of the stellar population of galaxies, such as rest-frame (u − r)res colour, stellar mass, extinction, and mass-weighted age with a precision of 0.16 ± 0.07 mag, 0.14 ± 0.04 dex, 0.27 ± 0.1 mag, and 0.08 ± 0.04 dex, respectively. Using a dust-corrected colour–mass diagram, we divided the full sample into populations of red and blue galaxies and considered the dependencies between stellar mass and age. Throughout cosmic time, red sequence galaxies remain older and more massive than blue cloud galaxies. The star formation history of a complete subsample of galaxies selected in the redshift range 0.05 ≤ z ≤ 0.015 with <$\mathrm{log} M \mathrm{\gt }_\mathrm{[M_\odot ]}$>8.3 shows an increase in the SFRD up to z ∼ 3, under the results obtained in earlier studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Exploring Changes in Quasar Spectral Energy Distributions across C IV Parameter Space

We examine the UV/X-ray properties of 1378 quasars in order to link empirical correlations to theoretical models of the physical mechanisms dominating quasars as a function of mass and accretion rate. The clarity of these correlations is improved when (1) using C IV broad emission line equivalent width (EQW) and blueshift (relative to systemic) values calculated from high signal-to-noise ratio reconstructions of optical/UV spectra and (2) removing quasars expected to be absorbed based on their UV/X-ray spectral slopes. In addition to using the traditional C IV parameter space measures of C IV EQW and blueshift, we define a "C IV ∥ distance" along a best-fit polynomial curve that incorporates information from both C IV parameters. We find that the C IV ∥ distance is linearly correlated with both the optical-to-X-ray slope, α ox , and broad-line He II EQW, which are known spectral energy distribution indicators, but does not require X-ray or high spectral resolution UV observations to compute. The C IV ∥ distance may be a better indicator of the mass-weighted accretion rate, parameterized by L/L Edd , than the C iv EQW or blueshift alone, as those relationships are known to break down at the extrema. Conversely, there is only a weak correlation with the X-ray energy index (Γ), an alternate L/L Edd indicator. We find no X-ray or optical trends in the direction perpendicular to the C IV distance that could be used to reveal differences in accretion disk, wind, or corona structure that could be widening the C IV EQW–blueshift distribution. A different parameter (such as metallicity) not traced by these data must come into play.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Elastic and resonance structures of the nucleon from the hadronic tensor in lattice QCD: Implications for neutrino-nucleon scattering and hadron physics

We compute the Euclidean hadronic tensor from charge density operators and extract elastic and resonance structures by employing exponential fits to the four-point function correlator, as well as a Bayesian reconstruction inverse algorithm to obtain the corresponding spectral density for qualitative comparison. We present the determination of the nucleon’s Sachs electric form factor using the hadronic tensor formalism and verify that it is consistent with that from the conventional three-point function calculation. Beyond the elastic peak, we observe a structure located approximately 0.5–0.7 GeV above the nucleon mass in the Bayesian reconstruction. This structure is interpreted as a mixture of the Roper resonance [𝑁⁡(1440)], and states with both positive and negative parities in this mass region, as well as multihadron states. Assuming the observed structure is dominated by 𝐽 𝑃 = 1/2 ± states, we extract the transition electric form factor 𝐺$^*_𝐸$⁡(𝑄 2 ) and the corresponding longitudinal helicity amplitude 𝑆 1/2 ⁡(𝑄 2 ), and compare them with those determined from the CLAS experimental data of nucleon-to-Roper transition. Although fitting to the four-point correlation function or using the inverse algorithm does not resolve individual resonances, it nevertheless enables the determination of total inclusive lepton–nucleon scattering cross sections in appropriate energy bins. This lattice QCD calculation presents the first major step toward studying the inclusive 𝑁 → 𝑋 contributions within the hadronic tensor formalism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Optical polarization and spectral properties of the hydrogen-poor superluminous supernovae SN 2021bnw and SN 2021fpl

ABSTRACT New optical photometric, spectroscopic, and imaging polarimetry data are combined with publicly available data to study some of the physical properties of the two hydrogen-poor superluminous supernovae (SLSNe) SN 2021bnw and SN 2021fpl. For each SLSN, the best-fitting parameters obtained from the magnetar model with Modular Open-Source Fitter for Transients do not depart from the range of parameter obtained on other SLSNe discussed in the literature. A spectral analysis with SYN++ shows that SN 2021bnw is a W type, fast evolver, while SN 2021fpl is a 15bn type, slow evolver. The analysis of the polarimetry data obtained on SN 2021fpl at four epochs (+1.8, +20.6, +34.1, and +43.0 d, rest frame) shows >3σ polarization detections in the range of 0.8–1 per cent. A comparison of the spectroscopy data suggests that SN 2021fpl underwent a spectral transition a bit earlier than SN 2015bn, during which, similarly, it could have underwent a polarization transition. The analysis of the polarimetry data obtained on SN 2021bnw does not show any departure from symmetry of the photosphere at an empirical diffusion time-scale of ≈2 (+81.1 d rest frame). This result is consistent with those on the sample of W-type SLSN observed at empirical diffusion time-scale ≤ 1 with that technique, even though it is not clear the effect of limited spectral windows varying from one object to the other. Measurements at higher empirical diffusion time-scale may be needed to see any departure from symmetry as it is discussed in the literature for SN 2017egm.

Astronomy & Astrophysics↗
DOE OSTI2026

Improved Accuracy in Semi-Experimental Structure Determination by Resolving Problems Associated with Rotation of Principal Inertial Axes of Isotopologues: Structures of 1,3-Oxazole ( c -C 3 H 3 NO)

The rotational spectrum of the normal isotopologue of 1,3-oxazole (c-C 3 H 3 NO) was observed from 43 to 750 GHz. Over 3900 transitions for the ground vibrational state are measured, assigned, and least-squares fit to sextic centrifugally distorted-rotor Hamiltonians. The measured frequencies and resulting spectroscopic constants from this extended spectral range, combined with previous measurements of the nuclear quadrupole coupling constants, will facilitate astronomical searches for oxazole across the majority of the range of modern radiotelescopes. Spectra for a set of 30 oxazole isotopologues, which include multiple isotopic substitutions of each atom, are used to determine the first semi-experimental equilibrium ($r$$^{SE}_{e}$) structure and semi-experimental substitution structure ($r$$^{SE}_{e}$), each using CCSD(T) computed values for the vibration–rotation interaction and electron-mass corrections. The large number of isotopologues, including 21 isotopologues observed for the first time, and the redundant substitutions of each atom provide sufficient spectroscopic information to determine the $r$$^{SE}_{e}$ structure with the expected high level of accuracy and precision (0.0001 or 0.0002 Å in bond distances and 0.013 to 0.025° in bond angles). In the course of this study, we analyzed a known issue for some $r$$^{SE}_{e}$ structure determinations of near-oblate asymmetric tops in which inclusion of individual isotopologues degrades the structure determination. We demonstrate that this problem primarily arises from the difference in the values of the computed vibration–rotation interaction corrections as evaluated at the computed re geometry vs the $r$$^{SE}_{e}$ geometry of the “real” molecule. Our solution to this problem substantially improves the $r$$^{SE}_{e}$ structure of oxazole and likely can be generalized to many other molecules.

Chemical structure↗
DOE OSTI2022

Aging of galaxies along the morphological sequence, marked by bulge growth and disk quenching

Aims: We revisit the color bimodality of galaxies using the extensive EFIGI morphological classification of nearby galaxies. Methods: The galaxy profiles from the Sloan Digital Sky Survey (SDSS) gri images were decomposed as a bulge and a disk by controlled profile modeling with the Euclid SourceXtractor++ software. The spectral energy distributions from our resulting gri SDSS photometry complemented with Galaxy Evolution Explorer (GALEX) NUV photometry were fitted with the ZPEG software and PEGASE.2 templates in order to estimate the stellar masses and specific star formation rates (sSFR) of whole galaxies as well as their bulge and disk components. Results: The absolute NUV-r color versus stellar mass diagram shows a continuous relationship between the present sSFR of galaxies and their stellar mass, which spans all morphological types of the Hubble sequence monotonously. Irregular galaxies to intermediate-type Sab spirals make up the “Blue Cloud” across 4 orders of magnitude in stellar mass but a narrow range of sSFR. This mass build-up of spiral galaxies requires major mergers, in agreement with their frequently perturbed isophotes. At high mass, the Blue Cloud leads to the “Green Plain”, dominated by S0a and Sa early-type spirals. It was formerly called the “Green Valley”, due to its low density, but we rename it because of its wide stretch and nearly flat density over ~2 mag in NUV-r color (hence sSFR), despite a limited range of stellar mass (1 order of magnitude). The Green Plain links up the “Red Sequence”, containing all lenticular and elliptical galaxies with a 2 order of magnitude mass interval, and systematically higher masses for the ellipticals. We confirm that the Green Plain cannot be studied using u - r optical colors because it is overlayed by the Red Sequence, hence NUV data are necessary. Galaxies across the Green Plain undergo a marked growth by a factor 2 to 3 in their bulge-to-total mass ratio and a systematic profile change from pseudo to classical bulges, as well as a significant reddening due to star formation fading in their disks. The Green Plain is also characterized by a maximum stellar mass of 10 11.7 M ⊙ beyond which only elliptical galaxies exist, hence supporting the scenario of ellipticals partly forming by major mergers of massive disk galaxies. Conclusions: The EFIGI attributes indicate that dynamical processes (spiral arms and isophote distortions) contribute to the scatter of the Main Sequence of star-forming galaxies (Blue Cloud), via the enhancement of star formation (flocculence, HII regions). The significant bulge growth across the Green Plain confirms that it is a transition region, and excludes a predominantly quick transit due to rapid quenching. The high frequency of bars for all spirals as well as the stronger spiral arms and flocculence in the knee of the Green Plain suggest that internal dynamics, likely triggered by flybys or (mainly minor) mergers, may be the key to the bulge growth of massive disk galaxies, which is a marker of the aging of galaxies from star forming to quiescence. The Hubble sequence can then be considered as an inverse sequence of galaxy physical evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Construction of the damped Ly⁢𝛼 absorber catalog for DESI DR2 Ly⁢𝛼 BAO

We present the Damped Ly⁢𝛼 Toolkit for automated detection and characterization of damped Ly⁢𝛼 absorbers (DLAs) in quasar spectra. Our method uses quasar spectral templates with and without absorption from intervening DLAs to reconstruct observed quasar forest regions. The best-fitting model determines whether a DLA is present while estimating the redshift and HI column density. With an optimized quality cut on detection significance (Δ⁢𝜒$^{2}_{𝑟}$ >0.03), the technique achieves an estimated 80% purity and 79% completeness when evaluated on simulated spectra with S/N>2 that are free of broad absorption lines (BALs). We provide a catalog containing candidate DLAs from the DLA Toolkit detected in DESI DR1 quasar spectra, of which 21 719 were found in S/N>2 spectra with predicted log 10 ⁡(𝑁 𝙷𝙸 )>20.3 and detection significance Δ⁢𝜒$^{2}_{𝑟}$ >0.03. We compare the Damped Ly⁢𝛼 Toolkit to two alternative DLA finders based on a convolutional neural network and Gaussian process models. We present a strategy for combining these three techniques to produce a high-fidelity DLA catalog from DESI DR2 for the Ly⁢𝛼 forest baryon acoustic oscillation measurement. The combined catalog contains 41 152 candidate DLAs with log 10 ⁡(𝑁 𝙷𝙸 )>20.3 from quasar spectra with S/N>2. We estimate this sample to be approximately 85% pure and 79% complete when BAL quasars are excluded.

79 ASTRONOMY AND ASTROPHYSICS↗