Papers for Tuesday, Mar 17 2026

Quasi-periodic fast-propagating (QFP) wave trains are a distinctive form of magnetohydrodynamic disturbance frequently observed in the solar corona. Yet their excitation mechanism and propagation characteristics are not well understood. In this study, we investigate a well-observed QFP wave event associated with an M6.5-class flare and coronal mass ejection that occurred in active region (AR) 12371 on 2015 June 22 by combining multi-wavelength observations from SDO/AIA and HMI with data-inspired 3D MHD simulations. The QFP wave trains propagating at high speeds of 1140$-$1760 km~s$^{-1}$ are detected in the AIA 171 Å channel, following global EUV wave fronts visible at 171 and 193~Å traveling at considerably lower speeds of 300$-$510 km~s$^{-1}$. Wavelet analysis reveals consistent 2--4 minutes periodicities in both the QFPs and flare quasi-periodic pulsations (QPPs) observed in UV/EUV and hard X-ray emissions, suggesting a common origin likely linked to intermittent magnetic reconnection. Guided by these observations, we construct realistic 3D MHD models incorporating dense fan-loop structures and periodic drivers applied at different locations. The simulations reproduce the key characteristics of the observed wave trains. Comparison between cases with and without a coronal background (non-loop plasma emission) indicates that coronal density structuring significantly modifies the detected wave amplitudes and propagation patterns. Our results highlight the importance of realistic coronal magnetic configurations in modeling QFP dynamics and suggest that their observed association with fan loops in AIA 171 Å may represent a temperature-dependent visibility effect rather than a genuine confinement of the waves.

Transmission spectroscopy has greatly advanced the study of exoplanet atmospheres, but stellar surface heterogeneities can contaminate transit spectra. Characterising how stellar spectra vary across the stellar disc is therefore essential to disentangle stellar and planetary contributions. Transit observations can probe the local stellar spectra along the planet's transit chord. We study centre-to-limb variations of line profiles across the surface of HD 189733 using the ESPRESSO spectrograph. Building on previous work, we assess the feasibility of applying the Doppler shadow technique with ESPRESSO and compare the results with solar observations and numerical simulations. We analyse spectra obtained during two transits of HD 189733 b. Each spectrum was cross-correlated with two masks of selected Fe I lines, producing four sets of cross-correlation functions (CCFs). Using a Doppler shadow methodology, we retrieved local stellar profiles along the transit chord. These were compared with previous studies, with disc-resolved solar spectra from IAG ATLAS and with transit simulations generated using SOAPv4 and synthetic spectra from Turbospectrum based on MARCS stellar atmosphere models under LTE and NLTE conditions. For three Fe I CCF sets we detect a statistically significant increase in line depth from stellar limb to centre, consistent with Turbospectrum predictions, although solar data show a weaker gradient. For one CCF set we also find that line widths decrease from limb to centre, consistent with solar observations but not reproduced by the simulations. These results demonstrate the capability of ESPRESSO to measure centre-to-limb variations of spectral line profiles on other stars. While the local CCF profiles of HD 189733 agree with solar data, discrepancies in line widths suggest that additional physical processes are required to reproduce the observed profiles.

The discovery of galaxies with extremely strong nebular continuum emission at high redshifts provide novel, unique insights into the conditions under which the first super-massive stars formed. Here we identify a galaxy at redshift $z=5.124$ observed by the JWST CAPERS survey that exhibits a prominent turnover in the rest-frame UV continuum and a pronounced Balmer `jump'. We model the entire JWST/NIRSpec Prism spectrum from rest-frame UV to optical wavelength, finding that a dominant ($>95\%$) nebular continuum emission can accurately reproduce the spectral shape across all wavelengths. We tested an alternative model with strong damped Ly$\alpha$ absorption (DLA), but found that it is not able to match the shape of the turnover without invoking a large freedom in the redshift of the absorber. The nebular continuum emission model reveals a hot ($T = (5.3\pm 0.2)\times 10^{4}$ K) and dense ($n_e = (5.4\pm 0.8)\times 10^{3} {\rm cm^{-3}}$) nebular region powering the origin of the spectral shape. We also note the presence of a `blue' companion source at the same redshift, offset by 3 kpc to the main galaxy. Intriguingly, the spectrum of this source show several hints of hydrogen and helium lines, but no metal lines are detected. We theorize that this companion galaxy might be comprised mainly of Population III (PopIII) stellar remnants and potentially powers the nebular continuum emission seen in the main galaxy. These results have important implications for the presence of a potential dominant population of super-massive and PopIII stars and their consequent excess UV brightness for a significant fraction of galaxies at cosmic dawn.

Dissipative effects appear in many early-Universe scenarios, yet their universal observational signatures and systematic confrontation with data remain largely unexplored. We employ the Open Effective Field Theory of Inflation (Open EFToI) to consistently incorporate dissipative and stochastic effects while preserving scale invariance. Dissipation enhances specific interaction channels of the Goldstone mode, generating distinctive primordial non-Gaussian signatures, beyond those generically produced by standard EFToI. In the weak-dissipation regime, this includes folded bispectrum shapes observationally more favoured than both the equilateral and orthogonal templates. Using the Modal bispectrum pipeline with the Planck CMB data, we obtain the likelihood and derive the first model-independent bounds on early-Universe dissipation. We find a marginalised upper bound on the dissipation scale $\gamma \leq 384\,H$ and a lower bound on the sound speed $c_s \geq 0.38$ at $95\%$ confidence level. The maximum likelihood for best-fit models reveals a degeneracy between $\gamma$ and $c_s$. These results open a model-independent window for probing departures from minimal inflation and discriminating between early-Universe scenarios with stochastic noise and dissipative effects.

Hydrodynamic simulations of neutron star mergers rely on the clear separation between the strong-interaction, weak-interaction, and hydrodynamic timescales. In this effective framework, weak Urca interactions are typically the slowest microscopic processes, and therefore the Urca rate determines the bulk-viscous dissipation. This assumed hierarchy of dissipative mechanisms can be decisively altered, without invalidating hydrodynamics, if the trajectory of the matter in a neutron star merger passes through the vicinity of a possible low temperature QCD critical point. The enhanced density fluctuations lead to critical slowing down and rapid growth of transport coefficients including bulk viscosity. While this growth is regulated by finite-time effects, finite-size effects, and the breakdown of hydrodynamic scale separation, which bound the correlation length, we demonstrate that the QCD contribution to bulk viscosity can rival the electroweak contribution in realistic conditions. Thus, critical dynamics could leave observable imprints on the hydrodynamic evolution of neutron star mergers.

A significant neutrino event with an estimated energy between $72\,\mathrm{PeV}$ and $2.6\,\mathrm{EeV}$ was recently observed by the KM3NeT experiment (KM3-230213A). When interpreted as cosmogenic in origin, this event can provide constraints on several phenomenological parameters of UHE proton sources. In this study, we present the best fit to the spectrum and composition of UHECRs that is consistent with multi-messenger constraints, including the detection of a single neutrino event by the KM3NeT detector in the energy range of KM3-230213A. From the best fit, we obtain the 68\% CL constraints on the parameters of a two-population model of UHECRs, comprising a mixed-composition population and a subdominant UHE proton population. Our results indicate that the detection of a single neutrino event in the energy range of KM3-230213A solely with the KM3NeT exposure requires strongly evolving UHE proton sources, consistent with high-luminosity active galactic nuclei. On the other hand, including the null observations from the Pierre Auger and IceCube observatories disfavors such strong evolution. In both cases, the observed proton fraction of UHECRs is primarily constrained by the composition data to be $\sim 20\%$ at $20\,\mathrm{EeV}$.

Key absorption features in the spectra of massive quiescent galaxies reveal a wealth of information about their stellar populations and, in some cases, the properties of gas within and around them. With the spectroscopic capabilities of \textit{JWST}, we are now able to probe deeper and farther into the near-infrared than ever before. It is therefore crucial that we fully understand the origins of observed spectral absorption and emission features. The goal of this document is to collate important rest-frame optical to near-infrared (NIR) spectral features of quiescent galaxies in the context of their physical origin (e.g., stellar photospheres, the interstellar medium, or a combination thereof). We present a look-up table summarizing key information, including a ``diagnostic'' column indicating whether lines are most sensitive to age, metallicity, surface gravity (IMF sensitivity), $\alpha$-enhancement, electron temperature or gas density. This compilation is intended to serve as a practical reference for interpreting rest-optical-NIR spectral features in quiescent galaxies, particularly as \textit{JWST} spectroscopy enables increasingly detailed studies of their stellar populations and surrounding gas.

Kilonovae from compact binary mergers are most informative when accompanied by a gravitational-wave signal, which can help identify the source as a binary neutron star (BNS) or a neutron star-black hole (NSBH) merger. However, future events will also be discovered serendipitously or through follow-up of other transients, without a confident identification of the progenitor. Hence, we ask whether the kilonova light curve alone can distinguish between these two progenitor channels. Using simulated BNS and NSBH populations together with semi-analytic light curve models, we compare their post-peak evolution across the optical $ugrizy$ bands. The strongest contrast appears in the blue $u$ band 2 days after peak and in the redder $i$ band 10 days after peak. In the $u$ band, typical BNS kilonovae decline by only $\sim 1$ mag within 2 days of peak, whereas NSBH kilonovae typically decline by $\gtrsim 3$ mag over the same interval. In the $i$ band, the trend reverses for most of the population, with NSBH kilonovae evolving more slowly than BNS kilonovae. We attribute this behavior to differences in ejecta mass, opacity, and diffusion timescale between the two merger classes. Although the quantitative overlap is model-dependent, the qualitative distinction persists across model variations, identifying post-peak decline to be a viable diagnostic for inferring whether the source was a BNS or an NSBH merger.

Multi-field models of inflation typically assume that interactions between particles can be treated perturbatively. Strongly-coupled models provide an intriguing alternative and may offer novel inflationary phenomenology. We study the "unparticle" scenario, where the inflaton is weakly mixed with a strongly-coupled sector, specified by a (gapless) conformal field theory. For certain choices of conformal scaling dimension, $\Delta$, the exchange of unparticles leads to distinctive non-Gaussian features in the primordial curvature distribution, including bispectra with enhanced squeezed limits and oscillations close to the equilateral regime. Efficiently analyzing these models using Cosmic Microwave Background (CMB) data is a challenge since the shapes are non-factorizable in momenta and often highly degenerate with single-field self-interactions. Here, we overcome these limitations using a library of tools, including neural-network factorization schemes and optimal CMB estimators. Our pipeline condenses 161 non-separable templates into just 7 factorizable forms, with negligible loss of signal-to-noise. We apply the model to the \textit{Planck} data, asking two key questions: (1) can we detect unparticles? (2) can we distinguish them from single-field self-interactions? Across $1\leq \Delta\leq 9$, we find a maximal signal-to-noise of $1.2\sigma$, implying no evidence for new physics. We also place the first CMB constraints on the modified consistency-condition-satisfying orthogonal bispectrum with $f^{\rm orth^*}_{\rm NL} = -12\pm12$. While many unparticle models are degenerate with single-field shapes, values of $\Delta$ close to half-integers have very different shapes, offering an intriguing future discovery channel. The methods developed herein can be directly applied to other classes of templates, motivating the exploration of models beyond the standard weakly-coupled paradigm.

We present gas-phase radial metallicity profiles for 20 main-sequence galaxies at $4<z<6$, primarily based on JWST NIRSpec IFU observations obtained as part of the JWST-ALPINE-CRISTAL programme. Our study aims to connect the metallicity gradients of these galaxies with their kinematic properties from [CII]158$\mu$m ALMA observations. We map the radial profiles of oxygen abundance using the strong-line method leveraging the rich set of rest-frame optical emission lines. Linear fits to the annular-binned radial profiles show that, on average, the metallicity gradients are slightly positive with a median of $+0.039 \pm 0.010{\rm dexkpc^{-1}}$. There are no substantial systematic offsets in gradients when using different line diagnostics. However, only three galaxies show a gradient $>0.05{\rm dexkpc^{-1}}$ at $1\sigma$, and none have a significant negative gradient. We investigate the correlation between the metallicity gradients and the intrinsic gas velocity dispersion $\sigma_0$, as well as the ratio $V_{\rm rot}/\sigma_0$ of the disks. Combining our sample with mass-matched literature samples at $3<z<7$, we find a negative shallow correlation between $V_{\rm rot}/\sigma_0$ and the metallicity gradients, but no strong relationships with $\sigma_0$. As $V_{\rm rot}/\sigma_0$ increases towards later cosmic times, the observed negative trend with $V_{\rm rot}/\sigma_0$ is consistent with the overall cosmic evolution of metallicity gradients from high to low redshifts. This suggests that disk maturity plays a crucial role in shaping the radial metallicity gradients. [Abridged abstract]

In the past few years, ALMA unveiled a variety of substructures (rings, spirals, crescents) in the continuum emission of most protoplanetary disks imaged at high spatial resolution. While the majority of disks presents axisymmetric ring-like structures in the dust brightness distribution, some sources display asymmetric morphologies (blobs, crescents) that have been often associated to vortices and/or mechanisms generated by the presence of one or more embedded planets. In this brief research report we present the analysis of the arc structure observed in the dust continuum emission of the disk around HD~163296, using high resolution ($\sim$8~au) matched continuum data from ALMA at four wavelengths. We characterize in detail the arc structures and present a kinematic signature observed in the CS(3-2) emission at the same location. Our results indicate that the crescent is caused by differential dust trapping in a local pressure maxima, for which plausible mechanisms can be the presence of a vortex or trapping in a Lagrangian point of the planet-star system.

We investigate how reliably the global properties of Milky Way-mass dark matter haloes can be recovered from dynamical data over a limited radial range, particularly $\lesssim 30~\mathrm{kpc}$ where observations are most sensitive but baryonic processes modify the halo structure. Using the ARTEMIS simulations, which produce varying degrees of baryon-induced contraction, we fit dark matter profiles over restricted radial ranges using commonly adopted parametric models. Assuming negligible observational uncertainties allows the systematic errors from these choices to be isolated. When fits are confined to inner radii, an NFW profile underestimates the virial mass by a factor of $\approx 2$ on average ($\approx 4$ for some systems), and the concentration by a factor of $\approx 2$. Einasto and generalised-NFW models provide excellent local fits but retain similar global biases. In contrast, the contracted halo prescription from Cautun et al. (2020) yields stable extrapolations and recovers unbiased halo mass estimates over all radii. The inferred mass improves systematically with increasing radial coverage, and tracers beyond $\gtrsim 50~\mathrm{kpc}$ largely eliminate the mean bias for all models. The local dark matter density at the Solar radius is recovered to within $\lesssim 5\%$ for all profiles other than NFW. These biases are sufficient to reconcile recent low Milky Way mass estimates derived from inner rotation-curve analyses with the canonical $\approx 10^{12}~\mathrm{M}_\odot$. We additionally find a halo-to-halo scatter of $\gtrsim 0.1$ dex ($\approx 25\%$) persists even under idealised conditions, setting a likely lower limit for the precision of halo mass estimates.

The Southern Galactic Disk Survey (GDS) monitored a mosaic of 268 fields along a $6^\circ$-wide stripe in the southern Galactic disk with simultaneous observations in $r'$ and $i'$ ($7^\mathrm{m} \lesssim r', i' \lesssim 18^\mathrm{m}$) from September 2010 to September 2019. The survey design and data characteristics, as well as first results in $r'i'$, were presented by Haas et al. (2012; Paper I). Hackstein et al. (2015a; Paper II) extended the photometry and analysis process, and introduced the first catalogue including photometry of all 268 fields in $UBVr'i'z'$ and $r'i'$ light curves comprising up to 272 observations per field made between September 2010 and May 2015. Here we describe our custom-made observational scheduler and conclude the GDS with $r'i'$ light curves of up to 407 observations per field until September 2019 and $UBVz'$ light curves for a fraction of the fields. $113\,449$ distinct sources are identified as variables. Together with Paper II, we identified $77\,592$ variables that are not listed in either the International Variable Star Index (VSX) or the cross-match catalogue by Gavras et al. (2023). All emerging catalogues, comprising light curves, photometry, and reduced images, are made publicly available via the German Astrophysical Virtual Observatory (GAVO).

