This study investigates the impact of spectral filtering on color-matching functions (CMFs) and its implications for observer variability modeling. We conducted color matching experiments with a single observer, both with and without a spectral filter in front of a bipartite field. Using a novel computational approach, we estimated the filter transmittance and transformation matrix necessary to convert unfiltered CMFs to filtered CMFs. Statistical analysis revealed good agreement between estimated and measured filter characteristics, particularly in central wavelength regions. Applying this methodology to compare between Stiles and Burch 1955 (SB1955) mean observer CMFs and our previously published "ICVIO" mean observer CMFs, we identified a "yellow" (short-wavelength suppressing) filter that effectively transforms between these datasets. This finding aligns with our hypothesis that observed differences between the CMF sets are attributable to age-related lens yellowing (average observer age: 49 years in ICVIO versus 30 years in SB1955). Our approach enables efficient representation of observer variability through a single filter rather than three separate functions, offering potentially reduced experimental overhead while maintaining accuracy in characterizing individual color vision differences.

Transmission spectroscopy is a key avenue for the near-term study of small-planet atmospheres and the most promising method when it comes to searching for atmospheres on temperate rocky worlds, which are often too cold for planetary emission to be detectable. At the same time, the small planets that are most amenable for such atmospheric probes orbit small and cool M dwarf stars. As the field becomes increasingly ambitious in the search for signs of even thin atmospheres on small exoplanets, the transit light source effect (TLSE), caused by unocculted stellar surface heterogeneities, is becoming a limiting factor: it is imperative to develop robust inference methods to disentangle planetary and stellar contributions to the observed spectra. Here, I present STCTM, the STellar ConTamination Modeling framework, a flexible Bayesian retrieval framework to model the impact of the TLSE on any exoplanet transmission spectrum, and infer the range of stellar surface parameters that are compatible with the observations in the absence of any planetary contribution. With the "exotune" sub-module, users can also perform retrievals directly on out-of-transit stellar spectra in order to place data-driven priors on the extent to which the TLSE can impact any planet's transmission spectrum. The input data formats, stellar models, and fitted parameters are easily tunable using human-readable files and the code is fully parallelized to enable fast inferences. [shortened for arxiv; see full summary in the PDF]

With the advancement of third-generation gravitational wave detectors, the identification of strongly lensed gravitational wave (GW) events is expected to play an increasingly vital role in cosmology and fundamental physics. However, traditional Bayesian inference methods suffer from combinatorial computational overhead as the number of events grows, making real-time analysis infeasible. To address this, we propose a deep learning model named Squeeze-and-Excitation Multilayer Perceptron Data-efficient Image Transformer (SEMD), based on Vision Transformers, which classifies strongly lensed GW events by modeling morphological similarity between time-frequency spectrogram pairs. By integrating Squeeze-and-Excitation attention mechanisms and multilayer perceptrons , SEMD achieves strong feature extraction and discrimination. Trained and evaluated on simulated datasets using Advanced LIGO and Einstein Telescope noise, the model demonstrates robustness and generalization across different detector sensitivities and physical conditions, highlighting the promise of deep learning for rapid identification of strongly lensed GW signals.

The temporal and spatial coincidence between FRB 20200428 and hard peaks of the X-ray burst from SGR 1935+2154 suggests their potential association. We attributed them to the plasma synchrotron maser emission and synchrotron radiation of electrons in weakly magnetized, relativistically moving plasma blobs, and Monte Carlo simulation analysis shows that our model can predict observable fast radio burst outbursts and associated hard X-ray bursts with current telescopes. We constrained the properties of the blobs, including the Lorentz factor $\Gamma=5-30$, the magnetization factor $\sigma=6\times10^{-5}\sim 2\times 10^{-4}$, the electron Lorentz factor $\gamma_{\rm e,s}=(1.8-3.3)\times10^4$, and the plasma frequency $\nu_P=2.48 -42.61$ MHz. The inferred size of the blobs is $\sim 10^{9-10}$ cm, and it is located $\sim 10^{12-14}$ cm from the central engine. By adopting fine-tuned parameter sets, the observed spectra of both the FRB 20200428 outbursts and X-ray bursts can be well represented. The peak flux density ($F_{\rm\nu_{ pk}}$) of plasma maser emission is sensitive to $\sigma$ and $\nu_P$. Variation in $F_{\rm \nu_{pk}}$ can be more than 10 orders of magnitude, while the flux density of the synchrotron emission only varies by $1-2$ orders of magnitude. This can account for the observed sub-energetic radio bursts or giant radio pulses from SGR 1935+2154.

The detection of compact binary mergers with sub-solar masses at gravitational-wave observatories could mark the groundbreaking discovery of primordial black holes (PBHs). Concurrently, evidence for a nHz stochastic gravitational wave background observed by pulsar timing arrays (PTAs) could suggest a non-astrophysical origin, potentially arising from scalar-induced gravitational waves (SIGW). In this work, we analyze the connection between the two phenomena in the case where they share a common origin: the collapse of large primordial curvature perturbations in the early universe. We focus on sub-solar PBH populations within reach of upcoming experiments, including the current and future runs of LIGO-Virgo-KAGRA as well as the third generation observatories such as the Einstein Telescope and Cosmic Explorer. Using a Bayesian framework with physically motivated priors, we perform a consistent model comparison that incorporates existing astrophysical bounds together with the discovery potential of future detectors. Our analysis lends stronger support for the SIGW interpretation over the astrophysical one, as the narrowed priors place greater weight on the region of highest likelihood. Ultimately, we illustrate that combining PTA data with interferometer searches can deliver correlated evidence for new physics across multiple gravitational-wave bands.

NGC 6302 is a spectacular bipolar planetary nebula (PN) whose spectrum exhibits fast outflows and highly ionized emission lines, indicating the presence of a very hot central star (~220,000 K). Its infrared spectrum reveals a mixed oxygen and carbon dust chemistry, displaying both silicate and polycyclic aromatic hydrocarbon (PAH) features. Using the JWST Mid-Infrared Instrument (MIRI) and Medium Resolution Spectrometer, a mosaic map was obtained over the core of NGC 6302, covering the wavelength range of 5--28 micron and spanning an area of ~18.5 arcsec x15 arcsec. The spatially resolved spectrum reveals ~200 molecular and ionized lines from species requiring ionisation potentials of up to 205 eV. The spatial distributions highlight a complex structure at the nebula's centre. Highly ionized species such as Mg VII and Si VII show compact structures, while lower-ionization species such as H^+ extend much farther outwards, forming filament-defined rims that delineate a bubble. Within the bubble, the H^+ and H_2 emission coincide, while the PAH emission appears farther out, indicating an ionization structure distinct from typical photodissociation regions, such as the Orion Bar. This may be the first identification of a PAH formation site in a PN. This PN appears to be shaped not by a steady, continuous outflow, but by a series of dynamic, impulsive bubble ejections, creating local conditions conducive to PAH formation. A dusty torus surrounds the core, primarily composed of large (micron-sized) silicate grains with crystalline components. The long-lived torus contains a substantial mass of material, which could support an equilibrium chemistry and a slow dust-formation process.

We present GammaPBHPlotter, a public Python code for calculating and plotting the Hawking radiation gamma-ray spectra of primordial black holes in the mass range of $10^{14}$ to $10^{18}$ grams. This tool allows users to compute the monochromatic and mass-averaged spectra of black holes over a range of parameters. We include the primary/direct Hawking emission, the secondary emission from the decay and hadronization of unstable particles, the final state radiation, and the in-flight annihilation gamma-ray emission components.

The growing number of gravitational-wave detections from binary black holes enables increasingly precise measurements of their population properties. The observed population is most likely drawn from multiple formation channels. Population-synthesis simulations allow detailed modelling of each of these channels, and comparing population-synthesis models with the observations allows us to constrain the uncertain physics of binary black hole formation and evolution. However, the most detailed population-synthesis codes are computationally expensive. We demonstrate the use of normalising flows to emulate five different population synthesis models, reducing the computational expense, and allowing interpolation between the populations predicted for different simulation inputs. With the trained normalising flows, we measure the branching ratios of different formation channels and details of binary stellar evolution, using the current catalogue of gravitational-wave observations.

Reliable, autonomously, deployment of telescopes enables a wide range of possible science cases. In this paper, we present a method for multi-stage telescope alignment with a simple commercial imaging sensor. For these studies, we use a design of an example three mirror anastigmat telescope and consider how the average spot size across the detector changes as a function of primary and secondary mirror positioning. This multi-stage alignment procedure will consist of three subprocesses, starting with a coarse alignment and converging down to a finer alignment before moving on to a stage where the telescope will refine its misalignments for data acquisition. This alignment strategy has been tested and meets diffraction limited requirements on a subset of misalignment cases from a statistical Monte Carlo simulation given misalignment tolerances on the telescope.

Context. We investigate the binding energies of atoms to interstellar dust particles, which play a key role in their growth and evolution, as well as for the chemical reactions on their surfaces. Aims. We aim to compute the binding energies of abundant atoms in the interstellar medium (C, N, O, Mg, Al, Si, S, Ca, Fe, and Ni) to silicate dust. Methods. We used the Geometries, Frequencies, and Non-covalent Interactions Tight Binding (GFN1-xTB) method to compute the binding energies. An FeMgSiO$_4$ periodic surface model, containing 81 local minima on the surface, was used. Results. A range of binding energies was found for each element. The median of the binding energies follows the order Si (14.8 eV) > Al (12.8 eV) > Ca (12.7 eV) > C (9.5 eV) > O (8.1 eV) > N (6.2 eV) > Fe (6.0 eV) > S (5.2 eV) > Mg (2.4 eV). The probability distribution of binding energies for each element except Ca is statistically consistent with a log-normal distribution. Conclusions. In general, Si, Ca, and Al atoms have large binding energies. Thus, these atoms can stay on the silicate dust particles at high temperatures. The binding energies of the other atoms, C, N, O, Mg, S, Fe and Ni, are relatively weak. However, the computed binding energies for these elements are still far stronger than the energies associated with dust temperatures typical of the ambient interstellar medium (ISM), suggesting that silicate grains are generally stable against sublimation. We estimate sublimation temperatures for silicate grains to range from 1600 K to 3000K depending on assumed grain size and lifetime. These binding energies on silicate dust grains, estimated from first principles for the first time, provide invaluable input to models of dust evolution and dust-catalyzed chemical reactions in the interstellar medium.

Reconstructing the star formation history (SFH) of disk galaxies is central to understanding their growth and evolution, yet such estimates can be strongly biased by stellar radial migration - the redistribution of angular momenta over cosmic time. Using 186 Milky Way (MW) and Andromeda (M31) analogues from the TNG50 cosmological simulation, we compare star formation rates (SFRs) inferred from present-day stellar positions with those based on stellar birth radii to quantify the magnitude, spatial structure, and temporal evolution of migration-induced biases. We find that radial migration systematically produces artificial star formation in regions that had not yet formed stars. Notably, about 80% of galaxies exhibit outer-disk stars older than 10 Gyr, implying that these populations must have migrated from the inner disk. Similar effects appear in about 45% of galaxies at intermediate radii during early epochs, and in 30% of quenched inner disks within the past 4 Gyr. Migration also smooths SFHs, washing out localized bursts and periods of suppressed star formation by dispersing stars across neighboring radii. The strength and imprint of these distortions depend sensitively on galactic structure and evolutionary history: strong bars drive mean SFR overestimates of up to 75% in the inner disk and 150% in the outskirts; thinner, dynamically cold disks suffer average outer-disk biases up to 160%; while thick disks exhibit typical inner-disk biases up to 125%. Merger timing further modulates these patterns. Our results demonstrate that failing to account for stellar migration can lead to severe misinterpretations of when and where stars formed, with direct implications for the chemical and evolutionary histories of the MW and external galaxies.

We perform a series of simulations of magnetised Binary Neutron Star mergers, with varying magnetic field topologies in the initial data, as well as varying Equations of State, and mass ratios. In this paper, a companion paper to arXiv:2506.18995, we analyse the impact of the initial field configuration on the gravitational wave signal, the amplification of the magnetic field, and the ejected material. We investigate the dependence of the phase evolution of the gravitational wave in the post-merger on the initial magnetic field, finding that dephasing between the $(\ell=2,m=2)$ mode of the gravitational wave, and the $(2,1)$ and $(3,3)$ modes may be strongly impacted by the numerical reconstruction scheme. The magnetic field amplification during the Kelvin-Helmholtz dominated phase may be considerably enhanced by anti-aligned fields, or suppressed by toroidal fields. The post-merger amplification of the field due to winding may be suppressed by toroidal fields, and enhanced by asymmetries or mixtures of poloidal and toroidal fields. The field strength in the ejecta may be impacted by the initial magnetic field, with configurations which lead to large amplifications and those with mixtures of poloidal and toroidal fields preferentially emitting highly magnetised material in the polar regions, showing a weaker dependence of the magnetic field on the density of the ejecta than in cases that amplify the magnetic field less. We find that the magnetic field is largely randomly oriented in the ejected material, supporting such models used to estimate thermalisation timescales of ejected material. We find that configurations which begin with an initial bitant symmetry break this symmetry uniformly, independent of the initial configuration, when evolved without an enforced symmetry. This behaviour suggests the presence of a spontaneous symmetry breaking bifurcation in the solution.

As a benchmark for galaxy evolution and dark matter studies, the total mass of the Milky Way is a parameter of cosmological significance, and its value at large radii from the Galactic center remains highly uncertain. Following a hierarchical Bayesian inference approach, we measure the cumulative mass of the Milky Way using full 6D phase-space information of stars from the first data release of the Dark Energy Spectroscopic Instrument (DESI). We employ 330 blue horizontal-branch stars (BHBs) and 110 RR Lyrae stars (RRLs) in DESI covering Galactocentric distances in the range $\sim$50--100 kpc. Within 100 kpc from the Galactic center, we report an enclosed mass of $M(<100\ {\rm kpc}) = 0.57^{+0.08}_{-0.07}\times10^{12}$ M$_\odot$ and $M(<100\ {\rm kpc}) = 0.55^{+0.12}_{-0.10}\times10^{12}$ M$_\odot$ when using BHBs and RRLs, respectively. Extrapolating our mass profiles beyond the extent of our data, we find the virial mass of the Galaxy to be $M_{200}=0.85^{+0.16}_{-0.14}\times10^{12}$ M$_\odot$ and $M_{200}=0.78^{+0.19}_{-0.15}\times10^{12}$ M$_\odot$, respectively. We validate the effectiveness and limitations of our method using mock BHBs and RRLs from two AuriDESI halos. These tests show that the code recovers the enclosed mass of the mock galaxy with high precision and accuracy between 50 and 200 kpc, independent of the stellar tracer used and their spatial distribution. The tests also suggest an underestimation of the galaxy's cumulative mass at a level of up to $\sim20$\% if stars close to the Galactic center are used in the models. Our mass estimates lay the groundwork for future inference of the Galactic mass with upcoming DESI data releases and spectroscopic surveys mapping the halo.

We study the impact of the recent observation of the radio halo in the Coma galaxy cluster at 6.6 GHz with the Sardinia Radio Telescope on models based on turbulent reacceleration of electrons produced in dark matter annihilation processes. Observing at that frequency it is possible to obtain information on the electrons spectrum at energies where the effect of turbulent reacceleration becomes sub-dominant with respect to energy losses, and therefore to obtain information on the properties of seed electrons. Under the assumption that dark matter particles are neutralino-like particles annihilating at a rate close to the maximum allowed by Fermi-LAT upper limits in dwarf galaxies, we obtain some constraints on the intensity of the reacceleration and on the value of the neutralino mass. In particular, models with mass of the order of 10 GeV are generally disfavored, because they produce a high-frequency radio spectrum that can not reproduce the possible flattening observed between 5 and 6.6 GHz; on the other hand, models with mass of the order of 500 GeV, in order to reproduce the observed spectrum at frequencies below 100 MHz, require a reacceleration phase longer than $10^9$ yr, which would require more than one event responsible of the generation of turbulence in the cluster. The resulting optimal mass values are in the range 100--200 GeV, with a preference for the quark annihilation channel.

