Papers for Friday, May 23 2025

The Cosmic Microwave Background (CMB) data analysis and the map-making process rely heavily on the use of spherical harmonics. For suitable pixelizations of the sphere, the (forward and inverse) Fourier transform plays a crucial role in computing all-sky map from spherical harmonic expansion coefficients -- or from angular power spectrum -- and vice versa. While the Fast Fourier Transform (FFT) is traditionally employed in these computations, the Quantum Fourier Transform (QFT) offers a theoretical advantage in terms of computational efficiency and potential speedup. In this work, we study the potential advantage of using the QFT in this context by exploring the substitution of the FFT with the QFT within the \textit{healpy} package. Performance evaluations are conducted using the Aer simulator. Our results indicate that QFT exhibits potential advantages over FFT that are particularly relevant at high-resolution. However, classical-to-quantum data encoding overhead represents a limitation to current efficiency. In this work, we adopted amplitude encoding, due to its efficiency on encoding maximum data to minimum number of qubits. We identify data encoding as a potential significant bottleneck and discuss its impact on quantum speedup. Future improvements in quantum encoding strategies and algorithmic optimizations could further enhance the feasibility of QFT in CMB data analysis.

We present the large-scale spatial Lya profiles of galaxies from the Keck Baryonic Structure Survey (KBSS) at 2<z<3. This work also describes the Lya imaging for the KBSS-Lya survey for the first time. Our sample includes 734 Lya-selected galaxies and 119 continuum-selected galaxies with Lya narrow-band imaging, and we measure the spatial morphology of Lya and continuum emission for stacked subsamples of these two populations. These samples allow us to study the variation in Lya emission profiles over a broad range of UV continuum luminosities and Lya equivalent widths (EW_Lya), including systems with net Lya absorption in slit spectroscopy. We characterize the spatial profiles using two techniques: directly fitting an exponential function to the stacked profile, and a multi-component forward-modeling technique using the empirical large-scale PSF. We find that both methods yield similar results and that the forward-modeling technique self-consistently fits profiles exhibiting central Lya emission or Lya absorption, with the spatial scale of central Lya approximately matching that of the continuum emission. We also find extended Lya emission such that all our subsamples -- including central Lya absorbers -- are net Lya emitters on scales comparable to the circumgalactic medium (R > 50 kpc, theta > 6''). We find that the scale length of the Lya halo is not strongly dependent on the properties of the central galaxy, including its net continuum luminosity or EW_Lya, although we find a possible weak tendency of continuum-faint, high-EW_Lya galaxies to exhibit larger Lya halos in contrast with previous work.

We present a complete forward model of a realistic weak lensing galaxy catalogue based on the 740 Mpc hydrodynamical MillenniumTNG (MTNG) simulation. Starting with a complete particle and cell lightcone covering one octant of the sky with redshift range 0 < $z$ < 1.5, we apply a group and subhalo finder to generate the corresponding galaxy catalogue for a fiducial observer. For all galaxies, we compute both their intrinsic and lensing-induced shear. The intrinsic component is derived from the luminosity-weighted inertia tensor of stellar particles, while the extrinsic (gravitational) shear is obtained through full-sky ray-tracing on the same lightcone. This allows us to directly predict the impact of intrinsic alignment (IA) of galaxies on the shear correlation function and popular convergence statistics in a fully non-linear forward model. We find that IA modifies the convergence power spectrum at all angular scales by up to 20%, it significantly impacts the PDF, altering its tails by 10-20%, and distorts peak and minimum counts up to 30%, depending on redshift and scale. We also evaluate the impact of the IA signal on the shear correlation function finding that, along with a redshift dependence, the signal strongly increases for higher galaxy stellar mass cuts applied to the catalogue. Notably, with the highest stellar mass cut we apply, the intrinsic shear autocorrelation can become comparable to the gravitational shear component on small angular scales. Our results highlight the importance of accurately modeling IA for precision weak lensing cosmology with upcoming Stage IV surveys.

The effectiveness and stability of mass transfer in binaries system are crucial in determining its final product. Rapid binary population synthesis (BPS) codes simplify the complex physics of mass transfer by adopting parameterized prescriptions for the stability of mass transfer, accretion efficiency in stable mass transfer, and the efficiency of common-envelope ejection. We calibrate these uncertain parameters by comparing BPS models with observational data. White dwarf and main sequence binaries are an ideal population to study binary interaction, as they can be formed through stable or unstable mass transfer, or without interaction, which affect the orbital period and masses of the present-day population. The APOGEE-GALEX-Gaia catalog provides a homogeneous sample of over 500 systems with well measured radial velocities that can be used as a comparison baseline for BPS simulations of such binaries. We compare the distribution of observed maximum radial velocity variation ($\Delta RV_{\rm max}$) and estimated masses to BPS models simulated with COSMIC, varying the mass transfer and common-envelope ejection efficiency, and the criteria for mass transfer stability at key evolutionary stages. The $\Delta RV_{\rm max}$ comparison shows clear preference for a higher fraction of stable mass transfer during the first ascent giant branch, and for highly effective envelope ejection. For the systems with WD masses, there is a slight preference for non-conservative mass transfer. In COSMIC and similar codes, the envelope ejection efficiency and the envelope binding energy are degenerate parameters. Our result of high ejection efficiency may indicate that either additional sources of energy are required to eject the envelope, or that its binding energy is lower than traditionally assumed. Future comparisons to BPS simulations can be drawn for other datasets as they become available.

Binary stars containing a white dwarf and a main-sequence star, WDMS binaries, can be used to study a wide range of aspects of stellar astrophysics. We build a magnitude-limited sample of unresolved WDMS binaries from Gaia DR3 to enlarge these studies. We look for WDMS with available spectra whose location in the Gaia colour-magnitude diagram bridges between the evolutionary sequences of single white dwarfs and the main-sequence. To exclude spurious sources we apply quality cuts on the Gaia photometry and astrometry and we fit the SED (spectral energy distribution) of the objects with VOSA (Virtual Observatory SED Analyser) to exclude single sources. We further clean the sample via visual inspection of the Gaia spectra and publicly available images of the objects. We re-fit the SEDs of the finally selected WDMS with VOSA using composite models to measure their stellar parameters and we search for eclipsing systems by inspecting available ZTF and CRTS light curves. The catalogue consists of 1312 WDMS and we manage to derive stellar parameters for 435. This is because most WDMS are dominated by the main-sequence companions, making it hard to derive parameters for the white dwarfs. We also identify 67 eclipsing systems and estimate a lower limit to the completeness of the sample to be ~50% (~5% if we consider that not all WDMS in the studied region have Gaia spectra). Our catalogue increases by one order of magnitude the volume-limited sample we presented in our previous work. Despite the fact that the sample is incomplete and suffers from heavy observational biases, it is well characterised and can therefore be used to further constrain binary evolution by comparing the observed properties to those from synthetic samples obtained modeling the WDMS population in the Galaxy, taking into account all selection effects.

A stellar common envelope occurs in a binary system when the atmosphere of an evolving star expands to encompass an orbiting companion object. Such systems are predicted to evolve rapidly, ejecting the stellar envelope and leaving the companion in a tighter orbit around a stripped star. We used radio timing to identify a pulsar, PSR J1928+1815, with a spin period of 10.55 ms in a compact binary system with an orbital period of 3.60 hours. The companion star has 1.0 to 1.6 solar masses, eclipses the pulsar for about 17% of the orbit, and is undetected at other wavelengths, so it is most likely a stripped helium star. We interpret this system as having recently undergone a common envelope phase, producing a compact binary.

Recently, it was pointed out that invoking a large value of the CMB optical depth, $\tau_{\rm CMB} = 0.09$, could help resolve tensions between DESI DR2 BAO data and the CMB. This is larger than the value of $\tau_{\rm CMB} = 0.058$ measured from the Planck low-$\ell$ polarization data. Traditionally, $\tau_{\rm CMB}$ is thought of as a constraint on reionization's midpoint. However, recent observations and modeling of the Ly$\alpha$ forest of high-$z$ quasars at $5 < z < 6$ have tightly constrained the timing of the last $10-20\%$ of reionization, adding nuance to this interpretation. Here, we point out that fixing reionization's endpoint, in accordance with the latest Ly$\alpha$ forest constraints, renders $\tau_{\rm CMB}$ a sensitive probe of the duration of reionization, as well as its midpoint. We compare low and high values of $\tau_{\rm CMB}$ to upper limits on the patchy kinematic Sunyaev-Zeldovich (pkSZ) effect, another CMB observable that constrains reionization's duration, and find that a value of $\tau_{\rm CMB} = 0.09$ is in $\approx 2\sigma$ tension with existing limits on the pkSZ from the South Pole Telescope.

The high-luminosity blazar OQ 334 is a leading exponent of the intriguing rare phenomenon of alternating between a flat-spectrum radio quasar (FSRQ) and a BL Lac (BLL) states. Its two optical continuum outbursts observed during the $\sim$ 1.5-year long time span, starting Jan 2018, had been shown to coincide with transition from the FSRQ to BLL state, manifested by a sharp drop in the equivalent width of MgII broad emission-line. Recently, a continuous monitoring of its blazar state, over a much longer duration (past $\sim$ 5 years) has become possible by deploying the observed $\gamma-ray$ spectral slope ($\Gamma_{\gamma}$) as a diagnostic. This opens prospects of making a much less biased and statistically more robust check on the association of optical flaring with FSRQ $\rightarrow$ BLL transition. We find that all 4 such transitions ($\Gamma_{\gamma}$ becoming < 2.0), observed during the past $\sim$ 5 years, were accompanied by an optical flare. While this appears consistent with the transition to BL Lac state happening purely due to an enhanced optical continuum (flaring) swamping out the broad emission-lines, this simple scheme may need additional ingredients, considering the hint found for a day-like offset between the flaring and the state transition.

We show that the gravothermal collapse of self-interacting dark matter (SIDM) halos can deviate from local thermodynamic equilibrium. As a consequence, the self-similar evolution predicted by the commonly adopted conducting fluid model can be altered or broken. Our results are obtained using a novel, efficient kinetic solver called KiSS-SIDM for tracing the gravothermal evolution based on the Direct Simulation Monte Carlo (DSMC) framework. In the long mean free path stage, the code is a viable alternative to the fluid model, yet requires no calibration parameters. Further, this method enables a fully kinetic treatment well into the late, short mean free path, stage of the collapse. We apply the method to a canonical case with isotropic, velocity independent scattering. We find that although a fluid treatment is appropriate deep in the short mean free path core, departures from local thermodynamic equilibrium develop in the intermediate mean free path region bounding the core, which modify the late-time evolution.

