Papers for Thursday, Jul 10 2025

Protoclusters represent the most extreme environments in the very early Universe. They form from large-scale dark matter overdensities, harbouring an overabundance of galaxies fed by large gas reservoirs. Their early and accelerated evolution results in a distinct difference in the properties of galaxies resident in protoclusters versus the field, which is known to be in place by $z\sim 5-6$. We utilise JWST NIRCam observations of the A2744-z7p9OD protocluster at $z=7.88$ to constrain the properties of resident galaxies. We identify seven new protocluster members, bringing the total number to 23 and the total stellar mass of the protocluster to in excess of $10^{10}\ \rm{M_{\odot}}$. These galaxies are remarkably evolved just 650 Myr after the Big Bang, preferentially showing redder UV-slopes and stronger Balmer breaks than is typical of field galaxies. We use the PROSPECTOR spectral energy distribution fitting code to derive key galaxy properties, finding distinct populations in the core versus the outskirts of the protocluster. The core is largely composed of dusty, massive galaxies which can be characterised as undergoing a synchronised (mini)-quenched phase, while galaxies in the protocluster outskirts are undergoing recent bursts of star formation. Finally, a strong suppression of the continuum around the Ly$\alpha$-break evidences extreme neutral hydrogen column densities in many resident galaxies ($N_{\rm HI}\gtrsim10^{23}\ {\rm cm^{-2}}$). The A2744-z7p9OD system is the most extreme, evolved overdensity yet observed at $z>7$, with higher stellar masses, gas densities, and dust attenuation, revealing the intersection of local environment and high-redshift galaxy formation at their extremes.

Optical observations of gamma-ray bursts (GRBs) contemporaneous with their prompt high-energy emission are rare, but they can provide insights into the physical processes underlying these explosive events. The Transiting Exoplanet Survey Satellite's (TESS) large field of view and continuous observation capabilities make it uniquely positioned to detect and characterize prompt optical flashes from GRBs. In this work, we fit phenomenological models to the gamma-ray through optical spectral energy distributions (SEDs) of 24 bursts with arcsecond-level localizations that fell within the TESS field of view between 2018 July and 2024 December. In four cases, the extrapolation of the high-energy SED agrees with the observed optical flux to within 1-$\sigma$. In one case, there is a significant excess of optical flux relative to the extrapolation. In two cases, upper limits from TESS did not constrain the optical portion of the SED. In the remaining 17 cases, the optical flux is overpredicted by the extrapolation from high energies. This discrepancy could be explained by dust extinction in the host galaxy.

The internal structures of Uranus and Neptune remain unknown. In addition, sub-Neptunes are now thought to be the most common type of exoplanets. Understanding the physical processes that govern the interiors of such planets is therefore essential. Phase separation between hydrogen and water may occur in cold, water-rich intermediate-mass planets. We assess whether it could occur in Uranus, Neptune, K2-18\,b and TOI-270\,d, and investigate its effect on the planetary evolution and inferred internal structure. We couple planetary evolution models with recent \textit{ab initio} calculations of the hydrogen-water phase diagram, allowing for temperature shifts to account for uncertainties in miscibility gaps. We find that demixing may occur and could lead to a complete depletion of water in the outermost regions of Uranus and Neptune. Temperature offsets of up to 1100~K lead to a depleted region comprising as much as 16\% of the planet's mass, and an increase in planetary radius by nearly 20\%. For K2-18\,b, our models suggest that hydrogen-water demixing is ongoing and may explain the absence of water features in its JWST spectrum. A temperature offset of 500~K is required to get a complete depletion of water in the atmosphere of K2-18\,b. TOI-270\,d may also have experienced hydrogen-water demixing. When applying a similar temperature offset on the phase diagram as for K2-18\,b, we find a partial depletion of water in the atmosphere of TOI-270\,d, consistent with JWST's detection of water. Hydrogen-water immiscibility may play a key role in shaping the structure and evolution of both Solar System giant planets like Uranus and Neptune, and cold/temperate exoplanets such as K2-18\,b and TOI-270\,d. Accounting for such internal processes is crucial to accurately interpret atmospheric observations from current (e.g., JWST) and upcoming (e.g., ARIEL) missions.

We present a 152 ks XRISM/Resolve observation of the persistently accreting Z source GX 340+0. Simultaneous observations also occurred with NuSTAR and NICER for 22.47 ks and 2.7 ks, respectively. The source covered the normal branch to the flaring branching during the observations. The data from all three missions were modeled concurrently for each spectral branch. The superior energy resolution of XRISM/Resolve reveals structure within the iron emission line complex regardless of spectral state. We model the reprocessed Fe K line with a reflection model tailored for thermal illumination of the accretion disk by a neutron star. The currently available model encompasses the broad components, but narrow emission features remain at the ~5% level. These remaining features may be described by the presence of an ionized plasma in the system as has been observed in the Z source Cygnus X-2, but subsequent updates to the reflection model code may be able to explain these features.

Using the MOPTOP polarimeter on the 2-meter Liverpool Telescope, we conducted a monitoring campaign targeting three optically-discovered TDEs (AT2024bgz, AT2024pvu, and AT2024wsd) and two Bowen flares in AGN (AT2020afhd and AT2019aalc). The three thermal TDEs showed low intrinsic polarization levels ($\Delta P \approx 0$-$6\%$) with stable polarization angles. The Bowen flares also showed variable polarization degree ($ \Delta P \approx 0$-$8\%$), but significant polarization angle variability: AT2020afhd exhibited a $\Delta \theta = 83 \pm 8 \,^\circ$ shift at 150 days post-optical peak, while the AT2019aalc displayed quasi-periodic swings of $\Delta \theta \approx 40 \, ^\circ$ amplitude starting 190 days after peak brightness. The TDEs of this study are well described by models invoking rapid disk formation and reprocessed emission from optically thick outflows, whereas the Bowen flares reveal more complex reprocessing geometries, potentially consistent with TDEs occurring in AGN gas-rich environments. We find that moderate polarization is observed at later times for TDEs with low-Eddington ratios and highly extended photospheres. This implies that, as the accretion level declines, we expect more asymmetric reprocessing layers along a given viewing angle. Since the outflow density and velocity depend sensitively on the inclination angle, we expect TDEs with low-Eddington ratios and highly extended photospheres to exhibit varying levels of polarization. The polarization of AT2019aalc (Seyfert 1) hints at a clumpy, asymmetric outflow and the presence of a tilted, precessing accretion disk, while the polarization of AT2020afhd (AGN type 2) is consistent with the detection of a scattered light echo.

The magnetic field is known to play a crucial role in star formation. Dust polarization is an effective tool for probing the morphology of the field, yet it does not directly trace its strength. Several methods have been developed, combining polarization and spectroscopic data, to estimate the strength of the magnetic field, including the DCF method, which relates these quantities to the magnetic-field strength under the assumption of Alfvénic turbulence. Skalidis & Tassis (2021) (ST), relaxed this assumption to account for the compressible modes, deriving more accurate estimates of the field strength. We evaluate the accuracy of these methods in star-forming regions and propose a systematic approach for calculating the key observational parameters involved: the velocity dispersion (dv), the dispersion of polarization angles (d\theta), and the cloud density (rho). We use a 3D MHD chemodynamical simulation of a turbulent molecular cloud and generate synthetic observations, for seven different inclination angles. We employ various approaches for estimating the parameters dv, d\theta, and rho and find that the approach used to calculate these parameters plays a crucial role in estimating the magnetic field strength. We show that the value probed by both the DCF and ST methods corresponds to the median of the molecular-species-weighted POS component of the magnetic field. The ST outperforms DCF, accurately following the expected cosine trend with respect to the inclination angle, and remains within 1\sigma from the true strength of the field. Based on our analysis, we proposed that, in self-gravitating clouds, the intrinsic parameters (rho, dv, d\theta) should be calculated as follows: rho using radiative transfer analysis, dv using the second moment maps, and d\theta by fitting Gaussians to the polarization angle distributions to remove the contribution of the hourglass morphology.

Proto-planetary disks display the so-called size-luminosity relation, where their mm-wavelength fluxes scale linearly with their emitting areas. This suggests that these disks are optically thick in mm-band, an interpretation further supported by their near-black-body spectral indexes. Such characteristics are seen not only among disks in very young star-forming regions like Lupus (1-3 Myrs), but, as we demonstrate here, also among disks in the much older Upper Scorpius region (5-11 Myrs). How can disks shine brightly for so long, when grain growth and subsequent radial drift should have quickly depleted their solid reservoir? Here, we suggest that the "bouncing barrier" provides the answer. Even colliding at very low speeds (below 1cm/s), grains already fail to stick to each other but instead bounce off in-elastically. This barrier stalls grain growth at a near-universal size of $\sim 100 \mu m$. These small grains experience much reduced radial drift, and so are able to keep the disks bright for millions of years. They are also tightly coupled to gas, offering poor prospects for processes like streaming instability or pebble accretion. We speculate briefly on how planetesimals can arise in such a bath of 100-micron grains.

Photometric redshifts (photo-$z$'s) are crucial for the cosmology, galaxy evolution, and transient science drivers of next-generation imaging facilities like the Euclid Mission, the Rubin Observatory, and the Nancy Grace Roman Space Telescope. Previous work has shown that image-based deep learning photo-$z$ methods produce smaller scatter than photometry-based classical machine learning (ML) methods on the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample, a testbed photo-$z$ dataset. However, global assessments can obscure local trends. To explore this possibility, we used a self-organizing map (SOM) to cluster SDSS galaxies based on their $ugriz$ colors. Deep learning methods achieve lower photo-$z$ scatter than classical ML methods for all SOM cells. The fractional reduction in scatter is roughly constant across most of color space with the exception of the most bulge-dominated and reddest cells where it is smaller in magnitude. Interestingly, classical ML photo-$z$'s suffer from a significant color-dependent attenuation bias, where photo-$z$'s for galaxies within a SOM cell are systematically biased towards the cell's mean spectroscopic redshift and away from extreme values, which is not readily apparent when all objects are considered. In contrast, deep learning photo-$z$'s suffer from very little color-dependent attenuation bias. The increased attenuation bias for classical ML photo-$z$ methods is the primary reason why they exhibit larger scatter than deep learning methods. This difference can be explained by the deep learning methods weighting redshift information from the individual pixels of a galaxy image more optimally than integrated photometry.

We introduce NewCluster, a new high-resolution cluster simulation designed to serve as the massive halo counterpart of the modern cosmological galaxy evolution framework. The zoom-in simulation targets a volume of $4.1\sigma$ overdensity region, which is expected to evolve into a galaxy cluster with a virial mass of $5 \times 10^{14} M_\odot$, comparable to that of the Virgo Cluster. The zoom-in volume extends out to 3.5 virial radii from the central halo. The novelties of NewCluster are found in its resolutions. Its stellar mass resolution of $2 \times 10^{4} M_\odot$ is effective for tracing the early assembly of massive galaxies as well as the formation of dwarf galaxies. The spatial resolution of 68 parsecs in the best-resolved regions in the adaptive-mesh-refinement approach is powerful to study the detailed kinematic structure of galaxies. The time interval between snapshots is also exceptionally short-15 Myr-which is ideal for monitoring changes in the physical properties of galaxies, particularly during their orbital motion within a larger halo. The simulation has up-to-date feedback schemes for supernovae and active galactic nuclei. The chemical evolution is calculated for ten elements, along with dust calculation that includes the formation, size change, and destruction. To overcome the limitations of the Eulerian approach used for gas dynamics in this study, we employ Monte Carlo-based tracer particles in NewCluster, enabling a wide range of scientific investigations.

In this paper we address two major questions raised by recent James Webb Space Telescope observations of the young Universe, namely: 1) what are the seed initial masses, and how rapidly have supermassive black holes (BHs) with masses of 1e6-1e8Msun grown in active galactic nuclei (AGN) hosted by very young galaxies? 2) What are the plausible explanations for the super solar abundances of nitrogen in a fraction of young galaxies at high redshift, both with and without evidence of a massive central black hole? We focus mainly on the system GS3073. This system shows an exceptionally large log(N/O)=+0.42(+0.13/-0.10) in the gas close to the AGN. We show here that this abundance is consistent with the composition of gas ejected from massive asymptotic giant branch stars. Moreover, this system shows chemical properties matching those expected at a specific point of the evolution of the abundances in the extreme populations of the former nuclear star cluster wCentauri (wCen). This analogy, along with the N/O, C/O and Fe/O abundances in GS3073, lead to an estimate of an age range of 270-440 Myr for this object, much smaller than the redshift (z=5.5) age of about 1 Gyr. We also adopt the same criteria to estimate an age for GNz11. These two determinations constrain the BH mass versus age relation: accretion on the BH must proceed at intermittent superEddington rates in the first phases, and at a much lower rate after the first half gigayear of life of the Universe. The intermittency of accretion is also a fundamental requirement to allow the formation of the extreme (N rich, O depleted, He rich) populations today observed in wCen for a large range of metallicities.

