We propose a machine learning-assisted analytical reconstruction of the Milky Way rotation curve and discuss its implications in a relativistic spacetime context. The rotation curve is reconstructed using 73 observational data points over the range 0.1-95.56 kpc, and we compare the performances of Ridge regression, LASSO regression, and feed-forward neural networks using a physically motivated functional basis. Ridge regression yields the most stable prediction with R2 = 0.9824 +/- 0.0064 and RMSE = 3.75 km/s, while retaining analytical interpretability. We embed the reconstructed velocity profile into a static, spherically symmetric spacetime, enabling the determination of the redshift function and the mass function through Einstein's equations. We verify that all energy conditions are satisfied, the sound speed remains subluminal, the circular orbits are stable, and the gravitational energy is negative, confirming the attractive nature of gravity. This framework provides a statistically validated, data-driven alternative to conventional dark matter halo models and establishes a direct connection between kinematical observables and relativistic spacetime geometry.

We show that a regular solution for primordial tensor fluctuations can arise from the higher-order multipoles of the collisionless neutrino distribution after neutrino decoupling. Focusing on the leading case with the neutrino octupole mode (tensor-octupole mode), we derive the initial conditions for the Einstein-Boltzmann system and calculate angular power spectra for the cosmic microwave background (CMB) anisotropies. Compared with the standard gravitational-wave mode, the tensor-octupole mode has a weaker large-scale metric response. It therefore gives a suppressed reionization bump and a different oscillation phase in the tensor CMB power spectra, providing a way to separate it from the standard tensor mode.

Large-scale structure provides a powerful probe of inflationary physics through primordial non-Gaussianity (PNG): the galaxy power spectrum depends on local PNG through scale-dependent bias, while the galaxy bispectrum depends sensitively on both equilateral and orthogonal PNG. In this paper, we study whether galaxy surveys can also probe folded PNG, whose shape is enhanced near $k_1+k_2-k_3\rightarrow0$. We consider three inflationary models with folded PNG, including excited initial states, imaginary speeds of sound, and dissipative inflation. These models fall into two classes: cutoff-regulated cases in which the folded-enhanced region has a power-law width, and dissipation-regulated cases in which the enhanced region is exponentially narrow. We develop a numerical pipeline for computing the corresponding PNG contributions to the redshift-space galaxy power spectrum multipoles and bispectrum monopole within the EFTofLSS. Using Fisher forecasts, we show that most of the constraining power on folded PNG comes from the galaxy bispectrum. For the cutoff-regulated models, nuisance parameter marginalization causes only a mild loss of information at large folded enhancement, but finite Fourier-space binning degrades constraints once the folded region becomes narrower than the bin width. For the dissipation-regulated models, the exponentially narrow folded enhancement is hard to resolve and the observable signal instead comes from the broader support of the template, leading to weaker binning dependence but larger overlap with the equilateral template. Our results show that folded PNG is a distinctive and promising target for galaxy bispectrum analyses, and the detectability depends on the width and morphology of the folded enhancement. The numerical pipeline developed in this work is general and can be used to study a wide class of non-separable primordial bispectra.

Measuring the total mass of the Milky Way and nearby galaxy groups is difficult because classical dynamical estimators rely on assumptions about satellite orbital geometry that are rarely satisfied in practice, and because only a handful of satellite galaxies are typically available as kinematic tracers. We present a new framework that corrects the well-known Projected Mass Estimator (PME) using a residual neural network trained on thousands of simulated galaxy groups from the IllustrisTNG cosmological simulation. Separate networks are trained for each satellite sample size, from as few as 5 satellites up to 50, so that the correction automatically accounts for the statistical noise that dominates when only a small number of tracers is available. In tests on simulated halos, the classical PME systematically overestimates halo masses by factors of $M_{\rm proj}/M_{\rm true} = 1.30^{+0.72}_{-0.62}$ (using the 2D distance) and $1.46^{+0.97}_{-0.72}$ (using the 3D distance), with RMSE of 0.29 and 0.32 dex respectively. The neural-network correction reduces this to $M_{\rm proj}/M_{\rm true} = 1.02^{+0.30}_{-0.26}$ with an RMSE of 0.13 dex. Applied to the Milky Way, the method yields a total mass of $M_{\rm MW} = 1.144^{+0.399}_{-0.296}\times10^{12}\,M_\odot$, with estimates based on the brightest 5-10 satellites favoring a somewhat lower range of $(0.8$-$0.95)\times10^{12}\,M_\odot$. The modified PME gives a tighter constraint on the virial masses and the dark matter rate prediction in galaxy groups and clusters.

The maximum physical extent that galaxies can reach is poorly understood. In this regard, IC 1101, one of the most extended and massive galaxies known, provides a valuable opportunity to constrain the upper limit of galaxy sizes at the present epoch. Previous deep imaging of the system confirmed its enormous extension, but did not indicate whether it has an edge. We explore this issue using the deepest images ever taken of this galaxy; ultra-deep g- and r- band imaging from the INT/WFC, reaching {\mu} = 30 mag arcsec^-2 (3{\sigma} in an area equivalent to 10 x 10 arcsec^2). We model and subtract the scattered light from both stars and the galaxy itself using an extended PSF characterization and a hybrid wavelet-based deconvolution. Using a combination of surface brightness, colour, and stellar mass density profiles oriented at different position angles, we find that the main body of IC 1101 extends to Redge = 260 kpc along the semi-major axis (assuming the redshift of Abell 2029, z = 0.077), enclosing 3.4 x 10^12 M_sun in stars. This Redge is among the largest edge radii measured for any galaxy to date, placing IC 1101 at the extreme upper end of the mass-size relation. In addition, we report a large number of asymmetrical, very low surface brightness features around the galaxy that are spatially consistent with the large-scale disturbances observed in the intracluster medium through X-ray studies of the Abell 2029 cluster, in which IC 1101 is embedded. With a confirmed diameter of around 520 kpc, IC 1101 stands as the largest galaxy known to date; yet, its outskirts show clear signatures of ongoing mass assembly, indicating that its spatial extent is still growing.

Neutral regions of the intergalactic medium (IGM) during the Epoch of Reionization (EoR) remain largely unexplored due to the limitations of existing probes. Owing to discoveries of numerous high-redshift radio-bright sources, the 21-cm forest, a series of absorption features imprinted by the neutral IGM in the spectra of such sources, now offers an attractive probe of the thermal and ionization state of the predominantly neutral IGM at $z\gtrsim5.5$. We analyse archival upgraded Giant Metrewave Radio Telescope (uGMRT) observations of J352-15, the brightest known radio-loud quasar in the EoR ($z=5.82$), to measure the one-dimensional (1D) power spectrum of the 21-cm forest. By comparing the observed power spectrum with forward-modelled synthetic spectra generated from cosmological simulations spanning a wide range of ionization and X-ray pre-heating scenarios, we perform Bayesian inference even in the absence of a statistical detection. We also present an independent Murchison Widefield Array measurement, although its lower sensitivity prevents competitive constraints. Using uGMRT, we achieve a sensitivity of $3.62\,\rm mJy\,beam^{-1}$ per $6.1\,\rm kHz$ channel. While we do not detect the 21-cm forest statistically, the null detection jointly constrains the mean neutral hydrogen fraction, $\langle x_{\rm HI}\rangle$, and the mean temperature of the neutral IGM, $\langle T_{\rm HI}\rangle$. At the $68\%$ credible level, our analysis disfavours cold and substantially neutral IGM models at $z\approx5.6$, including models with $\langle T_{\rm HI}\rangle \lesssim 27\,\rm K$ for $\langle x_{\rm HI}\rangle=0.1$. These limits probe parameter space allowed by existing Ly$\alpha$ and 21-cm observations, indicating substantial pre-heating of the neutral IGM above the adiabatic cooling floor. This demonstrates that the 21-cm forest has entered the regime of observationally informative statistics.

Type Ibn supernovae are stripped-envelope explosions whose spectra indicate interaction with dense, helium-rich and hydrogen-poor circumstellar material (CSM), making them important probes of late-stage mass loss and progenitor stripping. We present extensive ultraviolet-to-near-infrared spectrophotometry of the Type Ibn SN 2023tsz, including two epochs of HST/STIS ultraviolet (UV) spectroscopy and ground-based optical and near-infrared follow-up observations. The spectra are dominated by intermediate-width emission lines at all phases after maximum light, suggesting that much of the luminosity originates in a cold dense shell (CDS) formed by interaction between the ejecta and CSM. We compare the observations to one-dimensional non-local-thermodynamic-equilibrium radiative-transfer models of a helium-star explosion with a mass of $4 M_{\odot}$ at the onset of helium burning. The models reproduce the strong optical and near-infrared He I lines and require an added X-ray irradiation field to match the highly ionized UV features. The spectra are best reproduced by models with an X-ray irradiation power of $L_X \approx 10^8 L_{\odot}$, with the preferred models favoring CDS radii of order $(1.5$--$2) \times 10^{15}$ cm, velocities of $\sim 5 \times 10^7$ cm s$^{-1}$, and interaction powers of a few times $10^{42}$ erg s$^{-1}$. In the optical, the preferred models shift from higher interaction power and smaller radii at early times to lower power and larger radii at later times. These results add to the growing evidence that at least some SNe Ibn arise from lower-mass helium stars whose final evolution is shaped by binary interaction.

JWST observations, which have pushed the discovery and characterization of galaxies to cosmic dawn, have revealed significant deficiencies in state-of-the-art galaxy formation simulations, motivating the need for novel, physically grounded models. We present the first results from the Azahar suite of ten high-resolution ($\sim$20 pc), large-volume zoom-in cosmological simulations, which follow the formation of thousands of galaxies and progressively incorporate radiative transfer (RT), cosmic rays (CRs), and magnetohydrodynamics in addition to the 'standard' baryonic physics. Our Non-Thermal model, which simultaneously includes RT and CRs on the fly, reproduces the observed UV luminosity function from high redshifts ($z \sim 14$) to cosmic noon ($z \sim 3$), as well as the evolution of the stellar mass function, the galaxy main sequence, and observed gas metallicities. It does so through a combination of burstier star formation, with enhanced variability on both long ($\gtrsim 50$ Myr) and short ($\lesssim 10$ Myr) timescales, and a transition from high to low outflow mass-loading factors over cosmic time as CR pressure builds up in the ISM. We find that simple (boosted) thermal SN feedback models fail to capture the cosmic evolution of the very first galaxies because they lack realistic feedback channels that operate on different timescales and have distinct thermodynamical properties, thereby either overproducing stellar mass or driving overly explosive outflows. Our results indicate that non-thermal galaxy formation physics is crucial for providing a robust theoretical framework with which to interpret the high-redshift galaxy populations now being uncovered by JWST.

We introduce a geometric framework that organizes the second-order statistics of the Elasser fields into polarization states on generalized Poincaré spheres. In this representation, energy, cross-helicity, residual energy, and the phase lag between counter-propagating wave packets emerge as complementary polarization parameters. We derive a Bloch-analogue equation governing the evolution of polarization states and show that distinct polarization geometries are associated with different cascade dynamics. The framework predicts that transitions in the turbulent spectral scaling coincide with changes in the polarization state of the interacting modes.

The thermal state of interstellar gas controls whether it remains warm and diffuse or cools into dense clouds. Twenty-one centimetre observations have established a statistical picture of neutral-gas phases, but their three-dimensional architecture has not previously been mapped. Here we present $\mathcal{P}_{\rm 3D}$, a three-dimensional reconstruction of thermal phase in the local interstellar medium (ISM), covering a 1 kpc-diameter region centred on the Sun and sampled on a 2 pc Cartesian grid. $\mathcal{P}_{\rm 3D}$ combines high-resolution 3D maps of dust extinction and interstellar far-ultraviolet radiation with a neutral-ISM thermochemical model. The reconstruction reveals cold clouds surrounded by thermally unstable envelopes and embedded in a pervasive warm phase. Within $|z| \leq 150$ pc of the Galactic midplane, $\sim$41 % of the dust-traced neutral mass is thermally unstable, implying phase cycling on a timescale of $\sim$2.9-6 Myr. Yet the cold phase has a narrow density distribution, consistent with internal Mach numbers that are transonic at most ($\mathcal{M}_s \lesssim 1.5$). Together, these findings favour a dynamically cycling multiphase interstellar medium and motivate reassessing star-formation models that assume a single-phase, strongly supersonic cold-gas density field.

We propose a hybrid realization of Starobinsky inflation in which the inflationary epoch ends through vacuum decay. The model consists of an effective two-field system with a metastable Starobinsky branch shifted with respect to the true one. During the observable stage, the inflaton slow-rolls along the false branch, until a first-order phase transition in an orthogonal direction connects it to the true branch and ends inflation abruptly. This old-inflation-like exit skips the last part of the would-be Starobinsky trajectory. As a result, the Cosmic Microwave Background pivot scale exits the Hubble radius further from the minimum of the false branch than in ordinary Starobinsky inflation, raising the scalar spectral tilt $n_s$ while preserving the characteristic small tensor-to-scalar ratio. This provides a simple way of moving Starobinsky inflation towards the high-$n_s$ region favoured by recent ACT-related combinations. The same vacuum transition leaves a stochastic gravitational-wave relic whose peak frequency is controlled by the tunnelling timescale and the subsequent reheating history.

3I/ATLAS is only the third known interstellar object to traverse the Solar System. We obtained high-resolution ($R\simeq 140000$) spectroscopy of 3I/ATLAS on 23 nights with VLT/ESPRESSO ($r_{h}=1.67-2.45$au). CN production rates were derived by fitting the $B^2\Sigma^+-X^2\Sigma^+$ band, using CometSpec, a publicly released Python package for flexible fluorescence modeling. The three forbidden [OI] lines were modeled to derive the green-to-red (G/R) flux ratio and a CO$_2$/H$_2$O proxy. We fit $Q_{CN}=a(r_{h}/1\text{au})^b$ to the production rate versus heliocentric distance relation. The pre-perihelion CN power-law index is $b_{pre}=-4.62^{+1.25}_{-1.22}$, within the range observed for solar-system comets. The CN residuals about the smooth model show substantial scatter that we attribute to a combination of systematics and possible physical drivers that the present cadence cannot disentangle. The G/R ratio decreases from $\sim 0.42\pm 0.14$ at $r_{h}\simeq 2.4$au to $\sim 0.114 \pm 0.004$ near perihelion. The asymmetric piecewise oxygen fit is strongly preferred over the symmetric one, consistent with thermal-inertia models, although denser post-perihelion sampling would be needed to confirm this. The CO$_2$/H$_2$O proxy spans $\sim 0.06-1.06$ and is comparable to the values reported by Subaru/HDS and SPHEREx, and to those measured for 2I/Borisov. We additionally identify NiI and FeI, and report post-perihelion production rates for C$_2$ and CH. Our high-resolution monitoring of 3I/ATLAS reveals a CO$_2$-rich coma that becomes progressively H$_2$O-dominated near perihelion, with a G/R asymmetry consistent with delayed H$_2$O sublimation predicted by thermal-inertia models. Together with the public release of CometSpec and the reduced data, these results provide both an empirical reference and a methodological framework for the next generation of interstellar-object monitoring campaigns.

James Webb Space Telescope (JWST) has discovered unexpectedly bright, rapidly growing galaxies in the early universe, which were not predicted by most previously existing galaxy formation models. Using synthetic JWST observations of the Azahar simulation suite, we show that comprehensive non-thermal physics ("Full-Physics") produces compact, self-regulated galaxies that match observations from $z = 12$ to $z = 3$. This model also produces broad surface brightness distributions, where the bright end is dominated by compact sizes and bursty star formation. This compact starburst scenario naturally explains the detection of bright $z > 10$ galaxies in flux-limited surveys. By contrast, a model with standard hydrodynamics yields systems that are smaller and more concentrated than current data, while a model with calibrated supernova feedback produces unphysically large systems nearly twice the size of those observed. At lower masses ($\mathrm{M}_{*} < 10^{8}\,\mathrm{M}_{\odot}$), the Full-Physics model predicts sizes that are smaller than can be resolved with JWST, consistent with extrapolations from observations of higher-mass systems. Future observations with higher resolution could resolve this population and elucidate the physics driving the formation of the first galaxies.

Quasi-periodic eruptions (QPEs) are a novel class of repeating nuclear transients, discovered exclusively in the X-ray band to date. Since their initial discovery in 2019, the QPE sample has grown to 13 sources, exhibiting large amplitude, quasi-regular eruptive variability patterns that are distinct from previously known modes of massive black hole variability. In this chapter, we provide a comprehensive overview of their observational characteristics. We review the X-ray spectral and timing properties of QPE eruptions, their long-term evolution, as well as the underlying quiescent emission, which is well described by thermally dominated, compact accretion disks. We discuss population-level emerging trends and selection biases, and present an updated census of their host galaxy properties. We also highlight the growing body of evidence pointing to strong connections between QPEs and tidal disruption events. Finally, we briefly summarize the key observational constraints on proposed QPE model interpretations, before looking ahead to the observational challenges and opportunities that will shape future progress in understanding this emerging population of nuclear transients.

Quasi-Periodic Eruptions (QPEs) are recurring soft X-ray transients that may be generated by inspirals of stellar-mass objects spiraling into supermassive black holes, known as extreme mass ratio inspirals (EMRIs). Independently, EMRIs and the gravitational-wave signals they generate are one of the key targets for the Laser Interferometer Space Antenna (LISA). What is the potential of a coincident detection of EMRIs both as a QPE and by LISA? Electromagnetic counterparts to LISA events would provide sky localization, enable standard siren measurements of the Hubble constant and constrain formation mechanisms of the corresponding inspirals. Combined observations would link the accurate measurements of black hole masses and spins to their galactic nuclear environments and would thus enable lasting synergies with various observations across the electromagnetic spectrum. However, most of the currently known QPEs imply EMRI orbital periods that place the frequencies of the corresponding gravitational-wave signal out of the LISA sensitivity band. Additionally, the selection biases of QPE detections and the LISA instrument may preclude a coincident detection. Future searches should focus on expanding the QPE catalogs and ultimately hunt for ``golden'' short-period QPEs that correspond to EMRIs that fall within the LISA band.