Context. Ultra-diffuse galaxies (UDGs) are an intriguing population of galaxies. Despite their dwarf-like stellar masses and low surface brightness, they have large half-light radii and exhibit a diverse range of globular cluster (GC) populations. Some UDGs host many GCs while others have none, raising questions about the conditions under which star clusters (SC) form in dwarfs. GAMA526784, an isolated UDG with both an old stellar body and an extended star-forming (SF) front, including many young SCs, provides an exceptional case to explore the link between UDG evolution and star cluster formation. Aims. This study investigates the stellar populations, SCs, ionised gas, and kinematics of GAMA526784, focusing on its potential to form massive GCs and its connection to broader UDG formation scenarios. Methods. Imaging from HST and Subaru/HSC, alongside MUSE spectroscopy, were used to analyse the galaxy's morphology, chemical composition, and kinematics. A combination of SED fitting and full spectral fitting was applied. Results. GAMA526784's central stellar body exhibits a low-metallicity ([M/H] ~ -1.0 dex) and an old age (~9.9 Gyr). The outskirts are much younger (~0.9 Gyr), but slightly more metal-poor ([M/H] ~ -1.2 dex). The stellar kinematics show low velocity dispersions (~10 km/s) and a coherent rotational field, while the ionised gas exhibits higher dispersions (~50 km/s), a misaligned rotation axis (~20 deg) and localised SF, what could be suggestive of a recent interaction. The young SCs span ages of 8-11 Myr and masses of log(M*/Mo)~5.0, while the old GCs have ~9 Gyr and stellar masses of log(M*/Mo)~5.5. Conclusions. GAMA526784's properties point to interactions that triggered localised SF, leading to the formation of young SCs. Future observations of its molecular and neutral gas will help assess its environment, and the trigger of this SF episode.

Luminous quasars at $z>4$ provide key insights into the early Universe. Their rarity necessitates wide-field multi-band surveys to efficiently separate them from the main astrophysical contaminants (i.e., ultracool dwarfs). To expand the sample of high-$z$ quasars, we conducted targeted selections using optical, infrared, and radio surveys, complemented by literature-based quasar candidate catalogs. In this paper, we report the discovery of \nqsos\ new quasars at $4.6<z<6.9$ (six at $z\geq6.5$), with $M_{1450}$ between $-$25.4 and $-$27.0. We also present new spectra of six $z>6.5$ quasars we selected, but whose independent discovery has already been published in the literature. Three of the newly discovered quasars are strong radio emitters (L$_{1.4~\rm GHz}$$=0.09-1.0\times$10$^{34}$erg s$^{-1}$ Hz$^{-1}$). Among them, one source at $z=4.71$ exhibits typical blazar-like properties, including a flat radio spectrum, radio-loudness $\sim$1000, and multi-frequency variability. It is also detected by SRG/eROSITA X-ray telescope (f$_{\rm 0.2-2.3keV} \sim 1.3\times10^{-13}$erg s$^{-1}$ cm$^{-2}$). In addition, for seven $6.3<z<6.9$ quasars we present near-infrared spectroscopy and estimate the central black hole mass from their C$\rm IV$ and Mg$\rm II$ broad emission this http URL masses (log[M$_{\rm BH,MgII}$]$=8.58-9.14~\rm M_{\odot}$) and Eddington ratios ($\lambda_{\rm Edd,MgII}=0.74-2.2$) are consistent with other $z>6$ quasars reported in the literature. A $z = 6.3$ quasar exhibits a velocity difference of approximately $9000$ km s$^{-1}$ between the C$\rm IV$ and Mg$\rm II$ emission lines, making it one of the most extreme C$\rm IV$ outflows currently known. Additionally, the sample includes three high-ionization broad absorption line quasars. One of these quasars shows potential evidence of an extremely fast outflow feature, reaching $48000$ km s$^{-1}$.

The Nuclear Stellar Disc (NSD) is a rotating, disc-like structure in the Galactic Center, believed to have a distinct star formation. However, its formation history and evolutionary links to other structures in the Galactic Center remain uncertain. This study aims to present the first comprehensive chemical census of the NSD by deriving abundance trends for 18 elements in 9 M giants in the metallicity range of -1.0 <[Fe/H]< 0.5. By comparing these trends with those of other Galactic populations we seek to understand the chemical relationships between these structures. We obtained high-resolution H- and K-band spectra of NSD stars using the IGRINS spectrometer mounted on the Gemini South telescope. The giants were analyzed consistently with stars from a comparison populations to minimize systematic uncertainties. The abundance trends of NSD stars exhibit similarities with those of the inner-bulge and Nuclear Star Cluster (NSC) populations across a broad range of elements with different chemical evolution histories. The trends for 17 elements align closely with the local thick-disc behaviour at subsolar metallicities. At super-solar metallicities, most elements follow the NSC and inner-bulge trends. Sodium is the only element exhibiting a distinct trend, with enhanced abundances in the NSD and NSC compared to both thin-disc and inner-bulge stars. The chemical similarity suggests that the NSD likely shares an evolutionary history with the NSC and possibly the inner-disc sequence. Further studies are required to determine potential evolutionary links to Liller 1 and metal-rich globular clusters. We find no evidence of typical globular cluster abundance signatures in our NSD stars with subsolar metallicities. Our study demonstrates the feasibility of obtaining high-quality abundance data even in highly dust-obscured regions of the Milky Way, paving the way for future surveys.

Currently, 20 co-orbital asteroids of Venus are known, with only one with an eccentricity below 0.38. This is most likely caused by observational biases since asteroids with larger eccentricities may approach the Earth and are easier to detect. We aim to assess the possible threat that the yet undetected population of Venus co-orbitals may pose to Earth, and investigate their detectability from Earth and space observatories. We used semi-analytical models of the 1:1 mean-motion resonance with Venus and numerical simulations to monitor close encounters with Earth on several co-orbital cycles. We analyzed observability windows and brightness variations for potential Venus co-orbitals as viewed from ground-based telescopes to assess their future detection feasibility with next-generation survey capabilities. There is a range of orbits with e < 0.38, larger at lower inclinations, for which Venus' co-orbitals can pose a collisional hazard to Earth. Current ground-based observations are constrained by periodic observing windows and solar elongation limitations, though the Rubin Observatory may detect some of these objects during favorable configurations. Space missions based on Venus' orbits may be instrumental in detecting Venus' co-orbitals at low eccentricities.

The detailed study of gas flows in local Active Galactic Nuclei (AGN) is essential for understanding the regulation of star formation and black hole growth, which are fundamental to galaxy evolution. One such AGN case study is NGC 5643, a nearby ($D_{L}\sim17.3$ Mpc) star-forming, late-type, Seyfert galaxy, where inflows and outflows have been observed in detail. NGC 5643 has been studied at multiple wavelengths, however, a key missing component is sensitive, high-resolution neutral hydrogen ($\mathrm{H\,I}$) observations. We present 21-cm observations of NGC 5643 with MeerKAT, revealing six low-$\mathrm{H\,I}$ mass ($M_{\text{$\mathrm{H\,I}$}}\sim10^{7} M_\odot$) sources surrounding NGC 5643 and $\mathrm{H\,I}$ in IC 4444, $\sim230$ kpc north of NGC 5643. In NGC 5643, $\mathrm{H\,I}$ extends beyond the stellar disk with several morphological and kinematical asymmetries. North of the disk is an extended 30 kpc tail with counter-rotating velocities. This is $\mathrm{H\,I}$ gas accreting onto the regularly rotating disk of NGC 5643 from the environment. Within the spiral arms of the disk, we identify extraplanar gas components, tracing galactic fountains driven by star formation regions. These fountains have a molecular gas component and show an increased $\mathrm{H}_2$/$\mathrm{H\,I}$ ratio. In the circum-nuclear region, we observe spatially unresolved $\mathrm{H\,I}$ absorption that is slightly blue-shifted ($\sim72$ \kms) with an $\mathrm{H\,I}$ emission counterpart at redshifted velocities. These MeerKAT observations provide a complete census of the $\mathrm{H\,I}$ in and around this nearby Seyfert galaxy, providing missing information on the cold gas flows fuelling the star formation and nuclear activity.

NIRCam and MIRI coronagraphy have successfully demonstrated the ability to directly image young sub-Jupiter mass and mature gas-giant exoplanets. However, these modes struggle to reach the sensitivities needed to find the population of cold giant planets that are similar to our own Solar System's giant planets ($T_{\rm eff} = 60 - 125$ K; $a=5 - 30$ AU). For the first time, we explore the high-contrast imaging capabilities of MIRI imaging rather than coronagraphy. Using data from the JWST GO 6122: Cool Kids on the Block program which targets nearby ($<6$ pc) M-dwarfs with NIRCam coronagraphy and MIRI imaging, we demonstrate that 21$\mu$m MIRI imaging can detect planets with the same temperature, mass, age, and orbital separations as Saturn and Jupiter. For systems within 3pc, 21$\mu$m MIRI imaging reaches the sensitivity needed to detect planets colder than Saturn ($<95$ K). NIRCam coronagraphy can achieve similar results only in the unlikely case that a cold giant planet is cloud-free. Motivated by these compelling findings, we extend our analysis to evaluate the measured performance of MIRI F2100W imaging versus NIRCam F444W coronagraphy to 70 pc and conclude that MIRI imaging offers the advantage for systems within 20pc. Microlensing surveys predict an occurrence rate as high as 1 - 2 low-mass giant exoplanets per star, suggesting that JWST MIRI imaging surveys of nearby systems may be poised to uncover a substantial population. This breakthrough enables a path towards the first direct characterization of cold giant exoplanets that are analogous to the solar system giant planets.

High-velocity clouds (HVCs) are interstellar gas clouds whose velocities are incompatible with Galactic rotation. Since the first discovery of HVCs in 1963, their origins have been debated for decades but are still not settled down, because of the lack of vital parameters of HVCs, e.g., the distance. In this work, we determined the distance to the high-velocity cloud, namely the Anti-Center Shell (ACS). We trace the ACS with extinction derived from K-giant stars with known distances and with the diffuse interstellar band (DIB) feature at 5780 $\text{Å}$ fitted on spectra of O- and B-type stars with distance. As a result, we provide a lower limit distance of ACS as $\sim$8 kpc, which extends the lower limit outward by approximately 4 kpc compared to previous work. A byproduct of the DIB method is that we detected a bar-shaped structure with a unusually high positive line-of-sight velocity. Its shape extends along the (l,b)=(155,-5)$^{\circ}$ sight-line and shows a slightly increasing trend in equivalent width and velocity as the distance increases.

