Papers for Thursday, May 29 2025

Next generation photometric and spectroscopic surveys will enable unprecedented tests of the concordance cosmological model and of galaxy formation and evolution. Fully exploiting their potential requires a precise understanding of the selection effects on galaxies and biases on measurements of their properties, required, above all, for accurate estimates of redshift distributions n(z). Forward-modelling offers a powerful framework to simultaneously recover galaxy n(z)s and characterise the observed galaxy population. We present GalSBI-SPS, a new SPS-based galaxy population model that generates realistic galaxy catalogues, which we use to forward-model HSC data in the COSMOS field. GalSBI-SPS samples galaxy physical properties, computes magnitudes with ProSpect, and simulates HSC images in the COSMOS field with UFig. We measure photometric properties consistently in real data and simulations. We compare redshift distributions, photometric and physical properties to observations and to GalSBI. GalSBI-SPS reproduces the observed grizy magnitude, colour, and size distributions down to i<23. Median differences in magnitudes and colours remain below 0.14 mag, with the model covering the full colour space spanned by HSC. Galaxy sizes are offset by 0.2 arcsec on average and some tension exists in the g-r colour, but the latter is comparable to that seen in GalSBI. Redshift distributions show a mild positive offset (0.01-0.08) in the mean. GalSBI-SPS qualitatively reproduces the stellar mass-SFR and size-stellar mass relations seen in COSMOS2020. GalSBI-SPS provides a realistic, survey-independent galaxy population description at a Stage-III depth using only parameters from the literature. Its predictive power will improve significantly when constrained against observed data using SBI, providing accurate redshift distributions satisfying the stringent requirements set by Stage IV surveys.

In the course of studying the 3C 220.3 lensing system, spectra were obtained with the Binospec instrument on the MMT for 511 additional objects in 3C 220.3's vicinity. These gave 146 good-quality galaxy redshifts and identified 126 Galactic stars. The galaxy redshift histogram shows a peak near 3C 220.3's redshift, but there is no evidence for or against a galaxy group within 2 Mpc of 3C 220.3 itself. The spectra revealed 12 AGN candidates including a likely $z\approx4.64$ broad-line QSO. Visible and near-infrared imaging with HST allowed morphological classifications of 14 galaxies. One system is a potential analog of the Milky Way-LMC system with stellar mass ratio $\sim$0.6.

The measurement of the temporal evolution in the redshift of distant objects, the redshift drift, is a probe of universal expansion and cosmology. We perform the first steps towards a measurement of such effect using the Lyman-$\alpha$ forest in the spectra of bright quasars as a tracer of cosmological expansion. Our goal is to determine to which precision a velocity shift measurement can be carried out with the signal-to-noise (S/N) level currently available and whether this precision aligns with previous theoretical expectations. A precise assessment of the achievable measurement precision is fundamental for estimating the time required to carry out the whole project. We acquire 12 hours of ESPRESSO observations distributed over 0.875 years of the brightest quasar known, J052915.80-435152.0 (z=3.962), to obtain high-resolution spectra of the Lyman-$\alpha$ forest, with median S/N of ~86 per 1 km/s pixel at the continuum. We divide the observations into two epochs and analyse them using both a pixel-by-pixel method and a model-based approach. This comparison allows us to estimate the velocity shift between the epochs, as well as the velocity precision that can be achieved at this S/N. The model-based method is calibrated using high-resolution simulations of the intergalactic medium, and it provides greater accuracy compared to the pixel-by-pixel approach. We measure a velocity drift of the Lyman-$\alpha$ forest consistent with zero: $\Delta v = -1.25\pm 4.45 {\rm ms^{-1}}$, equivalent to a cosmological drift of $\dot{v}=-1.43\pm 5.09 {\rm ms^{-1}yr^{-1}}$ or $\dot{z}= (-2.19\pm7.77) \times 10^{-8}{\rm yr^{-1}}$. The measurement uncertainties are on par with the expected precision. We estimate that reaching a 99% detection of the cosmic drift requires a monitoring campaign of 5400 hours of integration time over 54 years with an ELT and an ANDES-like high-resolution spectrograph.

We present a parameter survey of fragmentation in collapsar disks, using a revised version of the Chen & Beloborodov (2007) model that determines the structure of steady state hyperaccretion disks in a general relativistic and neutrino cooled framework. We map out the range of disk conditions leading to gravitational instability alongside an exploration of the dimensionless cooling time $\beta$, which together determine whether fragmentation is likely to occur. We estimate the initial mass and density of fragments, finding that they occupy a unique region in the space of self-gravitating compact objects, with masses $M_{\rm f} \sim 10^{-3} M_\odot -10^{-1} M_\odot$ and densities $\rho_{\rm f}\sim 10^8-10^{11}~{\rm g~cm}^{-3}$. We then calculate their migration and mass growth (via Bondi-Hoyle accretion) as they inspiral through the collapsar disk. During a fragment's migration to the central black hole, it can grow its mass up to a range $M_{\rm f}\sim 10^{-1} M_\odot - 1 M_\odot$. In most cases, the final fragment mass is larger than the minimum cold stable neutron star mass but much smaller than any observed neutron star. The fragment briefly achieves peak accretion rates comparable to (or even larger than) that of the central engine. We propose that these bound fragments may give rise to observable astrophysical phenomena, and we approximately model two of these: (i) gamma ray burst variability produced by a secondary, fragment-launched jet; (ii) the generation of non-vacuum gravitational waveforms accompanied by electromagnetic counterparts.

To better understand how galaxies assemble their structure and evolve over cosmic time, we present a new catalog of morphological measurements for over 340,000 sources spanning $0 < z < 12$, derived from deep JWST NIRCam imaging across four major extragalactic fields (CEERS, PRIMER-UDS, PRIMER-COSMOS, GOODS) compiled in the DAWN JWST Archive (DJA). We perform two-dimensional surface brightness fitting for all galaxies in a uniform, flux-limited sample. Each galaxy is modeled with both a Sérsic profile and a two-component (bulge and disk) decomposition, yielding consistent structural parameters - including effective radius, Sérsic index, axis ratio, and bulge-to-total ratio ($B/T$). To demonstrate the scientific application of our morphology catalogs, we combined these measurements with DJA photometric redshifts, physical parameters and rest-frame colors, and investigated the relation between total, bulge and disk sizes, Sérsic index ($n_S$), star formation activity, and redshift. Bulge-dominated galaxies (high $n_S$ and $B/T$) predominantly occupy the quiescent region of the $UVJ$ diagram, while disk-dominated galaxies are mostly star-forming. A significant bimodality persists, with quiescent disks and compact, bulge-dominated star-forming galaxies observed out to $z > 3$. Quiescent galaxies also show significantly higher stellar mass surface densities, nearly an order of magnitude greater at $z \sim 4$ than at $z \sim 1$. Our results confirm a strong and evolving link between morphology and star formation activity, and support a scenario in which bulge growth and quenching are closely connected. This work is a highly valuable addition to the DJA, adding a morphological dimension to this rich dataset and thus enabling a wider scientific application.

The baryon fraction of galaxy clusters is a powerful tool to inform on the cosmological parameters while the hot-gas fraction provides indications on the physics of the intracluster plasma and its interplay with the processes driving galaxy formation. Using cosmological hydrodynamical simulations from The Three Hundred collaboration of about 300 simulated massive galaxy clusters with median mass $M_{500}\approx7 \times 10^{14}$M$_{\odot}$ at $z=0$, we model the relations between total mass and either baryon fraction or the hot gas fractions at overdensities $\Delta = 2500$, $500$, and $200$ with respect to the cosmic critical density, and their evolution from $z\sim 0$ to $z\sim 1.3$. We fit the simulation results for such scaling relations against three analytic forms (linear, quadratic, and logarithmic in a logarithmic plane) and three forms for the redshift dependence, considering as a variable both the inverse of cosmic scale factor, $(1+z)$, and the Hubble expansion rate, $E(z)$. We show that power-law dependencies on cluster mass poorly describe the investigated relations. A power-law fails to simultaneously capture the flattening of the total baryon and gas fractions at high masses, their drop at the low masses, and the transition between these two regimes. The other two functional forms provide a more accurate description of the curvature in mass scaling. The fractions measured within smaller radii exhibit a stronger evolution than those measured within larger radii. From the analysis of these simulations, we conclude that as long as we include systems in the mass range herein investigated, the baryon or gas fraction can be accurately related to the total mass through either a parabola or a logarithm in the logarithmic plane. The trends are common to all modern hydro simulations, although the amplitude of the drop at low masses might differ [Abridged].

The young $\beta$ Pictoris system has been monitored with high-resolution optical spectrographs for decades. These observations have revealed strongly variable absorption in the Ca II H\&K lines attributed to in-falling cometary bodies. Since 2003, over 9000 HARPS observations of $\beta$ Pictoris have been taken and many of these have not yet been used for exocomet studies. We search these spectra for new exocomet phenomenology enabled by the long time coverage and large volume of this dataset. We systematically carry out telluric correction of the HARPS spectra using molecfit, compare multi-year observations of the Ca II and Na I lines and use a Bayesian fitting algorithm to extract exocomet line parameters. We explore the usage of an unbiased reference spectrum with which to calibrate the continuum, and investigate Keplerian orbital solutions to observed exocomet acceleration. We find a general absence of exocometary sodium line absorption, with only two instances of clear ($\sim 2$ % deep) exocometary sodium out of 198 nights of observation, as well as a weaker ($\sim 1$ %) feature that persists over 13 nights in 2004. We find that these events occur during times of deep Ca II absorption at the same red-shift, implying that strongly Ca II-evaporating exocomets also exhibit detectable levels of Na I, in spite of the vast majority of Na I being rapidly photo-ionised in close proximity to the star. We find long-lived Ca II absorption in 2017 and 2018 that persists on a timescale of a year, which may be difficult to explain with the classical exocomet model. Finally, we investigate two strongly accelerating, blue-shifted exocomes observed in 2019 that show strong and sudden departures from Keplerian motion, suggesting rapid changes to the dynamics of the exocomet cloud. We hypothesize that this is caused by the destruction of the comet nuclei shortly after their periastron passages.

We investigate galaxy groups that reside in the field but have been previously processed by galaxy clusters. Observationally, they would appear to have the same properties as regular field groups at first glance. However, one would expect to find quantifiable differences in processed groups as dynamical interactions within clusters perturb them. We use IllustrisTNG300 simulation to statistically quantify that processed groups of galaxies show different properties compared to regular field groups. Our analysis encompasses a broad range of groups with total masses between $8 \times 10^{11} \Msun$ and $7 \times 10^{13} \Msun$. We distinguish between processed groups that passed through a galaxy cluster and capture more galaxies, referred to as thief groups, and groups that did not capture any new members, referred to as non-thief groups. The employed statistical tools show that thief groups are generally less compact and contain more members, while non-thief groups seem to have the same properties as the field groups which makes them indistinguishable.