Multimessenger astronomy seeks to uncover the origins of cosmic rays and neutrinos. The IceCube Neutrino Observatory plays a key role in monitoring the sky for revealing high energy neutrinos and neutrino time clusters possibly associated with astrophysical sources, issuing alerts to the astrophysical community for significant excesses. This enables joint observations with other astronomical facilities that could reveal the hidden mechanisms behind the most extreme environments in the Universe. In particular, since 2006 the Gamma-ray Follow-Up (GFU) program shares cluster alerts with partner Imaging Air Cherenkov Telescopes. The faint cosmic signals, searched against large atmospheric backgrounds, are widely masked by the statistical penalties that arise when scanning the full sky in an unbiased way. Hence, targeted analyses of pre-selected neutrino source candidates have proven to increase our search sensitivity. Our understanding of astrophysical environments has improved in recent years, with evidence of neutrino emission from the blazar TXS 0506+056 and the Seyfert galaxy NGC 1068. The aim of expanding observational possibilities and engaging the broader scientific community through public cluster alerts has motivated the creation of a new list of target sources to be monitored by IceCube. This contribution presents the systematic compilation of this list, which extends the well-established focus on gamma-ray bright active galactic nuclei (AGN) to include X-ray bright AGN and binary systems.

The recent detection of fast-moving stars in the core of Omega Centauri ($\omega$ Cen), the most massive globular cluster (GC) in the Milky Way, has provided strong evidence for the presence of an intermediate-mass black hole (IMBH). As $\omega$ Cen, is likely the accreted nucleus of a dwarf galaxy, this IMBH also represents a unique opportunity to study BH seeding mechanisms and their potential role in the formation of supermassive BHs. We present Monte Carlo $N$-body models of $\omega$ Cen with detailed treatments for the loss cone dynamics involving stars, binaries, and compact objects. Starting with BH seeds of $500-5000 \, M_{\odot}$ (consistent with runaway collisions of massive stars), our cluster models grow IMBHs with masses of $\sim50{,}000 \, M_{\odot}$ after 12 Gyr, while successfully reproducing the present-day surface brightness and velocity dispersion profiles of $\omega$ Cen. We find a population of fast stars similar to those observed in the core of $\omega$ Cen, with the fastest stars originating from binaries that were tidally disrupted by the IMBH. The IMBH growth is primarily driven by mergers with $30-40 \, M_{\odot}$ BHs, suggesting a present-day IMBH-BH merger rate of $\sim(4-8)\times10^{-8}~\rm{yr}^{-1}$ in $\omega$ Cen-like GCs. Our models also predict a similar rate of tidal disruption events ($\sim5\times10^{-8}~\rm{yr}^{-1}$) which, depending on the frequency of $\omega$ Cen-like GCs per galaxy, may represent anywhere from $0.1\%$ to $10\%$ of the observed TDE rate.

The IceCube Neutrino Observatory is a cubic-kilometer Cherenkov array deployed in the deep, glacial ice at the geographic South Pole. An important feature of the instrumented ice are undulations of layers of constant optical properties over the footprint of the detector. During detector construction, these layers were mapped using stratigraphy measurements obtained from a stand-alone laser dust logger. While this system is very precise, its cost does not scale to the instrumented volume envisioned for the proposed IceCube-Gen2 Observatory. Here, we explore the possibility of obtaining equivalent stratigraphy data from camera footage recorded during the deployment of IceCube more than a decade ago. If successful, this could be an alternative technique to be considered for IceCube-Gen2.

Luminous and Ultra-luminous IR galaxies ((U)LIRGs) are critical for investigating feedback mechanisms due to a combination of intense star formation (SF) episodes and active galactic nuclei (AGN), particularly in the context of complex galaxy interactions. We conduct a detailed analysis of the II ZW 096 merging system using the Multi-Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope (VLT), combining high-resolution Narrow Field Mode (NFM) and large-area Wide Field Mode (WFM) observations. We mapped the morphology, kinematics, and ionizing radiation of the system's gas by fitting atomic emission lines and the optical continuum. We identify three or more distinct galaxies within II ZW 096, revealing rotational patterns and complex interactions consistent with a collapsing small galaxy group. The kinematics and ionization structures suggest high star formation rates and shock-driven processes, which align with this proposed scenario. Focusing on the D1 compact region, which contributes 40-70% of the system's IR emission, and combining information from archival multi-wavelength observations, we find strong evidence of a heavily obscured AGN powering it. Our analysis of the internal structure, interactions, and merger state of II ZW 096 offers novel insights into the galaxy evolution processes in this dynamic and highly chaotic system

Understanding the nature and evolution of dark energy (DE) is a central challenge in modern cosmology. In this work, we explore the constraining power of bright standard sirens -- gravitational wave (GW) events with electromagnetic counterparts - for probing the DE equation of state as function of redshift. Focusing on future GW observations from next-generation ground-based GW detectors such as the Einstein Telescope and Cosmic Explorer, we perform a comprehensive analysis using simulated binary neutron star (BNS) and neutron star-black hole (NSBH) events over five years of observation with a $75\%$ duty cycle. We consider three broad classes of DE models: (i) phenomenological parametrizations, specifically the Barboza-Alcaniz extension to the Chevallier-Polarski-Linder model; (ii) physically motivated scalar field scenarios, specifically hilltop quintessence; and (iii) evolving dark matter setup in which the matter density evolves as $(1+z)^{3+\alpha}$. For each case, we jointly infer the Hubble constant $H_0$ and model-specific DE parameters from the observed GW luminosity distances and spectroscopic redshifts. Our results demonstrate that bright sirens alone can yield competitive and independent constraints on the time evolution of DE indicating that multi-messenger cosmology has the potential to test a wide range of DE theories, bridging phenomenological and physically motivated models, and paving the way for precision cosmology in the era of GW astronomy.

IceCube is a neutrino observatory located at the South Pole that uses Antarctic ice as a medium for detection of Cherenkov photons. As such, analysis of the data relies on our understanding of the properties of ice within and around the instrumented volume. Over the years we have made significant progress in understanding the glacial ice and now have a comprehensive model that covers many of the relevant aspects of the photon propagation in it. In this report we give a historical overview of the ice description within the IceCube detector, list some of the remaining issues, and assess how much more improvement is still needed. As the IceCube Upgrade is expected to be installed in less than a year, with several new types of calibration devices aiming to further our understanding of ice, this is the perfect time to review the current state of the ice model.

The central 5 parsecs of our Galactic Center is rich in various molecular tracers. The region is observed to be hot, highly ionized, and volatile, resulting in complex chemistry and kinematics. Countless molecular observations of the central region orbiting the central super-massive black hole Sgr A* have revealed a stable grouping of clouds with a disk-like structure known as the circumnuclear disk (CND) orbiting Sgr A* between 1.5 and 3 pc. The CND houses the densest population of gas outside the ionized streams in the minispiral falling onto Sgr A*. However, the small-scale kinematics of the clouds and gas are poorly understood in this region. In this paper, we report single dish SiO (1-0) observations of the central 5 pc surrounding Sgr A* taken with the Green Bank Telescope. Most of the SiO (1-0) traces the southern streamer and Sgr A East-neither of which are associated with the CND. However, one region in the CND, known as the NE Arm, shows high-intensity SiO emission unlike in the rest of the CND, indicating the NE Arm as a candidate region of either cloud-cloud collision or early star formation. Further high resolution observations of the region are needed to distinguish between the two scenarios.

We present Kurucz-a1, a physics-informed neural network (PINN) that emulates 1D stellar atmosphere models under Local Thermodynamic Equilibrium (LTE), addressing a critical bottleneck in differentiable stellar spectroscopy. By incorporating hydrostatic equilibrium as a physical constraint during training, Kurucz-a1 creates a differentiable atmospheric structure solver that maintains physical consistency while achieving computational efficiency. Kurucz-a1 can achieve superior hydrostatic equilibrium and more consistent with the solar observed spectra compared to ATLAS-12 itself, demonstrating the advantages of modern optimization techniques. Combined with modern differentiable radiative transfer codes, this approach enables data-driven optimization of universal physical parameters across diverse stellar populations-a capability essential for next-generation stellar astrophysics.

Studying galaxy evolution requires knowledge not only of the stellar properties, but also of the interstellar medium (in particular the molecular phase) out of which stars form, using a statistically significant and unbiased sample of galaxies. To this end, we introduce here the integrated Extragalactic Database for Galaxy Evolution (iEDGE), a collection of integrated stellar and nebular emission lines, and molecular gas properties from 643 galaxies in the local Universe. These galaxies are drawn from the CALIFA datasets, and are followed up in CO lines by the APEX, CARMA, and ACA telescopes. As this database is assembled from data coming from a heterogeneous set of telescopes (including IFU optical data and single-dish and interferometric CO data), we adopted a series of techniques (tapering, spatial and spectral smoothing, and aperture correction) to homogenise the data. Due to the application of these techniques, the database contains measurements from the inner regions of the galaxies and for the full galaxy extent. We used the database to study the fundamental star formation relationships between star formation rate (SFR), stellar mass ($M_*$), and molecular gas mass ($M_{\rm mol}$) across galaxies with different morphologies. We observed that the diagrams defined by these quantities are bi-modal, with early-type passive objects well separated from spiral star-forming galaxies. Additionally, while the molecular gas fraction ($f_{\rm mol}=M_{\rm mol}/M_*$) decreases homogeneously across these two types of galaxies, the star formation efficiency (SFE=SFR/$M_{\rm mol}$) in the inner regions of passive galaxies is almost two orders of magnitude lower compared to the global values. This indicates that inside-out quenching requires not only low $f_{\rm mol}$, but also strongly reduced SFE in the galactic centres.

The identification of X-ray and CMB sources as galaxy groups and clusters is a prerequisite for cluster cosmology. But the identification of groups, especially nearby ones, suffers from projection effects which in turn affect the purity of the sample. In X-rays, the position of the cluster can be given either by the peak of the emission, or by the full information content of the cluster image. Similarly, the optical center, or its member galaxies, can describe the optical counterpart. With the progress of numerical simulations, it is currently feasible to reproduce both the optical group membership assignment and the behavior of the group outskirts in X-rays, and therefore there is an opportunity to define a reproducible group identification procedure. We performed two-way matching between X-ray contours, drawn at a fixed surface brightness level corresponding to a baryonic overdensity of ~500, to the projected galaxy positions using a modified Hausdorff distance (MHD). We used the volume-limited SDSS group catalog to evaluate the purity and completeness of the procedure, maintaining the constant performance of the optical group finder with redshift. We find that an MHD of 0.631 Mpc provides 90% purity. This is a clear improvement over the methods that rely on the optical counterparts' distance and richness. We study the purity versus MHD and the completeness versus redshift and velocity dispersion. Over half of nearby groups have X-ray emission, even those with velocity dispersions as low as 200km/s; this has never previously been demonstrated. The bulk of these groups follow the same scaling relations as the groups with a small separation between the optical and X-ray centers, removing feedback as an explanation for the lack of matches in previous studies. Instead, the problem is caused by over-merging in the optical group catalog construction and source confusion in X-rays.