The modified Castor, Abbott, and Klein (m-CAK) theory predicts different wind regimes based on the line force parameter for changes in ionization ($\delta$) and the rotation parameter ($\Omega$). Stationary hydrodynamic studies have reported ''forbidden regions'' or gaps in this parameter space where no steady-state solution exists, suggesting physical instabilities. We investigate the stability of wind solutions within these gaps for B-supergiants to determine if they correspond to physical instabilities or numerical artifacts. We perform 1D time-dependent hydrodynamic simulations, systematically exploring the full $(\Omega, \delta)$ space for three B-supergiant models ($T_{\rm eff}=15-25$ kK), adopting a fixed density boundary condition. Our simulations reveal stable stationary solutions continuously across the entire parameter space, effectively filling the reported gaps. The transition from fast to slow regimes is smooth but structurally complex. Within the gap, the velocity profile develops a distinct ''kink'' or extended plateau in the supersonic flow, allowing the wind to reach a stable state. The mass-loss rate ($\dot{M}$) varies smoothly without artificial jumps. We find that the $\dot{M}$ gradient depends on the radiative driving strength ($k$): while $\dot{M}$ increases with $\delta$ for standard driving ($k \approx 0.32$), it decreases for the weak-driving regime ($k = 0.1$), consistent with stationary predictions. Moreover, in this regime, the final solution depends on the initial flow acceleration, confirming multiple hydrodynamic solutions. We conclude the m-CAK solution space is continuous; reported forbidden regions are artifacts of stationary methods. Time-dependent simulations effectively bridge the regimes, suggesting these transitions correspond to metastable states.

Paleo-detectors are a proposed experimental technique for direct detection (DD) of dark matter (DM) via the read-out of DM-induced nuclear recoil tracks in natural minerals. The large detector mass required for the sensitivity of conventional DD experiments to rare events is replaced by the exposure of paleo-detectors to DM-induced nuclear recoils over geological timescales. In this paper, we extend previous theoretical predictions for canonical spin-independent coherent and spin-dependent scattering (proportional to $A^2$ and the spin of the nucleus, respectively). We estimate the sensitivity of paleo-detectors to interactions between weakly interacting massive particle (WIMP) DM and nuclei within the framework of a Non-Relativistic Effective Field Theory (NREFT), considering isoscalar couplings to nucleons for both elastic and inelastic scattering. Taking into account cosmogenic, astrophysical and radiogenic backgrounds, we project the 90% confidence-level (CL) upper limits on the isoscalar NREFT coupling constants for both scattering types. We consider representative read-out scenarios and examine several target minerals. The projected sensitivities of paleo-detectors are compared with the 90% CL limits from the XENON100, LUX-ZEPLIN, and PandaX-II experiments, as well as with the 95% Bayesian credible region of the 2D marginalized posterior distribution from SuperCDMS. For DM masses from 1 GeV-10 GeV, paleo-detectors are projected to have sensitivity superior to that of conventional experiments for WIMP-nucleus interactions via all NREFT operators, largely independent of read-out scenario or target mineral. For DM masses from 10 GeV-5 TeV, we find that the sensitivity of paleo-detectors is projected to be comparable to or better than that of conventional experiments for WIMP-nucleus interactions via several NREFT operators, depending on the read-out scenario and target mineral.

Pulsar Timing Array (PTA) experiments have entered a new era with evidence for a nanoHertz gravitational wave background (GWB). This review describes the physics of detection, detailing the noise models and cross-correlation techniques required to isolate the Hellings-Downs curve. We discuss astrophysical implications, arguing that the perceived tension between current amplitudes and standard merger models is largely resolved by new insights into supermassive black hole binary populations. Beyond the stochastic background, we review the framework for multi-messenger continuous gravitational-wave searches, highlighting targeted search campaigns and rigorous detection protocols. We also examine the potential to probe New Physics, including cosmic strings and ultralight dark matter. Critical challenges are addressed, including small-scale leakage bias in anisotropy searches and the separation of deterministic signals from the GWB and pulsar noise. Finally, we outline the field's future, from rapid data combination strategies to the sensitivity gains expected from the Square Kilometre Array Observatory (SKAO) and DSA-2000.

We conduct an environmental analysis around 167 star-forming galaxies (SFGs) detected by the Star Formation Across Cosmic Time (SFACT) survey over the redshift range 0.129 $\leq$ z $\leq$ 0.500. We use three environmental estimators to characterize the local galactic environments around the SFACT SFGs, on the scales of 100 kpc to several Mpc. We categorize these environments based on the relative clustering strength with respect to a deep environment comparison redshift sample. The SFACT SFGs tend to be less clustered than the environment comparison sample (ECS), with no significant change in relative clustering strengths over our redshift range. We find that any trends with the star-formation rates (SFRs) of the SFACT galaxies and their environments are likely related to their absolute magnitudes, a proxy for mass. Mergers and interactions with other luminous galaxies do not appear to be the primary driver of the star-formation activity seen within the SFACT SFGs.

The Universe is now extensively populated by discy galaxies with coherent galaxy-wise stellar rotation. This disc prevalence has been deemed a late-time phenomenon because the penetrating cold gaseous streams in the early Universe ($z\gtrsim 2$) fuel the star formation in galaxies too intensively to allow for thin disc formation. However, recent images taken by the James Webb Space Telescope (JWST) unveiled a prominent population of low-mass galaxies at high redshifts with flattened shapes, widely interpreted as early significance of discs given the well-established connection between flattening and discy morphology seen in the local Universe. It is noticed, on the other hand, that these galaxies show far more flattened systems than can be accounted for by randomly oriented oblate discs, and the axial ratio distributions are better explained by elongated prolate ellipsoids, an extremely rare spindle-like configuration at low redshifts. The true morphological nature of these early low-mass galaxies is fundamental to understanding the structure evolution of their discy descendants we see today, including our Milky Way. In this work, we discriminate the oblate disc and prolate spindle scenario by a decisive experiment with stellar kinematics at its core. The result clearly supports the prolate spindle scenario, and evidences an early Universe widely inhabited by linear stellar systems contrasting the current era dominated by planar discy galaxies, which suggests a dimensional transition in galactic structure over cosmic time.

Globular clusters are complex stellar populations that provide unique opportunities to study stellar evolution -- as the second brightest cluster, 47 Tuc is a prime target. Asteroseismology can be used to measure precise masses of stars and has recently been applied to red giants in globular clusters, but so far not for 47 Tuc. Here, we present a search for solar-like oscillations in red giants of 47 Tuc using 8.3 days of high-cadence Hubble Space Telescope data. We detect oscillations in two out of the five giants falling in the field of view, 5 arcmin from the cluster center. One is on the horizontal branch (HB) while the other is on the red giant branch (RGB) at a similar brightness. From the seismic signal, we measure the stellar masses to be $0.78\pm0.13\,$M$_\odot$ (HB) and $0.94\pm0.15\,$M$_\odot$ (RGB), and hence an inferred integrated mass loss along the upper RGB of $0.16\pm0.20\,$M$_\odot$. A mass uncertainty of less than 0.05M$_\odot$ would be required to obtain a useful estimate of the mass loss, while an uncertainty below 0.01M$_\odot$ would be required to measure the mass difference between the cluster's multiple chemical populations. The former would be attainable with observations of about 100 times more stars to form ensemble-averaged values (assuming a similar length campaign), or alternatively a longer campaign observing fewer stars. Detecting mass differences between the chemical sub-populations, would require a 20-day campaign observing several hundreds of stars. Our clear detection of oscillations and the prospects presented here warrant dedicated high-cadence campaigns of 47 Tuc, which are possible with NASA's Roman mission and future missions like HAYDN.

We report Subaru/PFS spectroscopic follow-up of a radio-loud quasar at $z=1.715$ from the UNVEIL radio AGN catalog and with X-ray detections. The PFS spectrum displays a broad MgII emission line with an $\mathrm{FWHM}\gtrsim3400\ km/s$, accompanied by a narrow absorption feature. The spectrum reveals a characteristic $\Lambda$-shape over the rest-frame wavelength ranging $\sim1500-3500\ Å$. This underlying UV continuum is too curved to be reproduced by simply applying dust extinction to the spectrum of typical unobscured quasars. Alternatively, it is well described by a blackbody spectrum with a temperature of $T\approx10000\ K$. This result is in good agreement with its UV to MIR photometry that can be well modeled by three blackbody components representing the SMBH envelope ($\mathit{T}\approx9700\ K$), dust torus ($T\approx1500\ K$), and host galaxy dust ($T\approx80\ K$). The source is marginally detected in the GALEX NUV, revealing a potential V-shaped spectral energy distribution around $1400\ Å$, reminiscent of the spectral feature reported for recently discussed LRDs whose V-shapes occur around $3000-4000\ Å$. This wavelength shift is broadly consistent with the temperature contrast between our blackbody component, with $T\sim10^4\ K$, and the lower effective temperature of $T\sim5000\ K$ expected for an optically thick photosphere surrounding the SMBH in LRDs. These properties suggest that this source might be caught in a transient evolutionary phase in which the dense gas envelope characteristic of LRD has begun to fragment, allowing us to witness the emergence of a quasar from an LRD-like state.

X-ray observations of galaxy clusters allow us to estimate the gas fraction, and thus the baryon fraction, and its evolution over time. This offers an additional cosmological probe as well as a probe of the gas behaviour in massive halos at the end of structure formation. However, cosmological and astrophysical effects are degenerate, and both should be modeled in order to explain observations; otherwise, the chosen baryonic model can potentially bias the cosmological results. We propose to quantify this effect by adopting a model-independent framework. We utilize Type Ia Supernovae to reconstruct the cosmic expansion history and apply the iterative smoothing method to infer the mass and redshift evolution of the hydrostatic mass bias. Our results confirm previous findings and show that the bias should evolve with time to reproduce CMB cosmological constraints.

Ultra-compact X-ray binaries (UCXBs) are a subclass of low-mass X-ray binaries (LMXBs) characterized by ultra-short orbital periods, typically less than $60-80\,$min. They consist of a compact mass-accretor and a hydrogen-poor mass-donor, in which the mass-accretor could be a neuron star (NS) or even a black hole (BH). UCXBs play an important role in multiple areas of astrophysics. In particular, they are considered strong, continuous gravitational wave (GW) sources in the low-frequency band, making them key targets for future space-based GW observatories such as LISA, TianQin and Taiji. As the most compact binaries, the formation and evolution of UCXBs remain highly uncertain. In this article, we review four classic formation channels: the white dwarf donor channel, the He star donor channel, the evolved main-sequence donor channel, and the accretion-induced collapse channel. We also discuss recent progress in these channels, covering evolutionary scenarios, the initial parameter space for UCXB formation, and associated objects. A comparison between observed UCXBs and theoretical expectations is provided, along with a discussion on the observed BH-UCXB candidates. The origin of UCXBs can be constrained by the chemical composition of mass-donors and their locations in diagrams of mass-transfer rate and X-ray luminosity versus orbital period. We also examine the implications of UCXBs for several astrophysical fields, including GW astronomy, multi-messenger astronomy, binary evolution, and NS physics under extreme conditions. Further progress will depend on multi-wavelength observations, the discovery of more UCXB samples, and more detailed theoretical simulations.

We develop a hybrid framework to identify kilonovae (KNe), using single-epoch, medium-band spectral energy distributions from the 7-Dimensional Telescope (7DT). The framework integrates an unsupervised anomaly classifier (\texttt{Isolation Forest}) to flag unusual events with a supervised multi-class classifier (\texttt{XGBoost}) that characterizes eight common transient types. Trained on realistically simulated 7DT photometry accounting for per-filter sensitivity, the classifier achieves macro $F_{1}\sim0.80$ ($\sim0.82$) with 20 (40) filters across eight classes, Type~Ia/Ibc/II SNe, SLSNe, TDEs, AGN, stellar variables, and asteroids. Without direct training, the anomaly detector recovers $>$90\% of simulated and observed optically detectable KNe (AT~2017gfo) with a low contamination fraction, with a caveat of limitations of the training sample such as limited redshift range of SNe ($z < 0.15$), and a relatively small number of early non-KNe spectra. A SHAP-based feature analysis reveals that only $\sim$40--50\% of the most informative filters are sufficient to retain near-baseline performance, while red-end filters contribute little. Combining the top-ranked half of the 40 7DT filters with a single LSST band reproduces the full-model accuracy within 1--2\%, suggesting practical follow-up strategies. These results demonstrate that 7DT's medium-band system enables rapid, interpretable classifications and reliable anomaly alerts from single-epoch data -- promising for gravitational-wave follow-up, Rubin alert stream filtering, and serendipitous transient discovery in the 7DT survey.

We report the discovery and characterization of PSR J1810-0623, a fully recycled millisecond pulsar with a spin period of 4.55 ms, discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) and followed up with FAST and the Green Bank Telescope (GBT). A phase-connected timing solution spanning over 6.5 years reveals a 15.4-day binary orbit with extremely low eccentricity (about 1.5E-5). Assuming a neutron star mass of 1.4 Msun, the inferred companion median mass (about 0.64 Msun) is consistent with a carbon-oxygen white dwarf, indicating an evolutionary origin in an intermediate mass Xray binary. The system's properties closely resemble those of other massive white dwarf binaries thought to form via Case A Roche lobe overflow, suggesting a prolonged accretion phase during which the neutron star was efficiently recycled. Polarimetric analysis of FAST data yields a moderate degree of linear polarization and a rotation measure of 86.6 pm 0.6 rad/m^2, offering constraints on the Galactic magnetic field. The inferred characteristic age (about 32 Gyr) and low surface magnetic field (about 1E8 G) indicate a highly recycled pulsar. Proper-motion measurements imply a modest transverse velocity, consistent with those of recycled millisecond pulsars in the Galactic field. Although the proximity of the globular cluster Pal 7 raises the possibility of a dynamical origin, discrepancies in dispersion measure and proper motion argue against a physical association. PSR J1810-0623 adds to the rare class of long-orbital period MSP-COWD systems and provides a valuable laboratory for studying pulsar recycling, binary evolution, and Galactic structure.

We present the first three-dimensional reddening maps of the Large and Small Magellanic Clouds (LMC and SMC) constructed using fundamental-mode RR Lyrae stars from the OGLE-IV survey. By applying a period-amplitude-color relation and a period-luminosity-metallicity calibration in the OGLE photometric system, we derive intrinsic colors, color excess $E(V-I)$, and photometric distances for more than 20,000 RRab stars in the LMC and 3,000 in the SMC. Spatial variations in reddening are modeled using an adaptive quadtree scheme, where robust reddening-distance relations are fit within each partition and distances are iteratively updated to achieve self-consistency. The resulting maps reveal resolved dust structures across both galaxies, including steep reddening gradients in the central LMC and flatter profiles in the SMC. The construction of the three-dimensional reddening maps further reveals that high-extinction regions exhibit reddening behavior inconsistent with a uniform extinction law, implying localized variations in dust properties. The final maps comprise 205 partitions for the LMC and 67 partitions for the SMC, and are released together with a Python-based query tool and GeoJSON data products. These 3D maps provide a foundation for distance-dependent reddening corrections and for probing the structure and physical conditions of the Magellanic interstellar medium, and future high-precision and cadence RR Lyrae sample from Gaia DR4 will support higher-resolution mapping and deeper exploration of dust substructure.

To expand the homogeneous microlensing planetary sample of the Korea Microlensing Telescope Network (KMTNet), we investigate six planetary candidates identified by the AnomalyFinder search in the 2023 prime-field data, namely KMT-2023-BLG-1592, OGLE-2023-BLG-0766, KMT-2023-BLG-0332, KMT-2023-BLG-0486, KMT-2023-BLG-0792, and OGLE-2023-BLG-1043. Light-curve modeling indicates that the first two events have planetary mass ratios of $\log q \sim -3.0$ and $-2.6$, while the third exhibits a brown dwarf mass ratio of $\log q \sim -1.4$. The remaining three events show the well-known degeneracy between the binary-lens single-source (2L1S) and single-lens binary-source (1L2S) interpretations. A Bayesian analysis yields companion masses of about 0.6 and 1.2 Jupiter masses for the two planetary systems, likely orbiting beyond the snow lines of M- or K-dwarf hosts. A review of the KMTNet planetary sample shows that candidates discovered by AnomalyFinder are significantly more likely to exhibit the 2L1S/1L2S degeneracy, consistent with the tendency of AnomalyFinder to detect subtler planetary signals.