Absorption outflows in quasars play an important role in understanding active galactic nuclei (AGN) feedback and their influence on galaxy evolution. The unprecedented spectral data provided by the Dark Energy Spectroscopic Instrument (DESI) opens new avenues to explore these outflows. We analyze five low-ionization absorption outflow systems in four intermediate-redshift quasars (2 < z < 3) using the data obtained by DESI in order to characterize their physical properties and energetics, and also to assess their role in AGN feedback. We use the spectra from DESI's Early Data Release to determine the ionic column densities, total hydrogen column densities, electron number densities, and ionization parameters via photoionization modeling and absorption line analysis. We derive the outflows' distance from the AGN, and its kinematic properties: mass-flow rates, kinetic luminosity, and momentum flux. Our study identifies five distinct mini-broad absorption line outflow systems, hosted by four quasars. The identified outflows exhibit hydrogen column densities of \log (N_H) = 20.0-20.7[cm^{-2}], ionization parameters of log (U_H) = -2.7 to -2.1, and electron number densities of log (n_e) =1.45-2.85[cm^{-2}]. The distances of the outflows from the central source are between 4.5 to 31 kpc, and the kinetic luminosities range from 2\times10^{-6} to 5\times10^{-3} of the Eddington luminosity. The outflows in J1407 and J1032 show the largest kinetic luminosities, with momentum flux ratios (\dot{p} / \dot{p}_{rad}) of about 2 and 0.2, respectively. Our findings highlight the vital role of DESI data in uncovering the diversity and significance of quasar outflows in galaxy evolution.

Geminga is the third gamma-ray pulsar firmly detected by imaging atmospheric Cherenkov telescopes (IACTs) after the Crab and the Vela pulsars. Most of its emission is expected at tens of GeV, and, out of the planned telescopes of the upcoming Cherenkov Telescope Array Observatory (CTAO), the Large-Sized Telescopes (LSTs) are the only ones with optimised sensitivity at these energies. We aim to characterise the gamma-ray pulse shape and spectrum of Geminga as observed by the first LST (hereafter LST-1) of the CTAO-North. Furthermore, this study confirms the great performance and the improved energy threshold of the telescope, as low as 10 GeV for pulsar analysis, with respect to current-generation Cherenkov telescopes. We analysed 60 hours of good-quality data taken by the LST-1 at zenith angles below 50$^\circ$. Additionally, a new Fermi-LAT analysis of 16.6 years of data was carried out to extend the spectral analysis down to 100 MeV. Lastly, a detailed study of the systematic effects was performed. We report the detection of Geminga in the energy range between 20 and 65 GeV. Of the two peaks of the phaseogram, the second one, P2, is detected with a significance of 12.2$\sigma$, while the first (P1) reaches a significance level of 2.6$\sigma$. The best-fit model for the spectrum of P2 was found to be a power law with $\Gamma = (4.5 \pm 0.4_{stat})^{+0.2_{sys}}_{-0.6_{sys}}$, compatible with the previous results obtained by the MAGIC. No evidence of curvature is found in the LST-1 energy range. The joint fit with Fermi data confirms a preference for a sub-exponential cut-off over a pure exponential, even though both models fail to reproduce the data above several tens of GeV. The overall results presented in this paper prove that the LST-1 is an excellent telescope for the observation of pulsars, and improved sensitivity is expected to be achieved with the full CTAO-North.

Modified Newtonian Dynamics (MOND) is an alternative to the dark matter hypothesis that attempts to explain the "missing gravity" problem in astrophysics and cosmology through a modification to objects' dynamics. Since its conception in 1983, MOND has had a chequered history. Some phenomena difficult to understand in standard cosmology MOND explains remarkably well, most notably galaxies' radial dynamics encapsulated in the Radial Acceleration Relation. But for others it falls flat -- mass discrepancies in clusters are not fully accounted for, the Solar System imposes a constraint on the shape of the MOND modification seemingly incompatible with that from galaxies, and non-radial motions are poorly predicted. An experiment that promised to be decisive -- the wide binary test -- has produced mainly confusion. This article summarises the good, the bad and the ugly of MOND's observational existence. I argue that despite its imperfections it does possess ongoing relevance: there may yet be crucial insight to be gleaned from it.

Long gamma-ray bursts (GRBs) serve as powerful probes of distant galaxies. Their luminous afterglow pinpoints galaxies independent of luminosity, in contrast to most flux-limited surveys. Nevertheless, GRB-selected galaxy samples are not free from bias, instead tracing the conditions favoured by the progenitor stars. Characterising the galaxy populations traced by GRBs is therefore important both to effectively use GRBs as probes as well as to place stronger constraints on the progenitor stars capable of forming long GRBs. Spatially-resolved spectroscopic observations with integral field units (IFUs) provide valuable insights into the interstellar medium and stellar populations of GRB host galaxies. In this paper we present results of the first two GRB host galaxies observed with the JWST/NIRSpec IFU with a spatial resolution of ~ 1.6 kpc; the hosts of GRB 150403A and GRB 050820A at redshifts z ~ 2.06 and z ~ 2.61, respectively. The data reveal two complex galaxy environments made up of two or more star forming galaxies that are likely interacting given their small spatial separation (< 20 kpc) and line of sight velocity offsets (< 100 km/s). The measured gas-phase metallicity, star formation rates (SFRs), and key diagnostic line ratios for each of the detected galaxies are overall consistent with the properties of other star forming galaxies and GRB hosts at z > 2. However, differences in the SFR and metallicities of the interacting galaxies highlight the importance of spatially resolved observations in order to accurately characterise the galaxy properties traced by GRBs.

Scaling relations between galaxy cluster properties are crucial for understanding cosmology and baryonic physics. Rigorous calibration of the $M-L_X$ relation, employing weak lensing mass and consistent statistical methodology, is challenging due to heterogeneous cluster samples. The release of LEGACY imaging data introduced the possibility of unifying the cluster selection. We present the all-sky extension of the CODEX catalog based on LEGACY data and introduce a Bayesian framework for calibrating the X-ray luminosity-mass relation, derived for 100 clusters with weak lensing mass measurements. Using the X-ray luminosity estimates for those clusters from ROSAT All-Sky Survey (RASS) data, we perform a power-law fit to the $M-L_X$ relation. Furthermore, taking advantage of the recently released eROSITA data (eRASS1), we assess the impact of point source contamination on cluster fluxes for 42 clusters in the eRASS1 footprint. The RASS fit yields a slope of $\beta = 0.75 \pm 0.09$, 1.7$\sigma$ lower than the best self-similar prediction, with marginal evidence for the redshift evolution of the normalization ($\gamma = 0.65 \pm 0.43$). As for the eRASS1 analysis, the slope is substantially steeper, $\beta = 1.11 \pm 0.15$, and in further agreement with the prediction of self-similarity. No additional evolution is also seen ($\gamma = 0.004 \pm 0.790 $). While our results provide the practical means for cosmological studies of both RASS and eRASS data, the link to cluster physics is much cleaner after the cluster flux contamination is reduced. We also analyzed the impact of the selection function on calibration, finding that its full modeling is essential.

Post-merger gravitational-wave emission from a binary neutron star merger carries crucial information about the equation of state (EoS) of matter at high temperatures. Current gravitational wave detectors have limited sensitivities at post-merger frequencies in the range [1.5, 4] kHz. Therefore, valuable inferences can only be made after combining information from multiple (BNS) events. Criswell et al. [Phys. Rev. D 107, 043021 (2023)] carries out an injection study to infer the radius posterior for a $1.6 M_{\odot}$ NS by combining the information from injected BNS events via the hierarchical Bayesian inference (HBI) formulation. This formulation utilizes empirical relations that connect the peak frequency of the post-merger remnant with the chirp mass of the system, and the EoS proxy parameter $R_{1.6}$. In this work, we extend the HBI formulation to other EoS proxy parameters (i.e., the radius of (NS), $R_{1.X}$, with masses 1.2, 1.4, and 1.8 $M_{\odot}$) and combine the four $R_{1.X}$ posteriors through the piecewise polytropic EoS model to obtain measurable constraints on the EoS of the NS. We show that the NS radii can be constrained to within $\sim$1 km ( $\sim$0.55 km) assuming uniform (astrophysical) prior on $R_{1.X}$ for injections in the A+ noise. We also study systematic biases in the analysis coming from the limitations of empirical relations.

Satellite galaxies endure powerful environmental tidal forces that drive mass stripping of their outer regions. Consequently, satellites located in central regions of galaxy clusters or groups, where the tidal field is strongest, are expected to retain their central dense regions while losing their outskirts. This process produces a spatial segregation in the mean mass density with the cluster-centric distance (the $\bar{\rho}-r$ relation). To test this hypothesis, we combined semi-analytical satellite orbital models with cosmological galaxy simulations. We find that not only the mean total mass densities ($\bar{\rho}$), but also the mean stellar mass densities ($\bar{\rho}^{\star}$) of satellites exhibit this distance-dependent segregation ($\bar{\rho}^{\star}-r$). The correlation traces the host's tidal field out to a characteristic transition radius at $\Re_{\star}$ $\approx$ $0.5$ $R_{\rm vir}$, beyond which the satellite population's density profile can have a slight increase or remain flat, reflecting the weakened tidal influence in the outskirts of galaxy clusters and beyond. We compare these predictions with observational data from satellites in the Virgo and Fornax galaxy clusters, as well as the Andromeda and Milky Way systems. Consistent trends in the satellite mean stellar mass densities are observed across these environments. Furthermore, the transition radius serves as a photometric diagnostic tool: it identifies regions where the stellar components of satellites underwent significant tidal processing and probes the gravitational field strength of the host halo.

I present infrared spectroscopy of 37 brown dwarf candidates in the Upper Sco association, 35 of which are classified as young and cool, making them likely members. This sample includes many of the faintest spectroscopically confirmed members ($K=16$-17 mag), which should have masses down to $\sim0.007$-0.01 $M_\odot$ for the range of ages in Upper Sco (7-14 Myr). Using my updated membership catalog for Upper Sco, I have estimated the initial mass function (IMF) for a field in the center of the association that encompasses $\sim80$\% of the known members. I have derived IMFs in the same manner for previous membership samples in three other star-forming populations, consisting of IC 348, Taurus, and Chameleon I. When using logarithmic mass bins, the substellar IMFs for Upper Sco and the other young regions are roughly flat down to the completeness limits of $\sim$0.01 $M_\odot$. These IMFs are broadly similar to mass functions recently measured for the solar neighborhood. Finally, I have used W1$-$W2 colors to search for excess emission from circumstellar disks among the late-type objects in my new census of Upper Sco. I measure an excess fraction of 52/200 for members with spectral types of M6.25-M9.5, which is similar to results from previous membership catalogs. For the L-type members, it is difficult to detect the small W2 excess emission produced by typical disks around brown dwarfs because of the large uncertainties in spectral types, which preclude accurate estimates of the photospheric colors. Thus, W2 photometry provides poor constraints on the presence of disks for the L-type members of Upper Sco.

Aims. This study investigates the magnetic evolution of solar active regions (ARs), with a particular focus on understanding how the magnetic morphology of simple and complex ARs changes throughout their lifetime. Methods. To analyse the magnetic evolution of ARs, we developed a Magnetic Evolution Method (MEM) that segments each region's lifetime into three phases: growth, main, and recovery. The method was applied to ARs observed between January 1996 and December 2020. Results. We found that complex active regions (CARs) have a mean lifetime of approximately 24 days, which is 8 days longer than that of simple active regions (SARs). Most CARs (94%) first appear with a simple magnetic structure and remain in this configuration for about 3 days (growth phase), before transitioning into complex structures for around 5 days (main phase), after which they typically revert to a simple state (recovery phase). The average lifetimes of SARs and CARs show no significant difference between solar cycles 23 and 24, suggesting that active region lifetimes are independent of the solar cycle. Conclusions. By tracking the full magnetic evolution of ARs, our study reveals that CARs typically become magnetically complex around three days after emergence and remain in that state for a limited but critical period. This temporal structure, uncovered using a novel method that follows ARs throughout their full development, provides important context for identifying the magnetic conditions associated with increased eruptive potential. The results offer a foundation for improving the forecasting of solar flares and magnetic clouds, and suggest that the magnetic evolution of ARs is largely independent of the solar cycle.