The dwarf galaxy Peekaboo was recently identified as a Local Volume (LV) gas-rich and 'eXtremely Metal-Poor' (XMP) dIrr. Its gas metallicity is Z~Zsun/50, with +/-1$\sigma$ uncertainty range of [Zsun/72-Zsun/35]). Its the "Tip of Red Giant Stars" (TRGB)-based distance is of 6.8$\pm$0.7 Mpc. HST data for its individual stars reveal that its older RGB stars comprise a smaller part, while the majority of visible stars have ages of less than one to a few Gyr. Thus, Peekaboo dwarf can be considered as the nearest record-low-Z dwarf. As such, the galaxy deserves a deeper multi-method study, including properties of both, young massive stars and the fainter older population, and its ionised gas and the dominant baryonic component of HI gas. We use the direct (Te) method for the east HII region, in which the [OIII]4363A line is well detected, to estimate its parameter 12+log(O/H). Since in the west HII region the line [OIII]4363A is not detected, its O/H is estimated via the empirical "strong-line" method. The resulting value of O/H is very close to that in the east HII region. The new spectroscopy of Peekaboo dwarf allows us to improve substantially the accuracy of its direct O/H estimate, which appears of 12+log(O/H) = 6.99$\pm$0.06. The new data reveal that emission lines in the E region consist of two components with the velocity difference of ~65 km/s. The fainter, approaching, component can be related to a fast-moving WR star thrown from a cluster or a binary system. Using the HST $V$ magnitudes and colour $V-I$, we identify tentative O-type and very hot candidate WO stars, which are likely the ionising stars of the studied HII regions. With the new optical spectra, the Peekaboo galaxy is confirmed as the lowest-metallicity dwarf in the Local Volume and the valuable object for indeep multi-method studies.

Mantle convection simulations are an essential tool for understanding how rocky planets evolve. However, the poorly known input parameters to these simulations, the non-linear dependence of transport properties on pressure and temperature, and the long integration times in excess of several billion years all pose a computational challenge for numerical solvers. We propose a physics-based machine learning approach that predicts creeping flow velocities as a function of temperature while conserving mass, thereby bypassing the numerical solution of the Stokes problem. A finite-volume solver then uses the predicted velocities to advect and diffuse the temperature field to the next time-step, enabling autoregressive rollout at inference. For training, our model requires temperature-velocity snapshots from a handful of simulations (94). We consider mantle convection in a two-dimensional rectangular box with basal and internal heating, pressure- and temperature-dependent viscosity. Overall, our model is up to 89 times faster than the numerical solver. We also show the importance of different components in our convolutional neural network architecture such as mass conservation, learned paddings on the boundaries, and loss scaling for the overall rollout performance. Finally, we test our approach on unseen scenarios to demonstrate some of its strengths and weaknesses.

This exploratory study investigates whether seismic signals can be used to infer fragmentation during a fireball event. Re-entry objects, particularly sample return capsules (SRCs) such as the one from the Hayabusa2 mission, behave similarly to slow meteors during atmospheric entry and provide valuable insights into natural fireball events. In this study, we initially analyse seismic signals from the Hayabusa2 SRC re-entry, which took place on December 5, 2020, over South Australia. The SRC's signature was captured by a dense network of seismic stations (Eakin, 2018; O'Donnell et al., 2020), offering a unique opportunity to investigate the signals' characteristics and verify their connection to the re-entry event. The ballistic trajectory was confirmed as the source shock mechanism for this event. We isolate this signal and use it as a reference for a ballistic shock signature and compare it to three other fireball case studies, including a suborbital re-entry and two natural meteoroids. Although factors such as local geology and atmospheric conditions were not considered in this preliminary study, our results show promise, with high correlations for events with purely ballistic trajectories and lower correlations for those involving fragmentation or airbursts. This implies that seismic data may be able to disambiguate whether any particular fireball event underwent significant fragmentation or airburst, key phenomena for assessing body strengths.

We here report the results from spectroscopic observations of a sample of 26 globular cluster (GC) and 21 faint fuzzy (FF) candidates in the lenticular galaxy NGC 1023 using the 10.4-m Gran Telescopio Canarias. Using the recessional velocities and stellar absorption features, we determine that 18 and 9 of the observed candidates are bona fide GCs and FFs, respectively. The majority of the rejected FF candidates are background emission line galaxies for which we determine their redshifts. We used the spectroscopic data to determine velocity, age, metallicity and extinction of all bona fide clusters. We find that FFs are clearly younger (age = 7-9 Gyr) than GCs (age > 10 Gyr). Both kind of clusters in this galaxy are metal-rich ([Fe/H] = -0.58 $\pm$ 0.33). The ages and metallicities of individual FFs reported here are the first such measurements in any galaxy and agree with the previously-reported measurement on stacked spectrum. The kinematical analysis reaffirms that the FFs belong to the disk of the galaxy, suggesting that their progenitors are most likely massive, compact disk clusters that have been able to survive for long timescales. We propose that the fuzzy appearance of FFs as compared to the GCs is a consequence of the dynamical evolution of their progenitor super star clusters in the disks of low-mass galaxies.

The peak of the matter power spectrum, known as the turnover (TO) scale, is determined by the horizon size at the time of matter-radiation equality. This scale can serve as a standard ruler, independent of other features in the matter power spectrum, such as baryon acoustic oscillations (BAO). Here, we present the first detection of the turnover in the galaxy auto-power spectrum, utilising the distribution of quasars (QSO) and luminous red galaxies (LRG) measured by the Dark Energy Spectroscopic Instrument (DESI) during its first year of survey operations in a model-independent manner. To avoid confirmation bias, we first analyse the data using data blinding methods designed for the DESI baryon acoustic oscillation, redshift space distortion and scale-dependent bias signals. We measure the angle-averaged dilation distance $D_V(z = 1.651) = (38.1\pm 2.5)r_H$ from the quasars and $D_{V}(z = 0.733) = (21.8\pm 1.0)r_H$ from the LRG sample in units of the horizon $r_H$ at the matter-radiation-equality epoch. Combining these two constraints and assuming a flat $\Lambda$CDM model with three standard neutrino species, we can translate this into a constraint of $\Omega_{m}h^2 = 0.139^{+0.036}_{-0.046}$. We can break the $\Omega_m$-$H_0$ degeneracy with low-redshift distance measurements from type-Ia supernova (SN) data from Pantheon+, we obtain a sound-horizon free estimate of the Hubble-Lemaître parameter of $H_0=65.2^{+4.9}_{-6.2}$ km/s/Mpc, consistent with sound-horizon dependent DESI measurements. On the other hand, combining the DESI BAO and TO, we find a truly DESI-only measurement of $H_0=74.0^{+7.2}_{-3.5}$ km/s/Mpc, in line with DESI-only full-shape results where the sound-horizon scale is marginalised out. This discrepancy in $H_0$ can be reconciled in a $w_0w_a$CDM cosmology, where the combination of DESI BAO and TO data yields $H_0 = 66.5\pm 7.2\;\mathrm{km/s/Mpc}$.

Accurate redshift determinations of both lenses and sources are critical for confirming strong-lens systems and fully realizing their scientific value. However, the thousands of strong-lens candidates now routinely discovered in wide-field imaging surveys make one-by-one follow-up observations impractical. In this work, we investigate the capability and efficiency of large-scale spectroscopic surveys in confirming strong-lens systems. As a case study, we cross-match strong lens candidates identified from the Hyper Suprime-Cam Subaru Strategic Program with Data Release 1 of the Dark Energy Spectroscopic Instrument (DESI). We find that DESI DR1 serendipitously observed putative lenses and/or lensed images in approximately 50\% of these candidates. Analyzing the DESI spectra for $\approx 500$ matched candidates that meet our specific interests, we determine both lens and source redshifts for 27 systems. Additionally, 76 candidate systems feature lensing galaxies at $z > 0.8$, and one candidate system contains a quasar within its lensing galaxy. Applying this approach to other strong-lens candidates will yield numerous additional confirmations, with a further several-fold increase anticipated from the final DESI data release. Our results highlight the growing importance of large-scale spectroscopic surveys in advancing strong lensing discoveries and investigations.

The possibility of the emergence of a stratified region in the uppermost part of the Earth's outer core with long-term magnetic field generation is assessed, taking into account uncertainties in the thermal conductivity of the Earth's core and the present-day heat flow across the core-mantle boundary (CMB). The radial structures of the Earth's outer core are calculated for various values of thermal conductivity and CMB heat flow using a one-dimensional thermo-chemical model. The results show that there exist solutions that allow both emergence of stable stratification and long-term magnetic field generation although their thickness of stratified region is thinner than 100 km. In order to satisfy both emergence of stratified region and long-term magnetic field generation, possible value of the present-day CMB heat flow (13-15 TW) suggests a thermal conductivity of 77-121 W/m/K at CMB, which is in good agreement with the values estimated from the electrical conductivity measurements under the Earth' s core condition. The thickness of the stratified region in this case is about 50 km, which is also consistent with the thickness of the stratified region estimated from the geomagnetic secular variation. However, the proposed values of thermal conductivity obtained by this analysis could be smaller when the present-day CMB heat flow becomes smaller than the constraint used in this study.

The H I Lyman-alpha (Ly$\alpha$) emission, with a wavelength of 1216 Å, is the brightest solar ultraviolet (UV) line. However, comprehensive observations of the Ly$\alpha$ emission line across the full solar disk remain limited. As part of the ASO-S mission, the Solar Disk Imager (SDI) has successfully captured full-disk images in the Ly$\alpha$ band. Gaussian fitting of SDI's spectral response function (SRF) yields a full width at half maximum (FWHM) of approximately 85 Å, which is significantly broader than the distance of Si III line at 1206 Å and the Ly$\alpha$ line. Thus, the emission contribution of Si III to the SDI Ly$\alpha$ passband needs to be considered. For flares, in practice, we calculated the integrated intensity ratio $I$(Si III)/$I$(Ly$\alpha$) by analyzing spectral observations from the SOLSTICE instrument. It yields values between 1.7% and 14.6%. Empirically, the ratio is proportional to the SXR flux. Further analysis of spectral data from the SUMER instrument reveals that the ratio $I$(Si III)/$I$(Ly$\alpha$) is approximately 0.5% for prominences, 0.7%--0.9% for the inner disk, and 1.4%--1.9% close to the limb. These findings suggest that $I$(Si III)/$I$(Ly$\alpha$) is minimal for prominences and the inner disk, and the varying ratios across regions align with the center-to-limb variation of the Si III and Ly$\alpha$ lines. Additionally, we compared Ly$\alpha$ image intensity with 304 Å, 1600 Å, and 1700 Å observations from AIA, as well as H$\alpha$ from CHASE, in multiple regions (a prominence region, two active regions, and a quiet region). A relatively higher correlation of about 85% is found between Ly$\alpha$ and 304 Å in active regions, whereas in the quiet region and prominence, their correlation coefficients are about 55%.