Studying the emergence of magnetic fields is essential for understanding the physical mechanisms behind various phenomena in the solar atmosphere. Most importantly, the emerging fields offer valuable insights into how energy and mass are transferred to the upper solar atmosphere. As a result, they have garnered significant attention from both observational and theoretical perspectives. In this article, we present two models of quiet-Sun-like magnetic fields generated by the Bifrost code. We compare these models with observations from the Swedish 1-meter Solar Telescope (SST) and the Interface Region Imaging Spectrograph (IRIS). By tracking the magnetic features in both the SST and Bifrost data, we determine the similarities and differences between the fields identified in the models and those observed. We conduct a quantitative comparison of various properties, such as flux content, flux densities, horizontal and line-of-sight velocities, lifetimes, sizes, and surface interactions. Additionally, we identify and analyze the properties of the largest emerging bipoles in the SST and Bifrost data. Our findings indicate that the magnetic bipoles in the Bifrost simulations are generally stronger than those observed with the SST. However, a qualitative comparison of the chromospheric and transition region responses to the emerging fields in the Bifrost models, SST, and IRIS observations shows similar heating processes occurring above and around the emerging fields. Finally, we outline our plans for future work aimed at studying the emergence of magnetic fields in the quiet Sun, with a particular focus on the chromosphere and upper atmospheric layers.

Local helioseismology comprises of imaging and inversion techniques employed to reconstruct the dynamic and interior of the Sun from correlations of oscillations observed on the surface, all of which require modeling solar oscillations and computing Green's kernels. In this context, we implement and investigate the robustness of the Hybridizable Discontinuous Galerkin (HDG) method in solving the equation modeling stellar oscillations for realistic solar backgrounds containing gravity and differential rotation. While a common choice for modeling stellar oscillations is the Galbrun's equation, our working equations are derived from an equivalent variant, involving less regularity in its coefficients, working with Lagrangian displacement and pressure perturbation as unknowns. Under differential rotation and axisymmetric assumption, the system is solved in azimuthal decomposition with the HDG method. Compared to no-gravity approximations, the mathematical nature of the wave operator is now linked to the profile of the solar buoyancy frequency N which encodes gravity, and leads to distinction into regions of elliptic or hyperbolic behavior of the wave operator at zero attenuation. While small attenuation is systematically included to guarantee theoretical well-posedness, the above phenomenon affects the numerical solutions in terms of amplitude and oscillation pattern, and requires a judicious choice of stabilization. We investigate the stabilization of the HDG discretization scheme, and demonstrate its importance to ensure the accuracy of numerical results, which is shown to depend on frequencies relative to N, and on the position of the Dirac source. As validations, the numerical power spectra reproduce accurately the observed effects of the solar rotation on acoustic waves.

We present a high-resolution cosmological zoom-in simulation of a group of field dwarf galaxies which includes on-the-fly radiative transfer (RT) and is evolved to $z=0$. Emission from stars is included according to age-dependent spectra, and a redshift-dependent UV background. The inclusion of RT results in the formation of eight additional faint dwarf galaxies with stellar masses of $10^{4}$ M$_{\odot}$ to several $10^{5}$ M$_{\odot}$ and only old stellar populations, similar to the observed Ultra-Faint Dwarf galaxies. They formed before and during cosmic reionisation and were mostly quenched by $z \sim 3-4$. The simulated galaxies follow many observed scaling relations such as the stellar mass-halo mass relation, the mass-size relation, and the luminosity-velocity dispersion relation. For the more massive dwarf galaxies, radiative feedback suppresses star formation, making it less bursty and reducing explosive outflows. This consequently reduces the dark matter core sizes by a factor of 2-3, rendering the core sizes ($\sim$ 1 kpc) more consistent with observations. The distribution of HI in the circumgalactic medium (CGM) is ubiquitous with a covering fraction of unity within $R_{vir}$, in good agreement with observations. It is rather insensitive to radiative or SN feedback at $z=0$, but at $z>5$ it is much higher in the RT simulation. In contrast, the distribution of low ions like SiII is very compact and declines sharply beyond the ISM scale. CIV and OVI have a more extended distribution, but their column densities are generally below the detection limit. Radiative feedback leads to smaller column densities of the metal ions, partly due to the reduction of total metal production, and partly because hard photons from the stellar radiation escape the ISM and further ionise the CGM. The abundance of CIV is particularly sensitive to the latter effect.

The Habitable Worlds Observatory (HWO) is a flagship mission concept proposing to characterize earth-like exoplanets at high contrast with a coronagraph instrument. The most in-depth, validated contrast error budgets made to date have been at the $10^{-9}$ planet-to-star contrast levels for the Roman Space Telescope. To obtain the raw contrast levels of $\leq 10^{-10}$ contrast needed for HWO, more modeling is needed of the vector diffraction effects within a coronagraph instrument. Several coronagraph architectures are based on the Lyot coronagraph, which utilizes a small occulting spot to block most of the on-axis starlight at a focal plane. In this paper, we present a study of the polarization aberrations generated by internal occulting spots. First, we use finite-difference time-domain (FDTD) modeling to produce the complex transmission of the occulter in each polarization state, and then we separately run those through a closed-loop wavefront control model in FALCO\cite{rigss2018falco} to determine the degradation in achievable contrast. We perform sweeps of such parameters as spot diameter, spot thickness, and angle of incidence to inform the HWO instrument design.

