Papers for Thursday, Jun 05 2025

In this paper, we use genetic algorithms, a specific machine learning technique, to achieve a model-independent reconstruction of $f(T)$ gravity. By using $H(z)$ data derived from cosmic chronometers and radial Baryon Acoustic Oscillation method, including the latest Dark Energy Spectroscopic Instrument (DESI) data, we reconstruct the Hubble rate which is the basis parameter for reconstructing $f(T)$ gravity without any assumptions. In this reconstruction process, we use the current value of the Hubble rate, $H_0$, derived by genetic algorithms. The reconstructed $f(T)$ function is consistent with the standard $\Lambda$CDM cosmology within the 1$\sigma$ confidence level across a broad temporal range. The mean $f(T)$ curve, adopting a quadratic form, prompts us to parametrize it using a second degree polynomial. This quadratic deviation from the $\Lambda$CDM scenario is mildly favored by the data.

Every halo of collisionless dark matter contains a $\rho = A r^{-1.5}$ density cusp at its center. This prompt cusp is a relic of the halo's earliest moments and has a mass comparable to the cutoff scale in the spectrum of initial density perturbations. In this work, we provide a framework to predict for each halo the coefficient $A$ of its central cusp. We also present a "cusp-NFW" functional form that accurately describes the density profile of a halo with a prompt cusp at its center. Accurate characterization of each halo's central cusp is of particular importance in the study of warm dark matter models, for which the spectral cutoff is on an astrophysically relevant mass scale. To facilitate easy incorporation of prompt cusps into any halo modeling approach, we provide a code package that implements the cusp-halo relation and the cusp-NFW density profile.

We study the Ly$\alpha$ escape fraction, $f_{\mathrm{esc}}^{\mathrm{Ly\alpha}}$ of Ly$\alpha$ emitters (LAEs) identified by Subaru/HSC narrowband imaging at $z = 2.2-6.6$, using publicly available deep imaging data from HST and JWST. We perform SED fitting for 127 LAEs at $0.4 - 5.0\, \mathrm{\mu m}$ to estimate their physical properties robustly, and confirm that two distinct LAE populations exist: young LAEs (< 100 Myr) and old LAEs (> 100 Myr). Young LAEs are characterized by burst-like star formation activity and low dust content, significantly differing from Lyman-break galaxies (LBGs) at the same stellar mass, while old LAEs show similar star formation activity to LBGs, yet with lower dust content and more compact morphology in rest-UV/optical than LBGs. The $f_{\mathrm{esc}}^{\mathrm{Ly\alpha}}$ of LAEs is anticorrelated with stellar mass, and this correlation is found to depend on the age of LAEs, such that old LAEs show a weaker anticorrelation than young LAEs, and tend to exhibit higher $f_{\mathrm{esc}}^{\mathrm{Ly\alpha}}$ than young LAEs at a given stellar mass. This implies that Ly$\alpha$ photons escape more efficiently from old LAEs, possibly through the low-density channels of HI and dust created by outflows. The average $f_{\mathrm{esc}}^{\mathrm{Ly\alpha}}$ of young LAEs remains nearly constant at $\sim40$% at $z = 2.2-6.6$, suggesting that the previously observed evolution of global $f_{\mathrm{esc}}^{\mathrm{Ly\alpha}}$ of star-forming galaxies (SFGs) is due to the changes in the LAE fraction among the SFGs. Converting $f_{\mathrm{esc}}^{\mathrm{Ly\alpha}}$ to Lyman continuum escape fraction using empirical relations, we demonstrate that LAEs alone can supply the ionizing photons necessary for reionization at $z\sim6$, causing rapid and late reionization.

We present COSMOS2025, the COSMOS-Web catalog of photometry, morphology, photometric redshifts and physical parameters for more than 700,000 galaxies in the Cosmic Evolution Survey (COSMOS) field. This catalog is based on our \textit{James Webb Space Telescope} 255\,h COSMOS-Web program, which provides deep near-infrared imaging in four NIRCam (F115W, F150W, F277W, F444W) and one MIRI (F770W) filter over the central $\sim 0.54 {\, \rm deg}^2$ ($\sim 0.2 {\, \rm deg}^2$ for MIRI) in COSMOS. These data are combined with ground- and space-based data to derive photometric measurements of NIRCam-detected sources using both fixed-aperture photometry (on the space-based bands) and a profile-fitting technique on all 37 bands spanning 0.3-8 micron. We provide morphology for all sources from complementary techniques including profile fitting and machine-learning classification. We derive photometric redshifts, physical parameters and non-parametric star formation histories from spectral energy distribution (SED) fitting. The catalog has been extensively validated against previous COSMOS catalogs and other surveys. Photometric redshift accuracy measured using spectroscopically confirmed galaxies out to $z\sim9$ reaches $\sigma_{\rm MAD} = 0.012$ at $m_{\rm F444W}<28$ and remains at $\sigma_{\rm MAD} \lesssim 0.03$ as a function of magnitude, color, and galaxy type. This represents a factor of $\sim 2$ improvement at 26 AB mag compared to COSMOS2020. The catalog is approximately 80\% complete at $\log(M_{\star}/{\rm M}_{\odot}) \sim 9$ at $z \sim 10$ and at $\log(M_{\star}/{\rm M}_{\odot}) \sim 7$ at $z \sim 0.2$, representing a gain of 1\,dex compared to COSMOS2020. COSMOS2025 represents the definitive COSMOS-Web catalog. It is provided with complete documentation, together with redshift probability distributions, and it is ready for scientific exploitation today.

The Little Red Dots (LRDs) are high-redshift galaxies uncovered by JWST, characterized by small effective radii ($R_{\rm eff} \sim 80-300$ pc), number densities intermediate between typical galaxies and quasars, and a redshift distribution peaked at $z \sim 5$. We present a theoretical model in which the LRDs descend from dark matter halos in the extreme low-spin tail of the angular momentum distribution. Within this framework, we explain their three key observational signatures: (i) abundance, (ii) compactness, and (iii) redshift distribution. Our model focuses on observed, not modeled, properties; it is thus independent of whether they are powered primarily by a black hole or stars. We find that the assumption that the prototypical LRD at $z \sim 5$ originates from halos in the lowest $\sim 1\%$ of the spin distribution is sufficient to reproduce both their observed number densities and physical sizes. The redshift evolution of their observability is driven by the interplay between the evolving compact disk fraction and cosmological surface brightness dimming. This effect leads to a well-defined "LRDs Era" at $4<z<8$, during which the LRDs are common and detectable; at $z<4$, they are bright but rare, while at $z>8$, they are common but faint. Finally, we test the predicted redshift trend against observational data, finding excellent agreement. Additional observational support comes from their excess small-scale clustering and spectral signatures of extreme core densities, both of which are expected outcomes of galaxy formation in low-spin halos. These findings suggest that the LRDs are not a fundamentally distinct population but the natural manifestation of galaxies forming in the rarest, lowest angular momentum environments.

Dwarf galaxies provide powerful laboratories for studying galaxy formation physics. Their early assembly, shallow gravitational potentials, and bursty, clustered star formation histories make them especially sensitive to the processes that regulate baryons through multi-phase outflows. Using high-resolution, cosmological zoom-in simulations of a dwarf galaxy from \textit{the Pandora suite}, we explore the impact of stellar radiation, magnetic fields, and cosmic ray feedback on star formation, outflows, and metal retention. We find that our purely hydrodynamical model without non-thermal physics - in which supernova feedback is boosted to reproduce realistic stellar mass assembly - drives violent, overly enriched outflows that suppress the metal content of the host galaxy. Including radiation reduces the clustering of star formation and weakens feedback. However, the additional incorporation of cosmic rays produces fast, mass-loaded, multi-phase outflows consisting of both ionized and neutral gas components, in better agreement with observations. These outflows, which entrain a denser, more temperate ISM, exhibit broad metallicity distributions while preserving metals within the galaxy. Furthermore, the star formation history becomes more bursty, in agreement with recent JWST findings. These results highlight the essential role of non-thermal physics in galaxy evolution and the need to incorporate it in future galaxy formation models.

The Kepler field hosts the best studied sample of field star rotation periods. However, due to Kepler's large 4" pixels, many of its light curves are at high risk of contamination from background sources. The new Kepler Bonus Background light curves are de-blended using a PSF algorithm, providing light curves of over 400,000 new background sources in addition to over 200,000 re-analyzed Kepler prime targets. These light curves provide the opportunity to search for new rotation periods. Here we apply a convolutional neural network trained on synthetic spot-modulated light curves to regress rotation periods from the Kepler Bonus light curves. We obtained periods for 32,159 total sources, 19,650 of which had previously been measured and 9,811 of which are new periods for both Kepler prime and background sources. Our method also detected 608 pulsation frequencies from asteroseismic oscillations in red giants. We validate our Kepler prime periods against literature values and present the full period sample. We find excellent agreement with previously-known literature periods, validating deep learning as a viable class of period determination methods. Comparing the periods and light curves of foreground-background pairs, we find that as many as 63% of periodic background light curves are still blended with the foreground, highlighting limitations of the de-blending technique.

We present the data reduction methodology used for the COSMOS-Web survey JWST NIRCam data. Covering 0.54 deg^2 with four broadband filters (F115W, F150W, F277W, F444W) and a total exposure time of approximately 270 hours, COSMOS-Web represents the largest contiguous field surveyed during JWST Cycle 1, posing unique data reduction challenges due to its extensive scale. By combining the official JWST Calibration Pipeline with custom improvements for noise removal, background subtraction, and astrometric alignment, we achieve high fidelity science-ready mosaics. We detail the systematic approach employed in the three stages of the JWST Calibration Pipeline. The data, collected in three epochs from January 2023 to January 2024, encompass 152 visits and have been processed into 20 mosaic tiles to optimize computational efficiency and data processing. The final data products achieve 5 sigma depths of 26.7-28.3 AB mag in 0.15" apertures. The processed and calibrated datasets are made available to the public.

We study a sample of nearby (z~0.2) low-luminosity dwarf (10^7 MSun < M* < 10^8 MSun) galaxies which have extreme (0.1 - 3 MSun/yr) star formation rates (SFRs) for this mass regime, making them plausible analogues of dwarfs at z~5.5. We compare the properties of these analogues to control samples of 'normal' dwarfs, which reside on the star formation main sequence (SFMS) at z~0.2 and are matched in their stellar mass and redshift distributions to the analogue population. The analogue and normal populations do not show differences, either in their half-light radii or the projected distances to nodes, filaments and massive galaxies. This suggests that the comparatively extreme SFRs in the analogues are not driven by them being anomalously compact or because they reside in specific environments which might provide a larger gas supply. However, the fractions of interacting galaxies and those that have early-type morphology are significantly elevated (by factors of ~5.6 and ~9 respectively) in the analogues compared to the normal population. Extrapolation of the redshift evolution of the star formation main sequence into our mass range of interest appears to underestimate the SFRs of observed dwarfs at z~5.5. Since current SFMS measurements remain dominated by low and intermediate redshift data (especially at low stellar masses), our study suggests that this underestimation may be driven by interactions (which are more frequent at earlier epochs) boosting the SFRs in the high-redshift dwarf population. Our results are consistent with a picture where higher gas availability, augmented by interactions, drives much of the stellar mass assembly of dwarf galaxies in the early Universe.

