Papers for Friday, Oct 20 2023

SCALES (Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy) is a 2 - 5 micron high-contrast lenslet-based integral field spectrograph (IFS) designed to characterize exoplanets and their atmospheres. The SCALES medium-spectral-resolution mode uses a lenslet subarray with a 0.34 x 0.36 arcsecond field of view which allows for exoplanet characterization at increased spectral resolution. We explore the sensitivity limitations of this mode by simulating planet detections in the presence of realistic noise sources. We use the SCALES simulator scalessim to generate high-fidelity mock observations of planets that include speckle noise from their host stars, as well as other atmospheric and instrumental noise effects. We employ both angular and reference differential imaging as methods of disentangling speckle noise from the injected planet signals. These simulations allow us to assess the feasibility of speckle deconvolution for SCALES medium resolution data, and to test whether one approach outperforms another based on planet angular separations and contrasts.

We present a dataset for investigating the impact of stellar activity on astrometric measurements using NASA's Solar Dynamics Observatory (SDO) images of the Sun. The sensitivity of astrometry for detecting exoplanets is limited by stellar activity (e.g. starspots), which causes the measured "center of flux" of the star to deviate from the true, geometric, center, producing false positive detections. We analyze Helioseismic and Magnetic Imager continuum image data obtained from SDO between July 2015 and December 2022 to examine this "astrometric jitter" phenomenon for the Sun. We employ data processing procedures to clean the images and compute the time series of the sunspot-induced shift between the center of flux and the geometric center. The resulting time series show quasiperiodic variations up to 0.05% of the Sun's radius at its rotation period.

The overabundance of super-early (redshift $z>10$), luminous ($M_{\rm UV} < -20$), and blue galaxies detected by JWST has been explained (Ferrara et al. 2023) as due to negligible dust attenuation in these systems. We show that such model correctly reproduces the UV luminosity function at $z>10$, and the star formation rate (SFR) density evolution. The model also predicts, in agreement with data, that the cosmic specific SFR grows as ${\rm sSFR} \propto (1+z)^{3/2}$. At $z \simeq 10$ the cosmic sSFR crosses the critical value $\rm sSFR^\star = 25\, \rm Gyr^{-1}$ and $\approx 45$% of the galaxies become super-Eddington driving outflows reaching velocities of $\approx 830 \,(\epsilon_\star/f_M)^{1/2}\, {\rm km\, s}^{-1}$, where $\epsilon_\star$ and $f_M$ are the SF efficiency and fraction of the halo gas expelled in the outflow, respectively. This prediction is consistent with the outflow velocities measured in 12 super-Eddington galaxies of the JWST/JADES sample. Such outflows clear the dust, thus boosting the galaxy luminosity. They also dramatically enhance the visibility of the Ly$\alpha$ line from $z>10$ galaxies, by introducing a velocity offset. The observed Ly$\alpha$ properties in GN-z11 ($z=10.6$) are simultaneously recovered by the outflow model if $\log N_{\rm HI} \simeq 20.1$, implying that the outflow is largely ionized. We make analogous predictions for the Ly$\alpha$ visibility of other super-early galaxies, and compare the model with Ly$\alpha$ surveys at $z>7$, finding that essentially all super-Eddington (sub-Eddington) galaxies are (not) detected in Ly$\alpha$. Finally, the sSFR positively correlates with the LyC escape fraction as outflows carve ionized, transparent channels through which LyC photons leak.

We present the stellar mass - stellar metallicity relation for 3491 star-forming galaxies at $2 \lesssim z \lesssim 3$ using rest-frame far-ultraviolet (FUV) spectra from the Ly$\alpha$ Tomography IMACS Survey (LATIS). We fit stellar population synthesis models from the Binary Population And Spectral Synthesis code (BPASS v$2.2.1$) to medium resolution (R $\sim 1000$) and high signal-to-noise ($>30$ per 100 km/s over a wavelength range of 1221 - 1800 \r{A}) composite spectra of galaxies in bins of stellar mass to determine their stellar metallicity, primarily tracing $\rm Fe/H$. We find a strong correlation between stellar mass and stellar metallicity, with stellar metallicity monotonically increasing with stellar mass at low masses and flattening at high masses ($M_* \gtrsim 10^{10.3} M_\odot$). Additionally, we compare our stellar metallicity measurements with the gas-phase oxygen abundance of galaxies at similar redshift and estimate the average $\rm [\alpha/Fe] \sim 0.6$. Such high $\alpha$-enhancement indicates that high-redshift galaxies have not yet undergone significant iron enrichment through Type Ia supernovae. Moreover, we utilize an analytic chemical evolution model to constrain the mass loading parameter of galactic winds as a function of stellar mass. We find that as the stellar mass increases, the mass loading parameter decreases. The parameter then flattens or reaches a turning point at around $M_* \sim 10^{10.5} M_\odot$. Our findings may signal the onset of black hole-driven outflows at $z \sim 2.5$ for galaxies with $M_* \gtrsim 10^{10.5} M_\odot$.

We present a population of 11 of the faintest ($> 25.5$ AB mag) short gamma-ray burst (GRB) host galaxies. We model their sparse available observations using the stellar population inference code Prospector-$\beta$ and develop a novel implementation to incorporate the galaxy mass-radius relation. Assuming these hosts are randomly drawn from the galaxy population and conditioning this draw on their observed flux and size in few photometric bands, we determine that these hosts have dwarf galaxy stellar masses of $7.0\lesssim\log(M_*/M_\odot)\lesssim9.1$. This is striking as only $14\%$ of short GRB hosts with previous inferred stellar masses had $M_* \lesssim 10^{9}\,M_{\odot}$. We further show these short GRBs have smaller physical and host-normalized offsets than the rest of the population, suggesting that the majority of their neutron star (NS) merger progenitors were retained within their hosts. The presumably shallow potentials of these hosts translate to small escape velocities of $\sim5.5-80$ km/s, indicative of either low post-supernova systemic velocities or short inspiral times. While short GRBs with identified dwarf host galaxies now comprise $\approx 14\%$ of the total Swift-detected population, a number are likely missing in the current population, as larger systemic velocities (observed from Galactic NS population) would result in highly offset short GRBs and less secure host associations. However, the revelation of a population of short GRBs retained in low-mass host galaxies offers a natural explanation for observed $r$-process enrichment via NS mergers in Local Group dwarf galaxies, and has implications for gravitational wave follow-up strategies.

Interstellar dust at high Galactic latitudes can influence astronomical foreground subtraction, produce diffuse scattered light, and result in harder (de-reddened) ultraviolet spectra of quasars. In a sample of 94 quasars at high latitude and low extinction, we evaluate the interstellar "gas-to-dust ratio", $N_{\rm H}/E(B-V)$, using hydrogen column densities (H I and H$_2$) and far-infrared estimates of dust reddening. In the Galactic plane, this ratio is $6.0\pm0.2$ (in units of $10^{21}~{\rm cm}^{-2}~{\rm mag}^{-1}$). In a sub-sample of 51 quasars with measurements of both H I and H$_2$ and $0.01 \leq E(B-V) \lesssim 0.1$, we find mean ratios $10.3\pm0.4$ (gas at all velocities) and $9.2\pm0.3$ (low velocity only). Including H$_2$ increases $N_{\rm H}$ by $\sim10$%, even at high latitude. On average, recent Planck estimates of $E(B-V)$ in low reddening sight lines are 12% higher than those from Schlafly & Finkbeiner (2011), and $N_{\rm HI}$ exhibits significant variations when measured at different radio telescopes. We also show that $E(B-V)$ is sensitive to far-infrared modeled dust temperature $T_d$ and emissivity index $\beta$. Consequently, gas-to-dust ratios have large, asymmetric errors at low $E(B-V)$. The ratios are elevated in sight lines with high-velocity clouds, which contribute $N_{\rm H}$ but little reddening. In Complex C, the ratio decreases by 40% when high velocity gas is excluded. Decreases in dust content are expected in low-metallicity gas above the Galactic plane, resulting from grain destruction in shocks, settling to the disk, and thermal sputtering in hot halo gas.

In 2016 the Compton Spectrometer and Imager (COSI) had a successful 46-day flight onboard NASA's Super Pressure Balloon platform. In this work we report measurements of the Galactic diffuse continuum emission (GDCE) observed towards the inner Galaxy during the flight, which in the COSI energy band (0.2 - 5 MeV) is primarily generated from inverse Compton radiation. Within uncertainties we find overall good agreement with previous measurements from INTEGRAL/SPI and COMPTEL. Based on these initial findings, we discuss the potential for further probing the GDCE with the 2016 COSI balloon data, as well as prospects for the upcoming satellite mission.

Pulsar timing arrays (PTAs) perform Bayesian posterior inference with expensive MCMC methods. Given a dataset of ~10-100 pulsars and O(10^3) timing residuals each, producing a posterior distribution for the stochastic gravitational wave background (SGWB) can take days to a week. The computational bottleneck arises because the likelihood evaluation required for MCMC is extremely costly when considering the dimensionality of the search space. Fortunately, generating simulated data is fast, so modern simulation-based inference techniques can be brought to bear on the problem. In this paper, we demonstrate how conditional normalizing flows trained on simulated data can be used for extremely fast and accurate estimation of the SGWB posteriors, reducing the sampling time from weeks to a matter of seconds.