Robust exoplanet characterization studies are underway, and the community is looking ahead toward developing observational strategies to search for life beyond our solar system. With the development of life detection approaches like searching for atmospheric chemical species indicative of life, chemical disequilibrium has also been proposed as a potentially key signature for life. Chemical disequilibrium can arise from the production of waste gases due to biological processes and can be quantified using a metric known as the available Gibbs free energy. The main goal of this study was to explore the detectability of chemical disequilibrium for a modern Earth-like analog. Atmospheric retrievals coupled to a thermodynamics model were used to determine posterior distributions for the available Gibbs free energy given simulated observations at various noise levels. In reflected light, chemical disequilibrium signals were difficult to detect and limited by the constraints on the CH4 abundance, which was challenging to constrain for a modern Earth case with simulated observations spanning ultraviolet through near-infrared wavelengths with V-band SNRs of 10, 20, and 40. For a modern Earth analog orbiting a late-type M dwarf, we simulated transit observations with the James Webb Space Telescope Mid-Infrared Instrument (MIRI) and found that tight constraints on the available Gibbs free energy can be achieved, but only at extremely low noise on the order of several ppm. This study serves as further proof of concept for remotely inferring chemical disequilibrium biosignatures and should be included in continuing to build life detection strategies for future exoplanet characterization missions.

Formaldehyde (H$_2$CO) and methanol (CH$_3$OH) have served as traditional tracers of the star formation process for decades. Studies of the environments which produce these species, though, have pointed to significant differences in the physical environments within which each molecule resides. In this paper we investigate the physical and chemical conditions which give rise to formaldehyde and methanol emission in the nearby starburst galaxy NGC 253. We employ high spatial (1.$''$6 or $\sim28$ pc) and spectral ($\sim10$ km/s) imaging of the NGC 253 central molecular zone (CMZ) from the ALCHEMI Large Program to constrain radiative transfer models of the dense gas volume density, and temperature, molecular species column density, and source filling factor within eight giant molecular clouds (GMCs). We also measure the relative abundances of the two nuclear spin isomers of CH$_3$OH to investigate its formation history. The physical and chemical conditions derived clearly indicate that H$_2$CO and CH$_3$OH originate from distinct physical environments. H$_2$CO traces low volume density and high kinetic temperatures, while CH$_3$OH traces high volume density and low kinetic temperatures. The H$_2$CO abundances are constant, though poorly constrained, within the eight NGC 253 GMCs analyzed, while the CH$_3$OH abundance shows a radial gradient from low to high values within the NGC 253 CMZ. Our findings highlight the complex chemical and physical differentiation of CH$_3$OH and H$_2$CO in the starburst environment of NGC 253. Methanol formation appears to be influenced by warm, dynamic processes rather than cold cloud chemistry, while formaldehyde primarily forms via gas-phase reactions. These results challenge the assumption of a direct chemical link between CH$_3$OH and H$_2$CO and underscores the impact of starburst-driven shocks, turbulence, and cosmic rays on molecular gas chemistry.

Eclipsing millisecond pulsars (MSPs) are a type of pulsar binaries with close orbits ($\lesssim1.0\,$d). They are important objects for studying the accretion history of neutron stars (NSs), pulsar winds, and the origin of isolated MSPs, etc. Recently, a new eclipsing MSP, PSR J$1928+1815$, was discovered by the Five-hundred-meter Aperture Spherical radio Telescope. It is the first known pulsar with a He star companion, as suggested in Yang et al. The system features a short orbital period of $\sim0.15\,$d and a relatively massive companion $\gtrsim1.0\,M_\odot$. However, the origin of PSR J$1928+1815$ remains highly uncertain. In this paper, we investigated the formation of the new subclass of eclipsing MSPs containing (evolved) He star companions through NS + He star channel. We found that if a NS binary undergoes subsequent mass-transfer phases following Case BA or Case BB, it may appear as an eclipsing MSP during the detached phase. Additionally, we obtained the initial parameter space for producing eclipsing MSPs with He star companions. Using binary population synthesis approach, we estimated their birth rate to be $\sim2.1-4.7\times10^{-4}\rm\,yr^{-1}$, corresponding to a total number of $\sim55-150$ systems in the Galaxy. Moreover, we concluded that PSR J$1928+1815$ may originate from the evolution of an NS+He star system with an initial orbital period of $\sim0.1\,$d, which can undergo the Case BB mass transfer.

Sedna-like objects (a.k.a. sednoids) are transneptunian objects (TNOs) characterized by large semimajor axes and exceptionally high perihelia. Their high-$q$ orbits are detached from the influence of the four giant planets and need extra perturbation to form. One hypothesis posits that close stellar flybys could have perturbed objects from the primordial scattering disk, generating the sednoid population. In this study, we run N-body simulations with different stellar encounter configurations to explore whether such a close stellar flyby can satisfy new constraints identified from sednoid (and detached extreme TNO) observation, including the low-inclination ($i<30^\circ$) profile and primordial orbital alignment. Our results suggest that flybys with field stars are unable to generate a sufficient population, whereas flybys within the birth cluster fail to produce the primordial orbital alignment. To meet the inclination constraint of detached extreme TNOs, flybys have to be either coplanar ($i_\star \sim 0^\circ$) or symmetric about the ecliptic plane ($\omega_\star \sim 0^\circ, i_\star \sim 90^\circ$). After taking into account their occurrence rate at the early stage of the Solar System, we conclude that close-in stellar flybys ($q_\star \le 1000$~au) that satisfy all constraints are unlikely to happen ($\lesssim$5\%). Future discoveries of additional sednoids with precise orbital determinations are crucial to confirm the existence of the low-inclination tendency and the primordial alignment, and to further constrain the early dynamical evolution of the Solar System.

Gaia Data Release 3 (DR3) published for the first time epoch photometry, BP/RP (XP) low-resolution mean spectra, and supervised classification results for millions of variable sources. This extensive dataset offers a unique opportunity to study their variability by combining multiple Gaia data products. In preparation for DR4, we propose and evaluate a machine learning methodology capable of ingesting multiple Gaia data products to achieve an unsupervised classification of stellar and quasar variability. A dataset of 4 million Gaia DR3 sources is used to train three variational autoencoders (VAE), which are artificial neural networks (ANNs) designed for data compression and generation. One VAE is trained on Gaia XP low-resolution spectra, another on a novel approach based on the distribution of magnitude differences in the Gaia G band, and the third on folded Gaia G band light curves. Each Gaia source is compressed into 15 numbers, representing the coordinates in a 15-dimensional latent space generated by combining the outputs of these three models. The learned latent representation produced by the ANN effectively distinguishes between the main variability classes present in Gaia DR3, as demonstrated through both supervised and unsupervised classification analysis of the latent space. The results highlight a strong synergy between light curves and low-resolution spectral data, emphasising the benefits of combining the different Gaia data products. A two-dimensional projection of the latent variables reveals numerous overdensities, most of which strongly correlate with astrophysical properties, showing the potential of this latent space for astrophysical discovery. We show that the properties of our novel latent representation make it highly valuable for variability analysis tasks, including classification, clustering and outlier detection.

We investigate the impact of dark matter on the structure and deformation of magnetars. We assume a perturbative approach for the magnetic field deformation and that the dark matter only interacts gravitationally with hardonic matter. Assuming that the DM is significantly softer than HM, we find that the magnetic field can affect dark matter through the deformation of space-time. The number of stars having a dark matter halo outside the hadronic matter surface increases with an increase in dark matter fraction and the stiffness of the dark matter equation of state. As the magnetic field deforms the stars from sphericity, we can have a situation where we have a non-symmetric dark matter halo outside the star. This can have interesting observational gravitational signatures unique to magnetars with having DM halo.

Relativistic jets can be produced within the accretion disk of an active galactic nucleus (AGN), leading to distinct thermal emission as they propagate through a dense disk environment. In this paper, we present a comprehensive study of dynamical evolution of jets embedded in an AGN disk and their associated observational properties, focusing on scenarios in which jets either successfully break out of the disk or become choked. By modeling the jet-cocoon system propagation, we calculate the thermal emission contributions from the jet-head shock breakout, disk cocoon, and jet cocoon components. Our results reveal that soft X-ray flares are the most prominent observable signatures, with duration ranging from O(10^2) s to O(10^5) s, occasionally exhibiting double-peaked light curves, whereas UV/optical flares are detectable only for powerful jets, persisting for several days to tens of days. This thermal emission serves as a critical electromagnetic counterpart to jet-producing events and provide insights into jet dynamics and AGN disk properties. Our findings highlight the importance of multi-wavelength follow-up observations to establish a diagnostic paradigm for candidate electromagnetic counterpart identification to AGN-embedded events and to distinguish thermal flares from AGN background variability.

Weak lensing alters galaxy sizes and fluxes, influencing the clustering patterns of galaxies through cosmic magnification. This effect enables the reconstruction of weak lensing convergence $\hat{\kappa}$ maps for DES and DECaLS by linearly combining galaxy overdensities across magnitude bins in the $g$, $r$, and $z$ photometry bands \citep{Qin+,Qin2+}. In this study, we enhance the lensing reconstruction method by addressing biases in the magnification coefficient estimation, which arise from incomplete consideration of selection effects, especially those induced by photometric redshift (photo-$z$) selection. Using a Random Forest-based photo-$z$ estimation for DECaLS and DES galaxies, we quantify the impact of photo-$z$ induced selection on magnification coefficient estimation. Our results show that neglecting photo-$z$ selection introduces significant biases in the magnification coefficient, leading to deviations in the reconstructed convergence map amplitude $A$, with values ranging from 0.4 to 3.5 depending on the survey, redshift, and magnitude cuts. By incorporating an improved magnification coefficient estimation that accounts for photo-$z$ selection, these biases are significantly reduced, with $A$ converging to $\sim 1$ as the magnitude cuts approach optimal values. This improvement is consistently observed across DES and DECaLS datasets and redshift bins, despite differences in survey strategies and depths. Our findings highlight the importance of addressing photo-$z$ induced selection to achieve unbiased weak lensing reconstructions and accurate cosmic magnification measurements.