We present a ground-based time-series photometric study of stellar variability in four intermediate- to old-age open clusters NGC 2192, NGC 2266, NGC 2509, and IC 1369 based on high-cadence Cousins R-band observations obtained with the 0.6 m VASISTHA telescope at the IERCOO observatory. The monitoring campaign comprises more than 34 h of time-series data, providing sensitivity to short-period variability on timescales of 0.02-2 d. We identified between 190 and 290 probable members in each cluster using a Gaussian Mixture Model. Structural parameters were derived from radial density profiles fitted with King models. Fundamental parameters were further constrained using color-magnitude diagram analysis with PARSEC isochrones, yielding ages of 0.3-1.6 Gyr and distances of 2.5-3.9 kpc. From the time-series photometry, we identify four new variable stars and seven previously uncharacterized periodic variables, including {\delta} Scuti and {\gamma} Doradus pulsators, as well as rotational variables. The detected variables exhibit periods between 0.12-0.90 d, with R-band amplitudes ranging from 0.01 to 0.20 mag. Periods were determined using Lomb-Scargle analysis of calibrated light curves. For a subset of variables, spectral energy distribution fitting was performed to derive effective temperatures (4300-10 000 K), radii (1.3-46 R_{\odot}), and luminosities (2-100 L_{\odot}), enabling reliable placement on the Hertzsprung-Russell diagram. We present PHOEBE light-curve modelling of the W UMa-type eclipsing binary Gaia DR3 2164531610149292288 in IC 1369, deriving its physical parameters and providing the first detailed characterization beyond its previously reported variability. These results demonstrate that combining dense-cadence ground-based observations with Gaia astrometry provides a reliable approach for identifying and characterizing variable stars in OCs.

Context. Data-constrained models incorporate observed photospheric magnetic fields. However, due to the lack of magnetic field information in the rest of the solar atmosphere, models rely on extrapolations that, in most cases, neglect the Lorentz force. Nevertheless, this force is present in the lower atmosphere and may play a key role in destabilising the equilibrium configuration and triggering eruptions. Aims. This study seeks to understand and reproduce a solar eruption SOL2014-12-18T21:41 that occurred in active region NOAA 12241, preceded by an M6.9 flare, and to investigate the impact of relaxing the initial force-free assumption. Methods. The resistive and compressible magnetohydrodynamic simulation is initiated using a non-force-free magnetic field extrapolated from a photospheric vector magnetogram taken minutes before the flare. The simulation includes a stratified atmosphere and non-ideal effects such as thermal conduction and radiative cooling. Results. A flux rope forms and rises in the simulation, carrying away dense material from the lower solar atmosphere. Its formation results from the non-zero Lorentz force acting on the initial sheared arcade, without assuming pre-existing flux ropes or photospheric driving motions. The flux rope is then deflected toward regions of low magnetic pressure, escaping the domain at 350 km/s with approximately constant acceleration. Conclusions. A robust numerical framework for modelling flaring active regions was applied to the eruption of NOAA AR12241 as a case study, assuming a realistic non-force-free magnetic field near the flare onset. It exemplifies how an initial Lorentz force imbalance can successfully trigger a flux rope formation that later escapes the simulation domain. It also enables comparison with real observations through the addition of a stratified atmosphere spanning from the photosphere to the corona.

We report a novel feature of relic gravitational waves (GWs) in non-singular bounce cosmologies that is testable in light of GWs astronomy. In non-singular bounce cosmologies, the effective potential $M_p^2 a^{\prime \prime}/a$ that governs the evolution of primordial GWs contains two peaks due to the existence of contraction phase prior to the standard expansion phase. Accordingly, relic GWs interference between the two peaks, resulting in a distinctive oscillatory feature in the spectrum, analog to the resonant tunneling effect in quantum mechanics. As a result, the GWs spectrum exhibits an oscillatory patterns on high frequecy regime, distinctive to other cosmological scenarios such as inflation. We show that the amplitude of GWs spectrum is high enough to reach the sensitivity of current and forthcoming GWs instruments, making our predictions falsifiable. Hence, our finding offers a promising way to experimentally test the non-singular bounce scenarios and search for new physics in early universe cosmologies.

BL Lacertae is a blazar known for its high flux variability and occasional broadband flares of unknown origin. It was in an extended flaring state from July 2020 until the end of 2021, making it an ideal candidate to study spectral and temporal properties during different flux states. We analysed five XMM--Newton EPIC observations of BL Lacertae taken up to the end of 2021. Temporal properties were investigated using fractional variability, minimum variability timescale, and the discrete correlation function. Detailed spectral modeling was performed on the two most variable observations, including correlation analysis between the soft (0.3--2.0 keV) and hard (2.0--10.0 keV) bands. Two of five observations were found to be highly variable with $F_{\mathrm{var}} = 19.16 \pm 0.32$ and $6.27 \pm 0.43$. The 2021 observation corresponds to the highest flux state. The shortest variability timescale in the 0.3--10 keV band is 1.24 ks. Assuming synchrotron-dominated X-ray emission, this timescale constrains the emission region size. Under equipartition between the magnetic field and radiating particles, this implies $B \approx 0.4\,\mathrm{G}$. A softer-when-brighter spectral trend was found, as commonly seen in blazars. Spectra were modeled with single power-law, log-parabola, and broken power-law models; the broken power-law gave the best fit by Akaike Information Criterion in most cases, with a strong break energy--flux correlation. A thermal blackbody component showed a positive temperature--flux correlation in some observations. The spectral break, interpreted as the synchrotron cooling break, shifts to higher energies with increasing flux. The source consistently showed softer-when-brighter behavior. Only one observation showed significant soft--hard band correlation. The data suggest the synchrotron peak moves into or across the X-ray band as the source brightens.

HAT-P-32 b is the first exoplanet for which absorption in hydrogen and helium lines has been measured simultaneously. In addition to the relatively large maximum depth of ~5% in both lines, observations have revealed very long pre-transit signatures. In this paper, we apply a 3D aeronomy model to simulate the detailed spectrally resolved absorption profiles at mid-transit, formed in the upper atmosphere, as well as the integral transit curves formed by atmospheric material accumulating around the star. By fitting the observations, we derive key atmospheric characteristics, including the atmospheric composition and metallicity, and constrain the stellar XUV, hard X-ray, and Ly{\alpha} fluxes, as well as the intensity of the stellar plasma wind. We show that the presence of metals crucially affects the absorption in hydrogen and helium lines by accelerating the dissociation of H2. We also demonstrate that several processes not previously included in the modeling of HAT-P-32 b are important for the interpretation of the observations, such as the 3D asymmetric structure of the atmosphere, heating and cooling by excited states of hydrogen and metals, and the stellar VUV flux.

We present our shape and rotational characterization of four near-Earth asteroids observed as part of the Visible NEAs Observations Survey. This work includes 61 new light curves obtained between 2020 June and 2024 September for asteroids (21088) Chelyabinsk and (66146) 1998 TU3 and the potentially hazardous asteroids (4660) Nereus and (297418) 2000 SP43. We combine these observations with archival data to derive their shape models, spin parameters, and refined rotation periods using a light-curve inversion method. Our derived rotational periods are consistent with previously published values for all the objects. We present constant-period models for each asteroid, while for (66146) 1998 TU3 we detect a possible Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) acceleration of $\upsilon \simeq 2.05 \times 10^{-8} \text{ rad day}^{-2}$, potentially adding it to the short list of asteroids known to experience the YORP effect.

An important goal of the disk instability model is to explain the superhump phenomenon. Superhumps are features found in the light curves of binary systems, characterized by a period slightly different from the binary orbital period. In cases where the superhump period is longer than the orbital period (positive superhumps), they have been interpreted as arising from an eccentric, precessing disk. This paper reviews the theory and simulations that indicate that the disk's eccentricity originates from a dynamical instability at the 3:1 resonance. The instability is described by a mode-coupling process involving the interaction of the disk eccentricity with the binary tidal potential. This instability provides critical constraints on the nature of the disk turbulence that enables the disk to reach this resonance.

The R Coronae Borealis (RCB) Star, V854 Cen, was not discovered until the 1980's even though it is 7th magnitude at maximum light. This is because it was in a faint state due to the presence of thick circumstellar dust clouds for nearly a century. RCB stars are known for having deep declines of up to 9 magnitudes at irregular intervals. The declines are caused by dust clouds which block the light from the star. The historical lightcurve of V854 Cen before discovery was investigated by examining plates taken by Harvard College Observatory and other observatories beginning in 1889. These observations show that V854 Cen was well below 7th magnitude from 1890 to 1980, often as faint as 16th or 17th magnitude.

The origin of Saturn's rings has been debated for decades. Measurements from Voyager and Cassini have suggested that the rings could be as young as ~100 Myr and composed of nearly pure water ice. Several scenarios have been proposed to explain these properties. One hypothesis (Wisdom et al 2022) is that the rings formed through the recent tidal disruption of a pre-existing moon, Chrysalis, which experienced a close encounter with Saturn following its highly eccentric orbit. However, the mechanism by which this hypothesis would have formed the rings remains largely unexplored, in particular, whether Chrysalis could supply ring material of the desired mass and composition. To address these questions, we perform smoothed particle hydrodynamics simulations to investigate the tidal response of Chrysalis during close encounters with Saturn. Our results demonstrate that preferential tidal stripping of the ice mantle from a differentiated Chrysalis can produce rings with both mass and composition resembling the present rings -- provided that the closest encounter occurs between the parabolic Roche limits for ice ~1.53Rs and rock ~1.07Rs -- consistent with Wisdom et al 2022. Moreover, multiple close encounters can extend the effective disruption limit by spinning up the body, enhancing the tidal stripping efficiency. Following close encounters, the rocky remnant of Chrysalis would have been removed in less than few kyr, either by collision with Saturn or ejection onto a hyperbolic orbit. These findings support the hypothesis that Saturn's rings could originate from a recent lost moon, and imply a highly dynamical evolution of the Saturnian system over the past few hundred million years.

We investigate the rotation curve of the Milky Way using a multi-component mass model including a stellar disk, a gaseous disk, a bulge/bar component, and a dark-matter halo. The stellar and gas contributions are calibrated using recent observational determinations of the Galactic surface-density distribution, while the dark-matter halo is modelled with standard spherical profiles. We compute the circular-velocity contributions of the different components using a combination of spherical mass reconstruction for the bulge and halo, and thin-disk Hankel-transform methods for the disk and gas components. We first fit the stellar surface-density profile to determine a fiducial bulge-disk decomposition and then use this calibration to predict the Galactic rotation curve. We find that, although the resulting stellar mass model reproduces the observed surface-density profile reasonably well, it does not provide a fully satisfactory description of the rotation-curve data, with the largest discrepancies arising in the inner Galaxy. We then consider an alternative RC-first calibration strategy, in which the bulge and disk parameters are adjusted to improve the kinematic fit. While this significantly improves the agreement with the observed rotation curve, the corresponding stellar surface-density profile becomes inconsistent with the independently inferred baryonic distribution. Our results highlight a tension between photometric and kinematic constraints within simplified axisymmetric models and indicate that a fully consistent description of the Milky Way mass distribution likely requires a more realistic treatment of the bulge/bar structure and of baryonic systematic uncertainties.

We present evidence for atomic absorption lines in the high-resolution 4-30 A X-ray spectrum of the neutron star RX J0822-4300 in the supernova remnant Puppis A. Comparison with model atmosphere calculations shows that features in the observed spectrum can be uniquely associated with redshifted and pressure-broadened transitions in highly ionized oxygen and neon. We also spectroscopically confirm the previously estimated strength of the surface magnetic dipole field; we detect both the linear and the quadratic Zeeman effect. We derive values for both the gravitational redshift and the acceleration of gravity at the stellar surface, yielding the first purely spectroscopic estimates for the radius and mass of a neutron star.

Galactic foreground emission plays a key role in cosmic microwave background (CMB) science, particularly for detecting primordial gravitational waves. A well-known lesson is the ``dust wave'' identified by BICEP2 in 2014, which was ruled out through a more careful analysis of foreground emission. To date, most estimates of Galactic foreground emission have relied on the assumption that for each line of sight, only one component is considered per emission mechanism. However, the results in this work suggest that more complex modeling -- particularly involving multiple components arising from either line-of-sight complexity or pixel mixing -- may be necessary to fully account for Galactic foregrounds, including dust and other components. More interestingly, the only available two-component dust estimate also fails due to oversimplified emission parameters, although it is conceptually superior to single-component alternatives. These results yield three key conclusions: (1) Due to the intrinsic three-dimensional complexity of the Galactic environment, where physical conditions vary with both distance and direction, the actual radiation from Galactic foreground components cannot be accurately characterized by single-component models. (2) Consequently, CMB experiments require more frequency bands to resolve these components. (3) Spatial variations of foreground emission parameters should not be simplified, because in this work, all such simplifications are found to degrade the estimates significantly.

We extend the work of Roychowdhury (2026) on skewness variations of the logarithmic flux, driven by large GeV flares in FSRQs, to a sample of 18 FSRQs. We find that they can be categorized into three groups, one where the skewness attains a persistent lower value after a large flare, one where it does not, and those where change in skewness is not significant. To provide a theoretical ground for these results, we use the statistical plasmoid model of Fermo et al. (2010) that self-consistently produces large plasmoids through merging which, when gain energy from the reconnection event and are Doppler aligned, produce large flares. We find that a downsampling of our simulation of 1500 runs to 18 statistically reproduces the observed distribution in p-values for change in skewness. We further compute the ensemble Shannon entropy of the system and the skewness, where the entropy is found to decrease at a $3\sigma$ level in both the groups where skewness either increases or decreases, as a direct evidence of increase in order in the system caused by a flare. We find that the power spectral densities of the simulated light curves are broken-power-laws, resembling a white noise+red noise broken by the typical cooling timescale in our system, in accordance with known blazar variability. We find that our results are robust to a $200-300\%$ change in several fiducial parameters of the simulation. Our stochastic simulation of plasmoids inside a blazar jet self-consistently reproduces key observable statistical properties of blazar GeV light curves.

The stellar-to-halo mass (SMHM) relation of brightest cluster galaxies (BCGs) provides key insight into the connection between BCG growth and the assembly of their host halos. We analyze this relation using the spectroscopic SDSS GalWCat19 cluster catalog, selecting 996 systems with log(M200) >= 13.6, log(Mstar) >= 10.5, and 0.02 <= z <= 0.125 to limit evolutionary effects and ensure stellar-mass completeness. We fit lognormal scaling relations with a Markov Chain Monte Carlo (MCMC) framework that accounts for measurement uncertainties and intrinsic scatter. For the fiducial SMHM relation, <log Mstar | M200> = alpha + beta log(M200/Mpiv) with log(Mpiv) = 14.2, we find a shallow slope beta = 0.17 +/- 0.03, normalization alpha = 11.04 +/- 0.01, and intrinsic scatter sigma_int = 0.19 +/- 0.01 dex. Recasting the relation in normalized form reduces the scatter to 0.16 +/- 0.01 dex, while including the magnitude gap M14 further reduces it to 0.14 +/- 0.01 dex. Variations in richness, redshift, and mass thresholds produce systematic shifts that are small compared to the statistical uncertainties, indicating that our inferred relations are robust to plausible selection choices. The reduced scatter when including M14 supports a picture in which BCG stellar mass reflects both halo mass and halo assembly history.

Time-resolved optical photometry, complemented by TESS data and long-term survey light curves, reveals that the transient object CRTS SSS100505 J093417-174421 is an eclipsing novalike cataclysmic variable of the VY~Scl subtype, with an orbital period of 0.16329188(8) d. An analysis of the light curves with eclipse mapping techniques and an entropy landscape procedure indicates an orbital inclination of 81.5 degrees and a mass ratio of 0.45. Eclipse maps reveal two diametrically opposed asymmetric arcs of enhanced emission in the intermediate and outer regions of an accretion disk elongated in the direction perpendicular to the line joining both stars, interpreted as tidally-induced spiral shock arms. The accretion disk is 50 per cent larger in the longer wavelength TESS data than in the optical range, in line with the expected radial temperature gradient of an opaque steady-state disk. The combination of a small optical disk radius (of 21 per cent of the orbital separation) and high orbital inclination explains the relatively faint absolute magnitude of $M_g = 7.44$ for a novalike variable.