The Sunspot Solar Observatory Data Archive (SSODA) stores data acquired with the suite of instruments at the Richard B. Dunn Solar Telescope (DST) from February 2018 to the present. The instrumentation at the DST continues to provide high cadence imaging, spectroscopy, and polarimetry of the solar photosphere and chromosphere across a wavelength range from 3500Å to 11,000Å. At time of writing, the archive contains approximately 374 TiB of data across more than 520 observing days (starting on February 1, 2018). These numbers are approximate as the DST remains operational, and is actively adding new data to the archive. The SSODA includes both raw and calibrated data. A subset of the archive contains the results of photospheric and chromospheric spectropolarimetric inversions using the Hazel-2.0 code to obtain maps of magnetic fields, temperatures, and velocity flows. The SSODA represents a unique resource for the investigation of plasma processes throughout the solar atmosphere, the origin of space weather events, and the properties of active regions throughout the rise of Solar Cycle 25.

Throughout its mission, the Mars Reconnaissance Orbiter (MRO) has often rolled about its along-track axis by up to 28° to partially compensate for the suboptimal location of the Shallow Radar (SHARAD) antenna along an edge of the spacecraft that is opposite the imaging payload deck, thereby enhancing the signal-to-noise ratio (S/N) of echoes returned from the surface. After recent modeling work predicted that a much larger roll would improve the S/N by ~10 dB relative to nadir-pointed observing, MRO began a limited series of 120° roll maneuvers to test the effects on radar sounding. Three such SHARAD very-large-roll (VLR) observations have been acquired since May 2023, and they show dramatic improvements in signal clarity and depth of penetration, with S/N increasing by 9, 11, and 14 dB over that of nearly coincident observations at 0° roll angle. In low dielectric terrains, the first and second VLR observations enabled basal detections at depths previously unachievable, reaching depths of 800 m in Medusae Fossae materials and 1500 m through the ice of Ultimi Scopuli, respectively. The second VLR observation also obtained enhanced reflections throughout the ice stack. In the higher dielectric terrain of Amazonis Planitia, the third VLR observation improved continuity of a dipping subsurface interface, but it revealed neither an extension of the interface to greater depths nor any deeper interfaces. The MRO mission intends to obtain more SHARAD VLR observations of polar terrains and of midlatitude glacial and ground ices, sediments, and volcanics.

Pluto's argument of perihelion is known to librate around $90^\circ$. This libration is related to the secular phenomenon known as the von Zeipel-Lidov-Kozai (vZLK) oscillation. In this work, we make a quantitative assessment of the influence of Neptune's mean motion resonance and of the other giant planets' secular perturbations on the libration of Pluto's argument of perihelion. Here, a parameter $k^2 = (1-e^2) \cos^2 I$ is the key where $e$ is eccentricity and $I$ is inclination. When $k^2$ of a Pluto-like object is larger than a certain critical value, libration of its argument of perihelion would not occur. The secular effect of other disturbing planets (Jupiter, Saturn, Uranus) plays a significant role in determining the critical $k^2$. The non-zero oscillation amplitude of the critical resonant argument also plays a role, although not a dominant one.

We extend the analytic expressions for polytropic wind-driven bubbles and their shock structures, formulated initially in Koo and McKee 1992(a,b), focusing on spherically symmetric configurations in astrophysical environments with $\rho\propto r^{-2}$, which arises naturally in the star-forming environment and has applications to winds flowing into a preexisting bubble. Wind luminosity is assumed to be constant, and as a result the shock velocities of these bubbles are constant in time. The ratio of specific heats is assumed to be the same in the shocked ambient medium and the shocked wind. Numerical results are presented for one selected ratio of wind density to ambient density. Exact ODEs are written for the compressed wind region and approximate solutions are found by fitting the ODE solutions. By analyzing the interactions between stellar winds and ambient media in the strong compression limit, we model the formation and evolution of spherical bubbles, highlighting their shock fronts and contact discontinuities. Our analytic method provides an intuitive approach to calculating the thickness of bubble shells, which is crucial for understanding their dynamics and observational characteristics. A numerical method explores conditions without explicitly requiring the strong compression limit, and then we compare numerical to analytical results under various conditions.

A group of poorly studied eclipsing variables (the classification of which is marked as uncertain and/or the period of brightness changes is uncertain) has been studied with the using of the photometric observations of the TESS mission and NSVS, ASAS-SN sky surveys. We also obtained some observations covering the brightness minima of our variables by our group using the telescopes at Astronomical Observatory on Kolonica Saddle (Slovakia) and Observatory and Planetarium in Hlohovec (Slovakia) during the "Variable-2024" astrocamp. The periods and classification were corrected. For NSV 575 and NSV 014 the periods were found for the first time, but it is doubtful that NSV 014 is an eclipsing variable, because there are no eclipses but the asymmetric wave is present, which indicates that the variable star can be re-classified as a low-amplitude pulsating one. Different methods were used for approximation of the light curves and further calculation of stellar system's parameters such as eclipse depths and durations, values of reflection effect and effect of ellipticity of stars. The initial period was estimated using the periodogram based on the trigonometrical polynomial fit of high order (up to 10). For better approximation of the complete eclipsing phase curve, the "New Algol Variable" (NAV) software was used. The methods of "asymptotic parabolas" and "wall-supported asymptotic parabolas" were used for calculation of moments of eclipses, which use only near-eclipse part of the observations instead of a complete curve. These methods were implemented in the software MAVKA among a larger set of features. For the variables NSV 489 and NSV 1884, our moments of eclipses and the ones found in the literature, were used for the O-C curves. For NSV 489, the period was adjusted taking into account the slope of the O-C diagram.

Ultra High Energy Cosmic Rays, UHECR, are charged particles with energies between $\sim10^{18}\,{\rm eV}$ and $\sim3\times10^{20}\,{\rm eV}\sim50\,{\rm J}$. They exhibit fundamental physics at energies inaccessible to terrestrial accelerators, challenge experimental physics and connect strongly to astronomical observations through electromagnetic, neutrino and even gravitational wave channels. There has been much theoretical and observational progress in the sixty years that have elapsed since the discovery of UHECR, to divine their nature and identify their sources. The highest energy UHECR appear to be heavy nuclei with rigidity extending up to $\sim10\,{\rm EV}$; A significant ($6.9\sigma$) dipole anisotropy has been measured but our poor understanding of the Galactic magnetic fields makes this hard to interpret; The UHECR luminosity density is $\sim 10^{44}$ erg Mpc$^{-3}$ yr$^{-1}$ which constrains explanations of their origin; The most promising acceleration mechanisms involve diffusive shock acceleration and unipolar induction; The most promising sources include intergalactic accretion shocks, and relativistic jets from stellar-mass or supermassive black holes. We explore the prospects for using the highest energy events, combined with multimessenger astronomy, to help us solve the riddle of UHECR.

Context: We present results of time-resolved optical spectroscopy and photometry of the short-orbital period dwarf nova EI Psc. Aims: This study aims to determine fundamental system parameters of EI Psc, study properties of accretion structures in the system, and investigate its origin and current evolution state. Methods: We analyse newly obtained time-resolved spectroscopic and photometric observations as well as archival data. We used light curve modelling, Doppler tomography, and MESA evolutionary models to study the characteristics of EI Psc. Results: The system contains a relatively low temperature ($T_{\rm eff} = 6130\,\mathrm{K}$) white dwarf with mass of $M_{\mathrm WD} = 0.70(4)\,\mathrm{M}_{\odot}$. The mass of the warm ($T_2 = 4440\,\mathrm{K}$) secondary is $M_2 = 0.13\,\mathrm{M}_{\odot}$. The inclination of the system is $i= 44.5°(7)$. The mass accretion rate is $\approx$ $4\times10^{-13}\,\mathrm{M}_\odot\,\mathrm{year}^{-1}$. The long-term light curve of the system shows outbursts and superoutbursts. The quiescence light curve is double-humped and is formed by the combination of radiation from the Roche lobe filling the hot secondary and the hot spot. The radius of the outer disc is about two times smaller than the tidal truncation radius. Most of the disc's emission consists of emission lines and radiation from the hot spot at the stream-disc impact region. Conclusions: These types of systems are formed from progenitors with a low mass WD $\mathrm{M}_{\mathrm{WD}} \lesssim 0.6\,\mathrm{M}_\odot$ and relatively massive secondaries $1.1-1.5\,\mathrm{M}_\odot$ with initial orbital periods on a scale of days. The number of similar systems is expected to be significantly lower than the usual CVs due to a lower forming rate of their relatively massive progenitors.

This study investigates the projected, quasi-symmetric $\sim\rm46\,kpc$-scale diffuse radio lobes surrounding the giant elliptical galaxy M\,87, utilizing well-sampled wideband ($\rm 60\,MHz-10.55\,GHz$) observations from MWA and VLA, supplemented by data from LOFAR and Effelsberg. The observed structures feature sharp edges and filaments, with nearly uniform and moderately steep spectral indices ($\alpha$, mostly within $-1.2\leq\alpha\leq-0.8$), indicating turbulence. Well-sampled radio spectra for the lobes' diffuse region are derived using the continuous injection (CI) model (with $\alpha_{\rm inj}\simeq-0.86$ and $\nu_{\rm b}\simeq1.72\rm\,GHz$), and for its three localized regions using the impulsive injection model (e.g., JP model). From energy equipartition analysis, we estimate the typical magnetic field strength in the lobes' diffuse region to be $B_{\rm eq}\simeq10\,\mu\rm G$. The age of the lobes is estimated as $\sim30-50\,\rm~Myr$, based on lifetimes derived from the CI and JP models and sound crossing time. Outflow powers of $\sim(0.2-2)\times10^{44}\,\rm erg\,s^{-1}$ for the lobes' diffuse components and $\sim(1-11)\times10^{44}\,\rm erg\,s^{-1}$ for the whole source are calculated. With this power assessment, we conclude that the galactic stellar wind has a negligible effect, the active galactic nucleus (AGN)-driven jet can provide the necessary energy for the whole system. Furthermore, we argue that while the wind driven by current AGN activity is unlikely to power the lobes' diffuse components, an average enhancement of AGN activity by a factor of $\sim 10^2$ over the past $\sim 30-50$ Myr remains plausible.

Since the emission of water molecules cannot be observed from Earth, less abundant isotopologues, such as H$_2^{18}$O and HDO, are used to trace water in star-forming regions. The main aim of this study is to determine HDO abundance in the hot core RCW 120 S2. We performed observations of the hot core in the 200-255~GHz range using the nFLASH230 receiver on the APEX telescope. Two HDO lines were detected toward RCW 120 S2. Their intensities are described by excitation temperature $286\pm2$~K and gas number density $\geq 10^9$~cm$^{-3}$. The emission originates from the hot core rather than the warm dense envelope surrounding a central young stellar object. The HDO column density ranges from $(3.9-7.9)\times10^{13}$~cm$^{-3}$ with a most probable value $5.6\times10^{13}$~cm$^{-3}$. The HDO abundance relative to hydrogen nuclei is $1.5\times10^{-9}$. This HDO abundance value is among the lowest reported for hot cores. Combined with the non-detection of the H$_2^{18}$O line, we conclude that protostellar heating in RCW 120 S2 is still in its early stages.