The galaxy LEDA 1313424 has recently been discovered to have nine rings [1]. The authors of the discovery paper interpret the rings to be the result of the passing of a smaller galaxy through the center of the LEDA galaxy 56 Myr ago. We point out several difficulties with this interpretation and propose instead that the rings of LEDA are the imprint upon baryonic matter of caustic rings of dark matter. Caustic rings of dark matter, with the pattern of radii observed in the LEDA galaxy, form if the dark matter is axions. Because they move slowly (LEDA's 70 kpc ring, for example, has speed 6.6 km/s) caustic rings of dark matter are plausible sites for the star formation observed in LEDA.

Context: This research presents unsupervised machine learning and statistical methods to identify and analyze tidal tails in open star clusters using data from the Gaia-DR3 catalog. Aims: We aim to identify member stars, detect and analyze tidal tails in five open clusters, BH 164, Alessi 2, NGC 2281, NGC 2354, and M67, of ages between 60 Myr and 4 Gyr. These clusters were selected based on the previous evidence of extended tidal structures. Methods: We utilized machine learning algorithms such as Density-Based Spatial Clustering of Applications with Noise (DBSCAN) and Principal Component Analysis (PCA), along with statistical methods to analyze the kinematic, photometric, and astrometric properties of stars. Key characteristics of tidal tails, including radial velocity, color-magnitude diagram, and spatial projections in the tangent plane beyond the cluster's Jacobi radius ($r_J$), were used to detect them. We used N-body simulations to visualize and compare the observables with real data. Further analysis was done on the detected cluster and tail stars to study their internal dynamics and populations, including binary fraction. We also applied the residual velocity method to detect rotational patterns in the clusters and their tails. Results: We identified tidal tails in all five clusters, with detected tails extending farther in some clusters and containing significantly more stars than previously reported (tails ranging from 40 to 100 pc, one to four times the Jacobi radius, with 100 - 200 tail stars). The luminosity functions of the tails and their parent clusters were generally consistent, and tails lacked massive stars. In general, the binary fraction was found to be higher in the tidal tails. Significant rotation was detected in M67 and NGC 2281 for the first time.

Stellar population properties are crucial for understanding galaxy evolution. Their inference for statistically representative samples requires deep multi-object spectroscopy, typically obtained with fiber-fed spectrographs that integrate only a fraction of galaxy light. The most comprehensive local Universe dataset is the Sloan Digital Sky Survey (SDSS), whose fibers typically collected ~30% of total flux. Stellar population gradients, ubiquitously present in galaxies, systematically bias SDSS toward central properties, by amounts yet to be quantified. We leverage CALIFA integral-field spectroscopy to simulate fiber-fed observations at redshifts z=0.005-0.4, accounting for seeing effects. We analyze systematic aperture correction trends across galaxy morphologies and derive correction recipes based on: fiber-measured indices, global g-r color, absolute r-band magnitude Mr, and physical half-light radius R50. Corrections for absorption indices typically reach >~15% of their dynamical range at z~0.02, decreasing to ~7% at z~0.1 (median SDSS redshift) and becoming negligible above z~0.2. Spiral galaxies exhibit the largest aperture effects due to their strong internal gradients. Our correction recipes, applied to the SDSS-DR7 dataset, significantly reduce scatter in stellar population diagnostic planes and enhance bimodality in age-sensitive diagrams. Corrections reveal systematic overestimates of old galaxy fractions by up to 10% and an underestimate by >~0.2 mag of the transition luminosity at which old galaxies become dominant. Aperture corrections significantly impact observational tracers of stellar populations from fiber spectroscopy. Absorption indices corrections applied to SDSS-DR7 will provide a robust local benchmark for galaxy evolution studies.

Cosmological models of hierarchical structure formation predict the existence of a widespread population of dual accreting supermassive black holes (SMBHs) on kpc-scale separations, corresponding to projected distances < 0".8 at redshifts higher than 0.5. However, close companions to known active galactic nuclei (AGN) or quasars (QSOs) can also be multiple images of the object itself, strongly lensed by a foreground galaxy, as well as foreground stars in a chance superposition. Thanks to its large sky coverage, sensitivity, and high spatial resolution, Euclid offers a unique opportunity to obtain a large, homogeneous sample of dual/lensed AGN candidates with sub-arcsec projected separations. Here we present a machine learning approach, in particular a Convolutional Neural Network (CNN), to identify close companions to known QSOs down to separations of $\sim\,$0".15 comparable to the Euclid VIS point spread function (PSF). We studied the effectiveness of the CNN in identifying dual AGN and demonstrated that it outperforms traditional techniques. Applying our CNN to a sample of $\sim\,$6000 QSOs from the Q1 Euclid data release, we find a fraction of about 0.25% dual AGN candidates with separation $\sim\,$0".4 (corresponding to $\sim$3 kpc at z=1). Estimating the foreground contamination from stellar objects, we find that most of the pair candidates with separation higher than 0".5 are likely contaminants, while below this limit, contamination is expected to be less than 20%. For objects at higher separation (>0".5, i.e. 4 kpc at z=1), we performed PSF subtraction and used colour-colour diagrams to constrain their nature. We present a first set of dual/lensed AGN candidates detected in the Q1 Euclid data, providing a starting point for the analysis of future data releases.

Carbon-deficient giants (CDGs) are a rare and chemically peculiar class of stars whose origins remain under active investigation. We present an asteroseismic analysis of the entire known CDG population, selecting 129 stars observed by $Kepler$, K2, and TESS to obtain seismic constraints. We detect solar-like oscillations in 43 CDGs. By measuring $\nu_{\rm max}$ and applying seismic scaling relations, we determine precise masses for these stars, finding that 79\% are low-mass ($M \lesssim 2~M_\odot$). The luminosity distribution is bimodal, and the CDGs separate into three chemically and evolutionarily distinct groups, characterized by clear trends in sodium and CNO abundances, $\alpha$-element enhancement, and kinematics. We find that two of these groups are only distinguished by their initial $\alpha$-element abundances, thus effectively reducing the number of groups to two. Lithium enrichment is common across all groups, linking CDGs to lithium-rich giants and suggesting a shared evolutionary origin. We find that spectroscopic $\log g$ is systematically offset from seismic values. Group~1 CDG patterns are most consistent with formation through core He-flash mixing, while the more massive and more chemically processed Groups~2 and 2$\alpha$ likely formed through mergers involving helium white dwarfs, possibly in hierarchical triples. Pollution from AGB stars appears very unlikely, given the unchanged [C+N+O] abundance across all groups.

Radio relics, arc-like polarized sources with highly aligned magnetic fields, are typically found on the outskirts of merging galaxy clusters. The magneto-ionic media responsible for the significant coherence observed in radio relics remain poorly understood. Low-frequency measurements of radio relics are essential for constraining depolarization models, which provide crucial insights into the magnetic field distribution. However, these measurements are challenging due to the emission properties and interferometer systematics. We have detected polarization signals from the northwest radio relic in Abell 746 at 650 MHz with the upgraded Giant Meterwave Radio Telescope, marking the first-ever detection of polarisation from radio relics below 1 GHz. At this frequency, the average Rotation Measure (RM) corrected magnetic fields align well with shock radio emission, typical of radio relics. The fractional polarization at 650 MHz is $\sim 18\pm4 \%$. Our results indicate that a single internal depolarization model cannot explain the observed depolarization spectra, suggesting a non-uniform magnetic field distribution or complex contribution of different polarized regions in the radio relic. Our detection of polarization signals at 650 MHz reveals critical insights into radio relic magnetic field structures, offering a low-frequency approach to understanding ICM magnetic fields in merging galaxy clusters.

Coherent, periodic radio emission from pulsars has been widely interpreted as evidence of neutron stars as strongly magnetized compact objects. In recent years, radio pulses have also been detected from white dwarfs (WDs) in tight binary systems, raising the question of whether isolated WDs could similarly host pulsar-like emission. We conducted the most sensitive search to date for coherent radio signals from five isolated, rapidly rotating, and magnetized WDs, using the Five-hundred-meter Aperture Spherical radio Telescope (FAST), the Green Bank Telescope (GBT), and the Australia Telescope Compact Array (ATCA). No pulsed or continuum radio emission was detected down to $\mu$Jy levels. These non-detections place the most stringent observational constraints yet on the existence of isolated WD pulsars. Our results suggest that either such emission is intrinsically weak, narrowly beamed, or requires binary-induced magnetospheric interactions absent in solitary systems. Comparison with the known radio-emitting WDs highlights the critical role of companion interaction in enabling detectable emission. This work expands on prior surveys by targeting sources with the most favorable physical conditions for WD pulsar-like activity and employing highly sensitive, targeted observations. Future observations with next-generation facilities such as the SKA will be essential to explore fainter or sporadic emission from massive, magnetic WDs and to investigate their potential as compact radio transients further.

Recent observations with the Five-hundred-meter Aperture Spherical Telescope (FAST) have revealed abundant reservoirs of neutral hydrogen (HI) in low redshift post-starburst galaxies (PSBs), raising the question of why star formation ceases rapidly in these systems. In this study, we present a detailed analysis of the shape of the integrated HI spectra of 67 PSBs. We find that PSBs exhibit significantly higher HI spectral concentration values ($K$) compared to a matched sample from xGASS, and are more comparable to those of starburst galaxies. By extending our analysis to spatially resolved HI data from THINGS and ATLAS$^{\rm 3D}$, we show that both centrally concentrated HI distributions and dynamically unsettled HI can effectively increase $K$, while non-axisymmetric structures only contribute to the scatter of the $K$ distribution. Distinguishing between central concentration and dynamically unsettled gas as the origin of high $K$ can be achieved by measuring the spectral asymmetry ($A_{\rm F}$), making the $K$-$A_{\rm F}$ plane a powerful diagnostic tool for identifying galaxies with unsettled HI using integrated spectra alone. Based on their location in the $K$-$A_{\rm F}$ plane, we find that most PSBs are not dominated by unsettled HI, but rather exhibit elevated central gas concentration. Both modes of gas redistribution in PSBs may eventually contribute to their quenching.

Low-Earth Orbit (LEO) satellites are crucial for communications, navigation, and Earth observation. However, their operational lifetimes are strongly influenced by orbital decay due to atmospheric drag. This work presents a simplified analytical model of orbital decay in LEO, incorporating atmospheric density variations and satellite parameters. Using standard drag force equations and density models, we explore how altitude, satellite mass, and cross-sectional area affect decay timescales. Comparisons are made with existing NASA models and previous literature. The results highlight the importance of accurate atmospheric density representation and solar activity in predicting satellite lifetimes, especially relevant in the context of increasing space debris and mega-constellations.

We present a potential test of the origin of the $\gamma$-ray Galactic Centre Excess (GCE). We demonstrate how gravitational microlensing by stellar mass objects along the line of sight to the Galactic Bulge can distinguish between the possibility of extensive emission due to dark matter self-annihilation from more prosaic astrophysical sources, namely millisecond pulsars. Such an astrophysical origin would result in emission from a population of small, currently unresolved point-like sources - in contrast to the expected smoother emission resulting from dark matter annihilation. Given that the scale of gravitational microlensing, that is, the Einstein radius for stellar mass lenses, and hence, the degree of induced magnification, is sensitive to the size of the emitting region, such microlensing will induce time variability in the emission of astrophysical sources, whereas $\gamma$-ray emission from dark matter annihilation will effectively be immune to such influences. However, we find that detecting microlensing-induced variability requires significantly greater sensitivity than that of current or planned $\gamma$-ray detectors. For a small population of bright GCE sources, more than an order-of-magnitude increase in effective area over Fermi-LAT would be required, with events remaining extremely rare. For a large population of faint sources, events would occur multiple times a year, but would only be detectable with a four-order-of-magnitude improvement. Whilst microlensing might not be a definitive test of the origin of the GCE, in future observations, it may prove useful in determining the properties of any point-like source population.

We discovered a hyperluminous dust-obscured galaxy with mysterious blue-excess emission (BlueDOG) in rest-frame UV of its spectral energy distribution (SED) from a multi-wavelength survey in the AKARI Deep Field - South (ADF-S). We present the results of SED analysis with multiwavelength photometric data and spectroscopic analysis, observed with Gemini-S/GMOS, FLAMINGOS-2, to explore the origin of blue-excess emission of a hyperluminous BlueDOG, ADFS-KMTDOG-102, at z=2.6. The SED analysis shows that this BlueDOG is a highly massive system (log $M_{*}$/$M_\odot=12.3$) with substantial extinction. Additionally, the proportion of the old stellar population exceeds that of the young stellar population, which suggests stellar evolution cumulated from the early universe. The mass of supermassive black hole (SMBH) estimated using the extinction-corrected broad H$\rm\alpha$ emission line yields log $M_{\rm BH}$/$M_\odot$=10.2. We discuss the similarity between the BlueDOG and 'Little Red Dots' (LRDs), recently discovered with the James Webb Space Telescope, showing SED shapes remarkably similar to those of LRDs. The UV emission line ratios indicate that the emission lines are primarily powered by the central active galactic nuclei (AGN). In contrast, the origin of the blue-excess UV continuum remains ambiguous, since both recent star formation and AGN-induced scattered light are viable explanations, based on the results from the SED fitting and scattered light modeling.

We report the discovery of a vertical velocity anomaly (VVA) in the stellar component of the Galactic disk, consistent with the impact of a dark, ultra-compact dwarf galaxy (UCD)-sized object plunging into the Milky Way. The anomaly spatially coincides with a suite of gaseous disturbances -- including an \Hone\ void, a molecular shell (CO 16.134--0.553), and a vertical \Hone\ filament -- previously interpreted as signatures of a dark matter subhalo (DMSH) collision. Analysis of Gaia DR2 astrometry reveals a statistically significant vertical velocity dip co-located with these features, supporting a dynamical origin. The absence of a luminous source at the filament's terminus suggests a dark or failed UCD. These results provide rare observational evidence for a low-mass DMSH, with implications for the substructure of dark matter halos and constraints on the $\Lambda$CDM model.

We present a multi-wavelength study on the mini-outburst of the black hole X-ray binary MAXI J1348-630 to investigate the origin of its optical emission. Using the X-ray data from Insight-HXMT and the optical observations from the Las Cumbres Observatory, we find a powerlaw correlation between the optical and X-ray fluxes with slopes of about 0.4, indicating a hybrid contribution to the optical emission from both viscously heated disk and irradiation of the outer disk illuminated by the inner X-ray emission, but inconsistent with jet-dominated optical emission. Time delay analysis shows that the optical emission precedes the X-ray by about 8.5 days, which can be naturally interpreted by the disk instability model (DIM). Meanwhile, spectral energy distribution (SED) fitting with an irradiated disk model successfully reproduces the observed optical and X-ray emission, with negligible contribution from the jet. The color-magnitude diagram further supports the disk-dominated optical origin of MAXI J1348-630 during its mini-outburst. Our results suggest that the DIM plays a central role in understanding mini-outbursts of X-ray binaries.

Protoplanetary nebulae (PPNs) represent a critical evolutionary bridging stage between the asymptotic giant branch (AGB) and planetary nebula (PN) stages. Their dynamical structures provide key insights into late stellar evolution. Here, we report CO and $^{13}$CO imaging observations of the carbon-rich PPN IRAS 06530$-$0213, a source exhibiting the unidentified 21 micron emission band, conducted with the Northern Extended Millimeter Array (NOEMA). The CO maps reveal a spherical halo (diameter $\sim10''$) surrounding a central barrel-shaped torus, where the torus displays an inner diameter of $\sim1.5''$ and an outer diameter of $\sim2''$. Through three-dimensional morpho-kinematic modeling with the SHAPE software, we determine that IRAS 06530$-$0213 experienced its final thermal pulse during the AGB phase $\sim6500$ years ago, transitioning into the PPN phase $\sim4500$ years ago. Our analysis indicates that the blue-shifted CO emission, a feature also detected in several other PPNs and PNs previously misattributed to the interstellar clouds, actually originates from the obscuration of the central nebula by the remnant AGB halo. These findings are expected to deepen our understanding of the dynamical structures of PPNs, as well as their pivotal transitional role in the late-stage evolution of low-to-intermediate mass stars.

We investigate the structural and asteroseismic consequences of mass accretion in massive stars within close binary systems. Using MESA, we model the evolution of the 10 M$_{\odot}$ accretor through and after a Roche lobe overflow phase. In addition to changing the surface composition of the star, mass accretion also significantly modifies the internal structure by expanding the convective core and altering chemical stratification near the core-envelope boundary. This partial core rejuvenation creates a distinct mean molecular weight gradient and leaves a persistent local density modulation. In the late stages of mass transfer, changes in density and sound-speed profiles become apparent and influence stellar oscillations. We analyse the asteroseismic properties of the post-mass transfer models compared to single stars of the same mass and central hydrogen abundance. In the gravity mode regime, the altered Brunt-Väisälä frequency leads to period spacing patterns with larger amplitudes and phase shifts. For low and intermediate-order pressure modes, we find systematic frequency deviations linked to changes in the sound-speed profile. Weight function analyses confirm that these differences arise primarily from structural modifications near the convective core boundary. Furthermore, small frequency separations, sensitive to localised sound-speed gradients, reveal periodic variations attributable to the density discontinuity at the convective core edge. The accretor exhibits a larger sound-speed gradient integral and a longer acoustic radius ratio compared to the single star, consistent with its expanded core. Our results demonstrate that mass accretion imprints measurable asteroseismic signatures on both gravity and pressure modes. These signatures provide powerful diagnostics for identifying post-interaction stars, for refining stellar age and structure estimates in binary systems.

Accretion variability is ubiquitous in YSOs. While large outbursts may strongly affect the disk, the effects of moderate bursts are less understood. We aim to study the physical response of the disk around the eccentric binary system DQ Tau to its periodic accretion changes. We organized a multi-wavelength campaign centered on four JWST/MIRI spectra. We targeted three periastrons (high accretion state) and one apastron (quiescence). We used optical and near-infrared spectroscopy and photometry to measure how the accretion luminosity varies. We decomposed the multi-epoch SEDs into stellar, accretion, and rim components. We fitted the solid-state features using various opacity curves and the molecular features using slab models. We find the inner disk of DQ Tau to be highly dynamic. The temperature, luminosity, and location of the inner dust rim vary in response to the movement of stars and the L_acc variations. This causes variable shadowing of the outer disk, leading to an anti-correlation between the rim temperature and the strength of the silicate feature. The dust mineralogy remains constant, dominated by large amorphous olivine and pyroxene grains, with smaller fractions of crystalline forsterite. The excitation of CO (1550-2260 K), HCN (880-980 K), and hot H2O (740-860 K) molecules as well as the luminosity of the [NeII] line correlate with the accretion rate, while the warm (650 K) and cold (170-200 K) H2O components are mostly constant. CO emission, originating from a hot (>1500 K) region likely within the dust sublimation radius, is most sensitive to L_acc changes. In comparison with other T Tauri disks, DQ Tau is highly C-poor and displays moderately inefficient pebble drift. We conclude that even moderate accretion rate changes affect the thermal structure in the planet-forming disk regions on short timescales, providing a crucial benchmark for understanding disk evolution.

Small-pitch populated focal planes are a critical enabling technology for the next generation of highly multiplexed astronomical instruments. This work reports on ongoing prototyping activities for 6.2 mm-pitch alpha-beta (theta-phi) fiber positioner modules. Following a concise review of the underlying concept, the performance of trillium-based and independently actuated robotic designs is compared. Key operational parameters including XY positioning repeatability, non-linearity, backlash, and tilt are quantitatively assessed, with emphasis on their implications for overall instrument performance. Initial prototyping results are highly encouraging, indicating performance levels approaching the preliminary design specifications.

The behavior of lithium (Li) in Population I main sequence stars challenges standard stellar theory. Two phenomena stand out: the solar Li problem which extends to Li depletion in solar analogs and the Li dip observed in mid-F stars within open clusters. Building on the meridional circulation-driven radial mixing framework previously developed to explain Li-enriched red clump stars, we explore its relevance to Li depletion on the main sequence. First, our models reproduce the observed $A(\text{Li})$-Age correlation in solar analogs. Through detailed isochrone analysis, we find good agreement between the simulated and observed $A(\text{Li})$-$T_{\text{eff}}$ relationships within the solar analog parameter space. However, the predicted solar Li abundance ($\sim 1.5\,\text{dex}$) is still higher than current solar measurements. Second, our models partially explain the Li dip phenomenon in mid-F cluster stars. The models accurately reproduce Li distributions on the cool side of the Li dip in most clusters and capture the Li behaviors on the hot side observed in systems like the Hyades. However, we identify limitations in the models' ability to fully reproduce the dip morphology, particularly due to the rotation velocity distribution of sample stars in this temperature zone.

We present a novel analytical framework employing Physics-Informed Neural Networks (PINNs) to constrain the cosmological constant $\Lambda$ through the analysis of stellar orbits around the supermassive black hole (SMBH) Sgr A* at the Galactic center. Focusing on the well-observed S2 star, we use an inverse PINN (iPINN) architecture to infer orbital elements and estimate the total precession angle from astrometric data. By isolating the contribution from $\Lambda$, which is defined as the difference between the total precession and the Schwarzschild precession, we derive a stringent upper bound of $\Lambda \leq 5.67 \times 10^{-40}, \mathrm{m}^{-2}$, which is approximately two orders of magnitude tighter than previous estimates obtained using similar data-driven methods. Extension of our analysis to two additional long-period S-stars, S1 and S9, reveals that while the cosmological precession becomes relatively more prominent in such systems, limited orbital coverage introduces significant uncertainties in parameter estimation. Among the cases examined, the constraint derived from S2 remains the most robust. Our results highlight the potential of PINN-based approaches for extracting physical insights from sparse or noisy astronomical data. Future applications to next-generation observational data and further methodological improvements in machine learning are expected to refine the cosmological constraints and enable broader tests of gravitational theories.

Climate change is inducing a global atmospheric contraction above the tropopause (~10 km), leading to systematic decrease in neutral air density. The impact of climate change on small meteoroids has already been observed over the last two decades, with documented shifts in their ablation altitudes in the mesosphere (~50-85 km) and lower thermosphere (~85-120 km). This study evaluates the potential effect of these changes on meteorite-dropping fireballs, which typically penetrate the stratosphere (~10-50 km). As a case study, we simulate the atmospheric entry of the fragile Winchcombe carbonaceous chondrite under projected atmospheric conditions for the year 2100 assuming a moderate future emission scenario. Using a semi-empirical fragmentation and ablation model, we compare the meteoroid's light curve and deceleration under present and future atmospheric density profiles. The results indicate a modest variation of the ablation heights, with the catastrophic fragmentation occurring 300 m lower and the luminous flight terminating 190 m higher. The absolute magnitude peak remains unchanged, but the fireball would appear 0.5 dimmer above ~120 km. The surviving meteorite mass is reduced by only 0.1 g. Our findings indicate that century-scale variations in atmospheric density caused by climate change moderately influence bright fireballs and have a minimal impact on meteorite survival.

Many galactic globular clusters (GCs) contain at least two stellar populations. Recent observational studies found that the radial distributions of the first (P1) and second population (P2) differ in dynamically-young GCs. Since P2 is conventionally assumed to form more centrally concentrated, the rapid mixing (or even inversion) in some GCs but not others is puzzling. We investigate whether dynamical processes specific to certain GCs might cause this. Specifically, we evaluate the expansion of P2 by binary-single interactions in the cluster core and whether these can mix the P1/P2 radial distributions, using a set of toy-models with varying numbers and masses of primordial binaries. We find that even one massive binary star can push the central P2 outwards, but multiple binaries are required to fully mix P1 and P2 within a few relaxation times. We also compare our results to observed properties of mixed young GCs (NGC 4590, 5053, or 5904).

The star-forming complex of Corona Australis (CrA) is one of the closest and most isolated molecular clouds. It belongs to a chain of clusters that show age gradients with distance from the galactic plane. We aim to provide suggestions regarding its formation history by examining the stellar and disk populations, stellar multiplicity, and interstellar absorption. We made a census of stars and disks using Gaia DR3 and infrared data. Interstellar absorption in the direction of each star was derived by comparing SpTy from the literature and Gaia colors. Stellar multiplicity analysis accounts for both direct observation of visual companions (Gaia data and high-contrast imaging) and indirect detection of the presence of companions (eclipsing and spectroscopic binaries, and astrometry). The properties of the disks were obtained from the slopes of the spectral energy distributions. As found in previous studies, the CrA complex can be divided into two regions: a younger region (CrA-Main: 3+/-1Myr) and an older one (CrA-North: 6.7+/-0.3Myr), slightly younger than previously thought. While CrA-Main still appears bound to the gas, CrA-North is unbound and expanding. The stars that belong to CrA-North were in the most compact configuration 3.72Myr ago. At that time, CrA-Main and CrA-North were much closer to each other than they appear now. The fraction of disk-bearing stars is higher in Main than in North, as also expected due to the younger age of CrA-Main. We propose a formation history scenario for the CrA-complex. It started between 15 and 18 Myr ago with SNe explosions in the Upper Centaurus-Lupus complex, followed by a quiescent phase with little star formation. A star formation episode about 7Myr ago formed CrA-North stars. About 3.7Myr ago, a second SN explosion south of CrA-North triggered star formation in CrA-Main. This last SN might have been the origin of the pulsar RX J1856.5-3754.

Asteroseismology of solar-like oscillations in giant stars allow the derivation of their masses and radii. For members of open clusters, this provides an age of the cluster that should be identical to the one derived from the colour-magnitude diagram, but independent of the uncertainties that are present for that type of analysis. We aim to identify and measure the properties of giant members of the open cluster NGC6633, and combine these with asteroseismic measurements to derive a precise and self-consistent cluster age. Importantly, we wish to constrain the effects of rotational mixing on stellar evolution, since assumptions about internal mixing can have a significant impact on stellar age estimates. We identify five giant members of NGC6633 using photometry and Gaia data, supplemented by spectroscopic literature measurements. These are combined with asteroseismic measurements from CoRoT data and compared to stellar-model isochrones. This constrains the interior mixing to a low level and enables the most precise, accurate and self-consistent age estimate so far for this cluster. Asteroseismology of the giants and the cluster colour-magnitude diagram provide self-consistent masses of the giant members and their radii constrain the stellar interior mixing to a low level. The [C/N] ratios and Li abundances also suggest that rotation has had very little influence on the evolution of the stars in NGC 6633. This results in an age estimate of 0.55+0.05-0.10 Gyr for NGC6633. Four giant members appear to be in the helium-core burning evolutionary phase as also expected from evolutionary timescales. For the largest giant, the evidence remains inconclusive. A comparison to other age and mass estimates for the same stars in the literature uncovers biases for automated age estimates of helium-core burning stars.

The eight ultra-long period pulsars (ULPPs) in radio bands have been discovered recently, e.g., ASKAP J1935+2148 with a spin period of 53.8\,min, which are much longer than those of normal pulsars, spanning from 0.016\,s to 23.5\,s, however the origins, spin evolutions and emission mechanisms of these sources are still puzzling. We investigate how the ultra-long period of ASKAP J1935+2148 is evolved by the braking of relativistic particle wind, in a time scale of about 0.1 - 1 Myr, from a normal pulsar with local superstrong magnetic fields. In addition, it is noticed that the ULPPs in the period versus period derivative diagram are much below the ``death line", implying their different characteristics from the normal pulsars. Five sources (including ASKAP J1935+2148) in total eight ULPPs share the rotational energy loss rates to be lower than their respective radio emission luminosities, a phenomenon that can be accounted for by the sustainable radio bursts induced through the reconnection of locally concentrated magnetic field this http URL diversity and complexity of ULPP radio emissions should be closely related to the presence of magnetic reconnection rather than rotational powered discharges in the gaps. Furthermore, it is suggested that the coherent radio emissions of pulsars may have two origins, one from the rotation-powered electric voltage that accounts for the normal pulsar phenomena and the other from the magnetic reconnection-induced continual radio bursts that account for the ULPP observations.

The heating and structure of the solar chromosphere depends on the underlying magnetic field, among other parameters. The lowest magnetic flux of the solar atmosphere is found in the quiet Sun internetwork and is thought to be provided by the small-scale dynamo (SSD) process. We aim to understand the chromospheric structure and dynamics in a simulation with purely SSD generated magnetic fields. We perform a 3D radiation-magnetohydrodynamic (rMHD) simulation of the solar atmosphere, including the necessary physics to simulate the solar chromosphere. No magnetic field is imposed beyond that generated by an SSD process. We analyse the magnetic field in the chromosphere, and the resulting energy balance. Plasma at chromospheric temperatures reaches high into the atmosphere, with small, transient regions reaching coronal temperatures. An average Poynting flux of $5\times10^6~\mathrm{erg\;cm}^{-3}$\;s$^{-1}$ is found at the base of the chromosphere. The magnetic field in the chromosphere falls off more slowly with height than predicted by a potential field extrapolation from the radial component of the photospheric field. Starting in the middle chromosphere, the magnetic energy density is an order of magnitude larger than the kinetic energy density and, in the upper chromosphere, also larger than the thermal energy density. Nonetheless, even in the high chromosphere, the plasma beta in shock fronts and low-field regions can locally reach values above unity. The interactions between shocks and the magnetic field are essential to understanding the dynamics of the internetwork chromosphere. The SSD generated magnetic fields are strong enough to dominate the energy balance in the mid-to-upper chromosphere. The energy flux into the chromosphere is $8.16\times 10^{6}~\mathrm{erg\;cm^{-2}\;s^{-1}}$, larger than the canonical values required to heat the quiet sun chromosphere and corona.

We present the first detections of gas around a 2.1 Myr old debris disc-bearing intermediate-mass star, HD 44892. Through our analysis we deem the alternative age of 800 Myr unlikely. Gas is detected both in $^{12}$CO (2-1) emission through ALMA Band 6 observations and absorption in Ca II K and H seen with high-resolution UVES spectroscopy. The star exhibits a 12 $\mu$m fractional excess of $7.86^{+0.11}_{-2.27}$, placing it in the transition stage between protoplanetary and debris discs, at a level comparable to HD 141569. Our detection of 1.3 mm emission yields the dust mass of 0.019$\pm$0.009 $M_{\oplus}$ assuming 115 K temperature. Assuming a gas temperature of 20 K, the resulting CO gas mass is (7.86$\pm$2.05)$\times$10$^{-5}$ $M_{\oplus}$. For a warmer gas of 50 K it increases to (1.62$\pm$0.17)$\times$10$^{-4}$ $M_{\oplus}$. From the $^{13}$CO and C$^{18}$O non-detections we derive upper limits which suggest this CO gas mass estimate should be treated as a lower limit, because the CO emission might be optically thick. From the analysis of UVES spectroscopic data we find variability in the spectrum of HD 44892, in both calcium Ca II K and H lines, which is absent from sodium Na I D1 and D2 lines, and can be interpreted by transiting circumstellar gas. Through comparison to spectra of nearby stars we exclude the possibility of these features being caused by interstellar absorption. Both the dust mass, which is within an order of magnitude of HD 141569, and the gas mass derived here indicate a late gas dispersal stage of the protoplanetary disc.

Cosmic filaments are thought to host a substantial fraction of the missing baryons at redshifts $z<2$. In this study, we probe the baryonic content of these filaments using localized Fast Radio Bursts (FRBs). Filaments are identified from the galaxy distribution in the Dark Energy Spectroscopic Instrument (DESI) imaging surveys using the DisPerSe algorithm. We find tentative evidence ($\sim 3 \sigma$ significance) for a divergence in the relationship between the dispersion measure (DM) contributed by the intergalactic medium and redshift for FRBs whose signals intersect cosmic filaments compared to those that do not, suggesting excess baryons in the filamentary structures. Assuming an isothermal $\beta$-model gas profile with $\beta=2/3$, this discrepancy is best explained by a central baryon overdensity of $\delta = 20.78$, broadly consistent with previous simulation and observational results. The inferred baryon fraction residing in filaments decreases with redshift, from approximately $0.25-0.30\,\Omega_b$ at $z=0.02$ to $0.15-0.30\,\Omega_b$ at $z=0.5$, and $0.03-0.04\,\Omega_b$ at $z=0.8$. These estimates are likely lower bounds, particularly at $z>0.5$, due to the limited number of identified filaments and localized FRBs at higher redshifts. We also examine various factors that may affect the statistical significance of our results. Our method offers an independent approach to tracing baryons in cosmic filaments and underscores the importance of expanding localized FRB samples and deepening galaxy surveys, which are key steps toward refining these estimates and addressing the missing baryon problem.

We present the spectral characterization of Jupiter Trojans using data from the Gaia mission DR3 spectral catalog (320 spectra of good quality) and from the literature. Gaia spectrophotometry is available in 16 spectrophotometric points covering the 0.33-1.08 $\mu$m range. The final dataset includes 519 Trojans, 291 in the leading swarm (L4) and 228 in the trailing swarm (L5), which we classified using the Bus-Demeo and Mahlke classification schemes. The Trojan population is dominated by featureless asteroids with red spectral slopes belonging to the D-type, and an important fraction ($\sim$ 40\%) belong to the Z-class (in Mahlke taxonomy), characterized by very red slopes. The L4 swarm shows a higher spectral variability and a higher amount, by a factor of 2, of less spectrally red asteroids belonging to the C, P, and X classes, mainly associated with families members. Once excluding peculiar families members the two swarms have very similar average slope and albedo values (7.86 $\pm$ 0.15\% and 7.35 $\pm$0.15\% for the L4 and L5 swarms, respectively). This points toward a common origin of the bodies of the two swarms, likely from the Transneptunian region. We do not observe a spectral color bi-modality distribution, conversely to the results previously reported in the literature. The spectral slope distribution, peaked at 9-11 \%/1000 Å, it is very narrow compared to that of Transneptunian objects (TNOs) or cometary nuclei. Trojans visible spectral slopes match those of the less red TNOs, and lack of extremely red and organic rich bodies. We suggest that the Jupiter Trojans might have been captured by Jupiter from the Centaurs and scattered disk population, and that the lack of extremely red objects among Trojans is due to removal of the organic-rich crust through the sublimation of volatiles and collisions as TNOs migrated inward in the Solar System.

We present the first dataset of collisional (de)-excitation rate coefficients of HCN induced by CO, one of the main perturbing gases in cometary atmospheres. The dataset spans the temperature range of 5-50 K. It includes both state-to-state rate coefficients involving the lowest ten and nine rotational levels of HCN and CO, respectively, and the so-called "thermalized" rate coefficients over the rotational population of CO at each kinetic temperature. The derivation of these coefficients exploited the good performance of the statistical adiabatic channel model (SACM) on top of an accurate interaction potential computed at the CCSD(T)-F12b/CBS level of theory. The reliability of the SACM approach was validated by comparison with full quantum calculations restricted at the lowest total angular momentum of the system. These results provide essential input to accurately model the distribution among the rotational energy levels and the abundance of HCN in cometary atmospheres, accounting for deviations from local thermodynamic equilibrium that typically occurs in such environments.

In this work, we find constraints on the parameter space of the Ricci-Cubic Holographic dark energy (RCHDE) from various observational data sets like Hubble data, cosmic-chronometer data, Baryon-acoustic oscillation data, and also data from gamma-ray bursts. RCHDE is formed from the cubic invariant, which in turn is built from the cubic contractions of the Riemann and Ricci tensors. We have used the Markov chain Monte-Carlo (MCMC) sampling technique to find constraints on the model parameters via Bayesian inference. Contour plots have been obtained for the model parameters, showing their marginalized and joint probability distributions. The best-fit regression lines are found for the constrained model and compared with the standard $\Lambda$CDM model to verify and validate the model. To complement this data analysis mechanism, we have also performed an enhanced machine learning analysis using observational Hubble parameter data. This approach serves to validate the model's predictive power through independent, data-driven regression techniques. Different graphical illustrations of the machine learning techniques have been presented to understand the results. These illustrations reveal a strong agreement between the Hubble parameter predictions from the machine learning models, the theoretical RCHDE model, and observational data.

This study examines the communication of astrobiology and the Search for Life Elsewhere (SLE) in academic papers, press releases, and news articles over three decades. Through a quantitative content analysis, it investigates the prevalence of speculations and promises/expectations in these sources, aiming to understand how research results are portrayed and their potential impact on public perception and future research directions. Findings reveal that speculations and promises/expectations are more frequent in news articles and press releases compared to academic papers. Speculations about conditions for life and the existence of life beyond Earth are common, particularly in news articles covering exoplanet research, while promises of life detection are rare. Press releases tend to emphasize the significance of research findings and the progress of the field. Speculations and promises/expectations in news articles often occur without attribution to scientists and in quotes of authors of the studies, and slightly less so in quotes of outside experts. The study highlights the complex dynamics of science communication in astrobiology, where speculations and promises can generate public excitement and influence research funding, but also risk misrepresenting scientific uncertainty and creating unrealistic expectations. It underscores the need for responsible communication practices that acknowledge the speculative dimension of the field while fostering public engagement and informed decision-making.

During Juno's only flyby of Europa, the Jupiter Energetic Particle Detector Instrument (JEDI) measured complex dropouts in the energetic ion flux in Europa's wake. We investigate the causes of these dropouts, focusing specifically on energetic protons of ~100 keV and ~1 MeV, using back-tracking particle simulations, a prescribed description of Europa's atmosphere and a three-dimensional single fluid magnetohydrodynamics (MHD) model of the plasma-atmosphere interaction. We investigate the role of magnetic field perturbations resulting from the interaction between Jupiter's magnetospheric plasma and Europa's atmosphere and the presence of field-aligned electron beams in Europa's wake. We compare the simulated effect of the perturbed fields on the pitch angle distributions of the ion losses to Juno-JEDI measurements. We find that at ~100 keV, field perturbations are the dominant factor controlling the distribution of the losses along the flyby, while at ~1 MeV a combination of field perturbations and absorption by the surface due to short half bounce periods is required to explain the measured losses. We also find that the effect of charge-exchange with Europa's tenuous atmosphere is weak and absorption by dust in Europa's environment is negligible. Furthermore, we find that the perturbed magnetic fields which best represent the measurements are those that account for the plasma interaction with a sub-/anti-Jovian asymmetric atmosphere, non-uniform ionization of the atmosphere, and electron beams. This sensitivity to the specific field perturbation demonstrates that combining observations and modeling of proton depletions constitute an important tool to probe the electromagnetic field and atmospheric configurations of Europa.

Blue horizontal-branch (BHB) stars are evolved low-mass objects that have completed their core hydrogen burning main-sequence (MS) stage and have lost significant mass during the red giant phase culminating in the helium flash. We determine the fraction of BHBs in binary systems over a wide range of separations in the inner Galactic Halo to constrain mass-loss mechanisms and evolutionary pathways. Using a catalog of 22,336 BHB candidates from Gaia DR3, we analysed radial velocity variations found in spectra (263 spectra of 89 targets) acquired using the Ondrejov Echelle spectrograph at the Astronomical Institute of the Czech Academy of Sciences together with archival spectra from the Ultraviolet and Visual Echelle Spectrograph (UVES) and Fiber-fed Extended Range Optical Spectrograph (FEROS). We searched for wide common proper motion pairs, binary candidates with enhanced astrometric noise, and binaries with astrometric orbital solutions in Gaia DR3. Archival light curves from Gaia DR3 and the Zwicky Transient Facility (ZTF) were checked for binary induced variations. Synthetic SEDs and binary detection probabilities were modelled to account for selection effects. We find a binary fraction of <2.2%, far lower than the rates for their main-sequence (MS) and red-giant branch (RGB) progenitors (30-50%). This suggests that BHBs are either not descendants of binary systems, or that existing companions do not survive the BHB formation process. The negligible binary fraction implies single-star evolution could dominate BHB formation, contrasting with EHB stars where binarity is critical. Our results challenge models of mass loss on the RGB and highlight the need for alternative mechanisms.

A sub-parsec supermassive black hole binary at the center of the galaxy 3C 66B is a promising candidate for continuous gravitational wave searches with Pulsar Timing Arrays. In this work, we search for such a signal in the third data release of the Parkes Pulsar Timing Array. Matching our priors to estimates of binary parameters from electromagnetic observations, we find log Bayes factor of 0.22, highlighting that the source can be neither confirmed nor ruled out. We place upper limits at 95% credibility on the chirp mass $M < 6.05 \times 10^{8}\ M_{\odot}$, and on the characteristic strain amplitude $\textrm{log}_{10}(h_0)< -14.44$. This partially rules out the parameter space suggested by electromagnetic observations of 3C 66B. We also independently reproduce the calculation of the chirp mass with the 3-mm flux monitor data from the unresolved core of 3C 66B. Based on this, we outline a new methodology for constructing a joint likelihood of electromagnetic and gravitational-wave data from supermassive black hole binaries. Finally, we suggest that firmly established SMBHB candidates may be treated as standard sirens for complementary constraints on the universe expansion rate. With this, we obtain constraints on the Hubble constant with 3C 66B.

We investigate the escape of ionizing (Lyman-continuum; LyC) photons from 49 star-forming galaxies at redshifts $\sim 1-1.5$, using far-ultraviolet (FUV) imaging from the Ultra-Violet Imaging Telescope (UVIT) onboard AstroSat. The sample spans a wide range of stellar masses and UV luminosities. LyC emission is undetected in most galaxies (42/49), and stacking these galaxies yields only an upper limit on the observed LyC-to-nonionizing UV flux density ratio ($({F}_{\lambda,\rm{LyC}}/{F}_{\lambda,\rm{UV}})_{\rm{obs}}<0.12$). Including all galaxies (with 7 LyC-leaker candidates) produces a marginal $2.4\sigma$ detection, suggesting that the average LyC signal is driven by a small number of sources. To identify the conditions favorable for LyC escape, we perform stacking analyses in bins of stellar mass, UV slope, compactness, inclination, and star formation rate surface density. A stacked LyC signal is detected at a significance of $\sim3\sigma$ from a subset of galaxies that are characterized as being compact, have high star formation rate surface densities, and blue UV continuum slopes, despite each of these being individually undetected in LyC. This provides the first systematic evidence at $z\sim 1-1.5$ linking these properties to LyC escape, consistent with trends observed in the lower redshift universe. Additionally, LyC leakage appears more efficient in low-mass galaxies ($\log_{10}(M_{*}/M_{\odot})<9.5$), with their average absolute escape fraction ranging from $\langle f_{\text{esc,abs}} \rangle\sim0.1-0.2$ depending on stellar population assumptions. These results support the scenario that compact, low-mass starbursts were key contributors to the ionizing photon budget during cosmic reionization.

The day-side atmospheres of the hottest ultra-hot Jupiters (UHJs) have long been subject to speculation about cloud formation, often without direct observational evidence. Here, we present a detailed analysis of the panchromatic day-side emission spectrum of WASP-121b$\unicode{x2014}$one of the hottest known UHJs$\unicode{x2014}$covering a broad wavelength range of $\sim$0.6-5.1$\unicode{x00B5}$m, based on archival JWST observations from NIRISS and NIRSpec/G396H. We report statistically significant detections of several key molecular species, including H$_2$O (13.4 $\sigma$), CO (14.7 $\sigma$), SiO (4.9 $\sigma$), TiO (5.4 $\sigma$), and VO (6.6 $\sigma$), establishing WASP-121b as one of the most thoroughly characterized exoplanetary atmospheres to date. Additionally, we present the robust detection of Titanate (CaTiO$_3$) clouds at 6.7$\sigma$$\unicode{x2014}$the first such detection in any exoplanet atmosphere. Our analysis further reveals strong evidence of TiO depletion, likely due to sequestration into refractory condensates such as Titanate clouds. The precisely constrained molecular abundances yield a super-solar C/O ratio of 0.963$\pm$0.024, a sub-solar Si/O ratio of 0.034$\pm$0.024, and a metallicity of 4.7$_{-1.38}^{+1.99}$ $\times$solar. These findings offer a unique window into the atmospheric chemistry of an extreme UHJ, positioning WASP-121b as a key benchmark for next-generation models of atmospheric evolution and dynamics.

Though type-Ia supernovae (SNe Ia) are found in all types of galaxies, recent local Hubble constant measurements have disfavored using SNe Ia in early-type or quiescent galaxies, aiming instead for better consistency with SNe Ia in star-forming, late-type host galaxies calibrated by Cepheid distances. Here we investigate the feasibility of a parallel distance ladder using SNe Ia exclusively in quiescent, massive ($\log M_*/M_{\odot} \geq 10$) host galaxies, calibrated by tip of the red giant branch (TRGB) distances. We present TRGB measurements to four galaxies: three measured from the Hubble Space Telescope with the ACS F814W filter, and one measured from the JWST NIRCam F090W filter. Combined with literature measurements, we define a TRGB calibrator sample of five high-mass, early-type galaxies that hosted well-measured SNe Ia: NGC 1316 (SN 2006dd), NGC 1380 (SN 1992A), NGC 1404 (SN 2007on, SN 2011iv), NGC 4457 (SN 2020nvb), and NGC 4636 (SN 2020ue). We jointly standardize these calibrators with a fiducial sample of 124 Hubble-flow SNe Ia from the Zwicky Transient Facility that are matched in host-galaxy and light-curve properties. Our results with this homogenized subsample show a Hubble residual scatter of under 0.11 mag, lower than usually observed in cosmological samples of the full SN~Ia distribution. We obtain a measurement of the Hubble constant, $H_0 = 75.3 \pm 2.9$ km s$^{-1}$ Mpc$^{-1}$, including statistical and estimated systematic uncertainties, and discuss the potential to further improve the precision of this approach. As calibrator and supernova samples grow, we advocate that future cosmological applications of SNe Ia use subsamples matched in host-galaxy and supernova properties across redshift.

Titania doped with tantala is the high index material (high n) for the optical coatings used in LIGO and Virgo and its thermal noise limits LIGO/Virgo observations of astrophysical sources. In this paper, we study temperature induced changes to optical scatter of a multilayer highly reflective coating comprised of silica (low n) and titania doped with germania (high n) as a potential candidate to reduce coating thermal noise in ground-based observatories operating at room temperature. We observe that the scatter measured at 8 degree in a small region is low, with a median starting BRDF of $1.1 \times 10^{-7}\,\mathrm{str}^{-1}$ increasing to $1.2 \times 10^{-6}\,\mathrm{str}^{-1}$ through annealing. The results presented here show the potential of adopting titania doped with germania coatings for future upgrades to LIGO and Virgo and as a pathfinder coating for Cosmic Explorer, a next-generation detector.

We present very high-energy optical photometry and spectroscopic observations of SN 2024bch in the nearby galaxy NGC 3206 (\sim 20 Mpc). We used gamma-ray observations performed with the first Large-Sized Telescope (LST-1) of the Cherenkov Telescope Array Observatory (CTAO) and optical observations with the Liverpool Telescope (LT) combined with data from public repositories to evaluate the general properties of the event and the progenitor star. No significant emission above the LST-1 energy threshold for this observation (\sim 100 GeV) was detected in the direction of SN 2024bch, and we computed an integral upper limit on the photon flux of F_\gamma(>100 GeV) \le 3.61 \times 10^{-12} cm^{-2} s^{-1} based on six nonconsecutive nights of observations with the LST-1, between 16 and 38 days after the explosion. Employing a general model for the gamma-ray flux emission, we found an upper limit on the mass-loss-rate to wind-velocity ratio of \dot M/u_{w} \le 10^{-4} \frac{M_\odot}{yr}\frac{s}{km}, although gamma-gamma absorption could potentially have skewed this estimation, effectively weakening our constraint. From spectro-photometric observations we found progenitor parameters of M_{pr} = 11 - 20 M_\odot and R_{pr} = 531 \pm 125 R_\odot. Finally, using archival images from the Hubble Space Telescope, we constrained the luminosity of the progenitor star to log(L_{pr}/L_\odot) \le 4.82 and its effective temperature to T_{pr} \le 4000 K. Our results suggest that SN 2024bch is a type IIn-L supernova that originated from a progenitor star consistent with a red supergiant. We show how the correct estimation of the mass-loss history of a supernova will play a major role in future multiwavelength observations.

We present a joint analysis of HI absorption Zeeman measurements and the morphology of filamentary HI emission to investigate the three-dimensional structure of the magnetic field in the diffuse neutral interstellar medium (ISM). Our analysis is based on the Arecibo Millennium Survey and new data from the Five-hundred-meter Aperture Spherical radio Telescope (FAST) toward radio sources 3C 75, 3C 207, and 3C 409. Toward 3C 409, we make a 4$\sigma$ Zeeman detection and infer $B_{LOS}$ = 9.1 +/- 1.9$\mu$G, in agreement with Arecibo results. We quantify the dispersion of HI filaments at the locations and velocities of Zeeman components using GALFA-HI narrow-channel emission maps. Focusing on a subsample of 42 spectrally distinct components, we find a weak but statistically significant positive correlation (Spearman r = 0.3, $p = 0.01$) between $|B_{LOS}|$ and the circular variance of HI filament orientation angles. To examine its origin, we characterize the environments probed by HI absorption using dust emission, 3D dust maps, OH absorption, and CO emission. We find evidence that existing HI absorption Zeeman measurements trace magnetic fields that are coherent on parsec scales, probe primarily local gas ($100$-$500$ pc, often at distances consistent with the Local Bubble wall), and exhibit systematic differences in the magnitude of $B_{LOS}$. We attribute the correlation between Zeeman measurements and filamentary HI morphology to large-scale variations in magnetic field strength and/or inclination angle across different Galactic environments, which could arise due to the Local Bubble geometry or enhanced total field strength in denser regions.

In the incremental growth model, planetesimal formation constitutes the least understood step in the process of planetary formation. The two main difficulties in this regard are the collision/fragmentation and the drift barriers. Numerous solutions have been proposed to overcome these barriers, but often need a conjunction of processes to reach the conditions for planetesimal formation. We present numerical simulations, in which the protoplanetary disk turbulence is fully captured rather than modeled with a turbulent diffusion or turbulent viscosity parameter. When the turbulent cascade is taken into account, and in the case of weakly turbulent disks, not only can solid grains be highly concentrated in clusters, but their radial drift can also be slowed or even halted. These results open a unique path to planetesimal formation starting at disk canonical dust-to-gas ratio, namely Keplerian turbulence.

While the warm, ionized gas in the CGM at z>3 is now routinely observed around bright QSOs in Lya emission, little is known about the CGM hot phase due to its expected faintness in the X-ray band, often referred to as the ICM. Here, we report the analysis of 634 ks of Chandra X-ray observations in the MQN01 Cosmic Node, a region containing one of the brightest Lya nebulae and the largest galaxy overdensity discovered so far at z>3. We detect 66 net counts of X-ray emission in the 0.5-2 keV band extending to at least 30 kpc from the brightest QSO in MQN01. The morphology and spectrum are consistent with thermal emission from hot plasma in CIE. Photoionization is negligible, and IC is disfavored. A joint spatial and spectral MCMC analysis provides consistency with a beta-model with a steep density profile and a gas temperature kT~1.8 keV and virial halo mass Mvir~3e13 Mo. The inferred hot gas mass is Mhot(<Rvir)~2.6e12 Mo, which is ~8.3% of Mvir, or ~56% of the theoretical cosmological baryon budget of the halo. The hot gas also emits an exceptionally high Lx, with a measured L2-10~2.3e45 erg/s within the central 30 kpc. This system is a clear outlier in the Lx-Tx plane, indicating a thermodynamic state distinct from that of evolved lower-redshift hot halos. The cooling time in the inner 15-30 kpc is comparable to the local dynamical time, suggesting that the gas could become locally unstable in the absence of heating or feedback. Moreover, the thermal pressure associated with the detected CGM hot phase is large enough to confine the cold and dense clumps, which are required to reproduce the high Lya emission associated with the inner regions of the MQN01 structure. Although limited to a single system, our results provide unique information on the multi-phase properties of the CGM and a view of the nascent thermal hot gas phase observed in local galaxy clusters.

Photonic lanterns (PLs) are fiber-based waveguides that are capable of focal-plane wavefront sensing while simultaneously directing light to downstream science instruments. The optimal choice of wavefront reconstruction algorithm has yet to be determined, and likely depends on the particular observing scenario under consideration. Previous work in simulation suggests that PLs can be used for nonlinear wavefront sensing for several applications, including sensing the low-wind effect and correcting large-amplitude aberrations. We present the design of muirSEAL (miniature IR SEAL), a testbed designed to test PL wavefront reconstruction over Zernike modes and segmented-mirror offsets. We demonstrate throughput and linear wavefront reconstruction at multiple f-numbers. We further present initial laboratory imaging of a new photonic lantern fabricated at Lawrence Livermore National Laboratory.

We have conducted a search for strong gravitational lensing systems in the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys Data Release 10 (DR10). This paper is the fourth in a series of searches (following Huang et al. 2020; Huang et al. 2021; Storfer et al. 2024, Paper I, II, & III respectively). This is the first catalog of lens candidates covering nearly the entirety of the extragalactic sky south of declination $\delta\approx +32$ deg, all of it observed by the DECam, covering $\sim$14,000 $deg^2$. We impose a $z$-band magnitude cut of < 20 in AB magnitude. We deploy a Residual Neural Network and EfficientNet as an ensemble trained on a compilation of known lensing systems and high-grade candidates as well as nonlenses in the same footprint. The predictions from these two base models are aggregated using a meta-learner. After applying our ensemble to the survey data, we exclude known candidates and systems, and use our own visual inspection portal to rank images in the top 0.01 percentile of all neural network recommendations. We have found 811 new lens candidates. These include 484 new candidates in the Legacy Surveys DR9 footprint, all parts of which have been searched for strong lenses at least once before, either by our group or others. Combining the discoveries from this work with those from Paper I (335), II (1210), and III (1512), we have discovered a total of 3868 new candidates in the DESI Legacy Surveys.

The astrophysics of Long GRB (LGRB) progenitors as well as possible cosmological evolution in their properties still poses many open questions. Previous studies suggest that the LGRB rate density (LGRB-RD) follows the cosmic star formation rate density (SFRD) only at high-z and attribute this to the metallicity evolution of progenitor stars. For low z, opinions differ on whether the uptick in the LGRB RD is due to a distinct class of low-luminosity GRBs or perhaps even a different progenitor subclass. To investigate these questions, we utilize data from the Neil Gehrels Swift Observatory and ground-based observatories (redshift). To test the hypothesis that the observations can be mapped (with/without evolution) to the well-established cosmic SFRD, we consider three cases: no evolution, beaming angle evolution, and a simple power-law evolution. The comparison shows that the 'no evolution' case can be ruled out. Our study highlights that the beaming angle evolution or the simple power law evolution are also not sufficient to obtain a good match between the LGRB-RD and SFRD. Rather, the inclusion of multiple evolving properties of LGRBs in combination appears to be required to match the two rate densities in their entirety.

Astrophysical black holes are likely to be surrounded by various forms of matter in the form of disks or halos. While a number of studies have examined the impact of an environment on the lensing of light or gravitational waves from cosmological sources, these have, thus far, been carried out in either a Newtonian or post-Newtonian framework where the environment is superimposed on the black-hole spacetime. By using an exact solution in general relativity describing a black hole embedded within a realistic halo of Hernquist matter distribution, we study deflection angles and image amplification in a fully relativistic setup. It is shown that large ``bumps'' corresponding to the peak of the mass distribution can significantly adjust the inferences made for either the source or lens in various contexts. As an application, we consider ``echoes'' of gravitational waves, sourced by astrophysical lenses rather than being intrinsic to the compact object that produces the signal.

We present a novel theoretical framework for analysing the stability of rotating black hole spacetimes through the conformal properties of the Weyl tensor. By introducing a new conformal invariant constructed from the electric and magnetic parts of the Weyl tensor, we derive a master equation governing perturbations that unifies the Teukolsky and Regge-Wheeler- Zerilli formalisms. Our approach reveals previously unrecognised relationships between quasinormal mode frequencies and the conformal structure of the this http URL prove two fundamental theorems: (i) the conformal stability criterion, which relates mode stability to the sign-definiteness of our conformal invariant, and (ii) the isospectrality theorem for conformally related black hole spacetimes. Numerical calculations for Kerr black holes demonstrate that our formalism predicts new branches in the quasinormal mode spectrum, with frequencies differing from standard predictions by up to 3.7% in the near-extremal regime. These results have significant implications for gravitational wave astronomy and tests of general relativity in the strong-field regime.

Recent gravitational-wave (GW) observations offer a unique opportunity to probe the fundamental nature of compact objects. A growing body of research has focused on exploring the role of dark matter (DM) through the concept of DM-admixed neutron stars (NSs), where the presence of DM can significantly alter key physical properties of NSs, such as their mass, radius, and tidal deformability, ultimately affecting the predicted GW waveform emitted during binary coalescences. In this work, we present a novel observational test that, for the first time, places constraints on the influence of DM inside NSs using real GW data. By reanalyzing signals from events such as GW230529, GW200115, and GW200105, we derive new upper bounds on the DM fraction, $F_{\chi}$, and particle mass, $m_{\chi}$, under the assumption that DM is described by a scalar field with a self-interaction potential. We find that the upper bound on $F_{\chi}$ depends on the specific binary system under analysis, indicating that different DM configurations can be consistent with observations in different ways. In particular, the event GW190814 may be compatible with a DM halo configuration. In contrast, the other events analyzed (GW230529, GW200105 and GW200115) are consistent with DM forming a core inside the NS, yielding strong upper bounds on $F_{\chi}$. The corresponding values for the mass scale $m_{\chi}$ are also discussed in the text. This work offers a new approach to probing DM in the context of compact NS objects through GW observations.

Stratospheric High Altitude Balloons (HABs) have great potential as a remote sensing platform for Earth Observations that complements orbiting satellites and low flying drones. At altitudes between 20-35 kms, HABs operate significantly closer to ground than orbiting satellites, but significantly higher than most drones. HABs therefore offer a unique potential to deliver high spatial resolution imaging with large area coverage. Another two imaging parameters that are important for Earth Observation applications are spectral resolution and spectral range. In this paper, we therefore present the development and testing of a hyperspectral imaging system, able to record near-video-rate images in narrow contiguous spectral bands, from a HAB platform. In particular, we present the first stratospheric environmental tests and HAB flight of a snapshot hyperspectral camera, based on Computed Tomography Imaging Spectroscopy (CTIS), which is well suited to cope with the challenges posed by the motion of the HAB platform and the stratospheric environment. We have successfully acquired images with the system under both simulated stratospheric conditions in the Mars Simulation Laboratory at Aarhus University and during a 5 hour HAB flight mission named HEIMDAL from Kiruna in October 2024 as part of the REXUS/BEXUS 34/35 2024 campaign organized by DLR-SNSA. The study represents a step towards deploying the HAB platform for high quality land cover classification.

We present new insights into the phenomenology of post-Newtonian spin precession in black-hole binaries. Using multi-timescale methods, previous work has shown that the precession and nutation dynamics in such systems can be classified into so-called spin morphologies --mutually exclusive regions that partition the configuration space and characterize the motion of the black-hole spins relative to the binary's angular momentum. Radiation reaction can induce secular transitions between different morphology classes, which are generic occurrences during the inspiral of black-hole binaries. In this contribution, we systematically explore a more restrictive class of solutions in which multiple morphological transitions occur concurrently, i.e., within the same precession cycle. We find that all such cases can be mapped and characterized analytically, and we confirm these findings through numerical integrations. These coincident transitions correspond to extreme spin configurations in black-hole binaries with potential observational signatures in gravitational-wave astronomy.

The detection and classification of anomalies in gravitational wave data plays a critical role in improving the sensitivity of searches for signals of astrophysical origins. We present ABNORMAL (AI Based Nonstationarity Observer for Resectioning and Marking AnomaLies), a deep neural network (DNN) model for anomaly detection that is trained exclusively on simulated Gaussian noise. By removing dependence on real data for training, the method resolves a circular paradox in anomaly detection: training on real data implicitly involves prior segregation of stationary from non-stationary data but this is not possible unless all anomalies are detected first. ABNORMAL is an autoencoder-based DNN, commonly used in anomaly detection, with the key innovation that it is trained to predict statistical features of noise rather than reconstructing the noise time series themselves. The statistical features are obtained by applying Gabor and Wavelet filter banks, which implement time-frequency analysis, and are subsequently combined through multi-view fusion using a dual-path architecture. We quantify the performance of our method on simulated and real LIGO data. Application to data from the O1 to O3b observational runs uncovers a rich landscape of anomalies over different timescales, including many that do not fit within known classes.

We analyze the observational features of hot-spots orbiting parametrized black hole (BH) spacetimes. We select a total of four BH spacetimes, two from the Johanssen-Psaltis (JP) parametrization, and two from the Konoplya-Zhidenko (KZ) parametrization, corresponding to the most extreme configurations whose shadow sizes are within the $2\sigma$-constraints of the Event Horizon Telescope (EHT). We use the ray-tracing software GYOTO to simulate the orbit of a spherically symmetric hot-spot emitting synchrotron radiation close to a central parametrized BH object, in a vertical magnetic field configuration, and we extract the corresponding astrometric and polarimetric observables for the Stokes parameters I, Q and U, namely the time integrated fluxes, temporal fluxes and magnitudes, temporal centroid, temporal QU-loops, and temporal Electric Field Position Angle (EVPA). Our results indicate that at low inclination the astrometric observables extracted from the parametrized BH spacetimes considered are qualitatively similar to those extracted from the Schwarzschild one, with minor quantitative deviations caused by differences in the size and position of the secondary images. On the other hand, the polarimetric observables at high inclination present qualitative differences, but these are only visible for a short portion of the whole hot-spot orbit. Furthermore, the observables extracted from the JP parametrized BH models deviate more prominently from those of the Schwarzschild model than the ones extracted from the KZ parametrized BH models, with the JP model with a positive free parameter deviating the most among all models tested. Given the strong similarity among the observables extracted from all models tested, we point out that more precise observations are needed to successfully impose constraints on parametrized BH models via this method.

Gravitational redshift is an important prediction of General Relativity (GR). We propose a novel dual one-way frequency comparison scheme for gravitational redshift tests in satellite-ground clock experiments. Unlike conventional triple-link Doppler cancellation methods, our approach requires only two optical links, an uplink and a downlink, with Doppler data recorded independently at the satellite and ground station, respectively. Crucially, coherent combination of these measurements suppresses first-order Doppler frequency shift to second order while doubling the gravitational redshift signal. This scheme reduced engineering complexity (50$\%$ fewer links) and amplified redshift sensitivity with a factor of 2 for satellite-ground clock comparison. The scheme's viability hinges on high-precision clock synchronization ($\leq$ 0.1 ns for $10^{-18}$ clock comparison), which is potentially achievable via state-of-the-art time-transfer techniques. Systematic error analysis confirms residual Doppler, atmospheric, and Shapiro shifts remain below $10^{-17}$ for optimized orbits, establishing this scheme as a useful tool for probing fundamental physics beyond GR.

We investigate scalar-induced gravitational waves (SIGWs) in the framework of spatially covariant gravity (SCG), a broad class of Lorentz-violating modified gravity theories respecting only spatial diffeomorphism invariance. Extending earlier SCG formulations, we compute the general kernel function for SIGWs on a flat Friedmann-Lemaître-Robertson-Walker background, focusing on polynomial-type SCG Lagrangians up to $d=3$, where $d$ denotes the total number of derivatives in each monomial. We derive explicit expressions for the kernel in the case of power-law time evolution of the coefficients, and restrict attention to the subset of SCG operators whose tensor modes propagate at the speed of light, thereby avoiding late-time divergences in the fractional energy density of SIGWs. Instead of the usual Newtonian gauge, the breaking of time reparametrization symmetry in SCG necessitates a unitary gauge analysis. We compute the energy density of SIGWs for representative parameter combinations, finding distinctive deviations from general relativity (GR), including scale-dependent modifications to both the amplitude and the spectral shape. Our results highlight the potential of stochastic GW background measurements to probe spatially covariant gravity and other Lorentz-violating extensions of GR.

Magnetic reconnection can power spectacular high-energy astrophysical phenomena by producing non-thermal energy distributions in highly magnetized regions around compact objects. By means of two-dimensional fully kinetic particle-in-cell (PIC) simulations we investigate relativistic collisionless plasmoid-mediated reconnection in magnetically dominated pair plasmas with and without guide field. In X-points, where diverging flows result in a non-diagonal thermal pressure tensor, a finite residence time for particles gives rise to a localized collisionless effective resistivity. Here, for the first time for relativistic reconnection in a fully developed plasmoid chain we identify the mechanisms driving the non-ideal electric field using a full Ohm's law by means of a statistical analysis based on our PIC simulations. We show that the non-ideal electric field is predominantly driven by gradients of nongyrotropic thermal pressures. We propose a kinetic physics motivated non-uniform effective resistivity model, which is negligible on global scales and becomes significant only locally in X-points, that captures the properties of collisionless reconnection with the aim of mimicking its essentials in non-ideal magnetohydrodynamic descriptions. This effective resistivity model provides a viable opportunity to design physically grounded global models for reconnection-powered high-energy emission.

Small-scale, newly emerging internetwork (IN) magnetic fields are considered a viable source of energy and mass for the solar chromosphere and possibly the corona. Multiple studies show that single events of flux emergence can indeed locally heat the low solar atmosphere through interactions of the upward propagating magnetic loops and the preexisting ambient field lines. However, the global impact of the newly emerging IN fields on the solar atmosphere is still unknown. In this paper, we study the spatio-temporal evolution of IN bipolar flux features and analyze their impact on the energetics and dynamics of the quiet Sun atmosphere. We use high resolution, multi-wavelength, coordinated observations obtained with the Interface Region Imaging Spectrograph (IRIS), Hinode, and the Solar Dynamics Observatory (SDO) to identify emerging IN magnetic fields and follow their evolution. Our observational results suggest that only the largest IN bipoles are capable of heating locally the low solar atmosphere, while the global contribution of these bipoles appears to be marginal. However, the total number of bipoles detected and their impact estimated in this work is limited by the sensitivity level, spatial resolution, and duration of our observations. To detect smaller and weaker IN fields that would maintain the basal flux, and examine their contribution to the chromospheric heating, we will need higher resolution, higher sensitivity and longer time series obtained with current and next-generation ground- and space-based telescopes.

Neutrinos are believed to be the most abundant fermions in the Universe, but their masses are unknown, except for being non-zero but much smaller than other fermions. Cosmological relic neutrinos could also have non-zero chemical potentials (or asymmetries). Using neutrino-involved N-body simulations, we investigate the neutrino effects on the matter pairwise velocity, which itself is an interesting probe of cosmology. We find that for light-halo ($[10^{11},10^{13}]\ M_\odot$) mean pairwise velocity, in the transition range ($[4,15]\ \mathrm{Mpc}$), the effects of neutrino masses overwhelm the effects of neutrino asymmetries, while in the two-halo-group range ($[25,50]\ \mathrm{Mpc}$), for both light and heavy haloes ($[10^{13},10^{15}]\ M_\odot$), the effects of neutrino asymmetries dominate, making it possible to disentangle the two effects. We provide fitting formulae to quantify the effects of neutrino mass and asymmetry on halo-halo pairwise velocities.

Linking atmospheric measurements to the bulk planetary composition and ultimately the planetary origin is a key objective in planetary science. In this work, we identify the cases in which the atmospheric composition represents the bulk composition. We simulate the evolution of giant planets considering a wide range of planetary masses ($0.3-2~M_{\mathrm{J}}$), initial entropies ($8-11~k_{\mathrm{B}} m_\mathrm{u}^{-1}$), and primordial heavy-element profiles. We find that convective mixing is most efficient at early times (ages $\lesssim 10^7$ yr) and that primordial composition gradients can be eroded. In several cases, however, the atmospheric composition can differ widely from the planetary bulk composition, with the exact outcome depending on the details. We show that the efficiency of convection is primarily controlled by the underlying entropy profile: For low primordial entropies of $8-9~k_\mathrm{B} m_\mathrm{u}^{-1}$, convective mixing can be inhibited and composition gradients can persist over billions of years. The scaling of mixing efficiency with mass is governed by the primordial entropy. For the same primordial entropy, low-mass planets mix more efficiently than high-mass planets. If the primordial internal entropy would increase with mass, however, this trend could reverse. We also present a new analytical model that predicts convective mixing under the existence of composition (and entropy) gradients. Our results emphasize the complexity in the interpretation of atmospheric abundance measurements and show the great need to better understand the planetary formation process as it plays a key role in determining the planetary evolution and final structure.

Stellar haloes form through the disruption of satellite galaxies over time, making them a promising observable for constraining galaxy merger histories. We use a dynamical decomposition technique to isolate the stellar haloes of Milky Way-mass galaxies in the $100\,{\rm Mpc}$ EAGLE simulation and study their relationship to the merger histories of their hosts. We define the stellar halo as the stellar mass that is bound to the central subhalo but not associated with the disc or bulge components of a galaxy, and we quantify their merger histories using the most significant merger since $z=1$. Surprisingly, we find that the fraction of a galaxy's total stellar mass in the stellar halo, $f_{\rm SH}$, is not a reliable indicator of its merger activity. Contrary to common assumptions, disc galaxies with low $f_{\rm SH}$ do not necessarily have quiescent merger histories. In fact, roughly one quarter experienced a merger at $z \leq 1$ with a satellite whose stellar mass was at least 10 per cent of the host galaxy's stellar mass. These galaxies undergo mergers with satellites on circular orbits that are roughly co-planar with the pre-existing disc and thereby avoid contributing mass to the stellar halo. Instead, such mergers build thick, extended discs and supply fresh gas that often triggers a significant episode of star formation in the disc. Our results suggest that disc galaxies with low-mass stellar haloes, such as the Milky Way, can have varied and active merger histories, and that stellar haloes may not be a reliable tool for inferring galaxy merger histories.

We present the first spatially-resolved infrared images of supernova remnants (SNRs) in M33 with the unprecedented sensitivity and resolution of JWST. We analyze 40 SNRs in four JWST fields: two covering central and southern M33 with separate NIRCam (F335M, F444W) and MIRI (F560W, F2100W) observations, one $\sim$5 kpc-long radial strip observed with MIRI F770W, and one covering the giant HII region NGC 604 with multiple NIRCam and MIRI broad/narrowband filters. Of the 21 SNRs in the MIRI (F560W+F2100W) field, we found three clear detections (i.e., identical infrared and H$\alpha$ morphologies), and six partial-detections, implying a detection fraction of 43\% in these bands. One of the SNRs in this field, L10-080, is a potential candidate for having freshly-formed ejecta dust, based on its size and centrally-concentrated 21 $\mu$m emission. In contrast, only one SNR (out of 16) is detectable in the NIRCam F335M+F444W field. Two SNRs near NGC 604 have strong evidence of molecular (H$_2$) emission at 4.7 $\mu$m, making them the farthest known SNRs with visible molecular shocks. Five SNRs have F770W observations, with the smaller younger objects showing tentative signs of emission, while the older, larger ones have voids. Multi-wavelength data indicate that the clearly-detected SNRs are also among the smallest, brightest at other wavelengths (H$\alpha$, radio and X-ray), have the broadest line widths (H$\alpha$ FWHM$\sim$250-350 km/s), and the densest environments. No correlation between the JWST-detectability and local star-formation history of the SNRs is apparent.

The origin of Lyman Continuum (LyC) photons responsible for reionizing the universe remains a mystery, with the fraction of escaping LyC photons from galaxies at z$\sim$ 6 to 12 being highly uncertain. While direct detection of LyC photons from this epoch is hindered by absorption from the intergalactic medium, lower redshift analogs offer a promising avenue to study LyC leakage. We present Hubble Space Telescope Cosmic Origins Spectrograph (HST COS) observations of five low redshift (z$\sim$ 0.3) massive starburst galaxies, selected for their high stellar mass and weak [SII] nebular emission - an indirect tracer of LyC escape. Three of the five galaxies show LyC leakage, highlighting the reliability of weak [SII] as a tracer, especially in light of recent JWST discoveries of z $>$ 5 galaxies with similarly weak [SII] emission. The dust corrected LyC escape fractions, which represent the LyC photons that would escape in the absence of dust, range from 33% to 84%. However, the absolute escape fractions, which show the LyC photons escaping after passing through both neutral hydrogen absorption and dust attenuation, are significantly lower, ranging between 1% and 3%. This suggests that while the galaxies are nearly optically thin to HI, their high dust content significantly suppresses LyC photon escape. These [SII] weak, massive leakers are distinct from typical low-redshift LyC emitters, showing higher metallicity, lower ionization states, more dust extinction and higher star formation surface densities. This suggests that these galaxies constitute a distinct population, likely governed by a different mechanism facilitating LyC photon escape. We propose that the feedback-driven winds in these compact starbursts create ionized channels through which LyC photons escape, aligning with a picket-fence model.

Two notable anomalies in radio observations -- the excess radiation in the Rayleigh-Jeans tail of the cosmic microwave background, revealed by ARCADE2, and the twice-deeper absorption trough of the global 21cm line, identified by EDGES -- remain unresolved. These phenomena may have a shared origin, as the enhancement of the 21cm absorption trough could arise from excess heating. We investigate this scenario through the framework of axion-like particles (ALPs), showing that the resonant conversion of ALPs into photons can produce a photon abundance sufficient to resolve both anomalies simultaneously. Our model naturally explains the observed radio excess between 0.4 and 10GHz while also enhances the 21cm absorption feature at 78MHz. Furthermore, it predicts a novel power-law scaling of the radio spectrum above 0.5GHz and an additional absorption trough below 30MHz, which could be verified through cross-detection in upcoming experiments.

Fast radio bursts (FRBs), especially those from repeating sources, exhibit a rich variety of morphologies in their dynamic spectra (or waterfalls). Characterizing these morphologies and spectro-temporal properties is a key strategy in investigating the underlying unknown emission mechanism of FRBs. This type of analysis has been typically accomplished using two-dimensional Gaussian techniques and the autocorrelation function (ACF) of the waterfall. These techniques are effective and precise at all duration scales, but can be limited in the presence of scattered tails, complex morphologies, or recently observed microshot forests. Here, we present a technique that involves the tagging of per-channel arrival times of an FRB to perform spectro-temporal measurements using a Gaussian profile model for each channel. While scattering and dispersion remain important and often dominating sources of uncertainty in measurements, this technique provides an adaptable and firm foundation for obtaining spectro-temporal properties from all types of FRB morphologies. We present measurements using this technique of several hundred bursts across 12 repeating sources, including over 400 bursts from the repeating sources FRB 20121102A, FRB 20220912A, and FRB 20200120E, all of which exhibit recently observed microsecond-long ultra-FRBs, as well as 143 drift rates. In addition to retrieving the known relationship between sub-burst slope and duration, we explore other correlations between burst properties. We find that the sub-burst slope law extends smoothly to ultra-FRBs, and that ultra-FRBs appear to form a distinct population in the duration-frequency relation.

We explore transformations of the Friedman-Lemaître-Robertson-Walker (FLRW) metric and cosmological parameters that align with observational data, aiming to gain insights into potential extensions of standard cosmological models. We modify the FLRW metric by introducing a scaling factor, $e^{2\Theta(a)}$ (the cosmological scaling function, CSF), which alters the standard relationship between cosmological redshift and the cosmic scale factor without affecting angular measurements or Cosmic Microwave Background (CMB) anisotropies. Using data from DESI Year 1, Pantheon+ supernovae, and the Planck CMB temperature power spectrum, we constrain both the CSF and cosmological parameters through a Markov Chain Monte Carlo approach. Our results indicate that the CSF model fits observational data with a lower Hubble constant (although it is compatible with the value given by Planck 2018 within 1$\sigma$) and is predominantly dark-matter-dominated. Additionally, the CSF model produces temperature and lensing power spectra similar to those predicted by the standard model, though with lower values in the CSF model at large scales. We have also checked that when fitting a CSF model without dark energy to the data, we obtain a more negative conformal function. This suggests that the CSF model may offer hints about missing elements and opens up a new avenue for exploring physical interpretations of cosmic acceleration.

We present a numerical code that solves the Misner-Sharp system for a spherically symmetric cosmological model containing both a scalar field and a perfect fluid. While the code is capable of exploring general scenarios involving an minimally coupled scalar field and perfect fluid, we focus on the regime where the scalar field dominates the dynamics, particularly in the post-inflationary scalar field-dominated scenario, where the universe is governed by a rapidly oscillating scalar field for a period lasting a few $e$-folds. We analyse the threshold for PBH formation under quadratic and quartic potentials, evolving configurations from superhorizon scales. Our results confirm that a quartic potential behavior is similar to the radiation-dominated universe, resulting in a PBH formation threshold close to the well-established value in radiation backgrounds. Conversely, in the quadratic case, we observe a significant deviation from the expected dust-like behaviour, due to wave-like effects opposing the gravitational collapse. While numerical limitations prevent us from evolving a wide range of initial conditions to determine a precise threshold for PBH formation, our findings suggest that PBH formation may be suppressed with respect to the pure dust scenario, allowing the formation of stable solitonic structures instead. This study highlights the importance of properly accounting for wave dynamics in oscillating scalar fields when characterising PBH formation.

We aim to develop a new method to infer the sub-beam probability density function (PDF) of H2 column densities and the dense gas mass within molecular clouds using spatially unresolved observations of molecular emission lines in the 3 mm band. We model spatially unresolved line integrated intensity measurements as the average of an emission function weighted by the sub-beam column density PDF. The emission function, which expresses the line integrated intensity as a function of the gas column density, is an empirical fit to high resolution (< 0.05 pc) multi-line observations of the Orion B molecular cloud. The column density PDF is assumed to be parametric, composed of a lognormal distribution at moderate column densities and a power law distribution at higher column densities. To estimate the sub-beam column density PDF, the emission model is combined with a Bayesian inversion algorithm (the Beetroots code), which takes account of thermal noise and calibration errors. We validate our method by demonstrating that it recovers the true column density PDF of the Orion B cloud, reproducing the observed emission line integrated intensities. We apply the method to 12CO(J=1-0), 13CO(J=1-0), C18O(J=1-0), HCN(J=1-0), HCO+(J=1-0) and N2H+(J=1-0) observations of a 700 x 700 pc2 field of view (FoV) in the nearby galaxy M51. On average, the model reproduces the observed intensities within 30%. The column density PDFs obtained for the spiral arm region within our test FoV are dominated by a power-law tail at high column densities, with slopes that are consistent with gravitational collapse. Outside the spiral arm, the column density PDFs are predominantly lognormal, consistent with supersonic isothermal turbulence. We calculate the mass associated with the powerlaw tail of the column density PDFs and observe a strong, linear correlation between this mass and the 24$\mu$m surface brightness.

Self-interacting dark matter (SIDM) theories predict that dark matter halos experience core-collapse in late-stage evolution, a process where the halo's inner region rapidly increases in density and decreases in size. This process can be modeled by treating the dark matter as a gravothermal fluid, and solving the fluid equations to predict the density profile evolution. This model is incomplete without calibration to N-body simulations, through a constant factor $\beta$ included in the thermal conductivity for the long-mean-free-path limit. The value of $\beta$ employed in the gravothermal fluid formalism has varied between studies, with no clear universal value in the literature. In this work, we use the N-body code Arepo to conduct a series of isolated core-collapse simulations across a range of scattering cross-sections, halo concentrations, and halo masses to calibrate the heat transfer parameter $\beta$. We find that $\beta$ is independent of cross-section, halo concentration, and halo mass for velocity independent elastic scattering cross-sections. We present a model for an effective $\beta$ as a function of a dimensionless cross-section, to describe halo evolution in the long mean free path limit, and show that it accurately captures halo evolution as long as the cross section is not too large. This effective model facilitates comparisons between simulations and the gravothermal model, and enables fast predictions of the dark matter density profile at any given time without running N-body simulations.

The Zwicky Transient Facility (ZTF), with its extensive optical monitoring capabilities, has provided an unprecedented opportunity to study the long-term variability of active galactic nuclei (AGNs). In this work, we present a comparative analysis of optical colour and brightness variability for two $\mathrm{log(L_{Bol})}-z$ matched samples, consisting of 2095 Narrow-line Seyfert 1 (NLSy1) galaxies and a control sample consisting of 2380 Broad-line Seyfert 1 (BLSy1) galaxies. Using over six years of r-band and g-band light curves from the ZTF Data Release 22 (DR22), we characterize flux variability, fractional flux variability, and amplitude of temporal variability for each source in the samples. Our results indicate that BLSy1 galaxies exhibit significantly stronger variability compared to NLSy1s. To probe colour variability, we utilize quasi-simultaneous light curves, with half-hour epoch differences between $\mathrm{g}$- and $\mathrm{r}$-band measurements. We evaluated the colour index using both flux-flux space analysis and linear regression in magnitude-magnitude space. We find that large majorities of these sources -- 74% of NLSy1 and 79% of BLSy1 -- exhibit a clear ``bluer-when-brighter'' (BWB) trend, although part of this effect may arise from contamination by the non-varying, predominantly redder flux of the host galaxy. Furthermore, rest-frame structure function analysis reveals that BLSy1 galaxies are $1.44 \pm 0.06$ times more variable than NLSy1s. These results can provide valuable insights into the variability properties of AGN subclasses and their underlying physical drivers.

We present ultraviolet (UV) through near-infrared (NIR) photometric and spectroscopic observations of the nearby SN 2024pxl, the third Type Ia supernova (SN Ia) in NGC 6384. SN 2024pxl is a Type Iax supernova (SN Iax) with an intermediate luminosity ($M_r = -16.99\pm0.32$ mag) and an average SN Iax light curve decline rate. SN 2024pxl was discovered $\sim$3 days after first light, and the rising light curve follows a single power law that is inconsistent with significant interaction with a companion star or circumstellar material. Our extensive NIR photometric coverage is comparable to that of the well-observed SNe Iax 2005hk and 2012Z, and we demonstrate that the $J-H$ colors of SNe Iax differ from normal SNe Ia and appear to be more homogeneous as a class. Spectroscopically, we report the earliest-ever NIR spectrum of a SN Iax as measured from maximum light ($t\approx-9$ days): a featureless continuum with similarities to a $\sim$9,000 K blackbody, and the line velocities are consistent with a mixed-ejecta structure, with C, Si, and Fe having similar velocities and velocity evolutions. We find a tentative correlation between the $H$-band break Co II velocity $\sim$20 days post-peak and absolute magnitude, with more luminous SNe Iax showing faster Co II velocities. Our observations suggest that SN 2024pxl resulted from the thermonuclear disruption of a CO white dwarf star that undergoes deflagration burning.

Stellar atmospheric element abundance ratios of stars retain information about their birth conditions, helping elucidate their origin and nature. In this letter, we analyse and contrast the hydrostatic and explosive $\alpha$-element abundance ratios, and the ratio of the two (the hex ratio), for a large sample of Galactic globular clusters (GCs), halo substructures, satellite galaxies, and the Milky Way high-/low-$\alpha$ discs using data from the $APOGEE$ survey. Our results show that: $i$) Milky Way GCs and halo substructures appear to have qualitatively similar hex ratios across a broad range of [Fe/H], that are higher than that of dwarf satellite galaxies of similar [Fe/H]; $ii$) for all stellar populations studied, there is a trend in the hex ratio with [Fe/H]; $iii$) there is a weak trend in the hex ratio with respect to age for Galactic GCs, but not with initial or final GC mass; $iv$) there are no differences in the hex ratio between GCs formed $in$ $situ$ versus those labelled as accreted.