We first investigate an approximate structure of the top region (TR) of a jet, sandwiched by a front shock from which the surrounding matter (SM) inflows and a rear shock from which the jet matter (JM) inflows. Since pressure in the TR is higher than that in the laterally outer space, both JM and SM flowing in the TR are pressed out from the side of the TR. Supposing a steady flow of SM and JM there, we construct a simplified two dimensional model on a structure of the TR. With help of the model, we next infer what happens when precessing jets go through the surroundings in the SS433-W50 system presuming a supernova remnant (SNR) occupies W50. If we assume reasonable density distributions of the SNR and the interstellar matter in a 10 $\sim$ 100 pc distance range, the density of the surroundings is found to be much higher than that of the jet so that the jet is largely braked in the TR and that outflowing rate of the energy from the side of the TR becomes almost identical to the intrinsic energy flow rate through the jet. The outflowing energy could spread to the ambient space in a form of a bow shock but the situation of the shock propagation in the present case could be peculiar due to the presence of the precession. Particularly, all the mass and the energy outflowing from the inner side of the precession cone is considered to be concentrated around the axis of the precession cone. As the result, mass-compressed and energy-accumulated regions are expected to appear along the precession axis, which could be the origin of the lobes laterally extending from the main sphere of W50 observed in radio and X-rays.

We present a user-friendly, open-source Wolfram Language package that automates the construction of an effective homogeneous generalized Voigt rheology for a spherically symmetric, incompressible layered body with Maxwell solid layers. It provides a practical bridge between layered interior models and time-domain simulations of tidal evolution. The package combines three components: (i) a forward computation of the degree-2 tidal Love number based on the propagator-matrix formulation for incompressible stratified viscoelastic bodies; (ii) numerical identification of the secular relaxation poles and residues of the layered model; and (iii) inversion of the resulting response into the compliance of an equivalent homogeneous generalized Voigt body. The implementation is based on the equivalence established for multilayer Maxwell bodies and includes an optional dominant-mode selection procedure for obtaining reduced rheological models over a prescribed frequency range. The package returns the parameters of the equivalent homogeneous model, including elastic, gravitational, viscous, and Voigt-element contributions, in a format suitable for downstream numerical applications. As a case study, we apply the package to a five-layer lunar interior model and obtain its equivalent generalized Voigt representation, together with a reduced model that preserves the tidal response over the frequency interval relevant for orbital evolution while using fewer relaxation elements. This package makes the reduction from stratified viscoelastic interiors to effective homogeneous rheologies reproducible and accessible. It allows physical tidal dissipation models to be used in long-term orbital and spin-evolution studies without having to repeatedly solve the full layered boundary-value problem.

Although the James Webb Space Telescope (JWST) has received much attention for its ability to search deeper into the cosmos than ever before, it also enhances our capability to study objects closer to us in the Universe. We apply a methodology of subtracting intracluster light to the PLCK G165.7+67.0 (G165; $z$ = 0.35) cluster, revealing a population of unresolved point-like sources including globular clusters (GCs). By applying a fitting algorithm in color space used to select galaxy cluster members, we uncover over 900 globular cluster candidates from our point source sample. We also identify candidates by estimating the contribution of interlopers to the point source sample, yielding an estimate of 793$\pm$ 83 globular cluster candidates. We find the color-selected sources to be approximately correlated spatially with the intracluster light and lensing mass of the cluster. The observed luminosity function of the sources shows a turnover point fainter than the completeness limit, so we use fixed-parameter curve fitting models to predict a K-corrected turnover point between $-9.4 \leq M_{\rm F200W} \leq -10.7$ mag, although we predict the expected K-corrected turnover point should be closer to $-7.7 \leq M_{\rm F200W} \leq -8.4$ mag. We discuss the dynamical state of this disturbed galaxy cluster with a bimodal mass distribution using the spatial distribution of GC candidates and find that the radial profiles of our color-selected GC candidates are very consistent with the lensing-derived surface mass density at $>$50 kpc.

We utilise the kinetic Sunyaev-Zeldovich effect (kSZ) to measure the rotation of ionised gas within galaxy groups defined in the SDSS-DR7 galaxy sample, via their dipolar imprint on the cosmic microwave background (CMB). We estimate the direction of the projected angular momentum for each group by measuring the redshift dipole of satellite galaxies around their group centre. We find a clear redshift dipole in the stacked data for the SDSS groups. We then perform oriented stacking of the Planck CMB temperature map using the group centres and directions of angular momenta. We report a $2.3\sigma$ measurement of the coherent rotational kSZ effect (rkSZ) within the virial radii of SDSS groups with an average mass of $10^{14}h^{-1} \rm M_{\odot}$. We estimate the averaged rotational velocity of the sample to be $\sim 100-200 ~\rm km ~s^{-1}$, peaking at approximately half the virial radius. Our results are consistent within the errors with predictions based on the ELUCID constrained realisation simulation, with the predicted amplitude of the rkSZ signal being slightly lower near the centre. We also identify a systematic bias when estimating rotational velocities using the observed redshifts of galaxies, but find it to be subdominant for our analysis.

Photometric classification of Type Ia supernovae is essential for modern time-domain surveys, where spectroscopic confirmation is not always feasible for the full transient sample. In this work, we investigate a compact and physically interpretable feature representation derived from multi-band light curves and evaluate its performance using gradient-boosted decision trees on the Supernova Photometric Classification Challenge (SPCC) dataset. Starting from a reduced 16-feature model, we perform a systematic feature ablation study to determine which physical descriptors contribute most strongly to classification performance. The final compact model achieves an F1-score of approximately 0.844 and a precision--recall area under the curve (PR-AUC) of approximately 0.928. The ablation results show that temporal evolution provides the dominant classification signal, while brightness, color, and variability features supply complementary information. A reduced core of approximately ten physically meaningful features retains nearly the full performance of the compact model, indicating that reliable classification does not require large high-dimensional feature spaces. These results show that interpretable feature-based models can capture the essential astrophysical information needed for Type Ia photometric classification, with direct implications for survey cadence, filter coverage, and the design of transparent machine learning pipelines for future time-domain surveys.

We reconsider the late-time evolution of galaxy peculiar velocities in the 1+3 covariant approach to cosmological perturbation theory. It has recently been claimed that this approach predicts substantially stronger growth of peculiar velocities than standard metric-based perturbation theory -- on the grounds that the covariant treatment is fully relativistic whereas standard treatments are effectively Newtonian. We show that this is not the case. When the covariant equations are applied consistently, the $1+3$ approach reproduces exactly the standard perturbative result for peculiar-velocity growth. The stronger growth laws claimed in recent work arise from an inconsistent treatment of the coupled covariant system, in which terms constrained by the field equations are treated as if they were independent sources. Further claims are made that the stronger bulk flows can mimic accelerated expansion in a dust universe. We argue that these claims rest on a confusion between the kinematics of an arbitrarily chosen observer congruence and the physical expansion of the matter congruence traced by galaxies. We conclude that the standard treatment of peculiar velocities is correct and fully relativistic~-- and does not lead to anomalous bulk flows or to apparent accelerated expansion.

The phosphorus abundance distribution in field stars as a function of metallicity reveals a complex pattern. The LTE data for [P/Fe] in the low-metallicity range are sparse and scattered around [P/Fe]~ 0 dex. Near [Fe/H]~ -2 dex, the relative abundance [P/Fe] increases and reaches a maximum value of around [Fe/H]~ -1 dex. In this domain, P-rich stars and (super)phosphorus-rich stars are observed; the [P/Fe] value can exceed 1 dex. Until now, no attempts have been made to study the NLTE effects on the ultraviolet and infrared phosphorus lines in spectra of cool stars to test the robustness of the observed LTE phosphorus abundance distribution. We developed an atomic model of P I that can be used to analyze phosphorus lines in the spectra of cool dwarfs and giants in the NLTE approximation. The model was tested using the solar flux and intensity spectra, as well as the spectra of Procyon and sigma Boo. Profiles of 14 phosphorus lines in the infrared regions and equivalent widths were analyzed. Our NLTE phosphorus abundance in the Sun is (P/H)=5.35+/-0.04 dex. Using our NLTE model, we selected 12 ultraviolet and infrared phosphorus lines and calculated a grid of NLTE corrections for the following parameter ranges: Teff from 4000 to 6750 K, step 250 K; log g from 1 to 5 dex, step 1 dex; and Vt = 2 km/s, [Fe/H] from -3 to +0.5 dex, step 0.5 dex. The NLTE corrections were calculated for phosphorus abundance ratios of [P/Fe]=-0.4, 0.0, +0.4 dex. For the Sun, the NLTE correction is -0.08 dex. The grid of the NLTE corrections, as well as the direct line profile synthesis, were used to refine the literature data on the phosphorus abundance in metal-poor, intermediate-deficient, and solar-metallicity stars. NLTE corrections do not qualitatively alter the overall phosphorus abundance distribution over a wide metallicity range, and do not change the characteristic pattern of phosphorus-rich stars.

Accurate mass and radius measurements of small transiting exoplanets are essential for probing their compositions, formation histories, and potential habitability. We present a uniform analysis of six planetary systems (each hosting at least one small transiting planet): K2-79, K2-106, K2-111, K2-222, K2-263, and TOI-1634. Our study combines new CHEOPS transit observations with archival photometry from K2, TESS, and ground-based facilities, alongside new and archival radial velocity data from HARPS-N, HIRES, ESPRESSO, and others. For each system, we perform joint transit and RV modelling, achieving typical precisions better than 15% and 5% for mass and radius, respectively, and thus enabling precise bulk density determinations. These reveal a range of compositions, including rocky planets near the radius valley (e.g. K2-106 b, TOI-1634 b), intermediate-density planets requiring steam-rich or mixed volatile envelopes (e.g. K2-111 b, K2-263 b), and low-density regimes, consistent with gas dwarfs or water-worlds (e.g. K2-79 b, K2-222 b). Several systems show evidence of additional companions detectable via RVs but not seen in transit. The results highlight the value of coordinated CHEOPS and HARPS-N observations in delivering some of the most precise bulk densities for small planets to date and support the preparation for future atmospheric characterisation missions.

Context. For decades, reproducing the orbital period distribution of non-magnetic Cataclysmic Variables (CVs) seemed to require a drastic decrease, usually termed disruption, of angular momentum loss through magnetic braking at the fully convective boundary, which argued for a change in the dynamo mechanism operating in fully and partially convective stars. However, recent studies showed that the magnetic braking prescription traditionally used in CV evolution theory is clearly outdated as saturation, that is, a weak period dependence for rapidly rotating stars, is not included. Aims. Here we test an updated version of a saturated magnetic braking prescription that has been developed to explain the spin-down of single stars in the context of CV evolution. This prescription contains a boosting and a disruption parameter that represent the change in the strength of magnetic braking at the fully convective boundary. Methods. We performed state of the art MESA simulations for CVs with the revised saturated magnetic braking prescription. Results. As in previous studies, we found that magnetic braking needs to be stronger in close binaries than in single stars and that, in contrast to what is observed in single stars, magnetic braking needs to be reduced at the fully convective boundary. However, in contrast to previous studies of CV evolution, only a moderate disruption by a factor of 2 - 3 is sufficient to explain key features of the CV orbital period distribution and the measured mass-radius relation for CV donors. Conclusions. The relatively small decrease of the efficiency of magnetic braking at the fully convective boundary might have implications for our understanding of dynamo models for fully and partially convective stars.

The TOLIMAN project is engaged with the construction, launch and operation of a low-cost space telescope of unorthodox optical design. Its primary science goal targets an exhaustive search for temperate-orbit rocky planets around either star in the alpha Centauri AB binary within our nearest-neighbor star system. Despite their favorable proximity and brightness, the detection of terrestrial exoplanets around such nearby Sun-like stars remains problematic for contemporary instrumental approaches. By performing narrow-angle astrometric monitoring of binary stars at extreme precision, any exoplanets will betray their presence by way of gravitationally-induced perturbations on the binary orbit. Recovery of this signal is challenging for it amounts to only a few microarcseconds of angular deflection (at best), and so is normally thought to require a large (meter-class) instrument. By implementing an innovative optical and signal encoding architecture, the TOLIMAN space telescope aims to recover such signals with a telescope aperture of only 12.5cm. This paper gives an overview of key features of the mission; in particular the concepts underlying the optics to enable image registration at the extreme levels of precision required. An outline is also provided, sketching further mission components and systems incorporated into the 16U CubeSat spacecraft bus in which the science payload is housed - all of which are now under construction.

We present a unified \emph{background-level} framework that maps a prescribed late-time dark-energy density history $\rho_{\rm de}(z)$ onto an effective scalar-field description in a spatially flat FLRW universe. Working directly with $\rho_{\rm de}(z)$, we reconstruct the associated field trajectory $\phi(z)$, and field-space potential $V(\phi)$, together with a null energy condition (NEC) consistency check. We apply the method to three benchmark histories: (i) the Chevallier--Polarski--Linder (CPL) form; (ii) a smooth mirror AdS$\rightarrow$dS sign-switching profile in which $\rho_{\rm de}$ crosses zero at $z_\dagger$, interpolating between a positive late-time plateau and a negative high-$z$ plateau ($\Lambda_{\rm s}$CDM-like at the background level); and (iii) a shifted-$\tanh$ emergent profile that remains positive definite and approaches $\rho_{\rm de}\to 0^{+}$ at high redshift. Finally, treating the reconstructed potential, $V_{\rm tar}(\phi)$, as a target, we perform Bayesian model comparison directly in \emph{potential space} and rank representative analytic potential families by their Bayesian evidence. For CPL (restricting to the single-valued phantom branch for the potential-space comparison), the exponential potential has the highest evidence in the baseline analysis, while the shifted-$\tanh$ and hilltop quartic forms remain close competitors; for the sign-switching $\tanh$ target, the shifted-$\tanh$ potential is strongly preferred, and the emergent profile yields the same qualitative ranking. These results provide a practical dictionary between phenomenological expansion histories and the scalar-field potential shapes required to reproduce them at the background level.

We introduce an updated, time-dependent treatment to the split-tail ("croissant-like") heliosphere model with data-driven solar wind conditions at 1 au, to study the evolution of the heliosphere with solar-cycle variations in plasma speed, plasma density, and magnetic field intensity. The model produces a sub-Alfvénic and low beta region, not observed by the Voyagers, ~15 au ahead of the heliopause. The simulated magnetic field and radial flow depart from Voyager observations in this region, indicating that time-dependent effects alone are not sufficient to understand this regime of the heliosheath. We decompose fast and slow magnetosonic wave modes from time-dependent plasma pulse structures in the heliosheath, using a linear Riemann variable analysis, for the first time. Fast mode waves can both reflect at the heliopause and transmit into the interstellar medium, and their speeds are unaffected by the low beta plasma regime in front of the heliopause. The model reproduces the pf2 jump in the magnetic field at ~2020 in the interstellar medium and we attribute the source of pressure fronts observed by Voyager 1 in the interstellar medium, and pressure pulses observed by Voyager 2 in the heliosheath, to fast mode waves that are commonly recurring solar-cycle features. The presence of fast mode waves drive a highly variable termination shock, with average radial speeds of 6.05 au yr$^{-1}$ $\pm$ 5.37 au yr$^{-1}$ in the New Horizons direction. We find that the termination shock has a sinusoidal-like oscillatory motion in the rising phase of the solar cycle, and broad inward motions during the declining phase.

We present a toolkit, CosmoDS, designed to study cosmological models at the background level using dynamical system analysis within the Cobaya framework. Dynamical system analysis is a powerful mathematical approach for studying nonlinear systems and is widely used in cosmology to investigate the stability and evolution of different cosmological models, particularly those involving dark energy. In this code, we provide a framework for constraining cosmological models using a dynamical system formulation. Most importantly, the toolkit is directly integrated with the Cobaya interface, allowing users to take advantage of the sophisticated statistical and inference tools already implemented in Cobaya for cosmological parameter estimation and model analysis.

We report on the cosmic ray mass composition measured by the Telescope Array Low-energy Extension (TALE) hybrid detector. The TALE detector consists of a fluorescence detector (FD) station with 10 FD telescopes located at the Telescope Array (TA) Middle Drum FD Station (itself made up of 14 FD telescopes), and a surface detector (SD) array of scintillators. The array consists of 40 SDs with 400 m spacing and 40 SDs with 600 m spacing. In this paper, we present results on the measurement of the depth of shower maxima ($X_\mathrm{max}$) in the energy range from $10^{16.5}$ eV to $10^{18.5}$ eV collected over five years of the TALE hybrid detector. The $X_\mathrm{max}$ distributions were analyzed and compared with Monte Carlo simulations of proton, helium, nitrogen, and iron primaries, using the QGSJet II-04 hadronic interaction model. Our results indicate that the elongation rate of the mean $X_\mathrm{max}$, which is defined as the slope of $\langle X_\mathrm{max} \rangle$ versus cosmic ray energy, exhibits a break around $10^{17}$ eV. Up to this energy, the composition becomes increasingly heavy, characterized by a growing dominance of heavy nuclei and a steadily decreasing fraction of light primaries. Beyond this energy, the proton fraction increases significantly with energy. These findings suggest a transition from Galactic to extra-Galactic cosmic ray sources around the so-called second knee.

