adsabs/SciX_UAT_keywords huggingface dataset:
Kernel Selection for Gaussian Process in Cosmology: With Approximate Bayesian Computation Rejection and Nested Sampling

The Gaussian process (GP) has gained much attention in cosmology due to its ability to reconstruct cosmological data in a model-independent manner. In this study, we compare two methods for GP kernel selection: approximate Bayesian computation (ABC) rejection and nested sampling. We analyze three types of data: cosmic chronometer data, type Ia supernovae data, and gamma-ray burst data, using five kernel functions. To evaluate the differences between kernel functions, we assess the strength of evidence using Bayes factors. Our results show that, for ABC rejection, the Matérn kernel with ν = 5/2 (M52 kernel) outperformes the commonly used radial basis function (RBF) kernel in approximating all three data sets. Bayes factors indicate that the M52 kernel typically supports the observed data better than the RBF kernel but with no clear advantage over other alternatives. However, nested sampling gives different results, with the M52 kernel losing its advantage. Nevertheless, Bayes factors indicate no significant dependence of the data on each kernel.

Keywords: gaussian processes regression; astronomy data analysis; bayes factor; observational cosmology; astrostatistics strategies
The Age Dependence of Mid-infrared Emission around Young Star Clusters

Using the star cluster catalogs from the Hubble Space Telescope program Legacy Extragalactic UV survey (LEGUS) and 8 μm images from the IRAC camera on the Spitzer Space Telescope for five galaxies within 5 Mpc, we investigate how the 8 μm dust luminosity correlates with the stellar age on the 30-50 pc scale of star-forming regions. We construct a sample of 97 regions centered at local peaks of 8 μm emission, each containing one or more young star cluster candidates from the LEGUS catalogs. We find a tight anticorrelation with a Pearson correlation coefficient of r = -0.84 ± 0.05 between the mass-normalized dust-only 8 μm luminosity and the age of stellar clusters younger than 1 Gyr; the 8 μm luminosity decreases with increasing age of the stellar population. Simple assumptions on a combination of stellar and dust emission models reproduce the observed trend. We also explore how the scatter of the observed trend depends on assumptions of stellar metallicity, polycyclic aromatic hydrocarbon (PAH) abundance, fraction of stellar light absorbed by dust, and instantaneous versus continuous star formation models. We find that variations in stellar metallicity have little effect on the scatter, while PAH abundance and the fraction of dust-absorbed light bracket the full range of the data. We also find that the trend is better explained by continuous star formation, rather than instantaneous burst models. We ascribe this result to the presence of multiple star clusters with different ages in many of the regions. Upper limits of the dust-only 8 μm emission as a function of age are provided.

Keywords: interstellar medium; star formation; young star clusters; infrared sources; polycyclic aromatic hydrocarbons; stellar ages
Binding Energy Evaluation Platform: A Database of Quantum Chemical Binding Energy Distributions for the Astrochemical Community

The quality of astrochemical models is highly dependent on reliable binding energy (BE) values that consider the morphological and energetic variety of binding sites on the surface of ice-grain mantles. Here, we present the Binding Energy Evaluation Platform (BEEP) and database that, using quantum chemical methods, produces full BE distributions of molecules bound to an amorphous solid water (ASW) surface model. BEEP is highly automatized and allows one to sample binding sites on a set of water clusters and to compute accurate BEs. Using our protocol, we computed 21 BE distributions of interstellar molecules and radicals on an amorphized set of 15-18 water clusters of 22 molecules each. The distributions contain between 225 and 250 unique binding sites. We apply a Gaussian fit and report the mean and standard deviation for each distribution. We compare with existing experimental results and find that the low- and high-coverage experimental BEs coincide well with the high-BE tail and mean value of our distributions, respectively. Previously reported single BE theoretical values are broadly in line with ours, even though in some cases significant differences can be appreciated. We show how the use of different BE values impacts a typical problem in astrophysics, such as the computation of snow lines in protoplanetary disks. BEEP will be publicly released so that the database can be expanded to other molecules or ice models in a community effort.

Keywords: astrochemistry; interstellar dust processes; ice formation; surface ices; protoplanetary disks
A Sample Bias in Quasar Variability Studies