It has long been suggested that a fraction of the dark matter in the Universe could exist in the form of primordial black holes (PBHs) that have existed since the radiation dominated era. Recent studies have suggested that these PBHs may be the progenitors to the population of high-redshift, supermassive black holes (SMBHs) observed since the launch of JWST. For the first time, we have included PBHs in cosmological simulations, to test whether PBHs can sink to the center of collapsing halos, locate dense gaseous regions and experience significant growth. We tested PBH-to-DM mass ratios of $f_{\rm PBH}$ = $10^{-4}$ and $10^{-3}$, with an initial PBH mass of 1000 M$_\odot$, as inspired by recent observational constraints. We find that at $f_{\rm PBH} = 10^{-3}$, a number of PBHs were able to embed themselves in dense gas and grow to $10^{4}$-$10^{5}$ M$_\odot$ by $z=20$. These intermediate black holes (IMBHs) are possible progenitors to the highest redshift SMBH observations such as GNZ-11 ($10^{6}$ M$_\odot$ by $z=10$), outperforming light seed black hole (BH) growth seen in recent simulations without the need to invoke heavy seeding prescriptions. On the other hand, $f_{\rm PBH} = 10^{-4}$ resulted in no significant BH growth, emphasizing that the ability of PBHs to act as SMBH seeds is sensitive to the true value of $f_{\rm PBH}$ in the Universe, and showing that the $f_{\rm PBH} =10^{-4}-10^{-3}$ boundary marks the threshold above which SMBH seeding via 1000 M$_\odot$ PBHs becomes effective. This is the first step towards building a realistic PBH framework in cosmological simulations.