The GRAPES-3 large area muon telescope with its sixteen independent modules records the high energy (>1 GeV) muons continuously over 2.3 sr of the sky. However, the recorded muon rates are contaminated by instrumental effects and instabilities spanning both short- and long-timescales, such as variations in the efficiency of the detector. We present an automated, algorithmic method, which employs Bayesian blocks to discretize the data stream into periods and exploits the correlations among the sixteen independent modules of the muon telescope to separate the impact of these instrumental problems from those originating in physical effects of interest, allowing the Savitzky-Golay filter to be employed to mitigate the former. Compared to legacy methods, this method is less dependent on subjective input from experimental operators and provides a data stream free of all known instrumental effects over calendar years. The muon rate obtained with the new method shows a fairly better correlation with neutron monitor data, than that obtained with the legacy method.

Classical Wolf-Rayet (cWR) stars are evolved massive stars that have lost most of their H envelope and exhibit dense, extended atmospheres with strong, line-driven winds. Accurately modeling wind launching from optically thick layers remains a challenge. Two main approaches have advanced our understanding: 1D stationary atmosphere models with consistent hydrodynamics and time-dependent, multi-dimensional radiation-hydrodynamic simulations. Due to high computational demands, multi-dimensional models are limited in scope. Therefore, 1D hydrodynamically consistent models remain essential but must incorporate insights from 3D simulations. We compare averaged stratifications from recent multi-dimensional cWR models with 1D models computed using the hydrodynamically consistent PoWR$^{HD}$ code. We focus on winds driven by the hot iron opacity bump and explore how variations in 1D input parameters affect model outcomes. The 1D models reproduce the average 3D density structure well. While mass-loss rates are typically $\lesssim$0.2 dex higher in 1D models, small adjustments accounting for multi-dimensional dispersion reconcile the differences. 1D models tend to be more radially extended, with higher terminal velocities and lower effective temperatures. They reproduce the general velocity trends of 3D models but launch winds slightly further out and reach higher velocities during the hot iron bump. These differences also manifest in synthetic spectra computed from different 1D model approaches. Despite methodological variations, both 1D and averaged 3D models yield consistent stellar parameters when accounting for the variability seen in time-dependent simulations. For stars near the Eddington limit, reducing Doppler velocities in 1D models improves agreement in mass-loss rates, temperatures, and wind velocities. Matching temperature structures in optically thin layers remains an open challenge.

Identifying AGNs in dwarf galaxies is critical for understanding black hole formation but remains challenging due to their low luminosities, low metallicities, and star formation-driven emission that can obscure AGN signatures. Machine learning (ML) techniques, particularly unsupervised methods, offer new ways to address these challenges by uncovering patterns in complex data. In this study, we apply Self-Organizing Maps (SOMs) to explore the SED manifold of dwarf galaxies and evaluate AGN selection biases across diagnostics. We train a 51 by 51 SOM on 30,344 dwarf galaxies (redshift less than 0.055 and stellar mass below 10 to the 9.5 solar masses) from the NSA catalog using nine-band photometry from near-UV to mid-infrared. A set of 438 previously identified dwarf AGNs, selected via various methods, was mapped onto the SOM. AGNs identified by different methods occupy distinct and partially overlapping regions in SED space, reflecting selection biases tied to host properties. BPT selected AGNs cluster in higher-mass regions, while X-ray and variability-selected AGNs show broader distributions. WISE-selected AGNs are concentrated in lower-mass regions and form two clumps: one associated with bluer, starburst-like systems and the other with redder, more AGN-like SEDs. This separation may help distinguish true AGN hosts from starburst contaminants in WISE-selected samples. AGNs selected via traditional emission-line, broad-line, and WISE methods tend to avoid SOM regions linked to strong star formation. In contrast, a subset of AGNs in low-mass galaxies occupy regions indicative of high AGN luminosity relative to stellar content, highlighting luminous AGNs in faint hosts. These results demonstrate the utility of manifold learning for improving AGN selection in the low-mass regime.

We present a chemo-dynamical $N$-body/hydrodynamic simulation of an isolated Milky Way-like galaxy to investigate how bar formation influences star formation rates, stellar migration, and the resulting age and metallicity distributions of disk stars. Focusing on the transient epoch of bar formation, a phase that triggers gas inflows, enhances local star formation, and drives significant orbital migration, we find that the star formation rate in the inner disk exhibits a pronounced peak during this period. This behavior arises from the combined effect of vigorous star formation driven by strong spiral arms prior to bar formation and the subsequent suppression of star formation once the bar is established. In contrast, star formation in the outer disk persists after bar formation at modest levels, and enhanced outward migration of stars originally formed in the inner regions gives rise to a pronounced peak in the outer disk's stellar age distribution corresponding to the bar formation epoch. Moreover, stars formed during this epoch tend to exhibit higher gas-phase metallicities, reflecting their origin in more metal-rich inner regions. Although our model does not capture every detail of the Milky Way's complex evolution, our results highlight the dominant role of bar driven migration in segregating star formation activity and in shaping the long-term chemical and age structure of the Galactic disk. Recent observational studies suggest that the Milky Way's bar is approximately 8 Gyr old; therefore, our findings imply that the age distribution of stars in the solar circle and outer disk should show a corresponding peak around that age.

M dwarfs are the most abundant stars in the Galaxy and serve as key targets for stellar and exoplanetary studies. It is particularly challenging to determine their metallicities because their spectra are complex. For this reason, several authors have focused on photometric estimates of the M-dwarf metallicity. Although artificial neural networks have been used in the framework of modern astrophysics, their application to a photometric metallicity estimate for M dwarfs remains unexplored. We develop an accurate method for estimating the photometric metallicities of M dwarfs using artificial neural networks to address the limitations of traditional empirical approaches. We trained a neural network on a dataset of M dwarfs with spectroscopically derived metallicities. We used eight absolute magnitudes in the visible and infrared from Gaia, 2MASS, and WISE as input features. Batch normalization and dropout regularization stabilized the training and prevented overfitting. We applied the Monte Carlo dropout technique to obtain more robust predictions. The neural network demonstrated a strong performance in estimating photometric metallicities for M dwarfs in the range of -0.45<[Fe/H]<0.45dex and for spectral types as late as M5.0V. On the test sample, the predictions showed uncertainties down to 0.08dex. This surpasses the accuracy of previous methods. We further validated our results using an additional sample of 46 M dwarfs in wide binary systems with FGK-type primary stars with well-defined metallicities and achieved an excellent predictive performance that surpassed the 0.1dex error threshold. This study introduces a ML-based framework for estimating the photometric metallicities of M dwarfs and provides a scalable data-driven solution for analyzing large photometric surveys. The results outline the potential of artificial neural networks to enhance the determination of stellar parameters.

Ozone is a potential biosignature and disambuguator between Earth-like and Venus-like exoplanets due to its association on Earth with photosynthetically produced oxygen (O$_2$). However, the existence of ozone in Venus's observable atmosphere, a planet with no known life, raises the possibility of ozone biosignature false-positives on Venus-like exoplanets. We use a photochemical model of Venus's atmosphere to investigate the origin of its mesospheric ozone layer, and to predict how similar ozone layers would manifest for Venus-like exoplanets. For Venus, our model shows that the previously proposed fluxes of O atoms produced on the dayside and transported to the nightside cannot generate enough ozone to match the observed nightside ozone concentrations without also producing O$_2$ in excess of the observed upper limit. Nor can sufficient ozone be produced by varying the lower-atmosphere chemistry, atmospheric thermal structure, or received stellar flux in our model of Venus's atmosphere. These results imply that a presently unknown chemical pathway is responsible for the ozone production in Venus's nightside mesosphere. Ozone production rates from this pathway of 10$^5$--10$^7$ cm$^{-3}$s$^{-1}$ above the cloud layer on the nightside can re-produce the observed O$_3$ concentrations. Generalised to Venus-like exoplanets, known chemistry similarly fails to produce ozone in the abundance seen in the Venusian mesosphere. However, until the origin of Venus's ozone is understood, we cannot rule out that ozone production at concentrations observable with JWST will be common on abiotic Venus-like worlds, a possibility that limits the usefulness of ozone as a habsignature and as a biosignature.

The Gaia dataset has revealed many intricate Milky Way substructures in exquisite detail, including moving groups and the phase spiral. Precise characterisation of these features and detailed comparisons to theoretical models require engaging with Gaia's heteroscedastic noise model, particularly in more distant parts of the Galactic disc and halo. We propose a general, novel machine-learning approach using normalizing flows for denoising density estimation, with particular focus on density estimation from stellar survey data such as that from Gaia. Normalizing flows transform a simple base distribution into a complex target distribution through bijective transformations resulting in a highly expressive and flexible model. The denoising is performed using importance sampling. We demonstrate that this general procedure works excellently on Gaia data by reconstructing detailed local velocity distributions artificially corrupted with noise. For example, we show the multiple branches of the Hercules stream and the phase-space spiral can both be well captured by our model. We discuss hyperparameter choice to optimally recover substructure and compare our approach to extreme deconvolution. The model therefore promises to be a robust tool for studying the Milky Way's kinematics in Galactic locations where the noise from Gaia is significant.

We assess the accuracy of radial velocities of 671 stars located near Galactic clusters, measured by Robert Trumpler between 1924 and 1947 using the Lick Observatory 36 inch refractor equipped with prism spectrometers. We find that Trumpler's velocities share the same zero-point and scale as modern IAU radial velocity standards. Their accuracy ranges from 2 to 7 km/s. We also provide star identifications, B-band photometry, and notes for many stars. We provide a link to an online PDF that contains all 3,782 of Trumpler's individual stellar radial velocities and associated Julian dates. Trumpler's velocities provide a century-long baseline for detecting long-period stellar companions, gravitational perturbations from passing massive objects, and accelerations due to cluster and Galactic dynamics.

Magnetothermal (MHD) winds and jets originate in a wide range of regions of protoplanetary disks (1-30 au) and are thought to be the primary mechanisms driving accretion onto the central star. One indirect signature of these processes is the free-free emission from ionized gas close to the star. We analyze a sample of 31 Class II disks: 18 full disks (FD) and 13 transition disks (TD). All sources show evidence of excess free-free emission over the contribution of the thermal dust. We investigate the origin of this emission and whether it is associated with other observables. We first analyzed a sample of objects in Taurus, exploring correlations with the properties of the central star, the disk, and other disk-wind tracers. We compared our findings with a sample of TD for which free-free emission was shown to be likely associated with an MHD-wind/jet. We found no correlation between the detected free-free emission and either the X-ray or the [OI]6300A line properties. We found a strong correlation between the ionized mass loss rate, as inferred from the free-free emission, and the accretion rate, suggesting that free-free emission in FD is associated with an MHD-wind/jet. The detected free-free emission in both TD and FD is likely similarly associated with an ionized gas close to the star from an MHD-wind/jet. The free-free emission detected in TD shows hints of shallower correlations with accretion properties than in FD. Whereas the efficiency in transforming accretion into outflow might differ in TD and FD, considering the correlations between free-free emission and accretion properties, this difference could simply result from a bias toward strong accretors in the TD sample. Therefore, observations of a more complete and uniform sample are necessary to determine whether this change in correlations holds only for strong accretors or for TD in general.