Ultra-high-energy (UHE) cosmic neutrinos, with energies above 100~PeV, could be finally discovered in the near future. Measuring their flavor composition would reveal information about their production and propagation, but it is unclear if planned UHE neutrino telescopes will have the capabilities to do it. We circumvent this by proposing a new way to measure the UHE neutrino flavor composition that does not rely on individual telescopes having flavor-identification capabilities. We manufacture flavor sensitivity from the joint detection by one telescope sensitive to all flavors -- the radio array of IceCube-Gen2 -- and one mostly sensitive to $\nu_\tau$ -- GRAND. From this limited flavor sensitivity, predominantly to $\nu_\tau$, and even under conservative choices of neutrino flux and detector size, we extract important insight. For astrophysics, we forecast meaningful constraints on the neutrino production mechanism; for fundamental physics, vastly improved constraints on Lorentz-invariance violation. These are the first measurement forecasts of the UHE $\nu_\tau$ content.

Dark matter-dominated cores have long been claimed for the well-studied local group dwarf galaxies. More recently, extended stellar halos have been uncovered around several of these dwarfs through deeper imaging and spectroscopy. Such core-halo structures are not a feature of conventional cold dark matter (CDM), based on collisionless particles where smooth, scale-free profiles are predicted. In contrast, smooth and prominent dark matter cores are predicted for Warm and Fuzzy/Wave Dark Matter (WDM/$\psi$DM) respectively. The question arises to what extent the visible stellar profiles should reflect this dark matter core structure. Here we compare cosmological hydrodynamical simulations of CDM, WDM $\&$ $\psi$DM, aiming to predict the stellar profiles for these three DM scenarios. We show that cores surrounded by extended halos are distinguishable for WDM and $\psi$DM, with the most prominent cores in the case of $\psi$DM, where the stellar density is enhanced in the core due to the presence of the relatively dense soliton. Our analysis demonstrates that such behavior does not appear in CDM, implying that the small-scale cut-off in the power spectrum present for WDM and $\psi$DM provides a core-halo transition. Consequently, we estimate the mass of the $\psi$DM particle at this core-halo transition point. Furthermore, we observe the anticipated asymmetry for $\psi$DM due to the soliton's random walk, a distinctive characteristic not found in the symmetric distributions of stars in Warm and CDM models.

We study the tomographic clustering properties of the photometric cluster catalogue derived from the Third Data Release of the Kilo Degree Survey, focusing on the angular correlation function and its spherical harmonic counterpart, the angular power spectrum. We measure the angular correlation function and power spectrum from a sample of 5162 clusters, with an intrinsic richness $\lambda^*\geq 15$, in the photometric redshift range $z\in [0.1, 0.6]$, comparing our measurements with theoretical models, in the framework of the $\Lambda$-Cold Dark Matter cosmology. We perform a Monte Carlo Markov Chain analysis to constrain the cosmological parameters $\Omega_{\mathrm{m}}$, $\sigma_8$ and the structure growth parameter $S_8\equiv\sigma_8 \sqrt{\Omega_{\mathrm{m}}/0.3}$. We adopt Gaussian priors on the parameters of the mass-richness relation, based on the posterior distributions derived from a previous joint analysis of cluster counts and weak lensing mass measurements carried out with the same catalogue. From the angular correlation function, we obtain $\Omega_{\mathrm{m}}=0.32^{+0.05}_{-0.04}$, $\sigma_8=0.77^{+0.13}_{-0.09}$ and $S_8=0.80^{+0.08}_{-0.06}$, in agreement, within $1\sigma$, with 3D clustering result based on the same cluster sample and with existing complementary studies on other datasets. For the angular power spectrum, we derive statistically consistent results, in particular $\Omega_{\mathrm{m}}=0.24^{+0.05}_{-0.04}$ and $S_8=0.93^{+0.11}_{-0.12}$, while the constraint on $\sigma_8$ alone is weaker with respect to the one provided by the angular correlation function, $\sigma_8=1.01^{+0.25}_{-0.17}$. Our results show that the 2D clustering from photometric cluster surveys can provide competitive cosmological constraints with respect to the full 3D clustering statistics, and can be successfully applied to ongoing and forthcoming spectro/photometric surveys.

By utilizing the spatially-resolved photometry of galaxies at $0.2<z<3.0$ in the CEERS field, we estimate the resolved and unresolved stellar mass via spectral energy distribution (SED) fitting to study the discrepancy between them. We first compare $M_{\ast}$ derived from photometry with and without the JWST wavelength coverage and find that $M_{\ast}$ can be overestimated by up to 0.2 dex when lacking rest-frame NIR data. The SED fitting process tends to overestimate both stellar age and dust attenuation in the absence of rest-frame NIR data, consequently leading to a larger observed mass-to-light ratio and hence an elevated $M_{\ast}$. With the inclusion of the JWST NIR photometry, we find no significant disparity between the resolved and unresolved stellar mass estimates, providing a plausible solution to the conflict between them out to $z\sim 3$. Further investigation demonstrates that reliable $M_{\ast}$ estimates can be obtained, regardless of whether they are derived from spatially resolved or spatially unresolved photometry, so long as the reddest filter included in the SED fitting has a rest-frame wavelength larger than 10000 \AA.

Phylogenetic methods have long been used in biology, and more recently have been extended to other fields - for example, linguistics and technology - to study evolutionary histories. Galaxies also have an evolutionary history, and fall within this broad phylogenetic framework. Under the hypothesis that chemical abundances can be used as a proxy for interstellar medium's DNA, phylogenetic methods allow us to reconstruct hierarchical similarities and differences among stars - essentially a tree of evolutionary relationships and thus history. In this work, we apply phylogenetic methods to a simulated disc galaxy obtained with a chemo-dynamical code to test the approach. We found that at least 100 stellar particles are required to reliably portray the evolutionary history of a selected stellar population in this simulation, and that the overall evolutionary history is reliably preserved when the typical uncertainties in the chemical abundances are smaller than 0.08 dex. The results show that the shape of the trees are strongly affected by the age-metallicity relation, as well as the star formation history of the galaxy. We found that regions with low star formation rates produce shorter trees than regions with high star formation rates. Our analysis demonstrates that phylogenetic methods can shed light on the process of galaxy evolution.

Using cosmological simulations of galaxy cluster regions from The Three Hundred project we study the nature of gas in filaments feeding massive clusters. By stacking the diffuse material of filaments throughout the cluster sample, we measure average gas properties such as density, temperature, pressure, entropy and Mach number and construct one-dimensional profiles for a sample of larger, radially-oriented filaments to determine their characteristic features as cosmological objects. Despite the similarity in velocity space between the gas and dark matter accretion patterns onto filaments and their central clusters, we confirm some differences, especially concerning the more ordered radial velocity dispersion of dark matter around the cluster and the larger accretion velocity of gas relative to dark matter in filaments. We also study the distribution of shocked gas around filaments and galaxy clusters, showing that the surrounding shocks allow an efficient internal transport of material, suggesting a laminar infall. The stacked temperature profile of filaments is typically colder towards the spine, in line with the cosmological rarefaction of matter. Therefore, filaments are able to isolate their inner regions, maintaining lower gas temperatures and entropy. Finally, we study the evolution of the gas density-temperature phase diagram of our stacked filament, showing that filamentary gas does not behave fully adiabatically through time but it is subject to shocks during its evolution, establishing a characteristic z = 0, entropy-enhanced distribution at intermediate distances from the spine of about 1 - 2 $h^{-1}$ Mpc for a typical galaxy cluster in our sample.

This papers searches for evidence of mass concentrations along the path of radio pulses in the PPTA2 survey data release. Radio pulse travel times are influenced via gravitational fields along the path from the source to the observer. Transient time delays in transit are a useful measure of the matter distribution along the path. Many pulsars have very well understood timing solutions with predicable arrival times and can be used to sample the mass variation. Changes in the source, observer and mass concentration positions produce changes in arrival times which can be significant for precision pulsar times. Twelve candidates are reported from this search.

In this paper, we perform a dynamical study of the population of objects in the unstable quasi-Hilda region. The aim of this work is to make an update of the population of quasi-Hilda comets (QHCs) that have recently arrived from the Centaurs region. To achieve our goal, we have applied a dynamical criteria to constrain the unstable quasi-Hilda region that allowed us to select 828 potential candidates. The orbital data of the potential candidates was take from the ASTORB database and we apply backward integration to search by those that have recently arrived from the outer regions of the Solar System. Then we studied the dynamical evolution of the candidates from a statistical point of view by calculating the time-averaged distribution of a number of clones of each candidate as a function of aphelion and perihelion distances. We found that 47 objects could have been recently injected into the inner Solar System from the Centaur or transneptunian regions. These objects may have preserved volatile material and are candidates to exhibit cometary activity.

The complex chemistry that occurs in star-forming regions can provide insight into the formation of prebiotic molecules at various evolutionary stages of star formation. To study this process, we present millimeter-wave interferometric observations of the neighboring hot cores W3(H$_2$O) and W3(OH) carried out using the NOEMA interferometer. We have analyzed distributions of six molecules that account for most observed lines across both cores and have constructed physical parameter maps for rotational temperature, column density, and velocity field with corresponding uncertainties. We discuss the derived spatial distributions of these parameters in the context of the physical structure of the source. We propose the use of HCOOCH$_3$ as a new temperature tracer in W3(H$_2$O) and W3(OH) in addition to the more commonly used CH$_3$CN. By analyzing the physically-derived parameters for each molecule across both W3(H$_2$O) and W3(OH), the work presented herein further demonstrates the impact of physical environment on hot cores at different evolutionary stages.