We present a new 8.5 ks Chandra observation of Abell 1885, obtained as part of the Cluster Evolution Reference Ensemble At Low-z (CEREAL) survey of ~200 low-z galaxy groups and clusters. These data reveal that Abell 1885 is a strong cool core, with a central cooling time of 0.43 Gyr, and that the central galaxy hosts an X-ray luminous point source at its center (L=2.3x10^42 erg/s), indicative of a rapidly accreting supermassive black hole. In the context of the larger CEREAL sample, we constrain the fraction of clusters at z~0.15 with X-ray bright central AGN to be no more than 4.1%. Including radio data from LOFAR, GMRT, ASKAP, and the VLA and optical integral field unit data from SDSS MaNGA, we probe the details of cooling, feeding, and feedback in this system. These data reveal that cooling of the intracluster medium is highly suppressed on large (>10 kpc) scales despite a central supermassive black hole that is in the early stages of the self-regulation cycle (characterized by rapid accretion, physically small jets, and no large-scale low-frequency radio emission). To reconcile the large-scale quenching with a lack of visible large-scale feedback, we propose that the timescale on which energy is dissipated on large scales is significantly longer than the timescale on which black hole feeding operates on small (~pc) scales. These observations suggest that there are two separate timescales characterizing AGN feedback in clusters: the short timescale of small-scale feeding and feedback processes and a longer timescale by which energy is dissipated on large physical scales in the intracluster medium. This interpretation disfavors a model in which the energy is rapidly dissipated (e.g. shocks), which would synchronize the feeding and feedback timescales, and favors a model in which the heating effects of AGN feedback can linger long after the outburst has passed (e.g. turbulent mixing).

In this work, we present an equation of state formalism for hot Neutron Stars (NSs) which consistently includes the effects of finite temperature, hyperons as well as neutrino trapping, relevant for the study of proto-neutron stars, binary neutron star mergers and supernova explosions. Within a non-linear relativistic mean field description, the framework allows for a systematic variation of nuclear parameters within the range allowed by uncertainties in nuclear experimental data, ensuring compatibility with nuclear theory, astrophysical and heavy-ion data. We then investigate the role of nuclear and hypernuclear parameters on NS macroscopic properties and $f$-mode oscillations in hot neutron stars within Cowling approximation. Our results reveal that in hyperonic neutron stars with trapped neutrinos, the saturation nuclear density shows moderate to strong correlation with NS astrophysical observables. We also investigate the $C-Love$ universal relations for this model and provide fit relations for hot NS configurations in the neutrino-trapped regime.

The NASA Exoplanet Archive and the Exoplanet Follow-up Observing Program service are two widely used resources for the exoplanet community. The NASA Exoplanet Archive provides a complete and accurate accounting of exoplanetary systems published by NASA missions and by the community in the refereed literature. In anticipation of continued exponential growth in the number of exoplanetary systems, and the increasing complexity in our characterization of these systems, the NASA Exoplanet Archive has restructured its primary tables and interfaces, as well as extending and standardizing their modes of access. The Exoplanet Follow-up Observing Program service provides the exoplanet community with a venue for coordinating and sharing follow-up and precursor data for exoplanets, their host stars, and stars that might eventually be targets for future planet searches, and recently reached one million files uploaded by the community. In this paper we describe the updates to our data holdings, functionality, accessibility, and tools, as well as future priorities for these two services.

Energetic events on the Sun, particularly coronal mass ejections and high speed streams, regulate the near Earth space environment and give rise to space weather. A major terrestrial manifestation of such events are geomagnetic storms. A geomagnetic storm results in dissipation of energy from the solar wind into the atmosphere, leading to Joule heating and thermospheric expansion. This has serious consequences on Low Earth Orbit satellite lifetimes. Our work demonstrates the impact of different kinds of geomagnetic storms on satellite orbits. We also briefly discuss about some physical attributes of satellites that can make them prone to higher orbital decay. Our work highlights the importance of monitoring and predicting space weather, and assessing their impacts on space-based human technologies.

The COSMOS-Web survey is the largest JWST Cycle 1 General Observer program covering a contiguous ~0.54 deg$^2$ area with NIRCam imaging in four broad-band filters and a non-contiguous ~0.2 deg$^2$ with parallel MIRI imaging in a single broad-band filter, F770W. Here we present a comprehensive overview of the MIRI imaging observations, the data reduction procedure, the COSMOS-Web MIRI photometric catalog, and the first data release including the entire COSMOS-Web MIRI coverage. Data reduction is predominantly based on the JWST Science Calibration Pipeline with an additional step involving custom background subtraction to mitigate the presence of strong instrumental features and sky background in the MIRI images. We reach 5$\sigma$ (point source) limiting depths ($m_{F770W}$~25.51 based on $r$~0.3'' circular apertures) that are significantly better than initial expectations. We create a COSMOS-Web MIRI catalog based on the images presented in this release and compare the F770W flux densities with the Spitzer/IRAC CH4 measurements from the COSMOS2020 catalog for CH4 detections with S/N $>5$. We find that these are in reasonable agreement with a small median offset of $<0.05$ mag. We also derive robust 7.7$\mu$m number counts spanning five orders of magnitude in flux ($\sim$0.2-2300 $\mu$Jy) $\unicode{x2013}$ making COSMOS-Web the only JWST survey to date to efficiently sample such a large flux range $\unicode{x2013}$ which is in good agreement with estimates from other JWST and IRAC surveys.

We present dynamical mass measurements, obtained from multi-epoch Very Long Baseline Interferometry (VLBI) observations, for the young multiple stellar system EC\,95, located in the core of the Serpens star-forming region. Our dataset includes both archival data and more recent observations obtained as part of the \textit{Dynamical Masses of Young Stellar Multiple Systems with the VLBA project} (DYNAMO--VLBA), totalling 32 epochs over 12 years of observation. The observations span more than half of the EC\,95AB orbit, which has an estimated period of $21.6\pm0.1$ years, and yield masses of $2.15\pm0.10$ M$_\odot$ for EC\,95A and $2.00\pm0.12$ M$_\odot$ for EC\,95B. Furthermore, for the first time, we have estimated the mass of the third component (EC\,95C) using four available radio detections as well as an infrared observation. We find it to be 0.26 $^{+0.53}_{-0.46}$ M$_\odot$, suggesting that E\,95C is a low-mass T Tauri star. We conclude that EC\,95 is a hierarchical triple system comprised of a tight central binary proto-Herbig\,AeBe system and a lower mass companion on a wider orbit. These results contribute to improve our understanding of the stellar dynamics in this system and provide valuable insights into its triple nature.

We report on a pencil-beam survey of the Taurid Swarm, a possible concentration of bodies in the Taurid meteoroid stream associated with the 7:2 mean-motion resonance with Jupiter. Canada-France-Hawaii Telescope MegaCam observations reaching apparent magnitudes of 24.5 in the gri filter were taken over three nights. Rates of motion on the sky allowed for the quick elimination of main-belt objects from the over 1000 moving sources seen. Eight candidates with on-sky rates of motion consistent with Taurids were detected, but seven were subsequently shown to be non-Taurids (Hungarias, Mars-crossers, etc). One object might be a 60 m class Taurid but not enough data was collected and its orbit remains ambiguous. Our results are consistent with no Taurid Swarm members observed, and an upper limit of fewer than 3e3 - 3e4 objects down to H=25.6 +/- 0.3 (diameter of 34-76 m assuming a 2P/Encke-like albedo) at the 95% confidence level. While meteor observations have confirmed the Taurid Swarm's existence at meter and smaller sizes, our results indicate that the current mass budget of the swarm at 100 m sizes does not require an outsize parent to explain it.

The Hilbert-Huang transform (HHT) consists of empirical mode decomposition (EMD), which is a template-free method that represents the combination of different intrinsic modes on a time-frequency map (i.e., the Hilbert spectrum). The application of HHT involves introducing trials by imposing white noise on the signal and then calculating the ensemble mean process of the corresponding EMD to demonstrate its significance on the Hilbert spectrum. In this study, we develop a stacked Hilbert-Huang Transform (sHHT) method that generates the Hilbert spectrum for each trial and compiles all results to enhance the strength of the real instantaneous frequency of the main signal on the time-frequency map. This new approach is more sensitive to detecting/tracing the nonlinear and transient features of a signal embedded in astronomical databases than the conventional HHT, particularly when the signal experiences dramatic frequency changes in a short time. We analytically investigate the consistency of HHT and sHHT and perform numerical simulations to examine the dispersion of the instantaneous frequency obtained through sHHT and compare its advantages and effectiveness with those of conventional HHT. To confirm the feasibility of the sHHT, we demonstrate its application in verifying the signal of superorbital modulation in X-ray and binary black hole mergers in gravitational waves.

The presence of strong correlations between super-massive black hole (SMBH) masses and galaxy properties like stellar mass have been well-established in the local Universe, but how these scaling relations evolve with cosmic time is yet to be settled in both observations and theoretical models. Recent works have also highlighted the role of galaxy morphology on the scatter of the SMBH-galaxy mass scaling relations, while the impact of other galaxy properties remains poorly studied, like the role of galaxy environment. We use the state-of-the-art SHARK v2.0 semi-analytic model to explore the evolution of these galaxy-SMBH scaling relations to expand the available predictions from theoretical models to contrast with existing and upcoming observations. We find the relations between SMBH masses and both total and bulge stellar mass predicted by SHARK v2.0 to be in good overall agreement with observational measurements across a wide range of redshift and stellar masses. These scaling relations show a significant evolution as a function of cosmic time in SHARK v2.0, with SMBH masses $\sim1$ dex lower at $z=0$ compared to $z=9$ at fixed stellar mass and the scatter increasing by a factor of $\sim2-5$ towards low redshift. Both relations show a strong dependence with galaxy morphology and the main source for SMBH growth (gas accretion or mergers), with smaller trends with star formation rate, galaxy sizes, and environment. We find that galaxy morphology alone explains most of the scatter around both scaling relations, with other galaxy properties tying to the SMBH scaling relations through their correlations with morphology.

Future space-based far infra-red astronomical observations require background limited detector sensitivities and scalable focal plane array solutions to realise their vast potential in observation speed. In this work, a focal plane array of lens absorber coupled Kinetic Inductance Detectors (KIDs) is proposed to fill this role. The figures of merit and design guidelines for the proposed detector concept are derived by employing a previously developed electromagnetic spectral modelling technique. Two designs operating at central frequencies of $6.98$ and $12$ THz are studied. A prototype array of the former is fabricated, and its performance is experimentally determined and validated. Specifically, the optical coupling of the detectors to incoherent distributed sources (i.e. normalised throughput) is quantified experimentally with good agreement with the estimations provided by the model. The coupling of the lens absorber prototypes to an incident plane wave, i.e. aperture efficiency, is also indirectly validated experimentally matching the expected value of $54\%$ averaged over two polarisation. The noise equivalent power of the KIDs are also measured with limiting value of $8\times10^{-20}$ $\mathrm{W/\sqrt{Hz}}$.