When a flux-limited quasar sample is observed at later times, there will be more dimmed quasars than brightened ones, due to a selection bias induced at the time of sample selection. Quasars are continuously varying and there are more fainter quasars than brighter ones. At the time of selection, even symmetrical variability will result in more quasars with their instantaneous fluxes scattered above the flux limit than those scattered below, leading to an asymmetry in flux changes over time. The same bias would lead to an asymmetry in the ensemble structure function (SF) of the sample such that the SF based on pairs with increasing fluxes will be slightly smaller than that based on pairs with decreasing fluxes. We use simulated time-symmetric quasar light curves based on the damped random walk prescription to illustrate the effects of this bias. The level of this bias depends on the sample, the threshold of magnitude changes, and the coverage of light curves, but the general behaviors are consistent. In particular, the simulations matched to recent observational studies with decade-long light curves produce an asymmetry in the SF measurements at the few percent level, similar to the observed values. These results provide a cautionary note on the reported time asymmetry in some recent quasar variability studies.

Keywords: quasars; surveys; black hole physics
Kepler Data Analysis: Non-Gaussian Noise and Fourier Gaussian Process Analysis of Stellar Variability

We develop a statistical analysis model of Kepler stellar flux data in the presence of planet transits, non-Gaussian noise, and stellar variability. We first develop a model for the Kepler noise probability distribution in the presence of outliers, which make the noise probability distribution non-Gaussian. We develop a signal likelihood analysis based on this probability distribution, in which we model the signal as a sum of the star variability and planetary transits. We argue that these components need to be modeled together if optimal signal is to be extracted from the data. For the stellar variability model we develop an optimal Gaussian process analysis using a Fourier-based Wiener filter approach, where the power spectrum is non-parametric and learned from the data. We develop high dimensional optimization of the objective function, where we jointly optimize all the model parameters, including thousands of star variability modes, and planet transit parameters. We apply the method to Kepler-90 data and show that it gives a better match to the stellar variability than the existing methods, and robustly handles noise outliers. As a consequence, the planet radii have a higher value than what the existing methods give, including splines and celerite.

Keywords: exoplanet detection methods; non-gaussianity; astronomy data modeling; exoplanet astronomy; astronomy data analysis
Accurate Modeling of Lyα Profiles and Their Impact on Photolysis of Terrestrial Planet Atmospheres

Accurately measuring and modeling the Lyα (Lyα; λ1215.67 Å) emission line from low-mass stars is vital for our ability to build predictive high energy stellar spectra, yet interstellar medium (ISM) absorption of this line typically prevents model-measurement comparisons. Lyα also controls the photodissociation of important molecules, like water and methane, in exoplanet atmospheres such that any photochemical models assessing potential biosignatures or atmospheric abundances require accurate Lyα host star flux estimates. Recent observations of three early M and K stars (K3, M0, M1) with exceptionally high radial velocities (>100 km s-1) reveal the intrinsic profiles of these types of stars as most of their Lyα flux is shifted away from the geocoronal line core and contamination from the ISM. These observations indicate that previous stellar spectra computed with the PHOENIX atmosphere code have underpredicted the core of Lyα in these types of stars. With these observations, we have been able to better understand the microphysics in the upper atmosphere and improve the predictive capabilities of the PHOENIX atmosphere code. Since these wavelengths drive the photolysis of key molecular species, we also present results analyzing the impact of the resulting changes to the synthetic stellar spectra on observable chemistry in terrestrial planet atmospheres.

Keywords: stellar chromospheres; stellar activity; ultraviolet astronomy; low mass stars
Silicate Extinction Profile Based on the Stellar Spectrum by Spitzer/IRS