We present a spectral analysis of a black hole X-ray binary XTE J2012+381 during its 2022 outburst, using data from NICER and AstroSat. Combining data from NICER, LAXPC20, and SXT, we extract energy spectra covering the 0.7-10.0 keV range. We model the energy spectra using a series of physical models and find that a reflection-Comptonization model provides the best fit. Given the uncertainties in the black hole mass and source distance, we investigate the stability of the inferred spectral parameters by systematically varying the black hole mass (7.26, 11, and 16.5 M$_\odot$), source distance (3.3, 5.4, and 7.5 kpc), and spectral hardening factor (1.5, 1.7, and 1.9). We find that, across most combinations of these parameters, the spin solutions consistently lie in the high-spin regime, spanning values between $\sim$0.67 and $\sim$0.998, with only a limited subset of configurations favoring lower spins. In contrast, the disk inclination angle remains well constrained over the majority of the explored parameter space, typically ranging between $\sim$50° and $\sim$65°. Only a few parameter combinations yield higher inclination values.

Understanding the interplay between micrometeoroid impacts and solar wind irradiation is crucial for interpreting lunar surface evolution. Using reactive molecular dynamics simulations and surface binding energy (SBE) analyses, this study investigates the coupled effects of these two dominant space weathering processes on lunar regolith composed of Fe$_2$SiO$_4$. Our simulations reveal that micrometeoroid impacts significantly modify the lunar surface, creating structurally heterogeneous zones with varying SBEs across microcrater morphologies. Specifically, microcrater floors exhibit enhanced surface cohesion due to high-density compaction, whereas microcrater walls and ejecta show weakened structures. Applying Sigmund's sputtering theory with these SBEs indicates differential sputtering yields for Fe, Si, and O suggesting preferential retention of heavier elements like Fe. This selective sputtering mechanism supports the formation and growth of nanophase metallic iron (npFe$^0$) clusters, influencing the optical and compositional maturation of the lunar surface. These findings advance our understanding of lunar space weathering processes.

Context. Long-period Wolf-Rayet (WR) star binaries produced by mass transfer are predicted to be abundant, but are observationally rare. This yields constraints on the evolution of initially wide O star binaries, including those potentially leading to the formation of gravitational-wave sources through the Common Envelope Channel. Aims. We investigate this issue in the light of a new type of orbital evolution for initially wide O star binaries, which is driven by mass ejection at periastron passage during the Luminous Blue Variable (LBV) phase. Methods. The assumption that the mass ejection occurs instantly at periastron passage allows us to analytically describe the orbital evolution. This approach is motivated by our understanding of an Eddington-limit driven LBV phase. We perform population synthesis calculations for the WR stars in the Small Magellanic Cloud (SMC), and compare them to the observed SMC WR star population. Results. Different from mass transfer, our mass ejection scenario leads to increased orbital periods and eccentricities. The Galactic system WR 140 (orbital period 2895 d, eccentricity 0.9) could be a typical result of this evolution scenario. Our models predict measurable binary space velocities, and allow for the disruption of the binary. Our SMC population synthesis model predicts statistically 5.3 close, 3.7 long-period, and further 2 runaway single WR stars. With largely increased orbital periods and eccentricities, such WR+O star binaries may not be ruled out by past radial-velocity searches. Applying our scenario to the Gaia BH1 and BH2 systems, we find that it provides viable progenitor evolution models. Conclusions. The mass-ejection-driven orbital evolution could explain why so few wide WR binaries are observed, and why some of the apparently single WR stars have high space velocities. We discuss implications for gravitational-wave sources.

We present a comprehensive analysis of the 17-year Fermi-LAT observational data of Mrk 421 to investigate the spectral variations in the $\gamma$-ray bands. The light curve of the source in the 0.1--1000 GeV band with a 14-day time bin exhibits significant variability at a confidence level exceeding 5$\sigma$, which is accompanied by spectral variation, displaying a {\it harder-when-brighter} behavior. Moreover, its flux variation can reach up to one order of magnitude within one day, with a daily flux up to $(1.19\pm0.84)\times10^{-8}~{\rm erg~cm^{-2}~s^{-1}}$ on MJD 56152. The 17-year integrated spectrum of Mrk 421 necessitates a complex model for explanation, whereas its time-resolved spectra over one-day or several-day time intervals can be well fitted by a power-law model. We propose that the complex spectral shape of the 17-year integrated spectrum stems from the superposition of different spectral shapes in different flux states. By generating the GeV spectra that are simultaneously observed with the archived TeV observations and constructing the combined GeV--TeV spectra, we find that some combined GeV--TeV spectral shapes clearly imply different radiation origins for the GeV and TeV emissions, challenging the one-zone leptonic model. It is found that the flux follows a lognormal distribution, while the photon spectral index distributions can be well fitted by either a lognormal or a Gaussian functions. The possible nature of the $\gamma$-ray variability in Mrk 421 is discussed.

We investigate whether the first discovered fluorine-rich CEMP-no star, CS 29498$-$043, can be explained by a very metal-poor rotating massive star. We consider single rotating stellar models of 20 $M_{\odot}$ at a metallicity of $Z = 10^{-5}$, exploring initial rotation rates from $\upsilon_{\rm ini}/\upsilon_{\rm crit} = 0$ to $0.7$ in increments of $0.1$ ($0<\upsilon_{\rm ini}<644$ km s$^{-1}$). Rotational mixing enhances the production of light elements in the H--He layers, including fluorine. The ejected material can be nitrogen-rich without being fluorine-rich, whereas fluorine-rich ejecta are always predicted to be nitrogen-rich. The model providing the best fit to the abundances of CS 29498$-$043 is the $\upsilon_{\rm ini}/\upsilon_{\rm crit} = 0.6$ model ($\upsilon_{\rm ini} = 547$ km s$^{-1}$), which reproduces C, N, O, Na, Mg, and Al within the observational uncertainties. However, the predicted [F/Fe] $=2.8$ exceeds the observed value of [F/Fe] $=2.0 \pm 0.4$. By simultaneously varying the $^{15}$N($\alpha,\gamma$)$^{19}$F and $^{19}$F($\alpha,p$)$^{22}$Ne reaction rates within their acceptable ranges, the [F/Fe] ratio in the $\upsilon_{\rm ini}/\upsilon_{\rm crit} = 0.6$ model can be reduced to 2.2, providing a plausible solution to the abundance pattern of CS 29498$-$043. Our results support the hypothesis that fluorine-rich CEMP-no stars may originate from material enriched by a single, metal-poor, rotating massive star. A potential observational test of this scenario may be to check whether the nitrogen and fluorine abundances observed at the surface of CEMP-no stars are correlated.

We investigate the atomic-to-molecular gas transition in the molecular formation cloud HLCG 92-35. Using the ROHSA algorithm to decompose GALFA-H I data, we find the Lukewarm Neutral Medium (LNM) to be the dominant mass component, indicating a state driven out of thermal equilibrium by turbulence or past shocks. Spatial analysis reveals an inverse correlation between the phase distributions, with small-scale Cold Neutral Medium (CNM) structures embedded within an extended LNM this http URL Astrodendro, we identified 2,214 CNM clumps with sub-parsec scales. While the CNM mass spectrum steepens at high masses, its intermediate-mass slope matches that of CO clumps, suggesting that molecular clouds inherit the hierarchical structure of the CNM. Significant non-thermal linewidths and localized CNM-CO velocity offsets imply that the CNM consists of subsonic cloudlets moving collectively as aggregates. Our results show that these sub-parsec CNM structures are the fundamental building blocks of the cold interstellar medium, driven by thermal instability and turbulent compression.

The determination of the gravitational potential by the polyhedral method is revisited in the case where the surface of a body is composed of triangular facets. Based upon six test-shapes of astrophysical interest (sphere, spheroid, triaxial, lemon-shape, dumbell and torus) projected on nested grids, we verify that the convergence toward reference values is second-order in the step size of the grid, inside the body, at the surface and outside. We then show that the accuracy or computing time can be drastically enhanced by implementing the Repeated Richardson Extrapolation. This technique is especially efficient when the body's surface is smooth enough, and is therefore well adapted to the theory of figures (single and multi-layer fluids) and to dynamical studies (test-particle and mutual interactions), which require a large number of field evaluations. For real objects like asteroids that have very irregular terrains at small scales, the gain is modest. In that context, we estimate the discretization level beyond which the typical error in potential values due to altimetric uncertainties dominates over the contribution of sub-grid cavities and bumps. For bodies close to spherical, the criterion reads $T \gtrsim \frac{64 D}{3 \lambda},$ where $D$ is the diameter of the body, $\lambda$ the typical shape error and $T$ the number of triangular facets involved. The case of 433 Eros is considered as an example.

A hyperaccreting stellar-mass black hole (BH) has been proposed as the candidate central engine of gamma-ray bursts (GRBs). Comparing the predictions from the central engine models with the temporal behavior of GRBs is of great interest. In this paper, using the open-source GRMHD HARM-COOL code, we evolve several 2D magnetized hyperaccreting BH models with realistic equation of state in a fixed curved space-time background. We extend the code to include the calculation of neutrino annihilation power. We then study the time evolution of BH central engine parameters, i.e., the neutrino annihilation power, the Blandford-Znajke (BZ) power, and the initial magnetization $\sigma_0$. We find that the neutrino power is generally consistent with previous analytical results. Usually, the neutrino annihilation process tends to launch a thermal ``fireball'', while the BZ jet is Poynting-flux-dominated. Our results, especially the evolution characteristics of $\sigma_0$ may help to understand the complex GRB spectral behavior.

We investigate primordial non-Gaussianity in the Inflation without an Inflaton (IWI) framework, where scalar perturbations are generated at second order by primordial gravitational waves in Einstein gravity on an exact de Sitter (dS) background. Since scalar modes are produced nonlinearly from tensor modes, non-Gaussianity is an intrinsic prediction of the mechanism. We compute the corresponding scalar bispectrum, derive the relevant contribution to the three-point function of the scalar potential, and evaluate its shape numerically. We find that, unlike the scalar power spectrum, the bispectrum depends logarithmically on the ultraviolet cutoff set by the end of inflation, indicating a structural difference between the two- and three-point statistics in this scenario. Its shape is enhanced toward squeezed configurations, but its amplitude becomes strongly suppressed once the scalar power spectrum is normalized to the observed value. The resulting non-Gaussianity at CMB scales is therefore negligibly small, well below present observational sensitivity.

We present a study of the massive protocluster IRAS 15520$-$5234, which displays evidence of an explosive molecular outflow that unleashed a kinetic energy of at least 10$^{48}$ erg. The protocluster contains 16 dense cores detected in the ALMA band 6 continuum emission maps, having masses in the range from 0.2 to 11.0 M$_{\odot}$. Our analysis of CO $(2-1)$ emission reveals 28 well collimated outflow fingers, the majority of which follow a Hubble-Lemaître velocity law. The outflow fingers show no preferred orientation in the plane of sky and emerge from a common center of origin. We estimate the total mass, momentum, and kinetic energy of the outflow fingers and find that the values are at least one order of magnitude higher than the typical bipolar outflows associated with massive protostars. The morphology and kinematics of the outflow fingers suggest that the outflow associated with IRAS 15520$-$5234 is explosive in nature. We calculate the dynamical age of the explosive event to be approximately 6550 years. Additionally, we estimate the frequency of such explosive outflows in the Galaxy, which is one event every 83 years. Finally, we speculate that the rearrangement of masses within the massive protocluster and the dynamical interaction among the massive cores may result in the formation of such an energetic event.

Haumea is the third-largest of the five officially recognized dwarf planets and one of the four that reside in the trans-Neptunian region. It is among the most exotic bodies in the Solar System, with an exceptionally rapid rotation, a highly elongated triaxial shape, and a ring that orbits about three times more slowly than Haumea itself. Because of its large heliocentric distance, direct exploration by dedicated space missions is not feasible in the short term, so progress must rely on ground- and near-Earth facilities. Stellar occultations are among the most powerful tools to investigate trans-Neptunian objects. We present new predictions of stellar occultations by Haumea and its ring for stars down to Gaia G = 21, and assess their scientific potential, with special emphasis on the 4 May 2026 event. We computed occultation opportunities for the coming years and evaluated the 4 May 2026 geometry in detail, including Haumea's rotation phase, known 3D shape, pole orientation, and sky-plane motion, to estimate the expected shadow-path width. Because the target star has a very large Gaia RUWE, we also carried out a dedicated reliability analysis, including speckle observations. We identify eleven valuable events through 2030. For 4 May 2026, we derive an expected sky-plane shadow width of $2224 \pm 30$ km, substantially larger than conservative nominal assumptions and therefore highly favorable for observations. Speckle imaging reveals a companion at $\sim 0.12$ arcsec and $\Delta m \sim 3.1$; this companion is also expected to be occulted and shifts the nominal main-star path prediction on Earth by about 8 mas. These results confirm the strong scientific return expected from coordinated observations of upcoming Haumea occultations, especially the 4 May 2026 event, and provide an updated framework to improve constraints on Haumea's shape, density, ring properties, and environment.

The idea of the Chandrasekhar's conditions as equilibrium and stability of stars is revisited with a new universal three-parameter non-Maxwell distribution. We derive the maximum radiation pressures in the non-Maxwell distribution for a gas star and a centrally-condensed star, respectively, and thus we generalize the Chandrasekhar's conditions in a Maxwellian sense. By numerical analyses, we find that the non-Maxwellian distribution usually reduces the maximum radiation pressures in both a gas star and a central condensed star as compared with that cases if the gas is assumed to be a Maxwellian distribution.

Recent experiments in cosmology, particularly those aimed at detecting the faint, redshifted, global 21 cm hydrogen line (depth < ~200 mK, z > 7.5), have imposed stringent new requirements on radiometer calibration. In this work, we present a framework for circuit modeling and parameter inference to strengthen these calibration pipelines. This new approach enables in situ characterization of otherwise immeasurable systematics using physically motivated models. A combination of frequentist and Bayesian techniques are employed in a pipeline that supports iterative modeling, robust parameter estimation, and detailed uncertainty quantification. The framework is applied to the REACH telescope, where the precise correction of variations in the radio signal paths arising from component aging or environmental effects is critical. Circuit models of REACH's calibration sources are developed, with the goal of predicting source temperature corrections that are conventionally obtained from laboratory measurements. By fitting the models to measured data using a convolutional cost function, a strong agreement with RMS residuals no worse than -37 dB is obtained. However, Bayesian inference reveals that the resulting temperature corrections can have uncertainties on the order of 1 to 2 K, caused by reflection coefficient degeneracies, measurement noise, and errors in the models. To combat this, posteriors obtained from laboratory measurements are employed as updated priors, reducing correction uncertainties down to 75 mK. Ultimately, the framework provides a means of dynamically accounting for drift in system non-idealities over time, addressing the increasing precision demands of global 21 cm radio astronomy.

We present an investigation of the rest-frame optical/UV and X-ray properties for a sample of 3027 X-ray selected quasars between $1.5 \leq z \leq 3.5$ detected in the deepest Spectrum Roentgen Gamma/eROSITA data available and observed by the fifth iteration of the Sloan Digital Sky Survey (SDSS-V). We parametrize the CIV$\lambda1549$ emission line to infer the strength of accretion disc winds and perform X-ray spectral fitting. The X-ray spectral properties -- namely, the 2keV monochromatic luminosity (L$_\text{2keV}$) and spectral slope -- are not strongly correlated with wind strength. Despite this result, the X-ray selected sample is shifted towards lower CIV blueshifts and higher equivalent widths than the optically selected sample observed in previous SDSS surveys, and matching in optical luminosity, redshift, and Eddington ratio does not reduce these differences. We estimate the far-UV luminosity using the HeII$\lambda1640$ line luminosity and define the slopes between this and the 2500A monochromatic luminosity ($L_{2500}$) and L$_\text{2keV}$ ($\alpha_\text{ouv}$ and $\alpha_\text{uvx}$, respectively) in a similar manner to the familiar $\alpha_\text{ox}$ parameter, which tracks the spectral slope between $L_{2500}$ and L$_\text{2keV}$. The quantity $\alpha_\text{ouv}$ is more strongly correlated with wind strength in our sample than $\alpha_\text{ox}$. We show that the correlation between $\alpha_\text{ox}$ and wind strength is driven by the relationship between the optical luminosity and wind strength. Our results are consistent with a radiation line-driven wind, whereby the ionising far-UV photons must not over-ionise the gas. The hard X-ray photons are few enough in number to have a negligible effect on the ionisation state of the material.