During the Extreme Universe Space Observatory on a Super Pressure Balloon 2 (EUSO-SPB2) mission, we planned Target of Opportunity (ToO) operations to follow up on possible sources of $\gtrsim 10 \, {\rm PeV}$ neutrinos. The original plan before flight was to point the onboard Cherenkov Telescope (CT) to catch the source's path on the sky just below Earth's horizon. By using the Earth as a tau-neutrino to tau-lepton converter, the CT would then be able to look for optical extensive air shower signals induced by tau-lepton decays in the atmosphere. The CT had a field of view of $6.4^\circ$ vertical $\times$ $12.8^\circ$ horizontal. Possible neutrino source candidates include gamma ray bursts, tidal disruption events and other bursting or flaring sources. In addition, follow-ups of binary neutron star mergers would have been possible after the start of the O4 observation run from LIGO-Virgo-KAGRA. The resulting exposure is modeled using the NuSpaceSim framework in ToO mode. With the launch of the EUSO-SPB2 payload on the 13th May 2023, this summarizes the ToO program status and preliminary data, as available.

New reference wavelengths for atomic transitions of Mg VII and Si VII in the 272-281 A wavelength range are derived using measurements from the Extreme ultraviolet Imaging Spectrometer (EIS) on board the Hinode spacecraft. Mg VII and Si VII are important ions for measuring plasma properties in the solar transition region at around 0.6 MK. The six Si VII wavelengths are 13--21 mA and 7--11 mA longer than the values in the NIST Atomic Spectra Database (ASD) and the compilations of B. Edlen, respectively. The four Mg VII wavelengths are shorter than the values in the ASD by 8-12 mA but show reasonable agreement with the Edlen values. The new wavelengths will lead to more accurate Doppler shift measurements from the EIS instrument, and will be valuable for spectral disambiguation modeling for the upcoming Multi-Slit Solar Explorer mission.

Understanding the co-evolution of supermassive black holes (SMBHs) and their host systems requires a comprehensive census of active galactic nuclei (AGN) behavior across a wide range of redshift, luminosity, obscuration level and galaxy properties. We report significant progress with JWST towards this goal from the Systematic Mid-infrared Instrument Legacy Extragalactic Survey (SMILES). Based on comprehensive SED analysis of 3273 MIRI-detected sources, we identify 217 AGN candidates over a survey area of $\sim$34 arcmin$^2$, including a primary sample of 111 AGNs in normal massive galaxies ($M_{*}>10^{9.5}~M_\odot$) at $z\sim$0--4, an extended sample of 86 AGN {\it candidates} in low-mass galaxies ($M_{*}<10^{9.5}~M_\odot$) and a high-$z$ sample of 20 AGN {\it candidates} at $z\sim$4--8.4. Notably, about 80\% of our MIRI-selected AGN candidates are new discoveries despite the extensive pre-JWST AGN searches. Even among the massive galaxies where the previous AGN search is believed to be thorough, 34\% of the MIRI AGN identifications are new, highlighting the impact of obscuration on previous selections. By combining our results with the efforts at other wavelengths, we build the most complete AGN sample to date and examine the relative performance of different selection techniques. We find the obscured AGN fraction increases from $L_{\rm AGN, bol}\sim10^{10}~L_\odot$ to $10^{11}~L_\odot$ and then drops towards higher luminosity. Additionally, the obscured AGN fraction gradually increases from $z\sim0$ to $z\sim4$ with most high-$z$ AGNs obscured. We discuss how AGN obscuration, intrinsic SED variations, galaxy contamination, survey depth and selection techniques complicate the construction of a complete AGN sample.

RZ Piscium (RZ Psc) is well-known in the variable star field because of its numerous, irregular optical dips in the past five decades, but the nature of the system is heavily debated in the literature. We present multiyear infrared monitoring data from Spitzer and WISE to track the activities of the inner debris production, revealing stochastic infrared variability as short as weekly timescales that is consistent with destroying a 90-km-size asteroid every year. ALMA 1.3 mm data combined with spectral energy distribution modeling show that the disk is compact ($\sim$0.1--13 au radially) and lacks cold gas. The disk is found to be highly inclined and has a significant vertical scale height. These observations confirm that RZ Psc hosts a close to edge-on, highly perturbed debris disk possibly due to migration of recently formed giant planets which might be triggered by the low-mass companion RZ Psc B if the planets formed well beyond the snowlines.

The goal of this work is to probe the total mass distribution of early-type galaxies with globular clusters (GCs) as kinematic tracers, by constraining the parameters of the profile with a flexible modelling approach. To that end, we leverage the extended spatial distribution of GCs from the SLUGGS survey ($\langle R_{\rm GC,\ max} \rangle \sim 8R_{\rm e}$) in combination with discrete dynamical modelling. We use discrete Jeans anisotropic modelling in cylindrical coordinates to determine the velocity moments at the location of the GCs in our sample. We use a Bayesian framework to determine the best-fit parameters of the total mass density profile and orbital properties of the GC systems. We find that the orbital properties (anisotropy and rotation of the dispersion-dominated GC systems) minimally impact the measurements of the inner slope and enclosed mass, while a strong presence of dynamically-distinct subpopulations or low numbers of kinematic tracers can bias the results. Owing to the large spatial extent of the tracers our method is sensitive to the intrinsic inner slope of the total mass profile and we find $\overline{\alpha} = -1.88\pm 0.01$ for 12 galaxies with robust measurements. To compare our results with literature values we fit a single power-law profile to the resulting total mass density. In the radial range 0.1-4~$R_{\rm e}$ our measured slope has a value of $\langle \gamma_{\rm tot}\rangle = -2.22\pm0.14$ and is in good agreement with the literature.

We summarize the properties and initial data release of the JADES Origins Field (JOF), which will soon be the deepest imaging field yet observed with the James Webb Space Telescope (JWST). This field falls within the GOODS-S region about 8' south-west of the Hubble Ultra Deep Field (HUDF), where it was formed initially in Cycle 1 as a parallel field of HUDF spectroscopic observations within the JWST Advanced Deep Extragalactic Survey (JADES). This imaging will be greatly extended in Cycle 2 program 3215, which will observe the JOF for 5 days in six medium-band filters, seeking robust candidates for z>15 galaxies. This program will also include ultra-deep parallel NIRSpec spectroscopy (up to 104 hours on-source, summing over the dispersion modes) on the HUDF. Cycle 3 observations from program 4540 will add 20 hours of NIRCam slitless spectroscopy to the JOF. With these three campaigns, the JOF will be observed for 380 open-shutter hours with NIRCam using 15 imaging filters and 2 grism bandpasses. Further, parts of the JOF have deep 43 hr MIRI observations in F770W. Taken together, the JOF will soon be one of the most compelling deep fields available with JWST and a powerful window into the early Universe. This paper presents the second data release from JADES, featuring the imaging and catalogs from the year 1 JOF observations.

The size of the optics used in observatories is often limited by fabrication, metrology, and handling technology, but having a large primary mirror provides significant benefits for scientific research. The evolution of rocket launch options enables heavy payload carrying on orbit and outstretching the telescope's form-factor choices. Moreover, cost per launch is lower than the traditional flight method, which is obviously advantageous for various novel space observatory concepts. The University of Arizona has successfully fabricated many large-scale primary optics for ground-based observatories including the Large Binocular Telescope (LBT, 8.4 meter diameter two primary mirrors), Large Synoptic Survey Telescope (now renamed to Vera C. Rubin Observatory, 8.4 meter diameter monolithic primary and tertiary mirror), and the Giant Magellan Telescope (GMT, 8.4 meter diameter primary mirror seven segments). Launching a monolithic primary mirror into space could bypass many of the difficulties encountered during the assembly and deployment of the segmented primary mirrors. However, it might bring up unprecedented challenges and hurdles, also. We explore and foresee the expected challenges and evaluate them. To estimate the tolerance and optical error budget of a large optical system in space such as three mirror anastigmat telescope, we have developed a methodology that considers various errors from design, fabrication, assembly, and environmental factors.

Much of the exoplanet discovery efforts over the next several years are largely tasked with finding candidates for the upcoming Ariel mission. This is a role that TESS is well-suited for. Radial velocity follow-up is needed to confirm the planetary nature of these systems, as many of its planet candidates turn out to be eclipsing binary systems with small stellar secondaries. Focused Doppler follow-up to obtain these radial velocities is expensive. The Gaia mission's radial velocity measurements on its target stars are not adequately precise for a general search for orbiting planetary companions. While the RV data has yielded thousands of spectroscopic binary systems, only an extreme minority of these reach into the planetary-mass regime. Nonetheless, they do detect eclipsing binary systems that can masquerade as transiting hot Jupiters. In this work, we compare the Gaia DR3 Non-Single Stars catalogue to the current planet candidate list from TESS and determine several of these candidate planetary systems are actually eclipsing binaries. We find $\sim130$ eclipsing binaries among the TOI and EPIC lists that do not appear to be in the literature, including the previously statistically validated planet K2-256 b. This work illustrates the usefulness of \textit{Gaia} to exoplanet validation efforts and will help guide follow-up efforts for discovering new transiting planets for the Ariel mission by avoiding misdirecting valuable telescope time.

The distribution of Io's volcanic activity likely reflects the position and magnitude of internal tidal heating. We use new observations of Io's polar regions by the Juno spacecraft Jovian Infrared Auroral Mapper (JIRAM) to complete near-infrared global coverage, revealing the global distribution and magnitude of thermal emission from Io's currently erupting volcanoes. We show that the distribution of volcanic heat flow from 266 active hot spots is consistent with the presence of a global magma ocean, and/or shallow asthenospheric heating. We find that Io's polar volcanoes are less energetic but about the same in number per unit area than at lower latitudes. We also find that volcanic heat flow in the north polar cap is greater than that in the south. The low volcanic advection seen at Io's poles is therefore at odds with measurements of background temperature showing Io's poles are anomalously warm. We suggest that the differences in volcanic thermal emission from Io's poles compared to that at lower latitudes is indicative of lithospheric dichotomies that inhibit volcanic advection towards Io's poles, particularly in the south polar region.