In a recent paper we parameterised the evolution of the star-formation rate dispersion ($\sigma_{SFR}$) across the specific star-formation rate - stellar mass plane (sSFR-M$_{\star}$) using the Deep Extragalactic VIsible Legacy Survey (DEVILS) - suggesting that the point at which the minimum in the dispersion occurs (M$^{*}_{\sigma-min}$) defines a boundary between different physical mechanisms affecting galaxy evolution. Here we expand upon that work to determine the movement of galaxies through the sSFR-M$_{\star}$ plane using their recent star-formation histories (SFHs) and explore how this leads to the observed $\sigma_{SFR}$-M$_{\star}$ relation. We find that galaxies in sub-regions of the sSFR-M$_{\star}$ plane show distinctly different SFHs, leading to a complex evolution of the sSFR-M$_{\star}$ plane and star-forming sequence (SFS). However, we find that selecting galaxies based on stellar mass and position relative to SFS alone (as is traditionally the case), may not identify sources with common recent SFHs, and therefore propose a new selection methodology. We then use the recent SFH of galaxies to measure the evolution of the SFS, showing that it has varying contributions from galaxies with different SFHs that lead to the observed changes in slope, normalisation and turnover stellar mass. Finally, we determine the overall evolution of the sSFR-M$_{\star}$ plane from $z\sim1$ to today. In the second paper in this series we will discuss physical properties of galaxies with common recent SFHs and how these lead to the observed $\sigma_{SFR}$-M$_{\star}$ relation and evolution of the sSFR-M$_{\star}$ plane.

In part I of this series we discussed the variation of star-formation histories (SFHs) across the specific star formation rate - stellar mass plane (sSFR-M$_{\star}$) using the Deep Extragalactic VIsible Legacy Survey (DEVILS). Here we explore the physical mechanisms that are likely driving these observational trends, by comparing the properties of galaxies with common recent SFH shapes. Overall, we find that the processes shaping the movement of galaxies through the sSFR-M$_{\star}$ plane can be be largely split into two stellar mass regimes, bounded by the minimum SFR dispersion ($\sigma_{SFR}$) point. At lower stellar masses we find that large $\sigma_{SFR}$ values are likely observed due to a combination of stochastic star-formation processes and a large variety in absolute sSFR values, but relatively constant/flat SFHs. While at higher stellar masses we see strong observational evidence that Active Galactic Nuclei (AGN) are associated with rapidly declining SFHs, and that these galaxies reside in the high $\sigma_{SFR}$ region of the plane. As such, we suggest that AGN feedback, leading to galaxy quenching, is the primary driver of the high $\sigma_{SFR}$ values. These results are consistent with previous theoretical interpretations of the $\sigma_{SFR}$-M$_{\star}$ relation.

We conducted a comprehensive spectral and timing analysis of NGC 6814 using AstroSat's 2019 and XMM-Newton's 2021 observations. Cross-correlation analysis revealed a significant correlation between FUV (1541 Å)/X-ray and UVW1 (2910 Å)/X-ray variations, with delays of $\sim 15~\rm{ks}$ and $30~\rm{ks}$, respectively. We constructed four broadband SEDs after applying aperture correction (for the UVIT filter), subtracting host galaxy and emission line contributions from UV flux, and using mean X-ray spectra alongside selected UV data points. First, we fitted the SEDs with KYNSED model assuming various combinations of inclination, $\theta$, color correction factors, $f_{\rm col}$, and BH spins. Best-fit models were achieved for $\theta=70^{\circ}$ (consistent with past estimates for this source) and for spin $\leq 0.5$, while $f_{\rm col}$ is not constrained. KYNSED provided satisfactory fit to all SEDs in the case when the corona is powered by the accretion process, with $\sim 10-20$% of the accretion power transferred to the corona, $\dot{m}/\dot{m}_{\rm Edd}\sim 0.1$, corona radius of $\sim 6-10~r_g$, and height of $\sim7.5-35~r_g$. Model time-lags computed using the SED best-fit results are aligned well with the observed time-lags. Although some of the model parameters are not constrained, the important result of our work is that both the broadband X-ray/UV spectra and the X-ray/UV time-lags in NGC 6814 are consistent with the hypothesis of X-ray illumination of the disc in a lamp-post geometry framework. Within this model framework, we do not need to assume an outer or inner truncated disc.

Whether the presence of planets affects the lithium (Li) abundance of their host stars is still an open question. To investigate the difference of the Li abundance between planet-host stars (HS) and isolated stars (IS) with no detected planets, we analyze a large sample of stars with temperatures ranging from 4600 to 6600 K and metallicity ranging from -0.55 to +0.50. The sample consists of 279 HS whose spectra were taken from the California-Kepler Survey (CKS), which followed up planets detected by Kepler, and 171 IS whose spectra were taken from the Keck archive. The non-local thermodynamic equilibrium (non-LTE) effects were taken into consideration. It is found that the distribution of Li abundances in both the HS and IS groups are generally consistent with each other. This suggests that the presence of Kepler-like planets does not have a significant impact on Li depletion. We also found that the non-LTE corrections can not be neglected for stars with A(Li) over ~ 2.5 dex.

We study gas mixing in a simulated Milky Way-mass galaxy's circumgalactic medium (CGM) using cosmological zoom-in simulations. We insert tracer dyes in the CGM with different gas flows (shearing, coherent, and static) and diverse physical properties to track gas mixing and find that shearing flows generally exhibit a higher dye spread. We then correlate the extent and shape of dye spread with local gas properties for a generalised understanding of gas mixing. The velocity dispersion and magnitude of traceless symmetric shear (shear deformation) tensors in small regions (<5 kpc) around the dye injection locations most accurately predict the extent of dye spread after 200 Myr (Spearman rank 0.87). We use this to determine diffusivity calibration constants for subgrid-scale mixing models. While the dye shape after 200 Myr aligns well with velocity and magnetic field dispersion, the best alignment occurs with the dispersion of eigenvectors representing stretching (from traceless symmetric shear tensor) and plane-of-rotation (from antisymmetric shear or vorticity tensor) in large regions (>10 kpc) around the dye injection locations. Therefore, shear-related statistics and velocity dispersion are the properties best able to predict both the extent and shape of mixed gas in our simulations. Our dyes spread roughly linearly with time, suggesting superdiffusive behaviour in the CGM, potentially due to turbulent and large-scale coherent flows. We calculate physical Reynolds numbers of 10^2-10^3 using gas properties in the CGM of our simulations, which a resolution of a few hundred pc could resolve, but our 1 kpc resolution leads to significant numerical mixing. However, our correlation results are robust to changes in resolution. These results can be used to predict diffusion coefficients from local gas properties to model processes like magnetic field diffusion, heat transport, and metal mixing.

During the fall of late-type galaxies into clusters, they can experiment a variety of evolutionary mechanisms according their local environment. Consequently, studying the UV emission and the cold gas of late-type galaxies provide key insights in the evolution of short-lived starburst and galaxy quenching. In this work, we conduted a study of two 28' fields observed with UVIT-AstroSat in the central region of the Abell cluster A496 ($z=0.033$), including HI, data from NRAO VLA. We reported 22 cluster members detected in FUV; all of them are detected in HI, or have upper limits for the HI-mass. We find our FUV detected galaxies generally have higher specific star formation rates than other star forming galaxies. Most of the FUV galaxies with masses above 10$^9 \mathrm{M}_{\odot}$,and showing high sSFR have no close neighbors, pointing at RPS as the dominant mechanism affecting them. In contrast, most of the low-mass FUV objects present at least one companion, suggesting that tidal interactions also play an important role in the triggering of infalling galaxies. Combining the FUV-SFR with the HI properties of the observed galaxies in A496 we identify an evolutionary sequence consisting of five stages: (1) Pre-triggering, (2) Initial SF-triggering, (3) Peak of star-formation, (4) SF-fading, and (5) SF-quenching. During this path, normal gas-rich objects reach a gas-deficiency phase with SFR well below the main sequence. This process, prior to becoming a full passive galaxy, can be accomplished within a few 10$^{8}$ yr.

TOI-2285 b is a sub-Neptune-sized planet orbiting a nearby M dwarf, discovered through the TESS photometric survey and ground-based follow-up observations. The planet was initially reported to have an orbital period of 27.27 d, making it one of the lowest temperature sub-Neptunes transiting a bright M dwarf. However, additional TESS data reveal that its true orbital period is 13.64 d, half the original value, resulting in a warmer equilibrium temperature (358 K) than previously estimated (284 K). The misidentification likely resulted from the low signal-to-noise ratio of individual transit signals and the limited number of transits observed by TESS at that time. This case highlights the importance of carefully considering harmonic solutions for similar cases. The additional TESS data also reveal another planetary candidate with an orbital period of 9.67 d and a radius of 1.5 $R_\oplus$, requiring validation in future studies.

Mass drainage is frequently observed in solar filaments. During filament eruptions, falling material most likely flows along magnetic field lines, which may provide important clues for the magnetic structures of filaments. Here we study three filament eruptions exhibiting significant mass draining, often manifested as falling threads at a constant speed ranging between 30--300 km/s. We found that most of the falling material lands onto the hooked segments of flare ribbons, only a small fraction lands inside the hooks, and almost none lands onto the straight segments of ribbons. Based on these observations we surmise that before eruptions most of the filament mass is entrained by field lines threading the quasi-separatrix layers (QSLs), which wrap around the filament field and whose footpoints are mapped by the hooked ribbons, and that the magnetic reconnection involving these field lines is the major cause of the mass drainage during eruptions. Additionally, the light curves of the hooked ribbons suggest that during eruptions the earlier (later) QSL boundary of filaments is threaded by mass-loaded (depleted) field lines. By assuming that the constant-speed motion is due to a drag force balancing the gravity, we proposed a simplified analytical model to estimate the density contrast of the falling material. The estimated density contrast is then fed into a numerical model, in which filament threads fall along vertical magnetic field lines through a gravitationally stratified atmosphere. The resultant falling speeds, however, are short of observed values, which calls for further investigations.

In Limongi et al. (2024) we presented and discussed the main evolutionary properties and final fate of stars in the mass range 7-15 Msun. The evolutions of those models were computed by means of a medium size nuclear network that guaranteed a proper calculation of the nuclear energy generation and hence a good modeling of the physical evolution of these stars. In the present paper, we extend this study by computing the detailed chemical yields of stars in the mass range 9-15 Msun, i.e., those stars that explode as core collapse supernovae (CCSNe). The explosive nucleosynthesis is then computed in the framework of the thermal bomb induced explosion by means of the HYPERION code (Limongi and Chieffi 2020). We find that: (1) the yields of the intermediate mass elements (i.e., O to P) show a steep decrease as the inital mass decreases; (2) the yields of s-weak component, i.e., those produced by the slow neutron captures from Ga to to the first neutron closure shell, decrease almost linearly as a function of the initial mass with respect to the ones produced by the more massive stars; (3) the global contribution of the stars in the mass range 9.22-13 Msun to the yields of a generation of massive stars averaged over a standard initial mass function is negligible for essentially all the isotopes. In spite of this, however, the models of stars in this mass range can be fundamental to interpret the observations of specific supernovae.