We investigate the luminosity-temperature ($L$-$T$) relation of galaxy clusters as a probe for testing modified gravity (MG) theories, focusing on $f(R)$ gravity and symmetron models. Using an improved semi-analytic framework that incorporates angular momentum acquisition, dynamical friction, and shock heating within the modified punctuated equilibrium model, we compare predictions against hydrodynamical simulations and observational data. While massive clusters remain largely screened and follow standard $\Lambda$CDM predictions, low-mass systems ($kT \lesssim 1-2$ keV) exhibit systematic deviations characterized by steeper $L$-$T$ slopes in MG scenarios. Crucially, we demonstrate that these signatures cannot be mimicked by conventional astrophysical processes such as feedback or angular momentum effects, which primarily affect normalization rather than curvature. Our results establish the $L$-$T$ relation as a robust diagnostic tool for distinguishing general relativity from screened MG theories, with the strongest discriminatory power emerging at group scales accessible to current and future X-ray surveys. Moreover, a normalized reduced $\chi^2$ analysis of the $L$-$T$ relation shows that MG models provide significantly better agreement with observational data than $\Lambda$CDM, with several realizations achieving excellent fits while the $\Lambda$CDM model consistently performs worst.

When working with astronomical data, metadata is also important. A general-purpose file format for transmission, processing and archiving large datasets should facilitate, among other things, both efficient processing of bulk data and encoding of rich semantic metadata. When choosing a format for a particular purpose sometimes no existing format satisfies both these requirements adequately, but combining one data-efficient and one metadata-rich format can be made to do so. This paper discusses considerations for designing such hybrid data/metadata formats, and reviews some examples such as VOParquet, FITS-plus and ECSV. We focus on tabular data, but some of the considerations may apply to other datatypes such as arrays as well.

The anomalous rate of molecules ionization observed at the Central Molecular zone (CMZ) challenges known mechanisms of ionization observed in molecular clouds across the Galaxy, due to the exceptionally high levels of ionization measured (orders of magnitude above what cosmic rays can explain) and its uniform spatial distribution within the CMZ. Recent studies suggest that the source of the $511$~keV excess can be correlated with this anomalous ionization rate or contribute significantly to the ionization in the Galactic Centre (GC). One of the leading hypotheses attributes the $511$~keV signal to positron injection from radionuclides or pulsars distributed following the stellar bulge, which is rather flat around the GC and, hence, could help explaining the uniform ionization profile. In this work, we investigate whether such a population of sources, injecting MeV positrons at rates consistent with the $511$~keV observations, can account for the ionization levels and distribution observed in the CMZ. Our results indicate that positron injection alone falls short at explaining the anomaly, although their expected ionization is larger than expected from any previously studied candidates.

Likelihood-based forward modeling is standard in galaxy formation, but most implementations are formulated as forward maps rather than explicit trajectory-level likelihoods conditioned jointly on assembly history and environment. We introduce a Graph Path Likelihood Model (GPLM) on layered halo graphs, where temporal edges encode causal transport and coeval host edges encode environmental conditioning. On a fixed layered graph, the graph-conditioned path measure is written as $P(\mathbf{x}\mid G)\propto p_{\rm attach}(\mathbf{x}\mid G)\exp[-S(\mathbf{x}; G)]$, where $S$ is an effective action for dynamical increments and $p_{\rm attach}$ is a boundary measure for node entry. We also discuss a minimal preferential attachment-detachment prescription for the graph probability $P(G)$, which facilitates placing the likelihood within a cosmological ensemble of layered graphs. Trained on layered graphs reconstructed from TNG50-1, GPLM improves stellar- and gas-mass predictions over transport-only baselines. As fixed-graph applications, we evaluate dark-matter-deficient-galaxy operator averages, compute gas-channel response under controlled deformations, and compare full and host-ablated path measures through likelihood-ratio diagnostics. In these examples, higher-order satellites show a higher incidence of dark-matter deficiency and broader graph-to-graph variation, while the gas-rich response indicates more diverse environmental processing histories. GPLM thus provides a proof-of-principle likelihood framework in which trajectory likelihood ratios, operator averages, and response diagnostics become explicit statistical observables, with connections to astrophysical forward models, machine-learning emulators, and field-theoretic diagnostics.

We present GalRSG, a new catalogue of 227 Galactic Red Supergiants in the Scutum-Crux region, which provides a nearly coeval and co-distant sample, target of a long-term, panchromatic, high-cadence photometric monitoring campaign aimed at detecting pre-supernova variability and luminous eruptively driven mass-loss events.

Motivated by the value of CN-bearing molecules as tracers of interstellar physical conditions, we investigate the reactions of adsorbed CN radicals with acetylene and ethylene (C2H2 and C2H4) on interstellar dust-grain analogues using quantum chemical calculations. We find that reactivity is strongly controlled by the relative orientation of the reactants, with specific geometries either promoting or inhibiting reaction. We further show that, on ice, these reactions differ qualitatively from their gas-phase counterparts, stalling at the formation of the adduct complexes C2H2CN and C2H4CN and exhibiting newly emerged kinetic barriers for the neutral-radical association. We contextualize our calculations in the same reaction-diffusion framework that would be employed in astrochemical models, finding that, depending on the diffusion energy of the hydrocarbons, these reactions can be either negligible or efficient, highlighting the importance of the local ice structure in interstellar grain chemistry. These findings caution against the use of CN-based tracers that assume barrierless, bimolecular surface reactions involving CN radicals.

Although Gaia has identified the central stars of planetary nebulae (CSPNe) for about 70% of known Galactic planetary nebulae (PNe), reliable distance estimates remain incomplete, with fewer than one quarter having accurate parallaxes. Meanwhile, classical extinction-distance samples include only about 70 objects, corresponding to 1.8% of the Galactic PN population. We aim to construct a large and homogeneous catalogue of PN distances by refining extinction-distance measurements with Gaia DR3, providing an independent complement to CSPN parallax-based distances. We develop a Gaia-based extinction-distance method by combining an improved blue-edge approach with an extinction-jump model. PN distances are derived from stellar extinction jumps in line-of-sight extinction-distance profiles and are further constrained by comparisons with published distances, the spatial distribution of stars relative to the PN centre, and the PN radius-distance relation. We obtain distances for 1,066 PNe with a median relative uncertainty of 13%, with about 87% of the sample having uncertainties below 20%. The catalogue includes 765 PNe whose CSPN parallaxes have uncertainties greater than 20% and 128 PNe without CSPN parallaxes. This method complements CSPN parallax-based approaches and extends the traditional extinction-based method to higher Galactic latitudes. For PNe with discrepant literature distances, it helps identify the more reliable estimates and assess CSPN identifications. We find a likely misidentification of the reported CSPN for Fr2-36 and analyse 33 PNe with two CSPN candidates, suggesting improved identifications for 15 objects. This catalogue represents the largest homogeneous set of extinction-based PN distances to date and provides a robust benchmark for studies of Galactic structure, PN populations, and interstellar extinction.

Recent analyses of Fermi-LAT data have identified a nearly spherical, halo-like excess of gamma rays peaking at E_gamma ~ 20 GeV. If interpreted as dark matter annihilation, the excess directly probes the thermal freeze-out epoch and therefore any non-standard corrections to the early-Universe expansion rate. In this work we examine the implications of this tentative signal for CET Omega, an informational and modular extension of relativistic quantum field theory and cosmology. CET Omega predicts a universal state-dependent modification to the radiation energy density of the early Universe, characterized by a doubly logarithmic correction originating from renormalized modular fluctuations in the spectral triple of the theory. The correction is negligible during Big Bang nucleosynthesis and recombination but becomes relevant during thermal WIMP freeze-out. We derive the correction from the modular two-point function, justify the onset scale associated with the informational sector, and compute its quantitative impact on freeze-out through numerical solutions. We also analyze the evolution of the informational field Phi_Omega(x) and show that it freezes in before the freeze-out epoch and survives to the present time under gravitational advection. The resulting modification induces percent-level shifts in the relic abundance and sub-percent morphological corrections to the annihilation gamma-ray flux. We compare the scenario with Early Dark Energy, kination, and varying N_eff models, and show that the parameter range 10^{-4} < alpha_log < 10^{-2} remains consistent with Planck, BBN, and BAO constraints while predicting potentially observable deviations in the gamma-ray morphology accessible to next-generation MeV-GeV missions.

Active galactic nuclei (AGN) feedback plays a significant role in many aspects of galaxy formation and evolution and has become a key ingredient in cosmological simulations. However, the subgrid models of AGN feedback in cosmological simulations such as IllustrisTNG (hereafter TNG) often overlook recent progress in the small-scale modelling of black hole (BH) accretion and AGN physics. In this study, we improve on this by incorporating central aspects of the MACER model, a framework that treats AGN physics in greater detail, into the TNG feedback implementation. Specifically, we adopt MACER-prescriptions for feedback output for high and low accretion rates in a new model while the estimation of the accretion rate remains unchanged. We test this updated scenario both for idealized elliptical galaxies and for a cosmological box. Compared to the original TNG model, the MACER-based simulation leads to a higher star formation rate (SFR) and BH accretion rate in ellipticals, yielding a gas density profile in better agreement with observations. In the cosmological simulations, the time evolution of the SFR density, galaxy stellar mass function at $z=0$, and $M_{\star}-M_{\rm BH}$ relation at $M_{\star}>10^{10.5}\,{\rm M_{\odot}}$ are similar in both models. The MACER model better reproduces low-mass BHs in low-mass galaxies, and yields milder quenching in massive galaxies, although this is accompanied by the absence of a pronounced colour bimodality. Still, the similarity of the outcomes underlines the self-regulated nature of BH feedback: for different feedback energetics, the accretion rate tends to adjust such that a similar total AGN feedback energy is released.

Several studies have reported a nearly linear correlation between the molecular gas and star formation rate surface density, the so-called Kennicutt-Schmidt (KS) law. We aim to retrieve the KS relation for a sample of four star-forming galaxies located in nearby clusters, disturbed by the effects of the ram pressure stripping, as testing this law in galaxies subject to different environmental conditions can provide key information on the physics of star formation. To perform our analysis, we used ALMA band 6 and band 3 data coupled with MUSE data at spatial resolution of ~1 kpc. Moreover, we analyzed data of star-forming complexes detected through their H\alpha ionized gas emission. We also derived the star formation efficiencies of the star-forming regions nested in these big complexes using the star formation rates derived from spatially resolved HST images and various recipes for the corresponding cold gas phase. We find that ram-pressure-stripped galaxies show normal-to-low star formation efficiencies, depending on the position within the galaxy and on the local gas density: the inner dense regions in the disk show higher efficiencies with respect to the outer regions, including the gaseous tails. The global relation between the star formation rate density and the molecular gas surface density is superlinear, likely suggesting the shortening of the depletion times at high gas mass densities caused by the ram pressure. Within the star-forming complexes, the star formation efficiency is very similar to the one observed at 1 kpc scale in undisturbed star-forming disks. Interestingly, this result holds also for the star-forming complexes located in the stripped gas tails. The analysis of HST resolved clumps suggests that the molecular gas is not uniformly distributed within the star-forming complexes, but its density distribution follows a steeper profile.

Cosmic rays are charged energetic particles that permeate the interstellar medium. Their sizeable energy share and penetration power makes them essential players in the dynamical and chemical processes that rule Galactic evolution, such as the launching of outflows and the formation of star and planets. For these processes low-energy (MeV-GeV) CRs are particularly important, both because they are the most abundant and because they have the largest cross-section for ionization. The study of cosmic rays naturally connects with gamma-ray astronomy, as high-energy photons are the principal products of their interaction with the interstellar plasma. In this article, after reviewing our current understanding of Galactic cosmic rays as derived from direct measurements, we present the state of the art regarding Galactic cosmic rays covering their direct observables, their acceleration processes and models for their propagation in the Galactic Disk. We present then an excursus on the current state of gamma-ray observations, and propose new prospects for investigating the physical properties of Galactic cosmic rays, exploiting the observational capability of future MeV missions.

Very low-luminosity objects in nearby star-forming regions have been identified as promising proto-brown dwarf candidates. The study of their multiplicity can shed light on the dominant formation mechanism of these substellar objects. We aim at studying the multiplicity of the very low luminosity object IRAM 04191+1522. To do so, we have obtained 0.89mm ALMA observations with a very extended configuration, achieving an angular resolution of ~0.04 arcsec (6 au at 140 pc). We have complemented our data with new VLA observations, and ALMA archival data at 1.3mm. As a result, we resolve IRAM04191+1522 into a close binary candidate for the first time. The binary is detected in the ALMA continuum data with a projected separation of ~80 mas, or 11 au at a distance of 140 pc. The two sources are oriented in the East-West direction, with the eastern component being brighter and more extended than the western one, which is marginally resolved. The analysis of C18O(2-1) archival data reveals gaseous material in rotation around the binary, presumably from a circumbinary disk with ~27 au of radius centered on the faintest ALMA component. A fit of the position-velocity diagram allows us to estimate a total dynamical mass for the system of 50+-40 MJup. Therefore, we classify IRAM04191 as a tight proto-brown dwarf binary candidate. The VLA data reveals the detection of a single object closer to the western ALMA source, and with a spectral index consistent with a radio jet.

Binary black hole (BBH) mergers detected through Gravitational Waves (GWs) are a promising probe for the cosmic expansion. These sources are standard sirens for which we can directly measure the luminosity distance, but their redshift is degenerate with the determination of their source masses. In analogy to standard candles, the redshift of standard sirens can be obtained using a calibration based on the source mass spectrum, but without the need for a cosmological ladder. It has been recently shown that a mismodeling of the BBH mass spectrum is very likely to introduce a bias in the determination of the Hubble constant. To tackle this issue, we develop a BBH population model based on Normalizing Flows, trained on synthetic BBH catalogs generated from astrophysical prescriptions, including binaries formed through both isolated stellar evolution and dynamical environments. We validate this approach with a mock BBH dataset, demonstrating that the Normalizing Flow framework faithfully recovers the true distribution and eliminates systematic biases in the Hubble constant inference. By using this model on GWTC-4.0 data, we obtain $H_0 = 71.62^{+4.04}_{-4.00}\; km \; s^{-1} Mpc^{-1}$ at 68.3% credible interval. Assuming the astrophysical prescriptions present in B-POP, we also show that the determination of $H_0$ is degenerate with the fraction of binaries born in the dynamical and isolated formation channel, with a Planck cosmology favouring $\sim 35\%$ binaries formed in the dynamical environment while a SH0ES cosmology favouring a value of $\sim 25\%$.

The formation of cosmological large-scale structure is usually described in terms of the evolution of density fluctuations and their statistical measures, such as the power spectrum and correlation function. However, these statistics characterize the amplitude structure of density fluctuations and do not directly describe the spatial redistribution of matter that occurs during structure formation. In this work we formulate cosmological structure formation as a transport problem of mass distributions using the Wasserstein distance from optimal transport theory. The generative process from the initial linear density field to the observed galaxy catalog is treated as a hierarchical mapping from a continuous density field to a galaxy point process, and an approximate expression for the Wasserstein distance between them is derived under the small-fluctuation approximation. We show that this distance naturally decomposes into contributions associated with three physical processes: mass transport by gravitational evolution, galaxy formation bias, and shot noise arising from the discrete sampling of galaxies. The gravitational transport term is expressed as an integral of the matter power spectrum, while the galaxy formation contribution appears as a weighted integral of the galaxy correlation function. The sampling term corresponds to Poisson shot noise originating from the discreteness of the galaxy catalog. These results provide a unified framework for describing cosmological large-scale structure formation from the perspective of transport geometry and suggest that the Wasserstein distance may serve as a new statistical quantity linking continuous density fields with observed galaxy catalogs.