Radiative processes play a pivotal role in shaping the thermal and chemical states of gas across diverse astrophysical environments, from the interstellar medium (ISM) to the intergalactic medium. We present a hybrid cooling model for cosmological simulations that incorporates a comprehensive treatment of radiative processes, including parameterizations of the interstellar radiation field, cosmic ray rates, and dust physics. The model uses the chimes chemical network and combines on-the-fly non-equilibrium calculations with quasi-equilibrium cooling rates. The quasi-equilibrium rates account for the time-dependent free electron fractions of elements tracked in non-equilibrium, balancing computational efficiency with physical accuracy. We evaluate the performance under various conditions, including the thermal evolution of primordial gas at the cosmic mean density, the properties of the warm and cold neutral media in Milky Way-like galaxies, and the atomic-to-molecular hydrogen transition. We demonstrate that thermal equilibrium predictions for the neutral phases of the ISM underestimate the median gas pressures in simulations of isolated galaxies by up to 0.5 dex. Finally we find that the atomic-to-molecular hydrogen transition is shifted to lower densities by up to 1 dex if oxygen is not included in the chemical network. Our work provides a robust framework for studying the multi-phase ISM and its role in galaxy formation and evolution.

We investigate the consequences of non-ideal mixing between silicate, iron metal, and hydrogen for the structures of the cores of sub-Neptunes with implications for super-Earths, warm Neptunes, and ice giants. A method of extrapolating what we know about the miscibility in the three bounding binary systems MgSiO$_3$-H$_2$, MgSiO$_3$-Fe, and Fe-H$_2$ to the ternary composition space is used to deduce the phase equilibria of this system at relevant temperature and pressure conditions. We find that while separate silicate and metal phases can exist at shallow depths, the phases become entirely miscible deeper in the cores, thus altering the density structure of the cores. The assumption that the interiors of large rocky planets, either with extant magma oceans beneath H$_2$-rich envelopes, or evolved from such bodies, are composed of a differentiated metal core overlain by a silicate mantle is inconsistent with our understanding of the phase equilibria of these bodies.

