Papers for Friday, Jul 04 2025

Magnetic fields play a crucial role in planetary evolution and habitability. While the intrinsic magnetic fields of solar system planets are relatively well understood, the magnetic properties of exoplanets remain largely unconstrained, despite their potential ubiquity. Detecting exoplanetary magnetic fields is essential to advancing our understanding of planetary habitability beyond the solar system. This paper focuses on two promising spectropolarimetric techniques for detecting magnetic fields in hot exoplanets: direct detection through polarization signatures in the He I 1083 nm triplet and indirect detection via star-planet magnetic interactions manifesting as stellar hot spots. The direct method is particularly suited to close-in gas giants, leveraging the Hanle and Zeeman effects to detect low-amplitude magnetic fields. The indirect method can apply to both giant and low-mass planets by identifying magnetic connectivity-induced features in the stellar atmosphere. Although the interpretation of current detections remain tentative, upcoming high-resolution spectropolarimetric capabilities in the UV and near-infrared, particularly with future missions like HWO, promise to enable definitive measurements of exoplanetary magnetic fields. These advancements will open new avenues for probing the magnetic environments of exoplanets and their implications for atmospheric retention and habitability.

Light pollution, a typically underrecognized environmental issue, has gained attention in recent years. While controlling light pollution requires sustained efforts, Earth Hour offers a unique natural experimental setting to assess temporary lights-out measures. Using photometric and spectroscopic sensors, we observed up to 50% night sky darkening during Earth Hour from 2011 to 2024 in Hong Kong, primarily as a result of a small but critical number of lights-out instances in central business districts, as evidenced by crowd-sourced photography records. Weekend lighting pattern in the city during Earth Hour remained unaffected. The emission reductions mostly occurred in the 445-500, 500-540, and 615-650 nm spectral ranges-corresponding to peak emissions from LED billboard screens- and in the 585-595 nm range, associated with metal halide floodlights used for facades and billboards. Our study identifies these sources as major contributors to urban light pollution. By combining multimodal observations, we offer a comprehensive assessment of light pollution sources and the potential benefits of sustainable lighting practices in urban environments. This research highlights the importance of targeted light pollution mitigation efforts and provides critical insights for policymakers to enhance urban sustainability and human well-being.

We present deep $K_s$-band imaging of the extreme T subdwarf CWISE J221706.28$-$145437.6. Using the new photometry, we construct its spectral energy distribution and find this object exhibits exceptionally strong collision-induced absorption in the $H$ and $K$ band. The comparison with the nearest benchmark extreme T subdwarf WISEA J181006.18$-$101000.5 suggests the object would be cooler and more metal-poor than the benchmark.

The neutrino fast flavor instability dominates the evolution of neutrino flavor within the engines of core-collapse supernovae and neutron star mergers. However, theoretical models of neutrino flavor change that include many-body quantum correlations can differ starkly from similar mean-field calculations. We demonstrate for the first time that the inhomogeneous fast flavor instability is disrupted by many-body correlations using a novel tensor network framework that allows a continuous transition between mean-field and many-body results by tuning the singular value decomposition cutoff value. Generalizing the forward-scattering Hamiltonian to spatially varying conditions, we demonstrate that the timescale of flavor transformation scales logarithmically with system size, suggesting that many-body effects could outpace mean-field instabilities even for realistic neutrino numbers. Our results have significant implications for astrophysical explosion dynamics, nucleosynthesis, and observable neutrino signatures.

The XRISM Resolve X-ray spectrometer allows to gain detailed insight into gas motions of the intra cluster medium (ICM) of galaxy clusters. Current simulation studies focus mainly on statistical comparisons, making the comparison to the currently still small number of clusters difficult due to unknown selection effects. This study aims to bridge this gap, using simulated counterparts of Coma, Virgo, and Perseus from the SLOW constrained simulations. These clusters show excellent agreement in their properties and dynamical state with observations, thus providing an ideal testbed to understand the processes shaping the properties of the ICM. We find that the simulations match the order of the amount of turbulence for the three considered clusters, Coma being the most active, followed by Perseus, while Virgo is very relaxed. Typical turbulent velocities are a few $\approx100$ km s$^{-1}$, very close to observed values. The resulting turbulent pressure support is $\approx1\%$ for Virgo and $\approx 3-4\%$ for Perseus and Coma within the central $1-2\%$ of $R_{200}$. Compared to previous simulations and observations, measured velocities and turbulent pressure support are consistently lower, in line with XRISM findings, thus indicating the importance of selection effects.

One of the possible explanations for dark matter is that of compact dark objects of baryonic origin, such as black holes or even planets. Accumulating evidence, including the discovery of merging stellar mass black holes through gravitational waves, point to a population of such objects making up at least some fraction of dark matter. We revisit a historically heavily used probe, quasar spectra, from the new perspective of time variability and gravitational lensing. From a sample of 777 quasars selected from archival data we identify 19 that show decisive evidence of lensing by compact objects with masses measured in the range $5\times 10^{-5} < M/\mathrm{M}_{\odot} < 2\times 10^{-2}$ with 99\% confidence. This is much lower than what is hoped to be detected by even the most futuristic gravitational wave detectors and analysis strategies, but is crucial for theories of compact dark matter, such as primordial black holes predicted from quantum phase transitions in the early Universe.

Recent developments suggested that planet formation occurs in regions of the discs with low turbulent viscosity. There, the dynamical corotation torque is thought to play an important role by slowing down type I migration. We aim to provide a simple analytical prescription for the dynamical corotation torque for use in 1D global models of planet formation and evolution, and assess the importance of the dynamical corotation torque for the migration of low-mass planets in low-viscosity discs. We propose simple prescriptions for calculating in 1D the time evolution of the vortensities of the librating and orbit-crossing flows around a low-mass planet, which both enter the analytical expression for the dynamical corotation torque. One of our prescriptions involves a memory timescale for the librating flow, and 2D hydrodynamical simulations of disc-planet interactions are used to assess the memory timescale and validate our model. The orbital evolution of a low-mass planet is calculated by 1D simulations where the dynamical corotation torque features our prescriptions for the vortensities of the librating and orbit-crossing flows, and by 2D hydrodynamical simulations of disc-planet interactions, assuming locally isothermal discs. We find very good agreement between the 1D and 2D simulations for a wide parameter space, whether the dynamical corotation torque slows down or accelerates inward migration. We provide maps showing how much the dynamical corotation torque reduces the classical type I migration torque as a function of planet mass and orbital distance. The reduction is about 50\% for a 10 Earth-mass planet at 10 au in a young disc with surface density profile in $r^{-1/2}$ and alpha viscosity of $10^{-4}$. In discs with low turbulent viscosity, the dynamical corotation torque should be taken into account in global models as it can strongly slow down type I migration.

We performed a systematic analysis of how the orbital configuration of a merger can influence the structural formation of remnant galaxies using 531 merger pairs selected from IllustrisTNG-100. We comprehensively described the merger orbital configuration, considering the relative orbital motion of the merger pair and their disk orientations. We quantified the galaxy structures by dynamically defining four components: bulge, disk, warm component, and hot inner stellar halo. For mergers on spiral-in orbits, the disk planes of the two merging galaxies tend to be aligned with the orbital plane, leading to higher fractions for the disk and warm components, as well as lower fractions for the bulge and hot inner stellar halo components in the remnant galaxy. For mergers on direct collision orbits, the disk planes of the two galaxies tend to be perpendicular to the orbital plane, leading to lower fractions for disk and warm components, as well as higher fractions of the bulge and hot inner stellar halo in the remnant. Mergers can lead to either an increase or decrease in the disk and bulge mass fraction in the remnant compared to the progenitor galaxy, depending on the merger orbital configurations; however, in 93% of cases, mergers cause an increase in the hot inner stellar halo. As a result, the luminosity fraction of the hot inner stellar halo (but not the bulge) in galaxies at $z=0$ is highly correlated with its total ex situ stellar mass. In addition, we find that merger on spiral-in orbits is one of the possible reasons for the formation of recently discovered red but HI-rich (RR) galaxies.

JWST MIRI 15 micron time-series eclipse photometry presents a powerful way to probe for the presence of atmospheres on low-temperature rocky exoplanets orbiting nearby stars. Here, we introduce a novel technique, frame-normalized principal component analysis (FN-PCA) to analyze and detrend these MIRI time-series observations. Using the FN-PCA technique, we perform a uniform reanalysis of the published MIRI 15 micron observations of LHS 1478 b, TOI-1468 b, LHS 1140 c, TRAPPIST-1 b, and TRAPPIST-1 c using our new data reduction pipeline (Erebus) and compare them to different potential atmospheric and surface compositions. We also investigate additional public data sets with the sole purpose of understanding the instrument systematics affecting MIRI. We identify and categorize important detector-level systematics in the observations that are generally present across all 17 analyzed eclipse observations, which we illustrate as eigenimage/eigenvalue pairs in the FN-PCA. One of these eigenimage/eigenvalue pairs corresponds to the prominent ramp effect at the beginning of the time-series observations which has widely been reported for JWST and Spitzer photometry. For JWST/MIRI, we show that the detector settling time scales exponentially with the apparent magnitude of the target star $T_\mathrm{set} \mathrm{[hours]} = 0.063\exp^{0.427\cdot m_K} -0.657$. This uniform reanalysis and investigation of JWST/MIRI systematics is done in preparation for the 500 hour Rocky Worlds DDT survey, to demonstrate a data-driven systematic model usable across all MIRI 15 micron datasets.

Spatially resolved multi-wavelength analysis is essential to study galaxy formation and evolution. A UV-bright galaxy COS-2987030247 at $z = 6.81$ is one of the Rosetta Stones in the epoch of reionization for which JWST NIRSpec Integral Field Spectroscopy, NIRCam imaging, and ALMA data are available thanks to the RIOJA program. We identified the rest-frame optical emission lines from the ionized hydrogen, oxygen, and neon gas. The [OIII] 5008Å line emission and the NIRCam images show a complex kinematical and morphological structure where two bright main and three faint clumps are identified in a 10 kpc extent. The system is not a simple rotating disk, but instead consists of merger-induced star-forming clumps, as expected from recent cosmological simulations. The spatially resolved emission line fluxes show that dust attenuation, metal enrichment, and ionization parameter are preferentially enhanced in the star formation peaks. Our SED fitting suggests that the main clumps are in a moderately dust-attenuated star forming phase ($A_{\rm V} = 0.2$--$0.3$ and SFR(H$\alpha$) $\sim 10$\,M$_\odot$\,yr$^{-1}$) with almost zero escape fraction of ionizing photons. In contrast, the sub-clumps are dust-free and lying on or below the main sequence of star-forming galaxies. These sub-clumps may work as a perturber that triggers the clumpy starburst in the surrounding gas through the merger event.

We report the discovery of a complex, two-phase ionized gas structure in the star-forming galaxy COS-2987 at $z = 6.81$, revealed by the unprecedented synergy between JWST/NIRSpec IFS and ALMA observations. These observations detect key emission lines, including [OII]$\lambda\lambda$ 3727, 3730, [OIII] 4364, [OIII] $\lambda\lambda$ 4960, 5008, [OIII] 88 $\mu$m, as well as H$\alpha$ and H$\beta$. JWST spectroscopy alone indicates interstellar medium (ISM) properties that are typical for galaxies at $z\sim7$. These include low dust extinction ($A_{\rm V} \approx 0.14$ mag), moderate electron density ($n_{\rm e} \approx 500$ cm$^{-3}$), and low gas-phase metallicity ($\sim10\%$). However, the strong far-infrared [OIII] 88 $\mu$m emission detected by ALMA cannot be explained by a single-phase ionized medium with uniform electron density and temperature. Instead, a two-phase ISM model, comprising compact, high-temperature, high-density gas components ($T_e \approx 26,000$ K; $n_e \approx 600 \mathrm{cm}^{-3}$) and an extended, cooler, lower-density component ($T_e \approx 8,000$ K; $n_e \approx 50 \mathrm{cm}^{-3}$), successfully reproduces the observed line ratios of [OIII] 88 $\mu$m/[OIII] 5008 Å and [OIII] 4364 Å/[OIII] 5008 Å, with a volume ratio of 1 : 300 between the two phases. Our results demonstrate that JWST alone probes only a fraction of the ISM and highlight the critical importance of combining JWST and ALMA to reveal the density-stratified ISM of early galaxies.

