Papers

Ultra-fast outflows (UFOs), launched from the inner accretion disc in active galactic nuclei (AGN), are among the prime candidates for driving AGN feedback in their host galaxies. Determining their physical properties and energetics is therefore essential for assessing their impact on galaxy evolution. This review highlights recent observational progress, with particular emphasis on the groundbreaking insights enabled by high-resolution X-ray spectroscopy. The first X-Ray Imaging and Spectroscopy Mission (XRISM)/Resolve observations of AGN winds have revealed that the innermost UFOs are kinematically structured and inhomogeneous. In the archetypal quasar PDS 456, these observations have led to a substantial upward revision of the inferred mass outflow rate, momentum flux and kinetic power. Complementary UV, optical, IR and mm observations trace the ionised, atomic and molecular phases of AGN-driven outflows from parsec to galactic scales. However, the physical coupling between the accretion-disc wind and the galactic outflows remains one of the major unresolved questions in AGN feedback. Building on the pioneering XRISM results, NewAthena's X-ray Integral Field Unit (X-IFU) will characterise the physical properties and energetics of accretion-disc winds with unprecedented precision from the local Universe up to cosmic noon, while complementary multi-wavelength observations with facilities such as the Extremely Large Telescope (ELT), James Webb Space Telescope (JWST), Atacama Large Millimeter/submillimeter Array (ALMA), and Square Kilometre Array (SKA) will determine whether and how their energy and momentum are transferred to the host galaxy.