Protostellar disks are a ubiquitous part of the star formation process and the future sites of planet formation. As part of the Early Planet Formation in Embedded Disks (eDisk) large program, we present high-angular resolution dust continuum ($\sim40\,$mas) and molecular line ($\sim150\,$mas) observations of the Class 0 protostar, IRAS 15398-3359. The dust continuum is small, compact, and centrally peaked, while more extended dust structures are found in the outflow directions. We perform a 2D Gaussian fitting to find the deconvolved size and $2\sigma$ radius of the dust disk to be $4.5\times2.8\,\mathrm{au}$ and $3.8\,\mathrm{au}$, respectively. We estimate the gas+dust disk mass assuming optically thin continuum emission to be $0.6-1.8\,M_\mathrm{jup}$, indicating a very low-mass disk. The CO isotopologues trace components of the outflows and inner envelope, while SO traces a compact, rotating disk-like component. Using several rotation curve fittings on the PV diagram of the SO emission, the lower limits of the protostellar mass and gas disk radius are $0.022\,M_\odot$ and $31.2\,\mathrm{au}$ from our Modified 2 single power-law fitting. A conservative upper limit of the protostellar mass is inferred to be $0.1\,M_\odot$. The protostellar mass-accretion rate and the specific angular momentum at the protostellar disk edge are found to be between $1.3-6.1\times10^{-6}\,M_\odot\,\mathrm{yr^{-1}}$ and $1.2-3.8\times10^{-4}\,\mathrm{km\,s^{-1}\,pc}$, respectively, with an age estimated between $0.4-7.5\times10^{4}\,$yr. At this young age with no clear substructures in the disk, planet formation would likely not yet have started. This study highlights the importance of high-resolution observations and systematic fitting procedures when deriving dynamical properties of deeply embedded Class 0 protostars.

We report the confirmation of a sub-Saturn-size exoplanet, TOI-1194 b with a mass about $0.456_{-0.051}^{+0.055}$ $M_{J}$, and a very low mass companion star with a mass of about $96.5\pm1.5$ $M_J$, TOI-1251 B. Exoplanet candidates provided by the Transiting Exoplanet Survey Satellite (TESS) are suitable for further follow-up observations by ground-based telescopes with small and medium apertures. The analysis is performed based on data from several telescopes worldwide, including telescopes in the Sino-German multiband photometric campaign, which aimed at confirming TESS Objects of Interest (TOIs) using ground-based small-aperture and medium-aperture telescopes, especially for long-period targets. TOI-1194 b is confirmed based on the consistent periodic transits depths from the multiband photometric data. We measure an orbital period of $2.310644\pm0.000001$ d, and radius is $0.767_{-0.041}^{+0.045}$ $R_J$, and amplitude of RV curve is $69.4_{-7.3}^{+7.9}$ m/s. TOI-1251 B is confirmed based on the multiband photometric and high-resolution spectroscopic data, whose orbiting period is $5.963054_{-0.000001}^{+0.000002}$ d, the radius is $0.947_{-0.033}^{+0.035}$ $R_J$, and amplitude of RV curve is $9849_{-40}^{+42}$ m/s.

The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) instrument is a lenslet-based integral field spectrograph that will operate at 2 to 5 microns, imaging and characterizing colder (and thus older) planets than current high-contrast instruments. Its spatial resolution for distant science targets and/or close-in disks and companions could be improved via interferometric techniques such as sparse aperture masking. We introduce a nascent Python package, NRM-artist, that we use to design several SCALES masks to be non-redundant and to have uniform coverage in Fourier space. We generate high-fidelity mock SCALES data using the scalessim package for SCALES' low spectral resolution modes across its 2 to 5 micron bandpass. We include realistic noise from astrophysical and instrument sources, including Keck adaptive optics and Poisson noise. We inject planet and disk signals into the mock datasets and subsequently recover them to test the performance of SCALES sparse aperture masking and to determine the sensitivity of various mask designs to different science signals.

We present the results of observations of asteroid (248370) QN$_{173}$ obtained during July 2021 - January 2022 with three telescopes. Our analysis revealed the presence of the dust tail for about half of a year. The direct images of the asteroid were obtained with broad-band filters. No emissions were revealed in the spectra, and the spectrum of the asteroid closely matched that of a C-type asteroid. Created color and linear polarization variations along the tail were analyzed. The asteroid demonstrated a redder color compared to the Sun. Dramatic changes in dust productivity obtained in different filters were not detected. The $g-r$ color changes from $0.2^{m}$ to $0.7^{m}$ over the coma, and the linear polarization degree varies from about $1.2$\% to $0.2$\% and from $-0.2$\% to $-1.5$\% at the phase angle of $23.2^{\circ}$ and $8.16^{\circ}$. The total dust mass ejected until the latest observation on October 10 is $4.2 \times 10^7$~kg, with a maximum rate of 2.6~kg\,s$^{-1}$ based on the Monte Carlo modeling of the dust tail. The estimated asteroid size is 1.3~km. It is shown that large particles are concentrated around the nucleus, whereas smaller ones dominate in the tail. The evolution of (248370) QN$_{173}$ orbit and the orbits of the sample of the 464 short-periodic comets were followed. Ten of them approached the asteroid's orbit. These objects are not genetically related, despite the very close distance of their orbits for a relatively long time.

We investigate the impact of inelastic collisions between dark matter (DM) and heavy cosmic ray (CR) nuclei on CR propagation. We approximate the fragmentation cross-sections for DM-CR collisions using collider-measured proton-nuclei scattering cross-sections, allowing us to assess how these collisions affect the spectra of CR Boron and Carbon. We derive new CR spectra from DM-CR collisions by incorporating these DM-CR cross-sections into the source terms and solving the diffusion equation for the complete network of reactions involved in generating secondary species. Utilizing the latest data from AMS-02 and DAMPE on the Boron-to-Carbon ratio, we estimate a 95\% upper limit for the effective inelastic cross-section of DM-proton as a function of DM mass. Our findings reveal that at $m_\chi \simeq 2 \,\rm{MeV}$, the effective inelastic cross-section between DM and protons must be less than $\mathcal{O}(10^{-32})~{\rm cm}^2$.

In the absence of a sufficient amount of plasma injection into the black hole (BH) magnetosphere, the force-free state of the magnetosphere cannot be maintained, leading to the emergence of strong, time-dependent, longitudinal electric field (spark gap). Recent studies of supermassive BH magnetospheres by using analytical methods and particle-in-cell (PIC) simulations propose the possibility of the efficient particle acceleration and consequent gamma-ray emissions in the spark gap. In this work, we perform one-dimensional general relativistic PIC simulations to examine the gamma-ray emission from stellar-mass BH magnetospheres. We find that intermittent spark gaps emerge and particles are efficiently accelerated, in a similar manner to the supermassive BH case. We build a semi-analytic model of the plasma dynamics and radiative processes which reproduces the maximum electron energies and peak gamma-ray luminosities in the simulation results. Based on this model, we show that gamma-ray signals from stellar-mass BHs wandering through the interstellar medium could be detected by gamma-ray telescopes such as the Fermi Large Area Telescope, or the Cherenkov Telescope Array.

Observing 3D magnetic fields, including orientation and strength, within the interstellar medium is vital but notoriously difficult. However, recent advances in our understanding of anisotropic magnetohydrodynamic (MHD) turbulence demonstrate that MHD turbulence and 3D magnetic fields leave their imprints on the intensity features of spectroscopic observations. Leveraging these theoretical frameworks, we propose a novel Convolutional Neural Network (CNN) model to extract this embedded information, enabling the probe of 3D magnetic fields. This model examines not only the plane-of-the-sky magnetic field orientation ($\phi$), but also the magnetic field's inclination angle ($\gamma$) relative to the line-of-sight, and the total magnetization level (M$_A^{-1}$) of the cloud. We train the model using synthetic emission lines of $^{13}$CO (J = 1 - 0) and C$^{18}$O (J = 1 - 0), generated from 3D MHD simulations that span conditions from sub-Alfv\'enic to super-Alfv\'enic molecular clouds. Our tests confirm that the CNN model effectively reconstructs the 3D magnetic field topology and magnetization. The median uncertainties are under $5^\circ$ for both $\phi$ and $\gamma$, and less than 0.2 for M$_A$ in sub-Alfv\'enic conditions (M$_A\approx0.5$). In super-Alfv\'enic scenarios (M$_A\approx2.0$), they are under $15^\circ$ for $\phi$ and $\gamma$, and 1.5 for M$_A$. We applied this trained CNN model to the L1478 molecular cloud. Results show a strong agreement between the CNN-predicted magnetic field orientation and that derived from Planck 353 GHz polarization data. The CNN approach enabled us to construct the 3D magnetic field map for L1478, revealing a global inclination angle of $\approx76^\circ$ and a global M$_A$ of $\approx1.07$.

