Papers for Thursday, Jul 31 2025

In this study, we define the novel splashback depth $\mathcal{D}$ and width $\mathcal{W}$ to examine how the splashback features of dark matter haloes are affected by the physical properties of haloes themselves. We use the largest simulation run in the hydrodynamic MillenniumTNG project. By stacking haloes in bins of halo mass, redshift, mass-dependent properties such as peak height and concentration, and halo formation history, we measure the shape of the logarithmic slope of the density profile of dark matter haloes. Our results show that the splashback depth has a strong dependence on the halo mass which follows a power law $\mathcal{D}\propto\left(\log_{10}M\right)^{2.8}$. Properties with strong correlation with halo mass demonstrate similar dependence. The splashback width has the strongest dependence on halo peak height and follows a power law $\mathcal{W}\propto\nu^{-0.87}$. We provide the fitting functions of the splashback depth and width in terms of halo mass, redshift, peak height, concentrations and halo formation time. The depth and width are therefore considered to be a long term memory tracker of haloes since they depend more on accumulative physical properties, e.g., halo mass, peak height and halo formation time. They are shaped primarily by the halo's assembly history, which exerts a stronger influence on the inner density profile than short-term dynamical processes. In contrast, the splashback features have little dependence on the short term factors such as halo mass accretion rate and most recent major merger time. The splashback depth and width can therefore be used to complement information gained from quantities like the point of steepest slope or truncation radius to characterise the halo's history and inner structure.

Cosmological hydrodynamical simulations have become an indispensable tool to understand galaxies. However, computational constraints still severely limit their numerical resolution. This not only restricts the sampling of the stellar component and its direct comparison to detailed observations, but also the precision with which it is evolved. To overcome these problems we introduce the \emph{Superstars} method. This method increases the stellar mass resolution in cosmological galaxy simulations in a computationally inexpensive way for a fixed dark matter and gas resolution without altering any global properties of the simulated galaxies. We demonstrate the \emph{Superstars} method for a Milky Way-like galaxy of the Auriga project, improving the stellar mass resolution by factors of $8$ and $64$. We show and quantify that this improves the sampling of the stellar population in the disc and halo without changing the properties of the central galaxy or its satellites, unlike simulations that change the resolution of all components (gas, dark matter, stars). Moreover, the better stellar mass resolution reduces numerical heating of the stellar disc in its outskirts and keeps substructures in the stellar disc and inner halo more coherent. It also makes lower mass and lower surface brightness structures in the stellar halo more visible. The \emph{Superstars} method is straightforward to incorporate in any cosmological galaxy simulation that does not resolve individual stars.

While significant advances have been made in photometric classification ahead of the millions of transient events and hundreds of supernovae (SNe) each night that the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will discover, classifying SNe spectroscopically remains the best way to determine most subtypes of SNe. Traditional spectrum classification tools use template matching techniques (Blondin & Tonry 2007) and require significant human supervision. Two deep learning spectral classifiers, DASH (Muthukrishna et al. 2019) and SNIascore (Fremling et al. 2021) define the state of the art, but SNIascore is a binary classifier devoted to maximizing the purity of the SN Ia-norm sample, while DASH is no longer maintained and the original work suffers from contamination of multi-epoch spectra in the training and test sets. We have explored several neural network architectures in order to create a new automated method for classifying SN subtypes, settling on an attention-based model we call ABC-SN. We benchmark our results against an updated version of DASH, thus providing the community with an up-to-date general purpose SN classifier. Our dataset includes ten different SN subtypes including subtypes of SN Ia, core collapse and interacting SNe. We find that ABC-SN outperforms DASH, and we discuss the possibility that modern SN spectra datasets contain label noise which limit the performance of all classifiers.

The dynamics of star forming gas can be affected by many physical processes, such as turbulence, gravity, supernova explosions, and magnetic fields. In this paper, we investigate several nearby star forming regions (Orion, Upper Sco, Taurus, and Perseus) for kinematic imprints of these influences on the newly formed stars. Using Gaia DR3 astrometry and APOGEE DR17 radial velocities, we compute first-order velocity structure functions (VSFs) of young stars in galactic Cartesian coordinates in both 6D (3D positions and 3D velocities) and 4D (3D positions and each 1D velocity) to identify signatures of turbulence and anisotropic motion. We also construct 3D and 1D radial velocity profiles to identify coherent expansion trends, and compare stellar proper motions to plane-of-sky magnetic field orientations in Taurus and Perseus. We find that the VSFs are mildly anisotropic, with slightly different amplitudes, slopes, or features in different directions in several groups, but in general, they are all consistent with Larson's Relation at intermediate length scales, especially in less compact groups. In several cases, the VSFs exhibit features suggestive of local energy injection from supernovae. Radial velocity profiles reveal clear anisotropic expansion in multiple groups, with the most extreme cases corresponding to those with the most anisotropic VSFs. In Perseus, we find that the motions of young stars are preferentially perpendicular to the local magnetic field. We find multiple, overlapping causes in each group for the observed kinematics. Our findings support that young stars remember more than just the turbulent state of their natal clouds.

Cosmic birefringence, arising from a potential parity-violating interaction between cosmic microwave background (CMB) photons and evolving pseudo-scalar fields such as axion-like particles, can rotate the CMB polarization plane and induce an effective correlation between the CMB $E$- and $B$-mode polarization. In this work, we introduce a hybrid internal linear combination (ILC) method that combines both $E$- and $B$-mode frequency maps into the component separation pipeline, enabling the disentanglement of correlated and uncorrelated components of CMB polarization in the presence of cosmic birefringence and instrumental polarization angle miscalibration. We derive an analytic linear relation connecting the birefringence-induced correlated component of the CMB $E$- (or $B$-) mode field to the full CMB $B$- (or $E$-) mode field convolved with a modulating field. By performing linear regression between these fields across multiple sky patches, we directly estimate the birefringence angle at the field level. This allows us to distinguish cosmic birefringence from polarization angle miscalibration and foreground contamination, as the ILC responds differently to achromatic cosmic birefringence and chromatic systematic effects, with its weights projecting spatial or harmonic dependence only onto the latter. This non-parametric, field-level approach provides a novel way to probe cosmic birefringence directly in real space. When applied to realistic simulations of the forthcoming LiteBIRD satellite mission, our method yields constraints that are competitive with, and complementary to, existing power spectrum-based analyses. When applied to Planck Release 4 (PR4) data, we find a birefringence angle of $\beta = 0.32^\circ \pm 0.12^\circ$, a $2.7\sigma$ detection that remains robust against varying sky fractions.

We present an efficient and accurate pipeline for the analysis of the redshift-space galaxy bispectrum multipoles at one-loop order in effective field theory (EFT). We provide a systematic theory derivation based on power counting, which features the first comprehensive treatment of stochastic EFT contributions -- these are found to significantly improve the match to data. Our computational pipeline utilizes the COBRA technique that expands the linear matter power spectrum over a basis of principal components based on a singular value decomposition, allowing the cosmology dependence to be captured to sub-permille accuracy with just eight templates. This transforms the problem of computing the one-loop EFT bispectrum to a simple tensor multiplication, reducing the computation time to around a second per cosmology with negligible loss of accuracy. Using these tools, we study the cosmological information in the bispectrum by analyzing PTChallenge simulations, whose gigantic volume provides the most powerful test of the one-loop EFT bispectrum so far. We find that the one-loop prediction provides an excellent match to the bispectrum data up to $k_{\rm max}=0.15~h$ Mpc$^{-1}$, as evidenced by the precise recovery of the dark matter density $\omega_\text{cdm}$, Hubble constant $H_0$, and mass fluctuation amplitude $\sigma_8$ parameters, and the amplitude of equilateral primordial non-Gaussianity (PNG) $f_{\rm NL}^{\rm equil}$. Combined with the power spectrum, the COBRA-based one-loop bispectrum multipoles yield tighter constraints than the tree-level bispectrum monopole, with the posteriors on $\omega_{\text{cdm}}$, $H_0$, and $\sigma_8$ shrinking by 43\%, 31\%, and 4\%, respectively. This suggests that the COBRA-based bispectrum analysis will be an important tool in the interpretation of data from ongoing redshift surveys such as DESI and Euclid.

What is Dark Matter, and what is its concentration in the Milky Way remains an open question in physics. We show that if a significant fraction of dark matter is composed of sub-solar mass primordial black holes (PBHs), gravitational bremsstrahlung resulting from hyperbolic encounters between unbound PBHs within the galactic halos can generate distinctive chromatic gravitational-wave (GW) emission concentrated around the galactic dark matter halo, and it provides a direct window to discover such compact objects. We find that for both generalized NFW and Einasto dark matter profiles of Milky Way, the signal-to-noise ratio can be more than five in one year of observation for the upcoming ground based GW observatories Cosmic Explorer if PBH dark matter fraction $f_{\rm PBH} = 1$ over the mass range $10^{-14} M_\odot \lesssim m_{\rm PBH} \lesssim 10^{-8} M_\odot$. Our results show that the Galactic Center could appear as a GW-bright source, enabling new insights into dark matter and its distribution.

The advent of large aperture arrays, such as the ones currently under construction for the SKA project, allows for observing the Universe in the radio-spectrum at unprecedented resolution and sensitivity. To process the enormous amounts of data produced by these telescopes, scalable software pipelines are required. This paper helps address this by proposing a framework that allows for decentralized radio-interferometric image reconstruction, parallelizing by spatial frequency. This is achieved by creating pseudo-full-resolution problems for each node by using the local visibilities together with previous major cycle reconstructed images from the other nodes. We apply the proposed framework to both multiscale CLEAN and sparsity regularized convex reconstruction and compare them to their serial counterparts across four different data sets of varying properties in the context of two visibility partitions. We found that the parallelization framework allows for significantly improved reconstruction times for images of similar quality. This was especially the case for our larger datasets where we were able to achieve close to the optimal $2\times$ speedup.