The 9.7 and 18 μm interstellar spectral features, arising from the Si–O stretching and O–Si–O bending modes of amorphous silicate dust, are the strongest extinction features in the infrared. Here we use the "pair method" to determine the silicate extinction profile by comparing the Spitzer/IRS spectra of 49 target stars with obvious extinction with those of unreddened stars of the same spectral type. The 9.7 μm extinction profile is determined from all the 49 stars and the 18 μm profile is determined from six stars. It is found that the profile has a peak wavelength at ∼9.2–9.8 and ∼18–22 μm, respectively. The peak wavelength of the 9.7 μm feature seems to become shorter for stars of late spectral types; meanwhile the FWHM seems irrelevant to the spectral type, which may be related to circumstellar silicate emission. The silicate optical depth at 9.7 μm, Δτ 9.7, mostly increases with the color excess in J ‑ K S ( ${E}_{{{JK}}_{S}}$ ). The mean ratio of the visual extinction to the 9.7 μm silicate absorption optical depth is A V /Δτ 9.7 ≈ 17.8, in close agreement with that of the solar-neighborhood diffuse interstellar medium. When ${E}_{{{JK}}_{S}}\gt 4$ , this proportionality changes. The correlation coefficient between the peak wavelength and the FWHM of the 9.7 μm feature is 0.4, which indicates a positive correlation considering the uncertainties of the parameters. The method is compared with replacing the reference star by an atmospheric model spectral energy distribution and no significant difference is observed.

Keywords: interstellar dust extinction; interstellar dust; silicate grains; diffuse interstellar clouds
The Epoch of Giant Planet Migration Planet Search Program. II. A Young Hot Jupiter Candidate around the AB Dor Member HS Psc

We report the discovery of a hot Jupiter candidate orbiting HS Psc, a K7 (≈0.7 M ⊙) member of the ≈130 Myr AB Doradus moving group. Using radial velocities over 4 yr from the Habitable-zone Planet Finder spectrograph at the Hobby–Eberly Telescope, we find a periodic signal of ${P}_{b}={3.986}_{-0.003}^{+0.044}$ days. A joint Keplerian and Gaussian process stellar activity model fit to the radial velocities yields a minimum mass of ${m}_{p}\sin i={1.5}_{-0.4}^{+0.6}$ M Jup. The stellar rotation period is well constrained by the Transiting Exoplanet Survey Satellite light curve (P rot = 1.086 ± 0.003 days) and is not an integer harmonic nor alias of the orbital period, supporting the planetary nature of the observed periodicity. HS Psc b joins a small population of young, close-in giant planet candidates with robust age and mass constraints and demonstrates that giant planets can either migrate to their close-in orbital separations by 130 Myr or form in situ. Given its membership in a young moving group, HS Psc represents an excellent target for follow-up observations to characterize this young hot Jupiter further, refine its orbital properties, and search for additional planets in the system. *Based on observations obtained with the Hobby–Eberly Telescope (HET), which is a joint project of the University of Texas at Austin, the Pennsylvania State University, Ludwig-Maximillians-Universitaet Muenchen, and Georg-August Universitaet Goettingen. The HET is named in honor of its principal benefactors, William P. Hobby and Robert E. Eberly.

Keywords: exoplanet astronomy; exoplanet evolution; radial velocity
A Semi-empirical Approach to the Dynamic Coupling of CMEs and Solar Wind

Coronal mass ejections (CMEs) are one of the most relevant phenomena for space weather. Moreover, CMEs can negatively affect essential services and facilities. Therefore, to protect society, we require well-grounded knowledge of the physics that governs the propagation of CMEs from near the Sun to the orbit of Earth. In this work, we deduce expressions to approximate the main forces that affect the dynamic coupling between CMEs and the surrounding solar wind. Therefore, we explore the CME-solar wind dynamic coupling from a magnetohydrodynamic perspective, which, combined with a few reasonable assumptions, allows us to obtain expressions for the thermal and magnetic pressure forces, viscous and dynamic drag, and gravity. We simultaneously use our expressions to compute the trajectories of 34 Earth-directed CMEs. Our results, which are compared with in situ data, show significant quantitative consistency; our synthetic transits closely mimic their in situ observed counterparts. We conclude from our results that magnetic, thermal, and dynamic drag significantly surpass the other forces such as dynamic agents of CMEs in the interplanetary medium. In addition, we find that the initial relative speed of CMEs and solar wind is a determinant factor for the dynamic behavior of CMEs. In other words, subsonic CMEs are initially mostly affected by magnetic and thermal pressure forces, whereas super-magnetosonic CMEs are initially governed by inertial drag.