The discovery of the massive black hole (BH) system Gaia BH3 in pre-release Gaia DR4 data suggests that wide BH binaries with luminous companions may be significantly overrepresented at low metallicities. Motivated by this finding, we have initiated a spectroscopic survey of low-metallicity stars exhibiting elevated RUWE values in Gaia DR3, using the FEROS and APF spectrographs. We identify promising BH binary candidates as objects with instantaneously measured radial velocities (RVs) that are very different from their mean RVs reported in Gaia DR3. Thus far, we have observed over 500 targets, including a nearly complete sample of stars with $\text{[Fe/H]} < -1.5$, RUWE $> 2$, and $G < 15$. Our search has yielded one promising target exhibiting slow acceleration and an RV more than 98 km s$^{-1}$ different from its DR3 mean RV, as well as dozens of other candidates with smaller RV discrepancies. We quantify the sensitivity of our search using simulations, demonstrating that it recovers at least half of the BH companions within our selection criteria. We make all the spectra and RVs from our survey publicly available and encourage further follow-up.

We evaluate the reliability of CNEOS-derived ephemerides of fireball events given the absence of the underlying data. We analyzed 18 events that have both (i) sufficient satellite information to derive orbits and (ii) ground-based observational counterparts. We quantify the uncertainties on these calibrated events using the orbital similarity criterion $D_D$. We also examine the velocity components imbalance and identify discriminants that can indicate the accuracy of an event. We identify two groups in the CNEOS database. CNEOS data produces ephemeris determinations with $D_D$<0.1 for fireballs reported either (i) after late 2017 or (ii) with impact energies above 0.45 kt with 74-78% of events having $D_D$=0.03$\pm$0.02, while ~11% show $D_D$<0.008. Our statistical test confirms these two parameters as the only reliable discriminants that, when combined, explain the two accuracy groups. Daylight, z-velocity component, low altitude, long duration, and latitude might also indicate errors, although the limited dataset may obscure correlations. No clear discriminants are identified for more restrictive $D_D$ cut-offs. We provide estimates of orbital uncertainties for calibrated events. The hyperbolic fireball subset in the CNEOS database appears as an outlier in the velocity imbalance test. Our results confirm that the fidelity of CNEOS fireball data improved significantly from 2018, likely due to the deployment of next-generation space sensors, and show a growing number of high-velocity events. Hyperbolic candidates should be interpreted with caution, as their velocities and inclinations likely reflect measurement errors. Accuracy constraints remain limited by the dataset size, as evidenced by the lack of statistically significant dependence on duration, preventing strong conclusions from being drawn.