The computation of magnitudes and distances from low signal-to-noise observations is known to be problematic, in the sense that the magnitudes and distances tend to assume extreme values, or are even undefined or unphysical in the case of negative observed fluxes or parallaxes. In this work we show that magnitudes can be computed consistently at all signal-to-noise levels, and even for negative fluxes, if the prior information that the true flux or distance is non-negative is properly included. Furthermore, we derive an all-purpose estimator for distances from a prior implementing only the non-negativity of the true parallax. We apply our results to the case of combining magnitudes to colours, and magnitudes and distances to obtain absolute magnitudes. The resulting expressions are easy to compute and we show that the resulting distribution functions for magnitudes, colours, distances, and absolute magnitudes are not only consistent for all signal-to-noise levels and applicable to both, positive and negative observed fluxes and parallaxes, but also show no strong tails. While biases at very low signal-to-noise levels are unavoidable, the estimator for distances derived in this work is less biased than previously used estimators. We find that the magnitude, colour, distance, and absolute magnitude distributions for vanishing signals converge to limiting distributions, whose median values are important for assessing biases when working with data at low signal-to-noise levels.

It is still unknown how much water has escaped from Mars during its history. Hydrogen escape from Mars's atmosphere probably played a major role in drying the planet, but present-day Hloss rates (about 3x10^26 atoms per second on average) cannot explain the geological evidence for the large volumes of liquid water on ancient Mars. Here we used the three-dimensional Mars-Planetary Climate Model to show that H loss rates could have increased by more than one order of magnitude (6x10^27 atoms per second) during higher spin axis obliquity periods, notably in the last few million years when Mars's obliquity was about 35 deg on average. The resulting accumulated H escape over Mars's history translates into an approx. 80 m global equivalent layer, which is close to the lower limit of geological estimates, assessing the major role of atmospheric escape in drying Mars.

Long gamma-ray bursts (GRBs) are characterized by a brief gamma-ray flash followed by a longer-lasting multiwavelength afterglow. The basic mechanism is largely understood, and the early afterglow evolution often shows complex features that provide crucial insights into the transition between prompt and afterglow phases. We present a detailed analysis of GRB 210312B, detected by INTEGRAL, which exhibits both a precursor and a complex optical afterglow evolution. Through careful modeling using Markov chain Monte Carlo methods, we disentangled the contributions of an early optical flare and forward shock emission. Our analysis reveals a gamma-ray precursor 17 s before the main pulse with a significantly softer spectrum (hardness ratio 0.37 +/- 0.12 versus 1.9 +/- 0.4). The optical afterglow shows an early peak at 76.0^{+4.4}{-5.1} s characterized by a steep rise ({\alpha}{flare,1} = -4.1^{+1.6}{-2.3}) and decay ({\alpha}{flare,2} = 4.0^{+2.1}{-1.5}), followed by forward shock emission with a broad hydrodynamic peak at around 150 s. In the subsequent plateau phase, the afterglow initially has a complex structure before settling into a final power law decay consistent with an electron distribution index p = 2.36^{+0.18}{-0.15}. The negligible host extinction (A_{V,host} = -0.073^{+0.100}_{-0.078}) suggests we are observing the intrinsic afterglow spectrum. The host system consists of two luminous (M_B ~ -21.7) components separated by 11.5 kpc at z = 1.069, which are possibly an interacting galaxy pair. GRB 210312B provides a rare opportunity to study the prompt-to-afterglow transition in detail. The consistency of the forward shock component with standard afterglow theory supports our physical interpretation despite the lack of X-ray coverage.

The observed power-law spectra of relativistic jets from active galactic nuclei clearly indicate a synchrotron mechanism of radiation by particles that similarly possess a power-law energy spectrum. However, the issue of their angular anisotropy has not been given sufficient attention until recently, although the example of the solar wind (where a strongly magnetized wind is realized in a similar way) shows the importance of taking this circumstance into account. In this paper, we study the evolution of an initially isotropic power-law spectrum of radiating particles as they propagate along expanding relativistic jets. It is shown that for relativistic flows in which the electric field plays a crucial role, the preservation of the first adiabatic invariant does not lead to a decrease in the pitch angles of radiating particles as they enter the region of weak magnetic fields. This is due to the drift nature of the particle motion.

Herbig Ae/Be stars are young, pre-main-sequence stars that provide critical insights into the processes of stellar formation, early stellar evolution and protoplanetary this http URL of the key features of such stars are their circumstellar dusty disk and bipolar ionized outflows, which are key components for understanding planet formation processes and energy/matter deposition in the interstellar medium, respectively. In this context, imaging polarimetry is probably the sharpest tool to characterize the various structures and dynamics around the central star, due to the sensitivity of polarization to the morphology of the emitting, scattering and absorbing media. We take advantage of never published, near-ultraviolet polarimetric data of LK-H{\alpha}-233 taken by the Faint Object Camera aboard the Hubble Space Telescope in 1991, 1994 and 1995, which remained dormant in the archives despite their quality. Using the most recent and robust reduction pipeline for this instrument, we obtained high spatial resolution (0.0287 x 0.0287 arcsecond2) maps of this object at 4118 Å, together with polarimetric measurements. A dark lane, bisecting the approaching and receding polar outflows, suggests the presence of a circumstellar disk or dust torus, obscuring the pre-main sequence star and collimating the ejecta. Polarization reveals that the outflows have an X-shape structure with a significant centro-symmetric pattern in polarization angle, indicating that the outflows are both hollow and scattering the emission from the buried star. We constrain the half-opening angle of both the outflows and circumstellar disk, determine the inclination of the system and estimate the obscured star's intrinsic flux. This study highlights the importance of high-resolution polarimetric observations in understanding the complex environment around Herbig Ae/Be stars and advocates for future similar instruments.

The Fornax cluster, the second-largest galaxy cluster within 20 Mpc, presents an ideal environment for studying environmental effects on galaxy evolution. Utilizing data from the Southern Photometric Local Universe Survey (S-PLUS), this study explores the H$\alpha$+[NII] emission maps across an area of approximately 208 square degrees around NGC 1399. For such, a dedicated semi-automated pipeline, Pixel-to-Pixel Emission Line Estimate (PELE), was developed to generate emission line maps by processing S-PLUS images using the Three Filter Method. A morphological analysis was conducted using the ASTROMORPHLIB package to determine whether H$\alpha$+[NII] emitters exhibit perturbed features. The study successfully detected 77 H$\alpha$+[NII] emitters with $r<18$ mag, extending to four times the virial radius of the Fornax cluster. PELE demonstrated its ability to recover flux down to 2e-17 erg s$^{-1}$ cm$^{-2}$ when compared to H$\alpha$ maps from MUSE/VLT. Among the emitters, 25% are early-type galaxies (ETG) and 75% late-type galaxies (LTG). Signs of morphological perturbation or merger activity are observed in 44% of the LTG and in three ETG located beyond the cluster's virial radius. A significant fraction (91%) of the emitters are identified as recent infallers, primarily located in the northwestern region of the cluster, while others are associated with the infalling group Fornax A in the southwest. Disturbed, low-mass galaxies at larger cluster-centric distances provide evidence of galaxies begin transforming before entering the main cluster. This study demonstrates S-PLUS's effectiveness in detecting emitters, whose distribution reflects the Fornax cluster's assembly history, with LTG linked to recent infall from the field, possibly along a Fornax-Eridanus filament, and ETG may have evolved prior to entry.

We investigate the nonlinear evolution of cosmological perturbations in theories with scale-dependent perturbation growth, first in general and then focusing on Horndeski gravity. Within the framework of standard perturbation theory, we derive the second- and third-order kernels and show that they are fully determined by two effective functions, $h_1$ and $h_c$, which parametrize deviations from general relativity. Using the Wronskian method, we obtain solutions for the nonlinear growth functions and present explicit expressions for the resulting kernels, including bias and redshift space distortions. We show that the kernels are entirely dependent on the linear growing mode: once this is calculated, the kernels are analytic up to a time integral. Our approach provides a physically motivated framework for evaluating the one-loop galaxy power spectrum in scale-dependent theories, suitable for the forecasts and actual data analysis.

IceCube is a Cherenkov detector instrumenting over a cubic kilometer of glacial ice deep under the surface of the South Pole. The DeepCore sub-detector lowers the detection energy threshold to a few GeV, enabling the precise measurements of neutrino oscillation parameters with atmospheric neutrinos. The reconstruction of neutrino interactions inside the detector is essential in studying neutrino oscillations. It is particularly challenging to reconstruct sub-100 GeV events with the IceCube detectors due to the relatively sparse detection units and detection medium. Convolutional neural networks (CNNs) are broadly used in physics experiments for both classification and regression purposes. This paper discusses the CNNs developed and employed for the latest IceCube-DeepCore oscillation measurements. These CNNs estimate various properties of the detected neutrinos, such as their energy, direction of arrival, interaction vertex position, flavor-related signature, and are also used for background classification.

The scalar-induced gravitational waves (SIGWs) produced nonlinearly by the enhanced cosmological curvature perturbations can serve as a potentially powerful probe of primordial non-Gaussianity (PNG) in the early Universe. In this work, we comprehensively investigate the imprints of local-type PNG on the SIGW background beyond the widely used quadratic and cubic approximations. We develop a diagrammatic approach capable of analyzing SIGWs for PNG up to any order. Utilizing this approach, we derive semi-analytic formulas for the energy-density fraction spectrum of the isotropic background, the angular power spectrum of the energy-density anisotropies, as well as the angular bispectrum and trispectrum of the energy-density non-Gaussianity. Specializing to PNG up to quartic approximation (parameterized by $f_\mathrm{NL}$, $g_\mathrm{NL}$, and $h_\mathrm{NL}$), we numerically compute all contributions to these SIGW spectra. We find that PNG can significantly alter the amplitude of the energy-density spectrum and generate substantial anisotropies through the initial inhomogeneities. Furthermore, we observe that the angular bispectrum and trispectrum always vanish when the primordial curvature perturbations are Gaussian; otherwise, they do not, indicating their potential utility as probes of PNG. Therefore, we anticipate that both the isotropic background and the anisotropies of the SIGW background will provide essential information about the early Universe.