Keywords: solar coronal mass ejections; solar wind; magnetohydrodynamics; space weather; interplanetary medium
Unraveling the Complex Structure of AGN-driven Outflows. VI. Strong Ionized Outflows in Type 1 AGNs and the Outflow Size-Luminosity Relation

We present spatially resolved gas kinematics, ionization, and energetics of 11 type 1 and 5 type 2 active galactic nuclei (AGNs) with strong ionized gas outflows at z <0.3 using Gemini Multi-Object Spectrograph Integral Field Unit data. We find a strongly blueshifted region in [O III] velocity maps, representing an approaching cone in biconical outflows, and blueshifted and redshifted regions in Hα velocity maps, which show gravitationally rotating kinematics. AGN photoionization is dominant in the central region of most targets, and some of them also show ring-like structures of LINER or composite that surround the AGN-dominated center. Following our previous studies, we kinematically determine outflow sizes by the ratio between [O III] and stellar velocity dispersion. Outflow sizes of type 1 AGNs follow the same kinematic outflow size-[O III] luminosity relation obtained from the type 2 Integral Field Unit sample in Kang & Woo and Luo (updated slope 0.29 ± 0.04), while they are limited to the central kiloparsec scales, indicating the lack of global impact of outflows on the interstellar medium. Small mass outflow rates and large star formation rates of the combined sample support that there is no evidence of rapid star formation quenching by outflows, which is consistent with the delayed AGN feedback.

Keywords: agn host galaxies; star formation; quasars
The SN 2023ixf Progenitor in M101. I. Infrared Variability

Observational evidence points to a red supergiant (RSG) progenitor for SN 2023ixf. The progenitor candidate has been detected in archival images at wavelengths (≥0.6 μm) where RSGs typically emit profusely. This object is distinctly variable in the infrared (IR). We characterize the variability using pre-explosion mid-IR (3.6 and 4.5 μm) Spitzer and ground-based near-IR (JHK s ) archival data jointly covering 19 yr. The IR light curves exhibit significant variability with rms amplitudes in the range 0.2-0.4 mag, increasing with decreasing wavelength. From a robust period analysis of the more densely sampled Spitzer data, we measure a period of 1091 ± 71 days. We demonstrate using Gaussian process modeling that this periodicity is also present in the near-IR light curves, thus indicating a common physical origin, which is likely pulsational instability. We use a period-luminosity relation for RSGs to derive a value of M K = -11.58 ± 0.31 mag. Assuming a late M spectral type, this corresponds to $\mathrm{log}(L/{L}_{\odot })=5.27\pm 0.12$ at T eff = 3200 K and to $\mathrm{log}(L/{L}_{\odot })=5.37\pm 0.12$ at T eff = 3500 K. This gives an independent estimate of the progenitor's luminosity, unaffected by uncertainties in extinction and distance. Assuming the progenitor candidate underwent enhanced dust-driven mass loss during the time of these archival observations, and using an empirical period-luminosity-based mass-loss prescription, we obtain a mass-loss rate of around (2-4) × 10-4 M ⊙ yr-1. Comparing the above luminosity with stellar evolution models, we infer an initial mass for the progenitor candidate of 20 ± 4 M ⊙, making this one of the most massive progenitors for a Type II SN detected to date.

Keywords: supernovae; massive stars; circumstellar dust; variable stars
Photometric Analysis of the OGLE Heartbeat Stars

We present an analysis of 991 heartbeat stars (HBSs) from the OGLE Collection of Variable Stars. The sample consists of 512 objects located toward the Galactic bulge, 439 in the Large Magellanic Cloud, and 40 in the Small Magellanic Cloud. We model the I-band OGLE light curves using an analytical model of flux variations reflecting tidal deformations between stars. We present distributions of the model parameters that include the eccentricity, orbital inclination, and argument of the periastron but also the period-amplitude diagrams. On the Hertzsprung-Russell diagram, our HBS sample forms two separate groups of different evolutionary status. The first group, including about 90 systems with short orbital periods (P ≲ 50 days), consists of an early-type primary star lying on (or close to) the main sequence. The second group, including about 900 systems with long orbital periods (P ≳ 100 days), contains a red giant (RG). The position of the RG HBSs on the period-luminosity diagram strongly indicates their binary nature. They appear to be a natural extension of confirmed binary systems that include the OGLE ellipsoidal and long secondary period variables. We also present a time-series analysis leading to detection of tidally excited oscillations (TEOs). We identify such pulsations in about 5% of stars in the sample with a total of 78 different modes. This first relatively large homogeneous sample of TEOs allowed us to construct a diagram revealing the correlation between the TEO's orbital harmonic number and the eccentricity of the host binary system.