We investigate the diversity of dark matter (DM) density profiles in a large sample of late-type galaxies from the SPARC database, with the goal of testing whether a cusp-to-core transition occurs across galaxy mass scales. We perform Bayesian fits to high-quality rotation curves using flexible halo models that allow for variations in the inner slopes of DM density profiles. We quantify the central dark matter structure using the surface density within the inner region of the halo, defined as $\Sigma_{\rm DM}(<0.01r_{V_{\rm max}})$, and compare the SPARC galaxies with Milky Way dwarf satellites as well as galaxy groups and clusters. Our results reveal significant diversity in the inner density slopes of SPARC galaxies, ranging from steep cusps to shallow cores, and show that many of them lie below the cuspy profiles predicted by the cold dark matter model, consistent with core-like structures. In contrast, both lower-mass dwarf galaxies and higher-mass galaxy clusters tend to follow the cuspy DM halos. These findings suggest that baryonic feedback may induce a cusp-to-core transition in Milky Way-mass galaxies, as predicted by hydrodynamical simulations. However, observational limitations and modeling uncertainties still prevent a definitive conclusion. This study provides new empirical insights into the halo mass-dependent nature of DM inner structures and the role of baryonic processes in shaping them.

We present SLIDE, a pipeline that enables transient discovery in data from the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), using archival images from the Dark Energy Camera (DECam) as templates for difference imaging. We apply this pipeline to the recently released Data Preview 1 (DP1; the first public release of Rubin commissioning data) and search for transients in the resulting difference images. The image subtraction, photometry extraction, and transient detection are all performed on the Rubin Science Platform. We demonstrate that SLIDE effectively extracts clean photometry by circumventing poor or missing LSST templates. This is especially useful for transient analysis in the early years of LSST, when template coverage will be largely incomplete or when templates may be contaminated by transients present at the time of acquisition. We present multiband light curves for a sample of known transients, along with new transient candidates identified through our search. Finally, we discuss the prospects of applying this pipeline during the main LSST survey.

We consider the recently discovered \emph{pair-Alfvén shock wave} occurring in collisionless electron-positron plasmas. We perform a series of Particle-In-Cell studies in one and two dimensions in order to determine the stability conditions for such a shock and the mechanisms which sustain its growth. Building on our previous simulations, which established that these shocks are initially mediated by the Weibel instability before becoming Alfvénic, we demonstrate that the shock is sustained by self-generated Alfvén waves overtaking the shock in the upstream plasma. As a result growth is only possible when the guiding magnetic field strength, and hence Alfvén speed, is sufficiently small, $\omega_c \lesssim 0.4$ (in normalized units). Furthermore the production of the waves in the upstream is dependent on a resonance between the Alfvén wave mode and the thermal noise in the plasma, which is inhibited at high magnetization. This explains the conditional absence of this type of shock identified previously.

The High-Altitude Water Cherenkov (HAWC) observatory was designed to study gamma-ray sources in the energy range between a few hundred GeV up to few hundred TeV. It is composed of 300 Water Cherenkov Detectors (WCDs) that cover a surface of approximately 22000 m${}^2$, at 4100 m a.s.l. In this study, we use the HAWC WCDs as a very large horizontal particle tracker, searching for horizontal muon rate variations by season using 1.5 years of HAWC data. We look for a possible correlation between the effective temperature and the horizontal muon rate. In order to do this, we developed a method to calculate the effective temperature for the horizontal propagation of muons. This is the first time that a search for seasonal variations in the high-altitude horizontal muon rate is reported.

Dust-obscured galaxies (DOGs) with extremely red optical-to-infrared colors are often associated with intense starburst and AGN activity. Studying DOGs can provide insights into the processes that drive the growth of galaxies and their central supermassive black holes. However, the general DOG population is heterogeneous, spanning a wide range of evolutionary stages, and has X-ray obscuring column densities ($N_\mathrm{H}$) covering low-to-high levels. In this work, we focus on seven high Eddington ratio DOGs ($\log \lambda_\mathrm{Edd} \gtrsim -0.5$) to examine their X-ray obscuration properties using new and archival X-ray observations. We confirm that these systems are generally heavily obscured, with 6/7 having $N_\mathrm{H}\gtrsim10^{23}~\mathrm{cm^{-2}}$ and 3/7 having $N_\mathrm{H}\gtrsim10^{24}~\mathrm{cm^{-2}}$. Based on the observed similarity with the rare Hot DOG population, we argue that both high-$\lambda_\mathrm{Edd}$ DOGs and Hot DOGs likely trace the post-merger phase during which AGNs are enshrouded by large columns of dust-rich material.

We report improved measurements of the all flavor astrophysical neutrino spectrum with IceCube by combining complementary neutrino samples in two independent analyses. Both analyses show evidence of a harder spectrum at energies below $\sim$30~TeV compared to higher energies where the spectrum is well characterized by a power law. The spectrum is better described by a log parabola or a broken power law, the latter being the preferred model. Both, however, reject a single power law over an energy range 5~TeV-10~PeV with a significance $>4\sigma$, providing new constraints on properties of cosmic neutrino sources.

The IceCube South Pole Neutrino Observatory has discovered the presence of a diffuse astrophysical neutrino flux at energies of TeV and beyond using neutrino induced muon tracks and cascade events from neutrino interactions. We present two analyses sensitive to neutrino events in the energy range \SI{1}{TeV} to \SI{10}{PeV}, using more than 10 years of IceCube data. Both analyses consistently reject a neutrino spectrum following a single power-law with significance $>4\,\sigma$ in favor of a broken power law. We describe the methods implemented in the two analyses, the spectral constraints obtained, and the validation of the robustness of the results. Additionally, we report the detection of a muon neutrino in the MESE sample with an energy of $11.4^{+2.46}_{-2.53} $\,\si{PeV}, the highest energy neutrino observed by IceCube to date. The results presented here show insights into the spectral shape of astrophysical neutrinos, which has important implications for inferring their production processes in a multi-messenger picture.

As tracers of the underlying mass distributions, galaxies' peculiar velocities are valuable probes of the Universe, allowing us to measure the Hubble constant or to map the large scale structure and its dynamics. Yet, catalogs of peculiar velocities are noisy, scarce and prone to various interpretation biases. We aim to measure the bulk flow and the Hubble constant directly from the largest available sample of peculiar velocities, and without imposing a cosmological prior on the velocity field. To address these issues, we analyze the Cosmicflows-4 catalog, the most extensive catalog of galaxy peculiar velocities, reaching a redshift $z=0.1$. Specifically, we construct a forward modeling approach assuming only a flat Universe, which reconstructs the radial and bulk flows of the velocity field directly from measurements of peculiar velocities. Our method accurately recovers cosmological parameters within a radius of $120\ \mathrm{Mpc/h}$ that can then be compared with the predictions from $\Lambda{\rm CDM}$. Apart from a $3\sigma$ tension with $\Lambda{\rm CDM}$ on the magnitude around $120\ \mathrm{Mpc/h}$ associated with a $4\sigma$ tension on the supergalactic $X$ direction, we find a general agreement between the standard model and the observations. Lastly, our analysis suggests a Hubble constant value of approximately $75.8\pm0.4\ \mathrm{km/s/Mpc}$, exacerbating (or independently confirming) the existing ``Hubble tension", however, for the first time accomplished with the largest set of peculiar velocities in existence.

Intrinsically faint galaxies at $z\sim2-3$ offer critical insights into early galaxy formation, tracing low-metallicity, low-mass systems during Cosmic Noon and serving as analogs to reionization-era galaxies. We present ultra-deep JWST/NIRSpec spectroscopy of nine low-luminosity galaxies ($-17 \lesssim M_{\rm UV} \lesssim -20$, $M_\star \lesssim 10^9\,M_\odot$) at $z\sim2.5$ from the CECILIA program, with $\sim$29.5 hr in G235M/F170LP and 1 hr in G395M/F290LP. Our sample includes four LAEs, three rest-UV color-selected galaxies, and two serendipitous detections -- providing the most sensitive rest-optical spectra of individual faint galaxies at this epoch to date. Balmer-line measurements reveal low SFRs ($0.63 < \mathrm{SFR}/(M_\odot\,\mathrm{yr}^{-1}) < 5.43$) and a broad range of dust reddening ($0 < E(B-V) < 1$), with SFRs systematically below those of continuum-selected galaxies. Electron densities are low ($n_e \lesssim 200$cm$^{-3}$), and emission-line diagnostics indicate low [NII]/H$\alpha$, high [OIII]/H$\beta$, suggesting metallicities $12+\log({\rm O/H})\lesssim8.0$. We also present the first O1-BPT constraints in such faint high-redshift galaxies. Notably, two galaxies show low [OIII]/H$\beta$ despite high Ly$\alpha$ EWs and very low [NII]/H$\alpha$, consistent with the predicted turnover in this ratio at very low metallicities -- highlighting the need for complementary diagnostics (e.g., N2, O32) to identify metal-poor systems. Direct $T_e$-based abundances and expanded samples are needed to further trace metallicity and ionization trends in low-mass galaxies.