Galaxy-galaxy strong gravitational lenses can constrain dark matter models and the Lambda Cold Dark Matter cosmological paradigm at sub-galactic scales. Currently, there is a dearth of images of these rare systems with high signal-to-noise and angular resolution. The Nancy Grace Roman Space Telescope (hereafter, Roman), scheduled for launch in late 2026, will play a transformative role in strong lensing science with its planned wide-field surveys. With its remarkable 0.281 square degree field of view and diffraction-limited angular resolution of ~0.1 arcsec, Roman is uniquely suited to characterizing dark matter substructure from a robust population of strong lenses. We present a yield simulation of detectable strong lenses in Roman's planned High Latitude Wide Area Survey (HLWAS). We simulate a population of galaxy-galaxy strong lenses across cosmic time with Cold Dark Matter subhalo populations, select those detectable in the HLWAS, and generate simulated images accounting for realistic Wide Field Instrument detector effects. For a fiducial case of single 146-second exposures, we predict around 160,000 detectable strong lenses in the HLWAS, of which about 500 will have sufficient signal-to-noise to be amenable to detailed substructure characterization. We investigate the effect of the variation of the point-spread function across Roman's field of view on detecting individual subhalos and the suppression of the subhalo mass function at low masses. Our simulation products are available to support strong lens science with Roman, such as training neural networks and validating dark matter substructure analysis pipelines.

Exo-Kuiper belts have been observed for decades, but the recent detection of gas in some of them may change our view of the Solar System's youth. Late gas produced by the sublimation of CO (or CO$_2$) ices after the dissipation of the primordial gas could be the norm in young planetesimal belts. Hence, a gas-rich Kuiper belt could have been present in the Solar System. The high C/H ratios observed on Uranus and Neptune could be a clue to the existence of such late gas that could have been accreted onto young icy giants. The aim of this paper is to estimate the carbon enrichment of the atmospheres of Uranus and Neptune caused by the accretion of the gas released from a putative gas-rich Kuiper belt. We find that assuming a primordial Kuiper belt with a mass of tens of earth masses leads to significant CO gas accretion onto the giants, which can lead to high C/H ratios, especially for Uranus and Neptune. Our model shows that a relatively massive gas-rich Kuiper belt could have existed in the Solar System's youth, significantly enriching the atmospheres of Uranus and Neptune with carbon. Late gas accretion and its effect on outer giant planets metallicities could be a universal scenario, also occurring in extrasolar systems. Observations of sub-Jupiter exoplanets could provide very useful information to better constrain this scenario, with an enrichment in carbon and oxygen (for warm-enough planets) compared to other elements that should be inversely proportional to their envelope mass.

With over 1,000,000 images from more than 10,000 exposures using state-of-the-art high-contrast imagers (e.g., Gemini Planet Imager, VLT/SPHERE) in the search for exoplanets, can artificial intelligence (AI) serve as a transformative tool in imaging Earth-like exoplanets in the coming decade? In this paper, we introduce a benchmark and explore this question from a polarimetric image representation learning perspective. Despite extensive investments over the past decade, only a few new exoplanets have been directly imaged. Existing imaging approaches rely heavily on labor-intensive labeling of reference stars, which serve as background to extract circumstellar objects (disks or exoplanets) around target stars. With our POLARIS (POlarized Light dAta for total intensity Representation learning of direct Imaging of exoplanetary Systems) dataset, we classify reference star and circumstellar disk images using the full public SPHERE/IRDIS polarized-light archive since 2014, requiring less than 10 percent manual labeling. We evaluate a range of models including statistical, generative, and large vision-language models and provide baseline performance. We also propose an unsupervised generative representation learning framework that integrates these models, achieving superior performance and enhanced representational power. To our knowledge, this is the first uniformly reduced, high-quality exoplanet imaging dataset, rare in astrophysics and machine learning. By releasing this dataset and baselines, we aim to equip astrophysicists with new tools and engage data scientists in advancing direct exoplanet imaging, catalyzing major interdisciplinary breakthroughs.

Fast radio bursts (FRBs) are emerging as powerful cosmological probes for constraining the baryon fraction in the intergalactic medium (IGM), offering a promising avenue to address the missing baryon problem. In this paper, we analyze constraints on the IGM baryon fraction ($f_\mathrm{IGM}$) using 92 localized FRBs, incorporating a corrected probability distribution function for the IGM dispersion measure within three different cosmological models. We find that variations in the underlying cosmological model have a negligible impact on the inferred values of $f_\mathrm{IGM}$. While the NE2001 Galactic electron density model yields slightly higher $f_\mathrm{IGM}$ values compared to the YMW16 model, the results are consistent within the 1$\sigma$ confidence level. Additionally, there is no statistically significant evidence for redshift evolution in $f_\mathrm{IGM}$. Our analysis constrains $f_\mathrm{IGM}$ to the range $0.8 \sim 0.9$, providing strong support for the idea in which the majority of the missing baryons reside in the diffuse IGM.

Significant variability in broad emission line strengths of active galactic nuclei (AGN) over months to years has been observed, often accompanied by intrinsic continuum changes. Such spectral variability challenges the traditional AGN classification scheme, which attributes differences between Type 1 and Type 2 to geometrical effects, as transitions between these types occur on timescales shorter than viscous ones. In this work, using the {\sc cloudy} photo-ionization simulations, we investigated the response of the major emission line fluxes, in the optical/UV and hard X-ray bands, to changes in the intensity and shape of the continuum emission of the AGN under two scenarios: (i) changes in the X-ray power-law while keeping disc emission fixed, and (ii) broadband continuum variations. We demonstrate that BLR line fluxes are insensitive to X-ray power-law changes alone. Considering a well-studied case of the changing-look (CL) AGN Mrk 1018, which exhibits variations in the intrinsic disc emission, as well as the X-ray power-law, our simulations reproduce observed brightening and dimming trends of the BLR emission. Moreover, we show that the highly ionized Fe K$\alpha$ X-ray flux, primarily produced by the H-like and He-like ions of Fe, strongly depends on the X-ray strength of the intrinsic SED. These findings suggest that the origin of highly ionized Fe K$\alpha$ emission is in the coronal part of the accretion disk and that the CL phenomenon can be triggered by intrinsic changes in the accretion properties of AGN.

The low-frequency radio telescope of the Square Kilometre Array (SKA-Low), currently under construction in the remote Murchison shire in the Western Australia's outback, will observe the sky between 50 MHz and 350 MHz with unprecedented sensitivity and stringent requirements for polarization accuracy. In this work, we investigate the instrumental polarization purity of a SKA-Low prototype station by means of the Intrinsic Cross-Polarization Ratio (IXR) figure of merit. We derive all-sky experimental IXR maps using data from the Aperture Array Verification System 2 (AAVS2). The results are presented at three frequencies within the SKA-Low bandwidth (110, 160, and 230 MHz) with a quantitative comparison between observed and simulated all-sky IXR maps. Our findings show good agreement in IXR map distributions and promising consistency in their radial profiles, meeting SKA-Low's IXR specification overall. This study offers an empirical approach to verifying SKA-Low's polarization performance using all-sky observations from individual stations and will potentially support the telescope's early science commissioning phase.

The origin of fast radio bursts (FRBs), highly energetic, millisecond-duration radio pulses originating from beyond our galaxy, remains unknown. Observationally, FRBs are classified as non-repeating or repeating, however, this classification is complicated by limited observing time and sensitivity constraints, which may result in some repeating FRBs being misidentified as non-repeating. To address this issue, we adopt both empirical and machine-learning techniques from previous studies to identify candidates that may have been misclassified. We conducted follow-up observations of 36 such candidates, each observed for 10 minutes using the Five-hundred-meter Aperture Spherical Telescope (FAST). No radio bursts exceeding a signal-to-noise ratio of 7 were detected, with a typical 7 sigma fluence limit of ~0.013 Jy ms. We constrain the repetition rates of these sources using two statistical models of FRB occurrence. Combining our FAST non-detections with prior observations, we derive upper limits on the repetition rates of ~$10^{-2.6}$-$10^{-0.22}$ hr$^{-1}$ under a Poisson process, and ~$10^{-2.3}$-$10^{-0.25}$ hr$^{-1}$ under a Weibull process. This work presents one of the most stringent upper limits on FRB repetition rates to date, based on a sample size five times larger than those used in previous studies.

Core-collapse supernovae (CCSNe) are regularly observed electromagnetically, prompting targetted searches for their gravitational-wave emission. However, there are scenarios where these powerful explosions may not have any observable electromagnetic signal, but would still have strong detectable emission in gravitational waves and neutrinos. A regular CCSN explosion may be obscured by matter in the Galaxy. A star may undergo a failed CCSN explosion, where the stalled shockwave is not revived, and would eventually form a black hole. Higher mass progenitor stars may revive the shock, but form a black hole too quickly for the shockwave to reach the surface of the star and produce an electromagnetic signal. Previous work has shown that we can determine if a black hole forms from the CCSN neutrino emission if there are long duration sinusoidal modulations in the neutrino signal caused by the standing accretion shock instability (SASI). The SASI also produces an observable signature in the gravitational-wave emission. In this paper, we investigate if we can distinguish between different scenarios for electromagnetically dark CCSNe using the gravitational-wave emission alone. We find, using a reconstruction of the SASI mode, abrupt end times of the gravitational-wave emission, and the rate of change of frequency of the dominant mode, that we are able to accurately distinguish between an obscured CCSN, a failed CCSN, and an explosion with fast black hole formation.

Evidence for a low-frequency gravitational-wave background using pulsar timing arrays has generated recent interest into its underlying contributing sources. However, multiple investigations have seen that the significance of the evidence does not change with choice of pulsar modeling techniques but the resulting parameters from the gravitational wave searches do. PSR J1455-3330 is one of the longest-observed pulsars in the array monitored by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) but showed no evidence for long-timescale red noise, either intrinsic or the common signal found among many pulsars in the array. In this work, we argue that NANOGrav's piecewise-constant function used to model variations in radio-frequency-dependent dispersive delay should not be used for this pulsar, and a much simpler physical model of a fixed solar wind density plus a linear trend in dispersion measure is preferred. When the original model is replaced, (i) the pulsar's timing parallax signal changes from an upper limit to a significant detection, (ii) red noise becomes significant, and (iii) the red noise is consistent with the common signal found for the other pulsars. Neither of these signals are radio-frequency dependent. While the same physical motivation will not apply to many of the pulsars currently used in pulsar timing arrays, we argue for careful physically-motivated timing and noise modeling of pulsars used in precision timing experiments.

Protoplanetary disks observed in scattered light reveal essential insights into the disk's three-dimensional architecture and dust properties. These disks, which play a crucial role in planet formation, have complex structures where the visibility of the disk's backside can vary significantly based on several parameters. This study aims to explore the factors impacting backside visibility in protoplanetary disks, particularly under variations in inclination, dust distribution, grain characteristics, and outer disk morphology. Using RADMC-3D radiative transfer simulations, we investigate how these variables influence the appearance of the backside in scattered light images. Tapered disk models with exponential tapers, frequently obscure the backside, which supports the rarity of observed backside features. In cases where backside features are visible at lower inclinations, they likely indicate cut-off disks, as backside detection is challenging in standard tapered models at these inclinations. Additionally, factors like dust mass, grain distribution, and disk material stratification play crucial roles in backside observability, affecting its potential detection in real observations. This study contributes to understanding the detectability of the backside in protoplanetary disks, with implications for refining observational strategies and interpreting backside features in scattered light images. These findings help frame backside visibility as a critical aspect of assessing disk structure and evolution.