We present $\sim$500 ks XRISM observations covering the central and two northern regions of the Abell 2029 galaxy cluster. Resolve enables us to distinguish multiple emission lines from hydrogen-like and helium-like iron (Fe) ions. This study focuses on the multi-temperature structure of Abell 2029 using line-ratio diagnostics. Using a single-temperature collisionally ionized equilibrium model, we measure average plasma temperatures of 6.73 keV, 7.61 keV, and 8.14 keV in the central, inner northern, and outer northern regions, respectively, spanning a radial range up to 700 kpc. To further investigate thermal structure, we derive excitation and ionization temperatures by comparing observed emission-line flux ratios with atomic database predictions. Significant deviations from the single-temperature CIE model in the central and inner northern regions indicate the presence of multi-phase gas. The excitation and ionization temperatures range from 2.85 keV to 8.5 keV in the central region, 4.3 keV to 9.8 keV in the inner northern region, and 8.3 keV to 10.4 keV in the outer northern region. These temperature distributions are largely consistent with the previously observed temperature gradient of A2029. However, Resolve detects two notably cooler components--3.42 keV in the central region and $\sim$4.3 keV in the inner northern region--likely associated with displaced cool gas due to gas sloshing. Additionally, we thermally resolve a 2.85 keV gas component at the core of A2029--potentially a significant development in our understanding of gas cooling. We propose that this cooler gas is a direct product of ongoing cooling processes in A2029, having already cooled to its present temperature. If this temperature structure is stable and no heating mechanism is present, this reservoir is likely to cool to even lower temperatures and form stars.