As most massive stars are born in binary and other multiple-star systems, many are expected to exchange mass with a companion star or merge with it during their lives. This means that most supernovae (SNe) are from such binary products. Here, we focus on hydrogen-rich Type II SNe from accretors of binary mass transfer and stellar mergers. We compute various SN properties such as the explosion energies, nickel yields, and neutron star (NS) kick velocities, but also consider NS masses. We find tight correlations between these parameters and, e.g., the central specific entropy and core compactness. However, there is no obvious relation between these explosion properties and the evolutionary history of the pre-SN stars. We find linear relations between the nickel mass and the SN explosion energy and the NS remnant mass. We further group our models into progenitors of SNe IIP, SN 1987A-like and interacting SNe, predict their SN and SN-progenitor properties and compare to observations. Accretors of binary mass transfer and stellar mergers naturally produce SNe IIP with long plateau durations from progenitors with relatively small CO-cores but large envelope masses (c.f. SN 2015ba). Our models give rise to tight relations between the plateau luminosity and the nickel mass as well as the SN ejecta velocity as inferred observationally for SNe IIP. We speculate that cool/red supergiants at $\log\,L/L_\odot\,{\geq}\,5.5$ encounter enhanced mass loss due to envelope instabilities and could then explode in interacting SNe IIn. The rate of such SNe from our models seems compatible with observations. Some of our binary models explode as $10^6\,L_\odot$ blue supergiants that may have encountered enhanced and/or eruptive mass loss shortly before their SNe and could thus help understand interacting SNe such as SN 1961V and SN 2005gl but also superluminous Type II SNe such as SN 2010jl. [abridged]

Galaxy evolution is largely driven by star formation activity or by the cessation of it, also called star formation quenching. In this paper, we present star formation scaling relations for galaxies at different evolutionary stages. To do so, we used the integrated Extragalactic Database for Galaxy Evolution (iEDGE), which collects CO, optical continuum, and emission line information for 643 galaxies from the CALIFA IFU dataset. By considering the patterns described by star-forming and retired regions, we grouped the galaxies into quenching stages using the emission line classification scheme, QueStNA. We observed that the molecular gas mass ($M_{\rm mol}$) decreases from star-forming to retired systems and so does the molecular-to-stellar mass ratio ($f_{\rm mol}$). In contrast, star formation efficiency (SFE) is constant in the quenching stages dominated by star formation and rapidly declines afterwards. We observed that this rapid decline is more pronounced in the centre of the galaxies compared to the rest of the discs, reflecting the inside-out quenching displayed by nearby galaxies. We noticed that the relations between $M_{\rm mol}$ and the stellar mass ($M_*$) become increasingly shallow with the quenching stages; however, the relations between the star formation rate (SFR) and $M_{\rm mol}$ steepen when going from star-forming to retired systems. We observed that a three-dimensional relation between SFR, $M_*$, and $M_{\rm mol}$ exists for star-forming galaxies, while data points from other quenching groups are scattered across the parameter space. Taken together, these pieces of evidence indicate that the quenching of the galaxies cannot be explained solely by a depletion of the molecular gas and that a significant decrease in the SFE is necessary to retire the centre of the galaxies beyond the star formation green valley.

Thermonuclear supernovae (SNe) are the result of the nuclear transformation of carbon/oxygen (C/O) white dwarfs (WDs) to the radioactive element $^{56}\mathrm{Ni}$ and intermediate mass elements (IMEs) like Ca, Ar, etc. Most progenitor scenarios involve a companion star which donates matter to the exploding white dwarf, implying a fundamental prediction: the formation of a wake in the explosive ejecta as it runs into and moves past the companion star. This wake leaves an indelible imprint on the ejecta's density, velocity, and composition structure that remains fixed as the ejecta reaches homologous expansion. We simulate the interaction of the ejecta and Roche-lobe filling donor in a double degenerate double detonation Type Ia progenitor scenario and explore the detectability of this imprint in late-time nebular phase spectroscopy of Type Ia SNe under the assumption of local heating ($t > 200$ days). At these times, the velocity profiles of forbidden emission lines reflect the velocity distribution of all of the ejecta and the critical electron density for that forbidden line. We explicitly calculate line shapes for the [Co III] $11.89 \mu\mathrm{m}$ line that traces the initial $^{56}\mathrm{Ni}$ distribution and the [Ar III] $8.99 \mu\mathrm{m}$ line, which traces a typical intermediate mass element. We predict the viewing angle dependence of the line shape, present a tool to quickly calculate optically thin line shapes for various 3D density-velocity profiles and discuss JWST observations.

Only 40 exoplanetary systems with five or more planets are currently known. These systems are crucial for our understanding of planet formation and planet-planet interaction. The M dwarf L 98-59 has previously been found to show evidence of five planets, three of which are transiting. Our aim is to confirm the fifth planet in this system and to refine the system characteristics namely minimum masses, radii and the orbital parameters of the planets around L 98-59. We reanalysed RV and activity data from HARPS and ESPRESSO alongside TESS and HST transit data using a joint model. The parameter space was sampled using the dynesty nested sampler. We confirm the previously known fifth planet in the system's habitable zone with an orbital period of $23.07\pm0.08\,d$, a minimum mass of $3.0\pm0.5\,M_{\oplus}$ and an effective temperature of 289 K. We find an additional planet candidate in the RV data with an orbital period of $1.7361^{+0.0007}_{-0.0008}\,d$ and a minimum mass of $0.58\pm0.12\,M_{\oplus}$. This candidate (L 98-59.06) has a statistical significance between $2.9\sigma$ and $4.2\sigma$, details depending on the modelling of stellar variability. Moreover, we present evidence for a stellar rotation period of $76\pm4\,d$.

Gravitational lensing has transformed the field of gas tomography in the intergalactic medium (IGM) and circumgalactic medium (CGM). Here we use the brightest lensed galaxy identified to date, the Sunburst Arc ($z$$\approx$2.37), to constrain the physical size of foreground absorbers at $z$$\approx$2 in 2D. This galaxy is a confirmed Lyman continuum leaker, where its single leaking region is imaged 12 times over four separate arcs. The separations between the arcs allows for large scale tomography, while the distances between the images along an arc allow for small scale tomography. Using HST/WFC3 UVIS G280 grism observations, we extracted the spectra of the leaking region and fit for absorbers detected along these lines of sight using a binary population and spectral synthesis (BPASS) model for the galaxy. We identified two partial Lyman limit systems (pLLSs) and one Lyman limit system (LLS) across the different spectra and measured their physical sizes. We find consistent HI column densities across $\lesssim$2 kpc and an average HI mass of $\approx$10$^3$ ${\rm M}_\odot$ for the absorbers. Given the strong CIV lines associated with two of the absorbers, they are likely located within the CGM of foreground galaxies. The third absorber has no associated metal lines, so it is most likely within the IGM. This study provides the first tomography measurements of pLLSs/LLSs in the CGM and IGM at $z$$\approx$2.

Recent studies suggest that tidal disruption events (TDEs) with off-axis jets may manifest as optically overluminous events. To search for jet signatures at late times, we conducted radio observations of eight such optically overluminous ($M_{g, \rm peak} < -20.8$ mag) TDEs with the Very Large Array. We detect radio counterparts in four events. The observed radio luminosities ($L_{\rm 6 GHz} \sim 10^{38}$--$10^{39}$ erg s$^{-1}$) are two orders of magnitude lower than those of on-axis jetted TDEs, and we find no evidence for off-axis jets within rest-frame time of 3 yrs. Two of them (AT2022hvp and AT2021aeou) exhibit evolving radio emission, consistent with synchrotron emission from non-relativistic outflows launched near the time of first optical light. Two events (AT2020ysg and AT2020qhs) show no statistically significant variability, which can be attributed to either non-relativistic outflows or pre-existing active galactic nuclei. Compared to a control sample of fainter TDEs with $M_{g, \rm peak} > -20.5$ mag observed at similar rest-frame timescales ($t_{\rm rest} ~ 1.5$\,yr), our sample shows systematically more luminous radio emission, suggesting that optically overluminous TDEs may launch more powerful prompt non-relativistic outflows. We speculate that strong general relativistic effects near high-mass black holes ($M_{\rm BH} ~ 10^8\,M_\odot$) may play a key role. These findings motivate further investigation into the nature of relativistic disruptions around massive black holes and the physical conditions necessary for jet formation.

V1674 Her is one of the fastest and brightest novae, characterized by dense optical photometry in the pre-maximum phase, a rise from $g=17$ to 7 mag, in one-fourth of a day. We present a composite theoretical $V$ light curve model of its early rising phase starting from a quiescent brightness of $g=19.2$ mag. Our light curve model consists of a hot and bright white dwarf (WD) and irradiated accretion disk and companion star. We found that the earliest optical detection of ASAS-SN $g$ band brightness of $g=17.0$ at $t=0.014$ day from the onset of thermonuclear runaway can be explained with the irradiated accretion disk and companion star in the X-ray flash phase of a $1.35 ~M_\odot$ WD. This is the first detection in optical of an X-ray flash phase of a nova. Optically thick winds emerge from the WD photosphere at $t=0.04$ day, and optical flux is dominated by free-free emission from optically-thin ejecta just outside the WD photosphere. Our free-free emission model $V$ light curve reasonably reproduces the dense $g$ light curve of Evryscope that spans from $g=14.8$ (at 0.078 day) to $g=7.1$ (at 0.279 day), including a sudden change of slope in the $g$ light curve from slow to rapid rise at $g=14.3$ on day $0.1$. There is no indication of shocking power during the rising phase from $g=14.8$ to 7.1.

Classical Cepheids (CCs) are important probes for the large-scale warp structure of the Milky Way. Using Gaia DR3 CCs, we establish an optimal time-dependent warp model, where the warp height increases with radius following a power-law, the line of nodes (LONs) exhibit linear twisting with radius, following a leading spiral pattern, and the LONs undergo prograde evolution over time. Structurally, we identify significant warp features in the $5-9$ kpc region of the Galactic disk, where the warp model performs better than the flat model. Beyond 15 kpc, the model with the second Fourier term does not fit the observations well, whereas the model with twisted LONs better matches the data. Kinematically, we derived expressions for the vertical velocities using direct differentiation and then calculated the precession rates for each CC. Our results intuitively indicate a nearly uniform and low warp precession rate of $\omega = 4.86 \pm (0.88)_{stat} \pm (2.14)_{sys}$ km s$^{-1}$ kpc$^{-1}$ beyond 12.5 kpc, in agreement with classical kinematic estimates. Based on these findings, we propose a simple yet comprehensive time-dependent warp model, $Z_{w}(t) = 0.00019R^{3.08}\sin(\phi - (3.87R-41.79 + 4.86t))$, which provides a unified framework for describing both the geometric and kinematic evolution of the Galactic warp. We analyzed the impact of the adopted solar vertical velocity on the inferred warp precession rate and confirmed the reliability of the measured precession rate. In addition, we found that extinction treatment affects the warp amplitude in the inner disk, while its influence on the outer disk warp structure and the precession rate is negligible.

We investigate Power-Law Plateau (PLP) inflation in standard gravity and its consistency with ACT DR6 data. While many inflationary models, including the Starobinsky inflation, are disfavored by ACT observations, the PLP potential remains viable across a broad range of its parameters. Then, the dynamics of the reheating phase are investigated, where we mainly focus on the reheating temperature and its relationship with the inflationary phase and primordial gravitational waves. Incorporating the overproduction of the primordial gravitational waves can affect the effective number of relativistic species during the bounce. The constraint data on $\Delta N_{\rm eff}$ can impose a lower bound on the reheating temperature. This constraint will be more efficient for a stiff equation of state. It is determined that for $\omega_{re} > 0.58$, this constraint would be efficient. Combining the result of the reheating temperature and the inflationary phase, it is concluded that to have both a viable result standing in $1\sigma$ of ACT DR6 and also to satisfy the reheating lower bound, the total number of e-folds during the inflationary phase should be $N_k \lesssim 62$. Higher e-folds of expansion result in a reheating temperature below the bound, which is disfavored. Finally, for the constraint values of the reheating temperature, the energy spectrum of the gravitational waves has been explored. The results indicate that there is a higher chance of detection for lower reheating temperatures and higher reheating equation of state.