Understanding reflectance-related quantities for worlds enables effective comparative planetology and strengthens mission planning and execution. Measurements of these properties for Earth, especially its geometric albedo and phase function, have been difficult to achieve due to our Terrestrial situation -- it is challenging to obtain planetary-scale brightness measurements for the world we stand on. Using a curated dataset of visual phase-dependent, disk-averaged observations of Earth taken from the ground and spacecraft, alongside a physical-statistical model, this work arrives at a definitive value for the visual geometric albedo of our planet: 0.242 (+0.005/-0.004). This albedo constraint is up 30--40% smaller than earlier, widely-quoted values. The physical-statistical model enables retrieval-like inferences to be performed on phase curves, and includes contributions from optically thick clouds, optically thin aerosols, Rayleigh scattering, ocean glint, gas absorption, and Lambertian surface reflectance. Detailed application of this inverse model to Earth's phase curve quantifies contributions of these different processes to the phase-dependent brightness of the Pale Blue Dot. Model selection identifies a scenario where aerosol forward scattering results in a false negative for surface habitability detection. Observations of phase curves for Earth at redder-optical or near-infrared wavelengths could disentangle ocean glint effects from aerosol forward scattering and would help with understanding the utility of phase curve observations for the under-development Habitable Worlds Observatory.

Context. Phase curves of small bodies are useful tools to obtain their absolute magnitudes and phase coefficients. The former relates to the object's apparent brightness, while the latter relates to how the light interacts with the surface. Data from multi-wavelength photometric surveys, which usually serendipitously observe small bodies, are becoming the cornerstone of large statistical studies of the Solar System. Nevertheless, to our knowledge, all studies have been carried out in visible wavelengths. Aims. We aim to provide the first catalog of absolute magnitudes in near-infrared filters (Y, J, H, and K). We will study the applicability of a non-linear model to these data and compare it with a simple linear model. Methods. We compute the absolute magnitudes using two photometric models: the HG* 12 and the linear model. We employ a combination of Bayesian inference and Monte Carlo sampling to calculate the probability distributions of the absolute magnitudes and their corresponding phase coefficients. We use the combination of four near-infrared photometric catalogs to create our input database. Results. We produced the first catalog of near-infrared magnitudes. We obtained absolute magnitudes for over 10 000 objects (with at least one absolute magnitude measured), with about 180 objects having four absolute magnitudes. We confirmed that a linear model that fits the phase curves produces accurate results. Since a linear behavior well describes the curves, fitting to a restricted phase angle range (in particular, larger than 9.5 deg) does not substantially affect the results. Finally, we also detect a phase-coloring effect in the near-infrared, as observed in visible wavelengths for asteroids and trans-Neptunian objects.

We present Chandra observations of 63 sources from the Gemini Near Infrared Spectrograph$-$Distant Quasar Survey (GNIRS-DQS) of which 54 were targeted by snapshot observations in Cycle 24. A total of 55 sources are clearly detected in at least one X-ray band, and we set stringent upper limits on the X-ray fluxes of the remaining eight sources. In combination with rest-frame ultraviolet-optical spectroscopic data for these sources, we assess whether X-rays can provide a robust accretion-rate indicator for quasars, particularly at the highest accessible redshifts. We utilize a recently modified H$\beta$-based Eddington luminosity ratio estimator, as well as the C IV $\lambda$1549 emission-line parameter space to investigate trends and correlations with the optical-X-ray spectral slope ($\alpha_{\mathrm{ox}}$) and the effective hard-X-ray power-law photon index ($\Gamma$). We find that $\alpha_{\mathrm{ox}}$ does not improve current accretion-rate estimates based on H$\beta$ or C IV. Instead, within the limitations of our sample, we confirm previous findings that the C IV parameter space may be a better indicator of the accretion rate up to $z\sim3.5$. We also find that the average $\Gamma$ values for a small subset of our sources, as well as the average $\Gamma$ value in different groupings of our sources, are consistent with their respective relatively high Eddington luminosity ratios. Deeper X-ray observations of our X-ray-detected sources are needed for measuring $\Gamma$ accurately and testing whether this parameter can serve as a robust, un-biased accretion-rate diagnostic.

We investigate cool HI gas traced by Lyman series absorption in the vicinity of 256 galaxies (median stellar mass, log10(M*/Msun) = 8.7 and median redshift, z = 0.48) using 15 background quasars (median impact parameter, D ~ 140 physical kpc (pkpc)), as part of the MUSE Quasar-fields Blind Emitters Survey (MUSEQuBES). We find that the HI column density (N(HI)) profile around isolated, star-forming galaxies, which span ~ 3 dex in stellar mass, is well described by a power law with a slope of ~ -2.3 when expressed as a function of the normalized impact parameter, D/Rvir. The HI covering fraction (k) within the virial radius, for a threshold of log10(N(HI)/cm^{-2}) = 14, is significantly lower for high-mass passive galaxies compared to their isolated, star-forming counterparts. The k-profile of isolated, star-forming galaxies suggests a characteristic size of the HI-rich CGM of ~ 1.5 Rvir, consistently across the stellar mass range. The mean HI mass in the outer CGM (0.3-1 Rvir ) increases with stellar mass, ranging from ~ 10^{5} to 10^{6.6} Msun. The b-parameters of the strongest HI components correlate and anti-correlate with specific star-formation rate (sSFR) and mass, respectively, with > 2 sigma significance. Broad Lya absorbers (BLAs) with b > 60 km s^{-1} are predominantly associated with high-mass galaxies, likely tracing the warm-hot phase of the CGM. The velocity centroids of HI components indicate that the absorbers at D < Rvir are predominantly consistent with being gravitationally bound to their associated galaxy, regardless of stellar mass. Finally, leveraging ~3000 galaxies with spectroscopic redshifts from the wide-field Magellan follow-up of six MUSEQuBES fields, we find that non-isolated galaxies exhibit an HI-rich environment that extends approximately three times further than for their isolated counterparts.

We present primordial non-Gaussianity predictions from a new high-precision code for simulating axion-U(1) inflation on a discrete lattice. We measure the primordial scalar curvature power spectrum and bispectrum from our simulations, determining their dependence on both scale and axion-gauge coupling strength. Both the gauge-sourced power spectrum and the bispectrum exhibit a strong blue tilt due to our choice of an $\alpha$-attractor inflaton potential. We provide fitting functions for the power spectrum and bispectrum that accurately reproduce these statistics across a wide range of scales and coupling strengths. While our fitting function for the bispectrum has a separable form, results from high-resolution simulations demonstrate that the full shape is not separable. Thus, our simulations generate realizations of primordial curvature perturbations with nontrivial correlators that cannot be generated using standard techniques for primordial non-Gaussianity. We derive bounds on the axion-gauge coupling strength based on the bispectrum constraints from the cosmic microwave background, demonstrating a new method for constraining inflationary primordial non-Gaussianity by simulating the nonlinear dynamics.

In this work, we explore the generalized holographic dark energy (HDE) scenario. We relate the HDE density to the future-horizon scale via a function $f(a)$, which we reconstruct via spline-based nodal interpolation. We perform a Bayesian analysis to assess model's consistency with current observations, including baryon acoustic oscillations (BAO) from the Dark Energy Spectroscopic Instrument (DESI) DR1, Type Ia supernovae (SNe Ia) from the Union3 and Pantheon+ compilations, and local measurements of the Hubble constant, $H_0$, from SH0ES. We show that under specific conditions, the model reduces to $\Lambda$CDM with one node. We find strong statistical evidence against the standard HDE model $\Delta \chi^2_{\rm HDE, \; 3\text{-Node}} \sim 20$, and in contrast, the reconstructed HDE model, with three nodes, provides a better fit to the data than the $\Lambda$CDM model, specifically, $\Delta \chi^2_{\Lambda \rm CDM, \; 3\text{-Node}} \sim 12$.

Evidence suggests the existence of a large planet in the outer Solar System, Planet Nine, with a predicted mass of 6.6 +2.6 / -1.7 Earth masses (Brown et al., 2024). Based on mass radius composition models, planet formation theory, and confirmed exoplanets with low mass and radius uncertainty and equilibrium temperature less than 600 K, we determine the most likely composition for Planet Nine is a mini-Neptune with a radius in the range 2.0 to 2.6 Earth radii and a H-He envelope fraction in the range of 0.6 percent to 3.5 percent by mass. Using albedo estimates for a mini-Neptune extrapolated from V-band data for the Solar Systems giant planets gives albedo values for Planet Nine in the range of 0.47 to 0.33. Using the most likely orbit and aphelion estimates from the Planet Nine Reference Population 3.0, we estimate Planet Nines absolute magnitude in the range of -6.1 to -5.2 and apparent magnitude in the range of +21.9 to +22.7. Finally, we estimate that, if the hypothetical Planet Nine exists and is detected by upcoming surveys, it will have a resolvable disk using some higher resolution world class telescopes.

Cataclysmic variable (CVs) are close interacting binaries in which a white dwarf accretes materials from a low mass main sequence companion. CVs can experience nova eruptions due to low mass transfer rates. In the standard theory of CV evolution, the ejected materials during nova eruptions are assumed to leave the system in the form of fast, isotropic, optically thick winds, which predicts that novae only result in positive variation (expansion) of orbital period (i.e. positive $\Delta P$). In addition, the angular momentum losses (magnetic braking and gravitational radiation) only predicts a steady long-term decay in the orbital period of CVs, i.e. $\dot P$ is negative. Interestingly, an observation lasting over 30 years reveals positive and negative values for both $\Delta P$ and $\dot P$ in CVs, strongly conflicting with the standard evolutionary patterns. However, it cannot be excluded that these observations originate from short-term phenomena caused by nova eruptions because of a short timescale of observations. In this paper, we model the effect of instantaneous nova eruptions on the evolution of CVs, considering three mechanisms associated with mass loss in nova eruptions, including fast wind, Frank jet and binary-driven mass loss. By assuming that the observed $\Delta P$ and $\dot P$ are dominated by short-term phenomena, our results show that the binary-driven mass loss can explain almost all of the observations of normal CVs. However, the Frank jet may be needed for some of long-period CVs with evolved companions.