While galaxies generally live in clusters and groups, particularly small groups, or other structures such as walls or filaments, some galaxies appear relatively `isolated', i.e. lacking in close neighbours. Here we use the GAMA survey, and an objective measure of isolation, to explore the most isolated galaxies. We then determine whether these have particular physical characteristics. We find that, as expected, while grouped galaxies split between blue, disk-like and star forming and red, spheroidal and quiescent, isolated galaxies are preferentially the former. However, there is a broad overlap in properties, with some isolated galaxies sharing the latter characteristics. Considering only the star-forming galaxies, when controlled by mass we find no differences between the grouped and isolated galaxies in terms of their star formation and stellar populations. Given the lack of significant interactions for the isolated sample, this implies that the evolution of present-day star-forming galaxies in all environments (at least outside rich clusters) has been, and remains, dominated by internal processes.

Primordial black holes (PBHs) form from rare fluctuations in the high-density tail of the primordial distribution, making their abundance highly sensitive to primordial non-Gaussianity. Heavy fields during inflation can source such non-Gaussianity, producing characteristic signatures widely explored in cosmological collider physics, yet its impact on PBH formation remains largely unexamined. We investigate this question within quasi-single-field inflation (QSFI), the minimal realization that couples a heavy field to the inflaton. Unlike most PBH studies that consider only the bispectrum, the QSFI can generate a sizable trispectrum alongside the bispectrum, allowing us to investigate their combined effects on PBH formation. In our leading-order calculation, heavy-field non-Gaussianities can amplify the high-density tail and substantially enhance PBH production. Within the parameter scan, the largest enhancement occurs at the smallest heavy-field mass and the largest sound-speed ratio considered. The trispectrum contribution can be comparable to that of the bispectrum and may even dominate the total non-Gaussian correction. Our results establish PBH abundance as a novel observational channel for cosmological collider physics, demonstrating how heavy fields during inflation can significantly shape PBH formation.

Satellite galaxies in massive halos are seen to be preferentially quenched along the major axes of their central galaxies. In contrast, the AGN jet tends to point along the minor axis, which imposes interesting constraints on the role of AGN feedback in satellite galaxy quenching. Here we use the SIMBA, TNG100, EAGLE, and feedback-free SIMBA cosmological simulations at $z=0$ to show that this satellite anisotropy primarily originates from spatial segregation between quenched and unquenched populations. Quenched satellites are more centrally concentrated and more strongly aligned with the central galaxy's major axis, whereas unquenched satellites occupy larger radii and are more nearly isotropic. The persistence of this behaviour without stellar or active galactic nucleus feedback suggests that feedback is not important for satellite quenching. Beyond the halo boundary, the signal varies strongly among simulations and correlates with the excess number of galaxies along the major axis. This links large-scale anisotropic quenching to the direction-dependent distribution of groups, clusters and filamentary environments. Our results provide a unified geometrical interpretation of anisotropic satellite quenching across halo and extra-halo scales.

Measurements of the magnetic field strength in the solar corona at heliocentric distances $r\,{\gtrsim}\,1.1\,R_{\odot}$, using observations other than transient radio bursts, are rare. Circular polarization observations of thermal bremsstrahlung radio emission at frequencies $<$ 100 MHz, which typically originate at $r\,{>}\,1.1R_{\odot}$ in the solar corona, can be effectively used for this purpose. The thermal bremsstrahlung emission is randomly polarized, but in the presence of a magnetic field it propagates in two oppositely polarized circular modes whose differential absorption in the medium results in an observable circular polarization. We report a case study of two-dimensional imaging observations of the above mentioned circular polarization in the solar corona at 51\,MHz using the Gauribidanur RAdioheliograPH (GRAPH), and estimates of the associated coronal magnetic field strength. The peak degree of circular polarization ($dcp$), calculated from the observed total and circularly polarized intensity images with the GRAPH, is ${\approx}$2.7%. Using ray-tracing simulations we find that the observed peak $dcp$ corresponds to a radio source located at $r\,{\approx}\,1.65\,{\pm}\,0.15\,R_{\odot}$ at 51 MHz, and its field strength is ${\approx}\,1.55\,{\pm}\,0.3\,$Gauss.

It is now widely accepted that a supermassive black hole (SMBH) resides in the core of nearly every galaxy. Within the context of mergers of galaxies, the central SMBHs from each system should form a gravitationally bound binary system. Given the frequency of galaxy mergers throughout cosmic time, the existence of such binary SMBHs is likely to occur within active galactic nuclei. However, identifying them observationally remains highly challenging. The majority of reported candidates rely on indirect evidence rather than unambiguous detection. The quasar Markarian 231 (Mrk 231) stands out because of its distinctive ultraviolet to optical properties, which includes a strong wavelength dependent polarization and an unusual flux decrease at 2500A. In this work, we investigate the potential of spectropolarimetry as a powerful tool for identifying and characterizing SMBH binaries in AGNs. We have developed a simplified analytical model to reproduce the spectral energy distribution of the quasar Mrk 231 assuming a binary SMBH geometry, and we compute both its total and polarized emission. The main aim is to test whether two separate emitting sources could explain the unique UV to optical signatures of Mrk 231. Such investigations highlight the importance of a future high-sensitivity UV spectropolarimetric instrument. For instance, the proposed Pollux spectropolarimeter for the Habitable Worlds Observatory (HWO) has the potential to significantly contribute to our understanding of AGNs structure and the role of SMBH binaries in quasars.

The LUX-ZEPLIN (LZ) collaboration recently reported a singular anomalous nuclear recoil event at $248 \pm 23 \, \rm keV$, accompanied by a strict null result for unexplained excesses in their primary $5.4-50 \, \rm keV$ search window. The macroscopic momentum transfer required to generate this event is kinematically inaccessible to standard halo cold dark matter (CDM), naturally motivating models featuring light, relativistic relics. In this work, we investigate whether a contemporary flux of Hawking-boosted dark matter (DM) emitted by evaporating Primordial Black Holes (PBHs) in the $10^{11}-10^{13} \, \rm g $ mass range can source the anomaly. By parameterizing the nuclear scattering with heavy non-relativistic effective field theory (NREFT) operators and endothermic inelastic mass transitions, we formulate a strict kinematic exclusion. We demonstrate that the inherent thermal nature of the Hawking emission, coupled with the relativistic kinematics of the incident flux, inevitably overproduces low-energy recoils, aggressively violating the LZ background bounds. Utilizing this failure, we map the LZ low-energy null results into novel, highly stringent upper limits on the PBH abundance fraction, excluding $f_{\rm PBH} \gtrsim 10^{-6}$ for $M_{\rm PBH} \sim 2 \times 10^{11} \, \rm g $ at $ 90 \% $ confidence level.

Nearly all galaxies are believed to host a supermassive black hole at their center. In the context of galaxy mergers, the SMBHs from the progenitor galaxies are expected to form a bound binary system. Given the high frequency of such mergers throughout cosmic time, binary SMBHs should be relatively common in active galactic nuclei (AGNs). Yet, their direct detection remains a major observational challenge, with most candidates inferred only through indirect signatures. The quasar Markarian 231 presents unique ultraviolet-optical properties, including an atypical drop in the total flux around 2500 A and a strong wavelength-dependent polarization. In this work, we explore the potential of spectropolarimetry as a tool to detect and characterize SMBH binaries in AGNs. We construct a simple analytical model in order to reproduce the spectral energy distribution (SED) of Mrk 231 taking into account a binary SMBHs scenario and we compute the resulting total and polarized flux. Our goal is to test if the presence of two distinct sources can reproduce the unique ultraviolet-optical properties observed in Markarian 231. Such investigations support the case for future high-sensitivity UV spectropolarimeters in order to probe the inner structure of quasars. In particular, the next-generation instrument POLLUX, proposed for the Habitable Worlds Observatory, would be particularly suited to provide new insights into AGNs physics and supermassive black hole binaries in quasars.

The radial velocity (RV) technique enabled the discovery of the first exoplanets orbiting Sun-like stars. Since then, improvements in precision and sensitivity have enabled the detection of smaller planets, bringing us closer to the goal of discovering Earth analogs in the coming years. In this chapter, we provide an overview of the radial velocity method. We begin with a detailed description of the technique, including how precise radial velocities are measured and what limitations must be addressed to achieve reliable detection of small signals. We examine the unique challenges associated with the detection of circumbinary planets and outline the developments that have, for the first time, made their detection feasible using the radial velocity method. With over 2,000 planets detected to date, the radial velocity method remains one of the most reliable techniques, offering valuable insights into the demographics of exoplanets, especially at long orbital distances. As an example of its applications, we provide an analysis and interpretation of three publicly available radial velocity datasets: one from a single star, one for a circumbinary planet system, and another from a single star with transiting planets. Furthermore, we compile information on the spectrographs currently and soon-to-be available worldwide for radial velocity observations and a list of open-source tools designed for their analysis. In the Supplementary Material, we include a sample telescope proposal for obtaining radial velocity data, along with practical tips on writing a competitive proposal to aid early-career researchers entering the field. This chapter serves as a comprehensive guide for understanding the capabilities and applications of the radial velocity method, equipping readers with the knowledge and resources to address their own scientific questions using this technique.

This investigation aims to test the validity of general relativity (GR) on kiloparsec scales via galaxy-scale strong gravitational lensing (SGL) combined with multiple cosmological probes. To circumvent the circularity problem induced by the presumption of a cosmological model based on GR, we apply the distance sum rule in the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric to constrain the post-Newtonian parameter gamma PPN and spatial curvature Omega k independently of any cosmological model. The study introduces two methodological improvements: (i) parameterizing the dark-energy evolution factor with second-order Chebyshev orthogonal polynomials, which ensures uniformly high convergence across the full redshift baseline while retaining Omega k as an explicit free parameter; (ii) extending the original two-probe configuration (SGL + Type Ia supernovae) to a four-probe joint analysis that includes baryon acoustic oscillations and cosmic chronometers. Our results show that gamma PPN = 1.119 (+0.072/-0.072) and Omega k = 0.101 (+0.077/-0.076) at 68% confidence. The gamma PPN value is consistent with the GR prediction of unity at 1.7 sigma, while Omega k favors a mildly open universe, with flatness consistent at 2 sigma. Taking model-independent constraints as a benchmark, we quantify the systematic shift in gamma PPN induced by eight representative dark energy models, and find that all offsets are negligible, well below the baseline 1 sigma uncertainty, with no significant bias at current precision. These results confirm that conventional strong-lensing gravity tests introduce no measurable systematic bias, validating our model-independent framework.

We present JWST NIRSpec and MIRI observations of the metal-poor T dwarf WISEA J155349.96+693355.2. The combined NIRSpec/prism plus MIRI/LRS spectrum (R ~ 100) provides 0.6-12 $\mu$m spectroscopic coverage that reveals the presence of H$_2$O, CH$_4$, and NH$_3$ absorption features, enhanced collision-induced H$_2$ absorption, and weak or absent CO and CO$_2$ absorption, all indicators of a metal-poor low-temperature atmosphere. Comparison of the spectral data to a suite of atmosphere models yields an effective temperature $T_{\rm eff}$ = 960$^{+28}_{-32}$ K and a solar-scaled metallicity of [M/H] = -1.16$^{+0.28}_{-0.25}$ dex. We find indications that the uncertain trigonometric parallax, likely subject to Lutz-Kelker bias, may be overestimated. Adopting the photometric distance yields a bolometric luminosity log$_{10}L_{bol}/L_{\odot}$ = -5.29 $\pm$ 0.10, from which we infer $T_{\rm eff}$ = 974 $\pm$ 56 K, consistent with the temperature derived via our atmospheric model fits. From the moderate-resolution (R ~ 3,000) NIRSpec/G395H spectrum, we measure a radial velocity of -163 $\pm$ 5 km/s, which combined with the source's measured proper motion and estimated distance indicates a highly eccentric (e $\geq$ 0.9) Galactic orbit consistent with an accreted halo object. The UVW kinematics of WISEA J155349.96+693355.2, its bulk metallicity, and its modest alpha enrichment ([$\alpha$/Fe] = +0.15$^{+0.04}_{-0.06}$ dex) are all equivalent to stars identified in the Gaia-Enceladus Milky Way substructure, making it the first brown dwarf with kinematic and chemical evidence of association.

Monte Carlo sampling of strong-field quantum electrodynamics processes underpins simulations of high-intensity laser experiments and of astrophysical compact-object magnetospheres. Sampling an event requires the total rate of the process together with the cumulative probability that determines how energy is partitioned between the produced particles. Simulations typically tabulate both in advance and invert the tabulated probability numerically. Here we replace this procedure with elementary-function approximations for synchrotron radiation and the nonlinear Breit–Wheeler process. For each process, we approximate the auxiliary function that sets the total rate, as well as the cumulative probability, with Padé approximants chosen so that the inversion reduces to a quartic equation. This yields the sampled quantum parameter—electron $\chi_e$ or photon $\chi_\gamma$—in closed form. The approximations and the particle spectra sampled from them agree with the exact results to within $1\%$. The procedure requires no lookup tables, no interpolation, and no numerical root finding, and can be inserted directly into radiative particle-in-cell codes.

We present a comparative broadband spectro-temporal study of the ultra-compact X-ray binary (UCXB) candidates 1A 1246$-$588 and 4U 0614$+$091 using AstroSat observations. Spectral modelling shows that the spectrum of 1A 1246$-$588 is well described by thermal Comptonization and a soft thermal disk component. We identify a candidate narrow quasi-periodic oscillation (QPO) at $993\pm4$ Hz with a single-detector significance of $ ~ 3.7\sigma$, tentatively consistent with the upper kHz QPO previously reported from this source. An independent cross-instrument check does not confirm this feature, and we therefore treat it as a tentative, unconfirmed signature rather than a secure detection. In contrast, across the 4 observations, 4U 0614$+$091 evolves from a Comptonization-dominated hard state to a thermally dominated soft state, with the electron temperature decreasing from ~ 19 keV to ~2.3 keV, while the disk and blackbody temperatures increase. This evolution is accompanied by an increase in the disk contribution to the total flux from ~ 18% to ~ 50% and the appearance of an Fe emission line near 6.97 keV in the soft state. The inner disk radius remains approximately constant at ~13- 15 $R_g$, suggesting that the observed evolution is driven primarily by changes in the thermal properties of the accretion flow rather than substantial variations in disk truncation. When the source is in its hardest state, it displays the richest timing behaviour, with variability components spanning from ~ 12 Hz to nearly 1 kHz, whereas the softer observations are dominated by low frequency variability, along with a persistent component near ~140 Hz and a prominent ~ 100 Hz feature in O4. These results suggest that progressive cooling of the Comptonizing corona, along with enhanced thermal emission from the disk and boundary layer, play a key role in shaping the spectral and timing evolution of 4U 0614$+$091.

In this work, we present the General Interface for MPI-AMRVAC / Legolas Interconnection (or GIMLI, for short) as a toolkit for in-depth (magneto)hydrodynamic studies combining linear and non-linear methods. Acting as a bridge between the linear (M)HD spectroscopy code Legolas and the non-linear simulation code MPI-AMRVAC, GIMLI facilitates the initialization of non-linear MPI-AMRVAC simulations with linear solutions. We showcase this functionality on both Cartesian and cylindrical setups, covering representative plasma physical and astrophysical scenarios: Parker unstable magnetized atmospheres, interchange modes in a straight tokamak setup, resistive Harris sheet evolutions and thermally modified discrete Alfvén waves in coronal loops. We offer suggestions on how to combine linear and non-linear methods, and demonstrate how to leverage GIMLI in the process.

Gravitational-wave observations have revealed hundreds of stellar-mass black holes, yet only about two dozen are known in the Milky Way, almost all in binaries. We present POPKIN, a Python framework that couples single- and binary-star evolution with Galactic orbital dynamics to trace black-hole progenitors from the zero-age main sequence to the present-day isolated black-hole (IBH) population. Across ten models varying the supernova (SN) prescription, mass-transfer efficiency, and common-envelope ejection efficiency, the total IBH abundance is controlled primarily by the SN prescription. Our fiducial model, with a recently proposed metallicity- and stripping-history-dependent SN prescription, predicts $\sim4\times10^7$ IBHs in the Galaxy, including $\sim8\times10^4$ within $1\,\rm{kpc}$ of the Sun; alternative SN prescriptions predict $\sim(1-2)\times10^8$ IBHs. The fiducial model yields a bimodal mass distribution, peaking near $9\,M_{\odot}$ and $20\,M_{\odot}$, with a deficit at $13-17\,M_{\odot}$. This distinguishes it from alternative prescriptions, some of which populate the $2-5\,M_{\odot}$ mass-gap region. Non-kicked IBHs follow nearly circular orbits near the Galactic plane, with typical peculiar velocities of $20-30\,\rm{km\,s^{-1}}$, whereas kicked systems undergo stronger radial migration and span a broader velocity range. We estimate $\sim5\times10^3$ accreting IBHs with $F_{\rm X}>10^{-14}\,\rm{erg\,s^{-1}\,cm^{-2}}$, nearly all non-kicked; this estimate is sensitive to the adopted radiative-efficiency and hot-flow treatments. For a Roman-like bulge survey, our fiducial model predicts $\sim360$ intrinsic IBH microlensing events over five years in a $1.70\,\rm{deg^2}$ effective area, before survey-selection effects. We propose that long-timescale microlensing events from IBHs can strongly constrain the SN physics governing stellar-mass black hole formation.

We present a first background-level observational test of a thermogravity theory in which the Einstein equations are clipped to a trace-free version, the cosmological constant $\Lambda$ does not gravitate, and controlled violations of energy conservation obstruct the usual reinstatement of $\Lambda$ as an integration constant. We therefore set $\Lambda=0$ and ask whether late-time acceleration can instead be generated solely by energy non-conservation, with no extra background parameter relative to flat $\Lambda\mathrm{CDM}$. We contrast a minimal universal implementation with a model in which only CDM partakes in non-conservation. The universal model is strongly disfavored when combining Supernovae and BAO distances, with the inclusion of DESI BAO worsening the fit by $\Delta\chi^2\simeq 51$ relative to $\Lambda\mathrm{CDM}$, because matter creation ties the intermediate-redshift normalization of $H(z)$ too rigidly to the present acceleration. This conclusion, however, should be interpreted with caution, since universal non-conservation would modify the observational dictionary itself. By contrast, restricting non-conservation to CDM leaves baryonic and photon observables unaffected at the background level and therefore allows a self-contained analysis. The resulting $\xi\mathrm{CDM}$ model, which can be thought of as a one parameter extension of $\Lambda\mathrm{CDM}$, provides an improved fit relative to $\Lambda\mathrm{CDM}$ for the background dataset combinations considered, with improvements reaching $\Delta\chi^2=-4.85$ and a maximum Bayesian preference of $\ln\mathcal{B}=2.77$. However, when the BAO ruler is calibrated using BBN or CMB information, the reduced effective early-time CDM density increases the sound horizon and drives $H_0$ towards lower values, thereby increasing rather than alleviating the tension with the distance-ladder calibration.

Planet engulfment leads to planetary mass accretion and angular momentum deposition into the stellar envelope. This study uses MESA models to investigate the long-term impact of engulfment on stellar metallicity and rotational velocity, as well as the effect of rotational mixing on metal enrichment. Our models show that engulfment occurring during the late main-sequence (MS) stage yields significantly higher surface enrichment than at the ZAMS. Rotational mixing generally suppresses enrichment; however, metal-poor stars can sustain super-solar abundances for several Gyr, whereas high-metallicity stars remain largely unaffected. Although engulfment causes an instantaneous spin-up, this angular momentum injection is transient, and the star quickly re-converges to its standard rotational evolution. Thus, MS spin-up cannot serve as a long-term observational indicator. Because convective envelope thickness and MS lifetime dictate the magnitude and duration of enrichment, identifying engulfment signatures requires balancing the long-lived but diluted signals in 0.7 M_sun stars against the strong yet short-lived enrichment in 1.2 M_sun stars. Consequently, 1.0 M_sun stars offering both substantial enrichment and a relatively long MS lifetime represent optimal targets for observational searches.

Some active galactic nuclei (AGNs) are known to show extremely asymmetric broad Balmer lines, with their main peak redshifted or blueshifted by thousands of km s$^{-1}$, severely challenging our understanding of the standard AGN paradigm. We wanted to explore the causes of such asymmetric features by carrying out a spectropolarimetric study of a sample of bright AGNs with known extreme Balmer line profiles. We present optical spectropolarimetry of seven bright (V $<$ 17.5) Seyfert-1 galaxies obtained with the Very Large Telescope (VLT) FOcal Reducer/low dispersion Spectrograph 2 (FORS2) instrument from March 2012 to January 2013. The main broad Balmer lines (H$\alpha$ and H$\beta$, and H$\gamma$ occasionally) were sampled with a resolving power ranging from 350 to 500. We performed a detailed spectral decomposition of each H$\alpha$ profile to isolate kinematic substructures within the broad line region (BLR) and measured the polarization degree and angle of each component, providing constraints on the inner geometry and scattering properties. The broad H$\alpha$ components exhibit various scattering geometries. Their polarization as a function of full width at half maximum (FWHM) reveals two distinct behaviors, corresponding to either a BLR of roughly constant scale height or a stratified structure with decreasing height at higher velocities. This supports a disk-wind geometry linking equatorial and polar regions.

The acceleration of the Solar System barycenter (SSB) offers a unique precision probe of the local gravitational environment, enabling searches for otherwise invisible gravitating sources such as Planet Nine and nearby (primordial) black holes. Using the binary-pulsar timing measurements, we develop an analysis framework that combines the line-of-sight differential accelerations of 26 binary pulsars to jointly reconstruct the three-dimensional acceleration of the SSB and place directional upper limits. Across three smooth Galactic-potential baselines, the residual acceleration is consistent with zero. We obtain directional 95% upper limits of $0.198$, $0.262$, and $0.693~\mu\mathrm{as}\,\mathrm{yr}^{-1}$ over 50%, 75%, and all sampled directions, respectively. These limits are tighter by factors of 3.5, 3.6, and 1.6 than a matched Gaia EDR3 benchmark and improve previous pulsar constraints by more than an order of magnitude. We further extend the framework to the full SSB–pulsar two-endpoint response, enabling direct position-dependent mass constraints on unknown gravitational sources. A $10\,M_\odot$ object is excluded within 0.39, 0.34, and 0.21 pc over the same sky fractions. Our unknown-mass limits surpass the tidal-equivalent INPOP19a reference beyond 0.2 pc by about an order of magnitude over most of the sky at 1 pc.

Type IIP supernovae (SNe IIP) are the most common subclass of core-collapse SNe in observations. However, SNe IIP with short plateaus of the order of tens of days are rarely observed. The progenitors for this kind of SN can help to address the red supergiant issue in stellar evolution. In this article, we report optical photometry and spectroscopy of SN\,2025aedz, a rapidly post-peak declining SN IIP with a typical short plateau. It exhibits a peak absolute magnitude of $M_r=-17.16\pm0.03$\,mag. The $r$-band light curve shows a steep early post-peak decline of $\sim5\,\mathrm{mag}\,(100\,\mathrm{d})^{-1}$ followed by a relatively short plateau, with a plateau duration of $\sim50\pm3$\,d. The overall spectral evolution is consistent with that of normal SNe~IIP, showing a blue continuum with prominent Balmer P-Cygni profiles during the photospheric phase, followed by the gradual strengthening of hydrogen and metal lines as the ejecta cools down, although the metal lines remain weak and the expansion velocities decline rapidly. SN\,2025aedz is similar to the short-plateau SN\,2018gj in its overall evolution, whereas its pronounced early decline resembles that of SN\,2023ufx, which has the shortest plateau duration known so far. The radioactive tail of the bolometric light curve implies a synthesized $^{56}$Ni mass of $\sim0.03\pm0.01\,M_\odot$. Motivated by the steep early decline, we performed radiation hydrodynamic simulations by exploring different circumstellar material configurations to reproduce its early bolometric light curve. These simulations indicate that SN\,2025aedz originated from a progenitor with a relatively low-mass hydrogen envelope, possibly produced through enhanced mass loss or binary interaction, while the steep early decline is likely explained by additional luminosity from circumstellar interaction.

The centers of galaxies harbor the densest stellar environments, where a massive black hole (MBH) accelerates stars to such high velocities that direct collisions can result in high-energetic phenomena, such as gravitational wave sources, kilonovae, and supernova-like transients. These collisions can reshape the cluster's density profile and release gas that can be subsequently accreted by the MBH. However, the evolving rates of such phenomena from self-consistent dynamics around mass-growing MBHs remain largely unexplored. In this work, we simulate nuclear star clusters (NSCs) across a range of masses and density profiles by employing the GNC Monte Carlo code, that self-consistently models stellar dynamics and the subsequent accretion of released gas. We find that stellar collisions flatten the density cusp in the innermost regions ($r \lesssim 10^{-3}-10^{-2}$ pc) within $\sim 0.1-1$ Gyr. While high initial collision rates in steep cusps quickly decline due to stellar depletion, the interplay between collisions and MBH growth is important only in massive NSCs ($M_\star \sim 10^9 M_\odot$). As MBH grows, increased stellar velocities shift the balance toward destructive collisions of which relative velocities can be $\gtrsim 2500 {\rm \,km\,s^{-1}}$. Consequently, present-day destructive collision rates in massive clusters remain high ($10^{-4}\sim 10^{-3}{\rm yr}^{-1}$), whereas they are smaller in Milky Way-like NSCs or negligible in smaller NSCs. Our results highlight a crucial synergy between stellar dynamics and MBH growth, identifying massive galaxies as prime targets for observing transients from destructive stellar collisions.

Modeling the full inspiral-merger-ringdown signal requires combining analytically controlled inspiral physics with the nonlinear strong-field information supplied by numerical relativity. We represent the numerical relativity contribution beyond an analytic inspiral waveform as residual amplitude and phase corrections. Retaining the leading-order frequency-domain amplitude and the $3.5$ Post-Newtonian TaylorF2 phase, we use a Kolmogorov-Arnold network to learn these residual corrections from SXS waveforms. After training, the learned corrections are stored as explicit spline functions, so waveform evaluation no longer requires the network itself. On $75$ simulations excluded from training and model selection, the model achieves a median flat-noise mismatch of $2.7\times10^{-5}$. Our results demonstrate a machine learning driven waveform modeling strategy in which numerical relativity augments, rather than replaces, analytically known waveform structure.

Recent stellar-kinematic and neutral-hydrogen observations indicate a dark matter deficit within the central tens of kiloparsecs of nearby galaxies that grows systematically with stellar mass relative to hydrodynamical simulations. We show that this mass dependence can arise from exothermic inelastic dark matter with strongly velocity-dependent scattering. Two nearly degenerate dark matter states interact through vector and scalar mediators with opposite-sign contributions and unequal ranges. The resulting coupled-channel dynamics suppresses $s$-wave conversion at low velocity while retaining a $p$-wave enhancement at several hundred $\rm km/s$, where down-scattering injects kinetic energy and lowers central dark matter densities. At dwarf velocities, conversion remains suppressed while elastic scattering can still drive core formation consistent with the observed dwarf-clustering pattern. A late dark-sector phase transition, along with the low-velocity suppression, preserves a large excited-state fraction until halo formation. Using representative halos spanning the four observed stellar-mass bins, we find that the model accounts for the inferred dark matter deficit in all four bins, including its systematic growth with stellar mass, while its impact weakens toward cluster velocities. The growing discrepancy with hydrodynamical simulations may therefore offer a glimpse of previously hidden dark matter microphysics.

The D-region ionosphere is sluggish in nature while responding to the incoming solar radiation. We study the altitude (h) profile of mid-latitude D-region response time delay (Delta t) during three chosen solar flares, namely, C5.2, M5.2, and X2.2 classes. We solve "electron continuity equation" numerically at each and every h to obtain Delta t-h profile. We investigate the 1) latitudinal variation (over both Northern and Southern hemispheres) and 2) seasonal variation (throughout the year) of Delta t-h profiles of each of these solar flares separately. We go over noteworthy variations of Delta t-h profile for both 1) and 2) with reasonably different results over different hemispheres. Also, we study some contrasting behaviours of Delta t-h profile of X2.2 class for both 1) and 2) in comparison to C5.2 or M5.2 classes. We conclude the outcome with possible explanations.

Pulsar timing array (PTA) collaborations have reported evidence for a stochastic gravitational wave (GW) background in the nHz band. Should scalar-induced GWs, sourced at second order by enhanced primordial curvature perturbations, contribute to this signal, a coincidence of scales makes them directly sensitive to the softening of the equation of state around the QCD crossover. We solve the tensor and scalar equations of motion across the Standard Model (SM) thermal history and tabulate the transfer functions for direct use in present and future PTA analyses. We show that this SM effect modifies the height of the induced spectrum by up to $\approx55\%$ across the PTA band relative to the radiation-domination expectation, with either sign depending on whether the source modes cross the horizon before or after the crossover. Fitting the NANOGrav 15-year data with a broken-power-law curvature power spectrum, we find that including the crossover shifts the inferred peak amplitude and scale by an amount that could already be relevant for the comparison with primordial black hole overproduction bounds. The importance of this SM effect will grow as the statistical uncertainty on the amplitude and scale shrinks with future, more sensitive PTA datasets, at which point neglecting it could significantly bias the inference.

The exploration of the early Universe is being transformed by the James Webb Space Telescope (JWST), which delivers unprecedented angular resolution at infrared wavelengths and opens a unique window redward of the K band (2um). Thanks to NIRCam, NIRISS, NIRSpec, and MIRI instruments, which provide both imaging and spectroscopy with exquisite efficiency, new classes of sources have emerged within the first years of operations. The first half-Gyr of cosmic time is now routinely probed, revealing massive blue/red galaxies and a population of Active Galactic Nuclei (AGN) appearing as "little red dots" together with a rest-frame near-infrared/optical view of sources across the reionization and post-reionization epochs. Angular resolution will remain pivotal in the 2030s - 2040s, when extreme adaptive optics (AO) facilities will be deployed on both (8 - 10)m (e.g., the VLT Multi-Conjugate-Adaptive-Optics (MCAO) - Assisted Visible Imager and Spectrograph, MAVIS) and on extremely large telescopes, like the 39m ELT (e.g., Multi-conjugate adaptive Optics Relay For ELT Observation, MORFEO). Operating at the diffraction limit, these facilities will improve JWST's resolution, with ELT achieving a factor of ~ 6 smaller Point-Spread-Function (PSF). An ELT diffraction-limited PSF (with a Full Width Half Maximum, FWHM ~ 8-12 mas) in the near infrared will resolve spatial scales <100 pc at any redshift (z ~< 18), revealing abundant star-forming clumps with sufficient sensitivity. Leveraging gravitational lensing as a cosmic telescope, even with moderate magnification factors (mu ~ 4-8), diffraction-limited 8m and 39m telescopes will probe physical scales ~< 25 pc, enabling systematic studies of star formation down to star-cluster scale at cosmological distances. Such observations are poised to become routine in the 2030s - 2040s.

We investigate linear matter-density perturbations in full quadratic gravity. Working in the longitudinal gauge and adopting the sub-horizon and quasi-static approximations, we obtain an algebraic expression for the effective gravitational coupling in the theory and the equation for evolution of matter density perturbations. We identify the positive-real poles of this coupling and map their dependence on the quadratic couplings. In the formal infrared continuation, the reduced coupling recovers its general-relativistic value, whereas at large physical wavenumber it is suppressed when both quadratic couplings are nonzero; the single-coupling limits reproduce distinct behaviors. Numerical integrations of the growing mode are performed up to the first pole and show that pole-free regions of the reduced system can support regular matter-growth solutions, while poles obstruct a continuous late-time evolution. Our results show that quadratic gravity cannot generally reproduce the standard cosmological evolution of perturbations and in the regions where it can, it closely matches the behavior of $f(R)$ theories.

APOGEE delivers high-resolution, near-infrared spectroscopy and highly precise radial velocities (RVs), with extensive repeat observations spanning more than a decade. Long-term RV repeatability is nevertheless limited by instrumental zero-point systematics: exposure-to-exposure drift and fiber-dependent offsets that are only partially tracked by lamp-based wavelength calibration, and that can evolve as fibers are serviced or replaced. We infer and remove these zero points directly from the published per-visit RVs using an uber-calibration framework based on shared repeat observations. We solve sequentially for two sets of minimal calibration units: (plate, night) exposure-level offsets and time-dependent (fiber, observing year) offsets. We apply the method independently to APO and LCO DR17 and to APO DR19. After correction, the same-fiber repeatability within two days reaches a limiting precision of ~14 m/s, and the long-baseline single-visit precision improves by a factor of two from 98.2 to 44.7 m/s. We release per-visit, value-added RV catalogs for DR17 (APO+LCO) and DR19 (APO) containing 3,501,727 visits (3,336,141 fully corrected) of 930,519 unique stars, enabling more reliable RV variability, binary, and Galactic-dynamics studies.

Determining isotopic ratios on the surfaces of extraterrestrial objects can provide insight into their formation conditions and evolution. Recently, remote-sensing detections of the 4.14 {\mu}m HDO absorption band have enabled astronomers to estimate the D/H ratio on several icy objects in the outer solar system. Here, we quantify how ice phase and radiation processing affect the 4.14 {\mu}m absorption band. While the band is visible in H2O-ice samples with a crystalline component, we cannot detect it in amorphous H2O-ice. Furthermore, we find that radiation processing, using 10 keV electrons as a proxy, quickly makes this absorption feature undetectable at a rate consistent with the amorphization of an initially crystalline sample. Interestingly, we also find that recrystallizing the irradiated sample causes this band to reappear nearly to its original band depth. For the Saturnian satellite Mimas, we estimate that the HDO feature will decrease by a factor of two within ~1 x 10^4 yrs at typical depths probed by remote sensing spectroscopy. While we suspect that existing methods could still determine the D/H ratio from a surface composed entirely of crystalline H2O-ice, deriving it for surfaces with a significant amorphous fraction may lead to severe underestimation of the true D/H ratio. However, given the direct correlation between the amorphous fraction of the sample and the HDO band depth, we propose that one could still estimate the D/H ratio of a surface with mixed phases by determining the crystalline fraction of the surface H2O-ice using other absorption features.

Detecting quadratic quasi-normal modes in black hole ringdowns would provide evidence for nonlinear gravitational dynamics, while measuring their properties would enable tests of the corresponding predictions of general relativity. Specifically, second-order black hole perturbation theory predicts that their amplitudes scale with the product of the amplitudes of their parent linear modes, with a coupling coefficient that depends on the spin of the remnant black hole. This mode-specific coupling coefficient has not yet been directly measured from gravitational wave data. Here, we use Bayesian inference on the GW250114 ringdown, modeled with the $220$, $221$, and $220\times220$ modes, to infer the amplitude-coupling coefficient of the $220\times220$ mode. The inferred coefficient is consistent with predictions from numerical relativity fits and second-order perturbation theory; no comparison shows a deviation exceeding $1.3 \sigma$. Although the current data cannot distinguish among these spin-dependent predictions, this measurement establishes a basis for more stringent observational tests of nonlinear gravitational dynamics in black-hole ringdowns.