We probe the atomic hydrogen (HI) emission from the host galaxies of fast radio bursts (FRBs) to investigate the emerging trend of disturbance and asymmetry in the population. Quadrupling the sample size, we detect 13 of 14 new hosts in HI, with the only non-detection arising in a galaxy known to be transitioning towards quiescence. With respect to typical local Universe galaxies, FRB hosts are generally massive in HI ($M_{HI}>10^9 M_\odot$), which aligns with previous studies showing that FRB hosts also tend to have high stellar masses and are star-forming. However, they span a broad range of other HI derived properties. In our independent sample of repeater hosts, we observe a statistically insignificant preference towards lower HI masses compared to non-repeater hosts, similar to the low-significance trend toward lower stellar masses previously reported. Using visual inspection alongside various asymmetry metrics, we identify four unambiguously settled host galaxies, demonstrating for the first time that a disturbed HI morphology is not a universal feature of FRB host galaxies. However, we find another six that show clear signs of disturbance, and three which require deeper, more targeted observations to reach a conclusion; this brings the confirmed ratio of disturbed-to-settled FRB hosts to 11:4. Given that roughly a 1:1 ratio is expected for random background galaxies of similar type, our observed ratio yields a p-value of 0.065. Unlike earlier indications based on smaller samples, this no longer crosses the conventional threshold for statistical significance, though is still near enough to hint at a legitimate excess of disturbance among FRB hosts. Thus, an even larger sample size of FRB hosts observed in HI is required to fully clarify whether the trend is genuine or still a consequence of low-number statistics - a sample that upcoming data releases are well positioned to provide.

The detection of transient phenomena such as Gamma-Ray Bursts (GRBs), Fast Radio Bursts (FRBs), stellar flares, novae, and supernovae, alongside novel cosmic messengers like high-energy neutrinos and gravitational waves, has transformed astrophysics in recent years. Maximizing the discovery potential of multi-messenger and multi-wavelength follow-up observations, as well as serendipitous detections, requires a tool that rapidly compiles and contextualizes relevant information for each new event. We present Astro-COLIBRI, an advanced platform designed to meet this challenge. Astro-COLIBRI integrates a public RESTful API, real-time databases, a cloud-based alert system, and user-friendly clients (a website and mobile apps for iOS and Android). It processes astronomical alerts from multiple streams in real time, filtering them based on user-defined criteria and placing them in their multi-wavelength and multi-messenger context. The platform offers intuitive data visualization, a quick summary of relevant event properties, and an assessment of observing conditions at numerous observatories worldwide. We here describe its architecture, data resources, and main functionalities. We highlight the automatic collection of photometric data from a variety of large scale optical surveys, a recently added feature that significantly improves the capabilities of the Astro-COLIBRI platform.

The Dark MAtter Particle Explorer (DAMPE) instrument is a space-born cosmic-ray detector, capable of measuring ion fluxes up to $\sim$500 TeV/n. This energy scale is made accessible through its calorimeter, which is the deepest currently operating in orbit. Saturation of the calorimeter readout channels start occurring above $\sim$100 TeV of incident energy, and can significantly affect the primary energy reconstruction. Different techniques, analytical and machine-learning based, have been developed to tackle this issue, focusing on the recovery of single-bar deposits, up to several hundreds of TeV. In this work, a new machine-learning technique is presented, which profit of a unique model to correct the total deposited energy in DAMPE calorimeter. The described method is able to generalise its corrections for different ions and extend the maximum detectable incident energy to the PeV scale.

Observations of our Galaxy at very-high and ultra-high photon energies have revealed a rich population of sources, many of which have significant angular extension and/or are positionally coincident with energetic pulsars. We assessed whether the properties of these gamma-ray emitters are consistent with the predictions of a model for the non-thermal leptonic emission of a pulsar wind nebula expanding within a supernova remnant, including the possibility of particle escape across the different components of the system up to the ambient medium. We used a multi-zone model framework that describes the transport and radiation of a population of ultra-relativistic electron-positron pairs. The dynamics of the system is retrieved from numerical hydrodynamics simulations, which allows to follow the evolution through advanced stages. Model predictions for typical parameter setups are compared to source catalogs issued by the H.E.S.S. and LHAASO collaborations. If particles remain confined in the nebula, 1-100TeV emission is only partially consistent with the properties of observed sources in terms of surface brightness, angular extent, and photon index. In particular, the systems spend a long time in a bright and soft state that has no equivalent in observations. Conversely, including the possibility of energy-dependent particle escape across the object results in a much better agreement, both in the 1-10TeV range explored with H.E.S.S. or LHAASO-WCDA, and in the 20-100TeV range explored by LHAASO-KM2A. Emission components with intermediate to low brightnesses and large to very large sizes result from particles trapped in the remnant or escaped in the surrounding medium. Particle escape, clearly seen in the form of misaligned jets or halos around middle-aged pulsars, is a very important process in much younger systems with 1-100kyr ages, and shapes their appearance as TeV/PeV sources.

GRANDProto300 is one of the prototype experiments of the Giant Radio Array for Neutrino Detection. It will feature about 300 radio antenna detectors in Xiaodushang in Dunhuang, China, covering a total geometrical area of about 200 km$^2$. A main scientific goal of GRANDProto300 is the study of cosmic rays in the transition region ($10^{17}\, {\rm eV} < E < 10^{18.5}\, {\rm eV}$). Our study calculates the exposure of GRANDProto300 to cosmic rays and estimates the number of cosmic-ray events to be detected during a fixed observation period. The trigger efficiency reaches 50, 80, and 90% at $10^{17.5}$, $10^{17.9}$, and $10^{18.3}\, {\rm eV}$, respectively. The exposure of GRANDProto300 is 50 km$^2$ day sr at around $10^{17.5}$ eV, and the expected number of observed cosmic rays with energies above $10^{17}$ eV and zenith angles above 65$^{\circ}$ is about 130 events per day. GRANDProto300 will be able to measure the cosmic-ray energy spectrum in $10^{17.2}\, {\rm eV} < E < 10^{19.5}\, {\rm eV}$ through one-year observation, with a statistical precision about five times better than the previous spectral measurement by a mono-fluorescence detector of the Telescope Array Low-Energy Extension. The statistical uncertainty in the measurement of the mean depth of the air-shower maximum $X_{\rm max}$ is about five times better than the previous measurements using radio detectors at $10^{17.5}$ eV; systematic uncertainties should be a dominant contribution limiting our interpretation of the chemical composition of cosmic rays in the transition region.

Characterizing the astrophysical neutrino flux with the IceCube Neutrino Observatory traditionally relies on a binned forward-folding likelihood approach. Insufficient Monte Carlo (MC) statistics in each bin limits the granularity and dimensionality of the binning scheme. A neural network can be employed to optimize a summary statistic that serves as the input for data analysis, yielding the best possible outcomes. This end-to-end optimized summary statistic allows for the inclusion of more observables while maintaining adequate MC statistics per bin. This work will detail the application of end-to-end optimized summary statistics in analyzing and characterizing the galactic neutrino flux, achieving improved resolution in the likelihood contours for selected signal parameters and models.

The dehydrogenated cation of methylated benzene is known to exist in two isomeric forms: benzylium and tropylium. Structurally similar forms have been proposed for the -H cations of methylated polycyclic aromatic hydrocarbons, but their spectroscopic characterization remains limited, and their photophysical properties are still poorly understood. Previous studies identified 2-naphthylmethylium and benzyltropylium as specific long-lived isomers of the -H fragment of methylnaphthalene cations. Here, we investigate the photodissociation spectroscopy and photoprocessing of gas-phase C$_{11}$H$_9^+$ ions in the visible range. Experiments are conducted using the versatile laboratory astrophysics setup PIRENEA, which enables studies over long timescales ($\sim$1000 s) and allows photoprocessing to be combined with ion-molecule reaction experiments. We confirm the presence of 2-naphthylmethylium and benzyltropylium and additionally identify 1-naphthylmethylium, previously undetected in experiments. Moreover, we present the first complete quantitative analysis of the relative abundances of these isomers and provide clear evidence of interconversion among the three long-lived species. These isomerization processes occur below the dissociation threshold, a finding supported by molecular dynamics simulations. The photophysical properties of C$_{11}$H$_9^+$ isomers -- including isomerization and fluorescence -- make them intriguing candidates for consideration in astrophysical environments exposed to mild UV irradiation (h$\nu$ < 7 eV). Moreover, their detection via rotational spectroscopy in such regions is facilitated by their closed-shell electronic structure.

Models of astrophysical convection, such as mixing length theory, typically assume that the heat transport is independent of microphysical diffusivities. Such 'diffusion-free' behaviour is, however, not observed in numerical simulations employing standard fixed-flux or fixed-temperature boundary conditions, except possibly in extreme parameter regimes that are computationally expensive to achieve. Recent numerical and experimental work has suggested that internally heated and cooled convection can exhibit diffusion-free scalings in more numerically accessible regimes. Here, we present direct numerical simulations of 2.5D Cartesian rotating thermal convection driven by an internal heating and cooling function. The use of distributed heating and cooling functions alleviates sharp thermal boundary layers that would otherwise be present, allowing the flows to be simulated with modest computational resources. We show that for high Rossby numbers this set-up recovers mixing length theory scalings for the heat transport. The velocity amplitudes, in contrast, are observed to display diffusion-limited scalings. By comparing against boundary driven rotating convection, we show that internally heated cases have a larger fraction of their thermal dissipation occuring in the bulk of the fluid. We suggest this is connected to the increased convective efficiency observed in these cases. Our results indicate that 2.5D internally heated convection can be used as a computationally inexpensive test-bed to investigate some aspects of diffusion-free heat transport.

We explore the evolution of a giant planet that interacts with a circumbinary disc that orbits a misaligned binary by means of analytic models and hydrodynamical simulations. Planet-disc interactions lead to mutual tilt oscillations between the planet and the disc. Even if circumbinary gas discs form with an isotropic mutual misalignment to the binary, planet-disc interactions can cause giant planets to evolve towards coplanar or polar alignment. For a low-mass disc, the binary dominates the dynamical evolution of the planet leading to a wide range of circumbinary planet inclinations. For a high-mass disc, the disc dominates the dynamical evolution of the planet and planet inclinations move towards coplanar or polar alignment to the binary orbit, depending upon the initial disc inclination and the binary eccentricity. In addition, for a high-mass disc ($\sim 50\, M_{\rm p}$) and a high initial disc inclination, the planet can undergo Kozai-Lidov oscillations that can result in the planet being ejected from the system. For initially highly misaligned systems, the non-coplanarity of the planet and the disc can lead to long-lived inner misaligned disc rings that can become highly eccentric.

The Giant Radio Array for Neutrino Detection (GRAND) aims to detect radio signals from extensive air showers (EAS) caused by ultra-high-energy (UHE) cosmic particles. Galactic, hardware-like, and anthropogenic noise are expected to contaminate these signals. To address this problem, we propose training a supervised convolutional network known as an encoder-decoder. This network is used to learn a coded representation of the data and remove specific features from it. This denoiser is trained using high-fidelity air shower simulations specifically tailored to replicate the characteristics of signals detected by GRAND. In this contribution, we describe our machine-learning model and report initial results demonstrating the sensitivity enhancement resulting from our denoising algorithm when applied to realistically simulated GRAND signals with varying signal-to-noise ratios.

Thousands of small bodies, known as trans-Neptunian objects (TNOs), orbit the Sun beyond Neptune. TNOs are remnants of the planets' formation from a disc of gas and dust, so it is puzzling that they move mostly on eccentric orbits inclined to the planetary plane and show a complex red-to-grey colour distribution. A close stellar flyby can account for the TNOs' dynamics, but it is unclear if this can also explain the correlation between their colours and orbital characteristics. Assuming an initial red-to-grey colour gradient in the disc, our numerical study finds that the spiral arms induced by the stellar flyby simultaneously lead to the observed TNOs' colour patterns and orbital dynamics. The combined explanation of these TNO properties strengthens the evidence for a close flyby of another star to the young Solar System. Our study predicts that (1) small TNOs beyond 60 au will mostly be grey, and (2) retrograde TNOs will lack the colour most common to high-inclination TNOs. The anticipated TNO discoveries by the Vera Rubin telescope will be able to test these predictions. A confirmed flyby would allow us to reveal the chemical composition of the Solar System's primordial disc.