Swift J1357.2-0933 is a transient low-mass X-ray binary hosting a stellar-mass black hole. The source exhibits optical dips and very broad emission lines during both outburst and quiescence, which are thought to be the result of a high orbital inclination. We present phase-resolved spectroscopy obtained with the 10.4m Gran Telescopio Canarias (GTC). The spectra focus on the $\rm{H}\beta$ spectral region during X-ray quiescence. The emission line is exceptionally broad (full width at half maximum, FWHM > 4000 \AA), in agreement with previous studies focused on $\rm{H}\alpha$. A two-Gaussian fit to the prominent double-peaked profile reveals a periodic variability in the centroid position of the line. We also produced a diagnostic diagram aimed at constraining additional orbital parameters. Together, they allow us to independently confirm the orbital period of the system using a new dataset obtained five years after the previous outburst. However, our estimates for both the systemic velocity and the radial velocity semi-amplitude of the black hole reveal larger values than those found in previous studies. We argue that this could be explained by the precession of the disc and the presence of a hotspot. We found evidence of a narrow inner core in the double-peaked H$\beta$ emission profile. We studied its evolution across the orbit, finding that it is likely to result from the occultation of inner material by the outer rim bulge, further supporting the high orbital inclination hypothesis.

The standard approach to obtaining knowledge about the properties of the surface layer of a comet from observations of gas production consists of two stages. First, various thermophysical models are used to calculate gas production for a few sets of parameters. Second, a comparison of observations and theoretical predictions is performed. This approach is complicated because the values of many model characteristics are known only approximately. Therefore, it is necessary to investigate the sensitivity of the simulated outgassing to variations in the properties of the surface layer. This problem was recently considered by us for aggregates up to tens of microns in size. For millimetre-size aggregates, a qualitative extension of the method used to model the structural characteristics of the layer is required. It is also necessary to study the role of radiative thermal conductivity, which may play an important role for such large particles. We investigated layers constructed from large aggregates and having various thicknesses and porosity and evaluated the effective sublimation of water ice at different heliocentric distances. For radiative conductivity, approximate commonly used models and the complicated model based on the Dense Medium Radiative Transfer theory were compared. It was shown that for millimetre-size aggregates careful consideration of the radiative thermal conductivity is required since this mechanism of energy transfer may change the resulting gas productivity by several times. We demonstrate that our model is more realistic for an evolved comet than simple models parameterising the properties of the cometary surface layer, yet maintains comparable computational complexity.

The observation of interstellar objects 1I/'Oumuamua and 2I/Borisov suggests the existence of a larger population of smaller projectiles that impact our planet with unbound orbits. We analyze an asteroidal grazing meteor (FH1) recorded by the Finnish Fireball Network on October 23, 2022. FH1 displayed a likely hyperbolic orbit lying on the ecliptic plane with an estimated velocity excess of $\sim$0.7 km$\,$s$^{-1}$ at impact. FH1 may either be an interstellar object, indicating a high-strength bias in this population, or an Oort cloud object, which would reinforce migration-based solar system models. Furthermore, under the calculated uncertainties, FH1 could potentially be associated with the passage of Scholz's binary star system. Statistical evaluation of uncertainties in the CNEOS database and study of its hyperbolic fireballs reveals an anisotropic geocentric radiant distribution and low orbital inclinations, challenging the assumption of a randomly incoming interstellar population. Orbital integrations suggest that the event on March 9, 2017 (IM2) from CNEOS may have experienced gravitational perturbation during the Scholz fly-by, contingent upon velocity overestimation within the expected range. These findings suggest that apparent interstellar meteors may, in fact, be the result of accelerated meteoroid impacts caused by close encounters with massive objects within or passing through our solar system.

A complete and pure sample of quasars with accurate redshifts is crucial for quasar studies and cosmology. In this paper, we present CatNorth, an improved Gaia DR3 quasar candidate catalog with more than 1.5 million sources in the 3$\pi$ sky built with data from Gaia, Pan-STARRS1, and CatWISE2020. The XGBoost algorithm is used to reclassify the original Gaia DR3 quasar candidates as stars, galaxies, and quasars. To construct training/validation datasets for the classification, we carefully built two different master stellar samples in addition to the spectroscopic galaxy and quasar samples. An ensemble classification model is obtained by averaging two XGBoost classifiers trained with different master stellar samples. Using a probability threshold of $p_{\mathrm{QSO\_mean}}>0.95$ in our ensemble classification model and an additional cut on the logarithmic probability density of zero proper motion, we retrieved 1,545,514 reliable quasar candidates from the parent Gaia DR3 quasar candidate catalog. We provide photometric redshifts for all candidates with an ensemble regression model. For a subset of 89,100 candidates, accurate spectroscopic redshifts are estimated with the Convolutional Neural Network from the Gaia BP/RP spectra. The CatNorth catalog has a high purity of > 90% while maintaining high completeness, which is an ideal sample to understand the quasar population and its statistical properties. The CatNorth catalog is used as the main source of input catalog for the LAMOST phase III quasar survey, which is expected to build a highly complete sample of bright quasars with $i < 19.5$.

Observations of the Lyman-$\alpha$ forest in distant quasar spectra with upcoming surveys are expected to provide significantly larger and higher-quality datasets. To interpret these datasets, it is imperative to develop efficient simulations. One such approach is based on the assumption that baryonic densities in the intergalactic medium (IGM) follow a lognormal distribution. We extend our earlier work to assess the robustness of the lognormal model of the Lyman-$\alpha$ forest in recovering the parameters characterizing IGM state, namely, the mean-density IGM temperature ($T_0$), the slope of the temperature-density relation ($\gamma$), and the hydrogen photoionization rate ($\Gamma_{12}$), by comparing with high-resolution Sherwood SPH simulations across the redshift range $2 \leq z \leq 2.7$. These parameters are estimated through a Markov Chain Monte Carlo technique, using the mean and power spectrum of the transmitted flux. We find that the usual lognormal distribution of IGM densities cannot recover the parameters of the SPH simulations. This limitation arises from the fact that the SPH baryonic density distribution cannot be described by a simple lognormal form. To address this, we extend the model by scaling the linear density contrast by a parameter $\nu$. While the resulting baryonic density is still lognormal, the additional parameter gives us extra freedom in setting the variance of density fluctuations. With this extension, values of $T_0$ and $\gamma$ implied in the SPH simulations are recovered at $\sim 1-\sigma$ ($\lesssim$ 10%) of the median (best-fit) values for most redshifts bins. However, this extended lognormal model cannot recover $\Gamma_{12}$ reliably, with the best-fit value discrepant by $\gtrsim 3-\sigma$ for $z > 2.2$. Despite this limitation in the recovery of $\Gamma_{12}$, we argue that the model remains useful for constraining cosmological parameters.

The carbonaceous Bencubbin-like (CB), high-metal (CH), and Renazzo-like (CR) chondrites are metal-rich chondrites that have been suggested to be genetically linked and are sometimes grouped together as the CR chondrite clan. Of these, the CB and CH chondrites are thought to have formed in an impact-generated vapor-melt plume from material that may be isotopically akin to CR chondrites. We report Mo, Ti, Cr, and Hf-W isotopic data for CB and CH chondrites in order to determine their formation time, to assess whether these chondrites are genetically related, and to evaluate their potential link to CR chondrites. An internal Hf-W isochron for the CH chondrite Acfer 182 yields an age of 3.8 +- 1.2 Ma after the formation of Ca-Al-rich inclusions (CAIs), which is indistinguishable from the mean Hf-W model age for CB metal of 3.8 +- 1.3 Ma. The Mo isotopic data for CB and CH chondrites indicate that both contain some of the same metal and silicate components, which themselves are isotopically distinct. As such, the different Mo isotopic compositions of bulk CB and CH chondrites reflect their distinct metal-to-silicate ratios. CR metal exhibits the same Mo isotopic composition as CB and CH metal, but CR silicates have distinct Mo and Ti isotopic compositions compared to CB and CH silicates, indicating that CB/CH chondrites may be genetically related to CR metal, but not to CR silicates. Together, the new isotopic data are consistent with formation of CB and CH chondrites in different regions of a common impact-generated vapor- melt plume and suggest that the CB and CH metal may derive from a metal-rich precursor genetically linked to CR chondrites. The Hf-W systematics of CH and CB chondrites indicate that the impact occurred at 3.8 +- 0.8 Ma after the formation of Ca-Al-rich inclusions and, hence, up to ~1 Ma earlier than previously inferred based on Pb- Pb chronology.

Blazars display variable emission across the entire electromagnetic spectrum, with timescales that can range from a few minutes to several years. Our recent work has shown that a sample of five blazars exhibit hints of periodicity with a global significance $\gtrsim2\,\sigma$ at $\gamma$-ray energies, in the range of 0.1~GeV$<$E$<$800~GeV. In this work, we study their multiwavelength (MWL) emission, covering the X-ray, ultraviolet, optical, and radio bands. We show that three of these blazars present similar periodic patterns in the optical and radio bands. Additionally, fluxes in the different bands of the five blazars are correlated, suggesting a co-spatial origin. Moreover, we detect a long-term ($\approx$10 year) rising trend in the light curves of PG~1553+113, and we use it to infer possible constraints on the binary black hole hypothesis.

Context: Observations of the hot gas in distant clusters of galaxies, though challenging, are key to understand the role of intense galaxy activity, super-massive black hole feedback and chemical enrichment in the process of massive halos assembly. Aims: We assess the feasibility to retrieve, using X-ray hyperspectral data only, the thermodymamical hot gas properties and chemical abundances of a $z=2$ galaxy cluster of mass M500=7 x $10^{13} M_{\odot}$, extracted from the Hydrangea hydrodynamical simulation. Methods: We create mock X-ray observations of the future X-ray Integral Field Unit (X-IFU) onboard the Athena mission. By forward-modeling the measured 0.4-1 keV surface brightness, the projected gas temperature and abundance profiles, we reconstruct the three-dimensional distribution for the gas density, pressure, temperature and entropy. Results: Thanks to its large field-of-view, high throughput and exquisite spectral resolution, one X-IFU exposure lasting 100ks enables reconstructing density and pressure profiles with 20% precision out to a characteristic radius of R500, accounting for each quantity's intrinsic dispersion in the Hydrangea simulations. Reconstruction of abundance profiles requires both higher signal-to-noise ratios and specific binning schemes. We assess the enhancement brought by longer exposures and by observing the same object at later evolutionary stages ($z=1-1.5$). Conclusions: Our analysis highlights the importance of scatter in the radially binned gas properties, which induces significant effects on the observed projected quantities. The fidelity of the reconstruction of gas profiles is sensitive to the degree of gas components mixing along the line-of-sight. Future analyses should aim at involving dedicated hyper-spectral models and fitting methods that are able to grasp the complexity of such three-dimensional, multi-phase, diffuse gas structures.