Third-generation gravitational wave (GW) observatories such as the Einstein Telescope (ET) and Cosmic Explorer (CE) will detect hundreds of thousands of binary neutron star (BNS) mergers, reaching redshifts beyond $z\sim3$. To fully exploit joint GW and electromagnetic (EM) detections, dedicated strategies and adapted EM facilities are essential. We investigate the role of Integral Field and Multi-Object Spectroscopy (IFS and MOS) with the Wide-field Spectroscopic Telescope (WST) on next generation GW multi-messenger (MM) observations. We consider simulations of BNS populations, their GW detections with ET(+CE), and their EM counterparts: kilonovae (KNe) and gamma-ray bursts (GRBs). We consider two strategies: one in synergy with wide-field photometric surveys, and a galaxy-targeted one exploiting WST high multiplexing. We estimate the number of galaxies in GW error volumes, and identify observational challenges and mitigation strategies. We find that WST can detect KNe up to $z\sim0.4$ and $m_{\mathrm{AB}}\sim25$, and GRB afterglows beyond $z>1$ for $\Theta_{\mathrm{view}}\lesssim15^\circ$. KN observations are best scheduled 12-24 hours post-merger. For poorly localised GRBs, WST IFS can aid the identification. Mini-IFUs and galaxy catalogues complete to $z\leq0.5$ are key to EM counterpart detection. Even at low $z$, the number of galaxies can be huge-thousands at $z<0.1$, tens of thousands at $z<0.2$. Events at $z<0.3$ with localisation $<10$deg$^{2}$ are golden cases for WST, requiring few exposures to target all galaxies. Detecting and characterising EM counterparts of BNS detected in the extended volume explored by next-generation interferometers will be challenging. We show that high-sensitivity, wide-field, high-multiplex spectroscopic facilities are powerful instruments to fully exploit the new multi-messenger science opportunities enabled by next generation GW detectors.

The current state-of-the-art cosmic distance ladder requires three rungs--geometric distances, primary indicators, and Type Ia Supernovae--to achieve a 1% measurement of the Hubble constant $H_0$. The Habitable Worlds Observatory will have the sensitivity and resolution to reduce this to a two-step measurement, eliminating the third rung entirely and reaching into the Hubble flow with stellar distance indicators such as Cepheid variables and the tip of the red giant branch alone. We discuss the requirements for a program to measure $H_0$ to 1% with HWO here, including telescope and instrument design considerations. We also comment on the potential of HWO to measure distances to low-mass dwarf galaxies via their RR Lyrae stars.

The vast amount of energy released by active galactic nuclei (AGN) is increasingly recognized as a key driver of evolution not only in massive galaxies and clusters, but also in low-mass dwarf galaxies. Despite this, their role in the early stages of galaxy formation and in self-regulating the rapid growth of the first and abundant supermassive black holes (SMBHs) remains poorly understood. Through new high-resolution zoom-in cosmological simulations, we follow the co-evolution of $10^5 M_\odot$ black hole seeds with their host galaxy. The simulated suite progressively spans physics ranging from no AGN feedback and Eddington-limited thermal feedback, to more complex setups including non-Eddington-limited thermal, kinetic and radiative feedback. Across all our models, we find that black hole seeds efficiently reach masses of $\sim10^7 M_\odot$ by z=8. Although they exhibit notably different mass growth histories, these latter seem unimpeded by the presence of AGN feedback. The simulation including radiative feedback is the most distinct, with super-Eddington episodes driving fast and mass-loaded gas outflows (exceeding 2500 km $s^{-1}$) up to $\sim$50 kpc, along with minor stellar mass suppression in the host galaxy. Our measurements are in broad agreement with moderate luminosity quasars recently observed by JWST, producing overmassive black holes, dynamical masses of $\sim10^{9.5} M_\odot$, and high, though short-lived, Eddington fraction accretion rates. These results advocate for a scenario where AGN feedback allows for rapid SMBH growth during the reionisation era, while driving winds that extend deep into the intergalactic medium - shaping host galaxies as well as more distant surroundings.

Several observational phenomena suggest that the standard model of cosmology and particle physics requires revision. To address this, we consider the extension of general relativity known as massive gravity (MG). In this Letter, we explore the imprints of MG on the propagation of gravitational waves (GWs): their modified dispersion relation and their additional (two vector and one scalar) polarization modes on the stochastic GW background (SGWB) detected by pulsar timing arrays (PTAs). We analyze the effects of massive GWs on the Hellings-Downs curve induced by modification of the overlap reduction function. Our study consists of analyzing observational data from the NANOGrav 15-year dataset and the Chinese PTA Data Release I, and is independent of the origin of the SGWB (astrophysical or cosmological). By considering the bound on the graviton mass imposed through the dispersion relation, we scrutinize the possibility of detecting traces of MG in the PTA observational data. We find that massive GWs predict better fits for the observed pulsar correlations. Future PTA missions with more precise data will hopefully be able to detect the GW additional polarization modes and might be effectively used to constrain the graviton mass.

Most of the visible mass in a typical spiral galaxy is distributed in a thin disk, with a radial extent much larger than its thickness. While the planar disk structure, including non-axisymmetric features such as spiral structure, has been studied extensively, the vertical structure has not received comparable attention. This review aims to give a comprehensive, pedagogic introduction to the rich topic of vertical structure of a galactic disk in hydrostatic equilibrium and discuss the theoretical developments in this field in the context of recent observations. A realistic multi-component disk plus halo model of a galaxy has been developed and studied by us in detail. This takes account of both stars and interstellar gas, treated as isothermal components with different velocity dispersions, which are gravitationally coupled; further, the disk is in the gravitational field of the dark matter halo. This review focuses on this model and the results from it in different physical cases. The gas and halo crucially affect the resulting self-consistent stellar distribution such that it is vertically constrained to be closer to the mid-plane and has a steeper profile than in the standard one-component case, in agreement with modern observations. A typical stellar disk is shown to flare by a factor of few within the visible radial extent of the disk. These robust results question the sech^2 profile and a constant scale height, routinely used in the literature, for convenience. In an important application, the observed HI gas scale height is used as a constraint on the model which helps determine the shape and the density profile of the dark matter halo for galaxies. Finally, we outline some key, open questions which can be addressed in the near future using the above model, and new observational data -- for example, from IFU surveys and JWST -- for a better understanding of this topic.