GRANDProto300 (GP300) is a prototype array of the GRAND experiment, designed to validate the technique of autonomous radio-detection of astroparticles by detecting cosmic rays with energies between 10$^{17}$-10$^{18.5}$ eV. This observation will further enable the study of the Galactic-to-extragalactic source transition region. Since November 2024, 46 out of 300 antennas have been operational and collecting data stably. We present our cosmic-ray search pipeline, which involves several filtering steps: (1) coincidence search for signals triggering multiple antennas within a time window, (2) directional reconstruction of events, (3) exclusion of clustered (in time and space) noise events, (4) polarization cut, (5) selection based on the size of the footprint, and (6) other more arbitrary cuts in this preliminary stage, including visual cuts. The efficiency of the pipeline is evaluated and applied to the first batch of data, yielding a set of cosmic-ray candidate events, which we present.

We report the detection of a gravitationally lensed galaxy by the nearby spiral galaxy CGCG 012-070 ($z = 0.048$) using Integral Field Unit (IFU) observations with the Near-Infrared Spectrograph (NIRSpec) instrument on board the James Webb Space Telescope (JWST). The lensed galaxy is identified through the flux distributions of emission lines in the rest-frame optical, consistent with a source located at a redshift of $z\sim2.89$. The system is detected in [O III]$\lambda\lambda4959,5007$, H$\beta$, and H$\alpha$ emission lines, exhibiting line ratios typical of a star-forming galaxy. The emission-line flux distributions reveal three distinct components, which are modeled using an elliptical power-law (EPL) mass profile for the lens galaxy. This model provides a good characterization of the source and reveals a disturbed star-forming morphology consistent with those of galaxies at cosmic noon.

The radio emission of cosmic-ray air-showers changes significantly depending on parameters like signal frequency, magnetic field configuration and observing altitude. We use CoREAS simulations to adapt an existing signal model for the radio emission of inclined showers in the 30-80 MHz frequency band to the wide 50-200 MHz band. Our model uses a parametrisation of the charge excess fraction to isolate the geomagnetic emission component. We reconstruct the geomagnetic radiation energy by fitting a lateral distribution function, provided by the model, to the geomagnetic energy fluence distribution of the shower. After we correct for the shower geometry and air density, we correlate the radiation energy with the electromagnetic energy of the shower. We show that the method intrinsic energy resolutions < 5% for the sites of the Pierre Auger Observatory and GRANDProto300. For GRANDProto300, we test the reconstruction with simulations of a realistic, sparse antenna grid and with added noise, and find an energy resolution of < 10% with negligible bias. We do a similar study for a much larger array of 10, 000 km2 with 1 km antenna spacing. We find an intrinsic energy resolution of < 10%.

Galaxy quenching, the intricate process through which galaxies transition from active star-forming states to retired ones, remains a complex phenomenon that requires further investigation. This study investigates the role of active galactic nuclei (AGNs) in regulating star formation by analyzing a sample of 643 nearby galaxies with redshifts between 0.005 and 0.03 from the Calar Alto Legacy Integral Field Area (CALIFA) survey. Galaxies were classified according to the Quenching Stages and Nuclear Activity (QueStNA) scheme, which categorizes them based on their quenching stage and the presence of nuclear activity. We further utilized the integrated Extragalactic Database for Galaxy Evolution (iEDGE), which combined homogenized optical integral field unit and CO observations. This allowed us to examine how AGNs influence the molecular gas reservoirs of active galaxies compared to their non-active counterparts at similar evolutionary stages. Our Kolmogorov-Smirnov and chi-squared tests indicate that the star formation property distributions and scaling relations of AGN hosts are largely consistent with those of non-active galaxies. However, AGN hosts exhibit systematically higher molecular gas masses across all quenching stages except for the quiescent nuclear ring stage. We find that AGN hosts follow the expected trends of non-active quenching galaxies, characterized by a lower star formation efficiency and molecular gas fraction compared to star-forming galaxies. Our results suggest that signatures of instantaneous AGN feedback are not prominent in the global molecular gas and star formation properties of galaxies.

Eclipsing, short-period post-common envelope binaries (PCEBs) have been studied for several decades by eclipse timing variations (ETVs) which have been interpreted as being caused by circumbinary bodies. In this paper we report 355 new observations of 7 PECBs (HS0705+6700, NN Ser, NSVS 07826147, NSVS 14256825, NY Vir, QS Vir and RR Cae) and examine how the recent proposed models of these systems compare with our new observations. We find that none of the 18 recent models fit accurately with our data. We review alternative mechanisms of the ETVs, including magnetic effects, but conclude that they do not predict our results. Although we cannot exclude the presence of circumbinary bodies a combination of several mechanisms may be required to explain the observed ETVs.

We investigate the impact of Lorentz invariance violation (LIV) on radiation processes in astrophysical sources, focusing on synchrotron and inverse Compton interactions. We derive modified expressions for radiated power and photon energy under LIV assumptions and incorporate them into first-order Fermi acceleration models. Our analysis reveals energy thresholds beyond which LIV significantly alters particle dynamics and photon spectra, introducing non-physical divergences that highlight limitations in perturbative approaches. We model synchrotron self-Compton (SSC) emission in the presence of LIV and assess its consequences for photon fluxes from blazars, including Markarian 501 and the BL Lac population. LIV introduces distinct high-energy emission regions that deviate from standard expectations. Comparisons with observational data, particularly upper limits from the Pierre Auger Observatory, suggest that future multi-messenger observations could constrain LIV parameters through the non-detection of such excesses.

We present new UBV photometric observations (2009-2024) of the carbon-rich post-AGB star IRAS 22272+5435 (V354 Lac), combined with our data spanning over 30 years (1991-2024). The star shows pulsations with two closely spaced periods, 127 and 132 d, which cause brightness variations of variable amplitude. We also detect low-amplitude oscillations with a 661-day period, of uncertain origin. Color analysis reveals both temperature variations during pulsations and spectral anomalies. For the first time, we analyze long-term JHKLM photometry (1994-2024), showing increased KLM-band brightness from 1996 to 2004 - likely due to dust emission following a sudden mass ejection. The absence of strong optical variability suggests either a large angle between the ejection direction and the line of sight or a low optical depth of the ejected material. We discuss possible mechanisms and observational signatures of this dust-formation event.

We analyse the stellar mass-metallicity (M-Z) and stellar mass-star formation rate (M-SFR) relations for star-forming galaxies classified by their environment and compare them with matched control samples of field galaxies. Using data from the Galaxy And Mass Assembly (GAMA) survey and the filament catalogue, which categorises galaxies into filaments, tendrils, and voids, we correct emission lines for dust extinction and select star-forming galaxies based on the BPT diagram. Metallicity and star formation rate are estimated and used to fit the M-Z and M-SFR relations through both Bayesian and least-squares approaches. We find that metallicity increases in denser environments, while star formation rate decreases, with the most notable contrasts seen between filament/tendril galaxies and those in voids. Galaxies in filaments and tendrils that are not group members show little to no deviation from their control samples. Morphological analysis reveals no significant differences. Overall, galaxies in denser environments appear more chemically enriched with lower SFRs, likely due to processed material and reduced cold gas availability, while isolated void galaxies maintain higher SFRs and lower metallicities, possibly due to ongoing cold gas accretion. These results suggest that local environmental conditions, rather than large-scale structure alone, are the main drivers of the observed trends.

We present an efficient and accurate algorithm for solving the Poisson equation in spherical polar coordinates with a logarithmic radial grid and open boundary conditions. The method employs a divide-and-conquer strategy, decomposing the computational domain into hierarchical units with varying cell sizes. James's algorithm is used to compute the zero-boundary potentials of lower-level units, which are systematically integrated to reconstruct the zero-boundary potential over the entire domain. These calculations are performed efficiently via matrix-vector operations using various precomputed kernel matrices. The open-boundary potential is then obtained by applying a discrete Green's function to the effective screening density induced at the domain boundaries. The overall algorithm achieves a computational complexity of $\mathcal{O}(N^3 \log N)$, where $N$ denotes the number of cells in one dimension. We implement the solver in the FARGO3D magnetohydrodynamics code and demonstrate its performance and second-order accuracy through a series of test problems.

We present MIRI/JWST medium resolution spectroscopy (MRS) and imaging (MIRIM) of B14-65666, a Lyman-break and interacting galaxy at redshift $z$=7.15. We detect the H$\alpha$ line emission in this system, revealing a spatially-resolved structure of the H$\alpha$ emitting gas, which consists of two distinct galaxies, E and W, at a projected distance of 0.4". Galaxy E is very compact in the rest-frame UV, while W galaxy is more extended, showing a clumpy structure reminiscent of a tidal tail. The total H$\alpha$ luminosity implies that the system is forming stars at a Star Formation Rate (SFR) of 76$\pm$8 M$_{\odot}$ yr$^{-1}$ and 30$\pm$4 M$_{\odot}$ yr$^{-1}$ for E and W, respectively. The ionizing photon production efficiency is within the range measured in galaxies at similar redshifts. The high values derived for the H$\alpha$ equivalent widths (EW) and the distinct locations of the E and W galaxies in the $\log(\zeta_\mathrm{ion}$) $-$ EW (H$\alpha$) plane, indicate that the system is dominated by a young (less than 10 Myr) stellar population. The overall spectral energy distribution suggests that in addition to a young stellar population, the two galaxies may have mature stellar population and very different dust attenuation. The derived SFR and stellar masses identify the two galaxies as going through a starburst phase. The kinematics of the ionized gas traced by the H$\alpha$ line show a velocity difference of 175 $\pm$ 28 km s$^{-1}$ between the two components of B14-65666. The in-depth study of systems like B14-65666 reveal how galaxy mergers in the early Universe drive intense star formation, shape the interstellar medium, and influence the buildup of stellar mass, just 700 $-$ 800 Myr after the Big Bang.

Recently, IceCube has observed an excess of astrophysical neutrinos from the Galactic plane. Such a signal may indicate the presence of individual sources or a diffuse neutrino flux from the interactions of Galactic, hadronic cosmic rays. We consider the prospect of neutrino production in galactic X-ray binary systems. A model for neutrino production in the variable jets of these systems is discussed, and how such information may be incorporated within searches for astrophysical neutrinos from these sources. Two ongoing studies are presented on behalf of the IceCube Collaboration using 13 years of cascades and track-like events. In the first study, recently published model predictions are used to motivate the analysis of five selected black-hole X-ray binaries. Gamma-ray data from Fermi-LAT, as well as X-ray data from Swift/BAT and MAXI are used to search for temporal correlations between radiative cycles of the disk-jet system and potential neutrino emission. In the second study, a wider selection of variable X-ray binaries and their Swift/BAT light curves is considered and a stacking search is used to constrain the contribution of this source population to the Galactic neutrino flux.

We measure three-dimensional bispectra of halo intrinsic alignments (IA) and dark matter overdensities in real space from N-body simulations for halos of mass $10^{12}-10^{12.5} M_\odot /h$. We show that their multipoles with respect to the line of sight can be accurately described by a tree-level perturbation theory model on large scales ($k\lesssim 0.11\,h$/Mpc) at $z=0$. For these scales and in a simulation volume of 1 Gpc/$h$, we detect the bispectrum monopole $B_{\delta\delta E}^{00}$ at $\sim 30\sigma$ and the two quadrupoles $B_{\delta \delta E}^{11}$ and $B_{\delta \delta E}^{20}$ at $\sim 25\sigma$ and $\sim 15\sigma$, respectively. We also report similar detection significances for the lowest order multipoles of $B_{\delta EE}$ and $B_{EEE}$, although these are largely driven by stochastic contributions. We show that the first and second order EFT parameters are consistent with those obtained from fitting the IA power spectrum analysis at next-to-leading order, without requiring any priors to break degeneracies for the quadratic bias parameters. Moreover, the inclusion of higher multipole moments of $B_{\delta\delta E}$ greatly reduces the errors on second order bias parameters, by factors of 5 or more. The IA bispectrum thus provides an effective means of determining higher order shape bias parameters, thereby characterizing the scale dependence of the IA signal. We also detect parity-odd bispectra such as $B_{\delta \delta B}$ and $B_{\delta EB}$ at $\sim 10 \sigma$ significance or more for $k<0.15\,h$/Mpc and they are fully consistent with the parity-even sector. Furthermore, we check that the Gaussian covariance approximation works reasonably well on the scales we consider here. These results lay the groundwork for using the bispectrum of IA in cosmological analyses.