We present a novel method to differentiate stream-like and shell-like tidal remnants of stellar systems in galactic halos using the density-based approach of the clustering algorithm AstroLink. While previous studies lean on observation, phase-space, and action-space based criteria for stream and shell determination, we introduce AstroLink's ordered-density plot and cluster identification as a viable tool for classification. For a given data set, the AstroLink ordered-density plot reveals the density-based hierarchical clustering structure from which the resultant clusters are identified as being statistically significant overdensities. Using simulations of sub-halo disruptions in an external potential to generate samples of tidal structures, we find that the curvature of the ordered-density plot is positive for stream-like structures and negative for shell-like structures. Comparisons with more standard classification techniques reveal strong agreement on which structures typically fit into stream-like and shell-like categories. Furthermore, we investigate the properties of clustered stream and shell samples in radial phase space and energy-angle space. Given the sensitivity of stellar tidal structures to their host dark matter halos, the identification and subsequent classification of these structures provide exciting avenues of investigation in galactic evolution dynamics and dark matter structure formation.

As space traffic continues to increase in the cislunar region, accurately determining the trajectories of objects operating within this domain becomes critical. However, due to the combined gravitational influences of the Earth and Moon, orbital dynamics in this region are highly nonlinear and often exhibit chaotic behavior, posing significant challenges for trajectory determination. Many existing methods attempt to address this complexity using machine learning models, advanced optimization techniques, or sensors that directly measure distance to the target, approaches that often increase computational burden and system complexity. This work presents a novel initial orbit determination (IOD) algorithm for cislunar objects based solely on three angle-only measurements taken at three discrete times. The core methodology builds upon a differential corrections framework which iteratively refines the target's trajectory to satisfy line-of-sight constraints. Numerical simulations demonstrate that the algorithm enables accurate IOD using minimal observational data, making it particularly suited for onboard implementation in resource-constrained cislunar missions.