The distribution of heavy elements in stars is described using a phenomenological approach, in which Lagrange parameters related to temperature and the chemical potentials of protons and neutrons are introduced within a freeze-out concept. Slope parameters are considered which describe the gross behavior of the distribution of the heavy elements. Universality and deviations from universality are discussed, and various examples are provided. These slope parameters may be of interest for characterizing the conditions under which heavy elements form, but the astrophysical sites where heavy elements are produced remain to be determined.

MITESI is an optical testbed which simulates key characteristics of the ELT, ultimately producing an artificial natural guide star. It was designed primarily to enable testing of the METIS instrument's SCAO system in closed loop. In this contribution, we discuss the assembly, integration and verification process of the testbed, taking the project from design to assembled hardware.

Luminous red novae (LRNe) are a class of optical transients resulting from the mergers of binary stars or common envelope events. The population displays heterogeneous light curves with extended luminosity plateaus and/or secondary peaks, suggested to be powered by hydrogen recombination or shocks between pre- and post-merger ejecta. However, much of their spectroscopic behaviour remains unexplained. Using a sample of six LRNe, we identify telltale spectral evidence for shock processes across the LRNe luminosity range ($10^{38}$-$10^{41}\,\mathrm{erg/s}$): (i) fast ejecta sweeping up slower upstream material; (ii) composite epochs with two distinct components, namely a cool stellar-like continuum underneath a hot nebular recombination region; and (iii) extremely broad line wings extending to $1000$-$10\,000\,\mathrm{km/s}$. Such line profiles reveal photons scattering off hot electrons in an outflow ($T_e\sim5000$-$10\,000 \,\mathrm{K},\,v_{es}\!\sim\!300$-$500\, \mathrm{km/s}$), with recombination seen from both the upstream and downstream of the shock surface. Additionally, constraints on the density of the upstream medium are obtained from recombination-line ratios, where brighter LRNe show deviations from Case B, consistent with denser surrounding environments. Given the outflow's velocity and surrounding density constraints, the shock luminosity is sufficient to account for the total energetics of LRNe without requiring other energy-injection mechanisms. During the plateau phase, the scattering wings are hidden and line ratios are affected by radiative transport effects, which highlights the utility of post-plateau phase spectroscopy.

We present updated transit timing measurements for three hot Jupiters (Qatar-4 b, HAT-P-18 b, and CoRoT-1 b) by leveraging data collected from the MicroObservatory Telescope Network, a network of small, robotic ground-based telescopes, and the NASA Transiting Exoplanet Survey Satellite (TESS). By combining these data with archival published results, we present the most precise orbital solutions to date for all three systems, allowing for precise transit time predictions for future missions. We report an updated mid-transit time for Qatar-4 b of 2458919.5838 $\pm$ 0.000089 $\mathrm{BJD}_{\mathrm{TDB}}$ and an updated orbital period of 1.80536560 $\pm$ 0.00000021 days. For HAT-P-18 b, we find a mid-transit time of 2459743.85340 $\pm$ 0.000022 $\mathrm{BJD}_{\mathrm{TDB}}$ and an updated orbital period of 5.50802957 $\pm$ 0.00000012 days. For CoRoT-1 b, we report a mid-transit time of 2456268.99083 $\pm$ 0.000099 $\mathrm{BJD}_{\mathrm{TDB}}$ and an updated orbital period of 1.50896846 $\pm$ 0.000000071 days. Our results demonstrate improvements over recently published ephemerides, with reductions of 36.4%, 4.35%, and 17.5% in mid-transit time uncertainties and 65.0%, 77.4%, and 16.9% in orbital period uncertainties for Qatar-4 b, HAT-P-18 b, and CoRoT-1 b, respectively. The results of this study improve the precision of future transit predictions and demonstrate the value of coordinated small-telescope monitoring (and citizen science initiatives) when updating the orbital parameters of hot Jupiters.

Abridged abstract: Context: Gamma-ray binaries hosting a non-accreting neutron star and a massive star exhibit multi-wavelength emission on different spatial scales. The interaction between their winds produces an outflow that can inflate a bubble or form a bow shock as it interacts with the surrounding medium. LS 5039 shows extended (1 pc-scale) X-ray emission that may arise from one of these large-scale structures. Aims: We explain and predict the large-scale emission from LS 5039. Methods: We modelled the thermal and non-thermal emission from five scenarios, representing different evolutionary phases, assuming particle acceleration at the mixed-wind termination shock: three bubble scenarios, treated with a one-zone model, and two bow-shock scenarios, studied with a multi-zone approach. We also investigated the radiation from escaping particles. Results: The extended X-rays are best explained as synchrotron radiation. Some scenarios predict detectable radio emission, while escaping particles may provide a minor steady contribution to the gamma rays in the powerful bow-shock scenario. Escaping protons with energies of 0.1-1 PeV could also inject up to $\sim10^{36}$ erg s$^{-1}$ into Galactic cosmic rays for optimistic injection luminosities. Conclusions: Gamma-ray binaries can efficiently accelerate particles on large scales, producing broadband emission and 0.1-1 PeV cosmic rays. Our results can guide future multi-wavelength observations to constrain the large-scale interaction, age, and birthplace of LS 5039.

Pre-main sequence (PMS) early-type stars are rare due to their rapid evolution. Confirming PMS B stars is accordingly valuable for constraining higher mass star formation scenarios. Although the Orion Nebula (ON) region offers an ideal laboratory for such studies, its population of B stars remains poorly characterized. We aim to determine spectral types, masses, and ages of B and early A stars surrounding the ON, to identify robust PMS candidates. We combine optical and near-infrared spectroscopy to estimate spectral types using Brackett, He, Si, and Mg lines. Photometry from Gaia, 2MASS, and WISE is used to construct colour-magnitude diagrams to fit isochrones via chi-square minimisation, yielding stellar ages and masses. Dynamical masses from eclipsing binaries are employed as independent constraints on the mass range allowed for each spectral subtype. We derive stellar ages and masses for 53 stars and spectral classifications for 48 of them, identifying 30 PMS candidates, including some that were previously assigned to luminosity classes I-III. Several stars show spectroscopic variability, potentially linked to circumstellar material or binarity. Our combined spectroscopic and photometric approach identifies robust PMS B-star candidates and provides a validated framework for distinguishing them from evolved counterparts - in this instance, refining the census of early-types young stars in Orion. The age distribution of the PMS candidates offers preliminary clues about the star formation history of the ON region.

Enabled by improved instruments with increasing sensitivity, the ongoing gravitational wave census now contains 259 binary black holes, numerous enough to unveil trends, substructure, and subpopulations which may provide key clues to their underlying formation mechanisms. In this work, motivated by evidence for multiple formation channels including hierarchical formation, we build a natively multi component mixture model for the binary black hole population, in which each component has an independently recovered rate, mass, spin, and spin misalignment model. (The components share a common redshift distribution.) Using a model carefully tuned to avoid parameter degeneracies, a powerlaw model plus five successively higher mass gaussians, we recover overall merger rates versus mass and trends versus redshift which are consistent with previously published results. Too, we recover previously identified overall trends versus spin: preferential alignment and low spin at low mass; large spin and isotropic spins at high mass. Critically, however, our multi-component model disagrees with previously published results, finding all components except the lowest mass are consistent with isotropy. Too, our multi-component model has a roughly hierarchical spectrum of gaussian mass peaks, but without the expected correlations between spin and mass expected from naked hierarchical formation

We report the first observations of a long-duration very-high-energy (VHE; $E > 100$ GeV) flare from BL Lacertae (VER J2202+422), taken with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). On October 15, 2022, the Fermi-Large Area Telescope (LAT) detected elevated GeV activity originating from this blazar. This triggered a multiwavelength campaign, which includes observations from VERITAS, Swift, NuSTAR, and select optical and radio observatories. VERITAS observed the source for a total of $\sim 9.8$ hours between September 1, 2022 and December 1, 2022. An analysis of these data yields a $\sim 28 \sigma$ detection of the source. While previously observed VHE flares from BL Lacertae have lasted on time-scales of minutes to days, VERITAS continued to detect flaring activity from the source for over a month ($\sim 40$ days) after the original flaring activity was detected with Fermi-LAT. Broadband spectral modeling shows that a synchrotron self-Compton (SSC) model with an external inverse-Compton (EC) component is preferred over a one-zone SSC model.

Context. In recent years, clusters have been observed that host multi-component haloes, both in non-merging and merging systems. The existence of these multi-component haloes suggests that there is no clear distinction between the single components. Aims. Abell 2244 is an intermediate-mass cluster that hosts a double component diffuse radio emission. The aim of this paper is to carry out a in-depth study of the diffuse radio emission to constrain its origin and characterize its main radio properties. Methods. In this work we present LOFAR HBA, MeerKAT UHF, and L-band observations of the cluster Abell 2244. We investigated the nature of the diffuse radio emission, combining high sensitivity radio data with XMM-Newton deep X-ray observations. We also used mock LOFAR observations to investigate contamination of the emission from faint radio sources. Results. We find an integrated spectral index of $\alpha^{1279}_{144} = 0.9 \pm 0.1$ for both components, where only the radio halo shows spectral steepening at higher frequencies. These values are comparable with the spectral indices observed in disturbed massive clusters. The outer component does not follow the same radio X-ray correlation as the radio halo, which suggests a different physical origin. Conclusions. By analysing the physical and morphological properties of the diffuse emission, we find that the characteristics of the outer component of the emission are intermediate between those of radio haloes and of known megahaloes. Hence, we speculate that the source is either a morphologically disturbed radio halo, caused by a minor merger interaction, or a megahalo but we cannot reach a final classification. From the mock observations, we find that it is unlikely that the emission is caused by faint sources at low resolutions.

Radio flares from tidal disruption events (TDEs) are observed on two distinct timescales: prompt (around the time of optical peak) and delayed (hundreds to thousands of days after optical peak). A recent framework put forward to explain this diversity suggests that super-Eddington winds produce prompt flares while jets launched during a disk state transition at low accretion rates produce delayed flares. A unique observational prediction of this framework is that individual sources showing both types of flares may be expected to show a third flare. This third flare will occur when the (faster) jet catches up with the (slower) wind. The observational appearance of the third flare depends on the system parameters. As the two shock fronts collide the (brighter) jet encounters a change in density, modifying its synchrotron-self-absorption frequency and observed flux. Here we discuss the observational appearance of such a scenario, and then show, in new long-term monitoring radio observations of the TDE AT2020vwl, the discovery of a third radio flare. This flare is coincident with the time at which forward modeling of the first two flares in a wind-jet framework predicts a collision between the two shock fronts. We make predictions for the long-term radio evolution of two other TDEs showing double radio flares: ASASSN-15oi (which may show a weak flare or flattening of the lightcurve) and AT2024tvd (which is predicted to imminently undergo a third flare).

While masses and spins are routinely used to constrain compact binary formation channels, eccentricity provides an additional and potentially powerful diagnostic of binary origin, particularly for dynamically assembled systems. Recent advances in eccentric waveform modeling now make it possible to search for eccentric signatures in gravitational wave data; however, differences between waveform models can introduce systematic effects that may propagate into astrophysical population inference. In this work, we analyze 153 binary black holes, 2 binary neutron stars and 7 neutron star black hole binaries from the GWTC-4 catalog. We compare the source and population level inferences obtained with two eccentric waveform models, SEOBNRv5EHM and TEOBResumS-DALI, as well as with quasi circular waveform analyses. We find that the two eccentric models give broadly consistent source parameter estimates for most events, but some events exhibit subtle and coherent differences. These small, systematic offsets can accumulate in hierarchical population inference, leading to differences in inferred population properties, most notably in the redshift evolution and effective spin distribution. Because coherent event level biases can grow approximately as $\sqrt{N}$ for a catalog of N events, waveform systematics become increasingly important as gravitational wave catalogs expand. We also introduce a synthetic data framework that generates eccentric populations and corresponding RIFT posterior samples, enabling injection studies that test the recoverability of eccentric population properties.

AT2018cow, the prototypical luminous fast blue optical transient (LFBOT), shows a slowly fading plateau in the UV/optical light lasting at least 5 years after the explosion. The plateau SED is blackbody-like with temperature $\gtrsim$ 2e4 K and an emission radius ~40 Rsun, which has been attributed to a geometrically thin accretion disk around a compact object. Viscously powered disk models, however, require an unusually massive black hole $\gtrsim$ 200 Msun (for viscosity parameter $\alpha$ > 0.01). Moreover, accretion onto such a massive black hole produces X-ray emission that is likely in tension with the XMM-Newton constraint, unless the X-ray is highly anisotropic or absorbed. We propose instead that the UV/optical emission arises from the photosphere of a super-Eddington wind launched from the inner disk, together with irradiation and reprocessing by the geometrically thin outer disk. Fitting the model to four-band HST photometry at t ~ 4 yr, we find that all accretors from 1.4 Msun neutron star to ~100 Msun black hole produce similarly good fits. This degeneracy arises because the wind photospheric color temperature depends only weakly on the accretor mass. We also find that, at NIR-MIR wavelengths (1 to 10 um), wind free-free emission differs by an order of magnitude across the allowed mass range, providing a potential discriminant accessible to JWST. Another potential signature to differentiate the accretor mass and disk composition is the recombination line emission (HeII1640 or H$\alpha$) from the disk wind. The inferred disk parameters from our model are consistent with a micro-tidal disruption event, either tidal disruption or merger with a companion star, in which an initially geometrically thick disk transitions to a thin phase within the first year and subsequently the thin outer disk evolves on the observed ~10 yr viscous timescale.

We present an infrared, millimetre, and radio survey of 20 Class 0-III young stellar objects in the Ophiuchus A L1688 star-forming cluster, combining high-resolution (7-25 au) VLA and JWST observations with archival ALMA data. We implement physically motivated models to derive dust and ionised gas properties, spectral behaviour and their relative contributions through the millimetre-centimetre radio spectral energy distribution. Our models reveal circumstellar dust disks that are, on average, tens to hundreds of times more massive than millimetre-only estimates (subject to uncertainties arising from the choice of dust opacity) and contain millimetre-sized grains even at the Class 0 stage. Owing to the VLA's high resolution we are able to connect outflows to their origins, detecting protostellar jet emission on scales of 10s-1000s au. Our results represent a homogeneous characterisation of the dust and ionised gas properties in Ophiuchus and present a potential solution to the long-standing 'missing disk mass' problem. However, our understanding is still limited by resolution and sensitivity at frequencies <40 GHz. Future facilities like the SKA and ngVLA are needed to provide the necessary capabilities to fully spatially resolve this emission (<0.18") even in one of the closest star-forming regions.

Stars set the conditions for planet formation, planet evolution, and habitability -- all major topics in astronomy today. Stars are also important in their own right as the most visible component of galaxies. In cool stars, X-ray emission is powered by magnetic fields, and so far our Sun is the only system in which those fields are spatially resolved. High-resolution X-ray (HiReX) imaging can track the origin and evolution of those fields, see how they connect young stars to their disks and outflows, and measure the energy, mass, and momentum that radiation and coronal mass ejections (CMEs) carry into the circumstellar environment. This is crucial for understanding whether planets can form and survive in young stellar systems, whether life can develop on those planets, and how the star evolves over time. Intermediate-mass and high-mass stars blow winds and eventually evolve into degenerate objects such as white dwarfs, neutron stars, and black holes. Their evolution and death drive the chemical evolution of galaxies. High-resolution X-ray (HiReX) imaging can study the hottest components in those systems, such as the colliding winds of massive stars, accretion and nova explosions in CVs, and the shocks in outflows that form planetary nebulae. All these cases have in common that the X-ray emission is tracing the hottest, fastest, and most energetic components of the shocks. HiReX observations can reveal the temperature, spatial structures, and elemental abundances of different system components that no other wavelength can. While stars are physically small compared to more powerful objects such as accreting black holes and AGN, they are also much closer to us, allowing a HiReX mission to resolve a variety of physical phenomena fundamental to our understanding of how stellar systems form, evolve, and interact with their environment.