Cosmic strings, topological defects predicted by high-energy theories, may contribute to the late-time expansion of the Universe, effectively mimicking dynamical dark energy. We investigate four phenomenological extensions of the $\Lambda$CDM model involving a residual string network: (i) a non-relativistic component with positive energy density (Model~1), (ii) a velocity-dependent extension (Model~2), (iii) a non-relativistic string network with energy density allowed to take both positive and negative values (Model~3), and (iv) a general scenario with free energy and velocity parameters (Model~4). These models are constrained using \textit{Planck} CMB data, SDSS or DESI baryon acoustic oscillations, and Type Ia supernovae observations. Models~1 and~2 yield strong upper bounds on the string density, for example, $\Omega_{\mathrm{s}} < 0.00901$ at 95\% CL from the CMB+DESI+DESY5 combination for Model~2, and mildly shift the inferred value of $H_0$ upward, though they are not favored by Bayesian evidence. For the same combination, the bulk velocity is bound as $v_{\mathrm{s}} < 0.569$. Models~3 and~4 exhibit a consistent preference for slightly negative values of $\Omega_{\mathrm{s}}$, with CMB-only data yielding $\Omega_{\mathrm{s}} = -0.038^{+0.029}_{-0.022}$ and $v_{\mathrm{s}}< 0.574$ in Model~4, and a best-fit improvement of $\Delta \chi^2 = -6.07$. However, these improvements are not sufficient to overcome the Occam penalty, and the Bayesian evidence continues to favor $\Lambda$CDM. These findings demonstrate the power of current data to constrain exotic energy components and encourage further exploration of string-inspired extensions to $\Lambda$CDM, particularly those involving negative-tension networks.

The richest and largest structures in the cosmic web are galaxy superclusters, their complexes (associations of several almost connected very rich superclusters), and planes. Superclusters represent a special environment where the evolution of galaxies and galaxy groups and clusters differs from the evolution of these systems in a low-density environment. The richest galaxy clusters reside in superclusters. The richest superclusters in the nearby Universe form a quasiregular pattern with the characteristic distance between superclusters 120 - 140 $h^{-1}$Mpc. Moreover, superclusters in the nearby Universe lie in two huge perpendicular planes with the extent of several hundreds of megaparsecs, the Local Supercluster plane and the Dominant supercluster plane. The origin of these patterns in the supercluster distribution is not yet clear, and it is an open question whether the presence of such structures can be explained within the $\Lambda$CDM cosmological model. This review presents a brief story of superclusters, their discovery, definitions, main properties, and large-scale distribution.

Neutron stars (NSs) are thought to receive natal kicks at their formation in supernovae. In order to investigate the magnitude of these kicks, we analyze the proper motions and distance estimates -- either through parallax or dispersion measures -- of young isolated pulsars and infer their three-dimensional velocities relative to their local standard of rest. We find that the velocities based on parallax distances of pulsars younger than $10$ Myr follow a log-normal distribution with a mean of $\ln\left(v/\text{km}/\text{s}^{-1}\right)=5.60\pm0.12$ and a standard deviation of $0.68\pm0.10$, peaking at ${\sim}150$--$200$ km s$^{-1}$, which we adopt as our fiducial kick distribution. Using a previously established method that infers kick magnitudes through the eccentricity of Galactic trajectories, we also estimate the kick velocities of older pulsars, which we find to be consistent with our fiducial kick distribution. A log-normal fit to all pulsars with ages below $40$ Myr yields a more constraining (but possibly more prone to systematic errors) fit with mean and standard deviation of $5.67\pm0.10$ and $0.59\pm0.08$, respectively. Moreover, we (1) resolve the tension between our results and the Maxwellian distribution found by Hobbs et al. (2005), which has a ${\sim}50\%$ higher median velocity, by showing that their analysis is missing a Jacobian needed to correct for its logarithmic histogram bin sizes, and (2) argue that the bimodality found by others is not statistically significant and that previous results are consistent with our inferred kick distribution, effectively reconciling the literature on observed NS kicks.

Sulfur is essential for life, but its abundance and distribution in the interstellar medium remain uncertain, with over 90% of sulfur undetected in cold molecular clouds. Sulfur allotropes (S$_{\rm n}$) have been proposed as possible reservoirs, but the only detected interstellar molecule with a disulfide bond is S$_2$H in the Horsehead Nebula, making the estimation of sulfur chains abundances difficult. Here we present total-power ALMA images of H$_2$S, S$_2$H, SO$_2$, CO$^+$, and SH$^+$ towards the Horsehead nebula. These observations, with unprecedented sensitivity (rms $\sim$ 1.5 mK), provide the first detections of SH$^+$ and CO$^+$ in this region, together with the identification of a new S$_2$H line. The comparison of the spectroscopic images of H$_2$S, S$_2$H, SO$_2$, CO$^+$ and SH$^+$ shows that the S$_2$H emission originates from a warm gas layer adjacent to the photodissociation front. The emission peak of S$_2$H is offset from those of reactive ions such as SH$^+$, CO$^+$, and SO$^+$, suggesting that gas-phase reactions involving SH$^+$ and H$_2$S are not the dominant formation pathway of S$_2$H. Instead, we propose that S$_2$H is desorbed from irradiated grain surfaces by non-thermal processes. The SH$^+$ detection indicates that sulfur is not significantly depleted at the UV-irradiated edge of the molecular cloud, arguing against a major refractory sulfur reservoir in the interior of molecular clouds.

We report on optically selected stellar candidates of SARAO MeerKAT 1.3 GHz radio continuum survey sources of the Galactic plane. Stellar counterparts to radio sources are selected by cross-matching the MeerKAT source positions with \textit{Gaia} DR3, using two approaches. The first approach evaluated the probability of chance alignments between the radio survey and \textit{Gaia} sources and used AllWISE infrared colour-colour information to select potential stellar candidates. The second approach utilized a Monte Carlo method to evaluate the cross-matching reliability probability, based on populations of known radio-emitting stars. From the combined approaches, we found 629 potential stellar counterparts, of which 169 have existing SIMBAD classifications, making it the largest Galactic plane radio-optical crossmatch sample to date. A colour-magnitude analysis of the sample revealed a diverse population of stellar objects, ranging from massive OB stars, main-sequence stars, giants, young stellar objects, emission line stars, red dwarfs and white dwarfs. Some of the proposed optical counterparts include chromospherically/coronally active stars, for example RS CVn binaries, BY Dra systems, YSOs and flare stars, which typically exhibit radio emission. Based on Gaia's low-resolution spectroscopy, some of the stars show strong H$\alpha$ emission, indicating they are magnetically active, consistent with them being radio emitters. While MeerKAT's sensitivity and survey speed make it ideal for detecting faint radio sources, its angular resolution limits accurate counterpart identification for crowded fields such as the Galactic Plane. Higher frequency, and, thereby, better spatial resolution, radio observations plus circular polarization would be required to strengthen the associations.

We present two new radio continuum images obtained with Apertif at 1.4 GHz. The images, produced with a direction-dependent calibration pipeline, cover 136 square degrees of the Lockman Hole and 24 square degrees of the ELAIS-N fields, with an average resolution of 17x12" and residual noise of 33 uJy/beam. With the improved depth of the images we found in total 63692 radio sources, many of which are detected for the first time at this frequency. With the addition of the previously published Apertif catalog for the Bootes field, we cross-match with the LOFAR deep-fields value-added catalogs at 150 MHz, resulting in a homogeneous sample of 10196 common sources with spectral index estimates, one of the largest to date. We analyze and discuss the correlations between spectral index, redshift, linear sources size, and radio luminosity, taking into account biases of flux-density-limited surveys. Our results suggest that the observed correlation between spectral index and redshift of active galactic nuclei can be attributed to the Malmquist bias reflecting an intrinsic relation between radio luminosity and the spectral index. We also find a correlation between spectral index and linear source size with more compact sources having steeper spectra.

The circumstellar envelopes (CSEs) of asymptotic giant branch (AGB) stars harbour a rich variety of molecules and are sites of complex chemistry. Our current understanding of the circumstellar chemical processes of carbon-rich AGB stars is predominantly based on observations of a single star, IRC+10216, often regarded as an archetypical carbon star. We aim to estimate stellar and circumstellar properties for five carbon stars, and constrain their circumstellar CS abundances. This study compares the CS abundances among the sources, informs circumstellar chemical models, and helps to assess if IRC+10216 is a good representative of the physics and chemistry of carbon star CSEs. We modelled the spectral energy distributions (SEDs) and CO line emission to derive the stellar and outflow properties. Using these, we then retrieved CS abundance profiles with detailed radiative transfer modelling, imposing spatial and excitation constraints from ALMA and single-dish observations. We obtain good fits to the SEDs and CO lines for all sources and reproduce the CS line emission across various transitions and apertures, yielding robust estimates of the CS abundance profiles. Peak CS fractional abundances range from 1$\times$10$^{-6}$ - 4$\times$10$^{-6}$, with e-folding radii of 1.8$\times$10$^{16}$ - 6.8$\times$10$^{16}$ cm. We also derive reliable $^{12}$C/$^{13}$C and $^{32}$S/$^{34}$S ratios from CS isotopologue modelling. Our results refine previous single-dish CS abundance estimates and improve the relative uncertainty on the CS e-folding radius for IRAS 07454$-$7112 by a factor of $\sim$2.5. Chemical models reproduce our estimates of the CS radial extent, corroborating the CS photodissociation framework used therein. We find no significant differences between the derived CS abundance profiles for IRC+10216 and the rest of the sample, apart from the expected density-driven variations.

Time-series photometry at mid-infrared wavelengths is becoming a common technique to search for atmospheres around rocky exoplanets. This method constrains the brightness temperature of the planet to determine whether heat redistribution is taking place - indicative of an atmosphere - or whether the heat is reradiated from a low albedo bare rock. By observing at 15$\mu$m we are also highly sensitive to CO$_2$ absorption. We observed three eclipses of the rocky super-Earth LHS 1140c using MIRI/Imaging with the F1500W filter. We found significant variation in the initial settling ramp for these observations and identify a potential trend between detector settling and the previous filter used by MIRI. We analysed our data using aperture photometry but also developed a novel approach which joint-fits pixel light curves directly using a shared eclipse model and a flexible multi-dimensional Gaussian process which models changes in the PSF over time. We demonstrate using simulated data that our method has the ability to weight away from particular pixels which show increased systematics, allowing for the recovery of eclipse depths in a more robust and precise way. Both methods and an independent analysis detect the eclipse at $>5\sigma$ and are highly consistent with a low albedo bare rock. We recover a dayside brightness temperature of $T_\mathrm{day} = 561\pm44$ K, close to the theoretical maximum of $T_\text{day; max} = 537\pm9$ K. We rule out a wide range of atmospheric forward models to $>3\sigma$ including pure CO$_2$ atmospheres with surface pressure $\ge10$ mbar and pure H$_2$O atmospheres with surface pressure $\ge1$ bar. Our strict constraints on potential atmospheric composition, in combination with future observations of the exciting outer planet LHS 1140b, could provide a powerful benchmark to understand atmospheric escape around M dwarfs.

About half of the mass of all heavy elements with mass number A > 90 is formed through the slow neutron capture process (s-process), occurring in evolved asymptotic giant branch (AGB) stars with masses ~1-6 $\rm{M_{\odot}}$. The s-process can be studied by modeling the accretion of material from AGB stars onto binary barium (Ba), CH, and carbon-enhanced metal-poor (CEMP)-s stars. Comparing observationally derived surface parameters and 1D-LTE abundance patterns of s-process elements to theoretical binary accretion models, we aim to understand the formation of post-accretion systems. We explore the extent of dilution of the accreted material and describe the impact of convective mixing on the observed surface abundances. We compute a new grid of 2700 accretion models for low-mass post-accretion systems. A maximum-likelihood comparison determines the best fit models for observational samples of Ba, CH, and CEMP-s stars. We find consistent AGB donor masses in the mass range of 2-3 $\rm{M_{\odot}}$ across our sample of post-accretion binaries. We find the formation scenario for weak Ba stars is an AGB star transferring a moderate amount of mass ($\leq$0.5 $\rm{M_{\odot}}$) resulting in a ~2.0-2.5 $\rm{M_{\odot}}$ star. The strong Ba stars are best fit with lower final masses ~1.0-2.0 $\rm{M_{\odot}}$, and significant accreted mass ($\geq$0.5 $\rm{M_{\odot}}$). The CH and CEMP-s stars display lower final masses (~1.0 $\rm{M_{\odot}}$) and small amounts of transferred material (~0.1 $\rm{M_{\odot}}$). We find that Ba stars generally accrete more material than CEMP-s and CH stars. We also find that strong Ba stars must accrete more than 0.50 $\rm{M_{\odot}}$ to explain their abundance patterns, and in this limit we are unable to reproduce the observed mass distribution of strong Ba stars. The mass distributions of the weak Ba stars, CEMP-s, and CH stars are well reproduced in our modeling.