Gravitational-wave observations provide a unique window into the fundamental nature of massive objects. In particular, neutron star equations of state have been constrained due to the success of gravitational wave observatories. Recently, the possibility of detecting dark matter-admixed neutron stars via ground-based laser interferometry has been explored. Dark matter would impact the gravitational waveform of an inspiraling neutron star system through tidal parameters, namely the tidal deformability ($\lambda$, incurring a phase shift to the frequency evolution of the signal. This phase shift would depend both on the percentage of dark matter within the star and its particle nature, e.g., bosonic or fermionic. Indirect detection of dark matter through admixture within neutron stars can provide insight into the neutron equation of state, as well as constraints on the density of dark matter in the universe. In this work, we introduce \texttt{Darksuite}, a proposed extension of the \lal{} software framework, designed to model the gravitational wave signatures of dark-matter-admixed neutron stars. This framework employs simulations from the two-fluid, generally relativistic Tolman-Oppenheimer-Volkoff equations, wherein one fluid is ordinary nuclear matter and the other is dark matter. We demonstrate interpolation of values from a bank of simulations, enabling the study of binary systems where at least one component may be a dark-matter-admixed neutron star. By leveraging existing methodologies within \lal{} for tidal phase corrections and supplementing them with dark matter effects, \texttt{Darksuite} provides a means to generate and analyze gravitational waveforms for these exotic systems.

The exponential-tail behaviours of the probability density function (PDF) of the primordial curvature perturbation are confirmed in the mild-waterfall variants of hybrid inflation with the use of the stochastic formalism of inflation. On top of these tails, effective upper bounds on the curvature perturbation are also observed, corresponding to the exact hilltop trajectory during the waterfall phase. We find that in the model where the leading and higher-order terms in the expansion of the inflaton potential around the critical point are fine-tuned to balance, this upper bound can be significantly reduced, even smaller than the primordial black hole (PBH) threshold, as a novel perturbation-reduction mechanism than the one proposed in Ref. [1]. It makes PBH formation much difficult compared to the Gaussian or exponential-tail approximation. We also introduce Johnson's $S_U$-distribution as a useful fitting function for the PDF. It reveals a nonlinear mapping between the Gaussian field and the curvature perturbation, which enables us to apply the peak theory to estimate the PBH function.

The presence of Active Galaxy Nuclei (AGN) can affect the morphological classification of galaxies. This work aims to determine how the contribution of AGN affects the most used morphological parameters down to the redshift of z~2 in COSMOS-like conditions. We use a sample of > 2000 local non-active galaxies, with a well-known visual morphological classification, and add an AGN as an unresolved component that contributes to the total galaxy flux with 5%-75%. We moved all the galaxies to lower magnitudes (higher redshifts) to map the conditions in the COSMOS field, and we measured six morphological parameters. The greatest impact on morphology occurs when considering the combined effect of magnitude, redshift and AGN, with spiral galaxies being the most affected. In general, all the concentration parameters change significantly if the AGN contribution is > 25% and the magnitude > 23. We find that the GINI coefficient is the most stable in terms of AGN and magnitude/redshift, followed by the M20, Conselice-Bershady (CCON), and finally the Abraham (CABR) concentration indexes. We find that, when using morphological parameters, the combination of CABR, CCON and asymmetry is the most effective in classifying active galaxies at high-redshift, followed by the combination of CABR and GINI.

The interaction cross section of charged pions with air nuclei is a critical parameter for accurately simulating extensive air showers. Improving the modeling of high-energy pion interactions is essential for addressing the muon puzzle-the observed deficit of muons in simulations compared to indirect experimental estimates. As collider experiments cannot directly probe these interactions, we propose a novel measurement approach using muon bundles detected in deep-underground water Cherenkov detectors, such as IceCube and KM3NeT. This method aims to constrain the pion-air inelastic cross section, thereby reducing uncertainties in air shower simulations and advancing our understanding of cosmic ray interactions.

When dense cores in molecular clouds or filamentary structures collapse and form protostars, they may undergo fragmentation and form binary or multiple systems. In this paper, we investigated the key mechanisms influencing fragmentation by comparing the physical conditions of fragmented and unfragmented dense cores (~0.1 pc) in Orion A. Utilizing archival submillimeter continuum data from the James Clerk Maxwell Telescope (JCMT) and the Atacama Large Millimeter/Submillimeter Array survey of Class 0 and I protostars at a 0.''1 resolution, we identified 38 dense cores hosting single protostars and 15 cores hosting binary or multiple systems. We measured the dense cores properties with the Herschel dust temperature, Nobeyama 45m N$_2$H$^+$ J=1-0, and JCMT polarization data. Our results reveal that the dense cores hosting binary/multiple systems exhibit significantly higher density and Mach number compared to those hosting single protostars, while there are no correlations between the occurrence of fragmentation and the energy ratios of turbulence and magnetic field to gravity. Our results suggest that the higher density and supersonic turbulence of the dense cores can lead to local collapse and fragmentation to form binary/multiple systems, while the magnetic field has limited influence on fragmentation in the dense cores in Orion A.

Molecular emission is a powerful tool for studying the physical and chemical structures of dense cores. The distribution and abundance of different molecules provide information on the chemical composition and physical properties in these cores. We study the chemical segregation of three molecules (c-C$_3$H$_2$, CH$_3$OH, CH$_3$CCH) in the starless cores B68 and L1521E, and the prestellar core L1544. We applied the density-based clustering algorithms DBSCAN and HDBSCAN to identify chemical and physical structures within these cores. To enable cross-core comparisons, the input samples were characterised based on their physical environment, discarding the 2D spatial information. The clustering analysis showed significant chemical differentiation across the cores, successfully reproducing the known molecular segregation of c-C$_3$H$_2$ and CH$_3$OH in all three cores. Furthermore, it identifies a segregation between c-C$_3$H$_2$ and CH$_3$CCH, which is not apparent from the emission maps. Key features driving the clustering are integrated intensity, velocity offset, H$_2$ column density, and H$_2$ column density gradient. Different environmental conditions are reflected in the variations in the feature relevance across the cores. This study shows that density-based clustering provides valuable insights into chemical and physical structures of starless cores. It demonstrates that already small datasets of two or three molecules can yield meaningful results. This new approach revealed similarities in the clustering patterns of CH$_3$OH and CH$_3$CCH relative to c-C$_3$H$_2$, suggesting that c-C$_3$H$_2$ traces regions of lower density than to the other two molecules. This allowed for insight into the CH$_3$CCH peak in L1544, which appears to trace a landing point of chemically fresh gas that is accreted to the core, highlighting the impact of accretion processes on molecular distributions.

During a hypervelocity impact, both the impactor and target materials evaporate, generating an impact vapor plume with temperatures reaching several thousand K. As the plume cools through adiabatic expansion, chemical reactions are predicted to quench, leading to a non-equilibrium composition. However, it is still unclear how chemical reactions proceed during the cooling impact vapor plume and lead to the synthesis of organic molecules. In this study, to investigate the evolution of chemical composition within impact vapor plumes, we conducted a Monte Carlo chemical reaction simulation for complex organic synthesis, developed in our previous work. Our model does not rely on a predefined reaction network; instead, it utilizes imposed conditions for chemical changes and an approximate method for calculating reaction rates suited to our objectives. Additionally, we developed a new approach to couple these chemical reaction calculations with the rapid temperature and pressure decay in the vapor plume. Results show diverse organic molecule production depending on the impactor materials assumed in this study. These products include important precursors to biomolecules such as amino acids, sugars, and nucleobases. On the other hand, for all impactor compositions, the abundance of biomolecules themselves remains extremely low throughout the reactions from an impact to quenching. Therefore, our results suggest that biomolecules are not directly produced in impact vapor plumes but rather synthesized through reactions of these precursor molecules in aqueous solutions, following H2O condensation as the vapor plume cools. Many of the detected organic compounds, including the precursor molecules such as imine compounds and formamide, are not included in the reaction networks of previous kinetic model simulations, and their formation has not been predicted.

Understanding the dynamics of low angular momentum accretion flow around black holes (BHs) is essential for probing extreme plasma behavior in strong gravity, where shock formation can naturally produce variability signatures. In this paper, we perform general relativistic magnetohydrodynamic (GRMHD) simulations of low angular momentum accretion flows onto a BH with different BH spins to investigate the accretion dynamics near the central BH region. The simulation results show the standard and normal evolution (SANE) regime in all cases. In particular, we report the formation and persistence of standing shocks in low-angular-momentum accretion flows using multi (two and three)-dimensional GRMHD simulations for the first time. Previous studies did not detect such stable standing shock structures, making our findings a significant advancement in this field. The finding of shock dynamics can be further associated with some radiation features, such as flares observed in Sgr~A$^\ast$ and quasi-periodic oscillation (QPO) signals detected in some XRBs and AGNs.

Convective flow in Earth's iron-rich liquid core drives self-sustained dynamo action, generating Earth's magnetic field, which is strongest among all terrestrial planets of the solar system. Rock records show that this magnetic field has been operative in Earth for at least 3.4 billion years (b.y). However, advanced high pressure experiments have revised the value of the thermal conductivity of the outer core, which implies an age for the inner core of less than 1 b.y., when compositional convection begins. This creates a puzzle, with a gap between the observations of an early magnetic field on Earth and the young inner core. Previous work has suggested that the pre-inner core dynamo could have been generated in a magma ocean (MO) at the base of the mantle; however, the fluid dynamics of this scenario have received little attention. Here we numerically model the non-magnetic rotating flow in a MO above a convectively stable core in a configuration representing the pre-inner core days of Earth's evolution. Simulations here explore the importance of several dimensionless parameters on coupled core-MO convection -- the Rayleigh number, the ocean/core thermal diffusivity ratio, thermal expansion coefficient ratio, viscosity ratio, and layer thickness ratio. It is found that the MO can easily drive a flow of comparable magnitude in the core, and an approximately linear relationship is observed between the ratio of root-mean-square velocities in the core and the ocean, $(u_c^{RMS}/u_o^{RMS})$, and $(\Nu_o-1)$, where $\Nu_o$ is the Nusselt number for the MO, for the $\Nu_o$ of order 1 to 10 considered. Radial and azimuthal components of the core flow are of similar magnitude, so that, with comparable toroidal and poloidal components, we speculate that the MO-driven core flow could drive an early dynamo.

Wavelength calibration is a key factor for high-resolution spectroscopic measurements for precision radial velocities. Hollow-cathode lamps (e.g., ThAr), absorption cells (e.g., iodine cell), dielectric coated Fabry-Pérot etalons and laser frequency combs have been implemented over the years for precise wavelength calibration and wavelength drift measurements. However, due to their various impediments as wavelength calibrators, investigations of alternative methods remain of prime interest. In this paper, we examined the feasibility of low-cost (~ $1000) commercially available solid fused silica etalon with a broadband metallic coating as a calibrator. We studied the behaviour for two cavity spacings (free spectral range of 1/cm and 0.5/cm) with temperature from theoretical derivation and experimental data. Our setup had a temperature stability of 0.8 mK for a calibrator system using an off-the-shelf dewar flask with active stabilisation. Our result from radial velocity drift measurements demonstrated that such a calibration system is capable of providing higher signal-to-noise calibration and better nightly drift measurement relative to ThAr in the wavelength range between 470 nm and 780 nm. A similar result has been previously found for Fabry-Pérot etalons, and although the metalon solution lacks the efficiency of an etalon, it does offers a cost-effective broadband solution, which should be less prone to aging relative to complex dielectric mirror coatings. Nonetheless, long-term monitoring is required to understand the metalon behaviour in detail.

We present a scenario that could explain non-thermal particle acceleration in relativistic quasi-perpendicular electron-positron shocks, such as the termination shock of pulsar wind nebulae. The shock produces a strong electromagnetic precursor that propagates into the upstream plasma, which is initially threaded by a uniform background magnetic field. We show that the filamentation instability breaks the precursor into radiation filaments parallel to the shock normal. The transverse scale of the filaments is of the order of a few plasma skin depths. In the shock frame, the bulk Lorentz factor of the upstream plasma is significantly reduced inside the radiation filaments. Then, the instability produces a relativistic shear flow with strong velocity gradients on kinetic scales. The velocity gradients distort the background magnetic field lines, and generate a magnetic field component parallel to the shock normal that reverses across each radiation filament, a configuration that could trigger magnetic reconnection in the upstream plasma. These effects may accelerate particles before the plasma enters the shock.

Runaway stars are defined as stars that depart from their birth clusters at high peculiar velocities. There are two main mechanisms for the formation of runaway stars, i.e., the binary-supernova scenario (BSS) and the dynamical ejection scenario (DES). Investigating the binary fraction of runaway stars is an important step in further exploring the relative significance of the two mechanisms. We analyzed the binary fraction of 203 Galactic B-type runaway stars identified in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope Data Release 8 database. Our analysis of radial velocity variations in the runaway star sample reveals an observed spectroscopic binary fraction of $5.4\%\pm 1.6\%$, representing the proportion of objects that exhibit statistically significant variations in radial velocity with amplitudes larger than $\rm 16~km~s^{-1}$. We employed a Monte Carlo method to correct for observational biases and determined an intrinsic binary fraction of $27\%\pm 8\%$. The period and mass ratio distributions that best reproduce the observation are $f(P)\propto P^{-5.7}$ for $1\leq P\leq 1000$ days, and $f(q)\propto q^{-3.6}$ for $0.1\leq q\leq 1.0$, indicating a preference for binaries with shorter periods and less massive companions compared to a uniform distribution. The intrinsic binary fraction, in combination with previous studies on the binary fractions of runaway stars formed by the BSS and the DES, implies that both scenarios contribute comparably to the formation of Galactic B-type runaway stars, where the ratio of the BSS to the DES is 0.86.

We revisit the evolution of cosmic acceleration in a spatially flat $w_0w_a$CDM universe, in which the equation of state of dark energy takes the CPL parametrization, using the latest baryon acoustic oscillation (BAO) measurements from the Dark Energy Spectroscopic Instrument (DESI), in combination with Planck cosmic microwave background (CMB) data and several type Ia supernova datasets, including PantheonPlus, Union3, and DESY5. We analyze the deceleration parameter $q(z)$ and the jerk parameter $j(z)$ and further validate our results using the $Om(z)$ diagnostic. Our findings indicate significant deviations from the predictions of the $\Lambda$CDM model. Specifically, DESI BAO, DESI BAO + CMB, DESI BAO + CMB + Union3, and DESI BAO + CMB + DESY5 all provide strong evidence for a slowing down of cosmic acceleration at late times, as indicated by $j(0) < 0$ at more than 1$\sigma$ confidence level, within the framework of $w_0w_a$CDM model. These results suggest that in the $w_0w_a$CDM universe cosmic acceleration has already peaked and is now in a phase of decline.

The MAGIC telescopes, located at Observatorio El Roque de los Muchachos (La Palma, Spain) are two Imaging Air Cherenkov Telescopes observing the Very High Energy (VHE) gamma rays. They are run by an international collaboration composed of over 40 institutions from 12 countries. The first telescope was inaugurated in October 2003. The commissioning of the second finished in 2008. The MAGIC telescopes were designed to lower the energies to which ground based telescopes had access as well as to be able to point to any direction in the sky in less than 25 seconds. The former required the large reflective surface of 17 meters as well as an effort to optimise the mirror reflectivity and photo sensor sensitivity. The latter was achieved by minimising the weight of the full instrument using for instance carbon fibre reinforced plastic tubes for the mirror frame. The sensitivity of the MAGIC telescopes have been improving over the years thanks to hardware upgrades as well as new analysis techniques, which allowed the collaboration to keep a rich scientific program. The discovery of VHE emission from Gamma Ray Bursts and pulsars have called for a revision of the models that explain the production of gamma rays there. Both the observation of sources in flaring state as well as a systematic monitoring of sources have provided valuable data to better understand astrophysical sources both in our Galaxy and outside it. Relevant constraints on fundamental quantities like dark matter cross-section, quantum gravity scale and density of extragalactic background light have also been extracted from the observations.

We investigate the connection between stellar mass distribution, assembly history, and star formation timescales in low-redshift massive early-type galaxies (ETGs) by combining deep LegacySurvey imaging with MaNGA's spatially resolved spectroscopy. Focusing on stellar population properties, especially the [Mg/Fe] abundance ratio, we analyze stacked spectra using both absorption line indices and full-spectrum fitting. We find that, among massive ETGs with identical average stellar mass distributions beyond 5 kpc, those with higher central velocity dispersion ($\sigma_{cen}$) are older and more $\alpha$-enhanced, suggesting a connection between the in-situ star formation in the past and the central gravitational potential today for massive ETGs with a similar stellar accretion history. Conversely, at fixed $\sigma_{cen}$ and total stellar mass, galaxies with more extended stellar halos show lower [Fe/H], higher [Mg/Fe], and older ages, indicating an intriguing link between early starburst and quenching and later ex-situ assembly. These results demonstrate that the evolution of massive galaxies cannot be fully described by simple scaling relations alone, as the interplay between in-situ star formation and ex-situ accretion leaves distinct imprints in both their inner and outer stellar populations. Our findings highlight the importance of extending stellar population studies to large radii and underscore the scientific potential of next-generation IFU surveys and deep, high-resolution spectroscopy for probing the galaxy-halo connection.

Cosmic microwave background (CMB) photons are deflected by large-scale structure through gravitational lensing. This secondary effect introduces higher-order correlations in CMB anisotropies, which are used to reconstruct lensing deflections. This allows mapping of the integrated matter distribution along the line of sight, probing the growth of structure, and recovering an undistorted view of the last-scattering surface. Gravitational lensing has been measured by previous CMB experiments, with $\textit{Planck}$'s $42\,\sigma$ detection being the current best full-sky lensing map. We present an enhanced $\textit{LiteBIRD}$ lensing map by extending the CMB multipole range and including the minimum-variance estimation, leading to a $49$ to $58\,\sigma$ detection over $80\,\%$ of the sky, depending on the final complexity of polarized Galactic emission. The combination of $\textit{Planck}$ and $\textit{LiteBIRD}$ will be the best full-sky lensing map in the 2030s, providing a $72$ to $78\,\sigma$ detection over $80\,\%$ of the sky, almost doubling $\textit{Planck}$'s sensitivity. Finally, we explore different applications of the lensing map, including cosmological parameter estimation using a lensing-only likelihood and internal delensing, showing that the combination of both experiments leads to improved constraints. The combination of $\textit{Planck}$ + $\textit{LiteBIRD}$ will improve the $S_8$ constraint by a factor of 2 compared to $\textit{Planck}$, and $\textit{Planck}$ + $\textit{LiteBIRD}$ internal delensing will improve $\textit{LiteBIRD}$'s tensor-to-scalar ratio constraint by $6\,\%$. We have tested the robustness of our results against foreground models of different complexity, showing that a significant improvement remains even for the most complex foregrounds.

V Sge is a peculiar, highly luminous long-period (12.34h) binary star that can display a super-soft X-ray emitting component when in the faint phase of its V~ 10-13mag variability range. Apparently undergoing Eddington-limited accretion from its more massive secondary, it is in a very rare, short-lived evolutionary phase towards the double degenerate channel. Its complex and highly variable optical emission features, from Balmer and Heii to high-ionisation lines, including strong fluorescence features, have been challenging to interpret, especially given the absence of any absorption lines associated with photospheric features from either stellar component. With the detailed properties of V Sge, especially the donor, still controversial, we undertook a VLT/X-Shooter campaign over three months in 2023, obtaining high S/N, high resolution spectra that revealed multiple components in both high- and low-ionisation lines. This allows us to track V Sge's principal emitting regions via Doppler tomography, obtaining new insights into high accretion-rate dynamics. In particular, we identify a stationary, double-peaked emission core which we interpret as a circumbinary ring, analogous to SS433. This enables us to derive limits on the system masses. Furthermore, we find very broad emission-line wings whose mean velocity can vary over hundreds of kilometres per second on timescales of decades, yet ``flip'' between states in <1 week. We show that the super-soft X-ray source interpretation is able to account for these and other observational attributes significantly better than the hot binary model, concluding that V Sge could be one of the brightest known Galactic super-soft sources.

Understanding the provenance of CI and CM chondrites, the most primitive materials in our meteorite collections, is critical for shedding light on the Solar System's early evolution and contextualizing findings from recent sample return missions. Here we show that the parent bodies of CM chondrites originate from the Saturn formation region, whereas those of CI chondrites originate essentially from the primordial trans-Uranian disk. Using Nbody simulations to investigate the effect of giant planet growth and inward Type-I migration, along with the current observed distribution of CM, CI, and comet-like P types bodies in the asteroid belt, we demonstrate that CI- and CM-like bodies must have been implanted at different times in the belt. In contrast, CI and comet-like bodies were implanted at the same time. These different implantation periods are imposed by the fact that the gas disk profile entirely governs the radial distribution of bodies implanted by aerodynamic drag in the asteroid belt. A preferred location coincides with the inner edge of a gap opened by Jupiter. Saturn's growth likely drove the migration of CM-like bodies. In contrast, CI and comet-like bodies were transported at a later stage during the outward migration of Uranus and Neptune, driven by remaining planetesimals. Since CM chondrites are chondrule-rich, it follows that chondrule formation occurred mostly inward of the ice giants formation zone (under 10 au). A byproduct of our simulations is that only CM-like, not CI-like, bodies contributed to the water budget of the telluric planets.

Modern large scale cosmological hydrodynamic simulations require robust tools capable of analysing their data outputs in a parallel and efficient manner. We introduce SOAP (Spherical Overdensity and Aperture Processor), a Python package designed to compute halo and galaxy properties from SWIFT simulations after being post-processed with a subhalo finder. SOAP takes a subhalo catalogue as input and calculates a wide array of properties for each object. SOAP offers parallel processing capabilities via mpi4py for efficient handling of large datasets, and allows for consistent property calculation across multiple halo finders. SOAP supports various halo definitions, including spherical overdensities and fixed physical apertures, providing flexibility for diverse observational comparisons. The package is compatible with both dark matter-only and full hydrodynamic simulations, producing HDF5 catalogues that are integrated with the swiftsimio package for seamless unit handling.

For over 30 years, the Magneto-Rotational Instability has been accepted as the mechanism driving accretion disk turbulence. Its physical basis is well understood, where an interplay between centrifugal forces and magnetic tension transfers angular momentum between oppositely displaced fluid elements. In this work, we revisit this picture in global disk models and various magnetic field topologies and generalise it to non-axisymmetric instabilities like the Super-Alfvénic Rotational Instability (SARI). We use the open-source \texttt{Legolas} software to quantify all complex-valued linear eigenfunctions for the (near-)eigenmodes and visualise the resulting spatio-temporal variations in real space in a manner that can be compared to direct numerical simulations of disks. The field perturbations are fundamentally different between the (axisymmetric) MRI and the novel, ultra-localised SARI modes, which bear some resemblance to spiral modes in galaxies but differ in important ways. We use a combined numerical-analytical approach to study the polarization of the magnetic and velocity field perturbations. Finally, we compare disks of differing magnetic topology where many linear modes are merely superposed to recreate a visual impression of `turbulent' fields and quantify the resulting stresses. We find that even superposed, still linearly growing SARI modes can already provide the needed effective viscosity-related alpha-values invoked for angular momentum transport. 3D views on the magnetic field perturbations show that SARI modes of opposite azimuthal mode number may naturally introduce plasmoid and toroidal flux-tube like field deformations.

Primordial gravitational waves on small scales are not tightly constrained by current cosmological observations, which allows for the possibility of large amplitudes at small scales. We investigate second-order tensor induced gravitational waves (TIGWs) sourced by primordial gravitational waves and present the corresponding corrections to the total energy density spectrum of gravitational wave. We analyze primordial gravitational waves with large amplitudes generated by various models at small scales. Our results indicate that when primordial gravitational waves on small scales sufficiently dominate the current PTA observations, corrections to the total energy density spectrum from second-order TIGWs may become pronounced in certain frequency bands.

The Ricochet experiment aims to measure the coherent elastic neutrino-nucleus scattering process from antineutrinos emitted by a research nuclear reactor operated by the Institut Laue-Langevin (Grenoble, France). This article presents a description of the Ricochet experimental installation and the detector performance achieved during its commissioning with a mini-CryoCube module consisting of three 42-gram germanium cryogenic calorimeters. The baseline resolutions and background levels are reported both during reactor-on and reactor-off periods, and as noise mitigation techniques were improved. A baseline resolution of 40 eV electron equivalent was achieved for the ionization channel after setup improvements, and the phonon channel resolutions ranged from 50 to 80 eV of total phonon energy. In the energy region from 2 to 7 keV, a nuclear recoil rate of 15(2) events/(kg day keV) is measured during the reactor-off period selecting events in coincidence with muon veto signals. This rate is in agreement with the cosmogenic neutron rate calculated from GEANT4 simulations. After the rejection of events in coincidence with signals in the muon veto detectors, a combined 90% C.L. limit on the nuclear recoil background of < 9 events/(kg day keV) is obtained in that energy region during the reactor-on period, which is compatible with our GEANT4 model calculation corresponding to a total rate of 5 events/(kg day keV). The sensitivity of this analysis was however found to be limited by a surface event contamination which is currently being addressed by the Ricochet Collaboration with upgraded detectors.

We propose a model for dust formation in Type II supernovae (SNe) interacting with confined circumstellar material (CSM), motivated by recent time-domain surveys that have revealed a substantial fraction of SN progenitors to be surrounded by CSM ejected shortly before core-collapse. We simulate the pre-SN mass eruption and the resulting confined CSM using the open-source code CHIPS, and follow the subsequent evolution of the SN ejecta and its interaction with the CSM. We show that a cold dense shell (CDS) is formed at the radiative shock under a wide range of conditions and later undergoes rapid adiabatic cooling during free expansion, leading to efficient dust condensation. The resulting dust mass ranges from $\sim10^{-3}\,M_\odot$ to $0.1\,M_\odot$, depending on the mass and spatial extent of the CSM. We further calculate the infrared (IR) emission from the newly formed dust and find broad consistency with observations of SN~1998S. Notably, the IR light curve exhibits a rapid rise within $\lesssim10\,{\rm d}$, closely resembling that of kilonovae (KNe). This suggests that dust emission powered by confined CSM interaction may be also discovered in KN searches. Moreover, the high-density environment of the CDS may allow dust grains to grow to larger sizes, enhancing their survivability against destruction by reverse shocks propagating from the interstellar medium at later times.

The structure of spiral galaxies is essential to understanding the dynamics and evolution of disc galaxies; however, the precise nature of spiral arms remains uncertain. Two challenges in understanding the mechanisms driving spirals are how galactic environment impacts spiral morphology and how they evolve over time. We present a catalog characterizing the properties, dynamics, and evolution of m=2 spiral structure in 10 Milky Way-mass galaxies from the FIRE-2 cosmological zoom-in simulations. Consistent with previous literature, we find that FIRE-2 spirals are transient, recurring features simultaneously present in the disc at varying pattern speeds ($\Omega_p$) that broadly decrease with radius. These spirals persist on Gyr timescales (mean duration 1.90 Gyr), but fluctuate in amplitude on timescales of hundreds of Myr. Tidal interactions and bar episodes impact the resulting m=2 spiral structure; strong satellite interactions generally produce shorter-lived, stronger spirals with larger radial extent, and bars can increase $\Omega_p$. Galactic environment influences spiral structure; kinematically colder discs can support longer-lived, stronger spirals. The properties of identified spirals in FIRE-2 vary widely in radial extent (0.3-10.8 kpc), duration (1.00-6.00 Gyr), and amplitudes ($a_{2,\text{max}}$=0.018-0.192). We find the presence of spirals in all age populations, suggesting these are density wave-driven features. This work represents the first time that spiral structure has been cataloged in this manner in cosmological simulations; the catalog can be leveraged with current and forthcoming observational surveys, enabling systematic comparisons to further our understanding of galaxy evolution.

Context: SN 2024ggi is a Type II supernova (SN) discovered in the nearby galaxy NGC 3621 (D $\approx6.7\pm0.d$ Mpc) on 2024 April 03.21 UT. Its proximity enabled a detailed investigation of the SN's properties and its progenitor star. This work focuses on the optical evolution of SN 2024ggi at the nebular phase. Aims: We investigate the progenitor properties and possible asymmetries in the ejecta by studying the nebular phase evolution between days 287 and 400 after the explosion. Methods: We present optical photometry and spectroscopy of SN 2024ggi during the nebular phase, obtained with the Las Campanas and Gemini South Observatories. Four nebular spectra were taken at 287, 288, 360, and 396 days post-explosion, supplemented by late-time $uBVgri$-band photometry spanning $320-400$ days. The analysis of the nebular emission features is performed to probe ejecta asymmetries. Based on the [O I] flux and [O I]/[Ca II] ratio, and comparisons with spectra models from the literature, we arrive to an estimate of the progenitor mass. Additionally, we construct the bolometric light curve from optical photometry and near-infrared data to derive the synthesized nickel mass. Results: Our analysis suggests a progenitor zero-age-main-sequence mass between $12-15 M_\odot$. The late-time bolometric light curve is consistent with a synthesized $^{56}$Ni mass of $0.05-0.06 M_\odot$. The line profiles exhibit only minor changes over the observed period and suggest a roughly symmetrical ejecta, with a possible clump of oxygen-rich material moving towards the observer. No signatures of circumstellar material interaction are detected up to 400 days after the explosion.

Discrete symmetry-breaking phase transitions in the early universe may have caused the formation of networks of sheet-like topological defects, usually referred to as domain walls, which separate regions that have settled into different vacuum states. Field theory simulations predict the successive collapse of increasingly larger domains, which could potentially leave observable imprints in present-day large-scale structures. We use the parameter-free velocity-dependent one-scale model to provide an estimate of the final decay energy of these walls and their associated collapse rate, as a function of redshift. The energy released by collapsing walls can act as a seed for density perturbations in the background matter field, influencing structure formation. We estimate the dependence of the current mass of the resulting non-linear objects on the collapse redshift and wall tension, showing that domain walls can contribute to the formation of objects as massive as present-day galaxy clusters. Still, we confirm that the contribution of standard domain walls to structure formation is subdominant. In contrast, biased domain walls generally face much less stringent constraints on their tension, which allows for significantly higher collapse energies. Based on our analysis, we are able to show that the collapse of such biased wall networks can provide a significant contribution to structure formation, and, in particular, a mass excess at $z \gtrsim 7$ as suggested by JWST data.

We aim to characterize the multi-phase gas in the SPT2349-56 protocluster at z=4.3, known to host one of the most starbursting and AGN-rich high redshift this http URL this purpose we conducted APEX single dish observations of the [CII]158 micron (hereafter [CII]) line towards the Core and North components, previously imaged with the ALMA 12-m array. We also present the first [OIII]88 micron (hereafter [OIII]) line observations in such high redshift protocluster system. We obtain a [CII] line luminosity $\sim$1.7$\times$ more than the one recovered by ALMA towards the Core, while remarkably we recover 4$\times$ more [CII] line emission than the one found in deep ALMA images towards the North component, suggesting that the most massive gas reservoirs lie in the less extreme regions of this protocluster system. A minimum ionised gas mass of $\mathrm M_{\rm min}(H^+)$$ \sim$$3.7\times 10^{10}$\,\Msun\, is deduced from the [OIII] line, amounting to 30\% of the molecular gas mass in the same area. Finally we obtain star formation rate (SFR) estimates using the [OIII] line luminosity, and the corresponding ionised gas mass. These yield values that can surpass the far-IR continuum-derived SFR (under the assumption of a standard stellar IMF), which can be reconciled only if non-stellar ionising sources contribute to the [OIII] line luminosity, or a top-heavy stellar IMF produces a larger fraction of O stars per total stellar mass, a distinct possibility in High-Energy-Particle (HEP) rather than (UV-photon)-dominated environments in clusters. Future work using far-IR fine-structure and molecular/neutral-atomic lines is necessary for determining the thermal/ionisation states of the multi-phase medium and these line ratios must be measured over a wide range of spatial scales, which ultimately requires combining wide-field single-dish and high resolution interferometric observations.

We derive multiple constraints on dark energy and compare dynamical dark energy models with a time-varying equation of state ($w_0 w_a$CDM) versus a cosmological constant model ($\Lambda$CDM). We use Baryon Acoustic Oscillation (BAO) from DESI and DES, Cosmic Microwave Background (Planck) with and without lensing from Planck and ACT (noted CMB and CMBL, respectively), supernova (SN), and cross-correlations between galaxy positions and galaxy lensing from DES. First, we use pairs or trios of datasets where we exclude one type of datasets each time and categorize them as ``NO SN", ``NO CMB" and ``NO BAO" combinations. In all cases, we find that the combinations favor the $w_0 w_a$CDM model over $\Lambda$CDM, with significance ranging from 2.3$\sigma$ to 3.3$\sigma$. For example, DESI+DESY6BAO+CMB yields 3.2$\sigma$ without SN, DESI+DESY6BAO+DESY5SN yields 3.3$\sigma$ without CMB, and CMB+DESY5SN+DES3x2pts yields 2.6$\sigma$ without BAO. The persistence of this pattern across various dataset combinations even when any of the dataset is excluded supports an overall validation of this trending result regardless of any specific dataset. Next, we use larger combinations of these datasets after verifying their mutual consistency within the $w_0 w_a$CDM model. We find combinations that give significance levels $\sim$4$\sigma$, with DESI+DESY6BAO+CMBL+DESY5SN reaching 4.4$\sigma$. In sum, while we need to remain prudent, the combination of the first step that supports a validation of the pattern of these results beyond any single type of datasets and their associated systematics, together with the second step showing high-significance results when such datasets are combined, presents a compelling overall portrait in favor of a dynamical dark energy with a time-evolving equation of state over a cosmological constant, and constitutes a serious challenge to the $\Lambda$CDM model's reign. [Abridged]

JWST has shed light on galaxy formation and metal enrichment within 300 Myr of the Big Bang. While luminous galaxies at $z > 10$ often show significant [O III]$\lambda\lambda$4959, 5007 emission lines, it remains unclear whether such features are prevalent among fainter, more typical galaxies due to observational limits. We present deep imaging and spectroscopy of JADES-GS-z14-1 at $z_\mathrm{spec}=13.86^{+0.04}_{-0.05}$, currently the faintest spectroscopically confirmed galaxy at $z\approx 14$. It serendipitously received 70.7 hours of MIRI/F770W imaging in the JWST Advanced Deep Extragalactic Survey (JADES), the deepest MIRI exposure for any high-redshift galaxy to date. Nonetheless, we detect only tentative F770W emission of $7.9\pm2.8$ nJy at $2.8\sigma$ significance, constraining the total equivalent width of [O III]$\lambda\lambda$4959, 5007 + H$\beta$ to $520^{+400}_{-380}$ A, weaker than most $z > 10$ galaxies with MIRI detections. This source is unresolved across 16 NIRCam bands, implying a physical radius $\lesssim50$ pc. NIRSpec/PRISM spectroscopy totaling 56 hours reveals no rest-frame ultraviolet emission lines above $3 \sigma$. Stellar population synthesis suggests a stellar mass $\sim4\times 10^{7}$ $\mathrm{M_\odot}$ and a star formation rate $\sim 2$ $\mathrm{M_\odot yr^{-1}}$. The absence of strong metal emission lines despite intense star formation suggests a gas-phase metallicity below 10% solar and potentially a high escape fraction of ionizing photons. These deep observations provide rare constraints on faint, early galaxies, tracing the onset of chemical enrichment and ionization in the early Universe.

Recent observations with the James Webb Space Telescope (JWST) have suggested the existence of over-massive black holes (OMBHs) in high-redshift galaxies. In this paper, we propose a new mechanism for the formation of OMBHs, based on the accretion of globular clusters (GCs) in compact disk galaxies. We derive the conditions under which OMBHs can form, focusing on key parameters such as halo mass, redshift, and halo spin parameter. Our results show that at redshift $z = 10$, a halo with mass $10^{11}~M_{\odot}$ and a spin parameter of $0.02$ can form a black hole of $1.4 \times 10^{8}~M_{\odot}$ through GC migration and accretion via tidal disruption events (TDEs). The resulting black hole-to-stellar mass ratio can reach $\sim 0.1$, corresponding to the fraction of GC mass accreted onto the black hole. This mechanism thus provides a plausible explanation for the OMBHs observed by JWST. Furthermore, by combining our model with halo mass functions and the distribution of spin parameters, we construct black hole mass functions that successfully reproduce the number densities of massive BH candidates at $z \sim 5$ inferred from JWST observations, and UHZ1 and GHZ9 at $z \sim 10$.

The Vera C. Rubin Observatory will soon survey the southern sky, delivering a depth and sky coverage that is unprecedented in time domain astronomy. As part of commissioning, Data Preview 1 (DP1) has been released. It comprises a ComCam observing campaign between November and December 2024 with multi-band imaging of seven fields, covering roughly 0.4 square degree each, provides a first glimpse into the data products that will become available once the Legacy Survey of Space and Time begins. In this work, we search three fields for extragalactic transients. We identify six new extragalactic transients, and three known ones from a sample of 369,644 difference image analysis objects. Photometric classification using \texttt{Superphot+} indicates that this sample likely comprises six type Ia, two type II, two type Ibc and one type IIn supernovae. Our findings are in slight tension with supernova detection rate predictions from the literature of $12\pm3$ SN Ia and $3\pm1$ core-collapse supernovae likely due to the lack of suitable templates. Nevertheless, this work demonstrates the quality of the data products delivered in DP1 and indicates that Rubin Observatory Legacy Survey and Space and Time (LSST) is well placed to fulfill its discovery potential in time domain astronomy.

The unexpectedly high abundance of galaxies at $z > 11$ revealed by JWST has sparked a debate on the nature of early galaxies and the physical mechanisms regulating their formation. The Atacama Large Millimeter/submillimeter Array (ALMA) has begun to provide vital insights on their gas and dust content, but so far only for extreme 'blue monsters'. Here we present new, deep ALMA observations of JADES-GS-z11-0, a more typical (sub-$L^*$) $z > 11$ galaxy that bridges the discovery space of JWST and the Hubble Space Telescope. These data confirm the presence of the [O III] 88 $\mu$m line at $4.5\sigma$ significance, precisely at the redshift of several faint emission lines previously seen with JWST/NIRSpec, while the underlying dust continuum remains undetected ($F_\nu < 9.0 \, \mathrm{\mu Jy}$), implying an obscured star formation rate (SFR) of $\text{SFR}_\text{IR} \lesssim 6 \, \mathrm{M_\odot \, yr^{-1}}$ and dust mass of $M_\text{dust} \lesssim 1.0 \times 10^{6} \, \mathrm{M_\odot}$ (all $3\sigma$). The accurate ALMA redshift of $z_\text{[O III]} = 11.1221 \pm 0.0006$ ($\gtrsim \! 5\times$ refined over NIRSpec) helps confirm that redshifts measured purely from the Lyman-$\alpha$ break, even spectroscopically, should properly take into account the effects of potential damped Lyman-$\alpha$ absorption (DLA) systems to avoid systematic overestimates of up to $\Delta z \approx 0.5$. The [O III] 88 $\mu$m luminosity of $L_\text{[O III]} = (1.0 \pm 0.3) \times 10^{8} \, \mathrm{L_\odot}$, meanwhile, agrees well with the scaling relation for local metal-poor dwarfs given the SFR measured by NIRCam, NIRSpec, and MIRI. The spatially resolved MIRI and ALMA emission also underscores that JADES-GS-z11-0 is likely to consist of two low-mass components that are undergoing strong bursts of star formation yet are already pre-enriched in oxygen (~30% solar), only 400 Myr after the Big Bang.

We investigate the cosmological implications of a novel definition of field theory vacuum energy. The free field Hamiltonian represented as an ensemble of oscillators (in the Fourier space) usually implies the presence of mass scale for these oscillators, which in quantum field theory is of little importance since quantum energy spectrum of oscillator is mass independent. This mass scale, however, may be interesting due to its possible gravitational implications. Since black hole physics puts an upper limit on the total energy within a given region, one obtains constraint on the number of field oscillators. If the mass scale for field oscillators is set by the IR cutoff, then this number saturates the black hole entropy bound. Following this reasoning, one derives various kinds of dark energy models that maybe interesting for further study.

Supernova cooling provides a powerful probe of physics beyond the Standard Model (SM), in particular for new, light states interacting feebly with SM particles. In this work, we investigate for the first time the production of fermionic dark matter (DM) via the neutrino-devouring process inside a core-collapse supernova, which contributes to the excessive cooling. By incorporating state-of-the-art supernova simulation data and the full time evolution information, we derive stringent and robust limits on DM interactions. We exclude the cross sections down to $10^{-51}-10^{-58}$ cm$^2$ in the keV-MeV mass range for DM-electron scattering, and $10^{-49}-10^{-56}$ cm$^2$ in the 0.1-100 MeV mass range for DM-nucleon scattering, supplemented by complementary constraints from cosmology, astrophysics, LHC and direct detection experiments in the larger cross section regime. We also close almost the entire window in which fermionic DM constitutes $\mathcal{O}(1)$ fraction of DM for its coupling to electrons in the keV-MeV mass range.

We present a unified post-Newtonian framework for relativistic timing and coordinate transformations covering six time scales (TCB, TCG, TT, TDB, TCL, TL) and three reference systems (BCRS, GCRS, LCRS). Extending the IAU conventions, we define a Lunicentric Celestial Reference System (LCRS) metric that retains all contributions above a fractional threshold of 5x10^{-18} and timing terms above 0.1 ps by expanding the lunar gravity field to spherical-harmonic degree l=9 with Love number variations and including external tidal and inertial multipoles to the octupole. We derive closed-form mappings among TCB, TCG, TT, TCL and TL, yielding proper-to-coordinate time transformations and two-way time-transfer corrections at sub-picosecond accuracy. We evaluate secular rate constants and periodic perturbations arising from kinematic dilation, lunar monopole and multipoles, Earth tides and gravitomagnetic effects for clocks on the lunar surface, in low lunar orbits, at the Earth-Moon L1 point and in near-rectilinear halo orbits. Our analysis demonstrates that harmonics through l=9 and tides through l=8 are required to achieve 5x10^{-18} fractional stability, supporting sub-picosecond clock synchronization and centimeter-level navigation in cislunar space. This framework underpins high-precision time and frequency transfer, relativistic geodesy, quantum communication links and fundamental physics experiments beyond low Earth orbit.

In this work, we investigate geodesics and black hole shadows in the Kerr-Bertotti-Robinson spacetime. We show that the equations of motion for null geodesics are separable and admit analytical treatment, whereas timelike geodesics are generally non-separable. Approximate analytical expressions for the photon sphere and the innermost stable circular orbit are derived via perturbative expansions in the magnetic field strength. We further explore the black hole shadow using both numerical and analytical methods, examining the effects of the magnetic field, the observer's inclination angle and radial position. Deviations from the standard Kerr shadow are quantified, and a physical interpretation is provided by introducing asymptotic regimes defined relative to the magnetic field strength.

Inflationary tensor perturbations are treated as arising from a bath of gravitons produced by quantum particle creation at the end of inflation. We calculate the correlation function of the CMB temperature fluctuations produced by these gravitons in a model with an infrared cut off. The CMB photons are emitted from within a last scattering shell of finite thickness in redshift. We find the correlation function in terms of the separation of a pair of spacetime points of emission in both angle and redshift. In both variables, there is a significant amount of anti-correlation. The anti-correlation minimum has a relative magnitude compared to the central correlation maximum of about 20% in angle and 50% in redshift.

We have analyzed LVK gravitational wave events that show some evidence of eccentricity from TEOBResumS modeling parameter estimations and have confronted them independently with full numerical generated waveforms from our bank of nearly two thousand simulations of binary black holes. We have used RIFT for Bayesian parameter estimation and found that GW200208\_{22} KDE estimates favor eccentricities $e_{20} = 0.217_{-0.184}^{+0.076}$ upon entering the LVK band at $\sim20$Hz within a 90\% confidence limit. Within this event analysis we employed 39 new targeted full numerical relativity simulations and we have thus found a top improved likelihood $\ln \mathcal{L}$ matching waveform, compared to model-based analysis, with an estimated eccentricity at 20Hz, $e_{20}=0.200$, thus reinforcing the eccentric hypothesis of the binary. We have also used our full bank of numerical waveforms on GW190620 finding that it favors eccentricities in the range of {$0\leq e_{10}\leq0.3$}. New specifically targeted simulations will be required to narrow this eccentricity range.

Optical beamsplitters with similar properties for orthogonal, linear polarisation modes are required for realising polarisation-based speedmeter schemes to reduce back-action noise in gravitational-wave interferometers. In this paper, we investigate two beamsplitter coatings obtained from Laseroptik GmbH and Optoman on a best-effort basis that aim for a 50/50 power splitting ratio and equal overall phase shift for two orthogonal, linear polarisation modes interacting with the optic. We show that while Laseroptik GmbH opted for coating stack with 22 alternating layers of Ta2O5 and SiO2, Optoman produced a much thinner coating made of 5 SiO2 and SiOx (0 < x < 2) layers. With these strategies, the Laseroptik coating achieves an equal power reflectivity of 51% at 46 deg angle of incidence, and zero phase shift between both polarisations at 44.25 deg angle of incidence. The Optoman coating achieves power reflectivities of 49% for s-polarisation and 51% for p-polarisation with a differential phase shift around 5 deg largely independent of the angle of incidence.

Fermionic asymmetric dark matter (ADM) can be captured in neutron stars (NSs) via scatterings with the star constituents. The absence of dark matter annihilation due to its asymmetric nature leads to ADM accumulation in the NS core, potentially reaching densities sufficient to exceed the Chandrasekhar limit and trigger its gravitational collapse into a black hole (BH), eventually consuming the NS from within. Therefore, the existence and observation of old neutron stars provides a means to constrain the properties of ADM. We revisit previous constraints on the mass and scattering cross section off neutrons of fermionic ADM across a class of models. We critically examine common simplifying approximations used in the literature to derive these limits. Our analysis includes improved treatments of dark matter capture, thermalization, BH formation, accretion, and evaporation. We find that previous results can be relaxed by a few orders of magnitude once these effects are properly accounted for.