Context. Deriving the mass and large-scale asymmetries of radioactive isotopes offers valuable insights into the complex phases of a supernova explosion. Important examples are $^{56}$Ni, with its decay products $^{56}$Co and $^{56}$Fe, and $^{44}$Ti, which are studied through their X-rays emission lines and provide a powerful diagnostic tool to probe the explosive nucleosynthesis processes in the inner layers of the exploding star. Aims. In this framework, SN 1987A provides a privileged laboratory being the youngest supernova remnant from which the mass of Ti has been estimated. However, some tension exists in determining the initial mass of $^{44}$Ti. Previous analysis, relying on \textit{NuSTAR} and \textit{INTEGRAL} data, report $M_{44} = (1.5 \pm 0.3) \times 10^{-4}$ $M_\odot$ and $M_{44}=(3.1 \pm 0.8) \times 10^{-4} M_\odot$, respectively. In this paper we estimate the initial mass of $^{44}$Ti via its decay product, the $^{44}$Sc emission line at 4.09 keV, using \textit{Chandra} observations. Methods. We perform multi-epoch spectral analysis focusing on the inner part of the remnant, to minimize the contamination from the X-ray emission stemming from the shocked plasma. As a result, we provide the detection of $^{44}$Sc emission line in the central part of SN 1987A with a $\sim$99.7\% (3 $\sigma$) significance. Results. The simultaneous fit of the spectra extracted from observations between 2016 and 2021 provides a line flux of $6.8^{+2.2}_{-2.3}\times 10^{-7}$ photons s$^{-1}$ cm$^{-2}$ corresponding to a $^{44}$Ti mass $M_{44}=(1.6\pm0.5) \times 10^{-4} M_\odot$ (errors at the $90\%$ confidence level). The results obtained with our spectral analysis seem to align with those derived with NuSTAR.

Astrophysical neutrinos provide crucial insights into their sources and play a key role in multi-messenger astronomy. The neutrino flavor composition at Earth allows us to probe the mechanisms of neutrino production and cosmic ray acceleration, as well as the properties of the environments in which they originate. Understanding the flavor composition also offers a unique opportunity to test new physics in the neutrino sector. The IceCube Neutrino Observatory consists of 1 km$^3$ of ice instrumented with photomultipliers that detect neutrinos through Cherenkov radiation from their interaction products. Different neutrino interactions result in distinct event topologies, such as tracks, cascades, and double cascade events, which allow for the identification of the interacting neutrino type and measurement of the flavor composition of the astrophysical neutrino flux. In this contribution we present the results of the measurement of the flavor ratio of the High-Energy Starting Event Selection based on 12 years of data, a high-purity sample of neutrino interactions that occur inside the detector. In addition, we discuss various methods that could further improve the analysis in the future.

The dispersion measure (DM) analysis of extragalactic fast radio bursts (FRBs) has accounted for all cosmic baryons but remains limited by systematic uncertainties from various parameter degeneracies. We show that the prominent degeneracy between the baryon density ($\Omega_{\rm b}$) and the Hubble constant ($H_0$) can be disentangled through an independent $H_0$ value uniquely inferred from the absolute luminosity distance measurements of gravitational-wave (GW) standard sirens. By combining $104$ localized FRBs with $47$ GW events, we obtain a robust, late-Universe measurement of \ob $=0.0488\pm0.0064$ ($1\sigma$), in concordance with early-Universe observations of CMB + BBN. Notably, incorporating GW data helps not only avoid biases induced by the $H_0$ tension, but also mitigate those from the parameters of diffuse baryon fraction and DM distribution models. Although the current precision ($\sim 13\%$) is limited by sample size, the growing detections of both FRBs and GWs will make their synergy a powerful probe of low-redshift cosmology.

Studying the composition and origin of the inner region of our Galaxy -- the "Galactic bulge" -- is crucial for understanding the formation and evolution of the Milky Way and other galaxies. We present new observational constraints based on a sample of around 18,000 stars in the inner Galaxy, combining Gaia DR3 RVS and APOGEE DR17 spectroscopy. Gaia-RVS complements APOGEE by improving sampling of the metallicity, [Fe/H] in the -2.0 to -0.5 dex range. This work marks the first application of Gaia-RVS spectroscopy to the bulge region, enabled by a novel machine learning approach (hybrid-CNN) that derives stellar parameters from intermediate-resolution spectra with precision comparable to APOGEE's infrared data. We performed full orbit integrations using a barred Galactic potential and applied orbital frequency analysis to disentangle the stellar populations in the inner Milky Way. For the first time, traced by the field stars, we are able to robustly identify the long-sought pressure supported bulge. We show this stellar population to be chemically and kinematically distinct from the other main components co-existing in the same region. The spheroidal bulge has a metallicity distribution function (MDF) peak at around -0.70 dex extending to solar value, is dominated by a high-[alpha/Fe] population with almost no dependency on metallicity, consistent with very rapid early formation, predating the thick disc and the bar. We find evidence that the bar has influenced the dynamics of the spheroidal bulge, introducing a mild triaxiality and radial extension. We identify a group of stars on X4 orbits, likely native to the early spheroid, as this population mimics the chemistry of the spheroidal bulge, with a minor contamination from the more metal-poor ([Fe/H] < -1.0) halo. We find the inner-thick disc to be kinematically hotter (mean Vphi ~125 km/s) than the local-thick disc. ...

While cosmic voids are now recognized as a valuable cosmological probe, identifying them in a galaxy catalog is challenging for multiple reasons: observational effects such as holes in the mask or magnitude selection hinder the detection process; galaxies are biased tracers of the underlying dark matter distribution; and it is non-trivial to estimate the detection significance and parameter uncertainties for individual voids. Our goal is to extract a catalog of voids from constrained simulations of the large-scale structure that are consistent with the observed galaxy positions, effectively representing statistically independent realizations of the probability distribution of the cosmic web. This allows us to carry out a full Bayesian analysis of the structures emerging in the Universe. We use 50 posterior realizations of the large-scale structure in the Manticore-Local suite, obtained from the 2M++ galaxies. Running the VIDE void finder on each realization, we extract 50 independent void catalogs. We perform a posterior clustering analysis to identify high-significance voids at the 5$\sigma$ level, and we assess the probability distribution of their properties. We produce a catalog of 100 voids with high statistical significance, available at this https URL, including the probability distributions of the centers and radii of the voids. We characterize the morphology of these regions, effectively producing a template for density environments that can be used in astrophysical applications such as galaxy evolution studies. While providing the community with a detailed catalog of voids in the nearby Universe, this work also constitutes an approach to identifying cosmic voids from galaxy surveys that allows us to account rigorously for the observational systematics intrinsic to direct detection, and provide a Bayesian characterization of their properties.

Accurate reconstruction of the electric field produced by extensive air showers is essential for the radio-detection technique, as the key parameters of interest of the primary particles that generated the showers are the amplitude, polarization, frequency spectrum, and energy fluence carried by the electric field at each receiving radio antenna. Conventional electric-field reconstruction methods primarily focus on antennas with two horizontal polarizations. In this work, we introduce an analytic $\chi^2$ minimization method that is applicable to both two and three polarizations. This solution has been verified for simple and realistic antenna responses, with a particular focus on inclined air showers. Our method achieves standard deviations better than 4\% and 6\% for the estimation of the Hilbert peak envelope amplitude of the electric field and the energy fluence, respectively, with an antenna-response-dependent bias. Additionally, we have studied the dependence of the method with arrival direction showing that it has a good performance in the zenith range from 63$^\circ$ up to 80$^\circ$. This work also demonstrates that incorporating vertically polarized antennas enhances the precision of the reconstruction, leading to a more accurate and reliable electric-field estimation for inclined air showers.

The origin of the UHECR continues to puzzle, however, an excess of detection in the direction of the radio galaxy Centaurus A (Cen A) raises the possibility of this object being the first UHECR source identifiable. Cen A is known to be currently active, and also exhibits known past episodes of high activity. In this work, we investigate whether the known activity episodes in Cen A may be related to the excess events in the \textit{Centaurus region}. Analysing the energy of the events and the overall mass composition of UHECR, we report that an activity in the last $\sim30$ Myr is necessary to explain the excess of events. This period perfectly fits with the timescale where the transition regions and the Giant Lobes must be energized, as revealed by radio and $\gamma$ ray observations.

In this work, we derive upper limits on the physical energy-density fraction today of cosmological gravitational waves, denoted by $\Omega_{\rm{gw}}h^{2}$, by combining \emph{Planck} \& ACT \& SPT CMB data with DESI BAO data. In the standard cosmological model, we establish 95\% CL upper limits of $\Omega_{\rm{gw}}h^{2} < 1.0 \times 10^{-6}$ for adiabatic initial conditions and $\Omega_{\rm{gw}}h^{2} < 2.7 \times 10^{-7}$ for homogeneous initial conditions. In light of dynamical dark energy, we get $\Omega_{\rm{gw}}h^{2} < 7.2 \times 10^{-7}$ (adiabatic) and $\Omega_{\rm{gw}}h^{2} < 2.4 \times 10^{-7}$ (homogeneous). We also project the sensitivity achievable with LiteBIRD \& CMB Stage-IV measurements of CMB and CSST observations of BAO, forecasting 68\% CL uncertainties of $\sigma = 2.4 \times 10^{-7}$ (adiabatic) and $\sigma = 0.9 \times 10^{-7}$ (homogeneous) for ${\Omega_{\rm{gw}}h^{2}}$. The constraints provide critical benchmarks for exploring the cosmological origins of gravitational waves within the frequency band $f \gtrsim 10^{-15}$\,Hz and potentially enable joint analysis with direct gravitational-wave detection sensitive to this regime.

IceCube has detected a diffuse flux of high-energy neutrinos, whose origin still remains uncertain. Accreting supermassive black holes (SMBHs) have been proposed as plausible sources of neutrinos. Candidate sources include AT2019dsg, which is likely a stellar tidal disruption event (TDE), and AT2019dfr, an AGN flare. Both present delayed emission in the IR band with respect to the optical signal. This emission can be interpreted as the reprocessing of X-rays to optical light of the flare by dust located in a torus around the SMBH. An additional study using an optically detected sample of 63 accretion flares revealed another candidate as a potential high-energy neutrino counterpart: AT2019aalc, which is also accompanied by a dust echo. However, follow-up stacking analysis of the 63 nuclear flares using the full IceCube data sample did not show any significant excess over background. Motivated by these three suggested neutrino-TDE correlations, we analyze a more extensive catalog of IR flares, 823 dust-echo-like flares identified using WISE satellite data, against the IceCube 10-year sample of track events from the Northern Sky. Our analysis aims to perform sensitivity studies and assess the potential detectability of neutrino emission from these types of accretion flares. In addition, we carry out a correlation study of the 823 dust echo-like flares against a revised catalog of IceCube high-purity astrophysical alerts, and reevaluate the previous study of 63 nuclear flares against the same revised alerts sample.

In 2016, the IceCube Neutrino Observatory launched its realtime program. When a neutrino candidate of likely astrophysical origin is detected, a public alert is issued, typically within one minute. These alerts allow the astrophysical community to follow up on the region of the sky where the neutrino likely originated. Initially, the system issued around six track-signature alerts per year, with a highlight being IceCube-170922A, which was later associated with the flaring blazar TXS 0506+056. Since 2019, IceCube has expanded the selection criteria for neutrino candidates, increasing the track alert rate to around 30 per year with additional alerts for cascade signatures of probable astrophysical origin implemented in 2020. This work describes several improvements in the reconstruction of track alerts, which were introduced into the realtime stream in 2024, and details the two algorithms that are alternately used for reconstruction, depending on the reconstructed muon energy. The improvements result in more precise directional localizations, with a factor of 5 (4) reduction in the 50% (90%) contour area. Systematic errors affecting the reconstructions, such as the ice model or the geometry of the detector, have also been investigated to ensure statistical coverage. In addition to its application in the realtime stream, the improvements described here are also being applied to historical alerts with an updated catalog of track alerts forthcoming.