Next generation gravitational waves (GWs) observatories are expected to measure GW signals with unprecedented sensitivity, opening new, independent avenues to learn about our Universe. The distance-redshift relation is a fulcrum for cosmology and can be tested with GWs emitted by merging binaries of compact objects, called standard sirens, thanks to the fact that they provide the absolute distance from the source. On the other hand, fluctuations of the intervening matter density field induce modifications on the measurement of luminosity distance compared to that of a homogeneous universe. Assuming that the redshift information is obtained through the detection of an electromagnetic counterpart, we investigate the impact that lensing of GWs might have in the inference of cosmological parameters. We treat lensing as a systematic error and check for residual bias on the values of the cosmological parameters. We do so by means of mock catalogues of bright sirens events in different scenarios relevant to Einstein Telescope. For our fiducial scenario, the lensing bias can be comparable to or greater than the expected statistical uncertainty of the cosmological parameters, although non-negligible fluctuations in the bias values are observed for different realisations of the mock catalogue. We also discuss some mitigation strategies that can be adopted in the data analysis. Overall, our work highlights the need to model lensing effects when using standard sirens as probes of the distance-redshift relation.

18 mm-sized Orionid meteoroids were captured in 2019 and 2020 by the Canadian Automated Observatory's mirror tracking system. Meteor position measurements were made to an accuracy of $\sim1$ m and the meteors were tracked to a limiting magnitude of about $+7.5$ at the faintest point. The trajectory estimation shows the intrinsic physical dispersion of the Orionid radiant is $0.400^{\circ} \pm 0.062^{\circ}$. An erosion-based entry model was fit to the observations to reproduce ablation and fragmentation for each meteor, simultaneously reproducing the light curve, the dynamics, and the wake. Wake observations were found to directly inform the grain mass distribution released in the modelled erosion. A new luminous efficiency model was derived from simultaneous radar and optical observations and applied in the modelling to improve its accuracy. The results show that the apparent strength of Orionids varies with radiant location and time of appearance during the period of shower activity. The average differential grain mass distribution index was 2.15, higher than found from in-situ estimates, possibly due to the evolution of the physical properties of meteoroids since ejection. All Orionids showed leading fragment morphology which was best explained by stopping the erosion at the peak of the light curve, leaving a non-fragmenting meteoroid with $\sim10\%$ of the original mass. The inverted Orionid meteoroid average bulk density of $\sim300$ kg m$^{-3}$, corresponding to porosities of $\sim90\%$, is consistent with in-situ measurements of larger dust particles by Vega-2 at 1P/Halley and Rosetta at 67P.

The chemical fingerprint of a planet reveals information about its formation history regarding when and where it formed. The water content of a planet can help to constrain its formation pathway: If the planet formed in the outer regions of the disk and migrated inward, it will be water-rich due to the accretion of water-ice-rich solids. Conversely, formation in the inner disk, where water-ice is not available, will result in a smaller atmospheric water content due to the limited accretion of water vapor. This process complicates with the presence of gap-opening giant planets. A gas giant exerts a pressure bump exterior to its orbit, preventing the further influx of pebbles into the inner system, resulting in a water-poor environment. The different formation scenarios can help to constrain the formation of the HAT-P-11 system, which contains an inner sub-Neptune with a mass of 23.4 $\mathrm{M_{\oplus}}$ and substellar water abundances ($X_\mathrm{H_2O} \approx 0.11$), as well as an outer giant planet orbiting exterior to the water-ice line. Our planet formation model encompasses planetary growth through pebble and gas accretion, along with a pebble drift and evaporation module that enables us to track the chemical composition of the disk and the planets. We find that the presence of the gas giant is necessary to block water-ice-rich material, resulting in a substellar water content for the inner sub-Neptune. If the giant planet forms too early, not enough solid material can enter the inner disk regions, preventing the efficient growth of the inner planet. This highlights the importance of the timing of giant planet formation to explain the inner system structure. Our simulations predict a roughly stellar C/O ratio with superstellar C/H and O/H for HAT-P-11b, providing constraints for future observations that are essential for gaining a more detailed understanding of its formation.

It is commonly believed that accretion discs are truncated and their inner regions are described by advection dominated accretion flows (ADAFs) in the hard spectral state of black hole X-ray binaries. However, the increasing occurrence of a relativistically blurred Fe K$\alpha$ line together with a hard continuum points to the existence of a thin disc located near the innermost stable circular orbit (ISCO). Assuming the accretion in the hard state is via an ADAF extending to near 100 Schwarzschild radii, which is supplied by either a stellar wind from a companion star or resulting from an evaporated disc, we study the possible condensation of the hot gas during its accretion towards the black hole. It is found that a small fraction of the ADAF condenses into a cold disc as a consequence of efficient radiative cooling at small distances, forming a disc-corona configuration near the ISCO. This takes place for low accretion rates corresponding to luminosities ranging from $\sim 10^{-3}$ to a few per cent of the Eddington luminosity. The coexistence of the weak inner disc and the dominant hot accretion flow provides a natural explanation of the broad K$\alpha$ line in the hard state. Detailed computations demonstrate that such accretion flows produce a hard X-ray spectrum accompanied by a weak disc component with a negative correlation between the 2-10 keV photon index and the Eddington ratio. The predicted spectrum of Cygnus X-1 and the correlation between the photon index and the Eddington ratio are in good agreement with observations.

DQ Tau is a unique young high-eccentricity binary system that exhibits regular magnetic reconnection flares and pulsed accretion near periastron. We conducted NuSTAR, Swift, and Chandra observations during the July 30, 2022 periastron to characterize X-ray, near-ultraviolet (NUV), and optical flaring emissions. Our findings confirm the presence of X-ray super-flares accompanied by substantial NUV and optical flares, consistent with previous discoveries of periastron flares in 2010 and 2021. These observations, supported by new evidence, strongly establish the magnetosphere collision mechanism as the primary driver of magnetic energy release during DQ Tau's periastron flares. The energetics of the observed X-ray super-flares remain consistent across the three periastrons, indicating recurring energy sources during each passage, surpassing the capabilities of single stars. The observed flaring across multiple bands supports the Adams et al. model for magnetosphere interaction in eccentric binaries. Evidence from modeling and past and current observations suggests that both the mm/X-ray periastron flares and tentatively, the magnetic reconnection-related components of the optical/NUV emissions, conform to the classical solar/stellar non-thermal thick-target model, except for the distinctive magnetic energy source. However, our NuSTAR observations suffered from high background levels, hindering the detection of anticipated non-thermal hard X-rays. Furthermore, we report serendipitous discovery of X-ray super-flares occurring away from periastron, potentially associated with interacting magnetospheres. The current study is part of a broader multi-wavelength campaign, which is planned to investigate the influence of DQ Tau's stellar radiation on gas-phase ion chemistry within its circumbinary disk.

Extreme mass-ratio inspirals (EMRIs), namely binary systems composed of a massive black hole and a compact stellar-mass object, are anticipated to be among the gravitational wave (GW) sources detected by the Laser Interferometer Space Antenna (LISA). Similarly to compact binary mergers detected by current GW detectors, EMRIs can be used as cosmic rulers to probe the expansion of the Universe. Motivated by tensions in current cosmological observations as well as by alternative models of dark energy, modified gravity theories can affect the propagation of GWs across cosmological distances, with modifications commonly parametrised in terms of two phenomenological parameters, $\Xi_0$ and $n$. In this work we adopt a Bayesian approach to constrain for the first time parametrised deviations from General Relativity using the loudest simulated EMRIs detected by LISA as dark sirens with a simulated galaxy catalog. Assuming all the cosmological parameters except $\Xi_0$ are already tightly constrained, our forecasts show that $\Xi_0$ can be constrained to a few percent level (90% C.I.) with 4 years of LISA observations, unless EMRI detection rates turn out to be closer to current pessimistic expectations. These results quickly degrade if additional cosmological parameters are inferred simultaneously, but become more robust with an extended LISA observation period of 10 years. Overall, we find that EMRIs with LISA are better at constraining modified GW propagation than current second-generation ground-based GW detectors, but they will only be comparable to third-generation detectors in the most optimistic scenarios.

As IceCube surpasses a decade of operation in the full detector configuration, results that drive forward the fields of neutrino astronomy, cosmic ray physics, multi-messenger astronomy, particle physics, and beyond continue to emerge at an accelerated pace. IceCube data is dominated by background events, and thus teasing out the signal is the common challenge to most analyses. Statistical accumulation of data, along with better understanding of the background fluxes, the detector, and continued development of our analysis tools have produced many profound results that were presented at ICRC2023. Highlights covered here include the first neutrino observation of the Galactic Plane, the first observation of a steady emission neutrino point source NGC1068, new characterizations of the cosmic ray flux and its secondary particles, and a possible new era in measuring the energy spectrum of the diffuse astrophysical flux. IceCube is poised to make more discoveries and drive fields forward in the near future with many novel analyses coming online.