We investigate the effect of planetary corotation on energy dissipation within the magnetosphere-ionosphere system of exoplanets. Using MHD simulations, we find that tidally locked exoplanets have a higher cross-polar cap potential (CPCP) compared to fast-rotating planets with the same magnetic field strength, confirming previous studies. Our simulations show that for a given interplanetary magnetic field, an increase in corotation period leads to a higher CPCP. Notably, this difference in CPCP between tidally locked and rotating planets persists across a range of solar wind conditions, including extreme environments such as those experienced by hot Jupiters. Furthermore, we observe that variations in corotation have little impact on CPCP for Earth-sized planets. These results underscore the significance of both corotation dynamics and planetary size in understanding how exoplanets interact with their stellar environments.

We investigate primordial magnetic fields (PMFs) as a potential source of the parity-violating signatures recently implied by observations of the cosmic microwave background (CMB) anisotropies and the large-scale structure of the Universe. Among various possibilities, helical PMFs are of particular interest, as they inherently violate parity symmetry and can explain the observed magnetic fields, especially in void regions. PMFs, if generated in the early universe, can source curvature perturbations, which evolve into the present density fluctuations observed in galaxy surveys. Motivated by this, we focus on the imprint of helical PMFs on the trispectrum of the sourced primordial curvature perturbations, a leading-order statistics sensitive to parity-violating signals in three-dimensional space. We derive the analytic expressions for the trispectrum of the primordial curvature perturbations sourced by taking into account both the helical and non-helical PMFs and analytically reduce their expressions by using the pole approximation, whose validity is confirmed by comparison with the exact results. We find that, varying the ratio of the amplitude of the helical to non-helical power spectrum, the structure of the trispectrum qualitatively changes, particularly the ratio of the imaginary to real components of the trispectrum. Our findings highlight the primordial trispectrum as a promising probe of cosmological parity violation in the early universe and provide a theoretical basis for future precise observations of higher order statistics in the CMB anisotropies and the galaxy clustering.

In this work we consider a generalization of Starobinsky inflation obtained by combining power law ($R^\beta$), and $\alpha$-Starobinsky inflation ($E$-model). The Einstein frame potential for this model is that of power law Starobinsky inflation modified by a parameter $\alpha$ in the exponential. After computing power spectra for scalar and tensor perturbations numerically, we perform MCMC analysis to put constraints on the potential parameter $\alpha$, $\beta$ and $M$, and the number of e-foldings $N_{pivot}$ during inflation, using Planck-2018, BICEP/Keck (BK18) and other LSS observations. We find $\log_{10}\alpha= 0.37^{+0.82}_{-0.85}$, $\beta = 1.969^{+0.020}_{-0.023}$, $M=\left(3.54^{+2.62}_{-1.73}\right)\times 10^{-5}$ and $N_{pivot} = 47\pm{10}$. We compute the Bayesian evidences for our proposed model, power law Starobinsky inflation, $\alpha$-Starobinsky inflation and Starobinsky inflation. Considering the Starobinsky model as the base model, we calculate the Bayes factor and find that our proposed model is preferred by the CMB and LSS observations.

We introduce a non-parametric approach, the Stellar Population Synthesis with Equivalent Widths (SEW) method, to reconstruct spectral-resolution wavelength-dependent attenuation curves for 169,568 star-forming galaxies from the SDSS DR7. Composite attenuation curves, stacked across stellar mass and inclination bins, reveal systematic trends: higher stellar mass correlates with steeper attenuation slopes (lower $R_V$), while edge-on galaxies exhibit flatter curves due to geometric saturation effects. Radiative transfer modelling under a uniform dust-star mixture confirms that the observed slope evolution with inclination comes from the galaxy geometry; the slope evolution with stellar mass arises from intrinsic dust property variations, linked to mass-dependent grain processing mechanisms. Additionally, intermediate-scale structures (ISS) at 4870, 6370, and 7690 Å are tentatively detected. These findings underscore the interplay between dust geometry, grain evolution, and galactic environment, offering new insights into dust lifecycle models.

In this study we incorporate a new grid of kilonova simulations produced by the Monte Carlo radiative transfer code SuperNu in an inference pipeline for astrophysical transients, and evaluate their performance. These simulations contain four different two-component ejecta morphology classes. We analyze follow-up observational strategies by Vera Rubin Observatory in optical, and James Webb Space Telescope (JWST) in mid-infrared (MIR). Our analysis suggests that, within these strategies, it is possible to discriminate between different morphologies only when late-time JWST observations in MIR are available. We conclude that follow-ups by the new Vera Rubin Observatory alone are not sufficient to determine ejecta morphology. Additionally, we make comparisons between surrogate models based on radiative transfer simulation grids by SuperNu and POSSIS, by analyzing the historic kilonova AT2017gfo that accompanied the gravitational wave event GW170817. We show that both SuperNu and POSSIS models provide similar fits to photometric observations. Our results show a slight preference for SuperNu models, since the wind ejecta parameters recovered with these models are in better agreement with expectations from numerical simulations.

The study of quasi-periodic oscillations in X-ray binaries provides valuable insights into the physics of accretion around compact objects. The M82 galaxy hosts two ultraluminous X-ray sources (ULXs), one of which is suspected to harbor an intermediate-mass black hole. Using 39 NuSTAR observations acquired between 2014--2024, we investigate the aperiodic X-ray variability in M82. In particular, we study in detail the evolution of the QPO from M82 X-1 in the range 20--300 mHz. We do not find additional timing features in the data, besides a frequent broad noise component at lower frequencies. The QPO behaves similarly to other classes of low-frequency oscillations in accreting compact objects, both black holes and neutron stars.

The Dark matter Nanosatellite Equipped with Skipper Sensors (DarkNESS) deploys a recently developed skipper-CCD architecture with sub-electron readout noise in low Earth orbit (LEO) to investigate potential signatures of dark matter (DM). The mission addresses two interaction channels: electron recoils from strongly interacting sub-GeV DM and X-rays produced through decaying DM. Orbital observations avoid attenuation that limits ground-based measurements, extending sensitivity reach for both channels. The mission proceeds toward launch following laboratory validation of the instrument. A launch opportunity has been secured through Firefly Aerospace's DREAM 2.0 program, awarded to the University of Illinois Urbana-Champaign (UIUC). This will constitute the first use of skipper-CCDs in space and evaluate their suitability for low-noise X-ray and single-photon detection in future space observatories.

The main aim of this study is to examine the behaviour of physical parameters of an anisotropic compact star model demonstrating spherical symmetry in F(Q) modified gravity. To evaluate the behaviour and the stability of an anisotropic compact star model, we utilise the measured mass and radius of an anisotropic compact star model. This study obtained an anisotropic compact star model by solving Einstein field equations. The field equations have been simplified by an appropriate selection of the metric elements and the Karmarkar condition. By solving the field equation to develop a differential equation that establishes a relationship between two essential components of spacetime. A physical analysis of this model reveals that the resulting stellar structure for anisotropic matter distribution is a physically plausible representation of a compact star with an energy density of order $10^14 g/cm^3$. Using the Tolman-Oppenheimer-Volkoff equation, causality condition and Harrison-Zeldovich-Novikov Condition, we investigate the hydrostatic equilibrium and stability of the compact star Cen X-3. We further determined the mass-radius relation of this compact star for different values of delta}1.

We demonstrate that nuclear magnetic resonance based searches for dark matter (DM) have intrinsic and powerful sensitivity to dark photons and the axion-photon coupling. The reason is conceptually straightforward. An instrument such as CASPEr-Gradient begins with a large sample of nuclear spins polarised in a background magnetic field. In the presence of axion DM coupled to nucleons, the spin ensemble feels an effective magnetic field $\mathbf{B} \propto \nabla a$ that tilts the spins, generating a potentially observable precession. If the magnetic field is real rather than effective, the system responds identically. A real field can be generated by a kinetically mixed dark photon within the shielded region the sample is placed or an axion coupled to photons through its interaction with the background magnetic field. We show that all three signals are detectable and distinguishable. If CASPEr-Gradient were to reach the QCD axion prediction of the axion-nucleon coupling, it would simultaneously be sensitive to kinetic mixings of $\epsilon \simeq 3 \times 10^{-16}$ and axion-photon couplings of $g_{a\gamma\gamma} \simeq 2 \times 10^{-16}\,{\rm GeV}^{-1}$ for $m \simeq 1\,\mu{\rm eV}$.

Binary systems of compact objects in close orbit around a supermassive black hole (SMBH) may form in galactic nuclei, providing a unique environment to probe strong-gravity tidal effects on the binary's dynamics. In this work, we investigate precession resonances arising between the periastron precession frequency of a binary system and its orbital frequencies around the SMBH. By modeling the SMBH as a Schwarzschild black hole, we find that relativistic effects in the tidal field give rise to a significantly richer resonance spectrum compared to the Newtonian case. This result is supported by both perturbative and numerical analyses of the quadrupolar tidal interaction in the strong-gravity regime. Our results reveal new signatures for strong-gravity effects in such triple systems, with potential implications for gravitational-wave astronomy.

We study the implications of precision measurements of light element abundances in concordance with the Cosmic Microwave Background for scenarios of physics beyond the Standard Model that generate large inhomogeneities in the baryon-to-photon ratio. We show that precision Big Bang Nucleosynthesis (BBN) imposes strong constraints on any mechanism that produces large scale inhomogeneities at temperatures of the order or below a TeV. In particular, we find that inhomogeneities of the order of $20\%$ at comoving length scales larger than the comoving horizon at the temperature of $3~\mathrm{TeV}$ are in conflict with the measured light element abundances. This sensitivity to physics at such early times is because inhomogeneities in baryon number homogenize predominantly through diffusion, which is a slow process. BBN therefore acts as a novel probe of baryogenesis below the $\mathrm{TeV}$ scale, and readily rules out some of the proposed scenarios of baryogenesis in the literature. We discuss the implications for electroweak baryogenesis. In addition we show that precision BBN is a new probe of first-order phase transitions which produce a gravitational wave signal in the frequency range from pHz to mHz. This leads to constraints on the electroweak phase transition, as well as the first-order phase transitions that have been postulated to explain the pulsar timing array signal. We also discuss the future prospects of improvements in this probe.

Apparent horizon (AH) finders are essential for characterizing black holes and excising their interiors in numerical relativity (NR) simulations. However, open-source AH finders to date are tightly coupled to individual NR codes. We introduce BHaHAHA, the BlackHoles@Home Apparent Horizon Algorithm, the first open-source, infrastructure-agnostic library for AH finding in NR. BHaHAHA implements the first-ever hyperbolic flow-based approach, recasting the elliptic partial differential equation for a marginally outer trapped surface as a damped nonlinear wave equation. To enhance performance, BHaHAHA incorporates a multigrid-inspired refinement strategy, an over-relaxation technique, and OpenMP parallelization. When compared to a naïve hyperbolic relaxation implementation, these enhancements result in 64x speedups for difficult common-horizon finds on a single spacetime slice, enabling BHaHAHA to achieve runtimes within 10% of the widely used (single-core) AHFinderDirect and outperform it on multiple cores. For dynamic horizon tracking with typical core counts on a high-performance-computing cluster, BHaHAHA is approximately 2.1 times faster than AHFinderDirect at accuracies limited by interpolation of metric data from the host NR code. Implemented and tested in both the Einstein Toolkit and BlackHoles@Home, BHaHAHA demonstrates that hyperbolic relaxation can be a robust, versatile, and performant approach for AH finding.