Complete, accurate, and precise catalogues of exoplanet host star (EHS) properties are essential to deriving high-quality exoplanet parameters. This paper aims at homogeneously parameterising EHS and their exoplanets, using Gaia and GSP-spec data. For 2573 EHS, we computed their L*, R*, and M*, with no prior assumption from stellar evolution models. Their Galactic positions, kinematic and orbital properties were derived. We re-scaled Mp and Rp of 3556 exoplanets, fully consistently with the stellar data. These new stellar Teff, L*, and R* are in rather good agreement with literature values but are more precise. In particular, R* are derived with typically less than 3% uncertainty; this reduces the uncertainty on Rp significantly and allows for a finer analysis of the decrease in the number of planets around the evaporation valley. Larger differences, however, were found for M*. We note that the EHS population is rather diverse in terms of the chemical and Galactic properties, although they are all found in the Solar vicinity, close to the Local spiral arm. Most EHS belong to the thin disc, but some older thick disc and halo members have also been identified. For the less massive planets, the average Rp increases with the metallicity of the host star. For giant planets, a dichotomy between dense and inflated planets is found. Denser planets tend to be more massive as the metallicity of the host star increases, while inflated planets are more massive for less metallic hosts. If confirmed, this bimodality implies that the diversity of giant exoplanets depends on their Galactic birth locus, with dense giant planets being more numerous than inflated ones when [M/H]>~1.5 times Solar, as in the central Milky Way regions. The Gaia spectroscopic catalogue of exoplanets and their host stars is large, homogeneous, and precise. It would be a useful added-value for planetary studies.

Precise cosmological constraints from type Ia supernovae require adequately accurate corrections for host-galaxy extinction. Modelling these corrections is challenged by the problem of disentangling supernova intrinsic colours from host-galaxy interstellar this http URL latter is commonly modelled in a probabilistic way assuming an exponential distribution exp(-E(B-V)/\tau) as a universal prior which is applied across all types of supernova host galaxies. We test the robustness of the exponential model and its universality against predictions based on simulating dust and type Ia supernova distributions in host galaxies of different morphological types. We find substantial differences between predicted interstellar reddening in late- and early-type host galaxies, primarily driven by the stellar-to-dust mass ratios. The mean simulated reddening in late-type galaxies matches well those derived from type Ia supernova observations, but it is significantly lower for early-type host galaxies. The reddening distributions exhibit an excess of sight lines with vanishing reddening with respect to the exponential model, although the difference is quite mild for late-type galaxies. On the other hand, the distribution may peak at E(B-V)>0 when considering a population of young type Ia supernovae originating from lower heights within the dust disc. We demonstrate that assuming a universal reddening prior distribution for modeling peak magnitude-colour relation, which is currently a common practice, gives rise to a spurious scatter in the derived extinction properties. It may also bias relative distances between supernovae originating from different host-galaxy populations. The discrepancy between the simulated reddening in average early-type host galaxies and the observed occurrence of reddened supernovae suggests that reddening does not originate from interstellar dust expected in these galaxies.

We consider the 21\,cm brightness temperature as a probe of the Hubble tension in the framework of an inhomogeneous cosmological model. Employing Buchert's averaging formalism to study the effect of inhomogeneities on the background evolution, we consider scaling laws for the backreaction and curvature consistent with structure formation simulations. We calibrate the effective matter density using MCMC analysis using Union 2.1 Supernova Ia data. Our results show that a higher Hubble constant ($\sim73$\,km/s/Mpc) leads to a shallower absorption feature in the brightness temperature versus redshift curve. On the other hand, a lower value ($\sim67$\,km/s/Mpc) produces a remarkable dip in the brightness temperature $T_{21}$. Such a substantial difference is absent in the standard $\Lambda$CDM model. Our findings indicate that inhomogeneities could significantly affect the 21\,cm signal, and may shed further light on the different measurements of the Hubble constant.

Relativistic reflection features in the X-ray spectra of accreting black holes are considered to be generated by the illumination of the accretion disk by the hot corona. In this work, we present a numerical method for the emission line profile and the reflection spectrum produced by an off-axis X-ray source. The X-ray source is considered as a point source, as in the lamppost scenario, except that it is located off-axis and moves at arbitrary velocity. The observed flux for the distant observer is calculated directly without priority evaluation of the emissivity on the accretion disk, which allows our model to be applicable to the point source that deviates from the axis of the black hole spins and moves with a velocity. To study the impact of the off-axis geometry on the measurement of source properties, we simulate observations for a black hole binary with NuSTAR and eXTP. We compare the simulation with the observation of the phase-resolved spectra of the low-frequency quasiperiodic oscillation observed by the Insight Hard X-ray Modulation Telescope. Due to the nonaxisymmetric illumination on the accretion disk, parameters of the model are not reproduced by the lamppost model, including the corona height, radial velocity, and the reflection fraction. On the other hand, all the model parameters are recovered through the off-axis model.

In this work, we develop a novel formalism to include the effect of electron density fluctuations at ultra small scales (well below the sound horizon at last scattering) on the observed anisotropies of the Cosmic Microwave Background (CMB). We treat the electron field as an independent stochastic variable and obtain the required ensemble-averaged photon Boltzmann equations using Ito calculus. Beyond changes to the average recombination history (which can be incorporated in the standard approach) our work identifies two new effects caused by the clumpiness of the medium. The first is a correction to the Thomson visibility function caused by correlations of the electron fluctuations along the line of sight, leading to an additional broadening of the visibility towards higher redshifts which causes extra damping and smearing of the CMB anisotropies. The second effect is a reduction of the effective scattering rate in the (pre-)recombination era that affects the photon transfer functions in a non-trivial manner. These new effects are subdominant in LCDM but can be significant in cosmologies with an early onset of structure formation (e.g., due to generation of enhanced small-scale power) as suggested by a number of indicators (e.g., the abundance of high redshift galaxies observed by JWST). We discuss the relevance of these new effects to the Hubble tension, finding that corrections which cannot be captured by simple modifications to the average recombination history arise. This highlights how important an understanding of the recombination process is in cosmological inference, and that a coordinated simulation and analysis campaign is required as part of the search for the origin of the various tensions in cosmology.

Intermittency has been studied extensively in the fast and slow solar winds but to a far lesser extent in interplanetary coronal mass ejections (ICMEs). While ICMEs are often characterized by their relatively smooth, large-scale magnetic flux rope structures, a spectrum of fluctuations is nonetheless present at smaller scales. We have examined kurtosis and its scaling exponents at magnetohydrodynamic inertial scales in 49 ICMEs observed between 0.25 and 1 au by Parker Solar Probe and Solar Orbiter, and compared the results to those obtained for the ICME sheath regions and ambient solar wind intervals. Kurtosis behaves similarly in all intervals studied and presents a universal behavior typical of intermittent time series. The ICMEs displayed a radially invariant level of intermittency, suggesting that they are relatively static, well-developed turbulent environments. In the sheath regions, the level of intermittency increased with distance, indicating that the turbulence is not yet fully developed at small heliocentric distances. In addition to intermittent fluctuations related to turbulence, the sheath regions may possess a population of non-turbulent structures that increase the absolute value of kurtosis.

Bright SMGs contribute significantly to the star formation rate (SFR) density (20-50\%) and stellar mass density ($\sim$ 30-50\%) at $z=$ 2-4 with SFRs$\ge 1000$ M$_\odot$,yr$^{-1}$ and stellar masses of $\sim 10^{11}$-$10^{12}$ M$_\odot$. The number of bright SMGs with such high SFRs is hard to reconcile with the standard models of galaxy formation and evolution. In this paper we provide evidence that, in a small sample of 12 bright SMGs, the SFRs derived from spectral energy distribution (SED) fitting are significantly higher than those obtained using low-frequency radio emission as a proxy for star formation. Using the International LOFAR Telescope (ILT), which allows imaging at 144 MHz with sub-arcsecond angular resolution, we have produced deep images of a small sample of bright SMGs in the North Ecliptic Pole (NEP) field extracted from the NEPSC2 survey. For all 12 SMGs, we find radio-emitting mid-infrared galaxies at distances from a few arcseconds down to sub-arcsecond scales from the SMG and/or the presence of a radio-emitting AGN. The SFRs derived from the radio emission of the SMG, disentangled from the AGN-related radio emission, are systematically lower by a factor of $\sim 5$ (median value) than those derived from the multi-band SED fitting. We discuss whether our assumptions might be, at least in part, responsible for the observed discrepancy. We argue that the radio-derived SFRs are not systematically underestimated but can be affected by a significant dispersion ($0.3-0.5$ dex). Considering these new SFR estimates, the offset of the specific SFR of the 12 bright SMGs from the star-forming galaxy main sequence ($\Delta\mathrm{(SSFR)}$) is significantly reduced, with all 12 bright SMGs which are only a factor of 2 more star-forming than the main sequence galaxies.

We present a compilation of observed recurrence times ($t_{\rm rec}$) and infer the corresponding local mass-accretion rates ($\dot m$) for type I X-ray bursts, milliHertz quasi-periodic oscillating sources and recurrent novae eruptions. We construct models of the $t_{\rm rec}-\dot m$ relation for accreting white dwarfs and neutron stars and find that both are roughly consistent with a global inverse linear relation, connecting for the first time thermonuclear runaways on neutron stars and white dwarfs. We find that theoretical models of pure He bursts are in agreement with the best $t_{\rm rec}$ measurements in ultra-compact X-ray binaries at low $\dot m$ (4U~$0614+09$ and 2S~0918-549). We suggest that the transient Z source XTE~J1701-462 is a slow rotator, based on its mHz QPO properties. Finally, we discuss the implications for thermonuclear ignition and point out that the difference in eruption/burst energy ($E_{b_{WD}}/E_{b_{NS}}=2\times 10^4$) is consistent with the difference in area between neutron stars and white dwarfs $\left((R_{WD}/R_{NS})^2=4\times 10^4\right)$. We conclude that ignitions of thermonuclear shell flashes on neutron stars and white dwarfs depend primarily on the specific mass accretion rate and do not depend on the nature of the underlying compact object.

The chemical composition of exoplanets is thought to be influenced by the composition of the disks in which they form. JWST observations have unveiled a variety of species in numerous nearby disks, showing significant variations in the C/O abundance ratio. However, little is known about the composition and C/O ratio of disks around young stars in clusters exposed to strong ultraviolet (UV) radiation from nearby massive stars, which are representative of the environments where most planetary systems form, including ours. We present JWST spectroscopy of d203-504, a young 0.7 $\rm M_{\odot}$ star in the Orion Nebula with a 30 au disk irradiated by nearby massive stars. These observations reveal spectroscopic signatures of CO, H$_2$O, CH$_3^+$, and PAHs. Water and CO are detected in absorption in the inner disk ($r\lesssim 1$ au), where the estimated gas-phase C/O ratio is 0.48, consistent with the Solar value and that of the Orion Nebula. In contrast, \ch{CH3+} and PAHs are found in the extended surface layers of the disk. These results suggest that gas in the inner disk is chemically shielded from UV radiation while the surface layers of the disk experience UV-induced chemistry, potentially depleting their carbon content.