TOI-3884 b is a rare super-Neptune transiting a fully convective M dwarf that hosts a persistent giant polar spot. Because the planet occults this active region during every transit, the system offers a unique laboratory to directly probe the stellar surface and spot properties. We present seven James Webb Space Telescope (JWST) transits of TOI-3884 b observed with NIRISS and NIRSpec, spanning 0.5--5.3$\mu$m. While all visits show a recurring spot-crossing signature, each transit exhibits a distinct spot-crossing morphology, enabling us to infer a stellar rotation period of $P$=11.102$\pm$0.003d and tightly constrain the pole-on stellar orientation ($i_{*}$=40.8$\pm$0.3$^{\circ}$, $\lambda_{*}$=148.9$\pm$0.4$^{\circ}$) and spot properties ($R_{\rm{spot}}=0.576^{+0.006}_{-0.005}$R$_{*}$, $\phi_{\rm{spot}}$=84.69$\pm$0.12$^{\circ}$). We leverage this orbital configuration to measure the first empirical panchromatic spectrum of an M-dwarf starspot with JWST, establishing a direct observational benchmark for stellar atmosphere models in the fully convective regime. Comparison with 1D NewEra and SPHINX atmosphere models indicates that the spot is 185$\pm$2K cooler than the photosphere, consistent with previous ground-based measurements and expectations for mid-M-dwarf spot contrasts. While the models reproduce the observed contrasts at wavelengths longer than 1$\mu$m, they significantly underpredict the contrasts at shorter wavelengths. These results demonstrate that M-dwarf stellar atmosphere models alone may not fully capture the wavelength dependence of stellar contamination in transmission spectra and highlight the importance of empirical spot spectra for robust interpretation of planetary atmospheres, particularly in the optical.

A puzzling implication of current data interpreted in the $\Lambda$CDM cosmology is the preference for a negative sum of neutrino masses. Moving to $w_0w_a$CDM brings an appreciable fraction of the neutrino mass posterior back to positive values, while the constant equation-of-state dark energy case $w$CDM does not. We investigate a variety of one-parameter dark energy equations of state (DE EoS), each variation with particular physical properties, to understand whether a two-parameter DE EoS is required to bring the neutrino mass positive. The conclusion is that certain one-parameter DE EoS can suffice, implying that the data are pointing toward physical characteristics rather than a broad degeneracy. The required characteristics are identified as phantom dark energy at high redshift, crossing $w=-1$ at lower redshift.

The Einstein Telescope (ET) is a third-generation underground gravitational-wave observatory designed to extend the detection sensitivity down to a few Hertz. Newtonian noise is expected to limit the low-frequency sensitivity of ET, particularly in the 1.7-6 Hz band. Most existing estimates rely on analytical or semi-analytical models assuming homogeneous or layered media, neglecting geological heterogeneity and complex wave interactions. In this work, we present a numerical framework for Newtonian-noise estimation based on spectral-element simulations of a seismic wave field. As a proof of concept, we first benchmark the numerical results against analytical plane-wave predictions in a two-dimensional homogeneous medium with a single surface source, demonstrating excellent agreement for both bulk and cavern contributions. We then extend the model to an array of 30 stochastic surface sources to approximate stationary ambient seismic excitation. The P-wave fraction inferred from the simulated wave field is, in this simple homogeneous case, significantly lower than commonly assumed, indicating enhanced prospects for Newtonian-noise mitigation. The framework is readily applicable to three-dimensional simulations and to integration of detailed local seismic models and topography, offering strong potential for site-specific Newtonian-noise estimation.

The presence (and nature) of variations in the stellar initial mass function (IMF) at substantially sub-solar masses and metallicities ($m$$<$0.5M$_{\odot}$, [M/H]$\lesssim$$-$1) remains poorly constrained. Predictions from simulations vary widely, while observationally, resolved star studies of ultra-faint dwarf galaxies (UFDs) suffer from small sample sizes and background galaxy contamination due to low projected stellar densities. As an alternative metal-poor target, we measure the IMF from resolved stars towards a carefully selected field in the Small Magellanic Cloud (SMC), leveraging a plethora of independent constraints on the target field stellar population including distributions of distance, %extinction, age and metallicity. We resolve $>$15,000 stars down to 0.16M$_{\odot}$ within a single pointing of NIRCam onboard JWST, using an observing strategy that minimizes contamination from point-source-like background galaxies. We explore three different functional forms of the IMF, forward modeling observed color-magnitude diagrams (CMDs) and luminosity functions. We find a best-fit single power law IMF slope of $\alpha$=$-$1.61$^{+0.03}_{-0.03}$, consistent with UFDs probed down to similar limiting masses. Fitting a broken power law IMF, we find low- and high-mass slopes of $\alpha_{1}$=$-$1.44$^{+0.04}_{-0.04}$, $\alpha_{2}$=$-$2.17$^{+0.11}_{-0.11}$ respectively, consistent with solar neighborhood values. Assuming a lognormal IMF, we find a characteristic mass and lognormal width of $m_{c}$=0.12$^{+0.03}_{-0.03}$M$_{\odot}$, $\sigma$=0.61$^{+0.07}_{-0.06}$M$_{\odot}$, allowing for characteristic masses lower than local values as seen in some simulations as well as low-metallicity Galactic clusters. Lastly, we quantify the impact of assumptions required in our analysis and discuss potential future improvements.

As the third paper in a series regarding the estimation of luminosity functions (LFs) via kernel density estimation (KDE), we present a further generalization of our framework by extending its applicability to multiple flux-limited samples. While our previous works addressed single flux-limited datasets, many practical surveys collect data from several disjoint sky regions with varying observational depths and flux limits. We introduce a piecewise estimation framework that partitions the luminosity-redshift plane into disjoint regions according to the staggered flux limits of the sub-samples. Within each region, we integrate data from all surveys capable of detecting sources into a combined sample and apply the transformation-reflection KDE method using the corresponding local flux threshold as the truncation boundary. This strategy allows for the full utilization of all available sources while maintaining rigorous statistical consistency. The robustness of this approach is validated through Monte Carlo simulations. Furthermore, application to SDSS DR7 and 2SLAQ quasar data shows smooth transitions across flux limits and excellent agreement with parametric models. This approach is fully supported by our previously developed \texttt{kdeLF} Python package.

We study the cosmic evolution of radio-loud active galactic nuclei (AGNs) using a beaming-minimized sample of 4{,}555 steep-spectrum sources over $0<z\lesssim4$, compiled from the XXL survey, VLA-COSMOS, and other wide-field data sets. We model the rest-frame 1.4 GHz radio luminosity function (RLF) with a luminosity-and-density evolution (LADE; DE+LE) framework coupled to a flexible local LF family. Among the tested parameterizations, Model~C is statistically preferred and provides a globally consistent description of the binned RLFs while remaining compatible with local RLF measurements and Euclidean-normalized source counts. In the fiducial solution, the LE term rises toward cosmic noon ($z\sim2$--3) and then flattens or mildly declines, whereas the DE term decreases monotonically with redshift. This combined evolution naturally reproduces the observed luminosity-dependent turnover redshift $z_{\rm peak}(L)$ (often termed ``cosmic downsizing'') without imposing \emph{a priori} distinct evolutionary laws for low- and high-power sources. We further show that the same LADE functional family calibrated for star-forming galaxies also describes radio-loud AGNs when fitted independently, enabling a unified two-component (SFG+AGN) model consistent with both the local RLF and source-count statistics. Finally, converting the AGN RLF to a kinetic luminosity function yields a radio-mode black hole accretion rate density (BHAD) whose redshift dependence closely tracks the radio-based cosmic star formation rate density (after a conventional rescaling), with both histories peaking near $z\sim2$.

Context. Galactic globular clusters are a very important tool in explaining the characteristics of the Milky Way. Therefore it is essential to determine the kinematical and dynamical properties of the new star cluster candidates, especially at the low-latitude regions that suffer from heavy extinction and crowding. Aims. In this work, we report the first spectroscopic analysis for seven recently identified star cluster candidates: CWNU 4193, FSR 1700, Garro 02, Patchick 98, FSR 1767, Mercer 08, and BH 140. Our aim is to determine the kinematical properties, such as the mean cluster radial velocity, and the dynamical properties, such as the orbital parameters and the global dynamical mass, of these clusters in order to spectroscopically confirm the nature of these seven stellar systems. Methods. We collected the high-resolution infrared spectra of 33 candidate members of these clusters using the WINERED spectrograph mounted on the Magellan Clay 6.5 m telescope. Using the WINERED spectra, we measured the radial velocity of each individual star to confirm its membership in the clusters. From the confirmed members, we derived the mean cluster radial velocity of each cluster. In addition, for these clusters, we computed the orbital elements using the GravPot16 model and estimated their global dynamical masses based on the virial theorem. Results. As a result, we confirmed enough member stars (from three to seven stars per cluster) to reliably derive the mean cluster radial velocity and compute the orbital parameters of the clusters CWNU 4193, FSR 1700, Garro 02, and BH 140. For clusters CWNU 4193, FSR 1700, and BH 140, the number of confirmed members also allowed us to estimate their global dynamical masses. Therefore, we successfully derived key kinematical and dynamical properties for four of the most obscured star clusters in the Milky Way.