Extreme-ultraviolet (EUV) waves are one of the large-scale phenomena on the Sun. They are defined as large propagating fronts in the low corona with speeds ranging from a few tens km/s to a multiple of 1000 km/s. They are often associated with solar filament eruptions, flares, or coronal mass ejections (CMEs). EUV waves show different features, such as, wave and nonwave components, stationary fronts, reflection, refraction, and mode conversion. Apart from these, they can hit the nearby coronal loops and filaments/prominences during their propagation and trigger them to oscillate. These oscillating loops and filaments/prominences enable us to diagnose coronal parameters such as the coronal magnetic field strength. In this article, we present the different observed features of the EUV waves along with existing models.

(abridged) Context: Rotation is thought to influence the size of convective eddies and the efficiency of convective energy transport in the deep convection zones of stars. Rotationally constrained convection has been invoked to explain the lack of large-scale power in observations of solar flows. Aims: The main aims are to quantify the effects of rotation on the scale of convective eddies and velocity, the depths of convective overshoot, and the subadiabatic Deardorff layers. Methods: Three-dimensional hydrodynamic simulations of rotating convection in Cartesian domains were run. The results were compared with theoretical scaling results that assume a balance between Coriolis, inertial, and buoyancy (Archemedean) forces (CIA balance). Results: The scale of convective eddies decreases as rotation increases, and ultimately reaches a rotationally constrained regime consistent with the CIA balance. Using a new measure of the rotational influence on the system, it is shown that even the deep parts of the solar convection zone are not in the rotationally constrained regime. The simulations capture the slowly and rapidly rotating scaling laws predicted by theory, and the Sun appears to be in between these two regimes. Both, the overshooting depth and the extent of the Deardorff layer, decrease as rotation becomes more rapid. For sufficiently rapid rotation the Deardorff layer is absent. Conclusions: Relating the simulations with the Sun suggests that the convective scale even in the deep parts of the Sun is only mildly affected by rotation and that some other mechanism is needed to explain the lack of strong large-scale flows in the Sun. Taking the current results at face value, the overshoot and Deardorff layers are estimated to span roughly five per cent of the pressure scale height at the base of the convection zone in the Sun.

The TeV afterglow detected by LHAASO has been interpreted as arising from a narrow jet while the radio to X-ray afterglows are explained as arising from a wide structured jet. However, there is no model explaining the TeV and lower-energy multi-wavelength afterglows simultaneously. We here investigate a two-component jet model, including an inner narrow core and an outer wide wing with an angular structure, to explain both the early TeV afterglow and multi-wavelength afterglows that last up to 100 days. We find that the radio afterglow and the TeV upper limit imposed by H.E.S.S. observations combine to constrain the circum-burst density to be low at larger radii. Thus, a decreasing density profile with radius is favored. Considering that the rising TeV light curve during the afterglow onset favors a constant-density medium, we invoke a stratified density profile, including a constant-density profile at small radii and a wind density profile at large radii. We find that the two-component jet model with such a stratified density profile can explain the TeV, X-ray and optical afterglows of GRB 221009A, although the radio fluxes exceed the observed ones by a factor of two at later epochs. The discrepancy in the radio afterglow could be resolved by invoking some non-standard assumption about the microphysics of afterglow shocks, such as a decreasing fraction of accelerated particles with time. The total kinetic energy of the two components in our model is $\lesssim 10^{52}{\rm erg}$, significantly smaller than that in the single structured jet models.

The standard cosmological model predicts statistically isotropic cosmic microwave background (CMB) fluctuations. However, several summary statistics of CMB isotropy have anomalous values, including: the low level of large-angle temperature correlations, $S_{1/2}$; the excess power in odd versus even low-$\ell$ multipoles, $R^{TT}$; the (low) variance of large-scale temperature anisotropies in the ecliptic north, but not the south, $\sigma^2_{16}$; and the alignment and planarity of the quadrupole and octopole of temperature, $S_{QO}$. Individually, their low $p$-values are weak evidence for violation of statistical isotropy. The correlations of the tail values of these statistics have not to this point been studied. We show that the joint probability of all four of these happening by chance in $\Lambda$CDM is likely $\leq3\times10^{-8}$. This constitutes more than $5\sigma$ evidence for violation of statistical isotropy.

Evidence abounds that young stellar objects undergo luminous bursts of intense accretion that are short compared to the time it takes to form a star. It remains unclear how much these events contribute to the main-sequence masses of the stars. We demonstrate the power of time-series far-infrared (far-IR) photometry to answer this question compared to similar observations at shorter and longer wavelengths. We start with model spectral energy distributions that have been fit to 86 Class 0 protostars in the Orion molecular clouds. The protostars sample a broad range of envelope densities, cavity geometries, and viewing angles. We then increase the luminosity of each model by factors of 10, 50, and 100 and assess how these luminosity increases manifest in the form of flux increases over wavelength ranges of interest. We find that the fractional change in the far-IR luminosity during a burst more closely traces the change in the accretion rate than photometric diagnostics at mid-infrared and submillimeter wavelengths. We also show that observations at far-IR and longer wavelengths reliably track accretion changes without confusion from large, variable circumstellar and interstellar extinction that plague studies at shorter wavelengths. We close by discussing the ability of a proposed far-IR surveyor for the 2030s to enable improvements in our understanding of the role of accretion bursts in mass assembly.

The gamma-ray burst (GRB) GRB~211211A is believed to have occurred due to the merger of two neutron stars or a neutron star and a black hole, despite its duration of more than a minute. Subsequent analysis has revealed numerous interesting properties including the possible presence of a $\sim 22$~Hz quasiperiodic oscillation (QPO) during precursor emission. Here we perform timing analysis of Fermi and Swift gamma-ray data on GRB~211211A and, although we do not find a strong QPO during the precursor, we do find an extremely significant 19.5~Hz flux oscillation, which has higher fractional amplitude at higher energies, in a $\sim 0.2$~second segment beginning $\sim 1.6$~seconds after the start of the burst. After presenting our analysis we discuss possible mechanisms for the oscillation.

Circumnuclear star forming regions (CNSFR) are massive clusters found close to galactic nuclei. These entities give us an excellent opportunity to study star formation in environments with high metallicity and to relate it with active galactic nuclei. Our principal aim is to derive the physical properties and dynamical masses of the CNSFRs in the two rings of the spiral NGC 7469, categorized as a Luminous Infrared Galaxy (ULIRG) and hosting a Seyfert 1 nucleus. We have used archival data obtained with the MUSE spectrograph attached to one of the ESO VLT telescopes and we have applied the techniques shown in the first paper of the series. Regions in the studied galaxy show large sizes which can be explained by the stellar winds produced by WR stars. The inner ring regions seem to be more compact than the outer ones, showing higher electron densities and filling factors. The young stellar population of the clusters has contributions of ionising populations with ages around 5 Ma and its masses constitute less than a 1\% of the total mass of each cluster. The inner ring regions which are close to the active galactic nucleus probably are the only ones that have enough mass to survive the action of the AGN. They constitute the $\sim$ 90 \% of the total inner ring mass.

In this work we consider the potential of cometary impacts to deliver complex organic molecules and the prebiotic building blocks required for life to rocky exoplanets. Numerical experiments have demonstrated that for these molecules to survive, impacts at very low velocities are required. This work shows that for comets scattered from beyond the snow-line into the habitable zone, the minimum impact velocity is always lower for planets orbiting Solar-type stars than M-dwarfs. Using both an analytical model and numerical N-body simulations, we show that the lowest velocity impacts occur onto planets in tightly-packed planetary systems around high-mass (i.e. Solar-mass) stars, enabling the intact delivery of complex organic molecules. Impacts onto planets around low-mass stars are found to be very sensitive to the planetary architecture, with the survival of complex prebiotic molecules potentially impossible in loosely-packed systems. Rocky planets around M-dwarfs also suffer significantly more high velocity impacts, potentially posing unique challenges for life on these planets. In the scenario that cometary delivery is important for the origins of life, this study predicts the presence of biosignatures will be correlated with i) decreasing planetary mass (i.e. escape velocity), ii) increasing stellar-mass, and iii) decreasing planetary separation (i.e. exoplanets in tightly-packed systems).

The Near-Earth Object (NEO) Surveyor mission is a NASA observatory designed to discover and characterize near-Earth asteroids and comets. The mission's primary objective is to find the majority of objects large enough to cause severe regional impact damage ($>$140 m in effective spherical diameter) within its five-year baseline survey. Operating at the Sun-Earth L1 Lagrange point, the mission will survey to within 45 degrees of the Sun in an effort to find the objects in the most Earth-like orbits. The survey cadence is optimized to provide observational arcs long enough to reliably distinguish near-Earth objects from more distant small bodies that cannot pose an impact hazard. Over the course of its survey, NEO Surveyor will discover $\sim$200,000 - 300,000 new NEOs down to sizes as small as $\sim$10 m and thousands of comets, significantly improving our understanding of the probability of an Earth impact over the next century.

The Near Earth Object Surveyor mission has a requirement to find two-thirds of the potentially hazardous asteroids larger than 140 meters in size. In order to determine the mission's expected progress toward this goal during design and testing, as well as the actual progress during the survey, a simulation tool has been developed to act as a consistent and quantifiable yardstick. We test that the survey simulation software is correctly predicting on-sky positions and thermal infrared fluxes by using it to reproduce the published measurements of asteroids from the NEOWISE mission. We then extended this work to find previously unreported detections of known near Earth asteroids in the NEOWISE data archive, a search that resulted in 21,661 recovery detections, including 1,166 objects that had no previously reported NEOWISE observations. These efforts demonstrate the reliability of the NEOS Survey Simulator tool, and the perennial value of searchable image and source catalog archives for extending our knowledge of the small bodies of the Solar System.