We analyze 405 deep JWST/NIRSpec spectra of star-forming galaxies at $z=4.5$--$10.1$ from DREAMS and other lensing-cluster surveys to study chemical enrichment in intrinsically faint, low-mass galaxies. The sample covers $M_{\rm UV}\simeq-12$ to $-22$ and reaches $M_\star\simeq10^{5.7}\,M_\odot$, with 50\% of the sources at $M_{\rm UV}>-17.5$ magnified by $\mu>3$. From individual spectra and mass-binned stacks, we derive the gas-phase metallicity together with nitrogen and carbon abundances using the rest-frame UV and optical lines. \Add{We derive new empirical strong-line metallicity calibrations using direct-method measurements from stellar-mass-binned stacks, reaching a representative stellar mass of $M_\star\simeq10^{6.6}\,M_\odot$. Applying these calibrations, we trace the $z\sim6$ mass--metallicity relation down to $M_\star\simeq10^{6.6}\,M_\odot$, where it reaches $12+\log({\rm O/H})\simeq7.2$, with a low-mass slope slightly steeper than the local relation and in broad agreement with hydrodynamical simulations.} In the $M_\star\simeq10^{7.7}\,M_\odot$ stack, N/O from N\,{\sc iv}]\,$\lambda\lambda\,1483,1486$ exceeds that from [N\,{\sc ii}]\,$\lambda6583$ by $\simeq1.4$ dex. The UV--optical difference could indicate a localized, highly ionized N-rich component whose high N/O and subsolar C/O resemble nitrogen-rich globular-cluster populations with $M_\star\lesssim10^6\,M_\odot$. \Add{The combination of these abundance patterns and a He\,{\sc ii}\,\W4686/H$\beta$ ratio of $\simeq0.03$, well above BPASS predictions, suggests that WR stars may contribute both prompt CNO-cycle enrichment and hard ionizing radiation, with the inferred WR population potentially supplying enough nitrogen to account for the excess on globular-cluster scales.

We present the results of a survey program with JWST's MIRI LRS that obtained 5$-$14 $\mu$m R~100 spectra of twenty-three L0$-$T6 dwarfs. This spectral type range spans the formation, growth, and sedimentation of silicate condensate clouds, as well as the appearance of methane and ammonia in the atmospheres of brown dwarfs. The initial onset of silicate absorption at spectral type L0 as well as the evolution of this feature with spectral type are both examined in unprecedented detail. We detect 7.5$-$9.5 $\mu$m gas-phase SiO absorption for the first time in dwarf stars, with detections in M5.5$-$L2.5 dwarfs from combined Spitzer IRS and JWST spectra. This feature has previously been detected in K and M giants, with stronger absorption in these giants than in higher-gravity M dwarfs. In contrast, low-gravity L dwarfs show weaker SiO absorption than field-gravity L dwarfs. We also confirm water, ammonia, and methane absorption trends with spectral type previously seen in Spitzer IRS observations. We report the possible detection of CS2 in a T3.5 dwarf, although confirmation would require higher-dispersion spectroscopy. We find two newly resolved binaries, confirm two near-equal flux unresolved candidate binaries, and identify two additional possible non-equal flux binaries. After accounting for contamination by resolved and candidate binaries, we find that single field dwarfs are ~0.15$-$0.20 mag underluminous compared to existing empirical relations. The presented spectral library offers a spectrophotometrically calibrated basis for future interpretation of JWST MIRI spectra of brown dwarfs or of directly imaged or transiting exoplanets.

Stellar rotation plays a key role in the transfer of angular momentum, and a large sample of Be-type stars with reliable projected rotational velocities is crucial for understanding their formation and evolution. In this work, we derive the projected rotational velocities ($v$\,sin\,$i$) of 479 Be-type stars using the Fourier transform method, based on their LAMOST Medium-resolution Survey (MRS) spectra. Our results suggest that the Fourier transform method can provide reliable $v$\,sin\,$i$ values for Be-type stars by analyzing the \ion{He}{1}\,lines at 4922, 5015, 5047, and 6678 \,Åin their LAMOST MRS spectra. A K-S test indicates that Be-type stars with different H$\alpha$ emission line morphologies exhibit different $v$\,sin\,$i$ distributions, and Be-type stars with double-peaked emission have a higher fraction of rapid rotators than those with single-peak emission. The $v$\,sin\,$i$ distributions of our Be-type stars in the field, OB associations, and clusters show no significant differences. The deconvolved $v$\,sin\,$i$ distribution of our entire Be-type star sample does not exhibit a bimodal distribution but rather a single peak at $v\approx260$\,km$\cdot$s$^{-1}$. Based on the analysis of 105 stars in our sample, we find that the mean equatorial rotational velocity is 0.74 times the critical velocity. Furthermore, we investigate the relationship between $v$\,sin\,$i$ and the H$\alpha$ peak separation velocity for Be-type stars exhibiting double-peak H$\alpha$ emission lines, using Pearson and Spearman rank correlation coefficients.

Changing-look active galactic nuclei (CL AGNs) show the appearance or disappearance of broad emission lines on timescales of years. Among them, the repeating CL (RCL) AGNs may provide a clear clue for our understanding of the CL transitions because the same nucleus crosses some physical boundaries more than once. We search for RCL AGNs in known CL-AGN samples using long-term multi-band light curves, and selected 34 candidates for spectroscopic follow-up. We confirm 25 RCL AGNs, including 22 newly identified cases. Properties of these RCL AGNs are analyzed. The observed rest-frame intervals of the second transitions are mostly 3--4 yr, while the variations of the optical light curves suggest that some transitions might occur on timescales of several months. The latest spectra show that the on/off states correspond to higher/lower Eddington-ratio, in the expected direction relative to the parent CL-AGN samples. As seven RCL AGNs are well covered by nearly continuous single-band light curves, their on/off states can be found to follow multi-year optical excursions, and their Eddington ratios vary consistently with the photometric changes. We also find that the H$\beta$-only transitions occur at higher Eddington ratios than the transitions involving both H$\alpha$ and H$\beta$, suggesting a line-dependent Broad-Line-Region (BLR) visibility threshold. These results support a picture in which different accretion-flow processes drive reversible changes in the central ionizing emissions, while the observed RCL transitions are produced by the BLR breathing across line-dependent visibility thresholds.

We present a data-driven measurement of galaxy bias and scale-dependent relative bias for elliptical and spiral galaxies using angular auto- and cross-power spectra from DES and DESI Legacy Imaging Surveys DR8. We introduce the cross-tracer clustering ratio (CTCR), which uses the ratio of auto- to cross-power spectra to isolate the relative clustering of morphological tracers as a function of angular multipole $\ell$. Across both surveys, ellipticals are more strongly clustered than spirals; the cleanest CTCR constraints come from DESI, where the two morphological samples have better-matched redshift distributions. The difference is scale dependent: the relative bias is close to unity on large angular scales, $\ell \lesssim 50$, but increases toward smaller scales, reaching an average separation of $2.6\sigma$ at $\ell \sim 150-200$. A complementary linear-bias analysis confirms that ellipticals are positively biased relative to the matter field, while spirals are consistent with weak bias or anti-bias. Unlike narrowly selected LRG samples, our morphologically selected elliptical samples show little redshift evolution, consistent with a broader halo-mass distribution. The results are robust to two covariance estimators, contamination tests, luminosity splits, and comparisons with two N-body mock catalogs. These new findings provide empirical evidence that galaxy morphology imprints both the amplitude and scale dependence of galaxy bias, and establish CTCR as a useful observable for testing halo occupation and assembly-bias models with future imaging surveys.

The density dependence of the speed of sound in cold neutron star matter remains poorly constrained and is central to determining the high-density equation of state (EOS). While binary neutron star (BNS) merger simulations increasingly incorporate detailed microphysics, the direct impact of the sound-speed in the neutron star core on observable signatures has not been systematically explored. We address this by introducing a simplified parametrization that suppresses microphysics while allowing controlled variation of the sound speed at supranuclear densities through its derivative with respect to the baryon chemical potential. Using the WhiskyTHC code, we perform a suite of fully relativistic BNS merger simulations and identify correlations between the sound-speed slope parameter and key merger outcomes, including remnant properties and post-merger GW frequencies. These results demonstrate that multimessenger observables are sensitive to the behavior of matter at the highest densities reached in neutron star cores. We further analyze gravitational-wave signals from the CoRe database of binary neutron star merger simulations employing more realistic equations of state. Our analysis reveals approximately EOS-independent relations between the derivative of the sound speed and the peak post-merger gravitational-wave frequency. Although these relations cannot be considered truly quasi-universal, they nonetheless indicate that post-merger gravitational waves retain measurable information about the EOS at he highest densities. At the same time, the remaining EOS dependence highlights the difficulty of isolating the underlying high-density physics and motivates the development of targeted parameterized frameworks for interpreting future multimessenger observations.

We present the first blind search for OH 18-cm absorption with the Five-hundred-meter Aperture Spherical Telescope (FAST), conducted alongside the HI 21-cm absorption search. Our previous FAST blind HI absorption search identified 34 systems. In this work, we extend the search using 2024 and part of the 2025 CRAFTS and FASHI data (394.4 hr and 1622.1 deg$^{2}$) together with FATHOMER observations, yielding three known and four new HI absorbers, for a total of 41 HI absorption systems. We search for OH absorption in 19 HI absorption systems whose OH redshifted frequencies fall within the FAST band. The known OH absorber towards PKS 1413+135 was re-detected, making our survey the first blind survey to detect OH absorption. No new OH absorbers were identified. We examine the relationship between $N_{\rm{OH}}$ and $N_{\rm{HI}}$, applying survival analysis to account for upper limits. The analysis does not provide statistically significant evidence for either an $N_{\rm{OH}}$-$N_{\rm{HI}}$ correlation or redshift evolution of $N_{\rm{OH}}$/$N_{\rm{HI}}$. Finally, spectral stacking sets 3$\sigma$ OH column density upper limits of 4.93, 1.64, and 1.72 $T_{\rm{ex}}$/$c_{\rm{f,OH}}\times$10$^{12}$cm$^{-2}$K$^{-1}$ for associated, intervening, and combined samples, corresponding to [\rm{OH}]/[\rm{HI}] ratios of $<$1.66$\times$10$^{-8}$, $<$1.42$\times$10$^{-8}$, and $<$0.90$\times$10$^{-8}$, assuming $T_{\rm{ex}}$=10K for OH and $T_{\rm{s}}$=100K for HI. These results place the strongest constraints to date on the OH content in radio-selected HI absorbers and establish a blind-survey benchmark for future studies of molecular gas in HI-selected systems. They also demonstrate that known HI 21-cm absorbers provide an effective parent sample for systematic OH absorption searches, paving the way for future larger surveys.

Sunspot oscillations probe magnetic-plasma interactions in the solar atmosphere. We investigate oscillation periods in six sunspot active regions observed in the 8-10 um band with the Accurate Infra-red Magnetic-field Solar Telescope (AIMS) at Lenghu Observatory. Light curves were extracted from umbra, penumbra, and quiet Sun regions. Wavelet analysis was applied using four methods: pixel-wise with mean/median aggregation, and spatial mean/median prior to analysis. The weighted mean period is more physically meaningful than the dominant period, better representing multi-mode oscillations. A consistent period sequence holds across all datasets: umbra (U) < penumbra (P) < quiet Sun (Q), with typical values of 260-313 s, 286-374 s, and 294-382 s. This pattern is recovered with 100% consistency by pixel-wise mean/median and spatial median methods. Spatial mean shows 1/6 full, 3/6 partial, and 2/6 inconsistency. Spatial smoothing increases periods by up to 21%, most strongly in the umbra. All peak ratios are below 0.3 (mean 0.14 +/- 0.03), confirming multi-mode oscillations. These results establish the 8-10 um band as a valuable diagnostic window for sunspot magneto-convection and wave propagation.

We present synthetic images of galaxies that are in the stages of star formation quenching for the Cosmological Advanced Survey Telescope for Optical and UV Research (CASTOR) and Nancy Grace Roman Space Telescope (NGRST), based on simulations coming from the IllustrisTNG suite, as processed using the stellar population synthesis library \textsc{galaxev}. We account for the effects of dust and various sources of noise to produce mock observations that should mirror real observations. Using these synthetic images, we fit photometric observations in binned circular annuli using \texttt{FAST++} and a flexible star formation history, and recover well the spatially resolved stellar mass and star formation rate. We thereby measure various indicators (morphological metrics) of spatially resolved star formation activity in the context of galaxy quenching. We find that we are able to distinguish quenching galaxies from a mass-matched control sample of normal star forming galaxies. We additionally find that we can distinguish various quenching mechanisms, where galaxies consistent with an inside-out quenching signature can be separated from galaxies that display an outside-in signature. Using machine learning techniques the accuracy of this classification is reliable, and the progress through the quenching episode can be estimated for the different populations of quenching galaxies. We make predictions for the abundance of the various quenching populations in proposed surveys for CASTOR and NGRST, and find that these surveys will enable the classifications of thousands of quenching galaxies out to intermediate redshifts, and more when considering higher redshifts.

The Sun was almost certainly born in a stellar cluster, implying some degree of external forcing from neighboring stars during the Solar System's infancy. Published estimates of the strongest relevant stellar encounter, however, span a broad range, from relatively gentle perturbations to violently disruptive flybys. The modest inclination dispersion of the distant trans-Neptunian population has previously been used to argue that strong primordial encounters were unlikely and that the outer Solar System was not violently stirred at birth. In this work, we independently examine this constraint and ask whether this conclusion holds in the presence of Planet Nine. Because Planet Nine can reshape the distant trans-Neptunian population over gigayear timescales, the extent to which the present-day inclination distribution preserves a direct record of primordial cluster excitation is a priori unclear. To address this question, we compare long-term N-body integrations in which distant test particles begin either in strongly cluster-perturbed configurations or in comparatively quiescent ones, and are then evolved under the influence of the known giant planets, the Galactic tide, the effects of passing stars, and Planet Nine. We find that Planet Nine does not transform a highly excited primordial population into one resembling the observed low-inclination distant sample. Instead, comparisons between the simulated and the observed distant TNO populations indicate that the observed distant TNOs are most consistent with a relatively mildly perturbed birth environment, even in the presence of Planet Nine. If this comparatively low inclination dispersion is confirmed to be an intrinsic feature of the distant TNO population, rather than an artifact of observational bias, it constitutes a robust constraint on the severity of primordial cluster perturbations.

We review the present status of observations of the Cosmic Microwave Background. This book chapter is an update to an earlier review "Physics of the Cosmic Microwave Background Anisotropy" that appeared a decade ago in the compilation "One Hundred Years of General Relativity: From Genesis and Empirical Foundations to Gravitational Waves, Cosmology and Quantum Gravity." The earlier contribution summarized the state of the field then as well as the underlying cosmology and gravitational physics. This contribution reports on what new has occurred in the field during the intervening time 2015-2025. We also review future prospects.

The lunar gravitational-wave antenna concept uses the Moon as a resonant detector instrumented with precision seismometers, targeting the decihertz band between ground- and space-based observatories. We propose a compact monolithic fused-silica torsional-X seismometer that re-engineers garden-gate acceleration-to-rotation transduction for this regime through a high-tension dual-fiber suspension. Its designed millihertz-scale resonance and ultra-low mechanical dissipation enable a nearly order-of-magnitude improvement around $0.1\,\mathrm{Hz}$ compared with existing lunar seismometer concepts. Achieving this performance requires room-temperature operation, where fused-silica exhibits low mechanical loss, together with subdominant actuation noise. We demonstrate a room-temperature vacuum prototype validating the operating principle and core mechanical design, and derive requirements for a future lunar implementation capable of approaching the target sensitivity.

With the rapid surge in lunar exploration and the planned expansion of cislunar infrastructure, cislunar space has become a strategic focal point for global aerospace activities. This proliferation of spacecraft heightens the risk of fragmentation events, such as unintended explosions or orbital collisions which serve as the primary source of hazardous orbital debris. Given the potential threat these fragments pose to mission safety and long-term orbital sustainability, it is imperative to investigate their dynamical behavior within the Earth-Moon system. This study evaluates the dispersion of debris clouds following potential breakup events on Distant Retrograde Orbits (DROs) over a 30-day propagation period. The Circular Restricted Three-Body Problem (CR3BP) model is used to construct the reference orbits, while the NASA Standard Breakup Model is applied to simulate fragment generation at multiple locations along three DROs of varying sizes. These fragments are then propagated using the Bicircular Restricted Four-Body Problem (BCR4BP) for 30 days. To account for the variability of these events, multiple initial positions along each orbit are analyzed to capture a comprehensive range of post-explosion scenarios. Our analysis quantifies the fate of fragments within this window, specifically focusing on the escape mechanisms and the percentages of debris that either depart from the Earth-Moon gravitational sphere of influence or impact the lunar surface. Furthermore, we introduce an analytical approach to assess the potential collision risk to resident space objects operating within the vicinity of the parent orbit. The results provide insights into debris evolution and offer a foundation for developing safety guidelines for future cislunar activities.

Solar flares are intense energy release events in the solar atmosphere that may pose significant space weather hazards, which makes developing reliable prediction models essential. Although deep learning methods, particularly convolutional neural networks (CNNs), demonstrate strong predictive performance when using solar magnetograms, their scientific credibility is undermined by a lack of physical interpretability. Explainable artificial intelligence (XAI) offers a potential solution. However, current XAI studies in solar flare prediction are largely qualitative and lack systematic, theory-based, quantitative validation. We present a quantitative XAI framework that can decipher the physical basis of CNN-based solar flare prediction models. Using gradient-weighted class activation mapping (Grad-CAM), we identify model-focused regions (MFRs) in solar magnetograms. Then, we perform two key analyses to evaluate the predictive capability of magnetic parameters derived from MFRs and to quantitatively characterize their magnetic complexity. Our results reveal a strong physical correlation between MFRs and flare occurrence. Specifically, magnetic features extracted from MFRs demonstrate high predictive power for flares. Flare-producing active regions are characterized by magnetically complex configurations that are dominated by a single polarity rather than by balanced or purely unipolar structures. This finding is consistent with established physical theories of magnetic systems prone to flares. Our results suggest that CNNs can learn physically meaningful representations when trained on large-scale observations. Integrating XAI with quantitative magnetic field analysis improves the physical interpretability of deep learning-based flare prediction models, making them useful tools for prediction and modeling investigation in solar physics.

Aims. We investigate the evolution of X-ray reverberation lags in the black hole X-ray binary Swift J1727.8-1613 during its 2023 outburst, with the aim of probing the inner accretion flow geometry across spectral states. Methods. We analyzed NICER observations covering the low-hard state (LHS) and hard-intermediate state (HIMS). The time lags were computed using Fourier-based techniques, and we constructed lag-frequency and lag-energy spectra. To obtain a robust estimate of the soft lag amplitude, we focused on a frequency range in which the reverberation signal dominates and remains stable, thereby minimizing contamination from hard lags and phase-wrapping effects. Results. The soft lag amplitude increases rapidly from the LHS to the early HIMS and then stays near 10 ms throughout the HIMS. In the frequency range in which reverberation lags prevail, the lag shows little dependence on Fourier frequency. On the other hand, the observed low-frequency lags change from hard-lag dominated to soft-lag dominated, with amplitudes comparable to those measured at higher frequencies. Conclusions. These results suggest that the reverberation lag varies little during the HIMS, consistent with a relatively stable inner accretion geometry during this state. The apparent evolution of the lag amplitude from the LHS to the HIMS can be largely explained by the diminishing effect of hard lags and does not necessarily require significant changes in the intrinsic light-travel timescale. Swift J1727.8-1613 therefore provides a case in which the reverberation signal can be studied with reduced contamination over a broad frequency range, offering new insight into the evolution of the accretion geometry in black hole X-ray binaries.

To study interstellar gas and grain chemistry, gas-phase astrochemical models have been developed since the early 1980s and gas-grain models since the 1990s. Each published model includes various assumptions mainly for surface processes. In this paper, we compare recently added mechanisms for grain surface chemistry, namely, non-diffusive chemistry, radiolysis, and chemical desorption. Several formalisms for these processes have been added to our astrochemical model Nautilus, and we tested them, comparing the predicted gas-phase and ice abundances. Our predictions are also compared to gas and ice observed compositions. Our main findings are that radiolysis itself does not influence the results. Non-diffusive chemistry can have an impact on the gas-phase and ice species, but it depends on the adopted formalism. In particular, the one of Shingledecker & Herbst (2018) changes the main reservoirs of the species in the ices, impacting the species in the gas-phase as well. The adopted formalism for chemical desorption can produce differences in the gas-phase by up to a factor of ten. Last, our standard model, without non-diffusive chemistry and with the chemical desorption from Fredon et al. (2021), produces the best results in relation to observed gas-phase abundances, while the ice observed agreement is unchanged. The formalism for some grain surface processes are important even for gas-phase abundances. More experiments are needed to constrain their efficiency, however. For the chemical desorption, each formalism relies on an uncertain parameter, which is the fraction of the energy actually delivered to the products, that can be adjusted to reproduce the experiments.

In this paper, we first devise a geometrical model, featuring a torus-like flux rope based on the shape of 3DCORE model. The global shape of the torus is an ellipse, while the cross sections are circular along the torus. The thinnest point is located between the Sun center and photosphere. Deflections and inclination are considered as well. Using multiwavelength observations from perspectives of Earth, Ahead-STEREO (STA), and Solar Orbiter, we apply the model to three-dimensional (3D) reconstructions and tracking of the filament eruption, which was associated with a flare and a coronal mass ejection (CME) on 2024 October 8. The morphology, direction, and true velocity ($\sim$433 km/s) of the eruptive filament are obtained. It is found that the filament propagates nonradially, deflecting slightly eastward by $\sim$10 degrees and significantly southward by $\sim$40 degrees. Trajectory of the filament in the ecliptic plane reveals that the filament moves toward STA. The true direction of the eruptive filament using imaging and spectral observations is mutually verified by 3D reconstructions. The heliocentric distance of the filament increases from $\sim$1.68 to $\sim$2.94 solar radii within 35 minutes. Based on the results of 3D reconstructions, the true speed of the CME leading edge is evaluated to be 1046$-$1145 km/s.

Universal Time (UT1) is a key parameter characterizing Earth's rotation, and very long baseline interferometry (VLBI) is the mainstream technique for measuring UT1. To address the limitations in the timeliness and accuracy of existing global ionospheric models for single-frequency VLBI UT1 measurements, we construct a single-station regional ionospheric model using GNSS data from the VLBI stations on the Jilin-Kashi baseline. We apply this model to VLBI observations and compare its correction performance with that of a global predictive model and a global post-processed model. The results show that the line-of-sight ionospheric delays and baseline corrections calculated with the single-station regional model have precision close to that of the global post-processed model and are substantially better than those of the global predictive model. After correction with the single-station regional model, the derived UT1 values differ from the US Naval Observatory (USNO) reference values by a mean bias of -15.6 us and an RMS deviation of 82.3 us, both better than the results obtained with the other two model classes. A single-station regional ionospheric model constructed independently from GNSS data available at VLBI stations can effectively correct single-frequency VLBI observations and support quasi-real-time high-precision UT1 measurements. It therefore has important value for improving the timeliness of independent UT1 products.

Context. The upcoming GaiaNIR mission is currently planning to add a near-infrared spectrograph to its payload in order to enhance its scientific return, particularly for mapping the dust-obscured regions of the Milky Way. Aims. This study aims to identify the optimal wavelength region between 800 and 2300 nm for the proposed GaiaNIR spectrograph to maximize the precision of radial velocities and atmospheric parameters. Methods. To find its spectral range, we generated 10000 synthetic spectra from the BOSZ library across a wide range of stellar parameters, with resolutions varying from 5000 to 20000. By cross-correlating these mock observations with ideal templates, we assessed the statistical scatter of velocity residuals to isolate six candidate windows for further atmospheric parameter testing. Results. Our analysis finds that the 1926 - 1968 nm window at R = 16100 - 20100 in the K-band is the preferred strategic choice, because it has the potential to reach radial velocity precision for the brightest FGKM stars of about 160 - 260 m/s depending on resolution, while providing precision of the atmospheric parameters close to Gaia's Radial Velocity Spectrometer. We also identified a second region between 1158 and 1202 nm (R = 9300 - 11600), that has slightly lower radial velocity precision, but at a wider temperature range than the K-band. Both regions make it possible to derive abundances of ten species at these resolutions: O, Na, Mg, Si, Ca, Ti, V, Cr, Mn, Ni in the K-band, and Mg, Si, K, Ca, Ti, V, Cr, Mn, Co, Ni between 1158 - 1202 nm. Conclusions. The K-band delivers sufficiently precise measurements for the mission's primary cool star targets while taking advantage of significantly lower interstellar extinction, enabling the mapping of the dust-obscured regions of the Milky Way.

Cosmic rays are widely held to drive outflows from star-forming galaxies and profoundly influence galaxy evolution. Direct evidence for cosmic-ray carrying outflows is however lacking. At the same time there is increasing awareness of the importance of massive star clusters in the acceleration of cosmic rays in galaxies. Here we report on the discovery of a nascent outflow driven by the massive star cluster Westerlund 1. Giga-electronvolt gamma-ray emission coincident with a cavity visible in atomic hydrogen traces the emergence of a population of relativistic electrons out of the Galactic Disc. The emission is offset from tera-electronvolt gamma-ray radiation surrounding the cluster, but connects to it smoothly spectrally and spatially. The implied energy density of co-accelerated protons and nuclei, assuming standard non-thermal electron/proton injection efficiencies, is at least an order of magnitude higher than that in the general interstellar medium. These particles therefore have the potential to dynamically influence the outflow. This discovery suggests that cosmic-ray loaded outflows may be a common feature of young massive star clusters, with implications for the transport of cosmic rays into the halo of the Galaxy.

In this paper I will describe the basic principles of chemical evolution of galaxies, its main ingredients and uncertainties. By means of chemical evolution we can perform the so-called galactic archaeology, which consists in reconstructing the history of star formation of galaxies and in particular of the Milky Way, starting from the observed stellar and gas abundances. Galactic archaeology is a powerful tool to predict also the behaviour of galaxies at high redshift. In particular, we adopt the "time-delay model" which is a way of interpreting the [X/Fe] vs. [Fe/H] relations (X is the abundance of a specific chemical element) in terms of different timescales of the stellar progenitors of the chemical elements relative to Fe, which is the indicator of "metallicity". I will then describe how we can reconstruct the star formation histories of galaxies of different morphological type (spirals, ellipticals) starting from the available observations. Particular attention will be paid to the study of the Milky Way which is the best studied galaxy at the moment. I will start describing the first chemical evolution models and then the most recent ones, the main difference between the old and new models being the observational data to compare with. Finally, I will foresee which could be the future improvements to chemical models and what constraints can we derive to better understand galaxy evolution.

We investigate convolutional neural network (CNN) methods for reconstructing the high-redshift density field from late-time large-scale structure, focusing on how the physical scale of the CNN input changes when standard first-order reconstruction is applied beforehand. Using dark-matter-only $N$-body simulations, we compare three approaches: a single-input CNN, a dual-input CNN combining two physical scales, and a single-input CNN applied to the density field after standard reconstruction. We vary the physical side length of the input sub-box over $L_\mathrm{sub}\sim38$-$380~h^{-1}\mathrm{Mpc}$ while keeping its numerical size fixed at $39^3$ voxels, allowing us to examine the trade-off between spatial context and resolution. For the CNN applied directly to the evolved density field, the reconstruction performs best at $L_\mathrm{sub}\sim150$-$200~h^{-1}\mathrm{Mpc}$. After standard reconstruction, however, the preferred scale shifts to $L_\mathrm{sub}\sim38$-$114~h^{-1}\mathrm{Mpc}$. The single-input CNN after standard reconstruction consistently outperforms both the single- and dual-input CNNs without standard reconstruction according to the normalized loss, density probability distribution, Kullback-Leibler divergence, residual field, and Fourier-space correlation. These results indicate that coherent large-scale displacements are more efficiently recovered by perturbative reconstruction, while the CNN is better suited to modelling the remaining quasi-linear and non-linear evolution on smaller scales. The preferred post-reconstruction input range includes the effective receptive scale of approximately $60~h^{-1}\mathrm{Mpc}$ adopted in previous hybrid reconstruction studies. Our findings therefore support a physically motivated separation of scales between analytic and data-driven reconstruction and demonstrate the advantage of combining the two approaches.

We use statistical tools to analyse data from the Solar Dynamics Observatory Helioseismic and Magnetic Imager to determine the distribution of the magnetic flux of photospheric magnetic features and its variation over a full solar cycle. In particular, we use statistical figures of merit to test how well different types of probability density function represent the magnetic flux distribution inferred from the data and how their shape changes over the solar cycle. Our analysis shows that a double power law provides the best representation of the data over the full solar cycle and we present the dependence of the power law exponents on the phase of the solar cycle. The nature of the observed flux distributions at different times during the solar cycle is significant because it could be used to try and infer information about solar magnetic field generation mechanisms. We discuss potential implications of a double power law distribution for solar magnetic field generation.

Next-generation radio interferometers will produce massive data volumes, making it impractical to store original visibility measurements and later combine observations in $uv$ spatial frequency space. Visibility measurements at similar $uv$ locations measure the same signal but different noise realizations. In principle, these measurements can therefore be compressed by storing only the inferred mean visibility and its uncertainty. We propose modeling the visibility with a sparse Gaussian process (GP) and storing the resulting compact probabilistic model rather than raw visibilities. Using simulated Atacama Large Millimeter/submillimeter Array (ALMA) observations, we demonstrate that the sparse GP is flexible enough to represent the visibilities and recover images with high fidelity. We estimate compression factors of $10^3-10^5$ for an 8-hour Square Kilometre Array (SKA)-Mid observation, with further gains expected by extending the GP input space to include the spectral axis. Beyond data compression, the model exploits correlations in $uv$ space, boosting the signal-to-noise ratio compared with independent grid averaging. Once trained, the model can predict visibility and its uncertainty at any desired $uv$ coordinates, allowing imaging with arbitrary fields of view and image resolutions. The model may also be incrementally updated with new observations while filtering outliers based on the prediction.

We present the 5th data release (DR5) of the Southern Photometric Local Universe Survey (S-PLUS), covering 4592 square degrees across 2491 fields. Observations were conducted with the T80-South, a Brazilian robotic telescope equipped with the Javalambre 12-filter system, containing five broad- and seven narrowband filters. Data products feature FITS images and extensive catalogs containing fluxes, magnitudes, and shape parameters for over 113 million detections. In addition, several value-added catalogs are provided, offering photometric redshifts (photo-zs), object classifications, masks, and extinction coefficients. For the first time, this release includes coverage of 110 square degrees along the Galactic disk, facilitating new research into Galactic structure and stellar populations. The release also provides full coverage of the Hydra Supercluster and numerous other nearby clusters with improved data reduction and calibration, enhancing photo-z accuracy, which is vital for large-scale structure studies. A preliminary analysis of the Hydra and Antlia galaxy clusters up to $5 \times R_{200}$ yields an updated catalog of 1706 cluster members based on both spectroscopic and high-quality photo-zs. Our photometric data shows that both clusters have a similar proportion of galaxies with an H$\alpha$ excess relative to their clustercentric distance, though Hydra has a higher fraction near its center. Additionally, the spatial distribution of all objects in our sample highlights a bridge connecting both clusters. We verify that S-PLUS DR5 provides a solid foundation for future scientific investigations, ranging from Solar System studies to Cosmology.

With nearly 400 Gravitational Wave (GW) events detected by the LIGO-Virgo-KAGRA network and many more expected, similarities between signals produced by different astrophysical effects can complicate template-based searches and parameter estimation. In particular, GW modulations from spin precession can resemble the beating pattern induced by microlensing from compact objects with masses of $10$-$10^5,M_\odot$. We investigate this degeneracy and assess whether machine learning can distinguish between these effects. We generate 20,000 simulated GW signals for each class with a network optimal signal-to-noise ratio above 20 and train a convolutional neural network on time-frequency (Q-transform) spectrograms. The classifier achieves up to 95% accuracy in Gaussian noise and 82% in real detector noise. We also study classification between microlensed (ML) and unlensed non-spinning (UN) signals, as well as between unlensed non-spinning (UN) and unlensed precessing (UP) signals. While distinguishing UN from UP remains difficult even in Gaussian noise, ML vs. UN classification reaches up to 80% accuracy in real noise. We identify the regions of parameter space where the classifier performs best and evaluate the ML-UN network on real GW events, finding that the model trained on Gaussian noise generalizes better than the one trained on real noise. This work presents the first low-latency machine-learning pipeline for distinguishing microlensed from unlensed non-spinning GW signals.

Constraining the total mass of protoplanetary disks is crucial to determine the availability of material for planet formation. Yet, providing accurate and precise measurements of the disk mass is challenging. Investigating the gas dynamics is a powerful, tracer-independent method to precisely characterize disk masses. By fitting the velocity rotation curves of different molecular tracers with an accurate model including the disk thermal stratification and self-gravity, we constrain the stellar masses, disk masses, and scale radii for the disks around HD 97048 and WaOph 6. We obtain $M_\star=2.226 ^{+0.054}_{-0.049}\ M_\odot$, $M_\mathrm{d}=0.3 ^{+0.055}_{-0.061}\ M_\odot$ and $R_\mathrm{c}=172 ^{+24}_{-14}$ au for HD 97048, and $M_\star=0.956\ ^{+0.006}_{-0.006}\ M_\odot$, $M_\mathrm{d}=0.21 ^{+0.045}_{-0.038}\ M_\odot$ and $R_\mathrm{c}=647 ^{+193}_{-155}$ au for WaOph 6. We also measure the corresponding gas-to-dust and disk-to-star mass ratios. We efficiently extend the dynamical method to characterize embedded sources exhibiting features of absorption, for which a careful analysis is required to avoid biases in the retrieved velocity profiles. We prove the importance of including a beam smearing correction to the curves: if not, this observational effect can systematically bias the velocity profiles, altering the disk mass estimates up to $\sim45\%$. We provide comprehensive estimates of the systematic uncertainties on the best-fit parameters by bootstrapping over both the retrieved geometry and 2D thermal structure of the two disks: the overall uncertainty on the disk masses is $\sim20\%$. Finally, we investigate the connection between disk stability and the appearance of spiral morphologies in the mm continuum emission, by comparing the Toomre parameter of all dynamically weighed disks to date, showing that disks with mm-dust spirals have systematically lower values of Q.

Visible spectra of the Enceladus particle plume and E ring contain evidence for a change in spectral slope around 0.5 micron. This feature can be seen in data obtained by both the Visual and Infrared Mapping Spectrometer (VIMS) and Imaging Science Subsystem (ISS) onboard the Cassini Spacecraft, and is consistent with the slope change seen in the surface spectra of Saturn's rings and moons that has been attributed to either organics or iron compounds. The observed spectral features in the plume and E ring could represent either a non-ice contaminant in the plume particles or a deficit of sub-micron particles, so this spectral feature provides a new tool for assessing variations in the plume particle's composition and/or size distribution with time and space. The observed strength of this feature is consistent with the plume particles having an organic fraction similar to that measured by in-situ measurements, so there are good reasons to expect that this feature can be used to quantify the organic content of the plume particles. There are also hints of a potential absorption band around 0.45 micron in these spectra. If this feature can be confirmed, it could provide further constraints on the plume particles' composition.

Magnetars are neutron stars with ultra-strong magnetic fields ($B \sim 10^{14}$-$10^{15}$ G) and are promising candidates for high-energy particle acceleration. We present a multiwavelength analysis of the region surrounding CXOU J171405.7-381031, a magnetar associated with the supernova remnant (SNR) CTB 37B. The broadband spectral energy distribution spanning radio to TeV energies is modeled using leptonic and lepto-hadronic scenarios, with particle populations constrained using Markov Chain Monte Carlo techniques. Within an SNR framework, both scenarios provide acceptable descriptions of the gamma-ray data. The purely leptonic model reproduces the overall spectral shape but slightly underestimates the highest-energy flux measured by H.E.S.S., whereas a lepto-hadronic interpretation offers an improved description above $\sim 10$ TeV, with inverse-Compton scattering dominating the GeV emission and neutral-pion decay contributing at the highest energies. The required proton energy ($W_{\rm p} \gtrsim 10^{51}$ erg) can be substantially reduced if the remnant interacts with a dense ambient medium. A magnetar wind nebula scenario can reproduce the broadband spectrum but is strongly disfavored by the observed source morphology. Simulated Cherenkov Telescope Array Observatory (CTAO) observations indicate that exposures of $\sim 50$ h will constrain the proton cut-off energy, enabling a decisive test of hadronic emission in this region.

Lya emitters (LAEs) are galaxies with strong Lya emission, tracing early star formation and ionizing radiation. Their connection to active galactic nuclei (AGNs) is key to understanding the mechanisms behind (extended) Lya emission. In this work, we measure the fraction of LAEs identified as radio-emitting AGN (fAGN,radio) and the fraction of radio sources that exhibit Lya emission (fLya) to investigate the connection between radio AGN activity and Lya emission at 1.88 < z < 3.52. We identify 928 sources detected in both the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the LOw Frequency ARray (LOFAR) surveys. These matches are drawn from 55,109 spectroscopically confirmed LAEs and 27,625 radio sources. After applying completeness corrections, we obtain fAGN,radio = 1.77 $\pm$ 0.04% and fLya = 18.15 $\pm$ 0.14%. The fraction fAGN,radio increases from 0.4 $\pm$ 0.1% to 9.7 $\pm$ 1.3% with increasing Lya luminosity, while fLya rises from 0.7 $\pm$ 0.1% to 55.8 $\pm$ 14.5% with radio luminosity. `LAEs with radio AGN' and `optical AGN with Lya emission' show similar radio luminosities above the AGN threshold, although optical AGN have higher Lya luminosities. We find no significant correlation between Lya luminosity and either radio luminosity or spectral index. Lya line width increases with Lya luminosity but shows no correlation with radio size. Our results show that most Lya emission at 1.88 < z < 3.52 is powered by star formation, with radio AGN activity confined to a small luminous subset (1.77 $\pm$ 0.04%). The absence of correlations between Lya and radio properties suggests that Lya emission is governed primarily by host-galaxy gas properties rather than direct AGN jet coupling.

Coronagraphic imaging of exoplanets is limited by residual speckles that mimic planets. Advanced post-processing is essential for current and future instruments on the ground or in space. Current techniques are time-intensive and limited. ADI requires long sequences and is limited at small separations. RDI is also time-consuming and sensitive to speckle evolution, leading to imperfect subtraction. Coherence Differential Imaging (CDI), which we successfully demonstrated on SPHERE, offers a faster alternative by using the light incoherence between speckles and planets. However, its reliance on accurate instrumental models limits its performance. In this work, we introduce EPICX, an enhanced CDI method using gradient-boosted decision trees. EPICX models the differential signal as a high-dimensional regression problem, optimizing the discrimination between coherent speckle noise and incoherent planet signal. We validate this enhanced CDI method using simulated data for different coronagraphs, including those aboard JWST and Roman.

Observations of the aromatic infrared emission band at 3.3 $\mu$m often reveal satellite emission features in the 3.4 - 3.6 $\mu$m range. While the 3.3 $\mu$m band is attributed to the CH stretching vibration of polycylic aromatic hydrocarbons (PAHs), the satellite bands - particularly its prominent 3.4 $\mu$m component - is assigned to aliphatic CH stretching vibrations in hydrogenated and methylated PAH-like species. Our aim is to derive state-of-the-art infrared emission spectra for aliphatic-containing pyrene derivatives and compare them with astronomical observations. This will help refine our understanding of the contribution of these species to the 3.4 $\mu$m emission band. Mid-infrared spectra (1.4-25 $\mu$m) of gas-phase dihydropyrene, tehtrahydropyrene, hexahydropyrene methylated pyrene, and pyrene were recorded at temperatures ranging from 373 to 673 K, depending on the species. The band profiles were analyzed using a multi-component fitting tool, and empirical anharmonicity laws were derived to quantify the evolution of the band positions and widths with temperature. The obtained spectral data was combined with the results of a Monte Carlo emission model to simulate the emission spectra following UV-photon absorption, up to the dissociation limit ($\lesssim$6 eV). The resulting synthetic spectra were compared with James Webb Space Telescope observations of the Orion Bar region (PDRs4All program). Based on these state-of-the-art simulated spectra, we propose 1,2,3,6,7,8-hexahydropyrene as the carrier of the red component of the 3.4 $\mu$m band observed at 3.403 $\mu$m. While 1-methylpyrene may also contribute to the underlying emission plateau, its lack of a strong infrared band complicates detection in observed spectra, unlike hexahydropyrene. All experimental spectra and their temperature-dependent analyses are available in the new cosmicPAH-IRDB database.

Cosmology is an observational precision science, and detailed and scientifically rich simulations are necessary components of many analyses. These simulations are often computationally expensive and produce massive amounts of complex data. Sharing that data widely to enable further explorations, comparison with observations, and communication with general audiences is crucial to advancing the scientific frontier and engaging a broader community. In this paper, we introduce the OpenCosmo project, designed to provide flexible access and analysis modes to large cosmological simulation datasets. We make available an easy-to-use web-based data portal for retrieving downloadable subsets of flagship-scale cosmological datasets, and a sophisticated analysis library for further analysis on the returned data. OpenCosmo supports multiple levels of interaction -- from simple search and download to interactive exploration and computation -- by integrating with existing HPC and data infrastructure. Its architecture, based on Globus Compute, provides a scalable and adaptable framework that can be extended to other scientific domains seeking to couple data sharing with computational capability.

The polarized Sunyaev Zel'dovich (pSZ) effect is a cosmic microwave background (CMB) polarization anisotropy induced by Thomson scattering from free-electrons in non-linear structure. The pSZ signal is determined by the distribution of ionized gas tracing the cosmic web and the CMB quadrupole at the location of free-electrons - the remote quadrupole field. Measuring the pSZ effect provides a consistency check of the optical depth to reionization and sheds light on the anomalous nature of the large-scale CMB temperature anisotropies, such as the low observed CMB temperature quadrupole. In this paper, we demonstrate that a CMB-CMB-galaxy bispectrum summarizes several existing pSZ statistics, and that in our observable Universe the ideal galaxy sample to detect pSZ is at $z \sim 1-2$. We evaluate the bispectrum using CMB data from Planck and ACT with galaxy density from the unWISE galaxy redshift catalog as well as Planck cosmic infrared background (CIB) maps. We do not make a statistically significant detection of the pSZ effect, which is consistent with the expected O$(1)$ signal-to-noise from this data combination. The measured amplitude of the pSZ bispectrum provides constraints on the optical depth bias associated with large-scale structure (the amplitude of the pSZ signal) of $b_q=1.02 \pm 2.64$, the optical depth to reionization of $\tau_{\rm rei} = -0.01 \pm 0.14$, and the tensor-to-scalar ratio $r$ of $\sigma_r \sim 150$ ($n_t = 0$) or $\sigma_r \sim 3$ ($n_t = -1$). We forecast that future measurements could tighten the constraints on these quantities by roughly a factor of 3, which is sufficient to provide independent confirmation of the low CMB quadrupole and the optical depth to reionization.

The 2021 Event Horizon Telescope observations resolve the innermost jet region of the blazar 3C279 with unprecedented detail. The reconstructed images consistently reveal a compact core elongated nearly orthogonal to the large-scale jet axis. This rarely observed morphology recurs across multiple epochs and from 22-230 GHz and is therefore intrinsic rather than an imaging artifact. Geometric model fitting identifies several components with apparent speeds up to 10c, requiring bulk Lorentz factors greater than 10.3 and constraining viewing angles to extremely small values (smaller than one degree). Rest-frame brightness temperatures are systematically low (between 10^9 and 10^10 K), consistent with optically thin emission at 230 GHz. These results suggest that the jet bends toward the observer on sub-parsec scales, producing strong relativistic beaming. Possible drivers of the observed jet bending and temporal evolution include the jet's interaction with the interstellar medium, kink or Kelvin--Helmholtz instabilities, magnetic reconnection near the horizon, or binary-induced precession. However, the current temporal coverage of VLBI data remains insufficient to distinguish between these mechanisms. Continued multifrequency VLBI monitoring will be essential to constraining the dynamics and geometry of the jet base in 3C279.

At the South African Astronomical Observatory (SAAO), we have developed FiberPol$-$ a spectropolarimetric front-end for the SpUpNIC spectrograph on the 1.9 m telescope fed by fibers. In conjunction with SpUpNIC, it combines two niche techniques: fiber-spectroscopy and polarimetry to enable a wide range of science cases, particularly for studies of the interstellar medium, and serves as a technology pathfinder for implementing polarimetry larger facilities such as the 10 m Southern African Large Telescope. Commissioned in early 2025, FiberPol has already demonstrated key performance benchmarks. Based on a preliminary data reduction pipeline, we find that the uncertainties in the measured Stokes parameters $q$ and $u$, and the degree of polarization $p$ are in the range of 0.2% to 0.3% in 5 nm spectral bins across the 400 to 700 nm wavelength range. With improved calibration and refined analysis, we expect to reach our design goal of 0.1% accuracy per bin. Several science programs are already underway following its successful commissioning. FiberPol employs a Wollaston prism and half wave plate system to analyze linear polarization and extract spectropolarimetric data comprising wavelength ($\lambda$) and the Stokes parameters ($I$, $q$ and $u$). Its compact and modular design uses primarily small, commercially available optical and optomechanical components, making it straightforward to adapt to other spectrographs, including those on larger telescopes. This paper presents the design, laboratory validation, and commissioning results of FiberPol, and outlines planned upgrades to further enhance its scientific impact.

Motivated by the recent discovery of PSR~J0435+3233, a millisecond pulsar with an exceptionally large period derivative of $\dot{P}=4.9\times10^{-17}\ {\rm s\,s^{-1}}$ and a high spin-down luminosity of $\dot{E}=5.89\times10^{37}\ {\rm erg\,s^{-1}}$, we analyze $\sim$17~yr of observations obtained with the \textit{Fermi} Large Area Telescope for the pulsar. We identify the cataloged source 4FGL~J0435.5+3232, located only $0.01^{\circ}$ from the radio timing position, as the $\gamma$-ray counterpart of PSR~J0435+3233. Using $\gamma$-ray events collected within the validity interval of the radio timing ephemeris, we detect its $\gamma$-ray pulsations at a $\sim6.8\sigma$ confidence level. In this interval, the phase-resolved analysis shows that the pulsed emission is concentrated predominantly within the rotational-phase of $\phi\sim0.44$--$0.69$. We also derive its $\gamma$-ray luminosity of $L_{\gamma}=6.26\times10^{32}\ {\rm erg\,s^{-1}}$, assuming a distance of 1.2~kpc and isotropic emission. This luminosity corresponds to an apparent $\gamma$-ray efficiency of $\eta_{\gamma}\sim1.1\times10^{-5}$, revealing an exceptionally low $\gamma$-ray output despite the pulsar's young-pulsar-like rotational-energy budget. Our detection establishes PSR~J0435+3233 as a $\gamma$-ray MSP, and the striking combination of its high spin-down power and low apparent $\gamma$-ray efficiency provides a new probe of particle acceleration, radiation beaming, and viewing geometry in the magnetospheres of millisecond pulsars with extreme rotational properties.

Detections of gravitational waves from neutron star-black hole systems provide avenues for studying extreme matter, constraining binary formation channels, and testing the nature of compact objects in strong gravity. Eccentric signatures in the signal further enhance this potential by improving parameter estimation and offering clues about binary formation. Because eccentricity is primarily imprinted during the inspiral phase, it is often weakly constrained or missed entirely in binary black hole observations; in contrast, neutron star-black hole systems produce longer in-band signals, enabling more precise measurements of eccentricity and leaving a distinct imprint on parameter inference. In this work, we present a systematic parameter-estimation study exploring the impact of eccentricity on inference using injections simulated with the state-of-the-art eccentric waveform model SEOBNRv5EHM. We find that for systems like GW200105_162426, the measurement precision of eccentricity and correlated parameters improves as eccentricity increases, yielding tighter constraints at larger eccentricities. For the highest eccentricity considered in this study, $e=0.25$, we recover eccentricity with $1\sigma$ uncertainty as low as $4\times10^{-4}$. In addition, the constraints on effective spin $\chi_\mathrm{eff}$ and mass ratio $q$ improve relative to the quasi-circular case by factors of $\sim13$ and $\sim20$, respectively. On the other hand, we find no significant improvement in extrinsic parameters such as luminosity distance and sky localization, suggesting that for systems like GW200105_162426, the additional information provided by eccentricity in this sector is either negligible or degenerate with the information provided by higher-order modes.

The formation history of compact binary systems remains one of the key open questions in astrophysics. Theoretical studies generally favor isolated binary evolution for neutron star-black hole (NSBH) systems, which tends to produce nearly circular orbits. However, recent analyses of the gravitational-wave event GW200105 indicate that its source has measurable eccentricity, suggesting that alternative formation channels may also contribute. It has been shown that the intrinsic parameters of eccentric NSBH mergers, such as the component masses and spins, are much better measured than circular mergers with LIGO-Virgo-KAGRA (LVK) observatories. We explore how such eccentricity-enhanced parameter measurements can affect the inference of NSBH formation channels. We find that sharper measurements of the effective spin parameter $\chi_{\rm eff}$ increase the fraction of systems for which negative values can be confidently identified, allowing for the clear measurement of a spin-orbit misaligned event every $\sim 2.5$ eccentric NSBH detections for an isotropically distributed population in the fifth LVK observing run (O5). Improved NS mass measurements provide better constraining power for NS mass distributions, potentially revealing structure and tightening the bounds that can be drawn on their upper and lower masses. Similarly, the recovery of a metallicity-dependent BH mass distribution is improved by eccentricity-enhanced measurements. Finally, we show that the proposed population-level eccentricity distribution for dynamical-formation channels can be tested by the end of O5.

IXPE has provided for the first time detailed energy- and time-resolved X-ray polarimetry of Z-type neutron star low-mass X-ray binaries (NS-LMXBs) as they move along their color-color diagrams (CCDs). These sources can reach the highest polarization observed for NS-LXMBs in the 2-8 keV range when they move along the horizontal branch. In a previous paper, we characterized the spectral state of a sample of Z-sources using the CCD and estimated the polarization with model-independent analysis. Here, we present detailed spectropolarimetric analysis for each source on each branch using data from IXPE, NICER, and NuSTAR. The continuum X-ray emission of all the sources is well described with a combination of thermal accretion disk emission plus a harder Comptonized component. In addition, reflection features, in particular the relativistically broadened Fe line, are observed for our sources, except GX 5-1. For most of the sources and branches, the main contribution to the X-ray emission and polarization is due to Comptonization: moving from the horizontal branch (HB) to the normal branch (NB), the polarization degree (PD) in the 2-8 keV band varies from about 6% to 3-4%, while the PD is loosely constrained in the flaring branch (FB), due to the shorter exposures. These PD values are significantly higher than theoretical expectations for typical spreading or boundary layer configurations. The polarization of the disk is generally lower (below 3%) but still higher than predictions for an electron scattering-dominated, plane-parallel atmosphere above the disk observed at the corresponding inclination. Moreover, the polarization angle (PA) of the disk seems to be significantly misaligned and not perpendicular to that of Comptonization. We find no correlation between the polarization signal and the inclination, nor with the contribution of reflected photons throughout the Z-track.

The excursion-set formalism provides a key connection between primordial density fluctuations and the abundance of cosmic structures such as dark matter halos and voids, traditionally assuming Gaussian random walks. In this work, we extend this framework to fluctuations whose distribution presents strongly non-Gaussian tails. Such tails are beyond the reach of perturbative approaches to primordial non-Gaussianity based on moment expansion. We address the problem with rigorous, analytical derivations relying on the large deviation principle, suited for the study of rare fluctuations. We derive new first-passage time distributions for random walks with non-Gaussian statistics and obtain updated predictions for the halo mass function. We also study the two-barrier problem relevant to cosmic void formation, leading to a new analytical prediction for the void size function, with improved accuracy on large scales. Our results demonstrate the potential of large deviation techniques as a bridge between inflationary scenarios, often leading to strongly non-Gaussian tails, and late-Universe observables.

Solar energetic particle (SEP) event forecasting often involves integrating heterogeneous observations that differ in cadence, temporal coverage, format, and historical availability, posing challenges for reproducible analysis of data-driven approaches. This paper presents SEP-PRISM Data, a curated multi-source dataset designed for 24-hour ahead forecasting of operational SEP events, defined by proton flux exceeding 10 pfu in the GOES > 10 MeV channel. SEP-PRISM Data integrates flare records, active-region magnetic field parameters, coronal mass ejection (CME) catalogue data, GOES soft X-ray flux, and historical proton flux into a common window-based representation spanning 3 February 1986 to 10 September 2025. To improve temporal coverage and cross-source consistency, SHARP and SMARP magnetic products were aligned into a unified SMHARP archive, and CME records from DONKI and CDAW were aligned into a unified CDAWDONKI event set. Predictor variables were summarized over fixed non-overlapping 24-hour historical windows using minimum, mean, and maximum statistics and paired with targets defined over the subsequent 24-hour window, forming a supervised learning dataset. The resulting SEP-PRISM Data contains 14,464 labeled samples, including 650 positive operational SEP cases, and is intended to support reproducible benchmarking, model development, feature analysis, and future studies of space weather forecasting.

In this work, we investigate the prospects for probing primordial non-Gaussianity and associated symmetry of parity through scalar-induced gravitational waves (SIGWs), with third-generation gravitational-wave detectors. We develop a framework that accounts for contributions to the energy density spectrum of GWs arising from the scalar non-Gaussianity quantified by bispectrum and trispectrum, and perform parameter inference using simulated data from Einstein Telescope and Cosmic Explorer. The parity-odd component of the scalar trispectrum induces circular polarization in the stochastic gravitational-wave background (SGWB), providing a direct probe of parity-violation in the primordial Universe. We show that future interferometers can place competitive constraints on the parity-odd scalar trispectrum, along with the bispectrum and parity-even trispectrum. Moreover, we include the astrophysical contribution, which could act as a foreground for the SIGWs. We show that, despite the addition of such a foreground, we are still able to effectively constrain the cosmological parameters related to SIGWs and the astrophysical parameters as well.

Surface brightness fluctuations (SBF), the pixel-to-pixel variations in flux that arise from the statistical distribution of stars in galaxy images, provide a powerful tool to measure redshift-independent distances from photometric data alone. This method is particularly important in the era of large imaging surveys, such as those carried out during the Euclid mission. Here we present the first application of SBF to Euclid data, using the FAST-SBF code to measure stellar fluctuation amplitudes in the $I_\mathrm{E}$ band for a sample of galaxies in the Early Release Observations (ERO) of the Fornax galaxy cluster. Although the Euclid data reduction pipeline is not optimized for SBF measurements, extensive testing suggests that we are able to extract robust results. We calibrate the absolute fluctuation magnitude $\overline{M}_\mathrm{IE}$ as a function of the $(I_\mathrm{E}{-}H_\mathrm{E})$ colour for 15 galaxies in Fornax, then test this relation on galaxies in two other ERO fields: the Perseus cluster and dwarf satellite candidates around NGC 6744. Overall, we find reasonable agreement with distances in the literature and good agreement in cases where FAST-SBF indicates the measurements are robust. Finally, we compare our results against Stellar Population Tools (SPoT) simple stellar population models, and discuss the possibility of using SBF with Euclid to probe the underlying stellar populations of the galaxies.

The optical morphologies of low-mass galaxies can be used to directly trace their assembly and constrain models of galaxy evolution. We select a sample of 6211 low-mass ($7\lesssim {\rm log}(M_{\ast}/M_{\odot})\lesssim 10$) star-forming galaxies from the SAGAbg catalog that has high-completeness at low-redshifts ($z<0.1$). We obtain their galaxy maps in the $griz$ - bands of the Legacy Surveys and apply STATMORPH to calculate non-parametric morphological measures including the Gini index, $M_{20}$ measure, and CAS parameters. We study how resolution and signal-to-noise affect the morphology measures and find that the bulge strength measurements are the most reliable. The sequence in Gini$-M_{20}$ space is directly linked to the star-forming sequence of galaxies and is dominated by Sb/Sc/Ir morphologies. The $g$-band light distributions are the least concentrated among all the bands. The systematic trends of $M_{20}$ with respect to stellar mass and GALEX NUV-derived specific star formation rate (sSFR) strongly indicate that the galaxies with flatter light profiles are less massive with higher sSFR and vice-versa. We statistically infer that star-forming low-mass galaxies predominantly have disk morphologies with bulges becoming more prominent in quenched systems at higher masses (${\rm log}(M_{\ast}/M_{\odot})\gtrsim 9$).

Interactions between dark matter and dark energy can lead to an apparent phantom-crossing behavior that mimics the expansion history preferred by the latest cosmological observations from DESI baryon acoustic oscillations (BAO), Cosmic Microwave Background (CMB), and Type Ia supernovae (SNe Ia) data. In a previous paper [Khoury, Lin, and Trodden 2025 arXiv:2503.16415], we proposed a concrete particle physics realization of this idea, consisting of a strongly coupled dark sector in which a dark axion is coupled to dark baryons. In this paper, we investigate this idea further by comparing its predictions to the latest cosmological data. We implement the dark axion-dark baryon interaction model in a Boltzmann code and confront it with CMB, DESI DR2 BAO, and SNe Ia data. For the CMB+DESI DR2+DES-Dovekie combination, the best-fit model improves the fit relative to $\Lambda$CDM by $\Delta\chi^2=-14.48$. The preferred solution exhibits a non-monotonic dark-matter mass evolution: the mass decreases between matter-radiation equality and recombination, while increasing over the BAO/SNe-sensitive epoch, leading to an apparent phantom crossing in an effective dark-energy description. Interestingly, the same dynamics produces an Early Dark Energy-like energy injection near matter-radiation equality, but in the data-preferred region this component is too small to raise $H_0$ enough to substantially reduce the current tension.

Direct evidence for collective neutrino oscillations in core-collapse supernovae remains elusive. We show that this quantum phenomenon leaves a footprint on the abundance pattern of proton-rich nuclides in the solar system. Modeling the $\nu p$-process using a $20\,M_\odot$ progenitor, we map out the dependence of the total yields on the starting radius of the oscillations, self-consistently coupling hydrodynamics and nucleosynthesis. The oscillations boost key $p$-nuclides ($^{92,94}\text{Mo}$, $^{96,98}\text{Ru}$) and long-lived $^{92}\text{Nb}$ by up to two orders of magnitude, bringing their abundances into agreement with the observations. The best match is found when oscillations commence within $10\text{ km}$ of the proto-neutron star surface, indicating fast collective oscillations.

The collapse of false-vacuum domains during first-order phase transitions in the early Universe may lead to primordial black hole (PBH) formation whose signatures form a multimessenger complement to gravitational wave (GW) production. We focus on PBH formation through the gravitational collapse of false-vacuum domains, described using a junction condition formalism. This formalism develops the Schwarzschild collapse criterion dynamically, avoiding the usage of critical overdensity thresholds in a post-inflationary Universe, and is driven solely by the vacuum energy enclosed within shrinking false-vacuum domains without the assistance of particle or domain wall interactions in the false vacuum. We study the parameter space of phase transitions and identify regions producing observable GWs, observable PBHs, or both simultaneously. We investigate this phenomenology in polynomial and classically conformal scalar field potentials as model benchmarks. We find that the scalar fields with vacuum expectation values in the range of 1-100 MeV have the largest model parameter space available for these multi-messenger signals, which are testable with upcoming GW observatories, searches for Hawking radiation, and gravitational lensing surveys.

We investigate the relativistic scaling and electromagnetic feedback of Stern--Gerlach (SG) force driven spin transport in pair-plasma reconnection with SpinPIC2D. The model advances relativistic proper momentum and magnetic BMT spin precession, applies the SG force, deposits spin magnetization, and includes $\mathbf J_M=\nabla\times\mathbf M$ in Ampere's law. In a weak-seed scan at fixed $\gamma_{\rm tr}=2$, the normalized global magnetic flux-growth remains near the classical control-run for $\Xi\le0.1$, is $0.016$ at $\Xi=0.4$, and reaches approximately $0.12$ and $0.25$ at $\Xi=0.7$ and 1, respectively over $3\le t/\tau_{\rm sp}\le7$. Because $\partial/\partial y=0$ in the 2.5-D geometry, the direct $y$-directed SG term vanishes and the enhancement is indirect: sheet-normal SG sorting restructures branch-resolved electron and positron velocity distributions which results in changes in pressure moments and generates a layered magnetization current. Along a fixed-$\chi_{\rm sim}$ family, increasing $\gamma_{\rm tr}$ reduces $\Xi$ and suppresses branch sorting, whereas the matched-control flux-growth enhancement remains positive for $\gamma_{\rm tr}=2,3,5$. Retuning $\chi_{\rm sim}$ to hold $\Xi=1$ does not preserve nonlinear similarity: both the coupling and $J_{M,y}$ increase with $\gamma_{\rm tr}$, and the $\gamma_{\rm tr}=5$ case develops a multi-X-line state. Thus $\Xi$ orders the onset of SG-modified reconnection, while the nonlinear response also depends on the absolute spin coupling and magnetization-current amplitude.

Unlike the gravitational waves classically generated by moving macroscopic masses, the diffuse backgrounds of gravitons originate most likely from zero-point fluctuations of the gravitational field amplified by the evolution of the space-time curvature. The resulting entangled states lead to specific second-order correlation effects that could be eventually detected. For a quantitative analysis of this empirical expectation we scrutinize the interactions between the cosmic gravitons and the fundamental mode of a quantized electromagnetic field confined inside a closed optical resonator with perfectly-reflecting walls. We show that the Hanbury-Brown Twiss correlations of the photons are insensitive to the degrees of second-order coherence of the gravitons even barring for the exceedingly small couplings of the problem. Since the degree of second-order coherence of the photons does not reflect the correlation properties of the gravitons, the statistical properties of the gravitons (and their super-Poissonian statistics) cannot be inferred, even in principle, from the intensity correlations of the cavity modes.

We investigate the dynamics of vector excitations localized on cosmic strings in the Abelian-Higgs model. These confined gauge field fluctuations behave as massive vector degrees of freedom propagating along the string. We show that they remain bounded for all values of the scalar self coupling $\lambda$, including parameter regimes in which the previously studied conventional shape mode is no longer present. After determining their spectra and spatial profiles, we analyze their decay through nonlinear coupling to bulk radiation. Their amplitudes exhibit the characteristic power law relaxation associated with non-linear radiation. For sufficiently large $\lambda$ the scalar radiation channel becomes kinematically inaccessible and the decay is instead mediated by quasinormal modes which can be interpreted as Feshbach resonances of the string. We then study the full $3+1$ dimensional dynamics of the vector-excited string. Field theory simulations reveal a parametric instability between the vector mode and the Goldstone sector. In contrast to the shape mode, the vector mode resonantly excites both transverse directions simultaneously. In addition, we build an effective model that captures the instability and identifies the nonlinear couplings responsible for it. Our results show that vector excitations constitute a long-lived massive degree of freedom with a non-trivial role in the dynamics of local strings.

When a star is described as a spectral class G2V, we know its approximate mass, temperature, age, and size. At more than 5,700 exoplanets discovered, it is a natural developmental step to establish a classification for them, such as for example, the Harvard classification for stars. This exoplanet classification has to be easily interpreted and present the most relevant information about them and divides them into groups based on certain characteristics. We propose an exoplanet classification, which using an easily readable code, may inform you about a exoplanet's main characteristics. The suggested classification code contains four parameters by which we can quickly determine the range of temperature, mass, density and their eccentricity. The first parameter concerns the mass of an exoplanet in the form of the units of the mass of other known planets, where e.g. M represents the mass of Mercury, E that of Earth, N Neptune, or J Jupiter. The second parameter is the mean Dyson temperature of the extoplanet's orbit, for which we established four main classes: F represents the Frozen class, W the Water class, G the Gaseous class, and R the Roaster class. The third parameter is eccentricity and the fourth parameter is surface attribute which is defined as the bulk density of the exoplanet, where g represents a gaseous planet, w - water planet, t - terrestrial planet, i - iron planet and s - super dense planet. The classification code for Venus, could be EG0t (E - mass in the range of the mass of the Earth, G - Gaseous class, temperature in the range from 450 to 1000 K, 0 - circular or nearly circular orbit, t - terrestrial surface), for Earth it could be EW0t (W - Water class - a possible Habitable zone). This classification is very helpful in, for example, quickly delimiting if a planet can be found in the Habitable zone; if it is terrestrial or not.

In this study, we investigate the optical appearance of rotating Konoplya-Zhidenko non-Kerr black holes in the presence of thin accretion disks, with the aim of examining whether the information of deformation parameters manifest in observable signatures. By employing a fisheye camera model in conjunction with backward ray-tracing techniques, we simulate images for both prograde and retrograde accretion scenarios. The results indicate that the deformation parameter $\xi$ can partially mitigate the shadow deformation induced by the rotation parameter $a$. The inner shadow displays characteristic morphological transformations at varying observation angles, transitioning from an axisymmetric circular form at low angles to a hat-like configuration at higher angles. Furthermore, at high observational inclination angles, the direct image and the lensed image become distinctly discernible, and an increase in the deformation parameter $\xi$ enhances the observed intensity of the image. Interestingly, the motion behavior of the accretion flow influences the observed intensity distribution on the screen, a finding that is consistent with the redshift distribution. Hence, variations in the deformation parameter $\xi$, the observation angles, and the motion behavior of the accretion flow collectively influence the observable appearance of the black hole. We expect this work to provide valuable references for identifying observable signatures of spacetime deviations from general relativity.

Recently discovered asteroid 2024 YR$_4$ has an orbital period of almost exactly 4 years, a descending node located almost exactly 1 au from the Sun, with a perihelion slightly less than that. Its combination of semi-major axis and perihelion renders it a 'Potentially Hazardous Object' (PHO) in the Apollo class. It now has a low chance of colliding with the Earth on 22 December 2032, yet there is the potential for many further close encounters with Earth into the distant future. This paper investigates the feasibility of missions to this object in the short term, up to and including its close encounter in 2032, and exploits the preliminary mission design software known as 'Optimum Interplanetary Trajectory Software' (OITS). Many flyby opportunities are found with launch windows virtually throughout 2028. Sample Returns are also eminently feasible over this period. Rendezvous missions with 'New Horizons' spacecraft (adopted as a convenient reference mission) are available, although these require launches around late 2028 and early 2029, and with much longer flight durations. In summary it is found that 2024 YR$_4$ represents an 'opportunity rich environment' and bodes well for any future attempts by humanity to either examine this object close up, or even to deflect it if necessary.

Reconstructing the initial density field of the Universe from the late-time matter distribution is a nontrivial task with implications for understanding structure formation in cosmology, offering insights into early Universe conditions. Convolutional neural networks (CNNs) have shown promise in tackling this problem by learning the complex mapping from nonlinear evolved fields back to initial conditions. Here we investigate the effect of varying input sub-box size in single-input CNNs. We find that intermediate scales ($L_\mathrm{sub} \sim 152\,h^{-1}\,\mathrm{Mpc}$) strike the best balance between capturing local detail and global context, yielding the lowest validation loss and most accurate recovery across multiple statistical metrics. We then propose a dual-input model that combines two sub-boxes of different sizes from the same simulation volume. This model significantly improves reconstruction performance, especially on small scales over the best single-input case, despite utilizing the same parent simulation box. This demonstrates the advantage of explicitly incorporating multi-scale context into the network. Our results highlight the importance of input scale and network design in reconstruction tasks. The dual-input approach represents a simple yet powerful enhancement that leverages fixed input information more efficiently, paving the way for more accurate cosmological inference from large-scale structure surveys.

Periodic gravitational forcing by the Moon and Sun produces small oscillations in Earth's rotation known as nutations. Nutations are amplified by a resonance with a natural motion of the liquid core called the Free Core Nutation, whose amplitude is limited by friction-like processes at the core--mantle boundary. Previous studies have attributed this damping to the dissipation of electric currents induced in the lower conducting mantle, but, given current knowledge of the lower mantle, electromagnetic coupling appears insufficient to fully account for the observed lag. We show that additional dissipation arises from the interaction of the tidal flow inside the core with the topography of the core--mantle boundary, which excites internal waves that extract energy and momentum from the flow. Adapting a theory originally developed for tides over seafloor topography, we find that the observed damping can be fully accounted for by a topography of typical amplitude $\sim$5~km dominated by features of wavelength $\sim$1500~km. The dissipation is highest when the upper core is neutrally buoyant. Such amplitudes are larger than typical inferences from global seismic studies but are not ruled out, given regional seismic evidence for kilometer-scale features and the sparse constraints at these scales.

Fast Radio Bursts (FRBs) have become an indispensable tool for studying the ``missing baryons'', the Universe's ionisation properties, as well as the cosmological parameters. This is achieved by analysing the diffuse dispersion measure (${\rm DM}_{\rm diff}$) of FRBs as a function of redshift. However, the rapidly increasing data size requests more and more computational resources. In this work, we first develop an accelerated method for any cosmic dispersion measure by deriving an analytical approximation formula for flat, $\Lambda$CDM and $w$CDM universes. Focusing on FRBs, we show that our approximation works well for the ranges $0.01 \leq z \leq 2$, $0.2 \leq \Omega_m \leq 1.0$ and $-3.0 \leq w \leq -0.5$, with relative error to a numerically evaluated ${\rm DM}_{\rm diff}$ always smaller than $3.5 \%$ (in the worst case scenario). This error remains below observationally relevant ${\rm DM}$ scatter and is especially small near the concordance $\Lambda$CDM cosmology. Additionally, we perform a cosmological analysis of simulated FRB data and show that our approximation gives robust and unbiased results, even when applied in regions of parameter space where its relative error becomes larger than $3\%$. Finally, the approximation is more than $15$ ($2$) times faster than the numerical solution of $\Lambda$CDM ($w$CDM), and can reach a timing improvement of a factor of $25$ when used in an MCMC cosmological inference.

The observation of filamentary X-ray structures near bow-shock pulsar wind nebulae (PWNe) -- such as the Guitar, Lighthouse, and PSR J2030$+$4415 nebulae -- and of slow-diffusion regions around pulsars like Geminga, Monogem, and PSR J0622$+$3749, challenges the standard picture of cosmic-ray transport in the interstellar medium, implying a diffusion coefficient two orders of magnitude smaller than the Galactic average. The suppressed diffusion can be attributed to self-generated magnetic turbulence, driven -- via the non-resonant streaming instability -- by electron-positron pairs escaping the PWNe. This instability requires a net current, yet the beam of escaping pairs is expected to be charge-neutral. We show that a charge-neutral pair beam propagating through an electron-proton plasma can spontaneously generate a net current. Using fully kinetic two-dimensional particle-in-cell simulations with realistic mass ratio, we find that beam electrons get focused into self-generated magnetic filaments produced by the nonlinear evolution of the Weibel instability, while beam positrons remain unconfined. We show that in three-dimensional simulations the resulting net (positron) current drives the non-resonant streaming instability, further amplifying the magnetic field. This mechanism provides a pathway for the onset of charge asymmetries in initially charge-neutral pair beams and for the growth of magnetic fluctuations that efficiently scatter the beam particles, with implications for the formation of X-ray filaments and, potentially, for particle self-confinement in TeV halos around PWNe.

Understanding the thermodynamic properties of the circumgalactic medium (CGM) is key to uncovering the baryon cycle in galaxies. Here we present spatial and kinematic relationships between Mg II and O VI as representatives for low and high ion-bearing gas, in the cosmological zoom-in galaxy simulation suite FOGGIE, a set of Milky-way-like galaxy simulations with high CGM resolution. We find the O VI-bearing gas exists as a diffuse halo around the galactic disk, while the Mg II-bearing gas is more centrally located. We investigate the covering fraction, probability of co-observation, co-kinematic correspondence of these ions using two different analysis methods. We make both mock sightlines using two-dimensional projections of our simulations treating these cells as integrated lines of sight and we create one-dimensional ray objects and use the SALSA (Boyd et al. 2020) code to investigate individual gas structures that contribute most to the line of sight column densities, which we call mock absorbers. We explore the relative kinematics of these mock absorbers and find Mg II and O VI appear to have a co-kinematic relationship when looking at absorber pairs with the closest relative velocity like in Werk et al. (2016). However, this does not necessarily correspond with a close spatial separation meaning many O VI and Mg II absorber pairs only appear to be co-kinematic but are physically unrelated. Taking a more holistic look at Mg II and O VI absorber pairs reveals a much weaker correlation between these two ions.

We explore the formation of intermediate mass black holes (IMBHs), potential seeds for supermassive black holes (SMBHs), via runaway stellar collisions for a wide range of star cluster (surface) densities ($4\times10^3 M_\odot$ pc$^{-2} \lesssim \Sigma_\mathrm{h} \lesssim 4\times10^6 M_\odot$ pc$^{-2}$) and metallicities $(0.01 Z_\odot \lesssim Z \lesssim 1.0 Z_\odot)$. Our sample of isolated (>1400) and hierarchical (30) simulations of young, massive star clusters with up to $N=1.8\times10^6$ stars includes collisional stellar dynamics, stellar evolution, and post-Newtonian equations of motion for black holes using the BIFROST code. High stellar wind rates suppress IMBH formation at high metallicities ($Z\gtrsim0.2 Z_\odot$) and low collision rates prevent their formation at low densities ($\Sigma_\mathrm{h}\lesssim 3\times10^4 M_\odot$ pc$^{-2}$). The assumptions about stellar wind loss rates strongly affect the maximum final IMBH masses ($M_\bullet\sim 6000 M_\odot$ vs. $25000 M_\odot$). The total stellar mass loss from collisions and collisionally boosted winds before $t=3$ Myr can together reach up to $5$--$10\%$ of the final cluster mass. We present fitting formulae for IMBH masses as a function of host star cluster $\Sigma_\mathrm{h}$ and $Z$ which can be used to seed SMBHs in high resolution cosmological hydrodynamical simulations and in semi-analytic models for galaxy formation. Our results favour IMBH formation in dense low metallicity environments similar to $z\sim10$ James Webb Space Telescope (\textit{JWST}) proto globular clusters. IMBH formation is suppressed in the high metallicity and low density conditions of the local Universe.

Accretion discs that are tilted with respect to their compact hosts can warp out-of-plane through general relativistic frame-dragging. Warp influences disc dynamics in ways that have been studied extensively, especially as regards instabilities that might lead to rapid angular-momentum cancellation between neighboring rings of fluid and mass infall. We provide a review of warped-disc phenomena here, revisiting key hydrodynamical assumptions that impact calculations of the shear viscosity controlling instability thresholds. Relativistic effects at the level of gas-parcel orbits are included, as are external Lorentz forces applied by the compact primary's magnetic field. Semianalytic analysis reveals that intense magnetic fields can bring about new branches of warp modes and avoided crossings that significantly reduce the perpendicular viscosity at sub-Eddington accretion rates. Critical strengths required for misaligned torques to tear a thin disc may thus relax for systems like neutron star X-ray binaries or radio-loud active galactic nuclei.

The excitation and damping mechanisms for oscillation modes of gas giant planets are undetermined. We show that differential rotation may greatly enhance convective viscosity in giant planets, resulting in damping times of $t_{\rm damp} \sim 10^5-10^6 \, {\rm years}$ for f~modes and low-order p~modes. Radiative diffusion damps p~modes on time scales of $t_{\rm damp} \sim 10^3-10^7 \, {\rm years}$. While the lethargic convective motions cannot effectively excite f~mode or p~modes, storms driven by condensation of water and/or silicates may play a role. High-order p~modes are most effectively excited by cometary/asteroid impacts. Applying these calculations to solar system planets, water storms, rock storms, and impacts may all contribute to exciting the observed f~modes amplitudes of Saturn via ring seismology. Similar f~mode amplitudes with fractional gravitational perturbations of $\delta \Phi/\Phi \sim 10^{-10}-10^{-9}$ are expected for Jupiter and Uranus, apart from their lowest $\ell$ f~modes which could have larger gravitational perturbations of $\delta \Phi/\Phi \sim 10^{-7}$. Rock storms may contribute to mode driving in Jupiter, while water storms are more important for Uranus. The highest-amplitude p~modes are predicted to have periods of $\sim$10-30 minutes, with surface velocities of $\sim$10 {\rm cm/s} for Jupiter and Saturn, and $\sim$1 {\rm cm/s} for Uranus. These oscillation modes may be detectable with radial velocity measurements, ring seismology, or spacecraft Doppler tracking. However, both the damping and excitation physics are uncertain by orders of magnitude, so more careful examination of the relevant physics is required for robust estimates.

We investigate whether Little Red Dots (LRDs) are the dust-reddened, high-inclination counterparts of compact, blue broad-line AGNs (Little Blue Dots, LBDs) powered by super-Eddington accretion. We model the central engine as a geometrically thick, radiation-pressure supported accretion flow whose funnel produces strongly anisotropic, intrinsically blue ionizing continua, coupled to an equatorially concentrated BLR and dusty reprocessing clouds with modest covering factor. Using inclination-dependent spectral energy distributions (SEDs) as input to Cloudy, we show that the extreme broad Halpha EWs of JWST LRDs can be reproduced with global BLR covering factors of only 10%, fully consistent with standard Type~1 AGNs and far below unity. Large Balmer EWs arise because self-shadowing suppresses the high-inclination optical continuum while the BLR is illuminated by an ionizing-rich EUV SED. Weak high-ionization lines (e.g. HeII4686) follow from the orientation-dependent suppression of the XUV/soft X-ray continuum toward equatorial directions, without requiring a fully enclosing gaseous cocoon. Applying a gray dust attenuation law with AV~3 along high-inclination (LRD-selected) sightlines, our fiducial model reproduces the V-shaped UV-optical continua of LRDs and the large Balmer decrements; strong Balmer breaks arise only along the most obscured sightlines. A compact equatorial dust structure with modest global covering factor intercepts and reradiates only a small fraction of the bolometric luminosity, yielding a modest hot-dust bump and far-IR/sub-mm output consistent with current measurements and limits and implying small dust masses. This single-framework model links LRD and LBD observables through orientation, predicting correlated trends in Halpha EW, Balmer decrement, Balmer break, high-ionization line strengths, and IR emission.

Dual quasars separated at the kiloparsec scale are widely regarded as precursors to binary supermassive black holes and offer a key insight into the dynamical evolution of galaxy mergers. Our series of studies focuses on searching for dual quasars by using a selection strategy of zero proper motion and zero parallax to isolate candidate quasars near known ones and by follow-up spectroscopy of the candidates. This paper, the fourth in the series, reports the spectroscopic confirmations of our quasar pair candidates based on the spectroscopic data of the SDSS and DESI DR1. We newly identified 17 dual quasars and 143 projected quasars. The redshifts of the 17 dual quasars range from 0.573 to 2.758, with a median of 1.512. One notable system, J0023+0417, exhibits nearly identical spectral features in the two members and shows evidence of a potential foreground galaxy, making it a high-confidence strong gravitational lensing system. The redshifts of the 143 projected quasars are from 0.301 to 4.030, with a median of 1.596. Among them, four have projected distances below 30 kpc, offering valuable opportunities to probe the circumgalactic medium (CGM) of the foreground host galaxy through absorption lines.

Reliable prediction of the solar cycle is a formidable challenge, yet it is increasingly vital in our technology-dependent society as solar activity drives space weather. Various methods, including precursors, nonlinear curve fitting and extrapolation, statistical and Machine Learning (ML) models, and dynamo and surface flux transport (SFT) models, were implemented to predict past cycles. Analysing about 100 predictions for Solar Cycle 24 and over 130 for Solar Cycle 25, we find that most methods largely failed to predict the peak correctly: Cycle 24 was statistically predicted to be a strong cycle, whereas Cycle 25 was predicted to be a weak cycle. By and large, predictions made only after the cycle began became closer to reality. ML-based models also produced discouraging results. The polar field and its proxy-based predictions are the most physically supported approach to prediction; however, applying them much earlier, before the solar minimum, may yield inaccurate results. Dynamo models are progressively improving both in understanding and in forecasting; however, they need to improve by accurately assimilating the observed polar field data and additional physics, such as meridional flow variations. Solar dynamo theory, complemented by the SFT model and observations, demonstrates that the prediction of a cycle before the time of its previous cycle's maximum is meaningless. The current solar cycle is declining, and the community is now preparing for the prediction of the next cycle. Thus, this review will guide future studies.

HD 142527 is a well-studied intermediate-mass T Tauri star surrounded by a transitional disk with a large dust cavity, spiral structures, and an accreting low-mass companion. Despite extensive observations, the system's inner regions remain poorly understood, particularly regarding their influence on disk morphology and planet formation. This study aims to investigate the inner region of HD 142527 (<50 au) with high detection sensitivity thanks to dedicated postprocessing methods to search for undetected components and explore their potential role affecting the disk's structure and evolution. We analyze high-contrast imaging data obtained with VLT/SPHERE applying PACO and REXPACO algorithms, dedicated respectively to the detection of point-like sources and to the reconstruction of circumstellar disks with high reliability, while relying on both angular and spectral variations. We revisit the known companion HD 142527 B and update its photometry, astrometry and accretion rate estimates. Furthermore, we identify a new candidate companion (CC) at an angular separation of ~0.09" (~14 au), although it may also be a disk feature. Otherwise, it could be a young gas-giant planet or a brown dwarf with a mass of 12-50 $M_\rm{Jup}$. Additionally, we report the discovery of a tightly wound H$\alpha$ spiral feature in the inner disk, reconstructed for the first time by high contrast imaging. The spiral implies varying accretion dynamically linked to the known companion B and possibly to CC, suggesting ongoing interactions that influence the disk's structure. Our findings provide new insights into the complex interactions within the HD 142527 system, highlighting the role of multiple companions in driving disk asymmetries and facilitating planet formation. Future high-resolution observations and dynamical modeling will be essential to fully understand the system's architecture and evolution.

Little red dots (LRDs) are candidate high-redshift supermassive black holes accreting in dense gas. They remain undetected in X-rays. In previous work, we provided the first quantitative models that reproduce the optical and near-infrared spectra of LRDs with the \sirocco radiative transfer code, thereby constraining the properties of the surrounding gas. Here, we use these constraints to predict the X-ray attenuation produced by dense gas cocoons, and explore its dependence on Balmer-break strength, metallicity, intrinsic X-ray spectral energy distribution, and observed bandpass as a function of redshift. The X-ray constraints are very tight, requiring extinction by a Compton-thick gas column ($N_{\rm H}\sim10^{25}\,{\rm cm}^{-2}$) with moderate metallicity $0.05$-$0.1\,Z_\odot$ and intrinsically weak X-ray emission (bolometric to X-ray luminosity ratio, $k_{\rm bol,X}\gtrsim 30$), as observed in narrow-line active galactic nuclei with high accretion rates, to make LRDs sufficiently faint to evade detection. Intrinsically bright X-ray emitters as seen in typical broad-line active galactic nuclei would be detected even behind the typical Compton-thick gas columns with modest metallicity that were inferred from the optical spectra. Very low metallicity objects might be detected in X-rays even with low intrinsic X-ray luminosities, suggesting that LRDs are not (currently) chemically pristine.

The discovery of a population of massive, ancient quiescent galaxies within the first 2 Gyr of the Universe's history has led to significant tensions with models of galaxy formation. However, these analyses are often based on slit spectroscopy, which typically captures only the center-most region of these galaxies and, crucially, assumes these cores are representative of the entire galaxy. To illustrate the varying stellar populations present throughout these galaxies, we present an analysis of color gradients in four $z>3$, $\log(M_\star/M_\odot)>11$ quiescent galaxies which previous works have argued are in tension with models. Using medium-band photometry from MINERVA JWST observations, we measure resolved photometry in a series of elliptical annuli out to $0.7^{\prime\prime}$ ($\sim4~R_e$). We find negative color gradients in three galaxies, and for the most extreme color gradient ($\Delta(U-V)/\Delta R=-0.126\pm0.030~{\rm mag~kpc^{-1}}$), we find the stellar mass is 0.1 dex lower when compared to photometry measured within NIRSpec slits. In the limiting case where these color gradients are entirely driven by age, we find lessened tensions with extreme value statistics models out to $z\sim9.5$, though different stellar population modeling choices also contribute significantly. Ultimately, these findings highlight the need for integral field unit spectroscopy. Spatially-resolved spectra can provide the evidence needed to break the age-dust-metallicity degeneracy, and reliably separate the effects of the observed color gradients from the effects of different physical modeling assumptions on the formation histories of these galaxies.

Solar energetic particle (SEP) events, defined by proton flux exceeding 10 pfu in the > 10 MeV channel, pose major risks to spacecraft operations, astronaut safety, and high-latitude aviation. Due to the complexity and rarity of SEP events, reliable operational SEP forecasting remains an important challenge in space weather. Here we present a novel 24-hour-ahead realtime forecasting framework, SEPNET-PRISM, based on a multi-task learning structure and a thoroughly constructed list of features from multiple sources spanning multiple solar cycles, that jointly predicts SEP event occurrence and future proton and soft X-ray fluxes. SEPNET-PRISM extends the earlier-introduced SEPNET-based models by integrating a broader range of solar observations, including active-region magnetic parameters from SHARP and SMARP, solar-flare information, coronal mass ejections, soft X-ray flux, and historical > 10 MeV proton flux. As compared with SEPNET, the inclusion of SMARP data expands the temporal coverage of magnetic-field predictors to earlier solar cycles, while flux-based inputs provide additional precursor information. Evaluation on the CLEAR SEP benchmark dataset shows improved classification performance over the earlier SEPNET-O (operational version of SEPNET) on the newly aligned dataset. The best operational model is obtained when magnetic, radiative, and proton-flux predictors are combined, highlighting the value of expanded historical coverage and complementary precursor information for improving realtime SEP forecasting.

The 21-cm signal probes the intergalactic medium during the Epoch of Reionization (EoR) but is overwhelmed by astrophysical foregrounds orders of magnitude stronger than the cosmological signal. We evaluate two mitigation strategies: (i) Foreground Avoidance, restricting analysis to the EoR window in Fourier space, and (ii) Foreground Removal via Gaussian Process Regression, which exploits spectral smoothness to statistically separate contaminants and reclaim modes within the contaminated wedge. Both introduce systematic biases of up to $\approx$1$\sigma$ in astrophysical parameters such as the minimum star-forming halo mass and ionising escape fraction, with avoidance posteriors consistently broader than removal owing to its restricted visibility coverage. The global reionization history is recovered within the 95% credible interval, though the neutral fraction at late reionization epochs shows a persistent bias reflecting the difficulty of its inference from the power spectrum alone. Multi-redshift inference is susceptible to contamination from poorly mitigated bins. Excluding such bins significantly reduces parameter biases, but identifying them requires independent quality metrics. When restricted to identical length scales, both strategies recover similar power spectra, yielding posteriors in similar regions.

We present the EGIDE (The Edge-on Galaxies in the DESI survey) project - a catalogue of 149,215 edge-on galaxy candidates created using the data of the DESI Legacy Imaging Survey DR10 images. The catalogue size is ten times bigger than its predecessor and covers more than half of the sky. It is constructed in an automatic way utilizing the full power of manual annotations from the GalaxyZoo volunteers, implemented in the Zoobot neural model, which was fine-tuned to search for edge-on galaxies specifically. To ensure the credibility of the dataset, subsequent manual supervision was done. The EGIDE catalogue provides homogeneous SExtractor photometry in $griz$ bands, total stellar mass estimation, redshift values for 98% of the sample, star formation rates and other information. All of this is publicly available at The Edge-on Galaxy Database site. The preliminary analysis focused on differences between edge-on galaxies in the so-called blue sequence and red cloud populations. These galaxies demonstrate distinct properties: the number of redder galaxies drops with increasing $a/b$ ratio faster than for the bluer galaxies; galaxy thickness varies with galaxy colour: red sequence galaxies are thicker than blue cloud galaxies; the flattening ratio $q=b/a$ increases with total stellar mass $M_{\star}$ for the high-mass end. It is an intriguing result, that the same trend of $q$ increasing is detected from both the statistical models of figures of revolution and direct observations of edge-on galaxies in EGIDE independently. The full extent of the validity of this relationship can only be determined after correctly accounting for the contributions of the bulge and the PSF.

The population of Jupiter-sized exoplanets with orbital periods between 10 and 200 days (WJs) exhibits a broad range of orbital eccentricities and system architectures, suggesting a diversity of formation and migration pathways. In this work, we report the detection and characterization of two new eccentric WJs, TOI-2147 b and TOI-6019 b, initially identified as planet candidates by the Transiting Exoplanet Survey Satellite (TESS). We combined TESS photometry with ground-based follow-up observations, including multiband photometry from LCOGT and MuSCAT2, high-angular-resolution speckle imaging, and high-precision radial velocity measurements from the high-resolution Manfred Hirt Planet Finder Spectrograph (MaHPS). Using these data, we were able to confirm the planetary nature of both candidates. TOI-2147 b has a radius of $10.5 \pm 0.3\,\mathrm{R}_\oplus$ and a mass of $116 \pm 22\,\mathrm{M}_\oplus$. It orbits its slightly metal-poor ($\mathrm{[Fe/H]} = -0.29^{+0.07}_{-0.08}$) G-type host star on an eccentric orbit ($e = 0.29 \pm 0.07$) with a period of 26.2 days. TOI-6019 b has a radius of $12.3 \pm 0.3\,\mathrm{R}_\oplus$ and a mass of $149 \pm 15\,\mathrm{M}_\oplus$. It orbits a slightly evolved, solar-metallicity G-type sub-giant with a period of 14.5 days on a significantly eccentric orbit ($e = 0.48^{+0.05}_{-0.04}$). Both planets have bulk densities below that of Jupiter, indicating mildly inflated radii, with interior structure modeling using GASTLI. This suggests that tidal heating from the nonzero eccentricities likely contributes to this inflation and disfavors large atmospheric metal enrichment. No significant signals from additional companions were detected in the radial velocity time series or transit timing variations. Together with the elevated eccentricities, this is consistent with a high-eccentricity migration origin for both systems.

Discovery of dark-matter (DM) free dwarf galaxies in the NGC 1052 neighborhood has had a considerable impact on modern cosmology. They have been explained through a dwarf--dwarf head-on collision that is a rare event. We find that they could alternatively be associated with a head-on, 1:1 merger after it has been tuned to generate the E4 morphology of NGC 1052. Our simulations show that such mergers produce long-lived tidal features, associated with the remnant galaxy, and in the form of large tidal tails including tidal dwarf galaxies (TDGs). We underline that such tidal features are predicted by the hierarchical scenario in which massive galaxies are formed by galaxy mergers. The latter can reproduce both the tidal features in the NGC1052 outskirts and the observed dwarf galaxies. The simulated TDGs have similar sizes to those observed, while they are ten times smaller in the dwarf bullet scenario. However, we cannot reproduce the luminous globular cluster systems due to resolution limitations. Resolving the radial distance between the DM-free dwarfs is necessary to identify the scenario of their formation. We suggest that there should be many other examples of DM-free dwarf galaxies in the neighborhood of local massive galaxies and galaxy groups.

The present-day architecture of planetary systems contains information about their formation and migration histories. The origin of hot Jupiters (HJs, P $\lesssim$ 10 d, $R_{\rm p} > 8 R_\oplus$) has long been a matter of debate. While most of them are found to be ``lonely'', there is a rare population of HJs hosting small companions on inner orbits (eight known as of May 2026). Their peculiar architecture suggests a gentle disc-migration mechanism. In this study, we present the discovery and characterisation of the multi-planet system TOI-1533, comprising an inner sub-Neptune (TOI-1533 b, $P_{\rm orb} = 3.63$ d, $R_{\rm p} = 3.15 R_\oplus$) and an outer hot giant planet (TOI-1533 c, $P_{\rm orb} = 8.06$ d, $R_{\rm p} > 7.5 R_\oplus$) with substantial H/He by mass ($\rho_{\rm p} < 0.48$ g cm$^{-3}$), both transiting an active K-dwarf star ($T_{\rm eff} \approx$ 5130 K; $V$ (mag) $\approx$ 11). Our joint modelling of stellar activity and planetary signals from radial velocities (HARPS-N) and transits (TESS) allows us to detect their Keplerian signals (approximately $10~\sigma$) and to isolate the stellar modulation. The inclusion of simultaneous photometry in the multi-dimensional Gaussian processes formalism was a fundamental addition to the spectroscopic activity indicators, enabling the disentanglement of stellar activity from planetary signals. The mass ratio of the two confirmed planets ($M_{\rm b} / M_{\rm c}$ about 0.8), together with the super-Neptune mass of the large outer companion ($M_{\rm c} \approx 40 M_\oplus$), makes this system unusual compared to the other few HJs with low-mass inner companions.

Since times immemorial, total solar eclipses have inspired awe and wonder. In the modern scientific era, they have transformed into exclusive natural laboratories, offering fleeting but invaluable opportunities to study the Sun's faint outer atmosphere that is otherwise obscured by the intense glare of the photosphere. This unique vantage point has enabled revolutionary discoveries, from the identification of the element Helium and the first empirical validation of Einstein's General Relativity, to deciphering the corona's surprisingly high temperature. Today, ground-based eclipse experiments provide crucial data that complements and calibrates our space-based solar observatories, and offer high-resolution capabilities in the spatial, temporal as well as spectral domains. This chapter serves as a comprehensive guide detailing how to leverage modern observing equipments, detectors, and advanced computational techniques in image and data processing to conduct meaningful scientific investigations, bridging the gap between historical precedent and cutting-edge research.

Formation of primordial black holes (PBHs) requires a significant enhancement of curvature perturbations. This enhancement generates a twofold gravitational wave (GW) signature: a \emph{low-frequency} stochastic background of scalar-induced GWs (SIGWs) and a distinct \emph{high-frequency} signal from subsequent PBH binary mergers. Assuming a monochromatic PBH mass function, we use PBH abundance constraints on the primordial curvature power spectrum to evaluate the stochastic SIGW background. We also compute the stochastic GW background from mergers of the corresponding PBH binaries, incorporating merger-rate suppression effects to obtain realistic estimates. Furthermore, we derive a model-independent correspondence between the characteristic frequencies of these two signals. This unified framework links these otherwise distinct GW channels, enabling the same primordial fluctuations to be probed across widely separated frequency bands.

The SKA telescopes will bring unparalleled sensitivity across a broad radio band, a wide field of view across the Southern sky, and the capacity for sub-arraying, all of which make them the ideal instruments for studying the pulsar magnetosphere. This chapter describes the advances that have been made in pulsar magnetosphere physics over the last decade, and details how these have been made possible through the advances of modern radio telescopes, particularly SKA precursors and pathfinders. It explains how the SKA telescopes would transform the field of pulsar magnetosphere physics through a combination of large-scale monitoring surveys and in-depth follow-up observations of unique sources and new discoveries. Finally, it describes how the specific observing opportunities available with the AA* and AA4 configurations will achieve the advances necessary to solve the problem of pulsar radio emission physics in the coming years.

In May 2024, NOAA active region (AR) complex 13664/8 was one of the most productive regions of the current solar cycle, producing 12 X-class flares and over 20 coronal mass ejections (CMEs) and triggering the strongest geomagnetic storm since 2003. We investigate why this AR complex was so CME-productive. Using primarily SDO/AIA and SDO/HMI observations, we analyzed its photospheric magnetic evolution, eruption sources, eruption waiting times, and magnetic parameters in comparison with five other ARs. The region initially contained only AR 13664 and exhibited limited flare activity until AR 13668 emerged on May 4, after which clustered major flares and CMEs occurred. Rapid, complex flux emergence substantially increased the region's area, magnetic flux, complexity, and nonpotentiality. The increased complexity was manifested by at least 12 emerging bipoles and six collisional polarity inversion lines (cPILs) formed between nonconjugated polarities, all showing sustained collision and shearing and serving as eruption sources. Decay index distributions show systematically lower critical heights (<45 Mm) above CME source regions. The CME waiting time distribution exhibits two peaks, suggesting that multiple cPILs enhanced CME productivity by increasing both recurrent CMEs from the same source and disturbance-triggered CMEs from nearby sources. These results indicate that, in addition to sufficient nonpotentiality and rapid background-field decay, high magnetic complexity accompanied by dynamical collisional shearing at multiple cPILs was crucial for the extreme CME productivity of AR 13664/8.

We present angular correlation and power spectral analyses of the clustering of VST ATLAS galaxies, LRGs and QSOs and their cross-correlations with the ACT DR6 CMB lensing map with its $\approx3\times$ higher resolution than Planck DR3. We first use $\Lambda$CDM Halofit models at linear scales to estimate galaxy, LRG and QSO bias, $b$, from the ratio of the 3-D power spectra, $b= P_{gg}/P_{gK}$ and hence the mass clustering amplitudes, $\sigma_8$. For galaxies we find $b_G(z=0.15)=1.22\pm0.07$, $\sigma_8(z=0)=0.71\pm0.06$ and for LRGs $b_{LRG}(z=0.26)=2.76\pm0.11$, $\sigma_8(z=0)=0.63\pm0.15$, both consistent with $\Lambda$CDM predictions. But for QSOs, we estimate a bias of $b_Q(z=1.7)=4.44\pm0.24$ and $\sigma_8(z=0)=0.49\pm0.06$, the latter significantly lower than expected for $\Lambda$CDM and leads to an observed $0.15<z<1.7$ gravitational growth rate faster than predicted for $\Lambda$CDM. At smaller non-linear scales, HOD model fits imply a broad halo mass distribution for $z\approx0.15$ galaxies but much narrower mass distributions for LRGs and QSOs. More surprisingly, the HOD fit for QSOs implies $M_{eff}=5\times10^{14}M_\odot$, $\approx10\times$ higher than expected. Finally, while the LRG and QSO HOD models fit $w_{gg}$ and $w_{gK}$ consistently, the same is not true for the $z\approx0.15$ galaxies unless we hypothesise anti-bias with $b_G\approx0.55$. This result supports previous QSO lensing analyses that also found anti-bias at $r<5$h$^{-1}$Mpc. We suggest that while HOD + NFW halo models may well describe LRG mass profiles and clustering, at intermediate scales the dominant galaxy population (and QSOs) may follow the alternative hierarchical clustering model of Peebles(1974) where galaxies trace the mass more closely than in a halo model.

Plasma spectroscopy is a fundamental tool for diagnosing laboratory and astrophysical plasmas. Accurate interpretation of spectra depends upon precise modeling and comprehension of Stark broadening and other mechanisms affecting spectral lines. In this context, computer simulations have emerged as valuable tools, offering idealized experiments with well-defined conditions. Molecular dynamics simulations, in particular, excel at replicating the particle interactions within the plasma and their impact on the state of a radiating atom or ion. However, these simulations present challenges in tracking electron capture processes, since setting an unambiguous criterion to distinguish between bound and free electrons is not trivial. In this paper we introduce a new algorithm that, within a classical framework, precisely identifies the scenario in which an electron is captured by an ion and then follows a stable orbit around it. The algorithm's applicability extends to emitters with charges Z >= 1. The procedure enables the correct identification of valid time-histories of the electric microfield perturbing the emitting ion, which will be used for subsequent line shape calculations. The ionization balance results obtained from the application of this algorithm are compared with an additional method based on the potential energy of the particles in the simulations and with atomic kinetic simulations. For both methods, we find good agreement, therefore validating the use of this approach.

We analyze the Wheeler-DeWitt quantization of a spatially flat Friedmann-Lemaître-Robertson-Walker universe containing pressureless dust and a positive cosmological constant ($\Lambda > 0$). Following relational time framework, we establish a direct mathematical correspondence between the cosmological Hamiltonian and the radial Schrödinger equation for the scattering states of the non-relativistic hydrogen atom. This exact solvability allows us to rigorously construct the physical Hilbert space and ensure the self-adjointness of the Hamiltonian. As a concrete result, we show that the wave packets unitarily evolve depicting a non-singular quantum bounce, systematically replacing the classical Big Bang singularity. Finally, we discuss the physical relevance of this exact solution within the matter-bounce scenario. We demonstrate that this framework provides a robust quantum origin for a bounce during a dust-dominated contracting phase -- a necessary prerequisite for generating a scale-invariant spectrum of primordial perturbations -- derived from the unitary dynamics of the quantized background.

We propose a minimal, fully thermal mechanism that resolves the long-standing tension between achieving the observed dark-matter relic abundance and explaining the astrophysical signatures of self-interactions. The framework introduces two mediators: a light scalar $\phi$ (MeV scale) that yields the required, velocity-dependent self-interactions, and a heavy scalar resonance $\Phi_h$ (TeV scale) with mass $m_{\Phi_h}\!\approx\!2m_\chi$ that opens an $s$-channel resonant annihilation during freeze-out. This clearly decouples early-universe annihilation from late-time halo dynamics. A detailed numerical analysis identified a narrow predictive island of viability. A representative benchmark with $m_\chi\!=\!600$~GeV, $m_\phi\!=\!15$~MeV, and $m_{\Phi_h}\!\simeq\!1.2$~TeV reproduces the relic density and yields $\sigma_T/m_\chi\sim 0.1$--$1~\mathrm{cm}^2\!/\mathrm{g}$ at dwarf-galaxy velocities while satisfying cluster bounds. The model makes sharp, testable predictions: a narrow $t\bar t$ resonance near $1.2$~TeV within HL-LHC reach, and a spin-independent direct-detection signal $\sigma_{\rm SI}\!\sim\!7\times10^{-48}\,\mathrm{cm}^2$ within next-generation sensitivity. As an optional UV completion, we show that walking $\mathrm{SU}(3)_H$ gauge theory with $N_f=10$ naturally realizes the near-threshold relation $m_{\Phi_h}\!\approx\!2m_\chi$ and can furnish an effective anomalous dimension $\gamma\!\approx\!0.5$ which underlies a density-responsive dark-energy sector, suggesting a unified origin for the dark sector.

Axions that couple to electromagnetism are produced in the early Universe by, among other channels, freeze-in via the Primakoff process. For sufficiently large axion masses, the same coupling causes the axions to decay into two photons, which subsequently ionize the intergalactic medium. If this decay occurs in the redshift range $20 \lesssim z \lesssim 1100$, then the contribution to the cosmic microwave background optical depth $\tau_{\rm reio}$ can lead to a conflict with observations, excluding models with sufficiently strongly coupled, heavy axions and high reheating temperatures, $T_{\rm reh}$. Using large ensembles of explicit type IIB string theory models with up to $h^{1,1} = 100$ axions, we compute the full cosmic reionization history caused by the decays of multiple axions. We compare this to the posterior on the high-$z$ component of $\tau_{\rm reio}$ derived from parametric-independent constraints on the ionization state of the Universe, obtained in a full \textit{Planck} analysis presented in a companion paper. For $h^{1,1} = 20, 50, 100$, we find that approximately 15\%, 15\%, and 10\% of the models in the ensemble prefer $T_{\rm reh} \lesssim 10^{10}\,\text{GeV}$ at 95\% CL. We provide a publicly available code at:~\href{this https URL}{this http URL}, which computes the reionization history for arbitrary ensembles of decaying axions. Our analysis opens the door for future large-scale work studying the preference for low-temperature reheating in models with multiple axions.

Core-collapse supernovae serve as powerful probes of light, weakly coupled particles, such as dark photons. The conventional SN1987A cooling bound constrains the dark photon mass-mixing parameter space by requiring that the luminosity from the proto-neutron star core not exceed the observed neutrino emission. In this work, we revisit these limits by including the effect of dark matter (DM) captured inside the progenitor star before collapse. The trapped DM acts as an additional scattering target for dark photons, modifying their free-streaming length and, consequently, the supernova cooling rate. We perform a self-consistent analysis for both annihilating and asymmetric DM scenarios, incorporating light-mediator effects in the capture rate calculation. For annihilating DM, the equilibrium density remains too small to affect the bounds significantly. In contrast, asymmetric DM can accumulate to large densities, leading to the formation of a "dark photosphere" that suppresses the dark-photon luminosity and reopens regions of parameter space previously excluded by supernova cooling. The results presented are intended as a proof-of-principle demonstration of how astrophysical dark matter populations can alter supernova cooling constraints and do not provide precision exclusion limits on any given model.

We revisit the Starobinsky model of inflation in light of recent data from the Atacama Cosmology Telescope (ACT), which indicates a potential preference for a slightly larger scalar spectral index $n_s$ than predicted by the standard $R^2$ scenario. We demonstrate that a natural one-parameter generalization to a quadratic model $\sim L+L^2$ in the Lovelock invariant $L=R+\frac{\alpha}{4}{\cal G}$ ($\cal G$ is the Gauss--Bonnet term), can effectively resolve this minor tension. Scalar-tensor formulation of this theory yields an Einstein-frame Starobinsky-type scalar potential augmented by Gauss--Bonnet and derivative couplings, which modify the inflationary slow-roll dynamics. We show that a non-zero coupling $\alpha$ for the Gauss-Bonnet term can shift $(n_s, r)$ along a trajectory that brings the predictions into better agreement with the ACT likelihood. We also find that $L+L^2$ gravity, in its scalar-tensor formulation, is equivalent to Higgs inflation coupled to the Gauss--Bonnet term, and belongs to the Horndeski/galileon class of modified gravities. This work establishes the quadratic $f(L)$ gravity as a compelling and physically motivated extension that preserves the successes of Starobinsky inflation while improving its fit to modern precision cosmological data.

We study Higgs-like inflation in the framework of scalar-torsion gravity, focusing on the general class of $f(T,\phi)$ theories in which gravitation is mediated by torsion rather than curvature. Motivated by the increasing precision of cosmic microwave background and large-scale-structure observations, we examine whether Higgs-like inflation remains compatible with current data in this extended gravitational setting. Working within the slow-roll approximation, we analyze the inflationary dynamics both analytically and numerically. In the dominant-coupling regime we derive closed-form expressions for the scalar spectral index and the tensor-to-scalar ratio as functions of the number of e-folds, and we subsequently relax this assumption by numerically solving the slow-roll equations. Confrontation with the latest constraints from Planck 2018, ACT DR6, DESI DR1, and BICEP/Keck shows that Higgs-like inflation in $f(T,\phi)$ gravity is fully consistent with current bounds, naturally accommodating the preferred shift in the scalar spectral index and leading to distinctive tensor-sector signatures.

We discuss key aspects of the nature of radiation from global strings and its impact on the relic axion density. Using a simple model we demonstrate the dependence on the spectrum of radiation emitted by strings. We then study the radiation emitted by perturbed straight strings paying particular attention to the difference between the overall phase of the field and the small perturbations about the string solution which are the axions. We find that a significant correction is required to be sure that one is analyzing the axions and not the self-field of the string. Typically this requires one to excise a sizeable region around the string - something which is not usually done in the case of numerical field theory simulations of string networks. We have measured the spectrum of radiation from these strings and find that it is compatible with an exponential, as predicted by the Nambu-like Kalb-Ramond action, and in particular is not a ``hard'' spectrum often found in string network simulations. We conclude by attempting to assess the uncertainties on relic density and find that this leads to a range of possible axion masses when compared to the measured density from the Cosmic Microwave Background, albeit that they are typically higher than what is predicted by the Initial Misalignment Mechanism. If the decay is via a ``soft spectrum'' from loops produced close to the backreaction scale we find that $m_{\rm a}\approx 160\,\mu{\rm eV}$ and a detection frequency $f\approx 38\,{\rm GHz}$. If axions are emitted directly by the string network, and we use emission spectra reported in field theory simulations, then $m_{\rm a}\approx 4\,\mu{\rm eV}$ and $f\approx 1\,{\rm GHz}$, however this increases to $m_a \approx 125\,\mu{\rm eV}$ and $f\approx 30\,{\rm GHz}$ using our spectra for the case of an oscillating string. In all scenarios there are significant remaining uncertainties that we delineate.

General relativity has its successes at the local astrophysical level, however, it seems to be insufficient in describing the Universe at large scales. In this work we investigate how the most general field theories in the context of general relativity can accomodate a phantom-to-quintessence transition which may be essential element of realistic Dark Energy scenarios in the late Universe. As we demonstrate in a very detailed manner, this is impossible for a canonical and minimally coupled single scalar field theory, but it may be possible for ghost condensate theories like $k$-essence theories. We point out how the ghost instabilities may be eliminated, and we analyze the quantitative features of a $k$-essence theory that may realize a phantom-to-quintessence transition in the late Universe. We also qualitatively compare the difficulties and fine-tunings required for $k$-essence theories to realize a phantom-to-quintessence transition, and how such a transition is naturally realized in modified gravity, without unnecessary fine-tunings and ghost eliminations.

By exploiting small-scale structure formation probed by Lyman-$\alpha$ forest observations, we study constraints on a model of dark matter from dark sector decay. We compute the phase space distribution of the dark matter and the linear matter power spectrum. We map the non-thermal dark matter distribution in this dark matter model to an approximate thermal warm dark matter distribution, and use this approximation to obtain a constraint from the Lyman-$\alpha$ forest observation. We combine the latest Lyman-$\alpha$ forest bounds with the constraint from the Big Bang Nucleosynthesis. As these two probes offer highly complementary constraints, we impose strong limits on sub-GeV dark matter. Consequently, masses lighter than $\sim 10^{-1}$ GeV are excluded, thereby significantly limiting the allowed parameter space. More broadly, our findings demonstrate the utility of small-scale structure observations in testing non-thermal dark matter paradigms, offering valuable insights for exploring a wider class of late-time decay models.

Gravitational-wave events are interpreted in terms of Bayesian posteriors for their source properties inferred under unphysical reference priors. Though these parameter estimates are important intermediate data products for downstream analyses, we demonstrate that they are generically biased and therefore should not be used for astrophysical interpretation directly, as is common. Hierarchical parameter estimation is the solution, as joint analysis of the entire catalog of observations reduces statistical uncertainties and actually informs the correct prior, with population-informed event parameters now appropriate for astrophysical interpretation. As an example, we show how the most extreme measurements from a catalog can be derived and used to identify exceptional events from previous and ongoing observing runs, pointing out they are more informative about the population than any individual event. Using LIGO-Virgo-KAGRA data, we thus demonstrate that population inference is not optional to interpret gravitational-wave observations.

We derive a Stern--Gerlach control parameter, $\Xi$, comparing spin-driven cross-sheet displacement with the relativistic Larmor radius. It places heliospheric plasmas and most astrophysical jets in the negligible regime, some stellar-mass black-hole coronae in a transitional regime, and magnetar sheets in a strong, near-QED regime. Relativistic pair-plasma simulations at fixed $\gamma_{\rm tr}=2$ show that increasing $\Xi$ produces, within the coupled SG+$\mathbf J_M$ model, magnetic-moment sorting and magnetization-current feedback that enhance the normalized flux-growth rate through an additional spin-kinetic pathway beyond classical pressure- and geometry-controlled pair-plasma reconnection.

Triple systems have progressively been recognized as ubiquitous in our Universe and provide a good testing ground for wave generation and propagation in nontrivial environments. We study the dynamics of triple systems in a fully nonlinear setting. In particular, we analyze numerical relativity simulations of head-on collisions of black holes in the presence of a companion. We show evidence for Doppler and gravitational redshift in the ringdown, and clear signs of amplification by lensing. In certain cases, we also show the appearance of a second image, with hints of resonant mode excitation. Our results pave the way for the understanding of mergers in the vicinity of massive companions. Even in extreme setups we do not find collapse to black holes from lensed gravitational radiation.