Keywords: binary stars; tidal distortion; time domain astronomy; elliptical orbits; galactic bulge; magellanic clouds; stellar oscillations; astronomy data modeling; periodic variable stars; astronomical simulations
Survivor Bias: Divergent Fates of the Solar System's Ejected versus Persisting Planetesimals

The orbital architecture of the solar system is thought to have been sculpted by a dynamical instability among the giant planets. During the instability a primordial outer disk of planetesimals was destabilized and ended up on planet-crossing orbits. Most planetesimals were ejected into interstellar space, but a fraction were trapped on stable orbits in the Kuiper Belt and Oort cloud. We use a suite of N-body simulations to map out the diversity of planetesimals' dynamical pathways. We focus on two processes: tidal disruption from very close encounters with a giant planet, and loss of surface volatiles from repeated passages close to the Sun. We show that the rate of tidal disruption is more than a factor of 2 higher for ejected planetesimals than for surviving objects in the Kuiper Belt or Oort cloud. Ejected planetesimals are preferentially disrupted by Jupiter and surviving ones by Neptune. Given that the gas giants contracted significantly as they cooled but the ice giants did not, taking into account the thermal evolution of the giant planets decreases the disruption rate of ejected planetesimals. The frequency of volatile loss and extinction is far higher for ejected planetesimals than for surviving ones and is not affected by the giant planets' contraction. Even if all interstellar objects were ejected from solar system-like systems, our analysis suggests that their physical properties should be more diverse than those of solar system small bodies as a result of their divergent dynamical histories. This is consistent with the characteristics of the two currently known interstellar objects.

Keywords: comets; solar system formation; tidal disruption; comet volatiles
Halo Environment for Population III Star Formation

We analyze the environmental features and growth history of high-redshift halos from cosmological simulation data to determine the conditions that lead to Population III star formation. We use support-vector classification (SVC) to determine the separation in feature-space between Population III halos and halos that form no stars. We perform hyperparameter tuning but are unable to produce an SVC model that reliably classifies Population III halos. We perform feature selection and determine that among our included features, Lyman-Werner radiation and halo mass have the most significant impact on an SVC model's effectiveness.

Keywords: population iii stars; galaxy formation
CO(J = 1-0) Mapping Survey of 64 Galaxies in the Fornax Cluster with the ALMA Morita Array

We conduct a 12C16O(J = 1-0) (hereafter CO) mapping survey of 64 galaxies in the Fornax cluster using the Atacama Large Millimeter/submillimeter Array Morita array in cycle 5. CO emission is detected from 23 out of the 64 galaxies. Our sample includes dwarf, spiral, and elliptical galaxies with stellar masses of M star ~ 106.3-11.6 M ⊙. The achieved beam size and sensitivity are 15″ × 8″ and ~12 mJy beam-1 at the velocity resolution of ~10 km s-1, respectively. We study the cold gas (molecular and atomic gas) properties of 38 subsamples with M star > 109 M ⊙ combined with literature H I data. We find that (1) the low star formation (SF) activity in the Fornax galaxies is caused by the decrease in the cold gas mass fraction with respect to stellar mass (hereafter, gas fraction) rather than the decrease of the SF efficiency from the cold gas; (2) the atomic gas fraction is more heavily reduced than the molecular gas fraction of such galaxies with low SF activity. A comparison between the cold gas properties of the Fornax galaxies and their environmental properties suggests that the atomic gas is stripped tidally and by the ram pressure, which leads to the molecular gas depletion with an aid of the strangulation and consequently SF quenching. Preprocesses in the group environment would also play a role in reducing cold gas reservoirs in some Fornax galaxies.

Keywords: galaxy environments; galaxy clusters; molecular gas