We present a robust, apples-to-apples comparison between the photometric properties of the intracluster light (ICL) in the TNG300 magnetohydrodynamic cosmological simulation and those in Wendelstein Wide Field Imager (WWFI) observations. This is accomplished by generating synthetic $g'$-band images of 40 massive ($\log\left(M_{\rm 200, crit}/{\rm M}_{\odot}\right) > 14.5$) TNG300 clusters at $z \approx 0.06$, closely mimicking WWFI observations, and then performing identical photometric calculations on the synthetic and real images. Importantly, we apply the same observationally motivated satellite-masking procedure to both data-sets, which effectively removes any possible biases introduced by the halo finder. We first analyze the light distribution of the `smooth' stellar component of each cluster, composed of the brightest cluster galaxy (BCG) plus the ICL, and find that it tends to be about twice as extended in TNG300 than in observations, while also being approximately 1 $g'$ mag arcsec$^{-2}$ brighter. We then quantify $f_{\rm ICL}$, the ICL fraction relative to the BCG+ICL, by considering several ICL definitions: (i) the light dimmer than a surface brightness cut at 27 $g'$ mag arcsec$^{-2}$, (ii) the excess light over a de Vaucouleurs profile, (iii) the light beyond twice the half-light radius ($2 r_{\rm half}$), and (iv) the light beyond a fixed circular aperture of 30, 50, or 100 kpc. For most definitions, the median $f_{\rm ICL}$ is consistent between simulation and observations. However, the observations exhibit larger scatter in $f_{\rm ICL}$, which we attribute primarily to observational uncertainties in the total BCG+ICL luminosity rather than `true' cluster-to-cluster variation in the real Universe. We also find that most methods yield median $f_{\rm ICL}$ values near 0.3, which is consistent with a BCG/ICL transition radius around $2 r_{\rm half}$.

We review solar studies using AIA, HMI, and EVE data from the SDO spacecraft that revolutionized our physical understanding of the Sun. The relevant SDO studies cover the entire 15-year lifetime of SDO, from 2010 May 1 to 2025 May 1. The discussed phenomena and their physical interpretations include (in chronological order): (1) MHD Waves and Oscillations (AIA, HMI); (2) Propagating MHD Waves (AIA); (3) Coronal Loop Cross-Sectional Temperatures (AIA); (4) Size Distributions of Solar Flare Parameters (AIA); (5) Spatio-Temporal Evolution and Diffusion (AIA); (6) The Rosner-Tucker-Vaiana (RTV) Scaling Law (AIA); (7) The Fractal-Diffusive Self-Organized Criticality Model (AIA); (8) Automated Temperature and Emission Measure Maps (AIA); (9) Automated Pattern Recognition Codes (AIA); (10) Kelvin-Helmholtz Instability in Reconnetion Outflows (AIA); (11) Hydrodstatics of Coronal Loops (AIA); (12) Magnetic Energy Dissipation (HMI); (13) Global Energetics of Solar Flares (AIA).

We have analyzed the effects of rotation on mass-radius relationships for single-layer and two-layer planets having a core and an envelope made of pure materials among iron, perovskite and water in solid phase. The numerical surveys use the DROP code updated with a modified polytropic equation-of-state (EOS) and investigate flattening parameters $f$ up to $0.2$. In the mass range $0.1 M_\oplus < M < 10 M_\oplus$, we find that rotation systematically shifts the curves of composition towards larger radii and/or smaller masses. Relative to the spherical case, the equatorial radius $R_{eq}$ is increased by about $0.36f$ for single-layer planets, and by $0.30f$ to $0.55f$ for two-layer planets (depending on the core size fraction $q$ and planet mass $M$). Rotation is an additional source of confusion in deriving planetary structures, as the radius alterations are of the same order as i) current observational uncertainties for super-Earths, and ii) EOS variations. We have established a multivariate fit of the form $R_{eq}(M,f,q)$, which enables a fast characterization of the core size and rotational state of rocky planets and ocean worlds. We discuss how the observational data must be shifted in the diagrams to self-consistently account for an eventual planet spin, depending on the geometry of the transit (circular/oblate). A simple application to the recently characterized super-Earth candidate LHS1140b is discussed.

We present an analysis on the detection viability of refraction effects in Kepler's exoplanet atmospheres using binning techniques for their light curves in order to compare against simulated refraction effects. We split the Kepler exoplanets into sub-populations according to orbital period and planetary radius, then search for out-of-transit changes in the relative flux associated with atmospheric refraction of starlight. The presence of refraction effects - or lack thereof - may be used to measure and set limits on the bulk properties of an atmosphere, including mean molecular weight or the presence of hazes. In this work, we use the presence of refraction effects to test whether exoplanets above the period-radius valley have H/He atmospheres, which high levels of stellar radiation could evaporate away, in turn leaving rocky cores below the valley. We find strong observational evidence of refraction effects for exoplanets above the period-radius valley based on Kepler photometry, however those related to optically thin H/He atmospheres are not common in the observed planetary population. This result may be attributed to signal dampening caused by clouds and hazes, consistent with the optically thick and intrinsically hotter atmospheres of Kepler exoplanets caused by relatively close host star proximity.

Various ``cosmic shorelines" have been proposed to delineate which planets have atmospheres. The fates of individual planet atmospheres may be set by a complex sea of growth and loss processes, driven by unmeasurable environmental factors or unknown historical events. Yet, defining population-level boundaries helps illuminate which processes matter and identify high-priority targets for future atmospheric searches. Here, we provide a statistical framework for inferring the position, shape, and fuzziness of an instellation-based cosmic shoreline, defined in the three-dimensional space of planet escape velocity, planet bolometric flux received, and host star luminosity; explicitly including luminosity partially circumvents the need to estimate host stars' historical X-ray and extreme ultraviolet fluences. Using Solar System and exoplanet atmospheric constraints, under the restrictive assumption that one planar boundary applies across a wide parameter space, we find the critical flux threshold for atmospheres scales with escape velocity with a power-law index of $p=6.08^{+0.69}_{-0.48}$, steeper than the canonical literature slope of $p=4$, and scales with stellar luminosity with a power-law index of $q=1.25^{+0.31}_{-0.22}$, steep enough to disfavor atmospheres on Earth-sized planets out to the habitable zone for stars less luminous than $\log_{10} (L_\star/L_\odot) = -2.22 \pm 0.21$ (roughly spectral type M4.5V). If we relax the assumption that one power law must stretch from the hottest exoplanets to the coolest Solar System worlds, the narrower question of ``Which warm planets have thick CO$_2$ secondary atmospheres?" is still poorly constrained by data but should improve significantly with planned JWST observations.

Fast Radio Bursts (FRBs) are highly energetic, extremely short-lived bursts of radio flashes. Despite extensive research, the exact cause of these outbursts remains speculative. The high luminosity, short duration, and high dispersion measure of these events suggest they result from extreme, high-energy extragalactic sources, such as highly magnetized and rapidly spinning neutron stars known as magnetars. The number of detected FRBs, including repeating ones, has grown rapidly in recent years. Except for FRB 20200428D, and FRB-like radio burst that is associated to Galactic magnetar SGR 1935+2154, no multi-wavelength counterpart to any FRB has been detected yet. The High Energy Stereoscopic System (H.E.S.S.) telescope has developed a {program} to follow up FRBs searching for their gamma-ray counterparts, helping to uncover the nature of FRBs and FRB sources. This paper provides an overview of the searches for FRB sources conducted by H.E.S.S., including follow-up observations and simultaneous multi-wavelength campaigns with radio and X-ray observatories. Among the FRB sources observed by H.E.S.S., nine are localized with redshifts ranging between 0.11 and 0.492 from 2015 to 2022. No significant very high energy (VHE) emission was detected during these observations. We report constraints on the VHE luminosity ranging from $10^{44}$ erg s$^{-1}$ and $10^{48}$ erg s$^{-1}$, placing limits on the FRB's region persistent VHE emission and potential FRB afterglow emission across timescales from hours to years.

The existence of supermassive black hole binaries (SMBHBs) is predicted by various cosmological and evolutionary scenarios for active galactic nuclei. These objects are considered as contributors into the gravitational wave (GW) background and emitters of discrete GW bursts. Yet, SMBHBs remain a rather elusive class of extragalactic objects. Previously we have identified the quasar J2102+6015 as potential SMBHB system based on its oscillating astrometric pattern. We analysed the available VLBI astrometry data and identified another case of astrometric oscillations in the source J0204+1514. We assume these oscillations as manifestations of orbital motion in a binary systems. We estimated parameters of the suspected SMBHB in this source and applied basic theoretical models for projecting its evolution toward coalescence. We also develop a toy model of SMBHB consistent with the discovered astrometric oscillations and give quantitative predictions of GW emission of such the source using the case of J0204+1514 as an example. As an ad hoc result, we also provide a re-assessed estimate of the redshift of J2102+6015, z=1.42. A toy model of the object containing SMBHB with parameters consistent with the observed astrometric oscillations of the source J0204+1514 enabled us to consider GW emission as the cause of the system's orbital evolution. Astrometric VLBI monitoring has an appreciable potential for future detections of SMBHBs as multi-messenger targets for both electromagnetic (in radio domain) and gravitational wave astronomy. To outline the contours of a future physical model connecting SMBHB with detectable GW manifestations, we apply the toy model to the source J0204+1514. We also provide justification for aiming future space-borne VLBI missions toward direct imaging of SMBHBs as a synergistic contribution into future multi-messenger studies involving prospective GW facilities.

A study of [C II] 158 micron and HI 21-cm spectroscopic images plus high velocity resolution optical and ultraviolet spectra has shown the structure of the Orion Nebula to be different from that found from the study of those data separately. The [C II] features recently identified as the [C II] Shell is shown to be part of the Veil-B HI foreground layer. Jointly called the Outer Shell, it covers the bright Huygens Region and the Extended Orion Nebula. Its maximum expansion velocity is 15 km/s. Closer to ThetaOne OrionisC there is a second expanding shell, called the Inner Shell. It has an expansion velocity of 27 km/s and probably results from a more recent period of strong wind from one or more of the Trapezium stars. Even closer to ThetaOne OrionisC there is a central high ionization bubble, freely expanding towards the observer but slowed in the opposite direction by photo-ionized gas coming off the Main Ionization Front. Utilization of spectroscopic measures of the equivalent width of Hbeta shows that the enhanced emission in [C II] seen just outside the visual wavelength boundaries of the Orion Nebula is not caused by limb-brightening of the Outer Shell. This enhanced emission is due to the radiation field of the Trapezium stars being filtered by intervening residual neutral hydrogen. A velocity component near 30 km/s (Heliocentric) first seen in HI is also present in [C II] and may result from a foreground cloud of the ISM.

A wide range of scientific imaging datasets benefit from human inspection for purposes ranging from prosaic-such as fault identification and quality inspection-to profound, enabling the discovery of new phenomena. As such, these datasets come in a wide variety of forms, with diverse inspection needs. In this paper we present a software package, Image Marker, designed to help facilitate human categorization of images. The software allows for quick seeking through images and enables flexible marking and logging of up to 9 different classes of features and their locations in files of FITS, TIFF, PNG, and JPEG format. Additional tools are provided to add text-based comments to the marking logs and for displaying external mark datasets on images during the classification process. As our primary use case will be the identification of features in astronomical survey data, Image Marker will also utilize standard world coordinate system (WCS) headers embedded in FITS headers and TIFF metadata when available. The lightweight software, based on the Qt Framework to build the GUI application, enables efficient marking of thousands of images on personal-scale computers. We provide Image Marker as a Python package, and as Mac and Windows 11 executables. It is available on GitHub or via pip installation.

The Juno spacecraft provides a unique opportunity to explore the mechanisms generating Jupiter's aurorae. Past analyses of Juno data immensely advanced our understanding of its auroral acceleration processes, however, few studies utilized multiple instruments on Juno in a joint systematic analysis. This study uses measurements from the Juno Ultraviolet Spectrograph (UVS), the Jupiter Energetic particle Detector Instrument (JEDI), and the Juno Magnetometer (MAG) from the first 20 perijoves. On magnetic field lines associated with the diffuse aurora, we consistently find small-scale magnetic field fluctuations with amplitudes of up to 100 nT on time scales of seconds to 1 minute. On magnetic field lines directly linked to the main emission, the electron distribution is field-aligned, mostly broad-band in energy, and accompanied by large-scale magnetic field perturbations of several 100 nT on time scales of tens of min (except one case). These large-scale perturbations are generally associated with quasistatic field-aligned electric currents. Small-scale magnetic fields are not resolved over the main emission zone closer than radial distances 4 Jovian radii due to the digitization limit of the magnetometer. However, in all cases where Juno crosses the main auroral field lines beyond 4RJ, the digitization limit is significantly reduced and we detect small-scale magnetic field fluctuations of 2 nT to 10 nT consistent with a turbulent spectrum. Associated energy fluxes projected to Jupiter can exceed 1000 mW/m2. The general broad-band nature of the electron distributions and the consistent presence of small-scale magnetic field fluctuations over the main emission support that wave-particle interaction can dominantely contribute to power Jupiter's auroral processes.

Counting the number and brightness of ionizing radiation sources out to a redshift of z ~ 1.2 will revolutionize our understanding of how the ionizing background is created and sustained by the embedded growth of meta-galactic structures. The shear number of sparsely separated targets required to efficiently construct redshift binned luminosity functions is industrial in scale, driving the need for low spectral resolution multi-object spectroscopy (MOS) with a short wavelength cut-off ~ 1000 Å, a sensitivity in the far-UV to better than 30 abmag, and an instantaneous field-of-view ~ (2')$^2$. A MOS on Habitable Worlds Observatory is the only instrument that could conceivably carry out such an ambitious observing program.

We present MomentEmu, a general-purpose polynomial emulator for fast and interpretable mappings between theoretical parameters and observational features. The method constructs moment matrices to project simulation data onto polynomial bases, yielding symbolic expressions that approximate the target mapping. Compared to neural-network-based emulators, MomentEmu offers negligible training cost, millisecond-level evaluation, and transparent functional forms. As a demonstration, we develop two emulators: PolyCAMB-$D_\ell$, which maps six cosmological parameters to the CMB temperature power spectrum, and PolyCAMB-peak, which enables bidirectional mapping between parameters and acoustic peak features. PolyCAMB-$D_\ell$ achieves an accuracy of $0.03\%$over $\ell \leq 2510$, while PolyCAMB-peak also reaches sub-percent accuracy and produces symbolic forms consistent with known analytical approximations. The method is well suited for forward modelling, parameter inference, and uncertainty propagation, particularly when the parameter space is moderate in dimensionality and the mapping is smooth. MomentEmu offers a lightweight and portable alternative to regression-based or black-box emulators in cosmological analysis.

We analyze the high-resolution XRISM/Resolve spectrum of the Fe K$\alpha$ emission line of the nearest active galactic nucleus, in Centaurus A. The line features two narrow and resolved peaks of Fe K$\alpha_1$, and Fe K$\alpha_2$ with a FWHM of $(4.8\pm0.2)\times10^2$ km/s each. A broad line with a FWHM of $(4.3\pm0.3)\times10^3$ km/s, and with a flux similar to the two narrow line cores, is also required. This broad component is not observed in the optical or IR spectra of Cen A. The line shape requires the existence of an emission region that extends from $\sim10^{-3}$ pc to $\sim10^1$ pc. Assuming that the emissivity follows a radial power-law profile of $r^{-q}$, we find $q\approx2$. This may indicate an extended corona, an emitting region that bends towards the corona, or a non-uniform density. When assuming $q=3$, the line shape can only be reproduced by including three emitting components in the model. The measured best-fit inclination is $24^{+13}_{-7}$ degrees, but higher inclinations are only slightly disfavored. A single blurred MYTorusL line profile can describe the line shape, but requires a large relative normalization. This could be due to past variability, modified abundances, or differing geometries. The line shape can be reproduced from the radii measured by reverberation mapping, but only if an additional extended emitting region at small radii is included.

NGC 4527 is a nearby edge-on spiral galaxy with both starburst and AGN features, hosting a LINER nucleus. We present a radio study of the large-scale structure and nuclear region of this galaxy, based on new uGMRT observations at 700 and 1230 MHz. Our continuum maps reveal extended emission tracing the stellar disk, with no evidence of a radio halo. The spectral index distribution and the presence of PAHs across the disk are consistent with ongoing star formation. In the nuclear region we resolve three compact sources: one at the galactic centre and two symmetrically aligned with the major axis at $\sim400$ pc. The spectral index values and the destruction of PAHs in the central source, together with previously detected X-rays emission, suggest the presence of a low-luminosity AGN. The two off-centre sources are consistent with a star formation ring, coincident with a molecular gas ring previously reported. We explore a scenario where super-Eddington accretion onto the black hole drives a dense wind that falls back onto the disk, triggering star formation in a circumnuclear ring.

Planets are complicated. Understanding how they work requires connecting individual objects to the context of broader populations. Exoplanets are easier to picture next to their closest Solar System archetypes, and planets in the Solar System are richer when seen alongside a growing community of known exoplanets in the Milky Way. The `exoatlas` toolkit provides a friendly Python interface for retrieving and working with populations of planets, aiming to simplify the process of placing worlds in context.

W51-IRS2 is known to be one of the most prolific sources of interstellar ammonia (NH$_3$) maser lines. So far, however, many of these inversion lines have rarely been studied. Here we report spectrally resolved line profiles for the majority of detected features and provide information on the variability of these maser components between 2012 and 2023. This includes the first tentative detection of a ($J$,$K$) = (5,2) maser in the interstellar medium and the first tentative detection of a (6,4) maser in W51-IRS2. Furthermore, we report for the first time NH$_3$ (9,6) maser emission below Local Standard of Rest velocities of 50 km s$^{-1}$ in this source as well as double maser features occasionally seen in other transitions. The detected maser lines strongly indicate vibrational pumping by $\approx$10 $\mu$m photons, which must be abundant due to the high kinetic temperature ($\approx$300 K) of the ammonia emitting gas. The detection of vibrationally excited NH$_3$, suggesting a vibrational excitation temperature consistent with the kinetic one, and a comparison with measured SiO line profiles is also presented. For the (10,7) line, we find a tentative correlation between flux density and line width compatible with unsaturated maser emission. The velocity drift of the so-called 45 km s$^{-1}$ maser features, reported to be +0.2 km s$^{-1}$ yr$^{-1}$ between 1996 and 2012, has either slowed down to values $<$0.1 km s$^{-1}$ or has entirely disappeared. In 2023, the component is only seen in ammonia inversion lines that are located at least 800 K above the ground state. The other features have faded. Possible scenarios explaining this phenomenon are discussed.

The Oblate Schwarzschild (OS) approximation is a method often used to compute the flux of X-rays emitted from a rapidly rotating neutron star. In this approximation, the oblate shape of the rotating star is embedded in the Schwarzschild metric, which is used to compute the redshift of photon energies as they propagate from the star to the telescope. In this paper, we demonstrate that there are small errors introduced by the standard treatment of photon redshift in the OS approximation and provide a simple method to correct these errors. These errors are constant in phase, so this results in a constant absolute reduction in the flux. For PSR J0740+6620, the most rapidly spinning of the pulsars observed by NICER, we estimate the flux errors are less than 1\%, which is an order of magnitude smaller than the uncertainty in the distance, so this does not affect the mass and radius constraints found for this pulsar. The errors for the other pulsars observed by NICER are even smaller. However, this correction should be included when analyzing data for more rapidly rotating X-ray pulsars with spin frequencies near 600 Hz.

Few planetary systems have measured mutual inclinations, and even less are found to be non-coplanar. Observing the gravitational interactions between exoplanets is an effective tool to detect non-transiting companions to transiting planets. Evidence of these interactions can manifest in the light curve through transit timing variations (TTVs) and transit duration variations (TDVs). Through analysis of Kepler photometry and joint TTV-TDV modeling, we confirm the detection of KOI-134 b, a transiting planet with mass and size similar to Jupiter on a period of ~67 days, and find that it exhibits high TTVs (~20-hr amplitude) and significant TDVs. We explain these signals with the presence of an innermost non-transiting planet in 2:1 resonance with KOI-134 b. KOI-134 c has a mass $M = 0.220^{+0.010}_{-0.011} M_\text{Jup}$ and a moderately-high mutual inclination with KOI-134 b of $i_\text{mut} = 15.4_{-2.5}^{+2.8}{^\circ}$. Moreover, the inclination variations of KOI-134 b are so large that the planet is predicted to stop transiting in about 100 years. This system architecture cannot be easily explained by any one formation mechanism, with other dynamical effects needed to excite the planets' mutual inclination while still preserving their resonance.

POEMMA-Balloon with Radio (PBR) is a scaled-down version of the Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) design, optimized to be flown as a payload on one of NASA's sub-orbital super pressure balloons circling the Earth above the southern oceans for a mission duration of more than 20 days. One of the main science objectives of PBR is to follow up astrophysical event alerts in search of neutrinos with very high energy ($10^8 \lesssim E_\nu /{\rm GeV} \lesssim 10^{10}$). Of particular interest for anticipated PBR observations, the KM3NeT Collaboration has recently reported the detection of the neutrino KM3-230213A with $10^{8.1} \lesssim E_\nu/{\rm GeV} \lesssim 10^{8.9}$. Such an unprecedented event is in tension with upper limits on the cosmic neutrino flux from IceCube and the Pierre Auger Observatory: for a diffuse isotropic neutrino flux there is a $3.5\sigma$ tension between KM3NeT and IceCube measurements, and about $2.6\sigma$ if the neutrino flux originates in transient sources. Therefore, if KM3-230213A was not beginner's luck, it becomes compelling to consider beyond Standard Model (BSM) possibilities which could lead to a signal at KM3NeT-ARCA but not at IceCube/Auger. We calculate the PBR horizon-range sensitivity to probe BSM physics compatible with observation at KM3NeT-ARCA and non-observation at IceCube/Auger. As an illustration, we consider a particular class of BSM physics models which has been described in the literature as a possible explanation of KM3-230213A.

Massive stars play a decisive role in the evolution of the Universe. In order to constrain their current state and structure, we need sufficiently complex models, constrained by astrometric, interferometric, and spectroscopic observations. However, they are not available for distant stars. Instead, we focused on the nearest massive stars in the Orion Belt. We obtained VLTI interferometric observations of Orion Belt stars and calibrated visibility data from the GRAVITY and PIONIER instruments. Additionally, we obtained spectroscopic data from the CFHT and CTIO observatories. For modelling, we used a modified version of PHOEBE2, extended with new interferometric and spectroscopic modules. To describe non-spherical, rotating, or Roche-like stars, integrals over triangular meshes have to be computed, using extensive grids of synthetic spectra. For fitting, we used the simplex algorithm and chi2 mapping of the parameter space. In this paper, we present single-star models of the B0Ia supergiant eps Ori. Interferometric visibilities indicate that the star is not spherical but rotating close to its critical velocity. The preferred distance, d=(384+-8)pc, corresponds to the median of distances for the Orion OB1b association. We obtained the following parameters: m=(28.4+-2.0)Msol, R=(27.6+-1.5)Rsol, Teff=25000 K, i=45deg, longitude of the ascending node, Omega=300deg, and Prot=4.3+1.0d. This compromise model provides a reasonable fit to wind-free Balmer line profiles, but there is still some tension between interferometric and spectroscopic datasets, corresponding to a faster- vs. slower-rotating star. Our fast-rotating model implies that circumstellar matter should be naturally present, in the form of wind or disk, and contribute to continuum radiation. The fast rotation of eps Ori is compatible with a merger, formed from a multiple system of comparable mass, like del, zet or sig Ori.

The Milky Way Survey of the Dark Energy Spectroscopic Instrument (DESI) has so far observed three classical dwarf spheroidal galaxies (dSphs): Draco, Sextans and Ursa Minor. Based on the observed line-of-sight velocities and metallicities of their member stars, we apply the axisymmetric Jeans Anisotropic Multi-Gaussian Expansion modeling (JAM) approach to recover their inner dark matter distributions. In particular, both the traditional single-population Jeans model and the multiple population chemodynamical model are adopted. With the chemodynamical model, we divide member stars of each dSph into metal-rich and metal-poor populations. The metal-rich populations are more centrally concentrated and dynamically colder, featuring lower velocity dispersion profiles than the metal-poor populations. We find a diversity of the inner density slopes $\gamma$ of dark matter halos, with the best constraints by single-population or chemodynamical models consistent with each other. The inner density slopes are $0.71^{+0.34}_{-0.35}$, $0.26^{+0.22}_{-0.12}$ and $0.33^{+0.20}_{-0.16}$ for Draco, Sextans and Ursa Minor, respectively. We also present the measured astrophysical J and D factors of the three dSphs. Our results indicate that the study of the dark matter content of dSphs through stellar kinematics is still subject to uncertainties behind both the methodology and the observed data, through comparisons with previous measurements and data sets.

The Large High Altitude Air Shower Observatory has detected very high-energy (VHE) gamma rays from NGC 4278, which is known to host a low-luminosity active galactic nucleus (AGN). Having only very weak radio jets, the origin of its VHE gamma rays is unclear. In this paper we first show that NGC 4278 has a massive molecular cloud surrounding the nucleus by analyzing data taken with the Atacama Large Millimeter/submillimeter Array. We then assume that cosmic ray protons are accelerated in a radiatively inefficient accretion flow around the supermassive black hole, which diffuse into the molecular cloud and produce gamma rays and neutrinos via $pp$ interactions. We model the gamma-ray spectra and find that the observations can be explained by such hadronic processes if the AGN activity was higher in the past than at present, and the diffusion coefficient in the molecular cloud is appreciably smaller than in the Milky Way interstellar medium. However, we also show that the high-energy neutrinos co-produced with the gamma rays are unlikely to be detectable even with IceCube-Gen2.

The role of magnetic fields in regulating the formation of massive stars remains much debated. Here we present sub-millimeter polarimetric observations with JCMT-POL2 at $850\:\mu$m of 13 regions of massive star formation selected from the SOFIA Massive (SOMA) star formation survey, yielding a total of 29 massive protostars. Our investigation of the $p'-I$ relationship suggests that grain alignment persists up to the highest intensities. We examine the relative orientations between polarization-inferred magnetic field direction and source column density elongation direction on small and large scales. On small scales, we find a bimodal distribution of these relative orientations, i.e., with an excess of near-parallel and near-perpendicular orientations. By applying a one-sample Kuiper test and Monte Carlo simulations to compare to a relative orientation distribution drawn from a uniform distribution, we statistically confirm this bimodal distribution, independent of the methods to measure structural orientation. This bimodal distribution suggests that magnetic fields are dynamically important on the local scales ($\lesssim 0.6\:$pc) of massive protostellar cores. We also examine how basic polarization properties of overall degree of polarization and local dispersion in polarization vector orientations depend on intrinsic protostellar properties inferred from spectral energy distribution (SED) modeling. We find a statistically significant anti-correlation between the debiased polarized fraction and the luminosity to mass ratio, $L_{\rm bol}/M_{\rm env}$, which hints at a change in the dust properties for protostellar objects at different evolutionary stages.

Relativistic pair beams created in the intergalactic medium (IGM) by TeV gamma rays from blazars are expected to produce a detectable GeV-scale electromagnetic cascade, but the cascade component is absent in the spectra of many hard-spectrum TeV-emitting blazars. One common explanation is that weak intergalactic magnetic fields deflect the electron-positron pairs away from our line of sight. An alternative possibility is that electrostatic beam-plasma instabilities drain the energy of these pairs before a cascade can develop. Recent studies have shown that beam scattering by oblique electrostatic modes leads to minimal energy loss. But these modes might be suppressed by linear Landau damping (LLD) due to MeV-scale cosmic-ray electrons in the IGM. In this work, we explore the impact of LLD on the energy-loss efficiency of plasma instabilities in pair beams associated with 1ES 0229+200. We find that LLD effectively suppresses oblique electrostatic modes, while quasi-parallel ones grow to larger amplitudes. In this way, LLD enhances the energy-loss efficiency of the instability by more than an order of magnitude.

Nearby stars offer prime opportunities for exoplanet discovery and characterization through various detection methods. By combining HCI, RV, and astrometry, it is possible to better constrain the presence of substellar companions, as each method probes different regions of their parameter space. A detailed census of planets around nearby stars is essential to guide the selection of targets for future space missions seeking to identify Earth-like planets and potentially habitable worlds. In addition, the detection and characterisation of giant planets and brown dwarfs is crucial for understanding the formation and evolution of planetary systems. We aim to constrain the possible presence of substellar companions for 7 nearby M-dwarf stars using a combination of new SPHERE/H2 HCI and archival RV and astrometric data. We investigate how combining these techniques improves the detection constraints for giant planets and brown dwarfs compared to using each method individually. For each star and each data set, we compute the mass limits as a function of semi-major axis or projected separation using standard techniques. We then use a Monte Carlo approach to assess the completeness of the companion mass / semi-major axis parameter space probed by the combination of the three methods, as well as by the three methods independently. Our combined approach significantly increases the fraction of detectable companions. Although no new companion was detected, we could place stronger constraints on potential substellar companions. The combination of HCI, RV and astrometry provides significant improvements in the detection of substellar companions over a wider parameter space. Applying this approach to larger samples and lower-mass companions will help constraining the search space for future space missions aimed at finding potentially habitable or even inhabited planets.

Massive black hole binaries (MBHBs) and other sources within the frequency band of spaceborne gravitational wave observatories like the Laser Interferometer Space Antenna (LISA), Taiji and Tianqin pose unique challenges, as gaps and glitches during the years-long observation lead to both loss of information and spectral leakage. We propose a novel data imputation strategy based on Kalman filter and smoother to mitigate gap-induced biases in parameter estimation. Applied to a scenario where traditional windowing and smoothing technique introduce significant biases, our method mitigates the biases and demonstrates lower computational cost compared to existing data augmentation techniques such as noise inpainting. This framework presents a new gap treatment approach that balances robustness and efficiency for space-based gravitational wave data analysis.

The super-soft source (SSS) phase of a nova eruption, observed a few days after the outburst, usually displays an absorbed X-ray thermal continuum with absorption features, emitted by the white dwarf (WD) atmosphere. However, the X-ray spectra of many novae in this phase display additional emission lines which likely originate from shocks in the novae ejecta. When the shocked plasma interacts with cold gas, narrow radiative recombination continua (RRCs) and charge exchange (CX) emission are observed. We present the analysis of high-resolution ChandraLETG X-ray grating spectra of Nova Sco 2023, observed 128 and 183 days after the optical peak, on 2023 August and October. At both epochs, the absorbed X-ray thermal continuum is well described by a Non-Local Thermal Equilibrium atmosphere model with a temperature T=750,000 K (kT = 65 eV). On day 128, the atmosphere is found to be outflowing at v=-3500 km s^-1. On day 183, the atmosphere brightened by a factor of ~2 and slowed down to v=-1500$ km s^-1. The discrete emission features of the spectrum consist of the C^+5, N^+5, and N^+6 RRCs, indicating a cold electron temperature of kT_e=1 eV on day 128, and kT_e=20 eV on day 183. The observed line series of H-like and He-like C^+5, N^+5, N^+6, and O^+6 show enhanced intensities of high-n (principal quantum number) transitions, consistent with a CX model of hot ions at kT~100 eV. The velocity shift of the CX lines remained at v=+-3000 km s^-1, which can be explained by a bipolar outflow. After Nova Ret 2020 (YZ Ret), Nova Sco 2023 is yet another nova in which we have found exquisite evidence of CX in astrophysical ionized plasma.

The accretion disks of active galactic nuclei (AGN) are widely considered the ideal environments for binary black hole (BBH) mergers and the only plausible sites for their electromagnetic (EM) counterparts. Graham et al.(2023) identified seven AGN flares that are potentially associated with gravitational-wave (GW) events detected by the LIGO-Virgo-KAGRA (LVK) Collaboration during the third observing run. In this article, utilizing an additional three years of Zwicky Transient Facility (ZTF) public data after their discovery, we conduct an updated analysis and find that only three flares can be identified. By implementing a joint analysis of optical and GW data through a Bayesian framework, we find two flares exhibit a strong correlation with GW events, with no secondary flares observed in their host AGN up to 2024 October 31. Combining these two most robust associations, we derive a Hubble constant measurement of $H_{0}= 72.1^{+23.9}_{-23.1} \ \mathrm{km \ s^{-1} Mpc^{-1}}$ and incorporating the multi-messenger event GW170817 improves the precision to $H_{0}=73.5^{+9.8}_{-6.9} \ \mathrm{km \ s^{-1} Mpc^{-1}}$. Both results are consistent with existing measurements reported in the literature.

We investigated the widths and profiles of cosmic filaments using the IllustrisTNG simulations. Filaments were identified with DisPerSE, using galaxy samples in simulations as input. Since the width of an individual filament can vary significantly along its spine, we divided each filament into segments with lengths between $1.5\,h^{-1}\,\mathrm{Mpc}$ and $2.5\,h^{-1}\,\mathrm{Mpc}$ and measure their properties. The typical width of these filament segments increases gradually from approximately $0.3\,\mathrm{Mpc}$ at redshift $z = 2.0$ to about $1.0-1.5\,\mathrm{Mpc}$ at $z = 0.0$. We find that the segment width correlates nearly linearly with the linear halo mass density, consistent with previous studies. A similar linear relation is observed between the segment width and the linear stellar mass density, providing a potential estimator for filament width. Furthermore, the density profiles of filaments with different widths exhibit self-similarity and can be described by a unified formula akin to the isothermal $\beta$-model. For segments with a given width, the rescaled density profiles show only mild evolution from $z = 2.0$ to $z = 0.0$. Within the filament width, the gas temperature decreases slowly from the center to the boundary, with thicker filaments generally containing hotter gas than thinner ones. These trends in filament width, density, and thermal profiles are consistently observed across the TNG50, TNG100, and TNG300 simulations, and align well with results from earlier studies. We briefly discuss the potential implications and applications of our findings.

We present recent progress in our understanding of the physical interaction mechanisms at work in evolved binaries of low-to-intermediate initial mass, which are surrounded by a stable disc of gas and dust. These systems are known as post-asymptotic giant-branch (post-AGB) binaries, but recently, it has been shown that some systems are too low in luminosity and should be considered as post-red-giant branch (post-RGB) instead. While the systems are currently well within their Roche lobe, they still show signs of active ongoing interaction between the different building blocks. We end this contribution with some future research plans.

We investigate the clustering of Lyman-break galaxies at redshifts of 3 $\lesssim z \lesssim$ 5 within the COSMOS field by measuring the angular two-point correlation function. Our robust sample of $\sim$60,000 bright ($m_{\rm UV}\lesssim 27$) Lyman-break galaxies was selected based on spectral energy distribution fitting across 14 photometric bands spanning optical and near-infrared wavelengths. We constrained both the 1- and 2-halo terms at separations up to 300 arcsec, finding an excess in the correlation function at scales corresponding to $<20$ kpc, consistent with enhancement due to clumps in the same galaxy or interactions on this scale. We then performed Bayesian model fits on the correlation functions to infer the Halo Occupation Distribution parameters, star formation duty cycle, and galaxy bias in three redshift bins. We examined several cases where different combinations of parameters were varied, showing that our data can constrain the slope of the satellite occupation function, which previous studies have fixed. For an $M_{\rm{UV}}$-limited sub-sample, we found galaxy bias values of $b_g=3.18^{+0.14}_{-0.14}$ at $z\simeq3$, $b_g=3.58^{+0.27}_{-0.29}$ at $z\simeq4$, $b_g=4.27^{+0.25}_{-0.26}$ at $z\simeq5$. The duty cycle values are $0.62^{+0.25}_{-0.26}$, $0.40^{+0.34}_{-0.22}$, and $0.39^{+0.31}_{-0.20}$, respectively. These results suggest that, as the redshift increases, there is a slight decrease in the host halo masses and a shorter timescale for star formation in bright galaxies, at a fixed rest-frame UV luminosity threshold.

Context. Galaxies fly inside galaxy clusters and ram pressure by the ICM can remove a large amount of the ISM from the galaxy, and deposit the gas in the ICM. The ISM decoupled from the host galaxy leaves a long trail following the moving galaxy. Such long trails track the galaxy motion and can be detected with sensitive data in Halpha. Aims. We study the Halpha tail trailing NGC 4569 in the Virgo cluster. Methods. The initial discovery was made with the deep Halpha imaging data with CFHT, from the VESTIGE project. The follow-up spectroscopic observations were made with APO/DIS, MMT/Binospec and CFHT/SITELLE. Results. Besides the known 80 kpc Halpha tail downstream of NGC 4569, the deep Halpha imaging data allow the Halpha tail detected to at least 230 kpc from the galaxy. More importantly, the Halpha clumps implied from the imaging data are confirmed with the spectroscopic data. The Halpha clumps show a smooth radial velocity gradient across about 1300 km/s, eventually reaching the velocity of the cluster. We build a simple model to explain the deceleration of stripped clumps and constrain the age to about 0.9 Gyr. Conclusions. This discovery, for the first time, demonstrates the full deceleration process of the stripped ISM. This discovery also showcases the potential with wide-field Halpha survey on galaxy clusters to discover intracluster optical emission-line clouds originated from cluster galaxies. These clouds provide kinematic tracers to the infall history of cluster galaxies and the turbulence in the ICM. They are also excellent multi-phase objects to study the relevant important physical processes.

"Planetesimal or pebble" is one of the most fundamental open questions in planet formation theory. Similarly, "satellitesimal or pebble" remains unsettled regarding the formation of the Galilean satellites. I focus on a unique characteristic of Callisto--its interior is estimated to be only partially differentiated based on gravitational field measurements. I robustly demonstrate that such a state is not achievable through satellitesimal accretion, which inevitably leads to significant differentiation, but can be maintained with pebble accretion. Pebbles can release their impact energy at the surface of the satellite, allowing efficient radiative cooling, and their impact velocities can be reduced by aerodynamic drag from the circumplanetary gas disk. If future missions such as JUpiter ICy moons Explorer (JUICE) confirm that Callisto is indeed only partially differentiated, it could provide the first observed evidence for the pebble accretion mechanism--not only in the context of satellite formation, but also in the broader framework of planet formation.

Galaxies moving through the gas of the intracluster medium (ICM) experience ram pressure stripping, which can leave behind a gas tail. When a disk galaxy receives the wind edge-on, however, the characteristic signature is not a typical jellyfish tail, but rather an unwinding of the spiral arms. We aim to quantify such asymmetries both in the gas and in the stellar component of a simulated galaxy. To this end, we simulate a gas-rich star-forming spiral galaxy moving through a self-consistent ICM gas. The amplitude and location of the asymmetries were measured via Fourier decomposition. We found that the asymmetry is much more evident in the gas component, but it is also measurable in the stars. The amplitude tends to increase with time and the asymmetry radius migrates inwards. We found that, when considering the gas, the spiral arms extend much further and are more unwound than the corresponding stellar arms. Characterizing the unwinding via simulations should help inform the observational criteria used to classify ram pressure stripped galaxies, as opposed to asymmetries induced by other mechanisms.

We present a novel approach for assessing the muon content of air showers with large zenith angles on a combined analysis of their radio emission and particle footprint. We use the radiation energy reconstructed by the Auger Engineering Radio Array (AERA) as an energy estimator and determine the muon number independently with the water-Cherenkov detector array of the Pierre Auger Observatory, deployed on a 1500 m grid. We focus our analysis on air showers with primary energy above 4 EeV to ensure full detection efficiency. Over approximately ten years of accumulated data, we identify a set of 40 high-quality events that are used in the analysis. The estimated muon contents in data are compatible with those for iron primaries as predicted by current-generation hadronic interaction models. This result can be interpreted as a deficit of muons in simulations as a lighter mass composition has been established from Xmax measurements. This muon deficit was already observed in previous analyses of the Auger Collaboration and is confirmed using hybrid events that include radio measurements for the first time.

The formation of planetesimals from cm-sized pebbles in protoplanetary disks faces significant barriers, including fragmentation and radial drift. We identify a previously unaccounted screening force, arising from mutual shielding of thermal gas particles between pebbles when their separation falls below the gas mean free path. This force facilitates pebble binding, overcoming key growth barriers under turbulent disk conditions. Unlike conventional mechanisms, screening forces operate independently of surface adhesion and complement streaming instability and pressure traps by enhancing aggregation in high-density regions. Our analysis predicts that screening interactions are most effective in the {middle disk regions ($ \sim 0.3$ to few AU),} consistent with ALMA observations (e.g., TW Hya) of enhanced dust concentrations. {Furthermore, we find that screening-induced pebble growth from centimeter to kilometer scales can occur on timescales significantly shorter than the disk lifetime ($\sim 10^5$ years). Importantly, this growth naturally terminates when particles smaller than the local gas mean free path are depleted, thereby avoiding runaway accretion.} Beyond planetary science, the screening forces have {potential} implications for high-energy astrophysics, dusty plasmas, confined particle suspensions and other relevant areas, suggesting a broader fundamental significance.

X-ray reflection spectroscopy has evolved as one of the leading methods to measure black hole spins. However, the question is whether its measurements are subjected to systematic biases, especially considering the possible discrepancy between the spin measurements inferred with this technique and those from gravitational wave observations. In this work, we use general relativistic magnetohydrodynamic (GRMHD) simulations of thin accretion disks around spinning black holes for modeling the accretion process, and then we simulate NuSTAR observations to test the capability of modern reflection models in recovering the input spins. For the first time, we model the electron density and ionization profiles from GRMHD-simulated disks. Our study reveals that current reflection models work well only for fast-rotating black holes. We model the corona as the base of the jet and we find that reflection models with lamppost emissivity profiles fail to recover the correct black hole spins. Reflection models with broken power-law emissivity profiles perform better. As we increase the complexity of the simulated models, it is more difficult to recover the correct input spins, pointing towards the need to update our current reflection models with more advanced accretion disks and coronal geometries.

Magnetic activity effects are typically neglected in asteroseismic modelling of solar-type stars, presuming that these effects can be accounted for in the parametrisation of the surface effects. It was however demonstrated that magnetic activity can have a significant impact on the asteroseismic characterisation using both forward and inverse techniques. We investigated whether frequency separation ratios, which are commonly used to efficiently suppress surface effects, are also able to suppress magnetic activity effects. Based on GOLF and BiSON observations of the Sun-as-a-star, we performed asteroseismic characterisations using frequency separation ratios as constraints to measure the apparent temporal evolution of the stellar parameters and their correlation with the 10.7 cm radio flux. Frequency separation ratios do not suppress the effects of magnetic activity. Both $r_{01}$ and $r_{02}$ ratios exhibit a clear signature of the magnetic activity cycle. Consequently, when these ratios are employed as constraints in asteroseismic modelling, magnetic activity effects are propagated to the stellar characterisation. Additionally, most stellar parameters correlate with the activity cycle, unlike the direct fitting of individual frequencies. Magnetic activity effects significantly impact asteroseismic characterisation, regardless of whether forward modelling or inverse methods are used. Standard techniques to suppress surface effects have proven ineffective against magnetic activity influences and systematic uncertainties of 4.7%, 2.9%, and 1.0% should be considered for the stellar age, mass, and radius, respectively. In preparation for future space-based photometry missions, it is therefore essential to enhance our theoretical understanding of these effects and develop a modelling procedure capable of accounting for or efficiently suppressing them.

Spherically symmetric accretion incorporating self-gravity constitutes a three-point boundary value problem (TPBVP) governed by constraints at the outer boundary, sonic point, and accretor surface. Previous studies have two limitations: either employing an incorrect formula for self-gravity potential in analytical treatments, or introducing additional input parameters in numerical implementations to circumvent solving the full TPBVP. To address these issues, we present a self-consistent TPBVP formulation, solved using the relaxation method. We also derive approximate analytical formulae that enable rapid estimates of self-gravity effects. Our analysis identifies a dimensionless parameter $\beta \equiv 2G \bar{\rho} r_\mathrm{out}^2/a_\mathrm{out}^2$ that characterizes the strength of self-gravity, where $\bar{\rho}$ and $r_\mathrm{out}$ are the mean density and outer radius of the flow, respectively, and $a_\mathrm{out}$ is the adiabatic sound speed of the external medium. For practical estimation, $\bar{\rho}$ may be approximated by the external medium density $\rho_\mathrm{out}$. We identify an upper limit for $\beta$, beyond which steady accretion becomes unsustainable -- a behavior consistent with classical gravitational instability that previous studies failed to capture. The accretion rate enhancement decreases monotonically as the adiabatic index $\gamma$ increases. For $\gamma=5/3$, self-gravity ceases to augment the accretion rate. These theoretical predictions are validated by our numerical solutions. We further apply our results to two astrophysical scenarios: hyper-Eddington accretion onto supermassive black hole seeds in the early Universe, where self-gravity is significant; and accretion onto stellar-mass objects embedded in active galactic nuclei (AGN) disks, where self-gravity is non-negligible under certain conditions and should be evaluated using $\beta$.

Although the differential rotation rate on the solar surface has long been studied using optical and extreme ultraviolet (EUV) observations, associating these measurements to specific atmospheric heights remains challenging due to the temperature-dependent emission of tracers observed in EUV wavelengths. Radio observations, being primarily influenced by coherent plasma processes and/or thermal bremsstrahlung, offer a more height-stable diagnostic and thus provide an independent means to test and validate rotational trends observed at other EUV wavelengths. We aim to characterize the differential rotation profile of the upper chromosphere using cleaned solar full-disc 17 GHz radio imaging from the Nobeyama Radioheliograph (NoRH), spanning a little over two solar cycles (1992 - 2020). A tracer-independent method based on automated image correlation was employed on daily full-disc 17 GHz radio maps. Our results suggest that the upper chromosphere rotates significantly faster than the photosphere at all latitudes, with a relatively flatter latitudinal profile. A very weak anti-correlation between the equatorial rotation rate and solar activity is also observed. Our findings reaffirm the potential of radio observations to probe the dynamics of the solar chromosphere with reduced height ambiguity. The overlap of the equatorial rotation rate found in this study with that for $304$ Å in the EUV regime lends additional support to the view that the equatorial rotation rates increase with height above the photosphere. Future coordinated studies at wavelengths with better-constrained height formation will be crucial for further understanding the complex dynamics of the solar atmosphere.

Probing the nature of dark matter (DM) remains an outstanding problem in modern cosmology. The 21 cm signal, as a sensitive tracer of neutral hydrogen during cosmic dawn, provides a unique means to investigate DM nature during this critical epoch. Annihilation and decay of DM particles, as well as Hawking radiation of primordial black holes (PBHs), can modify the thermal and ionization histories of the early universe, leaving distinctive imprints on the 21 cm power spectrum. Therefore, the redshifted 21 cm power spectrum serves as a powerful tool to investigate such DM processes. In this work, we systematically assess the potential of the upcoming Square Kilometre Array (SKA) to constrain DM and PBH parameters using the 21 cm power spectrum. Assuming $10,000$ hours of integration time, the SKA is projected to reach sensitivities of $\langle\sigma v\rangle \leq 10^{-28}\,{\rm cm}^{3}\,{\rm s}^{-1}$ and $\tau\geq 10^{28}\,{\rm seconds}$, for $10\,{\rm GeV}$ DM particles. It can also probe PBHs with masses of $10^{16}\,\mathrm{g}$ and abundances $f_{\mathrm{PBH}} \leq 10^{-6}$. These results indicate that the SKA could place constraints on DM annihilation, decay, and PBH Hawking radiation that are up to two to three orders of magnitude stronger than current limits. Furthermore, the SKA is expected to exceed existing bounds on sub-GeV DM and to probe Hawking radiation from PBHs with masses above $10^{17}\,{\rm g}$, which are otherwise inaccessible by conventional cosmological probes. Overall, the SKA holds great promise for advancing our understanding of both DM particles and PBHs, potentially offering new insights into the fundamental nature of DM.

Small low-density exoplanets are sculpted by strong stellar irradiation but their origin is unknown. Two competing scenarios are that they formed either with rocky interiors and H2-He atmospheres (`gas-dwarfs'), or instead with bulk compositions dominated by H2O phases (`water-worlds'). Observations of L 98-59 d have revealed its unusually low density and spectral characteristics, opening a window for disentangling its origin. We constrain the possible range of evolutionary histories linking L 98-59 d's birth conditions to these observations, using a coupled atmosphere-interior modelling framework. Observations are explained by a chemically-reducing mantle and early substantial (>1.8 mass%) sulfur and hydrogen content, inconsistent with both the gas-dwarf and water-world origin scenarios. Observed spectral features are explained by in-situ photochemical production of SO2 in an H2 background. L 98-59 d's interior comprises a permanent magma ocean, allowing long-term retention of volatiles within its interior over billions of years, consistent with surveyed Kepler Survey trends. Modelling of the `cosmic shoreline' suggests that early volatile-metal budget is a determinant of atmospheric retention - here, we explicitly show that large initial budgets can explain JWST observations of some super-Earth planets (including DDT targets). In breaking with the supposed gas-dwarf/water-world dichotomy, we reveal an evolutionary pathway defined by high-molar-mass atmospheres overlying magma oceans, and invite a more nuanced taxonomy of small exoplanets.

We report an observation of a transit of the hot Jupiter (HJ) KELT-23A b with the Keck Planet Finder spectrograph and a measurement of the sky-projected obliquity ($\lambda$) of its Sun-like ($T_{\rm eff} \approx 5900$ K) host star. We measured a projected stellar obliquity of $\lambda \approx 180^\circ$, indicating that the orbit of the HJ is retrograde relative to the direction of the stellar spin. Due to the slow sky-projected rotational velocity of the host star ($v \sin{i_\star} \approx 0.5$ km s$^{-1}$), the true orbit of the HJ could be closer to polar. HJs around stars with effective temperatures below the Kraft break -- such as KELT-23A -- are generally found to have prograde orbits that are well-aligned with the equatorial planes of their host stars (i.e., $\lambda \sim 0^\circ$), most likely due to spin-orbit realignment driven by stellar tidal dissipation. This system is therefore a unique outlier that strains migration and tidal theories. The fact that the HJ has a highly misaligned orbit may suggest that the planet arrived at its close-in orbit relatively recently, possibly via interactions with the wide-separation (570 AU) M-dwarf companion in the system, or that it has stalled near an antialigned or polar orientation while realigning. Using Gaia DR3, we determined the orbit of the stellar companion to be moderately face-on ($\gamma = 60 \pm 4^\circ$). We show that the distribution of observed systems in the $\gamma - \lambda$ plane can be broadly reproduced using a toy model in which the orbits of the planetary and stellar companions begin aligned with the equatorial plane of the primary star and, upon migrating inwards, the planet preferentially obtains either an aligned or polar orbit.

We consider a model of inflation in which the inflaton field is a rolling axion with a potential which is flat enough to support an intermediate phase of USR inflation. Because of the Chern-Simons interaction, one polarization of the gauge field experiences the tachyonic growth during the first slow-roll stage, inducing large corrections in curvature perturbations via the inverse decay effect. A non-trivial feature of our setup is that once the system enters the USR phase, the instability parameter falls off rapidly, terminating the gauge field production efficiently. Consequently, the power spectrum involves two separate peaks, the first peak is induced by the gauge field particles production while the second peak is due to standard USR mechanism. We show that the power spectrum at the intermediate scales develops strong scale-dependence $\propto k^m$ with the index $m >4$. Calculating the bispectrum, we demonstrate that non-Gaussianities with non-trivial shapes and multiple peaks are generated in this setup.

We investigate cosmic structure formation in the framework of a path-integral formulation of an $N$-particle ensemble in phase space, dubbed Resummed Kinetic Field Theory (RKFT), up to one-loop perturbative order. In particular, we compute power spectra of the density contrast, the divergence and curl of the momentum density and arbitrary $n$-point cumulants of the stress tensor. In contrast to earlier works, we propose a different method of sampling initial conditions, with a Gaussian initial phase-space density. Doing so, we exactly reproduce the corresponding results from Eulerian standard perturbation theory (SPT) at one-loop order, showing that formerly found deviations can be fully attributed to inconsistencies in the previous sampling method. Since, in contrast to SPT, the full phase-space description does not assume a truncation of the Vlasov hierarchy, our findings suggest that non-perturbative techniques are required to accurately capture the physics of cosmic structure formation.

The potentially distinct formation pathways of exoplanets and brown dwarfs may be imprinted in their elemental and isotopic ratios. This work is part of the ESO SupJup Survey, which aims to disentangle the formation pathways of super-Jupiters and brown dwarfs through the analysis of their chemical and isotopic ratios. In this study, we characterize the atmospheres of two young L4 dwarfs, 2MASS J03552337+1133437 (2M0355) and 2MASS J14252798-3650229 (2M1425), in the AB Doradus Moving Group. This involved constraining their chemical composition, $^{12}$CO/$^{13}$CO ratio, pressure-temperature profile, surface gravity, and rotational velocity. We have obtained high-resolution CRIRES+ K-band spectra of these brown dwarfs, which we analyzed with an atmospheric retrieval pipeline. Atmospheric models were generated with the radiative transfer code petitRADTRANS, for which we employed a free and equilibrium chemistry approach, which we coupled with the PyMultiNest sampling algorithm to determine the best fit. We report robust detections of $^{13}$CO (13.4 & 8.0 $\sigma$) and HF (11.6 & 15.8 $\sigma$) in 2M0355 and 2M1425, respectively. The objects have similar overall atmospheric properties, including $^{12}$CO/$^{13}$CO isotope ratios of $95.5\substack{+6.8 \\ -6.4}$ for 2M0355 and $109.6\substack{+10.6 \\ -9.6}$ for 2M1425. In both brown dwarfs, we find only tentative hints of H$_2^{18}$O (1.1 & 3.0 $\sigma$), with lower limits of H$_2^{16}$O/H$_2^{18}$O $\sim$ 1000. Both objects appear to be close to chemical equilibrium, considering the main spectral contributors. Future studies will put them into the context of other objects observed as part of the ESO SupJup Survey.

Blazars are a subclass of active galactic nuclei (AGN) that emit non-thermal radiation through relativistic jets, characterized by rapid flux and polarization variability. Extreme high synchrotron-peaked blazars (EHSPs), with synchrotron peaks exceeding 10$^{17}$ Hz, are essential for understanding the full range of blazar phenomena and testing jet physics models. However, the number of known extreme blazars is small, so this class of objects remains poorly studied. This work aims to systematically identify and characterize the most extreme $\gamma$-ray blazars using data from the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope. The focus is on spectral hardening, where the $\gamma$-ray spectrum becomes harder at higher energies, particularly during flaring episodes. This represents the first dedicated analysis of spectral hardening across a population of EHSPs, as previous studies explored it only in individual sources. We analyze 138 blazars selected from the 4FGL-DR2 catalog with high synchrotron peak frequencies and well-sampled light curves. Flaring periods are automatically identified, and each flare is analyzed, with the significance of spectral hardening assessed through a test statistic based on the likelihood ratio of two spectral models. We identify two flaring episodes with indications of spectral hardening, in 4FGL J0238.4$-$3116 and PKS 2155$-$304, the latter detected independently by both methods but referring to the same period. These events are consistent with expectations from statistical fluctuations, suggesting that spectral hardening is a rare occurrence (< 0.1 %). These results constrain its frequency and support a smoothly varying power-law blazar emission model, motivating future multi-wavelength studies to clarify whether these rare flares reflect distinct physical processes within blazar jets.

Parity-violation leaves tell-tale trails in many cosmological observables. We illustrate parity-odd primordial scalar trispectra, that despite being of modest strength, impart detectable chirality to scalar-induced gravitational waves (SIGW). This allows us to impose strong bounds on the parity-odd part of trispectrum. Over certain scales, we find SIGW directly quantify parity-violation in primordial non-Gaussianity, unobscured by the Gaussian contribution. Our results call for treatment of SIGW and parity-odd trispectrum as complementary predictions of parity-violating theories.

We present a cautionary assessment of the extended X-ray emission around PSR B0656+14 in eROSITA Cal-PV data in response to the work of Niu et al. 2025 (arXiv:2501.17046). The eROSITA PSF model is known to underestimate emission in the wings beyond 1'. This prevents a reliable detection of faint nebular emission around PSR B0656+14 as claimed in arXiv:2501.17046. In addition, spectral analysis shows the surrounding diffuse X-rays can be fitted with the same 2BB+PL model as the pulsar's emission itself. This strongly invalidates the interpretation by the authors of arXiv:2501.17046 that the X-ray emission in the (4-10)' region is associated with the degree-scale gamma-ray halo recently found by the High-Altitude Water Cherenkov Observatory (HAWC), and shows that it originates from the pulsar due to the wings of the PSF.

We report initial observations aimed at the characterization of a third interstellar object candidate. This object, 3I/ATLAS -- also C/2025 N1 (ATLAS) -- , was discovered on 2025 July 1 UT and has an orbital eccentricity of $e\sim6.2$, perihelion of $q\sim 1.35$ au, inclination of $\sim175^\circ$, and hyperbolic velocity of $V_\infty\sim 60$ km s$^{-1}$. 3I/ATLAS has an absolute magnitude of $H_V\sim12$, which corresponds to a nuclear radius of $\sim10\text{ km}$, assuming an asteroid-like albedo of $p\sim0.05$. The discovery of this object implies a spatial number density of $n_0\sim10^{-3}$ au$^{-3}$ for objects with radii greater than or equal to that of 3I/ATLAS. We report deep stacked images obtained using the Canada-France-Hawaii Telescope that display faint activity. Using images obtained from the Las Cumbres Observatory 0.36 m telescopes at Haleakala and the 2.0 m Faulkes Telescope North, we find a small light curve variation of less than 0.2 mag for the object over a $\sim29$ h time span. The visible/near-infrared spectral slope of the object is red, comparable to 1I/`Oumuamua. The object will be observable until September 2025, unobservable near perihelion due to low solar elongation, and observable again in November. This limitation unfortunately prohibits detailed observations at perihelion when the activity of 3I/ATLAS is likely to peak. Based on the experience of 1I/`Oumuamua and 2I/Borisov, we ask the community to maintain a constant observational presence while possible with photometric, spectroscopic, and polarimetric methods. Such observational data would constrain the (i) light curve, (ii) onset and variation of activity, and (iii) nongravitational effects. It is essential that the community collaborate to rapidly and comprehensively characterize these properties of 3I/ATLAS.

Recently, the microquasar V4641 Sgr and its surrounding is detected at TeV-PeV gamma-ray band. Interestingly, the spectrum follows a power-law function continuing up to 0.8 PeV as reported by LHAASO, and the morphology of the emission appears a puzzling elongated structure. In this work, we propose that the elongated UHE emission from V4641 Sgr could originate from the inverse Compton radiation of electrons with a very hard spectrum, which may result from shear acceleration mechanism in the jets driven by V4641 Sgr. We also calculate the corresponding X-ray synchrotron emission from the same electron population, predicting the potential range of non-thermal X-ray flux of the source. The recent observation by XRISM toward the central part of the UHE source could pose a constraint on the model parameters. In the future, a full coverage of the source by sensitive X-ray instrument and high-resolution TeV-PeV gamma-ray instrument may provide a critical test of the model.

Astrophysical flares are one of the possible prominent source classes of ultra-high-energy (UHE, $E > 10^{17}$ eV) cosmic rays, which can be detected by recording clusters of extensive air showers in arrays of detectors. The search for sources of neutral particles offers distinct advantages over searching for sources of charged particles, as the former traverse cosmic distances undeflected by magnetic fields. While no cosmic-ray photons exceeding $10^{17}$ eV have been definitively detected, identifying the clustering of events in cosmic-ray data would provide compelling evidence for their existence. We compare two analysis methods for detecting direction-time clustering in UHE extensive air showers: an approach that examines multiplets, and the stacking method, which analyzes sets of doublets that are not necessarily consecutive, thus making it sensitive to multiple flares. Both techniques combine time-clustering algorithms with unbinned likelihood study. Background events (initiated by hadrons) can be more efficiently distinguished from photon-induced events (signals) by using a photon tag that employs probability distribution functions to classify each event as more likely to be initiated by either a photon or a hadron. We demonstrate that these methods can effectively distinguish between events initiated by photons and those initiated by hadrons (background), and can accurately reproduce both the number of photon events within flares and their duration. We calculate the discovery potentials, i.e., the number of events required to identify a photon flare. The methods discussed can be used to search for cosmic ray sources and/or improve limits on the fluxes of UHE photons.

We present a new framework to incorporate feedback from massive interacting binaries in simulations of star cluster formation. Our new feedback model adds binary stellar evolution to the cluster formation code Torch, and couples it in AMUSE to the pre-existing modules for collisional stellar dynamics, magnetohydrodynamics, and mechanical and radiative feedback. Our model accounts for the effects of mass transfer on the stars' mass loss rates, their radiation spectra, and the timing of core-collapse supernovae. It also injects mass lost through non-conservative mass transfer and common envelope ejection into the interstellar medium. We demonstrate the use of our feedback model through simulations of isolated binaries in a gaseous medium, and of embedded clusters of massive binaries. Feedback from interacting binaries efficiently couples with the surrounding interstellar medium. It increases the size of HII regions, increases the kinetic and thermal energy of the gas, and increases the pressure within HII regions compared to models that use single star stellar evolution. Those differences arise from the ionizing radiation, which increases by three orders of magnitude, resulting in HII regions that expand due to thermal pressure rather than radiation pressure. The effects of stellar dynamics and the gravitational potential of the background gas cause the evolution of individual binaries to deviate from the predictions made by secular evolution, impacting the subsequent feedback from the binary. We conclude that massive interacting binaries are an important source of feedback in cluster-forming regions, and must be considered when studying the emerging timescales of young star clusters.

Ultra-hot Jupiters (UHJs) offer exceptional opportunities for detailed atmospheric characterization via emission spectroscopy. Here we present a comprehensive analysis of the dayside atmosphere of WASP-19b$\unicode{x2014}$one of the shortest-period UHJs$\unicode{x2014}$using archival JWST NIRSpec/PRISM observations spanning 0.6-5.3 $\mu$m. We report robust detections of H$_2$O (16.44 $\sigma$), CO (5.47 $\sigma$), and CO$_2$ (10.72 $\sigma$), along with marginal detections of CH$_4$ (3.76 $\sigma$) and C$_2$H$_2$ (2.45 $\sigma$). The retrieved composition reveals a highly carbon-rich atmosphere with a tightly constrained super-solar C/O ratio of 0.94$\pm$0.03. Elevated abundances of carbon-bearing species provide strong evidence (11.69 $\sigma$) for disequilibrium chemistry. We also detect condensate clouds, likely Al$_2$O$_3$(c), at high significance (17.28 $\sigma$), and constrain the atmospheric metallicity to 1.7$_{-0.7}^{+1.2}$ $\times$ solar. These results establish a precise benchmark for modeling dayside conditions in extreme irradiated atmospheres and demonstrate JWST's transformative capabilities for exoplanet science.

Early multi-wavelength observations offer crucial insights into the nature of the relativistic jets responsible for gamma-ray bursts and their interaction with the surrounding this http URL present data of GRB 240825A from 17 space- and ground-based telescopes/instruments, covering wavelengths from NIR/optical to X-ray and GeV, and spanning from the prompt emission to the afterglow phase triggered by Swift and Fermi. The early afterglow observations were carried out by SVOM/C-GFT, and spectroscopic observations of the afterglow by GTC, VLT, and TNG determined the redshift of the burst ($z = 0.659$) later.A comprehensive analysis of the prompt emission spectrum observed by Swift-BAT and Fermi-GBM/LAT reveals a rare and significant high-energy cutoff at ~76 MeV. Assuming this cutoff is due to $\gamma\gamma$ absorption allows us to place an upper limit on the initial Lorentz factor, $\Gamma_0 < 245$. The optical/NIR and GeV afterglow light curves be described by the standard external shock model, with early-time emission dominated by a reverse shock (RS) and a subsequent transition to forward shock (FS) emission. Our afterglow modelling yields a consistent estimate of the initial Lorentz factor ($\Gamma_{\rm 0} \sim 234$). Furthermore, the RS-to-FS magnetic field ratio ($\mathcal{R}_B \sim 302$) indicates that the reverse shock region is significantly more magnetized than the FS region. An isotropic-equivalent kinetic energy of $E_{\text{k,iso}} = 5.25 \times 10^{54}$ erg is derived, and the corresponding $\gamma$-ray radiation efficiency is estimated to be $\eta_{\gamma}$ = 3.1%. On the other hand, the standard afterglow model can not reproduce the X-ray light curve of GRB 240825A, calling for improved models to characterize all multi-wavelength data.

Large aperture telescopes are pivotal for exploring the universe, yet even with state-of-the-art manufacturing and launch technology, their size is limited to several meters. As we aim to build larger telescopes - extending tens of meters - designs in which the main mirror is based on liquid deployment in space are emerging as promising candidates. However, alongside their enormous potential advantages, liquid-based surfaces present new challenges in material science, mechanics, and fluid dynamics. One of the fundamental questions is whether it is possible for such surfaces to maintain their precise optical shape over long durations, and in particular under the forces induced by the telescope's accelerations. In this paper, we present a model and a closed-form analytical solution for the non-self-adjoint problem of the dynamics of a thin liquid film pinned within a finite circular domain. We use the 50-meter Fluidic Telescope (FLUTE) concept as the case study, and examine the liquid dynamics of the telescope under both slewing actuation and relaxation regimes, elucidating the role of geometrical parameters and liquid properties. The solutions reveal a maneuvering 'budget' wherein the degradation of the mirror surface is directly linked to the choice of maneuvers and their sequence. By simulating ten years of typical operation, we show that, while the maximal deformation might reach several microns, the spatial distribution of the deformation and their propagation rate allows the telescope to maintain its optical functionality for years, with at least a substantial portion of the aperture remaining suitable for astronomical observations. The model provides valuable insights and guidelines into the performance of liquid-film space telescopes, marking a crucial step toward realizing the potential of this innovative concept.

We examine the joint effects of relativistic accretion and memory burdened evaporation on the evolution of primordial black holes (PBHs). The memory burden effect, which delays the evaporation by inducing a backreaction and making the evaporation rate scale as an inverse power law of the PBH entropy, opens up a new window that allows PBHs with $M \lesssim 10^{15}~\mathrm{g}$ to survive until the present epoch. Meanwhile, accretion increases the mass of PBHs, thereby enhancing their chances of survival for a given initial mass. We consider two main scenarios: one where PBHs evaporate completely before big bang nucleosynthesis, and another where PBHs persist until today. In the case of evaporation, we analyse the emission of dark matter (DM) and dark radiation (DR) during the process of evaporation. Conversely, in the other case, the surviving PBHs themselves can contribute as DM. We further investigate how relativistic and non-relativistic accretion, together with memory burdened evaporation, impact the parameter space of the emitted DM, the abundance of stable PBHs as DM, and the contribution of DR to the effective number of relativistic degrees of freedom, $\Delta N_{\mathrm{eff}}$.

The Cosmological Principle posits that the Universe is isotropic on the largest scales. While widely supported, this foundational assumption remains testable. We analyse the angular distribution of over one million quasars from the Gaia-unWISE catalogue using Renyi entropy, a multiscale statistical measure sensitive to higher-order clustering. Dividing the sample into three redshift bins, we find that both the low- and high-redshift distributions are statistically consistent with isotropy. However, at intermediate redshift ($1 \leq z < 2.2$), we detect a statistically significant and scale-dependent anisotropy that persists under stringent masking, suggesting a physical origin. We interpret this as evidence for a transitional epoch in cosmic history, during which large-scale structures such as superclusters became prominent before their growth was gradually damped by the onset of accelerated expansion. These findings position Renyi entropy as a powerful probe of cosmic evolution and highlight the potential thermodynamic links between structure formation, entropy dissipation, and the emergence of large-scale isotropy.

Magnetic holes are structures commonly observed in various space plasma environments throughout the solar system, including the solar wind. These structures are characterized by a localized decrease in magnetic field strength, coincident with an increase in plasma density. Previous observational studies in the solar wind link the presence of Langmuir waves to magnetic holes, suggesting a strong correlation between these phenomena. We develop a model based on magnetic-moment conservation and its violation to explain the excitation of Langmuir waves in magnetic holes. Our model illustrates that magnetic holes induce changes in the electron velocity distribution function that emit electrostatic Langmuir waves due to the bump-on-tail instability. Using data from the Solar Orbiter spacecraft, we provide a comprehensive analysis of this process and test our predictions with observations. The consistency between the model and observations indicates that our proposed process is a viable mechanism for producing Langmuir waves in magnetic holes in the solar wind.

Probing primordial curvature perturbations on small scales, beyond those accessible using cosmic microwave background (CMB) primary anisotropies and Lyman-$\alpha$ forest data, remains a major open challenge. Current constraints on the scalar power spectrum at these scales are either weak or rely heavily on model-dependent assumptions about small-scale structure. In this work, we propose a novel method to probe the small-scale primordial power spectrum using scalar-induced tensor perturbations, which are inevitably sourced by curvature perturbations at second order in cosmological perturbation theory. While induced tensor modes have traditionally been studied in the context of the stochastic gravitational wave background, we highlight a complementary observable: the distinctive pattern of B-mode polarization they imprint on the CMB. We compute the angular spectrum of these B-modes arising from enhanced scalar perturbations and show that the resulting signal can be competitive with inflationary predictions for values of the tensor-to-scalar ratio targeted in upcoming CMB experiments, most notably CMB-Stage 4. We map the region of the scalar power spectrum to which these future B-mode experiments will be sensitive and compare with existing constraints, finding it to exceed current sensitivities at $k \sim O(1-10)\, \mathrm{Mpc}^{-1}$. In addition to providing a new CMB-based probe of the small-scale power spectrum, this work also motivates dedicated B-mode searches at higher multipoles ($\ell \gtrsim 100$).

Among the few ways that allow or could allow us to probe the early Universe from the observation of a flux of primordial particles, there is one possibility which has been little studied: the observation today of high energy neutrinos which could have been emitted shortly after the Big Bang, from the decay or annihilation of early universe relics. We perform a general study of such a possibility. To this end we first emphasise that these neutrinos could display various kinds of sharp spectral features, resulting from the primary energy spectrum at emission, and from how this spectrum is smoothed by redshift and radiative correction effects. Next we determine the ranges of mass (from a fraction of eV all the way to the Planck scale) and lifetime of the source particles along which we do not/we do expect that the sharp spectral feature will be altered by interactions of the neutrinos on their way to the detector, mainly with the cosmic neutrino background or between themselves. We also study the theoretical (i.e. mainly BBN and CMB) and observational constraints which hold on such a possibility. This allows us to delineate the regions of parameter space (mass, lifetime and abundance) that are already excluded, hopeless for future observation or, instead, which could lead to the observation of such neutrinos in the near future.

As astronomy enters an era defined by global telescope networks, petabyte-scale surveys, and powerful computational tools, the longstanding goals of astronomy education, particularly introductory ``ASTRO101'', but equally encompassing both higher and lower level courses, warrant fresh examination. In June 2024, the AstroEdUNC meeting at UNC--Chapel Hill convened 100 astronomers, education researchers, and practitioners to synthesise community perspectives on the purpose, content, and delivery of astronomy education. Beginning with historical vignettes, the meeting's deliberations were organised into six interrelated themes: (1) Context, highlighting astronomy's evolution from classical charting to multi-messenger discovery and its role as a connective thread across STEM and the humanities; (2) Content, exploring how curricula can balance essential concepts with authentic investigations and leverage open-source and AI-augmented resources; (3) Skills, arguing that astronomy should foreground scientific literacy, computational fluency, and communication through genuine data-driven inquiry; (4) Engagement, advocating for active-learning strategies, formative assessment, and culturally inclusive narratives; (5) Beyond the Classroom, emphasising scaffolding, universal-design practices, and K--12/community partnerships; and (6) Astronomy Education Research, outlining priority areas for assessing knowledge, attitudes, and long-term outcomes. We provide concrete recommendations for future astronomy education research development, underscoring the need for approaches to education that are authentic while meeting the learning and life goal needs of the students, a vibrant community of practice and robust researcher--practitioner partnerships to ensure that introductory astronomy is pertinent, applicable and inspiring to a broad student population.

We develop a new class of cosmological attractors which are compatible with the recent ACT results. They are based on two types of fractional Kaehler potentials, K, for a gauge-singlet inflaton phi which reduce, along the inflationary path, to the form N/(1-phi^qM)^p with qM=1, 2 and 0.1< p<10. The combination of these K's with the chaotic potentials phi^n (where n=2, 4) within a non-linear sigma model leads to inflationary observables which are consistent with the current data and largely independent from qM and n. Endowing these K's with a shift symmetry we also offer a supergravity realization of our models introducing two chiral superfields and a monomial superpotential, linear with respect to the inflaton-accompanying field. The attainment of inflation with subplanckian inflaton values and the large values for the tensor-to-scalar ratio, which increases with N, are two additional attractive features of our proposal.

We propose a novel mechanism for realizing slow-roll inflation that is fully consistent with observational data, based on conformal transformations acting exclusively on a complex scalar field -- without coupling to the gravitational sector. These transformations generically produce a plateau in the inflaton potential, as guaranteed by the maximum modulus theorem, thereby naturally satisfying the slow-roll conditions. Our framework utilizes squeezing operations generated by the Virasoro algebra without central extension, as developed in our earlier work. The resulting inflationary potentials depend on the Virasoro mode $n$, the power $m$ of the original potential, and the squeezing parameter $\theta$. We present approximate analytical expressions at leading order for the special case $n=-2$, and perform numerical analyses for both $n=-2$ and other values of $n$. These reveal parameter regimes in which the predicted cosmological observables $(n_{s},r)$ align remarkably well with current CMB measurements.

Space-borne gravitational wave detection will open the observation window in the 0.1 mHz$-$1 Hz bandwidth, playing a crucial role in the development of cosmology and physics. Precise clock synchronization among satellites is essential for the accurate detection of gravitational wave signals. However, the independent clock counting mechanisms of each satellite pose a significant challenge. This work reports the mathematical model of clock asynchrony, which is mainly dominated by the constant term factor and the linear term factor. Moreover, it experimentally verifies the clock asynchronization technique based on a dual-phasemeter system. Through experimentation, the impacts of these two aspects of clock asynchrony were confirmed, and post-processing techniques were employed to reduce these impacts to as low as $\rm 2\pi \times 10^{-6} rad/Hz^{1/2}@ 3mHz$. Specifically, the constant term factor is measured by Time-delay Interferometry Ranging (TDIR), while the linear term factor can be gauged by clock transmission link. This study provides a reference for understanding the clock asynchrony mechanism and processing clock synchronization issues. Additionally, a low additional noise clock synchronization test system is introduced to support such measurements.

General relativity predicts that gravitational waves are described by two polarisation states: the plus $+$ state and cross $\times$ state. However, alternate theories of gravity allow up to six polarisations. We employ the gravitational-wave null stream, a linear combination of three or more detectors where the $+$ and $\times$ signals add to zero, leaving behind noise and potentially gravitational waves in non-standard polarisation states. We develop a Gaussian process model to search for extra polarisations beyond general relativity. Using data from 42 three-detector events from LIGO-Virgo-KAGRA's Third Gravitational-Wave Transient Catalogue, we find no evidence of non-standard polarisations. We set upper limits on the fractional deviation in gravitational-wave strain to be as low as 0.39 at 90% credibility for the event GW190602_175927.

In earlier companion papers, we showed that non-singular primordial black holes (PBHs) could account for all the dark matter (DM) over a significantly wider mass range compared to Schwarzschild PBHs. Those studies, mostly based on phenomenological metrics, are now extended by considering the quantum-corrected space-time recently proposed by Zhang, Lewandowski, Ma and Yang (ZLMY), derived from an effective canonical (loop) quantum gravity approach explicitly enforcing general covariance. Unlike the BHs considered earlier, ZLMY BHs are free from Cauchy horizons, and are hotter than their Schwarzschild counterparts. We show that this higher temperature boosts the evaporation spectra of ZLMY PBHs, tightening limits on their abundance relative to Schwarzschild PBHs and shrinking the asteroid mass window where they can constitute all the DM, a result which reverses the earlier trend, but rests on firmer theoretical ground. While stressing the potential key role of quantum gravity effects in addressing the singularity and DM problems, our study shows that working within a consistent theoretical framework can strongly affect observational predictions.

Low-temperature detectors are a powerful technology for dark matter search, offering excellent energy resolution and low energy thresholds. COSINUS is the only experiment that combines scintillating sodium iodide (NaI) crystals with an additional phonon readout at cryogenic temperatures, using superconducting sensors (remoTES), alongside the conventional scintillation light signal. Via the simultaneous phonon and scintillation light detection, a unique event-by-event particle identification is enabled. This dual-channel approach allows for a model-independent cross-check of the long-standing DAMA/LIBRA signal with a moderate exposure of a few hundred kg d, while completely avoiding key systematic uncertainties inherent to scintillation-only NaI-based searches. COSINUS built and commissioned a dedicated low-background cryogenic facility at the LNGS underground laboratories. Data taking with eight NaI detector modules (COSINUS1$\pi$ Run1) is planned to begin in late 2025.

The momentum-dependent interaction (MDI) model, which has been widely used in microscopic transport models for heavy-ion collisions (HICs), is extended to include three different momentum-dependent terms and three zero-range density-dependent terms, dubbed as MDI3Y model. Compared to the MDI model, the single-nucleon potential in the MDI3Y model exhibits more flexible momentum-dependent behaviors. Furthermore, the inclusion of three zero-range density-dependent interactions follows the idea of Fermi momentum expansion, allowing more flexible variation for the largely uncertain high-density behaviors of nuclear matter equation of state (EOS), especially the symmetry energy. Moreover, we also obtain the corresponding Skyrme-like energy density functional through density matrix expansion of the finite-range exchange interactions. Based on the MDI3Y model, we construct four interactions with the same symmetry energy slope parameter $L=35$ MeV but different momentum dependence of $U_{\mathrm{sym}}$, by fitting the empirical nucleon optical potential, the empirical properties of symmetric nuclear matter, the microscopic calculations of pure neutron matter EOS and the astrophysical constraints on neutron stars. In addition, two interactions with $L=55$ and $75$ MeV are also constructed for comparison. Using these MDI3Y interactions, we study the properties of nuclear matter and neutron stars. These MDI3Y interactions, especially those with non-monotonic momentum dependence of $U_{\mathrm{sym}}$, will be potentially useful in transport model analyses of HICs data to extract nuclear matter EOS and the isospin splitting of nucleon effective masses.

Astrophysical and cosmological observations suggest the existence of beyond standard model ingredient known as dark matter (DM). One of the most supported class of theories suggests that DM is composed of weakly interactive massive particles (WIMPs), possibly detectable via weak interaction with standard matter resulting in the recoil of the latter. The motion of the Sun and Earth with respect to the Galactic Centre is expected to induce a strong directional dependence in the recoil spectrum. Direct detection experiments capable of measuring the angular features of the recoils gain access to a wide range of advantages such as the possibility to positively claim a discovery of DM. The CYGNO project sets into this context, with the aim of deploying a large directional detector for rare event searches as DM. It exploits a gaseous time projection chamber filled with a He:CF4 gas mixture with a segmented amplification stage and granular optical readout. In this thesis, it is presented the work carried out with small CYGNO prototypes to maximise the light yield without degrading spatial and energy resolution, the addition of highly electronegative gases to induce a reduction of the electron cloud diffusion while drifting towards the amplification stage, very relevant to precisely measure the topological information of the recoil tracks. Moreover, the potential performances of directional detectors in the context of a direct DM search are analysed with the use of rigorous statistical tools both in the improvement in setting limits in the WIMP to nucleon sensitivity and in the capability of discerning two different DM models exploiting a directional detector.

The local primordial density fluctuations caused by quantum vacuum fluctuations during inflation grow into stars and galaxies in the late universe and, if they are large enough, also produce primordial black holes. We study the formation of the primordial black holes in $k$-essence inflation models with a potential characterized by an inflection point. The background and perturbation equations are integrated numerically for two specific models. Using the critical collapse and peaks formalism, we calculate the abundance of primordial black holes today.

Most theories that attempt to reproduce the Modified Newtonian Dynamics (MOND) phenomenology for dark matter at galactic scales rely on ad hoc free functions, preventing them from being regarded as fundamental. In this work, we present a new theory that reproduces MOND, built on a supposed to be fundamental Yang-Mills gauge field based on SU(2), with a gravitational coupling constant, and emerging in a low-acceleration regime, below the MOND acceleration scale. The gauge field plays the role of the internal force in the dipolar dark matter (DDM) model. We discuss how certain solutions of this theory recover the deep MOND regime without introducing arbitrary functions in the action. Within this framework, the MOND phenomenology appears to be due to the existence of a new sector of particle physics.

A fast-rolling axion can transfer its kinetic energy to a gauge field via the Chern-Simons coupling, leading to copious production of gauge quanta, which can act as a source of gravitational waves (GWs) with potentially observable amplitudes. In this work, we investigate GW production in a spectator axion model when strong backreaction is taken into account. We find that decreasing the decay constant of the axion enhances GW production. Since the initial value of the axion is larger than its quantum fluctuations, such a condition imposes a lower bound on the axion dacay constant, which sets an upper bound on the amplitude of the energy spectrum of GWs. As a result, the amplitude of the predicted GW energy spectrum is lower than $10^{-10}$ in the nHz to mHz frequency range.

We classify $f(R)$ theories using a mathematical analogy between slow-roll inflation and the renormalization-group flow. We derive the power spectra and spectral indices class by class and compare them with the latest data. The framework used for the classification allows us to determine the general structure of the $f(R)$ functions that belong to each class. Our main result is that only two classes survive. Moreover, we show that the Starobinsky model is the only polynomial $f(R)$ that can realize slow-roll inflation. In fact, all other polynomials belong to a special class that can only realize constant-roll inflation, at least far enough in the past. We point out some of the issues involved in considering a smooth transition between constant-roll and slow-roll inflation in this class of models. Finally, we derive the map that transforms the results from the Jordan frame to the Einstein frame.

Holography and entropy bounds suggest that the ultraviolet (UV) and infrared (IR) cutoffs of gravitational effective theories are related to one another as a form of UV/IR mixing. Motivated by this, we derive a bound on the allowed scalar field range in theories with cosmic horizons. We show how this bound challenges several inflationary scenarios, such as $\alpha$-attractors and modular-invariant inflation. Besides, we find a relation between the number of extra spatial dimensions and the tensor-to-scalar ratio.