We study a class of geometries in which nonmetricity is fully determined by a vectorial degree of freedom and three independent coefficients. Formulating the simplest linear action in this geometry, implemented through Lagrange multipliers, naturally leads to an equivalence with the purely kinetic $k$-essence models with quadratic kinetic terms. A detailed dynamical systems analysis reveals that the $\Lambda$CDM phenomenology is embedded within the model. Crucially, we find that if stability conditions such as a positive sound speed squared and non-negative energy density are not enforced, the model generically exhibits instabilities and divergent behaviour in the phase space. These physical viability criteria allow us to isolate stable regions of the parameter space and derive well-motivated priors for parameter inference. Using Markov Chain Monte Carlo methods and late-time observational data, including cosmic chronometers, Pantheon$^{+}$ Type Ia supernovae, and DESI baryon acoustic oscillations, we constrain the degrees of freedom associated with nonmetricity and demonstrate the viability of the model. Remarkably, the model is found to be statistically indistinguishable from $\Lambda$CDM at late times. We discuss the implications of these results in light of the recent cosmic tensions, and give a possible explanation as to why the equivalent $k$-essence models have been missed as serious competitors to $\Lambda$CDM in the past. Finally, we review the geometric foundations of the theory and show that the integrable Weyl, Schrödinger and completely symmetric geometries are embedded within our framework as special cases.

Cosmological correlation functions of inflaton and graviton perturbations are the fundamental observables of early universe cosmology and remain a primary target for observations. In this work, we ask the following question: are these observables independent of one another? We find that in the parity-odd sector of inflationary perturbation theory, the answer is a resounding no! In earlier work we derived a correlator-to-correlator factorisation formula which states that parity-odd correlators factorise into lower-point correlators under some mild assumptions on the underlying theory. In this work, we show that these assumptions are satisfied in dynamical Chern-Simons gravity where the action of minimal inflation is augmented by a coupling between the inflaton and the gravitational Chern-Simons term. Such a theory gives rise to a parity-odd trispectrum of curvature perturbations, and we show that such a trispectrum can be expressed solely in terms of the bispectrum that arises due to the minimal coupling between the inflaton and graviton, and the graviton power spectrum which receives a parity-odd correction in this theory. The trispectrum is quadratic in this mixed inflaton-graviton bispectrum and can therefore be interpreted as a ``double copy". Our final expression for the parity-odd trispectrum is a relatively simple function of the external momenta that is rational and factorised.

We investigate the energy spectra of secondary cosmic ray particles associated with two distinct solar events: the magnetospheric effect (ME) of 5 November 2023 and ground-level enhancement (GLE 74) of 11 May 2024. Using data from the SEVAN and Neutron Monitor networks and energy release histograms from particle spectrometers, we reconstruct spectra and identify key differences between ME and GLE. CORSIKA-based simulations reveal that MEs are caused by galactic protons below geomagnetic cutoff rigidities (Rc = 7.1 GV at Aragats) penetrating the magnetosphere during geomagnetic storms, leading to localized flux enhancements at mountain altitudes but not at sea level. In contrast, SEP events initiated by GLEs can involve high-energy solar protons (>10 GeV), producing secondaries that reach sea level at middle latitudes. We present integral energy spectra and spatial correlation of detector responses, demonstrating that SEVAN's energy-resolved data offer new diagnostic tools for identifying hard-spectrum SERs. Our results refine the definition of ME and suggest a strategy for early warning of hazardous solar particle events based on real-time ground-based observations.

In this work, we extend the scalar sector of the conventional hyperchargeless inert triplet model (ITM) to include a second dark matter (DM) candidate, which appears to be a pseudo-Nambu-Goldstone boson (pNGB). The usual ITM with an extended scalar sector offers a DM candidate along with novel signatures at different experiments, e.g., colliders, gravitational wave detectors, etc. Nevertheless, hitherto unseen experimental detections have placed stringent constraints on the ITM parameter space. Moreover, triplet masses lighter than $1.9$ TeV, consistent with the existing or upcoming collider sensitivity reach, are already excluded from the DM observable, as they yield an underabundant relic density due to a strong $SU(2)_L$ gauge annihilation. Inclusion of a pNGB DM, via a complex $SU(2)_L$ scalar singlet and through the soft-breaking of a $U(1)$ symmetry, helps to revive the sub-TeV regime of the triplet DM. This resurgence relies on a proficient conversion between the two DM species. Using this inter-conversion, with the triplet DM as the lighter one between the two, we show that it is possible to push the triplet DM contribution to $50\% - 60\%$ of the total relic density. This offers a significant improvement over the traditional ITM with a single DM candidate, where the same can at most reach $10\% - 20\%$. Besides, the concerned bipartite DM framework also offers the possibility of a first-order phase transition along various constituent field directions. Among these, the one along the real $SU(2)_L$ singlet direction can be a strong one which subsequently yields detectable gravitational wave signals at the upcoming space-based gravitational wave detectors such as LISA, BBO, DECIGO, etc., alongside distinctive and complementary signatures at the various DM and collider quests.

We study in detail the effect of quantum decoherence in neutrino oscillations. We adopt a phenomenological approach that allows us to parametrize the energy dependence of the decoherence effects resulting from the modification of the neutrino dispersion relation. Using the open quantum system framework we derive decoherence parameters, which are usually connected to quantum gravitational effects. Furthermore, we study the sensitivity of decoherence on high-energy astrophysical neutrinos among all possible initial source compositions. We find that variation in the flux composition at neutrino telescopes can be a good probe to test such effects. Additionally, we show that a simple extension with heavy sterile neutrino decoherence produces verifiable signatures.

Extreme mass-ratio inspirals pose a difficult challenge in terms of both search and parameter estimation for upcoming space-based gravitational-wave detectors such as LISA. Their signals are long and of complex morphology, meaning they carry a large amount of information about their source, but are also difficult to search for and analyse. We explore how sequential simulation-based inference methods, specifically truncated marginal neural ratio estimation, could offer solutions to some of the challenges surrounding extreme-mass-ratio inspiral data analysis. We show that this method can efficiently narrow down the volume of the complex 11-dimensional search parameter space by a factor of $10^6-10^7$ and provide 1-dimensional marginal proposal distributions for non-spinning extreme-mass-ratio inspirals. We discuss the current limitations of this approach and place it in the broader context of a global strategy for future space-based gravitational-wave data analysis.

We study the Wheeler-DeWitt quantization of a spatially flat Friedmann-Lemaître-Robertson-Walker (FLRW) universe with pressureless dust (modeled via the Brown-Kuchař formalism) and a dynamical cosmological constant $\Lambda$ treated in the unimodular gravity framework, where unimodular time serves as a relational clock. Remarkably, the quantum dynamics of this system exhibit a mathematical correspondence to a non-relativistic hydrogen atom -- $\Lambda$ maps to energy eigenvalues, the volume variable to the radial coordinate, and the dust energy density parameter to the Coulomb potential strength. This analogy yields a continuous spectrum for positive $\Lambda$, analogous to scattering states. For $\Lambda > 0$, we prove the self-adjointness of the unimodular Hamiltonian, guaranteeing unitary evolution in unimodular time. By constructing wave packets from normalized stationary states, we demonstrate a quantum bounce that resolves the classical Big Bang singularity. The dynamics transition from semiclassical behavior far from the bounce to quantum-dominated regions featuring characteristic "ringing" oscillations due to interference near the bounce. We quantify quantum effects through expectation values and fluctuations of cosmological observables, finding evidence for persistent quantum effects in the late universe. Thus our results suggest that quantum gravitational effects may leave imprints on late-time cosmology, even beyond the bounce.

In this paper, we explore parity violation in a scalar trispectrum from a dynamical Chern-Simons gravity theory. So far, a graviton-mediated diagram with two vertexes being of general relativity has been studied in this theory by taking into account the impact of a modified dispersion relation of gravitons on graviton's bulk propagators. We instead study a parity-odd trispectrum from both a graviton-mediated diagram, where one of the two vertexes originates from the Chern-Simons term, and a contact diagram by using the bulk propagators in general relativity. After computing the scalar-scalar-tensor cubic interactions and the scalar quartic ones originating from the Chern-Simons term, first we show that the resultant parity-odd trispectrum vanishes in the case of Bunch-Davies initial conditions, which is consistent with a no-go theorem for a non-vanishing parity-odd trispectrum. Then, we discuss a way to acquire a non-vanishing parity-odd trispectrum from the viewpoint of non-Bunch-Davies initial conditions.

We perform a Bayesian analysis of the neutron star (NS) equation of state (EoS) based on a wide set of Skyrme functionals, derived from previous nuclear physics inferences. The novelty of this approach lies in starting from the full multidimensional posterior distribution of nuclear matter parameters, consistent with a comprehensive set of static and dynamic nuclear structure observables. We construct unified EoSs for $npe\mu$ matter, where the inner crust of the NS is treated using an extended Thomas-Fermi method, providing for the first time a fully consistent Bayesian treatment of the correlation of bulk with surface as well as with spin-orbit and effective mass parameters. We then employ a standard Bayesian framework to identify those EoSs that satisfy astrophysical constraints from NS mass measurements, the tidal deformability from GW170817, and NICER mass-radius observations. We also examine NS observables, such as the crustal moment of inertia, which is crucial in understanding pulsar glitches. Compared to previous works, we observe an increase in both the NS surface thickness and the crustal moment of inertia.

Collisional flavor instabilities, driven by differing neutrino and antineutrino reaction rates, are expected to occur in dense astrophysical environments like supernovae and neutron star mergers, but have yet to be incorporated in large-scale simulations. We derive, for the first time, analytical expressions for the asymptotic state resulting from a homogeneous and isotropic instability, and apply these predictions to two representative conditions from a neutron star merger simulation. We emphasize the importance of using a collision term that allows for both damping of flavor coherence and relaxation back to thermal equilibrium, which leads to a "quantum" equilibrium with nonzero coherence. These results can be implemented in a subgrid model of collisional flavor transformation, an important step toward the inclusion of flavor oscillation physics into global simulations.