Recently, the baryon acoustic oscillations (BAO) measurements from the DESI survey have suggested hints of dynamical dark energy, challenging the standard $\Lambda $CDM model. In this work, we adopt an improved inverse distance ladder approach based on the latest cosmological data to provide a model-independent perspective, employing a global parametrization based on cosmic age (PAge). Our analysis incorporates DESI DR2 BAO measurements, cosmic chronometer (CC) data, and type Ia supernovae (SNe) observations from either the DESY5 or PantheonPlus datasets. For the DESY5+DESI DR2+CC datasets, we obtain $H_0 = 67.91 \pm 2.33~\mathrm{km~s^{-1}~Mpc^{-1}}$. This value is consistent with the Planck 2018 result, while shows $2.0 \sigma$ tension with the SH0ES measurement. Furthermore, by mapping specific cosmological models into PAge approximation parameter space $(p_{\mathrm{age}}, \eta)$, our model-independent analysis reveals a notable deviation from the $\Lambda \mathrm{CDM}$ model, as indicated by the DESY5 and DESI DR2 datasets. Finally, DESY5+DESI DR2+CC datasets provide nearly decisive evidence favoring the PAge model over the standard $\Lambda \mathrm{CDM}$ model. These findings highlight the need for further investigation into the expansion history to better understand the deviations from the $\Lambda \mathrm{CDM}$ model.

The solar wind (SW) is a vital component of space weather, providing a background for solar transients such as coronal mass ejections, stream interaction regions, and energetic particles propagating toward Earth. Accurate prediction of space weather events requires a precise description and thorough understanding of physical processes occurring in the ambient SW plasma. Ensemble simulations of the three-dimensional SW flow are performed using an empirically-driven magnetohydrodynamic heliosphere model implemented in the Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS). The effect of uncertainties in the photospheric boundary conditions on the simulation outcome is investigated. The results are in good overall agreement with the observations from the Parker Solar Probe, Solar Orbiter, Solar Terrestrial Relations Observatory, and OMNI data at Earth, specifically during 2020-2021. This makes it possible to shed more light on the properties of the SW propagating through the heliosphere and perspectives for improving space weather forecasts.

The impact of metallicity on the Cepheid Leavitt law (LL) and, in turn, the Hubble constant, has been the subject of much recent debate. Here, we present a comprehensive analysis of metallicity effects on Cepheid LLs based on synthetic Cepheid populations computed using Geneva models and the SYCLIST tool. We computed 296 co-eval populations in the age range of 5-300 Myr for metallicities representative of the Sun, the LMC, and the SMC ($Z \in [0.014, 0.006, 0.002]$). We computed LLs in fourteen optical-to-infrared passbands and five reddening-free Wesenheit magnitudes. All Cepheid populations take into account distributions of rotation rates and companion stars. We show excellent agreement between the predicted populations and key observational constraints from the literature. Our simulations predict a significant LL slope-metallicity dependence ($\beta_{\rm M} > 0$) that renders LLs steeper at lower metallicity at all wavelengths. Importantly, $\beta_{\rm M} \ne 0$ implies that the intercept-metallicity dependence, $\alpha_{\rm M}$, depends on pivot period; an issue not previously considered. Comparison with $\alpha_{\rm M}$ measurements in individual passbands reported in the literature yields acceptable agreement on the order of agreement found among different observational studies. The wavelength dependence and magnitude of the disagreement suggests a possible origin in reddening-related systematics. Conversely, we report excellent agreement between our $\alpha_{\rm M} = -0.20 \pm 0.03$ mag dex$^{-1}$ and the value determined by the SH0ES distance ladder in the reddening-free H-band Wesenheit magnitude ($-0.217 \pm 0.046$), the currently tightest and conceptually simplest empirical constraint.

In this paper, we investigate the gravitational wave (GW) emission from stars tidally disrupted by black holes (TDEs), using a semi-analytical approach. Contrary to previous works where this signal is modeled as a monochromatic burst, we here take into account all its harmonic components. On top of this, we also extend the analysis to a population of repeated-partial TDEs, where the star undergoes multiple passages around the black hole before complete disruption. For both populations, we estimate the rate of individual GW-detections considering future observatories like LISA and a potential deci-Hertz (dHz) mission, and derive the GW background from these sources. Our conclusions, even if more conservative, are consistent with previous results presented in literature. In fact, full disruptions of stars will not be seen by LISA but will be important targets for dHz observatories. In contrast, GWs from repeated disruptions will not be detectable in the near future.

The formation of primordial black holes or other dark matter relics from amplified density fluctuations in the early universe may also generate scalar-induced gravitational waves (GW), carrying vital information about the primordial power spectrum and the early expansion history of our universe. We present a Bayesian approach aimed at reconstructing both the shape of the scalar power spectrum and the universe's equation of state from GW observations, using interpolating splines to flexibly capture features in the GW data. The optimal number of spline nodes is chosen via Bayesian evidence, aiming at balancing complexity of the model and the fidelity of the reconstruction. We test our method using both representative mock data and recent Pulsar Timing Array measurements, demonstrating that it can accurately reconstruct the curvature power spectrum as well as the underlying equation of state, if different from radiation.

The recent discovery of a large number of massive black holes within the first two billion years after the Big Bang, as well as their peculiar properties, have been largely unexpected based on the extrapolation of the properties of luminous quasars. These findings have prompted the development of several theoretical models for the early formation and growth of black holes, which are, however, difficult to differentiate. We report the metallicity measurement around a gravitationally lensed massive black hole at redshift 7.04, hosted in a galaxy with very low dynamical mass. The weakness of the [OIII]5007 emission line relative to the narrow Hbeta emission indicates an extremely low chemical enrichment, less than 0.01 solar. We argue that such properties cannot be uncommon among accreting black holes around this early cosmic epoch. Explaining such a low chemical enrichment in a system that has developed a massive black hole is challenging for most theories. Models assuming heavy black hole seeds (such as Direct Collapse Black Holes) or super-Eddington accretion scenarios struggle to explain the observations, although they can potentially reproduce the observed properties in rare cases. Models invoking "primordial black holes" (i.e. putative black holes formed shortly after the Big Bang) may potentially explain the low chemical enrichment associated with this black hole.

Machine learning can accelerate cosmological inferences that involve many sequential evaluations of computationally expensive data vectors. Previous works in this series have examined how machine learning architectures impact emulator accuracy and training time for optical shear and galaxy clustering 2-point function. In this final manuscript, we explore neural network performance when emulating Cosmic Microwave Background temperature and polarization power spectra. We maximize the volume of applicability in the parameter space of our emulators within the standard $\Lambda$-cold-dark-matter model while ensuring that errors are below cosmic variance. Relative to standard multi-layer perceptron architectures, we find the dot-product-attention mechanism reduces the number of outliers among testing cosmologies, defined as the fraction of testing points with $\Delta \chi^2 > 0.2$ relative to \textsc{CAMB} outputs, for a wide range of training set sizes. Such precision enables attention-based emulators to be directly applied to real data without requiring any additional correction via importance sampling. Combined with pre-processing techniques and optimized activation and loss functions, attention-based models can meet the precision criteria set by current and future CMB and lensing experiments. For each of Planck, Simons Observatory, CMB S4, and CMB HD, we find the fraction of outlier points to be less than $10\%$ with around $2\times10^5$ to $4\times10^5$ training data vectors. We further explore the applications of these methods to supernova distance, weak lensing, and galaxy clustering, as well as alternative architectures and pre-processing techniques.

We present the ExoPhoto database (this https URL), an extension of the ExoMol database, specifically developed to address the growing need for high-accuracy, temperature-dependent photodissociation cross section data towards short-UV wavelengths. ExoPhoto combines theoretical models from three major computational databases (ExoMol, UGAMOP and PhoMol) and experimental datasets from two experimental groups, providing extensive wavelength and temperature coverage. ExoPhoto currently includes photodissociation data for 20 molecules: AlH, HCl, HF, MgH, OH, NaO, MgO, O2, AlCl, AlF, CS, HeH+, CO, CO2, H2O, SO2, C2H2, C2H4, H2CO, and NH3, derived from theoretical models and supported by experimental data from 5 databases. ExoPhoto also includes detailed data on branching ratios and quantum yields for selected datasets. The data structure of ExoPhoto follows the ExoMol framework, with a consistent naming convention and hierarchical JSON-based organization. Photodissociation cross sections are stored in a set of .photo files which provide data as a function of wavelength with one file for each target molecule temperature. Future developments aim to include more photodissociation cross section data and to provide data for molecules in non-local thermodynamic equilibrium (non-LTE). These will expand the utility of ExoPhoto for advanced astrophysical, planetary modeling and industrial applications.

Little red dots (LRDs), a population of active galactic nuclei (AGNs) recently identified by JWST, are characterized by their compact morphology and red optical continuum emission, which is often interpreted as evidence for significant dust extinction of $A_V \gtrsim 3$ mag. However, the dust-reddened AGN scenario is increasingly challenged by their faint near-to-far infrared emission and a potential "dust budget crisis" in cases when the host galaxy is either undetectably low-mass or absent. In this study, we re-evaluate the dust extinction level in LRDs by modeling the UV-to-infrared spectra for various extinction laws and a broad range of dusty distribution parameters. Comparing the predicted infrared fluxes with observational data from the JWST MIRI, Herschel, and ALMA, our analysis finds that the visual extinction is tightly constrained to $A_V \lesssim 1.0$ mag for A2744-45924 and $A_V \lesssim 1.5$ mag for RUBIES-BLAGN-1 under the SMC extinction laws, with slightly weaker constraints for those with gray extinction in the UV range. The revised $A_V$ values yield a radiative efficiencies of $10\%$ for the LRD population, easing the tension with the Soltan argument for the bulk AGN population at lower redshifts. Moreover, this moderate extinction (or dust-free) scenario, with reprocessed emission spectra testable by future far-infrared observatories, provides a paradigm shift in understanding their natures, environments, and evolutionary pathways of massive black holes in the early universe.

In contrast to the large-scale primordial power spectrum $\mathcal{P}_{\zeta}(k)$ and primordial non-Gaussianity $f_{\mathrm{NL}}$, which are strictly constrained, the small-scale $\mathcal{P}_{\zeta}(k)$ and $f_{\mathrm{NL}}$ remain less restricted. Considering local-type primordial non-Gaussianity, we study the PBH and SIGW caused by large-amplitude small-scale primordial power spectrum. By analyzing current observational data from PTA, CMB, BAO, and abundance of PBH, and combining them with the SNR analysis of LISA, we rigorously constrain the parameter space of $\mathcal{P}_{\zeta}(k)$ and $f_{\mathrm{NL}}$. Furthermore, we examine the effects of different shapes of the primordial power spectrum on these constraints and comprehensively calculate the Bayes factors for various models. Our results indicate that SIGW generated by a monochromatic primordial power spectrum are more likely to dominate current PTA observations, with the corresponding constraint on the primordial non-Gaussian parameter being $-10.0<f_{\mathrm{NL}}<1.2$.