In 2024, due to some operational changes, the Hubble Space Telescope began exhibiting undeclared loss of lock events. This loss of lock can result in the smearing out of light from the target field during an exposure, which leads to data degradation, which in turn may require data to be retaken. In this work, we investigate this `roll-drift' effect in ACS/WFC images. We quantify the impact of roll-drift on measurable parameters in data by using simulations and existing data reduction techniques. We identify a threshold of one such measurable parameter beyond which data may likely be affected by roll-drift, so that users can quickly and easily assess whether their data needs further attention.

I present a spatially resolved analysis of the polycyclic aromatic hydrocarbon (PAH) population in the Pillars of Creation within the Eagle Nebula (M16) using the James Webb Space Telescope's (JWST) Mid Infrared Instrument (MIRI) and Near Infrared Camera (NIRCam) imaging. By using mid infrared PAH sensitive bands, I derive resolved maps of PAH size and ionization state across the pillars and connect these directly to variations in the radiation field and gas structure. I present the first spatial maps of PAH ionization and size in the Pillars of Creation. The analysis reveals clear internal gradients that show the PAH population is strongly shaped by local conditions within the cloud, such as the local radiation intensity and orientation of the nebular structure. The intensely radiated regions show a neutral and large PAH population, possibly due to electron recombination in these regions. I measure a mean PAH size of 198 carbon atoms with an error bar of 1.22 for M16 and use the resolved emission structure to obtain a first-order estimate of the electron density in the molecular cloud. These results provide direct evidence that PAH properties in M16 are governed by the interplay between radiation and density on sub-cloud scales, demonstrating the power of JWST imaging to probe dust processing in star-forming regions.

We present very long baseline interferometric (VLBI) radio imaging and Fermi/LAT GeV $\gamma$-ray observations of the 2019 eruption of the symbiotic recurrent nova V3890 this http URL VLBI imaging spans 8 -- 51 days after eruption, synchronous with the detected $\gamma$-rays. VLBI imaging shows the eruption starts out asymmetric on day 8 with an eastern component brighter than a western component. By day 32 the blast is rather circularly symmetric, and on day 49, the nova shell is brighter along the north--south axis. This morphological evolution is explained by interaction with circumstellar material (CSM) comprised of a spherical wind plus an over-density in the orbital plane. Comparing radio images to optical line widths gives an expansion parallax distance of 6.8 kpc. In the first 32 days or eruption, VLBI images capture $>$80 per cent of the integrated flux (as measured by the VLA), implying that synchrotron emission dominates. A second peak in the VLA light curve is explained by an image on day 48 that reveals the nova shell surrounded by a diffuse halo, powered by synchrotron emission from particles that have diffused upstream of the shock. The $\gamma$-rays appear around optical maximum and remain detectable for 23 days; marginally significant $\gamma$-rays reappear around day 60, concurrent with the second radio peak. Modelling indicates radio and $\gamma$-ray emission arise in distinct shock regions: $\gamma$-rays from dense CSM in the orbital plane, radio from the more spherical CSM component. X-ray observations constrain the spherical CSM density, which is higher than in other symbiotic recurrent novae. Assuming equipartition, we estimate the fraction of the post-shock pressure in magnetic fields, $\epsilon_B = 3 \times 10^{-4} - 2 \times 10^{-3}$.

As modern astronomy confronts unprecedented data volumes, automated pipelines and machine-learning techniques have become essential for processing and analysis. As these workflows grow more complex, astronomers also require input and inspection tools that can keep pace. To address challenges in navigating multidimensional datasets for quality control and scientific interpretation, we present the immersive Data Visualisation Interactive Explorer (iDaVIE), a virtual reality (VR) software suite developed in collaboration with the astronomy community. iDaVIE enables users to import and render large 3D data cubes within a VR environment, offering real-time tools for selection, cropping, catalogue overlays, and exporting results back into existing pipelines. Built on the Unity engine and SteamVR, the system uses custom plug-ins for efficient data parsing, downsampling, and statistical calculations. The software has already been integrated into workflows such as verifying HI data cubes from MeerKAT, ASKAP, and APERTIF, refining detection masks, and identifying new sources. Its intuitive interface aims to reduce the cognitive load associated with higher-dimensional data, allowing researchers to focus more directly on scientific goals. As an open-source, scalable, and adaptable platform, iDaVIE supports continued development and integration with other tools. Version 1.0 marks a significant milestone, with planned enhancements including subcube loading, advanced rendering modes, video-generation scripts, and collaborative capabilities. By pairing immersive visualisation with robust interaction tools, iDaVIE seeks to transform how researchers engage with complex datasets and enhance productivity in the era of big data.

The implementation of the GRAVITY+ Adaptive Optics (GPAO) system at VLTI enables unprecedented sensitivity and stability in optical interferometry. This allows high-precision characterization of directly imaged exoplanets at medium spectral resolution, providing a new pathway for studying planetary atmospheres. We aim to quantify and characterize the short- and long-term stability of GRAVITY+ through a consecutive seven-hour observation of the bright and stable star beta Pictoris, providing a benchmark for future exoplanet observations. We developed SIMTERFERE, a data-driven simulation tool that reproduces GRAVITY+ on-star observations using ancillary instrument and telemetry data. By comparing the simulations with the measured coherent fluxes, we traced the origins of systematic flux variations and assessed their impact on exoplanet contrast measurements. We find that the approximately 10% variations are dominated by throughput changes driven by variable fiber coupling, which depends on wavefront stability, atmospheric dispersion, and residual fiber offsets. These variations appear as smooth continuum changes across wavelength and can be effectively mitigated using second-order polynomial corrections. After removing these instrumental effects, the remaining approximately 1% variations are almost purely of telluric origin, which we can reliably correct down to the photon-noise limit (0.1% precision) using a contrast spectrum approach with linear airmass interpolation. The GRAVITY+ inferometric instrument is highly stable: low-order continuum and telluric variations can be corrected with high precision, making it uniquely capable of high-fidelity characterization of directly imaged exoplanets.

Stochastic gravitational-wave backgrounds (SGWBs) of primordial origin offer a powerful probe of early-Universe physics and possible dark-sector dynamics. While most searches focus on the GW power spectrum, additional information is encoded in higher-order correlators that characterize the statistical properties of the signal. In this work we study non-Gaussian features of a cosmological SGWB generated at second order by vector fluctuations, a class of sources well motivated in early-Universe scenarios. Within this framework we develop tools to characterize higher-order GW correlators and compute representative four-point functions that generate a connected contribution to the GW trispectrum. We show that the trispectrum amplitude scales as the square of the GW power spectrum and peaks in characteristic folded momentum configurations, reflecting the structure of the nonlinear source. We then explore the observational implications. First, we demonstrate that the connected trispectrum contributes to the variance of two-point overlap reduction functions, including the Hellings-Downs curve relevant for pulsar timing arrays. We then construct the optimal estimator to measure the connected trispectrum with ground-based interferometers. Our results highlight how non-Gaussian SGWB statistics provide a complementary observable to probe the origin of GW backgrounds and to distinguish cosmological from astrophysical sources.

Context. The prompt-emission time profiles of GRB 230307A and other long-duration compact object merger (COM) candidates exhibit a unique set of temporal properties, characterised by a deterministic evolution of waiting times and pulse widths. Aims. We searched the Fermi/GBM catalogue for other unidentified long COM candidates exhibiting temporal properties similar to those observed in GRB 230307A. Methods. We examined the temporal and spectral prompt-emission properties of GRBs featuring at least eight light-curve peaks. For candidates, all with unknown redshifts, that exhibited properties similar to GRB 230307A, we analysed their trajectories in the Ep,i-Eiso plane as a function of redshift. We then evaluated the joint likelihood of their compatibility with the Ep,i-Eiso relation satisfied by the bulk of long GRBs. Furthermore, we calculated their minimum variability timescales (MVTs) for comparison against known COM and collapsar populations. Results. We identified 9 COM candidates with unknown redshifts and demonstrated that there are at least two outliers of the Ep,i-Eiso relation with 3.1 sigma (Gaussian) confidence level. Furthermore, their MVTs are more consistent with those of COM than with collapsar GRBs. Conclusions. These results indicate that this specific set of temporal properties can serve as a diagnostic tool to distinguish long-duration COMs from the broader collapsar population. Furthermore, our findings suggest that the fraction of unidentified COMs among long GRBs may be larger than previously assumed.

Pulsar wind nebulae (PWNe) are bubbles of relativistic particles, powered by the rotational energy loss of the central pulsars. The Crab Nebula, powered by the Milky Way's most energetic pulsar, was discovered by the Large High Altitude Air Shower Observatory (LHAASO) as a PeV gamma-ray emitter, thereby establishing it as an extreme particle accelerator along with multiwavelength observations. Here we report LHAASO's detection of a point-like ultrahigh-energy (UHE, photon energy $E>100\,$TeV) gamma-ray source associated with the PWN powered by PSR~J1849-0001, a pulsar of spindown power 50 times lower than the Crab pulsar. The measured gamma-ray spectrum extends to PeV energies following a power-law distribution, with the PeV luminosity a few times higher than that of the Crab Nebula. Combined X-ray observations constrain the average magnetic field within the PWN to about $3\mu\,$G, and reveal an extreme particle acceleration efficiency approaching or even exceeding unity. The result challenges the particle acceleration theory in PWN and implies non-ideal magnetohydrodynamics (MHD) conditions within the accelerator, potentially involving magnetic reconnection upstream of the termination shock.

T Coronae Borealis (T CrB) is a nearby recurrent nova expected to erupt in the near future, offering a unique opportunity to study particle acceleration and high-energy emission from novae in real time. We investigate the production of gamma-rays and neutrinos following the T CrB outburst by combining three-dimensional hydrodynamical simulations with a detailed diffusive shock acceleration model. Our simulations account for the complex circumbinary medium, including the red giant wind, equatorial density enhancement, and accretion disk. We compute spatially resolved spectra of accelerated protons and electrons at the forward shock, accounting for downstream velocity gradients and variations in shock properties. Using a multi-zone approach, we synthesize hadronic gamma- ray emission from proton-proton interactions, leptonic gamma-rays from inverse-Compton scattering, and the associated neutrino emission. We present predicted gamma-ray spectra, light curves, and images from our numerical models of T CrB, and assess their detectability with current gamma-ray and neutrino observatories. We find that the early high-energy emission is dominated by the ejecta, with the accretion disk significantly boosting the gamma-ray flux and particle normalization during the first hours after the outburst. By incorporating velocity gradients in the post-shock flow, we demonstrate that maximum particle energies can reach the PeV scale in high-energy explosion scenarios. We show that while GeV gamma-rays are prominent messengers, neutrino detection is feasible primarily in models with high explosion energy and high ambient density.

A finite, scheme-independent Casimir contribution to the stress-energy tensor arises naturally for quantum fields in universes with non-trivial spatial topology. We compute this Casimir stress-energy tensor contribution for a conformally coupled scalar field and for a minimally coupled scalar field. We show that, for the conformally coupled case, the backreaction of this contribution to the Einstein equations during an expanding de Sitter phase drives anisotropic expansion even when the Universe begins in a locally homogeneous and isotropic state. We conclude that quantum imprints of the underlying non-trivial topology inevitably give rise to local departures from homogeneity and isotropy.

We show that the Lagrangian of a Nambu-Goto $p$-brane satisfies the identity $\mathcal{L}_{\rm [\it p \rm]}=T_{\rm [\it p \rm]}/(p+1)$, with $T_{\rm [\it p \rm]}$ denoting the trace of the corresponding energy-momentum tensor, independently of the properties of the gravitational field. While for $p=0$ this reduces to the standard $\mathcal{L}_{\rm [0]}=T_{\rm [0]}$ relation, which determines the on-shell Lagrangian of point particles and their fluids, more generally it depends explicitly on the $p$-brane dimensionality. We explore the implications of this Lagrangian identity for the dynamics of non-self-intersecting cosmic string loops in a homogeneous and isotropic universe within $f(R,\mathcal{L}_{\rm m})$ gravity, showing that, unlike in general relativity, the proper mass of a cosmic string loop may evolve over cosmological timescales regardless of its small size or tension. Finally, we extend the analysis to the more general case of closed $p$-branes in $(N+1)$-dimensional Friedmann-Lemaître-Robertson-Walker spacetimes.

Ultralight dark matter is a hypothetical class of particle with a number of interesting theoretical and experimental properties, many of which are best understood by direct analogy with or application of undergraduate physics. We present a series of exercises and discussions which may inspire the reader to bring contemporary research on ultralight dark matter into the undergraduate classroom.

Blazars are unique astrophysical environments whose high-energy $\gamma$-ray spectra are susceptible to modulations in the presence of ultralight axions. We search for these modulations, induced by axion-photon mixing, in Fermi-LAT spectral data of previously unexplored blazar targets, focusing in particular on blazars 1ES 1959+650 and B2 1811+31, whose flare states provide a clean testbed for axion activity. In both cases, we find no evidence for axions, and set exclusion regions on the axion-photon coupling for masses between $10^{-9}$ eV $\lesssim$ $m_a$ $\lesssim$ $10^{-8}$ eV, with sensitivities typically reaching $g_{a \gamma \gamma} \sim 10^{-11} - 10^{-10}$ GeV$^{-1}$ depending on the assumed blazar modeling choices. We examine the broad impact of modeling uncertainties, finding that the resulting constraints can vary substantially across plausible configurations. We discuss the implications of these systematic effects and their relevance for similar blazar-like searches in the future.

Neutron-capture reactions provide essential nuclear-physics input for modeling the synthesis of heavy elements in stars. The growing precision of stellar spectroscopy and isotopic measurements in presolar SiC grains now demands cross sections with improved accuracy over the full energy range, and access to unstable nuclei relevant to slow (s-) process branchings and the intermediate (i-) process. This article reviews recent progress in direct neutron-capture measurements, focusing on time-of-flight (TOF) experiments at CERN n_TOF and complementary activation techniques. Substantial advances have been achieved for stable s-only and bottleneck isotopes, significantly improving constraints on s-process models. In parallel, the combination of high instantaneous neutron fluxes and advanced detector systems has facilitated first-time measurements of several radioactive branching-point nuclei. Feasibility studies, however, reveal current limitations related to sample availability, background conditions, and restricted energy coverage. In this context, the complementarity between TOF and activation emerges as a central strategy. Future developments, including high-flux facilities and novel inverse-kinematics experiments in ion storage rings, are expected to extend the boundaries of neutron-capture measurements, overcoming current limitation

Large-scale land enclosure for speculative mega-development constitutes a non-equilibrium spatial process whose velocity, topology, and irreversibility remain poorly quantified. We study the Pantai Indah Kapuk 2 (PIK2) coastal mega-development north of Jakarta, Indonesia, using eight years (2017--2024) of Sentinel-2 land-use/land-cover (LULC) data at 10-meter resolution. The landscape is projected onto a Marxian probability simplex partitioning terrestrial pixels into Commons, Agrarian, and Capital fractions. Fisher-Rao (FR) geodesic distances on this simplex identify a transformation pulse of $0.405$~rad/yr during 2019--2020, coinciding with major construction activity. Absorbing Markov chain analysis yields expected absorption times into the built environment of $46.0$~years for cropland and $38.1$~years for tree cover, with a pooled built-area self-retention rate of $96.4\%$. Percolation analysis reveals that a giant connected component containing $89$--$95\%$ of all built pixels persists at occupation probabilities $p \in [0.096, 0.162]$, far below the random percolation threshold $p_c \approx 0.593$, indicating planned rather than stochastic spatial growth. The box-counting fractal dimension of the urban boundary increases from $d_f = 1.316$ to $1.397$, consistent with increasingly irregular frontier expansion. These results suggest that information-geometric and statistical-mechanical tools can characterize the kinematic and topological signatures of capitalist spatial accumulation with quantitative precision.

We analyse screening in multi-field scalar-tensor theories, focusing on systems with a dilaton coupled to matter and an axion with a dilaton-dependent kinetic term, in the presence of both planetary and stellar density profiles. Using analytic arguments and fully coupled numerical solutions, we identify a regime in which full screening for a dark-energy-light, effectively unpinned string-dilaton, can occur without fine-tuning. The backreaction of the dilaton's partnered axion field can suppress the exterior scalar charge by selecting a minimum-energy configuration (the BBQ mechanism), yielding robust screening for generic axion gradients. In this regime screening is achieved by cancelling the dilaton's gradient rather than localising it. This reduces the exterior scalar charge and allows for gravity tests in the solar system to be passed. We then show that the more familiar thin-shell intuition need not apply in the multi-field setting. Axion surface gradients can drastically reshape the dilaton profile and drive a more localised transition without generically suppressing the fifth force. The exterior charge can remain essentially unchanged or even be enhanced as the shell is made thinner by a kinetically coupled field. Multi-field two-derivative dynamics therefore decouple localisation in thin shells from screening, evade single-field no-go arguments, and reopen viable parameter space for cosmologically light dilaton-like scalars with strong couplings to matter.

We obtain exact charged AdS black hole solutions in Einstein Lambda gravity including the effects of Born Infeld nonlinear electrodynamics and Perfect Fluid Dark Matter. The influence of the PFDM and BI parameters on the event horizon is analyzed. We compute the conserved and thermodynamic quantities and verify that they satisfy the first law of thermodynamics. Thermal stability is studied in the canonical ensemble using the heat capacity and Helmholtz free energy showing how PFDM and BI parameters affect local and global stability regions. We further investigate the thermodynamics in the extended phase space by treating the cosmological constant as thermodynamic pressure obtaining consistent conserved quantities and confirming the first law. The Ehrenfest equations are analytically verified demonstrating that the critical behavior corresponds to a second order phase transition. Heat engines associated with these black holes are also constructed to examine how PFDM and BI parameters influence their efficiency. Finally we analyze the geodesic structure through timelike and null trajectories using the effective potential determining conditions for stable and unstable circular orbits the innermost stable circular orbit and the photon sphere. PFDM significantly modifies the orbital structure while BI corrections are weaker.

Galaxy clusters are powerful probes of astrophysics and cosmology through gravitational lensing: the clusters' mass, dominated by 85% dark matter, distorts background light. Yet, mass reconstruction lacks the scalability and large-scale benchmarks to process the hundreds of thousands of clusters expected from forthcoming wide-field surveys. We introduce a fully automated method to reconstruct cluster surface mass density from photometry and gravitational lensing observables. Central to our approach is DarkClusters-15k, our new dataset of 15,000 simulated clusters with paired mass and photometry maps, the largest benchmark to date, spanning multiple redshifts and simulation frameworks. We train a plug-and-play diffusion prior on DarkClusters-15k that learns the statistical relationship between mass and light, and draw posterior samples constrained by weak- and strong-lensing observables; this yields principled reconstructions driven by explicit physics, alongside well-calibrated uncertainties. Our approach requires no expert tuning, runs in minutes rather than hours, achieves higher accuracy, and matches expertly-tuned reconstructions of the MACS 1206 cluster. We release our method and DarkClusters-15k to support development and benchmarking for upcoming wide-field cosmological surveys.

We investigate the properties of neutron stars when a massive scalar field, which could comprise all dark matter, is non-minimally coupled to the Ricci scalar. This coupling generates additional contributions to the field's effective mass, leading to tachyonic instabilities inside neutron stars and giving rise to rich phenomenology. Within this framework, we obtain neutron-star configurations with maximum masses exceeding 2 $M_\odot$, even when hyperons, which typically soften the equation of state and significantly lower the maximum mass, are included. Furthermore, we find that larger coupling strengths lead to multiple solutions for the scalar-field configuration. We analyze the structure of the corresponding effective potential responsible for this behavior. We also investigate how the inclusion of a scalar self-interaction term, in addition to the non-minimal coupling, modifies the resulting neutron-star properties.

The origin of cosmic angular momentum is a fundamental question in structure formation. We propose a novel mechanism that generates spatial angular momentum directly from quantum fluctuations during inflation. A spectator complex scalar field with global U(1) symmetry stores internal angular momentum via field-space rotation. Inflationary perturbations create spatial gradients that, upon horizon re-entry, couple to the background charge density and source a bulk momentum flow. During nonspherical gravitational collapse, this flow converts into net angular momentum. For primordial black holes forming from such collapse, the dimensionless spin can reach \(\chi\sim 0.1-1\) when the small-scale power spectrum is enhanced to produce detectable abundances-far exceeding tidal torque theory predictions. This establishes a testable link between inflation, primordial perturbations, and black hole spin distributions accessible to gravitational-wave observations.

We report a laboratory measurement of how impurities shift Coulomb crystallization in a strongly interacting ionic system. This is achieved by using laser cooled Ca$^+$ crystals doped with a controlled number of Xe$^{12+}$ highly charged ions. We find that the crystallization threshold is unchanged at low impurity concentration, but shows a clear crossover once the impurity content becomes sufficiently large, after which the shift grows approximately linearly. Complementary measurements reveal that this global effect originates from a local pinning of the crystal around the impurities. We further show how the measured shift could impact standard models of crystallization in white dwarfs and neutron stars. Our results provide an experimental route to incorporating impurity effects into models of multicomponent Coulomb matter, relevant to stellar crystallization and strongly coupled plasmas.

We investigate how a spatially modulated real scalar background $\phi(\vec{x})$ can modify phonon propagation in the context of Superfluid dark matter (SFDM). Using a simple toy model with quartic condensate and coupling $-g\phi^2|\Psi|^2$, we derive the local equation of state and the effective sound velocity $c_s(\vec{x})$. For $g>0$, modulation tends to increase the effective mass of the condensate and make the medium less rigid, suppressing $c_s^2\propto m_{\Psi,\mathrm{eff}}^{-4}$ up to a ``dust-like'' regime, $c_s^2\to 0$. We implement this modulation for the background scalar field by imposing rational profiles, through Padé radial profiles, and show the corresponding variation of $c_s^2(r)$ for different $g$, discussing implications for the structure of SFDM cores and the possible formation of inhomogeneous regions of dark matter.

The Einstein Telescope (ET) is a third-generation underground gravitational wave observatory designed to achieve an unprecedented sensitivity down to 3 Hz. Waves propagating in the soil due to anthropogenic or natural vibration sources generate density fluctuations which cause gravitational attraction, resulting in motion of the mirrors of the laser interferometer known as Newtonian noise. The latter is computed by integrating density fluctuations due to seismic wave fields over the soil domain surrounding the test mass. ANNA Newtonian Noise Analysis is a toolbox that computes Newtonian Noise from a seismic wave field defined on a finite element mesh, using Gaussian quadrature. 3D finite element meshes composed of linear and quadratic tetrahedral (4-node and 10-node) and brick (8-node and 20-node) elements are supported. The user computes (or interpolates) a seismic wave field on a finite element mesh and the toolbox computes the total Newtonian noise, as well as the bulk and surface contributions. ANNA runs in the MATLAB programming and numeric computing platform and is compatible with the open-source GNU Octave Scientific Programming Language; a Python version is also available. The toolbox is verified for plane P- and S-waves propagating in an elastic homogeneous full space with a mirror suspended in a spherical cavity, assuming that the wavelength is much larger than the radius of the cavity, so that wave scattering can be ignored. Excellent agreement with analytical solutions is obtained. Similar good agreement is reported for the Newtonian noise on a test mass suspended at a finite distance above the free surface of a homogeneous elastic halfspace in which a Rayleigh wave propagates. The proposed finite element framework provides a physically consistent and computationally efficient approach for computing gravitational-seismic coupling in heterogeneous media.

The nature of dark energy remains a central problem in cosmology. A compelling possibility is that dark matter is macroscopic, consisting of composite objects formed in the early Universe. We introduce the QCD-AQN framework, a well-motivated scenario in which dark matter is composed of dense aggregates of quarks and antiquarks matter stabilised by axion domain walls. The framework proposes a unified explanation for both dark matter and the observed matter-antimatter asymmetry. Particular emphasis is placed on existing observational constraints and on observational tests. Finally, we explore a possible QCD-based scenario for dark energy.