Observing the stars in our night sky tells us that giant, supergiant and hypergiant stars hold an unique importance in the understanding of stellar populations. Theoretical stellar models predict a rich tapestry of evolved stars. These evolved stars result in supernova explosions for massive stars and the shedding of the outer layers of low- and intermediate-mass stars to form white dwarf stars. Stars on these pathways provide elemental and kinematic feedback that shapes their host galaxies and synthesize the elements that we observe in the Universe today. Observational surveys in the Milky Way and the Local Group of galaxies shape our physical understanding of stellar evolution. Hypergiant stars represent stellar evolution at the extremes, displaying evidence of intense stellar winds combined with luminosities and radii that are unrivalled. Supergiant stars are used as elemental beacons in distant galaxies, and supergiants that cross the Cepheid variable instability strip can be used as rulers to measure the scale of our Universe. The overwhelming majority of stars in the Milky Way will become giant stars. Giants are hundreds of times larger and thousands of times more luminous than the Sun, which makes them excellent tracers of stellar structure. In this chapter, I describe the observational and theoretical manifestations of stellar evolution at various stages and masses. I highlight evolutionary connections between evolved products and comment on their eventual endpoints.

This study employs well-localized fast radio bursts (FRBs) to constrain cosmological parameters using two model-independent methods: the reconstruction of the Hubble parameter $H(z)$ with an artificial neural network and cosmography. By integrating FRB data with type Ia supernovae (SNe), BAO from DESI DR2, and cosmic chronometers (CC), we derive constraints on the Hubble constant ($H_0$), deceleration parameter ($q_0$), and jerk parameter ($j_0$). For the reconstruction method, our MCMC analysis with FRB-only provides $H_0 = 69.9 \pm 5.8 \, \text{km} \, \text{s}^{-1} \, \text{Mpc}^{-1}$, corresponding to a precision of $\sim 8\%$. A joint analysis with FRB+SNe+DESI gives $H_0 = 68.85_{-0.48}^{+0.47} \, \text{km} \, \text{s}^{-1} \, \text{Mpc}^{-1}$, reaching a precision below 1\%. The cosmographic approach with FRBs alone provides $H_0 = 65.83_{-4.87}^{+3.77} \, \text{km} \, \text{s}^{-1} \, \text{Mpc}^{-1}$, $q_0 = -0.45_{-0.31}^{+0.26}$, and $j_0 = 1.17_{-0.58}^{+0.70}$ with a precision for the Hubble constant of $\sim 6\%$. Besides, the DESI DR2 dataset yields $H_0 = 65.20_{-1.28}^{+1.29} \, \text{km} \, \text{s}^{-1} \, \text{Mpc}^{-1}$, $q_0 = -0.29\pm 0.07$, and $j_0 = 0.58_{-0.04}^{+0.03}$, indicating an accuracy of $\sim 2\%$ for the Hubble constant. Combining the FRB, SNe, DESI, and CC datasets provide tighter constraints, e.g., $H_0 = 67.88_{-0.53}^{+0.52} \, \text{km} \, \text{s}^{-1} \, \text{Mpc}^{-1}$, $q_0 = -0.42_{-0.03}^{+0.02}$, and $j_0 = 0.56 \pm 0.02$. Notably, these findings provide a statistically significant indication of deviation from the $\Lambda$CDM prediction of $j_0 = 1$, suggesting possible departures from standard cosmology at a $1\sigma$ confidence level. They also confirm a model-independent accelerated expansion ($q_0 < 0$), demonstrating the potential of FRBs as valuable cosmological probes.

While high-energy astrophysical neutrinos are well-established, their flavor composition remains relatively unconstrained. In IceCube, long muon tracks from $\nu_\mu$ charged-current interactions are easily identified but the detector geometry does not allow sufficient resolution to distinguish cascade-type events. The Neutron Echo - a delayed light signal in the detector from neutron capture and de-excitation - can probe the shower's hadronic content and thus the underlying interaction. A significant background arises from late PMT afterpulses, which are temporally coincident with the physics signal. The traditional IceCube data acquisition system has a limited readout window with significant deadtime between triggers, which is insufficient to capture the late pulses. A recently developed deadtime-free readout mode, with an extended window, enables their detection. An observed excess in the delayed time spectrum over the background would be compatible with the Neutron Echo hypothesis. In this contribution, we summarize the physics scope of delayed signals, discuss the timing spectrum of signal and PMT background, highlight the capabilities of the new system for recording late pulses, and emphasize the potential of IceCube for particle identification through delayed signals.

Properties of interstellar environment near the solar system have been probed by missions like IBEX, Voyager over the last two decades. Although it has been well recognized that properties of cosmic rays up to the PeV energy can be affected by the local interstellar environment, detailed modeling has not been done. We show that a three component model for the cosmic ray proton and helium spectra from GV to several PV can naturally explain the energy dependence of the dipole anisotropy of cosmic ray fluxes by considering effects of the local interstellar environment on cosmic ray transport, addressing the so-called cosmic ray anisotropy problem. In particular, it is shown that the dipole amplitude and position angle below $\sim 100$ TeV are very sensitive to the velocity of the heliosphere in the local interstellar cloud and the motion of the local interstellar cloud in the local standard of rest. Better measurement of cosmic ray flux anisotropy by experiments like LHAASO and properties of the local interstellar environment by future missions like IMAP will be able to test this model.

The IceCube Upgrade, currently under construction at the geographic South Pole, is the next development stage of the IceCube detector. It will consist of seven new columns of novel optical sensors and advanced calibration devices densely deployed at the centre of the existing array. The sensors are frozen into the ice in boreholes created by hot water drilling. The refreezing forms hole ice around modules with optical properties that differ from the surrounding glacial ice. A key objective of the IceCube Upgrade is to enhance our understanding of the optical properties of both bulk ice and refrozen hole ice. Precise ice modelling is crucial for the directional reconstruction of TeV-PeV neutrinos, as resolutions at such high energies can be strongly impacted by uncertainties in ice properties and optical sensor response. An improved directional reconstruction performance will translate to a boost in neutrino source sensitivities using IceCube data collected over the last 12 years. In this contribution, we present the expected improvements in reconstruction performance resulting from advances in hole-ice modelling and the resulting impact on IceCube's sensitivity to astrophysical neutrino sources across three distinct event samples.

We present preliminary results on an updated full-sky analysis of the cosmic-ray arrival direction distribution with data collected by the High-Altitude Water Cherenkov (HAWC) Observatory and IceCube Neutrino Observatory with complementary field of views covering a large fraction of the sky. This study extends the energy range to higher energies. The HAWC Observatory, located at 19$^{\circ}$N has analyzed 8 years of cosmic-ray data over an energy range between 3.0 TeV and 1.0 PeV and confirms an energy-dependent anisotropy in the arrival direction distribution of cosmic rays seen by other experiments. Combined with recently published results from IceCube with 12 years of data, the combined sky maps with 93\% coverage of the sky -- between 70$^{\circ}$N and 90$^{\circ}$S -- and the corresponding angular power spectra largely eliminate biases that result from partial sky coverage.

The LHAASO collaboration has recently released the spectrum and the angular distribution of the gamma-ray Galactic diffuse emission from 1 TeV to 1 PeV measured with the Kilometer-2 Array (KM2A) and the Water Cherenkov Detector Array (WCDA). We show that those data are in remarkably good agreement with a set of pre-existing models that assume the emission to be produced by the Galactic population of cosmic rays if its spectral shape traces that measured by CALET, DAMPE as well as KASCADE at higher energies. No extra-components besides the CR sea is needed to explain LHAASO results. Spatial dependent CR transport models, although not required to reproduce LHAASO results, are in better agreement with them respect to conventional ones and needed to consistently reproduce Fermi-LAT and neutrino data.

Supernova remnants (SNRs) have long been suspected to be the primary sources of Galactic cosmic rays. Over the past decades, great strides have been made in the modelling of particle acceleration, magnetic field amplification, and escape from SNRs. Yet while many SNRs have been observed in nonthermal emission in radio, X-rays, and gamma rays, there is no evidence for any individual object contributing to the locally observed flux. Here, we propose a particular spectral signature from individual remnants that is due to the energy-dependent escape from SNRs. For young and nearby sources, we predict fluxes enhanced by tens of percent in narrow rigidity intervals; given the percent-level flux uncertainties of contemporary cosmic-ray data, such features should be readily detectable. We model the spatial and temporal distribution of sources and the resulting distribution of fluxes with a Monte Carlo approach. The decision tree that we have trained on simulated data is able to discriminate with very high significance between the null hypothesis of a smooth distribution of sources and the scenario with a stochastic distribution of individual sources. We suggest that this cosmic-ray energy-dependent injection time (CREDIT) scenario be considered in experimental searches to identify individual SNRs as cosmic-ray sources.

The dispersion measures (DMs) of Fast Radio Bursts (FRBs) arise predominantly from free electrons in the large-scale structure of the Universe. The increasing number of FRB observations have started to empirically constrain the distribution of cosmic baryons, making it crucial to accurately forward model their propagation within cosmological simulations. In this work, we present a method for measuring FRB DMs in IllustrisTNG that continuously traces rays through the simulation while reconstructing all traversed line segments within the underlying Voronoi mesh. Leveraging this technique, we create over $20$ publicly available DM catalogs, including a full-sky DM map observed from a Milky Way-like environment. Our method addresses a problem in previous TNG-based studies, in which a sparse snapshot sampling in the path integral leads to a misestimation of the standard deviation and higher moments of the DM distribution $p(\rm{DM}|z)$ by over $50\%$. We show that our results are consistent with the most recent observational data from the DSA-110, ASKAP, and CHIME. We offer a functional form for $p(\rm{DM}| z)$ that provides very good fits across all redshifts from $0$ to $5.5$, showing that the previously proposed log-normal distribution is not well matched to the data. We find that using simulation box sizes smaller than $35\,h^{-1}\,$Mpc or resolutions with baryonic masses of less than $5\times10^{8}\,h^{-1}\,M_{\odot}$ can distort the derived DM signal by more than $8\%$. Our findings provide new insights into how the cosmic web shapes FRB signals, and highlight the importance of accurate methodology when comparing cosmological simulations against observations.

Mapping the local and distant Universe is key to our understanding of it. For decades, the Sloan Digital Sky Survey (SDSS) has made a concerted effort to map millions of celestial objects to constrain the physical processes that govern our Universe. The most recent and fifth generation of SDSS (SDSS-V) is organized into three scientific ``mappers". Milky Way Mapper (MWM) that aims to chart the various components of the Milky Way and constrain its formation and assembly, Black Hole Mapper (BHM), which focuses on understanding supermassive black holes in distant galaxies across the Universe, and Local Volume Mapper (LVM), which uses integral field spectroscopy to map the ionized interstellar medium in the local group. This paper describes and outlines the scope and content for the nineteenth data release (DR19) of SDSS and the most substantial to date in SDSS-V. DR19 is the first to contain data from all three mappers. Additionally, we also describe nine value added catalogs (VACs) that enhance the science that can be conducted with the SDSS-V data. Finally, we discuss how to access SDSS DR19 and provide illustrative examples and tutorials.

The widely anticipated outburst of recurrent nova T Coronae Borealis (T CrB), which is near the end of its 80-year cycle, provides an excellent opportunity to search for neutrinos from novae. Novae are an energetic class of transients, which have been studied for hundreds of years. Because many of them are located nearby, novae provide an excellent astrophysical laboratory to study shock-powered emission in our own backyard. Several recent novae have previously been detected in GeV gamma rays, and the 2021 outburst of RS Ophiuchi was detected up to TeV energies, with evidence for a hadronic origin of the observed emission. Previous searches for GeV-TeV neutrinos from novae, predicted to occur alongside their gamma-ray emission, have been performed using data from the IceCube Neutrino Observatory. However, no significant neutrino signals from novae have yet been observed. We present plans for follow-up of T CrB in real time with IceCube, using datasets spanning GeV to PeV neutrino energies. Due to its closer distance and higher optical flux, which has been well measured in two historical eruptions, the expected neutrino signal from T CrB is several times stronger than that from RS Ophiuchi. Furthermore, T CrB is located in the Northern sky at a declination where IceCube's sensitivity is an additional factor of a few better than at the location of RS Ophiuchi, which is beneficial to this search.