Highly-differential spectroscopic studies have revealed that the Sun is deficient in refractory elements relative to solar twins. To investigate the role of giant planets on this signature, we present a high precision abundance analysis of HARPS spectra for 50 F- and G-type stars spanning -0.4<[Fe/H]<+0.5. There are 29 stars in the sample which host planets of masses > 0.01 MJup. We derive abundances for 19 elements, and apply corrections to 14 of them for systematic errors associated with one dimensional (1D) model atmospheres, or the assumption of local thermodynamic equilibrium (LTE), or both. We find that, among the solar twins in our sample, the Sun is Li poor in comparison to other stars at similar age, in agreement to previous studies. The sample shows a variety of trends in elemental abundances as a function of condensation temperature. We find a strong correlation in these trends with [Fe/H], with a marginally-significant difference in the gradients for stars with and without giants planets detected, that increases after applying 3D and non-LTE corrections. Our overall results suggests that the peculiar composition of the Sun is primarily related to Galactic chemical evolution rather than the presence of giant planets.

I propose a scenario that allows white dwarfs (WDs) to launch relatively powerful jets when they enter a common envelope evolution (CEE) or experience a grazing envelope evolution (GEE) with a red giant branch star (RGB) or an asymptotic giant branch (AGB) star. In this, still a speculative scenario, the accretion for a time is mainly onto an accretion disk with a radius of ~ 1Ro that increases in mass. The accretion disk launches the powerful two opposite jets by releasing gravitational energy, up to several times super-Eddington, as its mass increases. The jets that the disk launches remove high-entropy gas from the disk's outskirts and the envelope that the WD inflates because of nuclear burning on its surface. The motivations to allow WDs to launch powerful jets are recent findings, from the morphologies of post-CEE planetary nebulae, that jets play a major role in the CEE and the accumulating evidence that jets power luminous red novae by jets, as their morphologies indicate. I strengthen my call to include jets in the simulation and modeling of the CEE, consider the GEE as a phase preceding the CEE in many (but not all) cases, and include jets as a major ingredient in modeling and stimulating all energetic luminous red novae.

Despite their somewhat-frequent appearance in EUV imaging of off-limb flares, the origins of Supra-Arcade Downflows (SADs) remain a mystery. Appearing as dark, tendril-like downflows above growing flare loop arcades, SADs themselves are yet to be tied into the standard model of solar flares. The uncertainty of their origin is, in part, due to a lack of spectral observations, with the last published SAD spectral observations dating back to the Solar and Heliospheric Observatory / Solar Ultraviolet Measurements of Emitted Radiation (SOHO/SUMER) era in 2003. In this work, we present new observations of SADs within an M-class solar flare on April 2nd, 2022, observed by the Hinode EUV Imaging Spectrometer (EIS) and NASA Solar Dynamics Observatory. We measure Fe XXIV 192.02 Angstrom Doppler downflows and non-thermal velocities in the low-intensity SAD features, exceeding values measured in the surrounding flare fan. The ratio of temperature-sensitive Fe XXIV 255.11 Angstrom and Fe XXIII 263.41 Angstrom lines also allow the measurement of electron temperature, revealing temperatures within the range of the surrounding flare fan. We compare EIS line-of-sight Doppler velocities with plane-of-sky velocities measured by AIA, to construct the 3D velocity profile of four prominent SADs, finding evidence for their divergence above the flare loop arcade - possibly related to the presence of a high altitude termination shock. Finally, we detect 'stealth' SADs, which produce SAD-like Doppler signals, yet with no change in intensity.

Much of our intuition about Effective Field Theories (EFTs) stems from their formulation in flat spacetime, yet EFTs have become indispensable tools in cosmology, where time-dependent backgrounds are the norm. In this work, we demonstrate that in spacetimes undergoing significant expansion-such as accelerated FLRW and de Sitter backgrounds-the contributions of operators with mass dimension $\Delta$ to physical observables grow factorially with $\Delta$ at fixed couplings. This behavior stands in stark contrast to expectations from flat spacetime. As a result, the cosmological EFT expansion is generally asymptotic rather than convergent, even at tree level. To illustrate this phenomenon, we analyze simple toy models involving a massless or conformally coupled scalar field interacting with a heavy scalar with zero or infinite sound speed. We demonstrate that meaningful EFT predictions can still be extracted via appropriate resummation techniques, performed in both Fourier and Mellin-momentum space. In the infinite sound speed limit, where the heavy field is effectively non-dynamical, the resummed EFT reproduces the exact result of the full theory. In other cases, the EFT captures only the local part of the dynamics, omitting nonlocal terms, which are exponentially suppressed in the large-mass limit for the Bunch-Davies state.

The first direct detection of gravitational waves by the LIGO collaboration, GW150914, marked the start of a new exciting era in astronomy, enabling the study of the Universe through a new messenger. Since then, the field has grown rapidly, with the development of increasingly more sophisticated techniques to detect, analyze and interpret the signals. In this paper we revisit GW150914, presenting updated estimates of its source parameters using a waveform model developed within the EOB formalism, able to describe gravitational-wave emission from generic non-circular, non-planar binaries. We provide a comprehensive analysis of the signal and its properties, considering and contrasting various scenarios for the source: from the simplest, aligned-spin quasi-circular binary black hole merger, to more complex scenarios, including precession, eccentricity or both. Unsurprisingly, we find that the signal is consistent with a quasi-circular ($e < 0.08$ at $15$ Hz), slowly spinning $(\chi_{\rm eff} = -0.03^{+0.12}_{-0.13})$ binary black hole merger, a-posteriori validating a considerable body of works. This is the first analysis performed with an inspiral-merger-ringdown model containing both eccentricity and precession.

The study of shear layer instability in compressible flows is key to understanding phenomena from aerodynamics to astrophysical jets. Blumen's seminal paper [``Shear layer instability of an inviscid compressible fluid," J. Fluid Mech. {\bf 40}, 769--781 (1970)] established a linear stability framework for inviscid compressible shear flows, emphasizing velocity gradients and compressibility effects. However, the nonlinear regime remains insufficiently explored. This research extends Blumen's framework by conducting a weakly nonlinear stability analysis using the method of multiple scales to derive amplitude equations, such as the Landau-Stuart and complex Landau equations. Perturbation variables are expanded in a power series to capture amplitude evolution beyond linear theory. Finite boundary conditions are incorporated to enhance physical applicability. The study analyzes how compressibility and Mach number influence nonlinear saturation, revealing Mach-dependent bifurcations in Kelvin-Helmholtz instability (KHI) with alternating stable and unstable regimes. Phase portraits and trajectories illustrate transitions, saturation, and spiral decay, which are relevant to astrophysical shear flows. Bifurcation analysis reveals both supercritical and subcritical Hopf behavior in compressible shear flows, underscoring the importance of nonlinear effects in the onset and evolution of flow instabilities. Qualitative and quantitative results of instability evolution have been shown from nonlinear stability analysis. This work bridges the critical gap between linear and fully nonlinear stability analyses by offering a systematic weakly nonlinear framework and the nonlinear dynamics for compressible shear layers. It generalizes the earlier linear results and provides new predictions about bifurcation behavior and long-time state selection in compressible flows.

We present a scenario in which primordial black holes (PBHs) form in a post-inflationary radiation-dominated (RD) phase and then experience significant accretion during a phase of early matter dominated (EMD). We show that PBH masses could grow by up to two orders of magnitude. Restricting to the linear perturbation regime, we compute the gravitational wave (GW) spectrum that features two peaks. The high-frequency peak is associated with the PBH formation in the RD phase, while the low-frequency peak is due to the sudden transition from EMD to the later, standard RD phase. We identify a PBH mass range where one or both peaks can be observed by a combination of different GW detectors. Finally, we show the signal-to-noise ratio of the total GW spectrum for PBHs in the asteroid mass window, where they could comprise the totality of dark matter.

The Laser Interferometer Space Antenna (LISA) will enable direct observations of low-frequency gravitational waves, offering unprecedented insight into astrophysical and cosmological phenomena. LISA's heterodyne interferometric measurement system requires phase-locking five of its six onboard lasers with tunable frequency offsets to ensure that all beatnotes remain within the metrology system's operational range, despite Doppler-induced frequency shifts. The selection of these offset frequencies -- collectively forming a frequency plan -- is a complex optimization problem constrained by the spacecraft's orbital dynamics and instrument limitations. While previous work established an algorithmic solution for deriving time-dependent frequency plans, this study takes a complementary approach by systematically analyzing and cataloging all possible laser locking configurations. We present an automated method to explore, validate, and classify viable locking schemes, identifying 36 unique non-frequency-swapping configurations and 72 additional frequency-swapping configurations for an arbitrary choice of primary laser. This exhaustive classification provides a foundation for frequency planning across the full range of operational scenarios.

Sterile neutrinos are compelling dark matter (DM) candidates, yet the minimal production mechanism solely based on active ($\nu_a$)-sterile ($\nu_s$) oscillations is excluded by astrophysical observations. Non-standard self-interactions in either active ($\nu_a-\nu_a$) or sterile ($\nu_s-\nu_s$) sector are known to alter the sterile neutrino DM production in the early Universe, which could alleviate the tension with astrophysical constraints to some extent. Here we propose a novel solution where scalar-mediated non-standard interactions between active and sterile neutrinos ($\nu_a-\nu_s$) generate new production channels for $\nu_s$, independent of the active-sterile mixing and without the need for any fine-tuned resonance or primordial lepton asymmetry. This framework enables efficient sterile neutrino DM production even at vanishingly small mixing angles and opens up new viable regions of parameter space that can be tested with future $X$-ray and gamma-ray observations.

The proposed Carbon Investigation (Carbon-I) Imaging Spectrometer is designed to measure variations of greenhouse gases in Earth's atmosphere. The instrument will survey the Earth from its own spacecraft at an altitude of approximately 610 km. It will use a coarse ground sampling distance (GSD) of <400 m in global mode for land and coastal monitoring and finer 35 m GSD in target mode to sample key regions. The identification and quantification of greenhouse gases require continuous spectral sampling over the 2040-2380 nm wavelength range with <1 nm spectral sampling. The proposed design builds upon Jet Propulsion Laboratory's (JPL) experience of spaceflight Dyson imaging spectrometers to achieve spectral sampling of 0.7 nm per pixel. This paper presents the proposed Carbon-I optical design comprised of a freeform three-mirror anastigmat telescope that couples to a F/2.2, highly uniform Dyson-inspired imaging spectrometer. The high uniformity and throughput enables Carbon-I to measure Earth's greenhouse gas concentrations with unprecedented precision and spatial sampling.

We formulate an approach to quantum gravity, called the ring paradigm. Gravity is mediated superluminally, and the graviton is described as a phonon on the grid of matter in the Universe. This theory has very interesting applications to cosmology and would ultimately solve the old problem of the cosmological constant. It further gives new impulses to the scalar field theories because the gravitational ring decays to some phantom field. As is obvious, we radically break the Lorentz invariance, which means that some generalization of the Haag-Lopuszanski-Sohnius theorem in quantum field theory is possible.

Many years of development have gone into producing instruments that meet the required noise performance of the LISA interferometric detection system. Concurrently, software simulations have been used to extensively develop the data analysis libraries to be used in the LISA pipeline, not least among which are the Time Delay Interferometry (TDI) algorithms. To bridge the gap between these two, we are developing a hardware-in-the-loop testbed to apply realistic, time-varying delays to signals traveling between phasemeters. We have shown that the testbed adds a sufficiently low amount of noise across the entire relevant LISA spectrum. We have also injected a realistic gravitational wave signal, generated via the LISA Data Challenge codebase, and successfully extracted it using TDI to remove the obscuring frequency noise of the carrier signal. Future efforts will expand the testbed to create a representative simulation of the entire LISA constellation, with an eye towards its use as a tool to aid in the development of the LISA data analysis pipelines.