Several phenomenological aspects of low-luminosity neutron star transients, such as atolls, remain poorly understood. One such source, MAXI J1807+132, entered its latest outburst in July 2023. To thoroughly characterize this outburst, we conducted an extensive observational campaign spanning radio to X-ray wavelengths. Here, we present the results of this campaign, which covered the period from before the outburst to the return to quiescence. We detected a delay between the X-ray and optical rise times, which is consistent with the predictions of the disk instability model with a truncated disk. The color evolution and optical/X-ray correlations, along with infrared and radio detections, support the presence of jet synchrotron emission during the gradual decay phase following the peak. We also report for the first time in an X-ray binary a near-orthogonal rotation of the optical polarization just before a small flare, after which the jet is thought to be quenched. The main outburst is followed by several high-amplitude, rapid reflares in the optical, ultraviolet, and X-ray bands, the origin of which remains difficult to constrain.

We present a catalog of 142 dark galaxy candidates in a region covered by the Arecibo Legacy Fast ALFA (ALFALFA) survey. We start with 344 ALFALFA HI sources without optical counterparts and remove those that do not seem to have dark galaxy origin. To do that, we first eliminate 83 sources that are known HI clouds probably formed from tidal interactions between galaxies and 13 sources that have optical counterparts. We then remove 56 sources located near other HI sources, which are likely to be HI clouds. We further exclude 10 sources that have nearby HI sources within the ALFALFA beam and 40 sources potentially associated with nearby galaxies. We perform visual inspection of optical images from DESI Legacy Imaging Survey with an $r$-band surface brightness limit of $\sim 28.5 \rm \ mag\ arcsec^{-2}$ as well as NUV images from GALEX to confirm the absence of stellar emission. We additionally inspect infrared images from WISE and AKARI for dust emission. As a result, we identify 142 dark galaxy candidates and analyze their physical properties by comparing with luminous galaxies. We find that the dark galaxy candidates generally have smaller dynamical masses, higher HI-to-dynamical mass ratios, and are located in less dense regions when compared to luminous galaxies, which is consistent with results from cosmological simulations. This sample provides an important testbed for studying the role of dark matter in galaxy formation and evolution.

The masses of white dwarfs (WDs) in intermediate polars (IPs) can be determined from the shape of their hard X-ray spectra. Here we study the importance of all possible systematic uncertainties in this X-ray spectroscopy method, including finite radii and rotation of magnetospheres, finite accretion column height and accretion-flow inclination relative to the WD surface. We also investigate the importance of accretion-heated envelopes on WD surfaces in IPs which are increasing WD radii. Their presence changes the commonly used mass-radius relation for cold white dwarfs. As a first approximation we use thick ($10^{-4}M_\odot$) hydrogen envelope models with a surface temperature of 30 kK. We present a new model grid of hard X-ray spectra of high-luminous IPs computed among other things with using a new mass-radius relation. This grid is used for fitting Swift/BAT spectra of 47 IPs. The average WD mass in this sample is 0.82 $M_\odot$ and coincides with the average WD mass in cataclysmic variables obtained by optical methods. This means that the calculated hard X-ray spectra and the assumptions made that the magnetospheric radii in IPs are close to the corotation radii, and the relative heights of the accretion columns are small are basically correct, because most IPs have high luminosities. But this universal grid (as well as previous universal grids) cannot give correct results for the low-luminous IPs with probably relatively tall accretion columns on the WD surfaces. Such IPs have to be investigated with individual accretion column models.

We present in this proceeding the results of the test phase of the GRAVITY+ adaptive optics. This extreme AO will enable both high-dynamic range observations of faint companions (including exoplanets) thanks to a 40x40 sub-apertures wavefront control, and sensitive observations (including AGNs) thanks to the addition of a laser guide star to each UT of the VLT. This leap forward is made thanks to a mostly automated setup of the AO, including calibration of the NCPAs, that we tested in Europe on the UT+atmosphere simulator we built in Nice. We managed to reproduce in laboratory the expected performances of all the modes of the AO, including under non-optimal atmospheric or telescope alignment conditions, giving us the green light to proceed with the Assembly, Integration and Verification phase in Paranal.

We present a comprehensive spectro-polarimetric and timing analysis of eleven black hole X-ray binaries, namely Cyg X-1, 4U 1630-47, Cyg X-3, LMC X-1, 4U 1957+115, LMC X-3, Swift J1727.8-1613, GX 339-4, Swift J151857.0-572147, IGR J17091-3624, and MAXI J1744-294, using quasi-simultaneous observations from {\it IXPE}, {\it NICER}, {\it NuSTAR}, and {\it AstroSat}. Timing analyses reveal characteristic type-B and type-C Quasi-periodic Oscillations (QPOs) across different spectral states, often associated with episodic radio ejections. Broadband ($0.7-60$ keV) spectral modelling, employing disc, Comptonization, and reflection components, reveals degeneracies in constraining disc-corona geometries. Polarimetric measurements in 2-8 keV band detect significant polarization degrees (PDs) ranging from $3-20.6\%$ ($1.2-21.4\%$) in harder (softer) states, with moderate to strong energy dependence, except for LMC X-1, Swift J151857.0-572147, and MAXI J1744-294, where no significant polarization is detected. We report the polarization detections of Cyg X-3 (PD $\sim 21.4\%$, SIMS), LMC X-3 (PD $\sim 2.4\%$, HSS) and IGR J17091-3624 (PD $\sim 9\%$, LHS) using the recent {\it IXPE} observations. A positive correlation is found between PD and the Comptonized photon fraction ($cov_{\rm frac}$), while an anti-correlation is observed with the disc-to-Comptonized flux ratio ($F_{\rm ratio}$) across spectral states. The combined timing, spectral, and polarimetric results, together with constraints from radio jet observations, suggest a radially extended corona within a truncated disc for Cyg X-1, Swift J1727.8-1613, IGR J17091-3624, and GX 339-4, whereas the disc-corona geometry remains unconstrained for 4U 1957+115, LMC X-3, and 4U 1630-47. We discuss the implications of these findings for understanding accretion geometries and highlight prospects for future X-ray polarimetric studies.

We have developed a novel methodology utilizing molecular dynamics simulations to evaluate the ionization yields of nuclear recoils in crystalline silicon. This approach enables analytical exploration of atomic-scale transport within the lattice without necessitating parameterization. The quenching factors across the nuclear recoil energy range from 20 eV to 10 keV have been thoroughly investigated. A remarkable agreement with experimental data is achieved, particularly for the minimal energy regime conducted to date, reaching the level of a single electron-hole pair. This work presents evidence of a crucial and fundamental distribution of the quenching factor, which can be associated to the collisional interactions underlying the transport phenomena. The region below 4 keV of the quenching factor, where discrepancies have been observed with the Lindhard's model, is found to be significantly attributed to the lattice binding effect and the specific crystal structure. In contrast, a gradual functional relationship is identified below approximately 100 eV, indicating that the quenching factor is influenced by the crystallographic orientation of the target material. From a distributional perspective, our analysis allows for the determination of the minimum exclusion mass for the dark matter nucleon elastic scattering channel at 0.29 $\mathrm{GeV}/c^2$, thereby significantly enhancing sensitivity for the sub-$\mathrm{GeV}/c^2$ mass region.

Quasi-periodic pulsations (QPPs) are a common feature of impulsive solar flare emissions, yet their driving mechanism(s) remain unresolved. We study long-period (3.4 min) QPPs in the M3.7 flare on Feb 24, 2023, using SDO/AIA and Solar Orbiter/STIX to pinpoint their source and driver. We analyzed UV and HXR emissions from four subregions covering the flare ribbons and footpoints of the erupting filament. We find the QPP characteristics varied significantly across these regions. The strongest UV pulsations, which correlated strongly with HXR emissions, originated from a sequence of flare kernels at the expanding front of the southern flare ribbon. In contrast, another expanding ribbon front showed no pulsations. The spatial coincidence of UV and HXR sources points to non-thermal electrons as a common driver, with the HXR spectrum showing a soft-hard-soft evolution matching the pulsation timescale. We also observed UV pulsations from plasma injections into a separate filament channel, distinct from the main flare kernels. Our results suggest that the spatiotemporal evolution of the pulsations is driven by the propagation of magnetic reconnection in an asymmetric geometry. We propose that slipping reconnection along quasi-separatrix layers is key to generating and moving the UV kernels, potentially modulated by another time-varying process.

Several cosmological observations (e.g., Cosmic Microwave Background (CMB), Supernovae Type Ia, and local distance ladder measurements such as Cepheids) have been used to measure the global expansion rate of the Universe, i.e., the Hubble constant, $H_{0}$. However, these precision measurements have revealed tensions between different probes that are proving difficult to solve. Independent, robust techniques must be exploited to validate results or mitigate systematic effects. We use the Cosmic Chronometer (CC) method, which leverages the differential age evolution of passive galaxies, to measure $H(z)$, without any assumption of the underlying cosmology. Unlike previous CC studies, we used only brightest cluster galaxies (BCGs), the oldest and most massive galaxies in the Universe, to construct a pure and homogeneous sample. In this work we used a sample of 53 BCGs in massive, Sunyaev-Zel'dovich selected galaxy clusters (0.3 $< z <$ 0.7) with Southern African Large Telescope (SALT) spectroscopic observations. We used optical spectra to measure D4000$_{\rm n}$ of the BCGs to obtain a new direct measurement of $H(z) = 72.1 \pm 33.9(\rm stat) \pm 7.3$(syst) km s$^{-1}$ Mpc$^{-1}$ at $z=0.5$. By using BCGs, we significantly reduced the systematic errors to 10% by minimising the stellar mass and metallicity dependence of the method. The dominant uncertainty, and limitation for our study, is statistical, and we need larger, homogeneous samples of the oldest, most massive galaxies. By using the $Planck$+BAO prior of $\Omega_{m}$ and $\Omega_{\Lambda}$, the projected Hubble constant is $H_{0}$ = $54.6 \pm 25.7(\rm stat) \pm 5.5$(syst) km s$^{-1}$ Mpc$^{-1}$, consistent with both CMB and Cepheid measurements.

Relativistic jets from AGN are an important driver of feedback in galaxies. They interact with their environments over a wide range of physical scales during their lifetime, and an understanding of these interactions is crucial for unraveling the role of supermassive blackholes in shaping galaxy evolution. The impact of such jets have been traditionally considered in the context of heating the large-scale environments. However, in the last few decades there has been additional focus on the immediate impact of jet feedback on the host galaxy itself. In this review we outline the development of various numerical simulations since the onset of studies of jets to the present day, where sophisticated numerical techniques have been employed to study jet feedback including a range of physical processes. The jets can act as an important agent of injecting energy in the host's ISM, as confirmed both in observations of multi-phase gas, as well as in simulations. Such interactions have the potential to impact the kinematics of the gas as well as its star formation. We summarize the recent results from simulations of jet feedback on kpc scales, and outline the broader implications for observations and galaxy evolution.