Recent attempts to fully resolve the Hubble tension from early dark energy models seem to favor a primordial Harrison-Zeldovich Universe with its scalar spectrum being extremely scale invariant. Restoring the Harrison-Zeldovich spectrum within the single-field inflationary paradigm appears to be infeasible, turning to the multi-field approach from either curvaton or waterfall models. In this Letter, we successfully align with the Harrison-Zeldovich spectrum within a single-field chaotic inflation by a non-minimal derivative coupling, and the previously disfavoured chaotic potential by Planck+BICEP/Keck data in the standard $\Lambda$CDM model now returns back to the scope of future polarization observations of the cosmic microwave background.

We compare the spatial clustering and physical properties of the active galactic nuclei (AGN) and star-forming galaxies (SFG) at fixed stellar mass using a volume limited sample from the Sloan Digital Sky Survey (SDSS). The analysis of the two-point correlation function shows that the AGN are more strongly clustered than the SFG. The closer proximity to the $5^{th}$ nearest neighbour for the AGN than that for the SFG indicates that AGN prefer the denser regions. We compare the distributions of the $(u-r)$ colour, star formation rate (SFR), D$4000$ and concentration index of the AGN and SFG after matching their stellar mass distributions. The null hypothesis can be rejected at $>99.99\%$ confidence level in each case. The comparisons are also carried out at different densities. The differences persist at the same significance level in both the high and low density regions, implying that such differences do not originate from the variations in the density. Alternatively stated, the AGN activity can be triggered in both the high and low density regions. An analysis of the correlations between the different physical properties at fixed stellar mass reveals that the anticorrelations of SFR with morphology, colour and recent star formation history are $2-3$ times stronger for the AGN than for the SFG. It suggests that the presence of a bulge and the availability of gas are the two most crucial requirements for AGN activity. We propose a picture where the galaxies at fixed stellar mass may have widely different assembly histories, leading to significant variations in bulge properties and cold gas content. Whether a galaxy of a given stellar mass can acquire the suitable conditions for AGN activity remains uncertain due to a broad diversity of assembly history. We conclude that AGN are stochastic phenomena owing to an underlying role of the assembly bias.

We introduce a non-perturbative method to constrain the amplitude of local-type primordial non-Gaussianity ($f_{\rm NL}$) using squeezed configurations of the CMB lensing convergence and cosmic shear bispectra. First, we use cosmological consistency relations to derive a model for the squeezed limit of angular auto- and cross-bispectra of lensing convergence fields in the presence of $f_{\rm NL}$. Using this model, we perform a Fisher forecast with specifications expected for upcoming CMB lensing measurements from the Simons Observatory and CMB-S4, as well as cosmic shear measurements from a Rubin LSST/Euclid-like experiment. Assuming a minimum multipole $\ell_{\rm min}=10$ and maximum multipole $\ell_{\rm max}=1400$, we forecast $\sigma_{f_{\rm NL}}=175$ ($95$) for Simons Observatory (CMB-S4). Our forecasts improve considerably for an LSST/Euclid-like cosmic shear experiment with three tomographic bins and $\ell_{\rm min}=10$ and $\ell_{\rm max}=1400$ ($5000$) with $\sigma_{f_{\rm NL}}=31$ ($16$). A joint analysis of CMB-S4 lensing and LSST/Euclid-like shear yields little gain over the shear-only forecasts; however, we show that a joint analysis could be useful if the CMB lensing convergence can be reliably reconstructed at larger angular scales than the shear field. The method presented in this work is a novel and robust technique to constrain local primordial non-Gaussianity from upcoming large-scale structure surveys that is completely independent of the galaxy field (and therefore any nuisance parameters such as $b_\phi$), thus complementing existing techniques to constrain $f_{\rm NL}$ using the scale-dependent halo bias.

Systemic racism is a ubiquitous theme in societies worldwide and plays a central role in shaping our economic, social, and academic institutions. The Vera C. Rubin Observatory is a major US ground-based facility based in Chile with international participation. The Observatory is an example of excellence and will deliver the largest survey of the sky ever attempted. Rubin's full scientific and social potential can not be attained without addressing systemic racism and associated barriers to equity, diversity, and inclusion (EDI). During Rubin's 2021 virtual Project and Community Workshop (PCW), the annual Rubin community-based meeting, an anti-Black racism workshop took place, facilitated by 'The BIPOC Project' organization. About 60 members from different parts of the Rubin ecosystem participated. We describe the motivation, organization, challenges, outcomes, and near- and long-term goals of this workshop.

Previous searches for a gravitational-wave signal from a possible neutron star remnant of the binary neutron star merger event GW170817 have focused on short ($<500$ s) and long duration (2.5 hr -- 8 day) signals. To date, no such post-merger signal has been detected. We introduce a new piecewise model which has the flexibility to accurately follow gravitational-wave signals which are rapidly evolving in frequency, such as those which may be emitted from young neutron stars born from binary neutron star mergers or supernovae. We investigate the sensitivity and computational cost of this piecewise model when used in a fully coherent 1800-second $\mathcal{F}$-statistic search on simulated data containing possible signals from the GW170817 remnant. The sensitivity of the search using the piecewise model is determined using simulated data, with noise consistent with the LIGO second observing run. Across a 100--2000 Hz frequency band, the model achieves a peak sensitivity of $h_{\text{rss}}^{50\%} = 4.4 \times 10^{-23} \text{Hz}^{-1/2}$ at 200 Hz, competitive with other methods. The computational cost of conducting the search, over a bank of $1.1 \times 10^{12}$ templates, is estimated at 10 days running on 100 CPU's.

A lot of Dark Matter (DM) models predict single-scattering Nuclear Recoil (NR) events as potential signals. Characterization of a detector's NR response is crucial for underground DM direct detection projects. The traditional calibrating methods characterize a detector with $\sim$ keV - MeV neutrons emitted from an accelerator, a neutron generator, or a radioactive source. These methods could be improved in several aspects: (a) the neutrons's energy could be more mono-energetic, (b) the to-be-calibrated NR energy should line up with the ROI (Research Of Interest) of the experiment, (c) the calibrating beam's flux should be appropriate. Accordingly, in the paper, we introduce a novel calibration method in which the beam of the target material (instead of neutrons) will be implemented to impinge the detector. For instance, calibrating a liquid helium detector with a helium beam. Since the helium beam can (i) be tuned precisely to have $\sim < $1\% energy jitter (so as the recoil energy's uncertainty), (ii) have an energy interval from $\sim$ 100 eV to tens of keV (depending on the configuration), and (iii) have a tunable flux down to 10$^{10}$/s. Since these unique advantages can also be applied to other noble gases; therefore, we believe the method can be implemented to calibrate other liquid noble gas detectors' NR response.

We extend the classical magnetohydrodynamics formalism to include nonlocal quantum behavior via the phenomenological Bohm potential. We then solve the quantum magnetohydrodynamics equations to obtain a new analytical form of the dynamic structure factor (DSF), a fundamental quantity linking theory and experiments. Our results show that the three-peak structure -- one central Rayleigh peak and two Brillouin peaks -- of the DSF arising from quantum hydrodynamic fluctuations becomes (in general) a five-peak structure -- one central Rayleigh peak and two pairs of peaks associated with fast and slow magnetosonic waves. The Bohm contribution influences the positions and characteristics (height, width, and intensity) of the peaks by introducing three significant modifications: (a) an increase in effective thermal pressure, (b) a reduction in the adiabatic index, and (c) an enhancement of effective thermal diffusivity. The multiple DSF peaks enable concurrent measurements of diverse plasma properties, transport coefficients, and thermodynamic parameters in magnetized dense plasmas. The potential for experimental validation of our theory looms large, particularly through future experiments conducted at state-of-the-art laser facilities.

Modifications inspired by quantum gravity in the kinematics of special relativity can manifest in various ways, including anomalies in the time of flight of massless particles and the emergence of decay channels for otherwise stable particles. Typically, these effects are studied independently; however, it may be necessary to combine both to perform a consistent analysis. In this work, we study the interplay between time-of-flight anomalies and neutrino instability in the context of a flavor-independent high-energy Lorentz-invariance violation (LIV) in the neutrino sector. Ensuring compatibility between both types of effects imposes strong constraints on the existence of early neutrinos with energies exceeding a maximum value determined by the scale of new physics. Such constraints depend on the specific LIV scenario and should be integrated into searches for high-energy neutrinos from gamma-ray bursts exhibiting LIV time shifts.

The backreaction of quantum fields in their vacuum state results in equilibrium structures that surpass the Buchdahl compactness limit. Such backreaction is encapsulated in the vacuum expectation value of the renormalized stress-energy tensor (RSET). In previous works we presented analytic approximations to the RSET, obtained by dimensional reduction, available in spherical symmetry, and showed that the backreaction-generated solutions described ultracompact fluid spheres with a negative mass interior. Here, we derive a novel approximation to the RSET that does not rely on dimensional reduction, but rather on a perturbative reduction of the differential order. This approximation also leads to regular stars surpassing the Buchdahl limit. We conclude that this is a consequence of the negative energies associated with the Boulware vacuum which, for sufficiently compact fluid spheres, make the Misner-Sharp mass negative near the centre of spherical symmetry. Our analysis provides further cumulative evidence that quantum vacuum polarization is capable of producing new forms of stellar equilibrium with robust properties accross different analytical approximations to the RSET.