We describe the second data release (DR2) of the FIRE-2 cosmological zoom-in simulations of galaxy formation, from the Feedback In Realistic Environments (FIRE) project, available at this http URL. DR2 includes all snapshots for most simulations, starting at z ~ 99, with all snapshot time spacings <~ 25 Myr. The Core suite -- comprising 14 Milky Way-mass galaxies, 5 SMC/LMC-mass galaxies, and 4 lower-mass galaxies -- includes 601 snapshots to z = 0. For the Core suite, we also release resimulations with physics variations: (1) dark-matter-only versions; (2) a modified ultraviolet background with later reionization at z = 7.8; (3) magnetohydrodynamics, anisotropic conduction, and viscosity in gas; and (4) a model for cosmic-ray injection, transport, and feedback (assuming a constant diffusion coefficient). The Massive Halo suite now includes 8 massive galaxies with 278 snapshots to z = 1. The High Redshift suite includes 34 simulations: in addition to the 22 simulations run to z = 5, we now include 12 additional simulations run to z = 7 and z = 9. We also release 4 dark-matter-only cosmological boxes used to generate zoom-in initial conditions for many FIRE simulations. Most simulations include catalogs of (sub)halos and galaxies at all available snapshots, and most Core simulations to z = 0 include full halo merger trees.

Strong gravitational lensing acts as a cosmic telescope, enabling the study of the high-redshift universe. Astronomical interferometers, such as the Atacama Large Millimeter/submillimeter Array (ALMA), have provided high-resolution images of strongly lensed sources at millimeter and submillimeter wavelengths. To model the mass and light distributions of lensing and source galaxies from strongly lensed images, strong lens modeling for interferometric observations is conventionally performed in the visibility space, which is computationally expensive. In this paper, we implement an image-plane lens modeling methodology for interferometric dirty images by accounting for noise correlations. We show that the image-plane likelihood function produces accurate model values when tested on simulated ALMA observations with an ensemble of noise realizations. We also apply our technique to ALMA observations of two sources selected from the South Pole Telescope survey, comparing our results with previous visibility-based models. Our model results are consistent with previous models for both parametric and pixelated source-plane reconstructions. We implement this methodology for interferometric lens modeling in the open-source software package lenstronomy.

Over several dedicated programs that include targets beyond the main cosmological samples, the Dark Energy Spectroscopic Instrument (DESI) collected spectra for 304,970 unique objects in two fields centered on the COSMOS and XMM-LSS fields. In this work, we develop spectroscopic redshift robustness criteria for those spectra, validate these criteria using visual inspection, and provide two custom Value-Added Catalogs with our redshift characterizations. With these criteria, we reliably classify 212,935 galaxies below z < 1.6, 9,713 quasars and 35,222 stars. As a critical element in characterizing the selection function, we provide the description of 70 different algorithms that were used to select these targets from imaging data. To facilitate joint imaging/spectroscopic analyses, we provide row-matched photometry from the Dark Energy Camera, Hyper-Suprime Cam, and public COSMOS2020 photometric catalogs. Finally, we demonstrate example applications of these large catalogs to photometric redshift estimation, cluster finding, and completeness studies.

As ancient star clusters, globular clusters (GCs) are regarded as powerful tracers of galaxy evolution and assembly. Due to their brightness and compact sizes, GCs are employed to probe the kinematics and stellar population properties of galaxies, from the central regions out into the halo where the underlying stellar light becomes too faint for spectroscopic studies. In this work, we present a comprehensive study of the GC system of M 104 (NGC 4594, also known as the Sombrero galaxy) based on literature spectroscopic catalogues and newly collected data from Very Large Telescope (VLT) MUSE integral-field spectroscopy combined with multi-object spectroscopy from VLT FLAMES and OSIRIS at the Gran Telescopio de Canarias (GTC). We present a new catalogue of 499 GCs with radial velocity measurements that span from the inner disc region out to $\sim$ 70 kpc (24$^{\prime}$). In addition to velocities, we measure metallicities from the MUSE, OSIRIS, and FLAMES spectra of 190 GCs. Together with literature values, we collected a sample of 278 metallicities. Comparing GCs observed with multiple instruments, we find a good agreement of velocity and metallicity measurements. Studying GC kinematics with a simple model confirms a decreasing velocity dispersion profile and low rotation velocities. The blue GCs appear to be more dispersion-dominated, while the red GCs follow the kinematics of the stars more closely. We find a large scatter of GC metallicities with distance from the centre and metal-rich GCs are found over all radii. We discuss that the GC metallicity distribution with a broad metal-poor component likely reflects the complex assembly history of M 104.

Some planetary nebulae (PNe) host X-ray-emitting hot bubbles shaped by stellar wind interactions and/or harbor X-ray-emitting central stars due to accretion, shocks within their fast stellar winds, or even chromospheric emission from binary companions. In both cases, the properties of the X-ray emission critically probe late-stages of stellar evolution for such low- and intermediate-mass stars. While extant Chandra and XMM-Newton observations have detected X-ray emission in PNe, the numbers known remain very small ($\sim40$) compared to the overall Galactic PNe population ($\sim4000$). We have initiated a project aimed at increasing the sample of known PNe with X-ray emission using both current and new space-based X-ray telescopes such as the Einstein probe. To further investigate their X-ray properties to elucidate what drives current X-ray PN detections, we have cross-searched the SRG {\it eROSITA-DE} eRASS1 source catalogue and Hong Kong (HASH) PNe Database. Five known X-ray PNe have been detected (Abell\,30, NGC\,2392, NGC\,3242, NGC\,5315, and LoTr\,5), two new X-ray PNe are revealed (IC\,1297 and NGC\,2867), one (K\,1-27) is removed from previous X-ray compilations, and another 11 previously detected X-ray emitting PNe are not recovered. A comparison of the X-ray flux of detected and undetected X-ray PNe reveals that eROSITA eRASS1 is sensitive to PNe with X-ray fluxes larger than $\approx2\times10^{-14}$ erg~cm$^{-2}$~s$^{-1}$. The frequency of occurrence is $\simeq$0.5\% among the 1430 HASH True PNe in the eRASS1 footprint.

Context. The Kepler Legacy sample is, to this day, the sample of solar-like oscillators with the most exquisite asteroseismic data. In this work, we carry out a detailed modelling of a subsample of these stars for which the surface lithium abundance has also been observed by the LAMOST survey and a photometric surface rotation as been measured. Aims. We aim at studying the impact of additional mixing processes on the asteroseismic modelling of Kepler Legacy G and F-type stars. We also investigate whether a single process can be invoked to reproduce the lithium depletion and asteroseismic constraints at the same time. Methods. We use detailed asteroseismic modelling techniques combining global and local minimization techniques. We start by using standard models and then aim at improving this solution by the addition of extra-mixing at the border of convective regions using either convective penetration or turbulence in radiative layers. Results. We find that lower mass models ($\sim 1M_\odot$) have no problem in reproducing the observed lithium depletion using only turbulence in the radiative zone, similarly to solar models. F-type stars, having a shallower convective envelope, are unaffected by additional turbulence at the BCZ, but require significant convective penetration values to actually reproduce the observed lithium depletion. The extent of this penetration is however incompatible with the frequency separation ratios. Conclusions. We conclude that the impact of extra-mixing is moderate for solar-type stars of the Kepler Legacy sample and well within the requirements of the PLATO mission. For more massive stars ($\sim 1.5M_\odot$), we conclude that the behaviour of the frequency separation ratios must be further investigated, as even models with large convective penetration at the base of their convective envelope are unable to reproduce them.

The characteristic rapid rise and decline at optical wavelengths of a kilonova is the product of the low ejecta mass ($\lesssim 0.05 M_\odot$) and high ejecta velocity ($\gtrsim 0.1$c). We show that, even at very early times ($\lesssim 2$ days), regions of ejecta fall below critical density and temperature thresholds at which non-local thermodynamic equilibrium (NLTE) effects become important. Here, we present an approximate method for calculating the ionization state of the ejecta that accounts for the NLTE impact of high-energy electrons produced in the beta decay of freshly synthesized $r$-process elements. We find that incorporating ionization from high-energy electrons produces an ``inverted" and ``blended" ionization structure, where the most highly ionized species are located in the fastest moving homologous ejecta and multiple ionization states coexist. In radiation transport calculations, the higher degree of ionization reduces line blanketing in optical bands, leading to improved agreement with the light curve properties of AT\,2017gfo such as the duration, decay rates, brightness, and colors. Our quasi-NLTE implementation helps to alleviate tensions in kilonova modeling: for high-velocity ($\sim 0.3c$) ejecta components our models require less mass for a given peak brightness in optical bands, by as much as a factor of 3; our models can explain the presence of observed features associated to Sr II, W III, Se III, and Te III under conditions where LTE models would predict only neutral species; and we naturally predict the coexistence of species like Sr II and Ce III without the need for fine-tuning of the ejecta properties.

We report the discovery of an exceptionally short burst recurrence time in the well-known clocked burster GS 1826$-$238, observed with the CubeSat X-ray observatory NinjaSat. The source had remained in a persistent soft spectral state since its hard-to-soft transition in 2015 July until a soft-to-hard transition occurred in 2025 May. On 2025 June 23, NinjaSat began monitoring GS 1826$-$238 in the hard state and continued until the source returned to a steady soft state. During this period, we detected 19 X-ray bursts: 14 during the hard state, 4 in the transitional state, and 1 in the soft state. In the hard state, we identified a new clocked bursting epoch, during which the burst recurrence time remained highly stable and unprecedentedly short among the clocked bursting phases of GS 1826$-$238, with $t_{\rm rec} = 1.603 \pm 0.040$ hr ($1\sigma$ error). Furthermore, this recurrence time deviates significantly by approximately 37% from the previously established empirical relation $t_{\rm rec} \propto F_{\rm bol}^{-1.05}$, where $F_{\rm bol}$ is the bolometric persistent flux. These results suggest that the ignition conditions for X-ray bursts differed from those in earlier epochs. We explore possible explanations for this discrepancy and propose a scenario in which a smaller fraction of the neutron star surface was involved in fuel accumulation and nuclear burning. This situation leads to a high local accretion rate compared to earlier epochs, which naturally results in the exceptionally short recurrence time, as well as the observed reductions of approximately 38% in blackbody normalization (proportional to the emitting area) and 29% in burst fluence.

Primordial Black Holes (PBHs) are compelling candidates for explaining the present-day relic abundance of cold dark matter (CDM), yet their formation typically requires finely tuned early-universe dynamics. In this work, we propose a novel PBH formation mechanism within a well-established magnetogenesis framework. This scenario simultaneously accounts for the large-scale magnetic fields observed today and generates an enhanced curvature power spectrum at intermediate scales, leading to PBH formation with masses that can survive until the present epoch. We identify a narrow reheating temperature range, $10^5\,\mathrm{GeV} \leq T_{re} \leq 3\times 10^5\,\mathrm{GeV}$, within which the resulting PBHs can constitute the entirety of the observed CDM abundance. Furthermore, our model predicts a stochastic gravitational wave (GW) background as a byproduct of the PBH formation process. Remarkably, the predicted GW signal lies within the sensitivity reach of upcoming space-based interferometers, such as the LISA, DECIGO, or SKA mission, offering a direct observational probe of this PBH generation mechanism.

We propose a new procedure for measuring the spin of a black hole with an unprecedented accuracy based on the gravitational lensing of millisecond pulsars. We derive the basic equations for lensing by a rotating black hole. We show that the frame dragging effect increases the deflection angle of a light ray co-rotating with the black hole. For the primary (secondary) images the angular positions are larger (smaller) for a rotating black hole by an amount on the order of tens of microarcseconds. The differential time delay of images for the case in which the lens is a rotating black hole is smaller than that in the case of non-rotating lens of the same mass, and it can be larger than a few milliseconds. We show that this quantity offers the possibility of reducing the error of spin measurement to less than one percent if we could measure the differential time delay with accuracy of microseconds. We also study relativistic images that are produced by light rays that rotate around the black hole before reaching the observer. The angular positions of relativistic images on the same side as the primary (secondary) image are a few microarcseconds larger (smaller) if the black hole is rotating. Furthermore, the differential time delay between relativistic images is about twelve orders of magnitude larger in the case of rotating lens.

A simple and natural extension of the standard Lambda cold dark matter ($\Lambda$CDM) model is to allow relic neutrinos to have finite chemical potentials. We confront this $\Lambda$CDM$\xi$ model, a $\Lambda$CDM with neutrino mass $M_\nu$ and degeneracy $\xi_3$ as additional parameters, with various cosmological data sets. We find that the $H_0$ and $S_8$ tensions become significant only in the presence of the cosmic microwave background (CMB) polarization data. Specifically, the global and local measurements agree to within 0.8$\sigma$ and 1.6$\sigma$ for the $H_0$ and $S_8$ tensions, respectively, when the CMB polarization data are not included. Therefore, the $H_0$ and $S_8$ tensions exist between CMB temperature and polarization data, both being global measurements. Fitting the $\Lambda$CDM$\xi$ model to the CMB temperature data, we find 3$\sigma$ evidence for nonzero neutrino mass ($M_\nu=0.57^{+0.17}_{-0.13}\,\mathrm{eV}$) and degeneracy ($\xi_3=1.13^{+0.41}_{-0.19}$), and the O(1) neutrino degeneracy parameter is compatible with Big Bang nucleosynthesis data. The scalar index $n_s$ exceeds 1 slightly, which is compatible with some hybrid inflation models. Furthermore, the recent DESI baryon acoustic oscillation data prefer the $\Lambda$CDM$\xi$ model to the Planck $\Lambda$CDM model. Similar results are obtained when including additional supernova data, while the inclusion of the Atacama Cosmology Telescope $\tau$ prior shifts the preferred $M_\nu$ and $\xi_3$ values closer to zero and brings $n_s$ back to the values favored when the polarization data are included.

We propose a revised formulation of General Relativity for cosmological settings, in which the Einstein constant varies with the energy density of the Universe. We demonstrate that this modification has only phenomenological impact of providing an effective dark energy density expression. Assuming a state close to vacuum, here defined by the vanishing product of the Einstein coupling constant and the Universe's energy density, we perform a Taylor expansion of the theory and hence extend it to the whole domain. In this framework, the (renormalized) vacuum energy problem is studied, and an additional constant pressure term, which induces a Chaplygin-like contribution to the dark energy sector, arises in the late-time dynamics. The correction to the late-time Hubble parameter is investigated by comparing theoretical predictions with the late Universe observational data. Our findings indicate that the current value of the stated vacuum energy is consistent with zero within 1$\sigma$. Implications of the modified $\Lambda$CDM model with respect to the Hubble tension are also discussed.

The environment surrounding a black hole or black hole binaries is generally expected to play an important role in understanding various astrophysical phenomena around them. In this paper, we study relativistic, low angular momentum, inviscid, and advective hot accretion flow onto a galactic supermassive black hole dressed with a cold dark matter halo. Focusing on different relativistic dark matter distributions with an inner density spike, we analyze the effect of the dark matter halo on the topology and properties of the accretion flow. Our results show enhancement of disk luminosity in the presence of dark matter, which depends on the nature and properties (halo mass and compactness) of the dark matter distribution. Under the assumptions of our accretion model, the dominant contribution to the disk luminosity for compact and massive halos arises from the inner regions of the accretion flow. Consequently, our analysis indicates that luminosity measurements can serve as an effective probe of the underlying dark matter density spike.

Neural simulation-based inference (SBI) is a popular set of methods for Bayesian inference when models are only available in the form of a simulator. These methods are widely used in the sciences and engineering, where writing down a likelihood can be significantly more challenging than constructing a simulator. However, the performance of neural SBI can suffer when simulators are computationally expensive, thereby limiting the number of simulations that can be performed. In this paper, we propose a novel approach to neural SBI which leverages multilevel Monte Carlo techniques for settings where several simulators of varying cost and fidelity are available. We demonstrate through both theoretical analysis and extensive experiments that our method can significantly enhance the accuracy of SBI methods given a fixed computational budget.

The most recent data release from the Atacama Cosmology Telescope (ACT) reveals a larger value of the scalar spectral tilt $n_s$, ruling out a broad class of inflationary attractors. In this paper, we consider inflationary models including the power law potential, the hilltop model, the polynomial $\alpha$-attractor and exponential $\alpha-$attractor, with non-minimally derivative coupling in the high friction limit, and show how the models can fit ACT data. We also derive constraints on the model parameters using the latest ACT data.