L-BASS is an instrument designed to make radiometric temperature measurements of the sky with an absolute accuracy of better than 0.1 K at 1.4 GHz. This will be achieved in two steps: first by measuring the sky temperature relative to that of the North Celestial Pole, using two horn-based antennas, and second with the sky antenna replaced with a calibrated cryogenic load to measure the absolute brightness temperature of the North Celesial Pole. Here we describe the design of the L-BASS two-antenna system and report on laboratory measurements to establish its performance at component and sub-system level.

We use analytical and $N$-body methods to study the evolution of dwarf spheroidal galaxies (dSphs) embedded in dark matter (DM) haloes that host a sizeable subhalo population. Dark subhaloes generate a fluctuating gravitational field that injects energy into stellar orbits, driving a gradual expansion of dSphs. Despite the overall expansion, the stellar density profile preserves its initial shape, suggesting that the evolution proceeds in a self-similar manner. Meanwhile, the velocity dispersion profile, initially flat, evolves as the galaxy expands: the inner regions heat up, while the outer regions cool down. Kinematically, this resembles gravothermal collapse but with an inverted evolution, instead of collapsing the stellar system expands within a fluctuating halo potential. As the half-light radius $r_{\rm half}$ approaches the halo peak velocity radius $r_{\rm max}$, the expansion slows, and the velocity dispersion peaks at $\sigma_{\rm max} \simeq 0.27 v_{\rm max}$. The stellar heat capacity remains positive for deeply embedded stars but diverges near $r_{\rm max}$, turning negative beyond this threshold, which indicates a phase transition in the dynamical response to energy injection. The relaxation time scales as $t_{\rm rel} \sim r_{\rm half}^{3/2}$, showing that orbital diffusion slows as the galaxy expands. Ultra-faint dSphs, having the smallest sizes and shortest relaxation times, are particularly sensitive to the presence of dark subhaloes. Some of our dSph models expand beyond the detection of current photometric surveys, becoming `stealth' galaxies with half-light radii and velocity dispersions similar to ultra-diffuse galaxies (UDGs), but orders of magnitude lower luminosities, representing a distinct, yet currently undetected, population of DM-dominated satellites.

As an alternative gravitational theory to General Relativity (GR), the Conformal Gravity (CG) has recently been successfully verified by observations of Type Ia supernovae (SN Ia) and the rotation curves of spiral galaxies. The observations of galaxies only pertain to the non-relativistic form of gravity. In this context, within the framework of the Newtonian theory of gravity (the non-relativistic form of GR), dark matter is postulated to account for the observations. On the other hand, the non-relativistic form of CG predicts an additional potential: besides the Newtonian potential, there is a so-called linear potential term, characterized by the parameter $\gamma^*$, as an alternative to dark matter in Newtonian gravity. To test CG in its non-relativistic form, much work has been done by fitting the predictions to the observations of circular velocity (rotation curves) for spiral galaxies. In this paper, we test CG with the observations from elliptical galaxies. Instead of the circular velocities for spiral galaxies, we use the velocity dispersion for elliptical galaxies. By replacing the Newtonian potential with that predicted by non-relativistic form of CG in Hamiltonian, we directly extend the Jeans equation derived in Newtonian theory to that for CG. By comparing the results derived from the ellipticals with that from spirals, we find that the extra potential predicted by CG is not sufficient to account for the observations of ellipticals. Furthermore, we discover a strong correlation between $\gamma^*$ and the stellar mass $M^*$ in dwarf spheroidal galaxies. This finding implies that the variation in $\gamma^*$ violates a fundamental prediction of Conformal Gravity (CG), which posits that $\gamma^*$ should be a universal constant.

It was the best of times, it was the worst of times, . . . - for the thermal/suprathermal particle populations in the vicinity of two traveling interplanetary shocks observed by Solar Orbiter on 2023-11-29 07:51:17 UTC and 2023-11-30 10:47:26 UTC at $\sim 0.83$ astronomical units from the Sun. We investigate these two very dissimilar shocks and elucidate their non-equilibrium features. We do not provide explanations of all observed features, our aim is to report them here for future reference. We use high-resolution data obtained with Solar Orbiter's Energetic Particle Detector (EPD), magnetometer (MAG), and Solar Wind Analyzer (SWA) to exhibit the very different natures of these two shocks and describe the detailed properties of suprathermal and energetic particles in their vicinity. We observe very different behavior of the energetic particle population because the two shocks are quite different. Solar wind protons and $\alpha$-particles are highly dynamic at the first, their beams appear to align well with rapid oscillations of the interplanetary magnetic field. Suprathermal particles associated with the second shock exhibit clear non-equilibrium and anisotropic features in their differential intensities at time scales comparable to the proton gyroperiod. The different geometries of the two shocks resulted in highly dissimilar populations of suprathermal and energetic particles in their vicinity. The first shock was associated with very interesting microphysics of the bulk plasma velocity distribution, the second resulted in similarly interesting microphysics of the accelerated particles. Both showed strong temporal variability of the particle populations at scales comparable to the proton gyroperiod.

A comprehensive analysis of the cosmological large-scale structure derived from galaxy surveys involves field-level inference, which requires a forward modelling framework that simultaneously accounts for structure formation and tracer bias. While structure formation models are well-understood, the development of an effective field-level bias model remains imprecise, particularly in the context of tracer perturbation theory within Bayesian reconstruction methods, which we address in this work. To bridge this gap, we have developed a differentiable model that integrates augmented Lagrangian perturbation theory, nonlinear, nonlocal, and stochastic biasing. At the core of our approach is the Hierarchical Cosmic-Web Biasing Nonlocal (HICOBIAN) model, which provides a positive definite description of tracer bias while incorporating a long- and short-range nonlocal framework via cosmic-web regions and deviations from Poissonity in the likelihood. A key insight of our model is that transitions between cosmic-web regions are inherently smooth, which we implement using sigmoid-based gradient operations. This enables a fuzzy, and, hence, differentiable hierarchical cosmic-web description, making the model well-suited for machine learning frameworks. Our approach accurately reproduces the primordial density field within associated error bars, as validated by two- and three-point statistics in Fourier space. Furthermore, we demonstrate that the methodology approaches the maximum encoded information consistent with Poisson noise. We introduce a Bayesian field-level inference algorithm that leverages the same forward modelling framework used in galaxy, quasar, and Lyman alpha forest mock catalog generation -- including nonlinear, nonlocal and stochastic bias with redshift space distortions -- providing a unified and consistent approach to the analysis of large-scale cosmic structure.

We report the discovery of a relativistic jet in Mrk~110, a narrow-line Seyfert 1 galaxy historically classified as a radio-quiet active galactic nucleus (AGN). Very Long Baseline Interferometry (VLBI) observations reveal intermittent jet activity during 2015--2016 and 2022--2024, with proper motion measurements yielding superluminal velocities of $\sim3.6\pm0.6\,c$ and $\sim2.1\pm0.2\,c$, respectively. The recent jet component decelerates to $\sim 1.5\pm0.2\,c$ at a projected distance of 1.1 parsec from the core, coinciding with the transition zone between broad-line and narrow-line regions. This deceleration accompanies dramatic spectral evolution from steep (the spectral index $\alpha \approx -0.63 \pm 0.04$) to inverted ($\alpha \approx +0.69 \pm 0.10$) as the 7.6 GHz flux density more than doubled. These episodic jet ejections and their evolutionary pattern match theoretical predictions from magnetically arrested disk (MAD) models for temporary jet formation in systems with Mrk 110's physical parameters on timescales of months to years. The observed jet deceleration distance matches expectations for relativistic outflows interacting with the circumnuclear environment. These findings demonstrate that the traditional radio-loud/quiet AGN dichotomy may reflect time-averaged states rather than intrinsic capabilities, suggesting that jets may form across the AGN population but become observable only during specific accretion phases when MAD conditions are temporarily established. Mrk 110 serves as a critical "missing link" between radio-loud and radio-quiet AGN, providing insight into jet formation mechanisms, environmental interactions, and physical processes that unify various AGN classifications.

Observations with the James Webb Space Telescope (JWST) reveal a previously unseen population of compact red objects, known as ``little red dots`` (LRDs). We study a new photometrically selected sample of 124 LRDs in the redshift range $z$ $\sim$ 3 - 10 selected from NIRCam coverage of the CEERS, NEP-TDF, JADES and JEMS surveys. For JADES, the NEP-TDF and CEERS, we compare SED models with and without AGN components and analyse the impact of an AGN component on the goodness of fit using the Bayesian information criterion (BIC). We find that whilst the $\chi^{2}$ of the majority of models containing AGN components is improved compared to models without AGN components, we show that the BIC suggests models without AGN are a more appropriate fit to LRD SEDs, especially when MIRI data is available. We also measure LRD clustering in the CEERS field, JADES field, and NEP-TDF, where we compare the spatial distribution of LRDs and galaxies with Kolmogorov-Smirnov tests of equality of distribution. We find that the neighbourhood of LRDs tends to be less dense compared to galaxies at all selections and masses and at similar redshifts. We further measure upper limit estimates for the halo masses of LRDs using abundance matching. Whilst the population of LRDs could be a mixture of several different inherent populations, as a whole it does appear that these systems are mostly hosting compact galaxies or star clusters in formation.

The timing analysis of transient events allows for investigating numerous still open areas of modern astrophysics. The article explores all the mathematical and physical tools required to estimate delays and associated errors between two Times of Arrival (ToA) lists, by exploiting Cross Correlation Function (CCF) techniques. The CCF permits the establishment of the delay between two observed signals and is defined on two continuous functions. A detector does not directly measure the intensity of the electromagnetic signal (interacting with its material) but rather detects each photon ToA through a probabilistic process. Since the CCF is defined on continuous functions, the crucial step is to obtain a continuous rate curve from a list of ToA. This step is treated in the article and the constructed rate functions are light curves that are continuous functions. This allows, in principle, the estimation of delays with any desired resolution. Due to the statistical nature of the measurement process, two independent detections of the same signal yield different photon times. Consequently, light curves derived from these lists differ due to Poisson fluctuations, leading the CCF between them to fluctuate around the true theoretical delay. This article describes a Monte Carlo technique that enables reliable delay estimation by providing a robust measure of the uncertainties induced by Poissonian fluctuations. GRB data are considered as they offer optimal test cases for the proposed techniques. The developed techniques provides a significant computational advantage and are useful analyzing of data characterized by low-count statistics (i.e., low photon count rates in c/s), as they allow overcoming the limitations associated with traditional fixed bin-size methods.

We present a comprehensive study of the structural evolution of Brightest Group Galaxies (BGGs) from redshift $z \simeq 0.08$ to $z = 3.7$ using the \textit{James Webb Space Telescope}'s 255h COSMOS-Web program. This survey provides deep NIRCam imaging in four filters (F115W, F150W, F277W, F444W) across $\sim 0.54~\mathrm{deg}^2$ and MIRI coverage in $\sim 0.2~\mathrm{deg}^2$ of the COSMOS field. High-resolution NIRCam imaging enables robust size and morphological measurements, while multiwavelength photometry yields stellar masses, SFRs, and Sérsic parameters. We classify BGGs as star-forming and quiescent using both rest-frame NUV--$r$--$J$ colors and a redshift-dependent specific star formation rate (sSFR) threshold. Our analysis reveals: (1) quiescent BGGs are systematically more compact than their star-forming counterparts and exhibit steeper size--mass slopes; (2) effective radii evolve as $R_e \propto (1+z)^{-\alpha}$, with $\alpha = 1.11 \pm 0.07$ (star-forming) and $1.40 \pm 0.09$ (quiescent); (3) star formation surface density ($\Sigma_{\mathrm{SFR}}$) increases with redshift and shows stronger evolution for massive BGGs ($\log_{10}(M_\ast/M_\odot) \geq 10.75$); (4) in the $\Sigma_e$--sSFR plane, a structural transition marks the quenching process, with bulge-dominated systems comprising over 80\% of the quiescent population. These results highlight the co-evolution of structure and star formation in BGGs, shaped by both internal and environmental processes, and establish BGGs as critical laboratories for studying the baryonic assembly and morphological transformation of central galaxies in group-scale halos.

Context: Europium (Eu) serves as a crucial tracer to understand the origin of rapid neutron-capture process (r-process) elements. An extensive effort was made in the last decade to model the chemical evolution of this element in the Galaxy. However, less attention was reserved to Eu in different galaxies of the Local Group (LG). Aims: By employing detailed and well-tested chemical evolution models, we investigate Eu enrichment across LG dwarf spheroidal galaxies, allowing for a direct comparison between model predictions for dwarf galaxies and the Milky Way (MW). Methods: Building upon an r-process enrichment framework that successfully reproduces the observed Eu abundance patterns as well as the supernova and compact binary merger rates in the MW, we build chemical evolution models for the Sagittarius, Fornax, and Sculptor dwarf spheroidal galaxies and test the enrichment scenario against the abundance patterns observed in these galaxies. Results: Models reproducing the Galactic Eu patterns significantly underestimate the [Eu/Fe] ratios observed in LG dwarfs. To address this "missing Eu" problem, we estimate the Eu production rate needed to match the observations and explore potential contributions either from prompt or delayed sources, assessing their compatibility with the MW observables. Conclusions: The same r-process enrichment frameworks cannot reproduce simultaneously the Eu patterns both in the MW and in dwarf galaxies. However, a scenario where additional Eu is provided by an increased production from delayed sources at low metallicity can theoretically reconcile the trends observed in the MW and in LG dwarfs, because of the small discrepancies (0.1 dex) between model predictions and observations found in this case. Further targeting and modelling of neutron-capture elements in LG galaxies are however needed to fill the gaps in our current understanding of the problem.

Observations of the Galactic Center using Imaging Atmospheric Cherenkov Telescopes (IACTs), such as H.E.S.S., MAGIC, and VERITAS, have revealed a very-high-energy (VHE, $\gtrsim 100$ GeV) gamma-ray source, HESS J1745$-$290, aligned with the dynamical center of the Milky Way. This source shows point-like emission ($\lesssim 0.1^\circ$) and a strong suppression in its energy-differential spectrum in the ten TeV energy regime, modeled well by a power-law with an exponential cutoff. The origin of this emission is debated, with candidate emission scenarios including dark matter annihilations, millisecond pulsars in the central stellar clusters, and hadronic interactions in the vicinity of Sagittarius A*. Deriving the sensitivity to these spectral models is key to discriminating the physical processes at work. We show that combining H.E.S.S., MAGIC, and VERITAS archival data can well described the observed emission by a power-law with an exponential energy cutoff within the present uncertainties. Given the near advent of the array of the Large-Sized Telescopes (LSTs) at CTAO-N, we timely simulate realistic upcoming observations of the central emission by the CTAO-N four-LST array, to derive the sensitivity to resolve the sharpness of the spectral energy cutoff. We find that 500 hours of four-LST observations taken at large zenith angles, possibly accumulated over several years, can significantly discriminate the dark-matter emission scenario from the leptonic and hadronic ones. Also, a preliminary 3$\sigma$ hint for such discrimination could emerge within the first year. We demonstrate, for the first time, that CTAO-N is able to provide new insights on differentiating among the above-mentioned emission scenarios in the next several years.

Ultra-light dark matter (ULDM) is a compelling candidate for cosmological dark matter. If ULDM interacts with ordinary matter, it can induce measurable, characteristic signals in pulsar-timing data because it causes the orbits of pulsars in binary systems to osculate. In this work, we investigate the potential of machine learning (ML) techniques to detect such ULDM signals. To this end, we construct three types of neural networks: an autoencoder, a binary classifier, and a multiclass classifier. We apply these methods to four theoretically well-motivated ULDM models: a linearly coupled scalar field, a quadratically coupled scalar field, a vector field and a tensor field. We show that the sensitivity achieved using ML methods is comparable to that of a semi-analytical Bayesian approach, which to date has only been applied to the linear scalar case. The ML approach is readily applicable to all four ULDM models and, in the case of the multiclass classifier, can distinguish between them. Our results, derived from simulated data, lay the foundation for future applications to real pulsar-timing observations.

We present the HI distribution, kinematics, mass modeling, and disk stability of the dwarf irregular galaxies ESO444-G084 and [KKS2000]23 using high-resolution, high-sensitivity MHONGOOSE survey data from MeerKAT. ESO444-G084 shows centrally concentrated HI emission, while [KKS2000]23 exhibits irregular high-density clumps. Total HI fluxes measured down to 10^19 and 10^18 cm^-2 are nearly identical, indicating that the increased HI diameter at lower column densities results mainly from the larger beam, with no significant extra emission detected. We derive total HI masses of (1.1 +/- 0.1) x 10^8 and (6.1 +/- 0.3) x 10^8 solar masses for ESO444-G084 and [KKS2000]23, respectively. Using PyFAT and TiRiFiC, we extract 3D rotation curves that reveal disk-like kinematics in both galaxies. ESO444-G084 shows a warp beyond ~1.8 kpc and a fast-rising curve consistent with a centrally concentrated dark matter distribution, while [KKS2000]23's more gradual rise suggests a more extended halo. Mass modeling with an isothermal halo and stellar mass-to-light ratios of 0.20 for ESO444-G084 and 0.18 for [KKS2000]23 yields consistent results. We analyze disk stability using spatially resolved Toomre Q and gas-to-critical surface density ratios, linking these with H-alpha and FUV-based star formation. ESO444-G084 supports localized star formation despite global stability, while [KKS2000]23 appears gravitationally unstable yet lacks H-alpha, suggesting that turbulence, gas depletion, or past feedback suppresses star formation. No inflows or outflows are detected, indicating internal processes regulate star formation. This study highlights the interplay between HI morphology, kinematics, dark matter distribution, and disk stability, showing how internal processes shape dwarf galaxy evolution.

The Burst Alert Telescope (BAT) onboard the Neil Gehrels Swift observatory has been serving as a survey instrument for the hard X-ray sky, and has detected thousands of X-ray sources (e.g., AGNs, X-ray binaries, etc). BAT monitors these X-ray sources and follows their light curves on time scales from minutes to years. In addition, BAT discovers hundreds of new X-ray sources in survey images stacked throughout the mission lifetime. We present the updated BAT survey catalog since the last published BAT 105 month survey catalog (Oh et al. 2018) with additional of 4.5 years of data until December 2017. Data since 2007 are reprocessed to include updated instrumental calibration. Analysis in this study shows that additional systematic noise can be seen in the 157-month mosaic images, resulting in decreases in the expected improvement in sensitivity and the number of new detections. The BAT 157-month survey reaches a sensitivity of $8.83 \times 10^{-12} \rm \ erg \ s^{-1} \ cm^{-2}$ for 90\% of the sky and $6.44 \times 10^{-12} \rm \ erg \ s^{-1} \ cm^{-2}$ for 10\% of the sky. This catalog includes spectra, monthly and snapshot light curves in eight energy bands (14-20, 20-24, 24-35, 35-50, 50-75, 75-100, 100-150, and 150-195 keV) for 1888 sources, including 256 new detections above the detection threshold of $4.8\sigma$. The light curves, spectra, and tables that summarize the information of the detected-sources are available in the online journal and in the catalog web page this https URL.

In this work, we propose a novel method to classify close binary stars, derived from the dynamical structure inherent in their light curves. We apply the technique to light curves of binaries from the revised Kepler Eclipsing binary catalog, selecting close binaries which have the standard morphology parameter, $c$, $\gt 0.5$ corresponding to semi-detached, over-contact and ellipsoidal systems. Using the method of time delay embedding, we recreate the non-linear dynamics underlying the data and quantify the patterns of recurrences in them. Using two recurrence measures, Determinism and Entropy, we define a new Dynamically Derived Morphology (DDM) parameter and compute its values for the Kepler objects. While as expected, this metric is somewhat inversely correlated with the existing morphology parameter (Spearman $\rho= -0.32$), the method offers an alternate classification scheme for close binary stars that captures their nonlinear dynamics, an aspect often overlooked in conventional methods. Hence, the DDM parameter is expected to distinguish between stars with similar folded light curves, but are dynamically dissimilar due to nonlinear effects. Moreover, since the method can be easily automated and is computationally efficient it can be effectively used for future sensitive large data sets.

The unexpected finding of a ring system around the Centaur (10199) Chariklo opened a new window for dynamical studies and posed many questions about the formation and evolutionary mechanisms of Centaurs as well as the relationship to satellites and outbursting activity. As minor planets that cross the orbits of the giant planets, Centaurs have short dynamical lifetimes: Centaurs are supplied from the trans-Neptunian region and some fraction migrates inward to become Jupiter-family comets. Given these dynamical pathways, a comparison of attributes across these classifications provides information to understand the source population(s) and the processes that have affected these minor planets throughout their lifetimes. In this chapter we review the current knowledge of satellites, rings, and debris around Centaur-like bodies, discuss the observational techniques involved, place the information into context with the trans-Neptunian Objects, and consider what the results tell us about the outer solar system. We also examine open questions and future prospects.

We report the delivery to the Mikulski Archive for Space Telescopes (MAST) of tables containing Root Mean Square (RMS) Combined Differential Photometric Precision (CDPP) values for all TESS 2-min cadence targets with Science Processing Operations Center (SPOC) light curves in Sectors 1-90. Each comma-separated values (CSV) file contains CDPP values for all 2-min light curves in the given sector. The tables include robust RMS CDPP values for the 15 trial transit pulse durations searched in the SPOC 2-min processing pipeline, ranging from 0.5-15.0 hr. For each pulse duration, CDPP is computed in the transit search for a trial transit centered on every cadence. The RMS value of the CDPP time series is a metric that may be employed to estimate signal-to-noise ratio for transits with the given duration and a specified depth. We will continue to deliver the RMS CDPP tables to MAST for each observing sector.

The COSMOS-Web survey, with its unparalleled combination of multiband data -- particularly near-infrared imaging from JWST's NIRCam (F115W, F150W, F277W, and F444W), provides a transformative dataset reaching down to approximately 28th magnitude (F444W) for studying galaxy evolution. In this work, we employ Self-Organizing Maps (SOMs), an unsupervised machine learning method, to estimate key physical parameters of galaxies, redshift, stellar mass, star formation rate (SFR), specific SFR (sSFR), and age, directly from photometric data out to redshift $z = 3.5$. SOMs efficiently project high-dimensional galaxy color information onto two-dimensional maps, revealing how physical properties vary among galaxies with similar spectral energy this http URL first validate our approach using mock galaxy catalogs from the HORIZON-AGN simulation, where the SOM accurately recovers the true parameters, demonstrating its robustness. Applying the method to COSMOS-Web observations, we find that the SOM delivers reliable estimates despite the increased complexity of real galaxy populations. Performance metrics (NMAD and RMSE typically between 0.1 and 0.3) confirm the method's precision, with about 70\% of predictions falling within one standard deviation (dex) of the reference values. While redshift estimation in COSMOS-Web presents challenges (median NMAD = 0.05), the overall success of the SOM highlights its potential as a powerful and interpretable tool for galaxy parameter estimation.

We propose a new method that reveal the nature of dark energy (DE) evolution. Specifically, the method consists of studying the evolving trend regarding the effective running Hubble constant: when it increases, it indicates a quintessence nature, and when it decreases, it reveals a phantom behavior. Within the framework of the dark energy models we analyze three parameterizations: the $w$CDM model, a reduced Chevallier-Polarski-Linder (CPL) model and a new theoretical model based on the possible creation of dark energy by the time-varying gravitational field of the expanding Universe. For each DE model, we construct a theoretical effective running Hubble constant, i.e. a function of the redshift, which highlights the difference between modified dynamics and the $\Lambda$CDM-one. Furthermore, these dark energy models are compared to the phenomenological model of a decreasing trend of the Hubble constant as a function of the redshift, called the power-law model (PL) and the $\Lambda$CDM one. These three theoretical functions are fitted against the binned SNe Ia data samples, i.e. the Pantheon and the Master samples, the latter being a collection of SNe Ia from 4 catalogs: Dark Energy Survey (DES), PantheonPlus, Pantheon and Joint Lightcurve Analysis (JLA), without duplicated SNe Ia, called the Master sample. The main result of our study is that the phenomenological PL model is statistically favored compared to the other proposed scenarios, both for the Pantheon and the Master samples. At this stage, the SNe Ia data do not indicate that the evolution of dark energy models among the studied ones is favored respect to the $\Lambda$CDM. Nevertheless, the binned Pantheon sample allows for a discrimination of the nature of dark energy at least at the $1\,\sigma$ level via the fit of the effective running Hubble constant.

In the past two decades, transit surveys have revealed a class of planets with thick atmospheres -- sub-Neptunes -- that must have completed their accretion in protoplanet disks. When planets form in the gaseous disk, the gravitational interaction with the disk gas drives their migration and results in the trapping of neighboring planets in mean motion resonances, though these resonances can later be broken when the damping effects of disk gas or planetesimals wane. It is widely accepted that the outer Solar System gas giant planets originally formed in a resonant chain, which was later disrupted by dynamical instabilities. Here, we explore whether the early formation of the terrestrial planets in a resonance chain (including Theia) can evolve to the present configuration. Using N-body simulations, we demonstrate that the giant planet instability would also have destabilized the terrestrial resonance chain, triggering moon-forming giant impacts in 20--50\% of our simulated systems, dependent on the initial resonance architecture. After the instability, the eccentricity and inclination of the simulated planets match their present-day values. Under the proposed scenario, the current period ratio of 3.05 between Mars and Venus -- devoid of any special significance in traditional late formation models -- naturally arises as a relic of the former resonance chain.

Fast magnetosonic waves are one of the two low-frequency plasma modes that can exist in a neutron star magnetosphere. It was recently realized that these waves may become nonlinear within the magnetosphere and steepen into some of the strongest shocks in the universe. These shocks, when in the appropriate parameter regime, may emit GHz radiation in the form of precursor waves. We present the first global Particle-in-Cell simulations of the nonlinear steepening of fast magnetosonic waves in a dipolar magnetosphere, and quantitatively demonstrate the strong plasma acceleration in the upstream of these shocks. In these simulations, we observe the production of precursor waves in a finite angular range. Using analytic scaling relations, we predict the expected frequency, power, and duration of this emission. Within a reasonable range of progenitor wave parameters, these precursor waves can reproduce many aspects of FRB observations.

The Vera Rubin Observatory's Legacy Survey of Space and Time (LSST) is expected to revolutionize time-domain optical astronomy as we know it. With its unprecedented depth, the LSST will survey the southern hemisphere sky, generating nearly 32 trillion observations over its nominal 10-year operation. Among these, approximately 10 million will be supernovae (SNe). These observations will uniquely characterize the SN population, enabling studies of known and rare SN types, detailed parameterization of their light curves, deep searches for new SN progenitor populations, the discovery of strongly lensed SNe, and the compilation of a large, well-characterized sample of superluminous SNe. We analyzed a sample of 22663 simulations of LSST light curves for core collapse SNe (CCSNe), modeled using the radiative transfer code STELLA. We analyzed this dataset with the software CASTOR, which enables the reconstruction of synthetic light curves and spectra via a machine learning technique that allows one to retrieve the complete parameter map of a SN. For each parameter we compared the observed and the true values, determining how LSST light curves alone will contribute to characterize the progenitor and the explosion. Our results indicate that LSST alone will not suffice for a comprehensive and precise characterization of progenitor properties and explosion parameters. The limited spectral coverage of LSST light curves (in most cases) does not allow for the accurate estimation of bolometric luminosity, and consequently, of the explosion energy and nickel yield. Additionally, the redshift-absorption degeneracy is difficult to resolve without supplementary information. These findings suggest that for the most interesting SNe, complementary follow-up observations using spectrographs and optical facilities (particularly in the infrared bands) will be essential for accurate parameter determination.

The distribution of baryons in the Universe remains a fundamental open question in astronomy, and the dispersion measure (DM) of Fast Radio Bursts (FRBs) serves as a valuable tool for probing this cosmic gas. We investigate the impact of the foreground cosmic web on FRB DMs using 61 localized FRBs and public galaxy catalogs. We test for the large-scale structure's impact on cosmological DM using two methods. First, we searched for a correlation between galaxy number density along the line of sight and extragalactic DM, and found a statistically significant positive correlation ($p$ = $1.76 \times 10^{-5}$). The shape of this correlation contains information about the cosmic baryon distribution, and can also be used to better constrain host galaxy DM by providing an estimate of the cosmic contribution on a per-source basis. We observe similar correlations in a mock FRB survey based on the IllustrisTNG cosmological simulation, where the DM is dominated by filaments in the IGM and not by halos. Next, we performed a stacking analysis that measures the average excess DM as a function of impact parameter of foreground galaxies to obtain spatial information about how ionized gas is distributed around galaxy halos. We report excess DM in the stacked signal for impact parameters up to Mpc scales ($\sim$3$\sigma$). Finally, we identified FRBs that do not appear to intersect intervening halos within $r_{vir}$, allowing us to estimate the fraction of baryons that reside in the IGM. We find $f_{\mathrm{IGM}} \geq 0.69$ at 95$\%$ confidence, indicating significant astrophysical feedback.

Gravitational lensing of gravitational waves (GWs) offers a novel observational channel that complements traditional electromagnetic approaches and provides unique insights into the astrophysical environments of GW sources. In this work, we investigate the repeated lensing of continuous gravitational wave (CW) sources in active galactic nucleus (AGN) disks by central supermassive black holes (SMBHs), focusing on the imprint of SMBH spin via the Lense-Thirring (LT) effect. Although typically weak and challenging to observe, the spin-induced precession of source orbits can accumulate over time, thereby modulating the lensing geometry. Such modulations influence the magnification, duration, and waveform structure of each repeated lensing event, and enhance the overall probability of lensing occurrences. Using matched filtering, we demonstrate that spin-dependent signatures may be detectable, suggesting that lensed CW signals could serve as an indirect probe of SMBH spin in AGNs.

This work presents a global solution for the internal and the external field of an axisymmetric rotating neutron star. It is shown that the twist of the internal field affects the external field, by increasing the number of open field lines and eventually the spin-down rate of the star. This effect is far more drastic if the toroidal field, and consequently the poloidal current flowing within the star, is allowed to populate the closed field lines of the magnetosphere, rather than if it remains confined in the star. We further remark that the internal field structure depends on the presence of a twisted magnetosphere: if the twist current is not allowed to flow in the magnetosphere it only occupied a narrow toroid at the interior of the star, whereas if the twist currents are allowed to flow in the magnetosphere the internal toroidal field may occupy a significant volume of the stellar interior. Strong magnetospheric currents may also impact the emission mechanisms, and lead to fluctuations in magnetar spin-down rates, moding and nulling of pulsars, a correlation between angular shear and twist, and the general morphology of the pulsar magnetic field leading to various observational manifestations. The magnetospheric toroidal fields may possibly dissipate, thus the system may switch from global twist to internal twist and consequently exhibit transient behavior.

Systematic abundance differences that depend on the condensation temperatures of elements have been observed, in particular for stars similar to the Sun; solar twins and solar analogs. Similar differences have also recently been shown to exist between solar abundances and abundances of refractory elements in primitive chondrites. Numerous mechanisms have been proposed to account for these effects, including also differences observed in binary systems. Rather than relying on specific mechanisms, this paper aims to show that the observed effects are a natural and unavoidable outcome of the star formation process itself, in which the associated outflows (winds and jets) carry away material, with efficiency varying with condensation temperature. By using analysis based on modeling results and scaling laws, the trends and magnitudes of the effects are investigated, in three contexts: 1) with respect to differences between the Sun and solar twins, 2) with respect to differences between the Sun and CI-chondrites, and 3) with respect to differences between members of binaries. It is shown that protoplanetary disk outflows indeed are expected to have differential abundance effects, with trends and magnitudes consistent with observed abundance effects. The qualitative as well as semi-quantitative character of the effects are reproduced, in all three contexts. The results demonstrate that the observed systematic differences are likely results of the disk outflows associated with the accretion process. In contrast to mechanisms relying on the tiny mass of planets leaving an observable signature, outflows carry away masses similar to the entire mass of the star, thus much more easily resulting in differential effects with the magnitudes observed, without for example having to assume that the abundance differences are limited to the convection zones of the stars.

The upcoming Vera Rubin Observatory Legacy Survey of Space and Time (LSST) will dramatically increase the number of strong gravitational lensing systems, requiring precise modeling of line-of-sight (LOS) effects to mitigate biases in lensing observations and cosmological inferences. We develop a method to construct joint distributions of external convergence ($\kappa_{\mathrm{ext}}$) and shear ($\gamma_{\mathrm{ext}}$) for strong lensing LOS by aggregating large-scale structure simulations with high-resolution halo renderings and non-linear correction. Our approach captures both smooth background matter and perturbations from halos, enabling accurate modeling of LOS effects. We apply non-linear LOS corrections to $\kappa_{\mathrm{ext}}$ and $\gamma_{\mathrm{ext}}$ that address the non-additive lensing effects caused by objects along the LOS in strong lensing. We find that, with a minimum image separation of $1.0^{\prime\prime}$, non-linear LOS correction due to the presence of a dominant deflector slightly increases the ratio of quadruple to double lenses; this non-linear LOS correction also introduces systematic biases of $\sim 0.1\%$ for galaxy-AGN lenses in the inferred Hubble constant ($H_0$) if not accounted for. We also observe a $0.66\%$ bias for galaxy-galaxy lenses on $H_0$, and even larger biases up to $1.02\%$ for galaxy-AGN systems if LOS effects are not accounted for. These results highlight the importance of LOS for precision cosmology. The publicly available code and datasets provide tools for incorporating LOS effects in future analyses.

Efficient turbulent acceleration of particles is indicated by recent astrophysical observations, but its mechanism is not well understood. Mirror acceleration has recently been proposed as an efficient mechanism for particle energization in turbulence-compressed magnetic fields. We employ a 3D particle-in-cell (PIC) simulation of pair plasma in magnetized and relativistic turbulence to study this new mechanism and its acceleration efficiency. By tracking individual particles, we see that reversal of a particle's moving direction and significant energy gain can happen during one mirror interaction and within one gyro-orbit. As expected for mirror acceleration, we statistically find that (1) energy gain is preferentially in the direction perpendicular to the local magnetic field and positively correlated with local magnetic field strengthening, and (2) the particle pitch angle distribution becomes increasingly anisotropic toward higher energies, with a concentration at large pitch angles. Our results demonstrate that the mirror acceleration causes a strong confinement of particles by stochastically increasing their pitch angles. This, in turn, facilitates repeated mirror acceleration with the mirroring condition well satisfied. We conclude that mirror acceleration is a promising mechanism accounting for efficient acceleration in magnetized and turbulent astrophysical plasmas.

We study the stellar distribution around supermassive black holes (SMBHs) in gas-rich nuclear star clusters (NSCs). NSCs could contain vast amounts of gas, which contribute significantly to shaping the stellar distribution, typically altering the stellar density cusp from the usual Bahcall \& Wolf 1976 solution and consequently affecting the dynamics in the NSC. The dense gaseous environment in NSCs gives rise to dynamical phenomena that are otherwise rare in other gas-free environments. Here we extend the derivation introduced in Bahcall \& Wolf 1976 to include an additional energy dissipation term associated with gas drag. We examine the effect of different forms of gas drag on the stellar density distribution. Finally, we discuss implications on the rates of tidal disruption events and other transients triggered by stellar interactions in gas-rich galactic nuclei.

The latest DESI DR2 results, when combined with other independent cosmological data on the Cosmic Microwave Background and supernovas, suggest a preference for dynamical dark energy. We propose a novel cosmological scenario, which features two distinct scalar fields. One governs the magnitude of the present-day dark energy density and is related to the size of extra-dimensions. Accounting for the observed smallness of this energy density requires the scalar to reside near the boundary of field space. The second field, responsible for the time evolution of dark energy and associated with the string coupling, must instead lie in the bulk to remain consistent with the non-observation of light string states. We show that a natural candidate for such dark energy dynamics is a quintessence modular-invariant potential, in which the second scalar field rolls down a negatively curved slope, starting from a self-dual critical point. We find that this scenario is in good agreement with the latest findings by DESI.

Miscidynamics describes coarse-grained neutrino transport under the assumption that flavor mixing is in local equilibrium. Here we introduce the concept of turbulent flavor-wave viscosity and develop techniques for including it in miscidynamics. This extension of the theory is necessary when neutrinos develop weak flavor instabilities as a result of astrophysical driving. The flavor-wave dispersion relation is obtained using a new and more general form of linear analysis, which does not require small flavor coherence. Flavor-wave transport is then approximated using wave kinetics and geometric optics. We hypothesize that the dynamical emergence of local oscillation invariants restricts the extent of flavor thermalization. In this view, flavor evolution is sculpted by both entropy-increasing and order-promoting factors.

We explore the end state of gravitational collapse under quantum gravity effects and propose that Planck Star Remnants (PSR), formed via nonsingular bounces, could serve as viable dark matter candidates. Within the framework of Loop Quantum Cosmology, we model the collapse of a homogeneous matter distribution and show that the classical singularity is replaced by a quantum bounce at the Planck density. By analytically matching the Friedmann Lemaitre Robertson Walker (FLRW) interior to an exterior Schwarzschild spacetime using the Israel junction conditions, we demonstrate that the bounce remains causally hidden from external observers, avoiding any observable re-expansion. This naturally leads to the formation of stable, non-radiating PSR, whose radius coincides with the Schwarzschild radius when the black hole mass approaches the Planck mass as a result of Hawking evaporation. We suggest that such remnants may originate from evaporating primordial black holes in the early universe, and estimate the relic abundance needed for PSR to account for the observed dark matter density. We also discuss some crucial differences between PSR and previous proposals of Planck mass relics. The scenario is shown to be consistent with existing astrophysical and cosmological constraints, offering a unified framework connecting quantum gravitational collapse, and the nature of dark matter.

The QCD axion has important connections to early universe cosmology. For example, it is often said that isocurvature limits rule out a combination of high axion decay constant, $f_a$, and high inflationary Hubble scale, $H_I$. High scales are theoretically motivated, so it is important to ask how robust this constraint is. We demonstrate that this constraint is naturally evaded when the quartic coupling of the complex $U(1)_\mathrm{PQ}$-breaking field is small. In this case, $f_a$ changes from a larger value during inflation to a smaller value in the later universe, suppressing isocurvature perturbations. Importantly, we show that in large parts of parameter space this solution is not jeopardised by overproduction of the axion through parametric resonance. The isocurvature bounds are thus dependent on UV physics. We have found that, even for the minimal QCD axion, large parts of UV parameter space at both high $f_a$ and high $H_I$ are in fact allowed, not ruled out by isocurvature constraints.

Analog gravity experiments are making remarkable strides in unveiling both the classical and quantum nature of black holes. By harnessing diverse states of matter, contemporary tabletop setups now replicate strong-field phenomena typically confined to the enigmatic regions surrounding black holes. Through these modern gravity simulators, physical processes once considered elusive may finally be brought into experimental reach. In this work, we investigate the bound spectrum of massless scalar excitations propagating within the effective geometry of a rotating acoustic metric. Specifically, we utilize an analog vortex endowed with a tunable parameter that emulates the spacetime of a rotating gravitational background. This model accommodates both the presence of a sonic horizon--characteristic of an acoustic black hole--for non-zero tuning parameters, and its absence when the parameter vanishes, yielding a horizonless, purely rotational vortex flow devoid of radial inflow. We focus on the latter case, where the vortex flow is purely rotational, and compute the spectral properties of the analog system. The resulting bound-state spectrum is found to be qualitatively consistent with that observed in recent experimental realizations of superfluid Helium giant quantum vortices featuring solid or hollow cores. This correspondence suggests that the analog spacetime geometry used here holds significant potential to replicate the phenomenology of cutting-edge laboratory experiments. In doing so, it offers new insight into the intricate landscape of analog black hole spectroscopy and, potentially, the physical topography of bounded, rotating astrophysical environments around black holes.

Initial orbital eccentricities of gravitational wave (GW) events associated with merging binary black holes (BBHs) should provide clues to their formation scenarios, mainly because various BBH formation channels predict distinct eccentricity distributions. However, searching for inspiral GWs from eccentric BBHs is computationally challenging due to sophisticated approaches to model such GW events. This ensures that Bayesian parameter estimation methods to characterize such events are computationally daunting. These considerations influenced us to propose a novel approach to identify and characterize eccentric BBH events in the LIGO-Virgo-KAGRA (LVK) collaboration data sets that leverages external attention transformer models. Employing simulated data that mimic LIGO O4 run, eccentric inspiral events modeled by an effective-one-body numerical-relativity waveform family, we show the effectiveness of our approach. By integrating this transformer-based framework with a convolutional neural network (CNN) architecture, we provide efficient way to identify eccentric BBH GW events and accurately characterize their source properties.

A quasi-optical (QO) test bench was designed, simulated, and calibrated for characterizing S-parameters of devices in the 220-330 GHz (WR-3.4) frequency range, from room temperature down to 4.8 K. The devices were measured through vacuum windows via focused beam radiation. A de-embedding method employing line-reflect-match (LRM) calibration was established to account for the effects of optical components and vacuum windows. The setup provides all four S-parameters with the reference plane located inside the cryostat, and achieves a return loss of 30 dB with an empty holder. System validation was performed with measurements of cryogenically cooled devices, such as bare silicon wafers and stainless-steel frequency-selective surface (FSS) bandpass filters, and superconducting bandpass FSS fabricated in niobium. A permittivity reduction of Si based on 4-GHz resonance shift was observed concomitant with a drop in temperature from 296 K to 4.8 K. The stainless steel FSS measurements revealed a relatively temperature invariant center frequency and return loss level of 263 GHz and 35 dB on average, respectively. Finally, a center frequency of 257 GHz was measured with the superconducting filters, with return loss improved by 7 dB on average at 4.8 K. To the best of our knowledge, this is the first reported attempt to scale LRM calibration to 330 GHz and use it to de-embed the impact of optics and cryostat from cryogenically cooled device S-parameters.

We discuss recent progress in the study of entanglement within cosmological frameworks, focusing on both momentum and position-space approaches and also reviewing the possibility to directly extract entanglement from quantum fields. Entanglement generation in expanding spacetimes can be traced back to the phenomenon of gravitational particle production, according to which the background gravitational field may transfer energy and momentum to quantum fields. The corresponding entanglement amount and its mode dependence are both sensitive to the field statistics and to the details of spacetime expansion, thus encoding information about the background. Gravitational production processes also play a key role in addressing the quantum-to-classical transition of cosmological perturbations. In order to directly extract entanglement from quantum fields, local interactions with additional quantum systems, working as detectors, have been suggested, leading to the formulation of the entanglement harvesting protocol. Despite harvesting procedures are currently unfeasible from an experimental point of view, various proposals for implementation exist and a proper modeling of detectors and local interactions is crucial to address entanglement extraction via realistic setups. In the final part of the work, we address entanglement characterization in position space, primarily focusing on black hole spacetimes. We first investigate a possible interpretation of Bekenstein-Hawking black hole entropy in terms of the entanglement entropy arising in discrete quantum field theories, on account of the area law. Then, we discuss the resolution of the black hole information paradox via the gravitational fine-grained entropy formula, which provides a new way to compute the entropy of Hawking radiation and allows to preserve unitarity in black hole evaporation processes.

We point out that inflationary superhorizon fluctuations can be effectively described by a set of equations analogous to those governing a superfluid. This is achieved through a functional Schrödinger approach to the evolution of the inflationary wavefunction, combined with a suitable coarse-graining procedure to capture large-scale dynamics. The irrotational fluid velocity is proportional to the gradient of the wavefunction phase. Marginalizing over short superhorizon modes introduces an external force acting on the fluid velocity. The quantum pressure characteristic of the superfluid plays a role in scenarios involving an ultra-slow-roll phase of inflation. Our superfluid framework is consistent with the standard Starobinsky approach to stochastic inflation while offering complementary insights, particularly by providing more precise information on the phase of the inflationary wavefunction. We also discuss a heuristic approach to include dissipative effects in this description.