Core-collapse supernovae are of particular interest in multi-messenger astronomy due to their potential to accelerate cosmic rays and produce high-energy neutrinos. One such supernova is the recent SN2023ixf located in M101 (the Pinwheel Galaxy). It is the closest (6.4 Mpc) and brightest (B band magnitude 10.8) core-collapse supernova in nearly a decade. This supernova likely had a progenitor surrounded by dense circumstellar material which, during the supernova, may have produced neutrinos when ejecta collided with the material. I will present results of a follow-up of this supernova using data collected from the IceCube Neutrino Observatory located at the South Pole. We obtain results consistent with background expectations with time-integrated energy flux ($E^2 dN/dE$) upper limits of 0.35 GeV/cm$^2$ for a 32-day time window and 0.44 GeV/cm$^2$ for a 4-day time window, both at 90% confidence level for an $E^{-2}$ power law. These correspond to values of $2.7 \times 10^{48}$ erg for the 32-day time window and $3.5 \times 10^{48}$ erg for the 4-day time window at the supernova.

Gamma-ray-bright active galactic nuclei (AGN) have been one of the most promising source classes of high-energy astrophysical neutrinos detected by IceCube. The first evidence of an IceCube point source was a blazar detected by the Fermi Large Area Telescope (LAT), TXS ~0506+056. Previous analyses have ruled out GeV-bright blazars as the predominant contributor to the high-energy astrophysical neutrino flux under simple correlation assumptions about the relationship between the fluxes of gamma rays and neutrinos. We present results from a more general and more sensitive search for correlation between neutrinos and GeV-selected AGN using improvements in the IceCube statistical methods and 13 years of data. We detect no correlation and set stringent constraints on neutrino emission by populations of GeV-detected AGN. These include constraints on the neutrino emission from subclasses of GeV-detected AGN, including BL Lacs, flat-spectrum radio quasars (FSRQ) and non-blazar AGN, using stacking analyses testing a variety of hypothesized relationships between neutrino and gamma-ray flux. We also present results from an analysis that is sensitive to a wider range of relationships between the gamma-ray and neutrino signal.

KM3NeT has recently reported an event where a muon of energy $120^{+110}_{-60}$ PeV was observed at its ARCA detector, which can stem from a very high-energy neutrino interaction in the vicinity of the detector. Besides revolutionizing our understanding of high-energy neutrino sources, this event can serve as a valuable probe for studying Beyond the Standard Model (BSM) interactions of neutrinos. In this work, we study the dark matter (DM)-neutrino interaction by assuming the neutrino for the event KM3-230213A is originated from a blazar. The flux of such neutrinos, traveling through DM distributed across astrophysical and cosmological scales, can get attenuated due to DM interactions. The detection of such event by KM3NeT allows us to place constraints on the interaction cross section at highest-ever neutrino energy. We derive both conservative constraints-neglecting flux attenuation from the host halo-and optimistic ones by including host halo contributions. Our results show that the energy-independent constraints are weaker than previous bounds. For energy-dependent case, the extreme energy of the event allows us to set some of the strongest limits on scattering cross sections. In future, more such neutrino events with well-understood origin will be essential in constraining or potentially discovering DM-neutrino interactions.

We present a semi analytic forecast for the detection of intermediate mass black hole (IMBH) binaries with the space based detectors LISA (millihertz band) and AMIGO (deci hertz band). A redshift dependent population model is built from extrapolated black hole mass functions and realistic pairing and merger time scales. Folding this population through the instrument sensitivities yields event rates of about 1e4 mergers for LISA in four observing years and between 1e2 and 1e3 mergers for AMIGO in three years, with complementary optimal detection bands in total mass and redshift. We then incorporate the minimal matter coupled modified gravity model f(R,T) = R + 2 lambda T. In this scenario the strain amplitude is rescaled by (1 - lambda / 8 pi)^(5/6) while the waveform phasing is unchanged at leading order. For Solar System compatible couplings abs(lambda) <= 2 x 10^-3, the peak detection redshift shifts by less than one percent and the total number of events changes by less than five percent. Combined statistics from LISA and AMIGO therefore provide a consistency check on f(R,T) gravity without spoiling standard IMBH forecasts.

We will present the way to derive a de Sitter inflationary solution within the so-called Starobinsky-Bel-Robinson gravity. Then, we will show by using the dynamical system method whether the obtained solution is stable or not. According to the stability of the de Sitter inflationary solution, we could judge which phase of our universe, among the two early and late-time phases, is more appropriate for this solution.

Multiple axions may emerge in the low-energy effective theory of Nature. Generically, the potentials describing these axion fields are non-diagonal, leading to mass mixing between axion states which can be temperature-dependent due to QCD instanton effects. As the temperature of the Universe drops, level crossing can occur, causing resonant conversion between axion states. In this work, we present an analytic study of the cosmological evolution of multi-axion systems including adiabatic and non-adiabatic resonant conversion from one axion state into another during the misalignment process. We show how the Landau-Zener formalism accurately captures the non-adiabatic resonant conversion, permitting an analytic description of the relic abundances of each axion field for nearly any arbitrary two-state axion mass matrix. As an application, we study the mixing of a QCD axion with an axion-like-particle for specific potentials to identify the predictions for haloscope experiments. We conclude that the detection of an axion off the expected QCD mass-coupling line predicts other haloscope targets if it mixes with the QCD axion.

We demonstrate Bayesian analyses of the complete gravitational-wave spectrum of binary neutron star mergers events with the next-generation detector Einstein Telescope. Our mock analyses are performed for 20 different signals using the TEOBResumSPA_NRPMw waveform that models gravitational-waves from the inspiral to the postmerger phase. They are employed to validate a pipeline for neutron star's extreme matter constraints with prospective detections and under minimal hypotheses on the equation of state. The proposed analysis stack delivers inferences for the mass-radius curve, the mass dependence of the quadrupolar tidal polarizability parameter, the neutron star's maximum density, the maximum mass and the relative radius, and the pressure-density relation itself. We show that a single event at a signal-to-noise ratio close to the minimum threshold for postmerger detection is sufficient to tightly constrain all the above relations as well as quantities like the maximum mass (maximum density) to precision of ${\sim}6$% (${\sim}10$%) at 90% credibility level. We also revisit inferences of prompt black hole formation with full spectrum signals and find that the latter can be robustly identified, even when the postmerger is not detectable due to a low signal-to-noise ratio. New results on the impact of the initial signal frequency and of the detector configuration (triangular vs. two-L) on the source's parameters estimation are also reported. An improvement of approximately one order of magnitude in the precision of the chirp mass and mass ratio can be achieved by lowering the initial frequency from 20 Hz to 2 Hz. The two-L configuration shows instead significant improvements on the inference of the source declination, due to geographical separation of the two detectors.

LISA data analysis represents one of the most challenging tasks ahead for the future of gravitational-wave (GW) astronomy. Characterizing the instrument's noise properties while fitting for all the other detectable sources is a key requirement of any robust inference pipeline. Noise estimation will also play a crucial role in searches and parameter estimation of cosmological and astrophysical stochastic signals. Previous studies have tackled this topic by assuming perfect knowledge of the spectral shape of the instrumental noise and of different possible types of GW Stochastic Backgrounds (SGWBs), usually resorting to parametrized templates. Recently, various works that employ template-agnostic methods have been presented. In this work, we take an additional step further, introducing flexible spectral shapes in both the instrumental noise and the stochastic signals. We account for the lack of knowledge of the exact shape of the individual contributions to the overall power spectral density by using splines to represent arbitrary perturbations of the noise and signal spectral densities. We implement a data-driven Reversible Jump MCMC algorithm to fit different components simultaneously and to infer the level of flexibility required under different scenarios. We test this approach on simulated LISA data produced under different assumptions. We investigate the impact of this increased flexibility on the reconstruction of both the injected signal and the noise level, and we discuss the prospects for claiming a successful SGWB detection.

We explore the phenomenology of the ``Dark Walker'' -- an $\text{SU}(3)$ theory with eight flavors of fundamental fermions in the dark sector. During inflation, its walking dynamics generate primordial non-Gaussianities through the exchange of unparticles, while accounting for the current dark matter relic abundance if the resulting massive pseudoscalars behave as strongly interacting massive particles (SIMPs). This provides a simple yet concrete example of a strongly coupled sector during inflation.

This work presents a geometric formulation for transforming nonconservative mechanical Hamiltonian systems and introduces a new method for regularizing and linearizing central force dynamics -- in particular, Kepler and Manev dynamics -- through a projective transformation. The transformation is formulated as a configuration space diffeomorphism (rather than a submersion) that is lifted to a cotangent bundle (phase space) symplectomorphism and used to pullback the original mechanical Hamiltonian system, Riemannian kinetic energy metric, and other key geometric objects. Full linearization of both Kepler and Manev dynamics (in any finite dimension) is achieved by a subsequent conformal scaling of the projectively-transformed Hamiltonian vector field. Two such conformal scalings are given, both achieving linearization. Arbitrary conservative and nonconservative perturbations are included, with closed-form solutions readily obtained in the unperturbed Kepler or Manev cases.

We show that a class of NJL-like models fails to reproduce the expected conformal limit of the speed of sound, making them unsuitable for analyzing the equation of state of dense matter. We then demonstrate how this issue can be resolved within a simple dynamical quark model.

A relativistic mean-field model with a crossing term of isovector-scalar and isoscalar mesons, the Cracow crossing terms (CCT) model, has previously been shown to be capable of explaining the light and compact object HESS J1731-347. Here, the model is supplemented with a correct treatment of the phase transition. Applying the thermodynamically consistent Gibbs conditions shows that a mixed-phase system occurs in the core of the neutron star over a wide range of densities, extending up to the neutron star crust, leading to an intriguing stellar structure with a thick and massive crust.

Early warning of gravitational waves (GWs) is essential for multi-messenger observations of binary neutron star and black hole-neutron star merger events. In this study, we investigate early warning prospects from eccentric compact binaries, whose mergers are expected to comprise a significant fraction of detected GW events in the future. Eccentric binaries exhibit oscillatory frequency evolution, causing GW frequencies to recur multiple times through their coalescence. Consequently, generating eccentric waveform templates for early warning requires specification of initial conditions. While the standard approach involves initiating waveform generation when the orbit-averaged frequency enters the detector band, we compare this with an alternative approach that uses the periastron frequency as the starting point. Our analysis shows that initializing at the periastron frequency yields an improved signal-to-noise ratio and sky localization. Additionally, including subdominant modes alongside the dominant $(2,2)$ mode leads to further improvements in sky localization. We explore the parameter space of primary mass $m_1 \in [1.4, 15] \, M_\odot$, spin $\chi_1 \in [0, 0.8]$, and eccentricity $e \leq 0.4$ across three detector configurations: O5, Voyager, and 3G. We find that in the O5 (Voyager) configuration, including eccentricity and subdominant modes, the sky localization area can be reduced by $2-80\% (2-85\%)$ at 1000 sq. deg. with increasing eccentricity from $e_5 = 0.1$ to $e_5 = 0.4$, yielding up to $41$ seconds (1 minute) of extra early warning time. For NSBH systems, subdominant modes contribute up to $70$ $(94)\%$ reduction for O5 (Voyager) scenario. In the 3G detector scenario, the sky area reduction due to eccentricity reaches $80\%$ (from $e_{2.5} = 0.1$ to $e_{2.5} = 0.4$) at 100 sq. deg., and subdominant modes enhance the reduction up to $98\%$ for NSBH systems.

We propose two novel solutions to the domain wall problem of the QCD axion by introducing a massless or light axion that also couples to gluons. The first solution applies when the new axion forms strings after inflation. Due to its mixing with the QCD axion, domain walls of the QCD axion are bounded by these strings and confined into cosmologically safe string bundles. This scenario predicts the existence of such string bundles, which may survive until today and leave observable signatures, such as gravitational waves, cosmic birefringence, and CMB anisotropies. The simultaneous detection of the QCD axion and any of these cosmological signatures would serve as a smoking-gun signal. The second solution assumes a homogeneous initial condition for the new axion. If it is sufficiently light, its potential temporarily induces a bias in the QCD axion potential before the onset of oscillations, rendering the domain walls unstable. In both scenarios, the Peccei-Quinn mechanism remains effective, and the strong CP problem is not reintroduced. We identify the viable parameter regions and discuss the resulting dark matter abundance.