Papers for Monday, Jun 09 2025

Astronomers, and in particular exoplaneteers, have a curious habit of expressing Bayes factors as frequentist sigma values. This is of course completely unnecessary and arguably rather ill-advised. Regardless, the practice is common - especially in the detection claims of chemical species within exoplanet atmospheres. The current canonical conversion strategy stems from a statistics paper from Sellke et al. (2001), who derived an upper bound on the Bayes factor between the test and null hypotheses, as a function of the $p$-value (or number of sigmas, $n_{\sigma}$). A common practice within the exoplanet atmosphere community is to numerically invert this formula, going from a Bayes factor to $n_\sigma$. This goes back to Benneke & Seager (2013) -- a highly cited paper that introduced Bayesian model comparison as a means of inferring the presence of specific chemical species -- in an attempt to calibrate the Bayes factors from their technique for a community that in 2013 was more familiar with frequentist sigma significances. However, as originally noted by Sellke et al. (2001), the conversion only provides an upper limit on $n_\sigma$, with the true value generally being lower. This can result in inflations of claimed detection significances, and this note strongly urges the community to stop converting to $n_\sigma$ at all and simply stick with Bayes factors.

Observations of Supermassive Black Hole (SMBH) at very high redshift hosted by massive galaxy provide a challenge to the SMBH and massive galaxy formation models which are predicated on the Lambda CDM model that predicts a Hubble constant currently in tension with late time observations. In a gravitational system with two asymptotic solutions such as a central mass in an expanding cosmological background or an expanding background with radiation and baryons, we find that a particle's free falling speed resulting from adding the free falling speeds due to the asymptotic solutions is also a theoretically viable solution [1]-[2]. In this work, we use this new model to examine the SMBH formation time-line and find that the change in the redshift-cosmological time relation provides sufficent cosmological time for the formation of SMBHs at high redshifts. We also find that the Hubble parameter-redshift relation in this model is consistent with late time observation at z<1.

The ram-pressure acceleration of cold gas by hot outflows plays a crucial role in the dynamics of multiphase galactic winds. Recent numerical studies incorporating radiative cooling have identified a size threshold for idealized cold clouds to survive within high-velocity outflows. This study extends the investigation to a more complex morphology of cold gas as observed in the interstellar medium. We conduct three-dimensional hydrodynamic simulations of ensembles of individual spherical clouds to systematically explore under which conditions the cold clouds can survive. We find that cloud ensembles can survive collectively -- even when individual clouds, if isolated, would be rapidly destroyed. Our results indicate that, besides the morphology, factors such as tight packing, small inter-cloud distance and higher fragmentation facilitate survival. We propose a novel multi-cloud survival criterion that accounts for collective properties of the cloud system, including total gas mass and the geometric configuration based on an effective volume filling fraction of the cold gas $F_V$. This fraction is computed by constructing a composite volume from individual enclosing conical boxes aligned with the wind, incorporating spatial overlap and cloud-tail spreading. The box dimensions scale with the critical survival radius $r_{\rm crit}$ from the single-cloud criterion. We find a universal threshold $F_{V,{\rm crit}}\approx 0.24$ that robustly separates surviving from destroyed systems across diverse geometric configurations. Our findings emphasize the critical importance of initial cloud distribution and fragmentation in governing the long-term evolution and survival of cold gas structures, providing insight into observed multiphase outflows and CGM dynamics.

W49B is a unique Galactic supernova remnant with centrally peaked, "bar"-like ejecta distribution, which was once considered evidence for a hypernova origin that resulted in a bipolar ejection of the stellar core. However, chemical abundance measurements contradict this interpretation. Closely connected to the morphology of the ejecta is its velocity distribution, which provides critical details for understanding the explosion mechanism. We report the first-ever observational constraint on the kinematics of the ejecta in W49B using the Resolve microcalorimeter spectrometer on the X-ray Imaging and Spectroscopy Mission (XRISM). Using XRISM/Resolve, we measured the line-of-sight velocity traced by the Fe He$\alpha$ emission, which is the brightest feature in the Resolve spectrum, to vary by $\pm$300 km s$^{-1}$ with a smooth east-to-west gradient of a few tens of km s$^{-1}$ pc$^{-1}$ along the major axis. Similar trends in the line-of-sight velocity structure were found for other Fe-group elements Cr and Mn, traced by the He$\alpha$ emission, and also for intermediate-mass elements Si, S, Ar, and Ca, traced by the Ly$\alpha$ emission. The discovery of the east-west gradient in the line-of-sight velocity, together with the absence of a twin-peaked line profile or enhanced broadening in the central region, clearly rejects the equatorially expanding disk model. In contrast, the observed velocity structure suggests bipolar flows reminiscent of a bipolar explosion scenario. An alternative scenario would be a collimation of the ejecta by an elongated cavity sculpted by bipolar stellar winds.

The abundance and mass distribution of galaxy clusters is a sensitive probe of cosmological parameters, through the sensitivity of the high-mass end of the halo mass function to $\Omega_m$ and $\sigma_8$. While galaxy cluster surveys have been used as cosmological probes for more than a decade, the accuracy of cluster count experiments is still hampered by systematic, such as the relation between observables and halo mass, the accuracy of the halo mass function, and the survey selection function. Here we show that these uncertainties can be alleviated by forward modeling the observed cluster population with simulation-based inference. We construct a pipeline that predicts the distribution of observables from cosmological parameters and scaling relations, and then train a neural network to learn the mapping between the input parameters and the measured distributions. We focus on fiducial X-ray surveys with available flux, temperature, and redshift measurements, although the method can be easily adapted to any available observable. We apply our method to mock samples extracted from the UNIT1i simulation and demonstrate the accuracy of our approach. We then study the impact of several systematic uncertainties on the recovered cosmological parameters. We show that sample variance and the choice of the halo mass function are subdominant sources of uncertainty. Conversely, the absolute mass scale is the leading source of systematic error and must be calibrated at the $<10\%$ level to recover accurate values of $\Omega_m$ and $\sigma_8$. However, the quantity $S_8=\sigma_8(\Omega_m/0.3)^{0.3}$ appears to be less sensitive to the accuracy of the mass calibration. We conclude that simulation-based inference is a promising avenue for future cosmological studies from galaxy cluster surveys such as eROSITA and Euclid as it allows to consider all the available observables in a straightforward manner.

We present a Field of Streams visualizing stellar structures in the Milky Way halo as viewed by the Dark Energy Camera (DECam). We use $g$- and $r$-band imaging from the Dark Energy Survey and the DECam Legacy Survey, covering $18{,}700 °^2$ across the sky. Using an isochrone-based matched filter in $g$ vs. $g-r$, we select old and metal-poor stars in three distance bins, and generate a false-color RGB image of the number density of selected stars. The DECam Field of Streams shows a variety of Milky Way halo structures, including dwarf galaxies, globular clusters, and an abundance of stellar streams, illustrating the significant progress that has been made in recent years in uncovering the building blocks of the Milky Way's stellar halo in deep, wide-area photometric surveys. This view of our Galaxy will be improved in the coming years as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) begins to collect data to greater depths and across a larger fraction of the sky than ever before.

Using deep JWST/NIRSpec spectra from the Blue Jay survey, we perform the first systematic investigation of neutral gas content in massive galaxies at Cosmic Noon based on the Ca II H, K absorption lines. We analyze a sample of 9 galaxies at 1.8 < z < 2.8 with stellar masses > 10.6, for which we detect neutral gas absorption both in Ca II and in Na I. After removing the stellar continuum using the best-fit model obtained with Prospector, we fit the excess absorption due to neutral gas in the Ca II H, K doublet and in the Na I D doublet, together with nearby emission lines produced by ionized gas. We measure covering fractions between 0.2 and 0.9 from the Ca II H and K lines, which are spectrally well resolved in the NIRSpec R ~ 1000 observations, unlike the absorption lines in the Na I D doublet. We measure the velocity shift, velocity dispersion, and column density separately for Ca II and Na I. About half of the galaxies present blueshifted Ca II, indicative of an outflow of neutral gas, consistent with previous results based on Na I. The velocity shift and the column density measured from Ca II are correlated with those measured from Na I, implying that these absorption lines trace gas in similar physical conditions. However, the column densities are not in a 1:1 relation, meaning that the relative amount of Ca II and Na I atoms along the line of sight varies with the gas column density. After discussing possible reasons for this behavior, we derive an empirical relation between the column density of Ca II and the column density of Na I and, in a more indirect way, of neutral hydrogen H I. This calibration offers a new way to estimate the outflow mass and the mass outflow rate for the neutral phase from current and future JWST observations of massive galaxies at Cosmic Noon and beyond

Type Ia Supernovae (SNe Ia) are crucial tools to measure the accelerating expansion of the universe, comprising thousands of SNe across multiple telescopes. Accurate measurements of cosmological parameters with SNe Ia require a robust understanding and cross-calibration of the telescopes and filters. A previous cross-calibration effort, 'Fragilistic', provided 25 photometric systems, but offered no public code or ability to add new surveys. We provide an open-source cross-calibration solution, available at this https URL . Using the Pan-STARRs (PS1) and Gaia all-sky telescopes, we characterise the measured filters from 11 photometric systems, including CfA, PS1, Foundation, DES, CSP, SDSS, and SNLS, using published observations of field stars. For the first time, we derive uncertainties on effective filter transmissions and modify filters to match the data. With the addition of direct observations of DA white dwarfs (Boyd et al. 2025), we simultaneously cross-calibrate our zeropoints across photometric systems and propagate to cosmology. With improved uncertainties from DA WDs, we find improvements to the calibration systematic uncertainty of x1.5 for the Pantheon+ (Brout et al. 2022) sample with a new systematic photometric uncertainty = 0.016 for FlatwCDM, and modest improvements to that of the DES5YR analysis. We find good agreement with previous calibration, and show that even these small calibration changes can be amplified by up to a factor of x6 in the inferred SN Ia distances, driven by calibration sensitivity in the colour-luminosity relations and SALT training. Initial results indicate that these changes cause dmu/dz = 0.025 and change the recovered value of Omega_M in LCDM by ~0.01. These may have a potentially larger impact in w0/wa space and inferences about evolving dark energy. We pursue this calculation in an ongoing full re-analysis of DES.

Jellyfish galaxies are extreme examples of how galaxies can transform due to dense environmental effects. These satellite galaxies suffer from ram-pressure stripping, leading to the formation of their distinctive gaseous tails. Some recent observational studies find that jellyfish galaxies are more likely to host active galactic nuclei (AGN) compared to central galaxies of the same mass, suggesting a link between ram pressure and supermassive black hole (SMBH) accretion. We use the IllustrisTNG cosmological-magnetohydrodynamical simulations, namely TNG50 and TNG100, to explore the presence of AGN in jellyfish galaxies with $M_{\rm{stellar}}\simeq10^{9.5-10.8}\,\rm{M}_\odot$ at redshift $z=0$ from the Zooniverse "Cosmological Jellyfish" citizen-science project. Compared to central galaxies, jellyfish are more likely to host an AGN ($L_{\rm AGN}\geq10^{44}\,\mathrm{erg\,s^{-1}}$) particularly at high stellar masses ($M_{\rm stellar}\gtrsim10^{10}\,\mathrm{M_\odot}$). Jellyfish are also more likely to host an AGN than satellites of the same mass, largely because many satellite galaxies are gas-poor and therefore have lower SMBH accretion rates. Compared to non-jellyfish satellites with similar gas content, jellyfish typically undergo stronger ram pressure and have higher central gas densities along with lower central gas sound speeds, although these effects are smaller at lower stellar masses ($M_{\rm stellar}\lesssim10^{10}\,\mathrm{M_\odot}$). Together with case studies of individual galaxies, our population analysis indicates that ram pressure can play a key role in fuelling AGN activity in a large fraction of jellyfish, where gas compression can lead to intense episodes of AGN feedback and star formation. Thus, it is essential to consider both environmental and secular processes for a more complete picture of satellite galaxy evolution.

Tidal disruption events (TDEs) have been proposed as valuable laboratories for studying dormant black holes. However, progress in this field has been hampered by the limited number of observed events. In this work, we present TDECat, a comprehensive catalogue of 134 confirmed TDEs (131 optical TDEs and 3 jetted TDEs) discovered up to the end of 2024, accompanied by multi-wavelength photometry (X-ray, UV, optical, and IR) and publicly available spectra. We also study the statistical properties, spectral classifications, and multi-band variability of these events. Using a Bayesian Blocks algorithm, we determine the duration, rise time , decay time, and their ratio for 103 flares in our sample. We find that these timescales follow a log-normal distribution. Furthermore, our spectral analysis shows that most optical TDEs belong to the TDE-H+He class, followed by the TDE-H, TDE-He, and TDE-featureless classes, which is consistent with expectations from main sequence star disruption. Using archival observations, we identify four new potentially repeating TDEs, namely AT 2024pvu, AT 2022exr, AT2021uvz, and AT 2019teq, increasing the number of known repeating events. In both newly identified and previously known cases, the secondary flares exhibit a similar shape to the primary. We also examine the IR and X-ray emission from the TDEs in our catalogue, and find that 14 out of the 18 IR events have associated X-ray emission, strongly suggesting a potential correlation. Finally, we find that for three subsamples (repeating flares, IR-, and X-ray-emitting events), the spectral classes are unlikely to be randomly distributed, suggesting a connection between spectral characteristics and multi-wavelength emission. TDEcat enables large-scale population studies across wavelengths and spectral classes, providing essential tools for navigating the data-rich era of upcoming surveys such as the LSST.

The hydrodynamical cosmological simulation \eagle{} is used to model the Halo Occupation Distribution (HOD) of moderate X-ray luminosity active galactic nuclei (mXAGN), extending previous work using only dark matter simulations and empirical relations. By examining mergers as a triggering mechanism, we focus on halos typical of galaxy groups and cluster-like systems with masses $\geq 10^{12.75}\,{\rm M}{\odot}\,h^{-1}$. We analyze simulation data to create catalogs of central and satellite galaxies. We study their merger history we quantify the percentage of minor and major mergers in the mXAGN sample. We obtain the HOD for central and satellite mXAGN across a redshift interval from \(z=2\) to the present epoch. Our results indicate that, across most redshifts, minor mergers slightly predominate as the primary mechanism for triggering mXAGN.

Quasi-Periodic Eruptions (QPEs) are recurrent X-ray eruptions discovered so far in the nuclei of low-mass galaxies. However, despite considerable observational progress, the origin of QPEs remains unclear. A variety of models have been proposed to explain their nature, but a definitive understanding has yet to be reached. Recently, chaotic mixtures of multiple overlapping eruptions with varying amplitudes have been observed in eRO-QPE1 obs1-features not reported in any other known QPE sources. This complex behavior presents a challenge to the existing QPE models. In this paper, we propose that the overlapping features may be the result of gravitational lensing. We analyze the light curve of eRO-QPE1 and compare its features to predictions from gravitational lensing scenarios. We discuss the implications for the trigger mechanism of QPEs in general. We show that the unique overlapping features observed in eRO-QPE1 may be naturally reproduced by gravitational lensing effects, without invoking a different physical origin from other known QPE sources.

In the identification of new planetary candidates in transit surveys, the employment of Deep Learning models proved to be essential to efficiently analyse a continuously growing volume of photometric observations. To further improve the robustness of these models, it is necessary to exploit the complementarity of data collected from different transit surveys such as NASA's Kepler, Transiting Exoplanet Survey Satellite (TESS), and, in the near future, the ESA PLAnetary Transits and Oscillation of stars (PLATO) mission. In this work, we present a Deep Learning model, named DART-Vetter, able to distinguish planetary candidates (PC) from false positives signals (NPC) detected by any potential transiting survey. DART-Vetter is a Convolutional Neural Network that processes only the light curves folded on the period of the relative signal, featuring a simpler and more compact architecture with respect to other triaging and/or vetting models available in the literature. We trained and tested DART-Vetter on several dataset of publicly available and homogeneously labelled TESS and Kepler light curves in order to prove the effectiveness of our model. Despite its simplicity, DART-Vetter achieves highly competitive triaging performance, with a recall rate of 91% on an ensemble of TESS and Kepler data, when compared to Exominer and Astronet-Triage. Its compact, open source and easy to replicate architecture makes DART-Vetter a particularly useful tool for automatizing triaging procedures or assisting human vetters, showing a discrete generalization on TCEs with Multiple Event Statistic (MES) > 20 and orbital period < 50 days.

For more than fifty years, astronomers have mapped the neutral hydrogen gas in the Galaxy assuming kinematically derived distances. We employ the distances of nearby young stars, which trace the gas from which they formed, in longitude-latitude-velocity space to map this gas without using kinematic distances. We denote this new method "pattern matching". Analysis of simulated spiral galaxies indicates that our pattern matching distances are 24% more accurate than kinematic distances for gas within 15 kpc of the Sun. The two methods provide similar agreement with parallaxes towards these masers, although the kinematic method shows a small systematic offset in the distance that is not present in the pattern matching distance. Using parallaxes and velocities for masers, we show that this novel method, when matched with nearby Cepheids, performs well compared to kinematics. This analysis is restricted to sources that have a reasonably good match with a member of our Cepheid sample. The distances derived here, and the associated map, have broad utility - from improving our understanding of star formation and the dynamical structure of the Galaxy, to informing 3-D dust maps.

The presence of dust in galaxies at redshifts $z>5$ is commonly connected with core collapse supernovae (SN). Galaxies with exceptionally large dust mass, of order of $1 - 3$\% of the galaxy stellar mass, have been detected during the last decade. The required SN dust yield is $\gtrsim 1~M_\odot$ per supernova, which is comparable to the theoretically predicted maximum. However, the reverse shock (RS) penetrating the SN ejecta significantly destroy the dust particles nucleating there through sputtering. The resulting net dust mass injected into the interstellar gas after processing by the RS turns out to be $\lesssim 0.1~M_\odot$ per SN. This makes the explanation of the existence of $z>5$ galaxies with dust masses as high as $M_d\gtrsim (0.01-0.03)M_\ast$ a challenging one. In this paper we present arguments in favor of an efficient inhibition of the sputtering behind the RS, caused by a strong radiation cooling from the dust in the supernova ejecta.

The quasar correlation function assesses the occurrence of quasar pairs as a function of separation, which is strongly influenced by quasar host halo masses. The empirical Trinity model recently inferred the redshift-dependent relationship between supermassive black hole (SMBH) mass, galaxy mass, and halo mass, using constraints other than correlation functions (e.g., quasar luminosity functions, active galactic nuclei occupation fractions, and SMBH mass-bulge mass relations). Hence, comparing the predicted quasar correlation functions from Trinity to real observations is an important test of Trinity's inferred SMBH -- halo relation. In this work, we use a compilation of observed two-point projected and redshift-space correlation functions from $0 \leq z \leq 3.5$. We find that Trinity accurately predicts quasar correlation functions within observed error bars, although observations do not have much constraining power at lower redshifts due to smaller observable volumes and lower quasar number densities. This finding is consistent with Trinity having the correct placement of quasars within their host galaxies and dark matter halos, without requiring quasar clustering constraints during model fitting. Using Trinity, we also predict the clustering as a function of quasar bolometric luminosity, finding that existing survey uncertainties are too large to show measurable differences ($\lesssim 0.3$ dex change in bias for $10^{42}$ erg s$^{-1}$ compared to $10^{46}$ erg s$^{-1}$ SMBHs across redshifts). This fact arises because most SMBH growth (and hence quasar luminosity) occurs in halos in a similar mass range ($10^{12}-10^{13} M_\odot$).

Helium-core white dwarfs (He WDs) formed through common envelope (CE) evolution offer valuable insight into binary interaction channels and compact remnant formation. Their cooling rates critically impact both detectability and age estimates in close binaries. Compared to He WDs formed via stable Roche-lobe overflow (SRLOF), those from the CE channel undergo markedly different mass-loss histories, resulting in distinct post-CE evolutionary behavior. We explore how the mass of the residual hydrogen envelope (Mh) shapes the cooling evolution of CE He WDs, focusing on the role of the bifurcation point in setting Mh and enabling residual hydrogen burning. Using the LPCODE stellar evolution code, we computed models of He WDs with masses from 0.20 to 0.42 solar masses, evolving from post-CE conditions to the white dwarf cooling track. Two evolutionary branches emerge: (i) non-flashing sequences, which cool rapidly with negligible hydrogen burning, and (ii) flashing sequences, where hydrogen shell flashes alter the envelope structure prior to cooling. Minimal-envelope models cool within 5-130 million years for effective temperatures between 12,000 and 27,000 K, and reach in approx. 300 million years at lower temperatures, remaining much younger than SRLOF counterparts. In contrast, models with more hydrogen retain active nuclear burning, delaying cooling and yielding ages of several billion years. Flashing sequences prolong the pre-white dwarf phase, though still shorter than in SRLOF evolution. The value of Mh also affects WD mass and surface gravity estimates, introducing systematic shifts with respect to SRLOF WDs. Our results show that CE He WDs follow distinct evolutionary paths, with important implications for interpreting the nature and fate of compact binaries hosting He WDs.

We present the X-ray polarization observation of G21.5-0.9, a young Galactic supernova remnant (SNR), conducted with the Imaging X-ray Polarimetry Explorer (IXPE) in October 2023, with a total livetime of approximately 837 ks. Using different analysis methods, such as a space-integrated study of the entire region of the PWN and a space-resolved polarization map, we detect significant polarization from the pulsar wind nebula (PWN) at the center of the SNR, with an average polarization degree of ~10% oriented at ~33° (north through east). No significant energy-dependent variation in polarization is observed across the IXPE band (2-8 keV). The polarization map, corrected for the effect of polarization leakage, reveals a consistent pattern in both degree and angle, with little change across the nebula. Our findings indicate the presence of a highly polarized central torus, suggesting low levels of turbulence at particle acceleration sites. Unlike Vela, but similar to the Crab Nebula, we observe substantial differences between radio and X-ray polarization maps. This suggests a clear separation in energy of the emitting particle populations and hints at an important, yet poorly understood, role of instabilities in the turbulence dynamics of PWNe.

The Transiting Exoplanet Survey Satellite (TESS) has surveyed nearly the entire sky in Full-Frame Image mode with a time resolution of 200 seconds to 30 minutes and a temporal baseline of at least 27 days. In addition to the primary goal of discovering new exoplanets, TESS is exceptionally capable at detecting variable stars, and in particular short-period eclipsing binaries which are relatively common, making up a few percent of all stars, and represent powerful astrophysical laboratories for deep investigations of stellar formation and evolution. We combed Sectors 1-82 of TESS Full-Frame Image data searching for eclipsing binary stars using a neural network that identified ~1.2 million stars with eclipse-like features. Of these, we have performed an in-depth analysis on ~60,000 targets using automated methods and manual inspection by citizen scientists. Here we present a catalog of 10001 uniformly-vetted and -validated eclipsing binary stars that passed all our ephemeris and photocenter tests, as well as complementary visual inspection. Of these, 7936 are new eclipsing binaries while the remaining 2065 are known systems for which we update the published ephemerides. We outline the detection and analysis of the targets, discuss the properties of the sample, and highlight potentially interesting systems. Finally, we also provide a list of ~900,000 unvetted and unvalidated targets for which the neural network found eclipse-like features with a score higher than 0.9, and for which there are no known eclipsing binaries within a sky-projected separation of a TESS pixel (~21 arcsec).

Population synthesis simulations of compact binary coalescences~(CBCs) play a crucial role in extracting astrophysical insights from an ensemble of gravitational wave~(GW) observations. However, realistic simulations are costly to implement for a dense grid of initial conditions. Normalizing flows can emulate the distribution functions of a simulated population of binary parameters and thereby enable empirical constraints on the astrophysical initial conditions and branching fractions of various formation channels given data from a catalog of GW observations. They can also be used for data amplification in sparse regions of the CBC parameter space to guide the development of phenomenological population models for rarely synthesizable systems with components in theorized mass gaps, without having to simulate a prohibitively large number of binaries. But flow predictions are wrought with uncertainties, especially for sparse training sets. In this work I develop a method for quantifying and marginalizing uncertainties in the emulators by introducing the Bayesian Normalizing flow, a conditional density estimator constructed from Bayesian neural networks. Using the exact likelihood function associated with density estimators I sample the posterior distribution of flow parameters with suitably chosen priors to quantify and marginalize over flow uncertainties. I demonstrate the accuracy, calibration, and data-amplification impacts of the estimated uncertainties for simulations of binary black hole populations formed through common envelope evolution. I outline applications of the methodology in simulation-based inference from growing GW catalogs and sketch other uses for general simulation-based approaches in GW astronomy.

We investigate the interaction between two filaments (F1 and F2) and their subsequent magnetic reconnection in active region (AR) NOAA 13296 and AR NOAA 13293 on May 9, 2023, utilizing high spatial and temporal resolution and multi-wavelength observational data from the Solar Dynamics Observatory, the New Vacuum Solar Telescope, and the Chinese H{\alpha} Solar Explorer. The movement of F1 from the southeast toward the northwest, driven by the motion of the positive magnetic polarity (P1), leads to a collision and reconnection with F2. This reconnection exchanges their footpoints, resulting in the formation of two new filaments (F3 and F4) consistent with "slingshot" type filament interaction. During the interaction, the current sheet moving due to the motion of F1 and the reconnection outflows moving along F3 and F4 were both observed. The current sheet is rarely observed in the slingshot type filament interaction, measuring approximately 2.17 Mm in length and 0.84 Mm in width. After the interaction, the F1 disappears whereas a portion of F2 remains, indicating that the interaction involves partial slingshot reconnection, due to the unequal magnetic flux between the filaments. The residual part of F2 will undergo another magnetic reconnection in the same interaction region with the magnetic loops connecting polarities N1 and P1. The material generated by the reconnection is continuously injected into F4, leading to its final morphology. The findings enhance our understanding of slingshot-type filament interactions, indicating that partial slingshot reconnections between filaments may be more common than full slingshot events.

We present a comparative observational data analysis of two DBI-type k-essence scalar field models for late-time cosmic acceleration: one treats dark energy with standard matter, while the other unifies dark energy and dark matter within a single scalar field, each governed by distinct non-canonical Lagrangians. The background dynamics are formulated via modified Friedmann and scalar field equations, constrained by Pantheon+ SNe Ia, Hubble, and BAO data, using Bayesian inference (NUTS in NumPyro) and a neural network-based emulator for efficiency. Introducing a nuisance parameter $\mu_0$ improves the robustness and interpretability of key cosmological parameters, particularly $H_0$, $r_d$, and $\Omega_{d0}$ by absorbing residual systematics, thereby mitigating the Hubble tension. The deceleration parameter $q(z)$ is computed for both models, revealing key differences in cosmic acceleration history: Model I yields a present value $q^{Bayes}_0=-0.589$ or $q^{ML}_0=-0.559$ with transition redshift $z^{Bayes}_{trans}=0.740$ or $z^{ML}_{trans}=0.686$, closely aligning with $\Lambda$CDM predictions, while Model II shows a stronger present acceleration $q^{Bayes}_0=-0.893$ or $q^{ML}_0=-0.819$ but with a later transition at $z^{Bayes}_{trans}=0.605$ or $z^{ML}_{trans}=0.564$, indicating a steeper late time evolution. Based on $\chi^2$, $\chi^2_\nu$, AIC, and BIC statics, Model II performs worse than Model I despite its conceptual appeal. Finally, Symbolic regression (GPlearn) is performed to reconstructs $\omega_{eff}(z)$. Our analysis shows that Model I is statistically favored across all metrics and offers a physically consistent alternative to $\Lambda$CDM, capable of addressing the Hubble tension.

We present two NuSTAR observations of the X-ray transient, MAXI J1752$-$457, following a superburst which was observed by MAXI/GSC in November, 2024. NuSTAR follow-up confirmed that MAXI J1752$-$457 is coincident with the previously observed Einstein Probe source, EP240809a. We performed a spectral analysis of the source during both NuSTAR observations, and we confirm that the hard X-ray spectra are consistent with the inclusion of a spherical blackbody component. At about 79 hours after the onset of the superburst, we find a blackbody temperature of $kT_\mathrm{bb}=0.60\pm0.1$ keV and $R_\mathrm{bb}/D_{8}=6.0^{+0.4}_{-0.3}$ km (not including corrections for scattering in the neutron star atmosphere), where $D_{8}$ is the source distance, which is not yet known, in units of 8 kpc. Furthermore, we report a hard X-ray excess which is fit well by a power law with photon index $\Gamma\approx4$, much steeper than those typically seen during accretion onto neutron stars at similar luminosities. We infer that the electron energy distribution in the Comptonizing medium which gives rise to the power law component differs significantly compared to purely accretion-powered neutron star transients even several days after the onset and rapid decay of superbursts. We also performed an energy-resolved timing analysis which showed that the source variability was dominated by red noise in the Comptonization component, suggesting coupling with an accretion disk, rather than being seeded by thermal emission from the neutron star surface.

We search for Lorentz invariance violation (LIV) using the spectral lag data for GRB 190114C using frequentist analysis, where we deal with the astrophysical nuisance parameters using profile likelihood. For this use case, we find a global minima for the $\chi^2$ as a function of energy scale of LIV ($E_{QG}$), well below the Planck scale. The best-fit $2\sigma$ central intervals for $E_{QG}$ are given by $2.81^{+0.96}_{-0.72}\times 10^{14}$ GeV and $9.10^{+2.59}_{-0.64}\times 10^{5}$ for linear and quadratic LIV, respectively and agree with the Bayesian estimates hitherto obtained in a previous work. Therefore, the results from frequentist analysis GRB 190114C agrees with Bayesian analysis and presents yet another proof of principles applications of profile likelihood in the search for LIV.

Depending on the astrophysical source and its environment, the accretion flows can exhibit a variety of behaviors and characteristics in accordance with the type of solutions. We study low-angular-momentum accretion flows onto black holes using two-dimensional general relativistic hydrodynamic (GRHD) simulations to find imprints of different types of accretion solutions. Such flows, relevant to X-ray binaries and wind-fed low-luminosity active galactic nuclei, often lack sufficient angular momentum to form standard accretion disks. We initialize simulations with semi-analytical transonic solutions defined by specific energy (${\cal E}_0$) and angular momentum ($\lambda_0$), allowing a systematic classification of flow types with: (i) an outer sonic point, (ii) an inner sonic point, and (iii) both, exhibiting shock-like transitions. Only solutions with two sonic points produce hot, thermally driven bipolar jets/outflows with Lorentz factors up to $\gamma\sim2$, despite the absence of magnetic fields. Using a general relativistic radiation transfer calculation, we compute broadband spectra and images at X-ray ($1 \, \rm keV$) from bremsstrahlung emission. Radiative properties depend strongly on the type of accretion solution. Solutions with inner sonic points produce the brightest and most extended X-ray emission, while outer-point solutions produce compact, fainter signals. These multidimensional models are thus essential for predicting radiative signatures and will enable the development of semi-analytical tools for interpreting X-ray binaries and possibly Sgr~A$^*$ in weak magnetic field regimes.

We report the timing results with Insight-HXMT observations of X-ray binary IGR J19294+1816 during its 2019 Type I outburst at the decline phase shortly following its peak. We analyze the light curves and power density spectrum (PDS) of the 2019 observations and reveal a peak at approximately $\nu_{NS} \sim 80.2$ mHz, corresponding to X-ray pulsations from the neutron star. In addition, a significant quasi-periodic oscillation (QPO) feature is observed at around $\nu_{QPO} \sim 30.2$ mHz from 10-50 keV, with the rms amplitude increasing with energy. Furthermore, we detect two QPOs at the frequency of $\sim 51.1$ mHz and $113.7$ mHz (corresponding to sidebands near $\nu_{NS} \pm \nu_{QPO}$) in 25-50 keV, exhibiting an rms amplitude of around 12$\%$. Wavelet analysis also shows multiple QPOs at the frequency of $\sim 30$ mHz, $50$ mHz and $ 110$ mHz and these QPO features show transient behaviors, the centroid frequencies of $\sim 30$ mHz remain nearly constant for different luminosities. Our research identifies IGR J19294+1816 as the second strong magnetic-field pulsar with significant sideband signals around the spin frequency. We explore various physical origins that could explain the presence of multiple QPOs.

Aims. The magnetized medium induces birefringence, splitting the light into two distinct wave modes. The differing propagation speeds of the two modes result in different trajectories. Strong gravitational lensing amplifies the birefringence and introduces an additional geometric rotation on top of the Faraday rotation. We compare the geometric rotation with the Faraday rotation. Methods. We construct the lens equation for massive objects in a magnetized plasma environment, and calculate the time delay difference between the two modes using two toy examples. We present that in the strong lensed radio sources, birefringence causes geometric rotation, which is a non-negligible effect, even with a weak magnetic field. Results. In both examples, the geometric delay causes a comparable or stronger rotation than the Faraday rotation and show a similar dependence on the wavelength of the signal. For a point lens with a strong magnetic field, the two wave modes exhibit distinct behaviours. The polarization of lensed sources can provide additional insights into the magnetic field and plasma environment.

Axions are an elegant solution to the strong CP problem for particle physics and a promising dark matter candidate. They can convert into photons under a strong magnetic field, while magnetars with extreme magnetic fields are natural labs for axion detection. Radio telescopes can detect the radio emission from axion-photon conversion near magnetars. In this study, we have refined the calculation of axion-photon conversion and developed the matched filtering integration method to largely improve the signal-to-noise ratio. We validate our method using end-to-end simulation and real observational data from TMRT. A new constraint is set with only 687 seconds of observations with TMRT. Using 10 hours of observation with the high-frequency receiver in FAST or SKA, we can reach the theoretical coupling constant prediction for the axion mass range from 1$\mu$eV to 100$\mu$eV. We validate the possibility of axion detection with radio telescopes and avoid spectrum confusion.

A statistical analysis of anisotropic quasiperiodic features of the spatial distribution of galaxy clusters obtained on the basis of spectroscopic and photometric redshifts in the interval $0.1 \leq z \leq 0.47$ has been carried out. Based on data from the SDSS~III catalog, we show that the preferential direction previously detected in the northern hemisphere (a narrow cone of directions: $\alpha_0=170^\circ \pm 5^\circ, \ \delta_0= 29^\circ \pm 5^\circ$), along which the one-dimensional distribution of projections of the Cartesian coordinates of clusters contains a significant ($\gtrsim (4 - 5) \sigma$) quasi-periodic component, can also be found using photometric redshifts, achieving a certain accuracy ($\Delta z \lesssim 0.013$). Based on data from the photometric DES$\times$unWISE catalog, we have analyzed the spatial distribution of clusters in the southern hemisphere, where a cone of close directions was also detected ($\alpha_0=346^\circ \pm 5^\circ,\ \delta_0= - 29^\circ \pm 5^\circ $), which are approximately an extention of the directions revealed in the northern hemisphere. The power spectra of one-dimensional distributions along these directions contain significant ($\gtrsim (4 - 5) \sigma$) features in the same interval of wave numbers $0.04 \lesssim k \lesssim 0.06~h$~Mpc$^{-1}$.

The KM3NeT research infrastructure comprises two neutrino telescopes located in the deep waters of the Mediterranean Sea, namely ORCA and ARCA. KM3NeT/ORCA is designed for the measurement of neutrino properties and KM3NeT/ARCA for the detection of high\nobreakdashes-energy neutrinos from the cosmos. Neutrinos are indirectly detected using three\nobreakdashes-dimensional arrays of photo\nobreakdashes-sensors which detect the Cherenkov light that is produced when relativistic charged particles emerge from a neutrino interaction. The analogue pulses from the photo\nobreakdashes-sensors are digitised offshore and all digital data are sent to a station on shore where they are processed in real time using a farm of commodity servers and custom software. In this paper, the design and performance of the software that is used to filter the data are presented. The performance of the data filter is evaluated in terms of its purity, capacity and efficiency. The purity is measured by a comparison of the event rate caused by muons produced by cosmic ray interactions in the Earth's atmosphere with the event rate caused by the background from decays of radioactive elements in the sea water and bioluminescence. The capacity is measured by the minimal number of servers that is needed to sustain the rate of incoming data. The efficiency is measured by the effective volumes of the sensor arrays.

The architectures of exoplanet systems are likely set during the initial planet-formation phase in the circumstellar disk. To understand this process, we have to study the earliest phases of planet formation. Complex sub-structures, believed to be driven by embedded planets, have been detected in a significant portion of disks observed at high angular resolution. We aim to extend the sample of such disks to low stellar masses and to connect the disk morphology to the expected proto-planet properties. We resolve the disk in the 2MASSJ16120668-3010270 system for the first time in scattered near-infrared light on scales of 10 au using VLT/SPHERE and reveal an exceptionally structured disk. We find an inner disk (inside 40 au) with two spiral arms, separated by a gap from an outer ring. By comparison with hydrodynamic models, we find that these structures are consistent with the presence of an embedded gas giant with a mass range between 0.1 and 5 MJup depending on the employed model. Our SPHERE observations find a tentative candidate point source within the disk gap, which may be consistent with this mass range if it indeed traces thermal emission by an embedded planet. This interpretation is somewhat strengthened by the proximity of this signal to compact mm continuum emission in the disk gap, which may trace circumplanetary material. It is, however, unclear if this tentative companion candidate could be responsible for the observed disk gap size, given its close proximity to the inner disk. The 2MASSJ16120668-3010270 system is one of only a few systems that shows this exceptional morphology of spiral arms located inside a scattered light gap and ring. We speculate that this may have to do with a higher disk viscosity compared with other systems such as PDS 70.

We analyze a wide set of historical magnetar burst observations detected with five different instruments, calibrating these to the energy range of Fermi-GBM observations for consistency. We find a striking correlation between a magnetar's characteristic age and both its typical burst energy and its burst activity level. Arguing that this bursting behaviour also correlates with true age, we interpret it as the result of a reducing high-stress volume of the crust in an aging magnetar: previous giant flares cause relaxation of large regions of its crust and inhibit burst clustering, whilst the reducing burst energy reflects the progressively shallower region of the crust where Hall drift can build stresses effectively, as the field decays through the range $\sim 10^{12}-10^{13}\,\mathrm{G}$. Low-energy bursts from very young magnetars may represent failures of weak regions of the crust that have only recently solidified.

X-ray spectral fitting in high-energy astrophysics can be reliably accelerated using Machine Learning. In particular, Simulation-based Inference (SBI) produces accurate posterior distributions in the Gaussian and Poisson regime for low-resolution spectra, much faster than other exact approaches such as Monte Carlo Markov Chains or Nested Sampling. We now aim to highlight the capabilities of SBI for high-resolution spectra, as what will be provided by the newAthena X-ray Integral Field Unit (X-IFU). The large number of channels encourages us to use compressed representations of the spectra, taking advantage of the likelihood-free inference aspect of SBI. Two compression schemes are explored, using either simple summary statistics, such as the counts in arbitrary bins or ratios between these bins. We benchmark the efficiency of these approaches using simulated X-IFU spectra with various spectral models, including smooth comptonised spectra, relativistic reflexion models and plasma emission models. We find that using simple and meaningful summary statistics is much more efficient than working directly with the full spectrum, and can derive posterior distributions comparable to those from exact computation using nested sampling. Multi-round inference converges quickly to the good solution. Amortized single round inference requires more simulations, hence longer training time, but can be used to infer model parameters from many observations afterwards. Information from the emission lines must be accounted for using dedicated summary statistics. SBI for X-ray spectral fitting is a robust technique that delivers well calibrated posteriors. This approach shows great promises for high-resolution spectra, offering its potential for the scientific exploitation of the X-IFU. We now plan to apply it to the current era of high-resolution telescopes, and further challenge this approach with real data.

As a soft X-ray imager with unprecedentedly large field of view, EP/WXT has detected many (fast) X-ray transients, whose nature is very intriguing. Whether there is gamma-ray counterpart for the X-ray transient provides important implications for its origin. Some of them have been reported to be associated with GRB, however, a systematic study on the gamma-ray emission of these X-ray transients is lacking. In this work, we implemented a comprehensive targeted search for gamma-ray counterparts to 63 X-ray transients reported by EP/WXT during its first year of operation, using the dedicated multiple-instrument search pipeline, ETJASMIN, with GECAM-B, GECAM-C, Fermi/GBM, and \textit{Insight}-HXMT data. We find that 14 out of 63 (22\%) EP/WXT X-ray transients have gamma-ray counterparts. For other transients, ETJASMIN pipeline provided upper limit of gamma-ray emission, which is more stringent than that given by individual instrument. Moreover, we investigated the properties of the X-ray transients and their gamma-ray counterparts, and explored the relation between the x-ray transient and gamma-ray counterpart.

Recently, Zhang & Liu (2024) proposed a turbulent convection model for multiscale anisotropies of cosmic rays (CRs), with an assumption of isotropic diffusion such that the anisotropies are statistically isotropic. However, this assumption may be unrealistic for TeV CRs, whose observations have revealed the significance of the local interstellar background magnetic field. To meet the difficulty, the turbulent convection scenario needs to be extended to cover anisotropic diffusion. In this paper, we focus on the parallel diffusion with isotropic pitch-angle scattering, which may be an approximation to the transport process driven by weak hydromagnetic waves in a magnetic flux tube, where fluctuations of the wave velocities lead to the turbulent convection. The consequence is the breaking of the statistical isotropy, while the overall shape of the angular power spectrum, $ \overline{C_\ell}\propto\ell ^{-\gamma -1} $ ($ \ell\gg 1 $), remains similar to that in the isotropic diffusion model, where $ \ell $ are degrees of spherical harmonics, and $ \gamma $ is the turbulence spectral index of the convection field. It is then expected that the power-law index of the TeV CR small-scale angular power spectrum can be explained with the Kolmogorov law $ \gamma =5/3 $, irrespective of the background magnetic field to some extent.

The compact mm emission ubiquitously found in radio-quiet active galactic nuclei (RQ AGN) exhibits properties consistent with synchrotron radiation from a small region ($\leq$1 light day) and undergoing self-absorption below $\sim$100 GHz. Several scenarios have been proposed for its origin, including an X-ray corona, a scaled-down jet, or outflow-driven shocks, which can be tested via mm polarimetry. In the optically thin regime, synchrotron emission is expected to show polarization up to $\sim$70\%, but disordered magnetic fields and Faraday rotation reduce this to a few percent for jets and outflows, while an X-ray corona is likely to result in complete depolarization. To investigate this, we conducted the first ALMA Band 3 full-polarization observations of three RQ AGN - NGC 3783, MCG 5-23-16, and NGC 4945. No polarized signal was detected in any of the AGN, with an upper limit of 0.5-1.5\%, supporting the X-ray corona scenario. However, we detected a compact source with 17\% polarization in NGC 3783, 20 pc away from the AGN, co-spatial with the mm and narrow-line outflow, likely linked to a shock propagating through the outflowing material. Additionally, combining our data with archival ALMA observations, we found typical mm variability in RQ AGN by a factor of 2; however, the sparsity of the data prevented a more detailed analysis of the total flux variability.

JUICE was launched in April 2023, and it is now in its cruise phase to Jupiter, where it is scheduled to arrive in July 2031. JUICE carries a radiation monitor, namely the RADiation hard Electron Monitor (RADEM) to measure protons, electrons, and ions, detecting particles coming mainly from the anti-Sun direction. On 2024 May 13, a large solar energetic particle (SEP) event took place in association with an eruption close to the western limb of the Sun as seen from Earth. Providentially, at that time JUICE was located very close to STEREO-A, being separated by 0.13 au in radial distance, 0.3 deg in latitude, and 1.6 deg in longitude. Methods. We analysed the IP context using in-situ measurements and studied the proton anisotropies measured by near-Earth spacecraft and STEREO-A. We focused on an isotropic period during the decay phase of the SEP event to compute the proton energy spectrum. We fit the STEREO-A spectrum and compared it to that measured by SOHO and JUICE. Conclusions. The RADEM instrument aboard JUICE is a valuable tool for measuring SEP events in the heliosphere, providing an excellent opportunity to study and characterise the energetic particle environment in the solar wind between 0.65 and 5.2 au. The intercalibration factors between the fluxes measured by STEREO-A and JUICE at the effective energies of 6.9 MeV, 13.3 MeV, 21.6 MeV, and 31.2 MeV are 1.02, 1.23, 1.12, and 0.95 respectively. These intercalibration factors are valid only until 2024 July 10, when the configuration of the RADEM instrument was changed.

We present new K-band spectroscopy for the giant elliptical galaxy M87 in the Virgo cluster, taken with the LUCI spectrograph at the Large Binocular Telescope (LBT). The new data are used to study line-strengths of K-band absorption features from different chemical species, namely Fe, Mg, Ca, Na, and CO, as a function of galactocentric distance, out to 40arcsec from the center (about half of the galaxy effective radius). The radial trends of spectral indices are compared to those for the bulge of M31, observed with the same instrument. For M87, most K-band indices exhibit flat radial profiles, with the exception of NaI2.21, that decreases outwards, with a negative radial gradient. Significant offsets are found between indices for M87 and those for the bulge of M31, the latter having weaker line-strengths for almost all features, but Fe and Ca, for which we find similar trends in both systems. We find that the behavior of CO features - most prominent in giant stars - is difficult to explain, consistent with previous results for the central regions of massive galaxies. In particular, the CO indices are stronger in M87 than M31, and do not exhibit significant radial gradients in M87, despite its IMF being bottom heavier than M31 especially in its central region. Predictions of state-of-the-art stellar population models, based on results from the optical spectral range, are able to match only the Na and Ca indices of M87, while a significant mismatch is found for all other indices. This shows that state-of-the-art stellar population models should be improved significantly in order to provide reliable constraints on the stellar population content of galaxies in the NIR spectral range.

The inner ~200 pc region of the Milky Way contains a nuclear stellar disc and a nuclear star cluster that are embedded in the larger Galactic bar. These stellar systems overlap spatially, which makes it challenging to separate stars that belong to the nuclear stellar systems, to deduce their internal dynamics, and to derive the central Galactic potential. Discrete stellar kinematics probe the mass distribution of a stellar system, and chemical tracers such as stellar metallicity can further separate multiple stellar populations that can have distinct kinematic properties. We took advantage of the information provided by discrete stellar kinematics and the metallicity in the Galactic centre using discrete chemo-dynamical modelling. We fitted axisymmetric Jeans models to discrete data of 4,600 stars. We fitted the stars as either one population plus a background component or as two populations plus a background that represents the bar. We tested the robustness of the inferred gravitational potential against a varying mass of the supermassive black hole, including dark matter, or a radially varying mass-to-light ratio. We obtained robust results on the fit with a single population and a background component. We obtained a supermassive black hole mass of (4.35$\pm 0.24) \times 10^6$ M$_\odot$, and we find that a dark matter component and radial variation in the mass-to-light ratio are negligible. We derived the enclosed mass profile of the inner ~60 pc and found a lower mass than reported in the literature in the region of ~5-30 pc. In our two-population fit, we found a high-[M/H] population that contributes more than 90% to the total stellar density. The properties of the high-[M/H] population are consistent with in situ formation after gas inflow from the Galactic disc via the bar. The distinct kinematic properties of the low-[M/H] population indicate a different origin. [abridged]

Using the stellar evolution code MESA, we mimic the negative jet feedback mechanism in common envelope evolution (CEE) of low-mass main sequence stars, M2=0.1-0.3Mo, spiraling inward inside the envelopes of asymptotic giant branch (AGB) or red giant branch (RGB) stars and find that the jets reduced the envelope density, therefore the jets' power, by a factor of ~0.5/(M2/0.1Mo). We mimic the energy that the jets deposit into the envelope by depositing energy into the outer envelope, a process that inflates the envelope, therefore reducing the density in the vicinity of the main sequence star, the accretion rate, and the jets' power. In deriving this expression for the negative jet feedback coefficient, we assume that the actual mass accretion rate is a fraction 0.2-0.5 of the classical Bondi-Hoyle-Lyttleton mass accretion rate and that the jets carry a fraction 0.25-0.5 of the accretion energy onto the main sequence star. Our study is another step in establishing the major role of jets in the onset and early phase of CEE, a possible grazing envelope evolution phase, and in transient events, such as luminous red novae, which these processes can power.

Leptonic one-zone radiation models are commonly used to describe multi-wavelength data and explore the physical properties of high-energy sources, such as active galactic nuclei. However, these models often require a large number of free parameters. In the context of possible parameter degeneracy and the complex landscape of the parameter space, we study how the choice of the fitting procedure impacts the characterization of the source properties. Furthermore, we examine how the data coverage and the uncertainties associated to the data influence the model parameter degeneracy. We generated simulated spectral energy distribution datasets with different properties, which we then fit with a numerical model. The model describes the relevant radiation processes with seven free parameters. We compare different optimization algorithms and study the parameter degeneracy using t-distributed stochastic neighbor embedding. Additionally, we applied the same fitting procedures to the observational data of two sources, Mrk 501 and PKS 0735+178. We demonstrate significant degeneracies in the seven-dimensional parameter space of the one-zone leptonic models caused by the incomplete wavelength coverage of the data. Given the same goodness-of-fit function, the best-fit result depends on the choice of the minimization algorithm. Source properties extracted from the best-fit solution to realistic datasets cannot be interpreted as the only solution due to significant degeneracies of the model parameters. Adding new energy ranges (e.g. MeV) and regular source monitoring would allow to reduce gaps in the data and significantly decrease the parameter degeneracy.

OH megamasers (OHMs) are extragalactic masers found primarily in gas-rich galaxy major mergers. To date, only $\sim$120 OHMs have been cataloged since their discovery in 1982, and efforts to identify distinct characteristics of OHM host galaxies have remained inconclusive. As radio astronomy advances with next-generation telescopes and extensive 21 cm HI surveys, precursors to the Square Kilometre Array (SKA) are expected to detect the 18 cm OH masing line with significantly increased frequency, potentially expanding the known OHM population tenfold. These detections, however, risk confusion with lower-redshift HI emitters unless accompanied by independent spectroscopic redshifts. Building on methods proposed by Roberts et al. (arXiv:2102.12486) for distinguishing these interloping OHMs via near- to mid-IR photometry and emission line frequencies, we apply these techniques to data from the Arecibo Legacy Fast ALFA [Arecibo L-band Feed Array] (ALFALFA) survey and a preliminary APERture Tile In Focus (Apertif) HI emission line catalog from the Westerbork Synthesis Radio Telescope. Our study, utilizing the Apache Point Observatory 3.5m telescope to obtain optical spectroscopic redshifts of 142 candidates (107 from ALFALFA and 35 from Apertif), confirms five new OHM host galaxies and reidentifies two previously catalogued OHMs misclassified as HI emitters in ALFALFA. These findings support the predictions from Roberts et al. (arXiv:2102.12486 [astro-ph.GA]) and underscore the evolving landscape of radio astronomy in the context of next-generation telescopes.

Recent studies suggest that the density structure of turbulent molecular clouds in the Milky Way and the Andromeda galaxy, M31, aligns with expectations from hydrostatic equilibrium (HE) and virial equilibrium (VE). This study extends the study of the M31 clouds by matching their observed surface densities to a spatially-averaged solution of the Lane-Emden equation. The results affirm that the M31 molecular clouds exhibit density profiles expected from HE, further supporting the earlier conclusion that VE holds at all radii within the clouds. This study also discusses HE in the context of the turbulent interstellar medium, considering energy balance, cloud lifetimes, and the pressure of the interstellar medium. A comparison between the Galactic and extragalactic clouds indicates that both share similar dynamical states. These findings contribute to a broader understanding of the interplay between turbulence, gravitational stability, and cloud evolution in the molecular interstellar medium.

Propagating slow magneto-acoustic waves are commonly observed in different coronal structures but are most prominent in active region fan loops. Their rapid damping with damping lengths of the order of a wavelength has been investigated in the past by several authors. Although different physical mechanisms have been proposed, significant discrepancies between the theory and observations remain. Recent high-resolution observations captured simultaneously by two different instruments reveal distinct damping lengths for slow magneto-acoustic waves although their passbands are similar. These results suggest a possible contribution of instrumental characteristics on the measurement of damping lengths. Here, we analyse the behavior of slow waves using a different pair of instruments in order to check the prevalence of such results. In particular, the cotemporal observations of active region NOAA AR12712 by the High-Resolution Coronal Imager (Hi-C 2.1) and the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) are utilised. The estimated oscillation periods of slow magneto-acoustic waves identified from these data are 2.7{\,}$\pm${\,}0.2{\,}min from SDO/AIA, and 2.8{\,}$\pm${\,}1.2{\,}min from Hi-C 2.1. The corresponding propagation speeds are found to be 46.0{\,}$\pm${\,}1.7{\,}km{\,}s$^{-1}$ and 48.1{\,}$\pm${\,}0.6{\,}km{\,}s$^{-1}$, respectively. Damping lengths were calculated by two different methods, the Phase Tracking Method (PTM) and the Amplitude Tracking Method (ATM). The obtained values from PTM are 4.0{\,}$\pm${\,}2.1{\,}Mm and 4.1{\,}$\pm${\,}0.3{\,}Mm while those from ATM are 3.4{\,}$\pm${\,}1.0{\,}Mm and 3.7{\,}$\pm${\,}0.1{\,}Mm, respectively, for the AIA and Hi-C data. Our results do not indicate any notable difference in damping lengths between these instruments.

Star formation and metal enrichment in galaxies are regulated by supernova (SN) explosions and metal yields from massive stars, which are sensitive to the high-mass end of the initial mass function (IMF). Recent JWST observations have reached a consensus on an invariant relation between stellar mass, metallicity, and star formation rate up to $z\sim 8$ and its breakdown at higher redshifts. It is crucial to understand the underlying physics, especially the role played by the IMF. We explore the impact of IMF on the chemical evolution of high-redshift galaxies and the interplay between IMF and galactic outflows. The ultimate goal is to constrain the high-mass end of the IMF by the cosmic star formation history and stellar mass-metallicity-star formation rate relation (MZSFR) inferred from observations at $z\sim 4-10$. Using the semi-analytical galaxy evolution code A-SLOTH, we follow galactic baryon cycles along merger trees built from cosmological simulations. Stellar feedback is modeled with up-to-date stellar evolution tracks covering the full metallicity range ($Z \sim 10^{-11} - 0.03$) and a broad stellar mass range ($m_\star\sim2 - 600\ \rm M_\odot$) including the metal yields from stellar winds, core-collapse SNe, (pulsational) pair-instability SNe, and Type Ia SNe. Assuming that the IMF follows a Kroupa-like shape with a varying upper mass limit $m_{\max}$, we find $m_{\max} \gtrsim 200\ \rm M_\odot$ is required to reproduce the observed MZSFR. Observational data at $z\gtrsim 6$ favor a galactic outflow model where the outflow mass is proportional to the ratio of supernova energy to halo binding energy. We conclude that very massive ($\gtrsim 200\ \rm M_\odot$) stars can play important roles in the star formation and chemical enrichment histories of high-$z$ galaxies. We also discuss their implications for transient sources of both electromagnetic waves and gravitational waves.

The interstellar hydride hydroxyl (OH) is a potential tracer of CO-dark molecular gas. We present new absorption line observations of OH at 18-cm wavelength towards four continuum sources. We compare these to the [CII] line at 1.9 THz obtained with SOFIA, observations of the neutral atomic hydrogen 21 cm line with the VLA, and CO lines obtained with APEX. We trace OH over a large range of molecular hydrogen column densities, and derive OH abundances with respect to molecular and total hydrogen column densities. Increased sensitivity and spectral resolution allowed us to detect weak and narrow features. We identify only one OH absorption component out of 23 without CO counterpart, yet several with intermediate molecular gas fractions. A potential association of [CII] 158 mu m emission with an OH absorption component is seen toward one sightline. Our results confirm that OH absorption traces molecular gas across diffuse and dense environments of the interstellar medium. At the sensitivity limits of the present observations our detection of only one CO-dark molecular gas feature appears in agreement with previous studies. We conclude that if OH absorption was to be used as a CO-dark molecular gas tracer, deeper observations or stronger background targets are necessary to unveil its full potential as a CO-dark molecular gas tracer, and yet it will never be an exclusive tracer of CO-dark molecular gas. For OH hyperfine-splitting transitions in the vicinity of photodissociation regions in W43-South, we detect a spectral and spatial offset between the peak of the inversion of the OH 1612 MHz line and the absorption of the OH 1720 MHz line on the one hand, and the absorption of the OH main lines on the other hand, which provides additional constraints on the interpretation of the OH 18 cm line signatures typical of HII regions.

Venus has regained attention on the international stage with the approval of three new missions by ESA and NASA. As the twin sister of Earth, Venus exhibits a distinct atmosphere, which casts a veil of mystery over the planetary evolution and is of great scientific significance. One of the most important components of Venus-the cloud-is believed to have significantly regulated its climate evolution and affect the environmental habitability. However, due to sparse in-situ measurements and the limitation of remote sensing, properties of these clouds remain largely unknown. Based on research conducted in past decades, this article reviews the observational structure of cloud properties, the progress of microphysical and simplified cloud model developments, and perspectives of future directions of this research field. Several possible solutions to the challenges associated with the coupling effect, ultraviolet absorption, and habitability are proposed and discussed in details, providing insights for future Venus' explorations.

Primordial black holes (PBHs) constitute a compelling dark matter candidate whose gravitational effects could significantly influence early cosmic structure formation. We investigate the impact of PBHs on Population III star formation through detailed $N$-body and hydrodynamic simulations, extending beyond previous semi-analytical approaches. Our results reveal a mass-dependent dichotomy in PBH effects: massive PBHs ($M_{\rm PBH} \gtrsim 10^2 M_\odot$) with sufficient abundance can accelerate structure formation and shift Pop III formation to higher redshifts, potentially conflicting with observational constraints from high-redshift galaxy surveys. Conversely, lower-mass PBHs can induce tidal disruption of gas-rich minihalos, suppressing star formation and delaying the cosmic dawn depending on their abundance. We quantify these competing effects to derive new constraints on the PBH mass function and their contribution to the total dark matter density, with implications for forthcoming observations with the James Webb Space Telescope and 21-cm cosmology experiments.

The detection of gravitational waves (GWs) by LIGO-Virgo and pulsar timing arrays (PTAs) has opened a new window into early universe cosmology. Yet, the origin of large-scale magnetic fields and the dynamics of the reheating epoch remain poorly understood. In this work, we study the generation of secondary GWs (SGWs) sourced by primordial magnetic fields produced via a Sawtooth-type coupling during reheating with a general background evolution. We show that the reheating equation of state significantly influences the spectral shape and amplitude of the magnetic fields. While a scale-invariant spectrum is typically needed to match observational bounds, this coupling naturally produces a strongly blue-tilted spectrum that remains consistent with current constraints. Crucially, the magnetic field continues to grow during reheating, leading to a GW signal with a broken power-law spectrum and a distinctive blue tilt on super-horizon scales. This SGW signal can fall within the sensitivity of upcoming detectors such as LISA, DECIGO, and BBO. The unique spectral features make this scenario distinguishable from other sources, offering a viable mechanism for cosmic magnetogenesis and a novel probe of the reheating era through GW observations.

We present the validation of TOI-2407 b, a warm Neptune-sized planet with a radius of 4.26 $\pm$ 0.26 R$_\oplus$, orbiting an early M-type star with a period of 2.7 days and an equilibrium temperature of 705 $\pm$ 12 K. The planet was identified by TESS photometry and validated in this work through multi-wavelength ground-based follow-up observations. We include an observation with the novel CMOS-based infrared instrument SPIRIT at the SPECULOOS Southern Observatory. The high-precision transit data enabled by CMOS detectors underscore their potential for improving the detection and characterisation of exoplanets orbiting M-dwarfs, particularly in the infrared, where these stars emit most of their radiation. TOI-2407 b lies within the boundaries of the period-radius Neptune desert, an apparent scarcity of Neptune-sized planets at short orbits. Further characterisation of TOI-2407 b, such as radial velocity measurements, will refine its position within planetary demographic trends. This system also provides a comparison case for the well-studied Neptune-sized planet Gliese 436 b, of similar radius, period and stellar type. Comparison studies could aid the understanding of the formation and evolution of Neptune-like planets around M-dwarfs.

New generation galaxy surveys targeting constraints on local primordial non-Gaussianity (PNG) demand $N$-body simulations that accurately reproduce its effects. In this work, we explore various prescriptions for the initial conditions of simulations with PNG, aiming to optimise accuracy and minimise numerical errors, particularly due to aliasing. We have used $186$ runs that vary the starting redshift, LPT order, and non-Gaussianities ($f^{\rm local}_{\rm NL}$ and $g^{\rm local}_{\rm NL}$). Starting with $3$LPT at a redshift as low as $z_{\rm ini}\simeq 11.5$ reproduces to $<1 \%$ the power spectrum, bispectrum and halo mass function of a high-resolution reference simulation. The aliasing induced by the PNG terms in the power spectrum produces a $ \leq 3 \%$ excess small-scale power at the initial conditions, dropping below $0.1\%$ by $z=0$. State-of-the-art initial condition generators show a sub-percent agreement. We show that initial conditions for simulations with PNG should be established at a lower redshift using higher-order LPT schemes. We also show that removing the PNG aliasing signal is unnecessary for current simulations. The methodology proposed here can accelerate the generation of simulations with PNG while enhancing their accuracy.

The properties of satellite halos provide a promising probe for dark matter (DM) physics. Observations motivate current efforts to explain surprisingly compact DM halos. If DM is not collisionless but has strong self-interactions, halos can undergo gravothermal collapse, leading to higher densities in the central region of the halo. However, it is challenging to model this collapse phase from first principles. To improve on this, we seek to better understand numerical challenges and convergence properties of self-interacting dark matter (SIDM) N-body simulations in the collapse phase. Especially we aim for a better understanding of the evolution of satellite halos. To do so, we run SIDM N-body simulations of a low mass halo in isolation and within an external gravitational potential. The simulation setup is motivated by the perturber of the stellar stream GD-1. We find that the halo evolution is very sensitive to energy conservation errors, and a too large SIDM kernel size can artificially speed up the collapse. Moreover, we demonstrate that the King model can describe the density profile at small radii for the late stages that we have simulated. Furthermore, for our highest-resolved simulation (N = 5x10^7) we make the data public. It can serve as a benchmark. Overall, we find that the current numerical methods do not suffer from convergence problems in the late collapse phase and provide guidance on how to choose numerical parameters, e.g. that the energy conservation error is better kept well below 1%. This allows to run simulations of halos becoming concentrated enough to explain observations of GD-1 like stellar streams or strong gravitational lensing systems.

Surprisingly compact substructures in galaxies and galaxy clusters, but also field halos, have been observed by gravitational lensing. They could be difficult to explain with collisionless dark matter (DM). To explain those objects, recent studies focused on the gravothermal collapse that halos consisting of self-interacting dark matter (SIDM) can undergo. However, simple models of elastic scattering could face problems explaining those compact objects during very later stages of the collapse and the post-collapse phase, where a black hole may have formed from DM. We aim to explain compact halos while avoiding the gravothermal catastrophe which typical SIDM models are subject to. Therefore, we investigate the evolution of a DM halo for an SIDM model consisting of two species with unequal masses, featuring only interactions between the different species but not within themselves. Employing $N$-body simulations, we study the effect of unequal-mass SIDM models on the evolution of an isolated DM halo. In particular, the late stages of its evolution with high central densities are simulated. We find that our two-species SIDM models can produce density cores with their size depending on the mass ratio of the two species. Moreover, the mass segregation caused by the unequal particle masses leads to a finite final density state or at least a slowly growing density, which depends on the mass ratio and the mass fraction of the two DM species. SIDM models consisting of two DM species can simultaneously explain DM halos with density cores, as well as systems that are denser in their centre than expected from collisionless DM, while avoiding the gravothermal catastrophe. They are a compelling alternative to single-species models, offering a rich phenomenology.

The Landau-Migdal parameter $G'_0$ characterizes the main part of the spin-isospin nucleon-nucleon interaction. Consequently, the $G'_0$ is closely related to the Gamow-Teller resonance (GTR), the beta and double-beta decay rates of finite nuclei, the spin response of hot and dense nucleonic matter that determines the neutrino-nucleon reaction rates in core-collapse supernovae (CCSNe) and binary neutron star (BNS) mergers, and finally the critical density for pion condensation in neutron stars. Historically, the $G'_0$ was obtained by fitting the peak location of experimental GTR spectra by using phenomenological pion exchange models, without strict uncertainty quantification. In this letter, for the first time, we report the Bayesian inference of $G'_0$ by using a self-consistent Skyrme Quasiparticle Random Phase Approximation (QRPA) model and joint constraints from experimental GTR measurements on $^{208}\mathrm{Pb}$, $^{132}\mathrm{Sn}$, $^{90}\mathrm{Zr}$. Our extracted $G_0'$ is $0.48\pm0.034$, which is close to the prediction of a few existing Skyrme models that consider the spin-isospin observables but is smaller than the extracted ones from pion-exchange models. We hint to possible reasons for this deviation, like the value of the nucleon effective mass $\frac{m^*}{m}$. Finally, we demonstrate the influence of $G'_0$ on neutrino opacities in CCSNe and BNS mergers. The new Skyrme parameterizations from our Bayesian study may also be used to study other spin-isospin-dependent phenomena.

Accurate modeling of the neutron star crust is essential for interpreting multimessenger observations and constraining the nuclear equation of state (EoS). However, standard phenomenological EoS models often rely on heuristic extrapolations in the low-density regime, which are inconsistent with microscopic predictions. In this work, we refine a unified meta-modeling framework for the EoS by incorporating low-density corrections based on energy density functionals constrained by ab initio neutron-matter calculations. Using Bayesian inference to combine information from astrophysical observations, nuclear theory, and experiments, we assess the impact of these corrections on key crustal properties, including the crust-core transition density and pressure, crustal composition, and moment of inertia. The improved model reduces uncertainties in the inner crust and emphasizes the importance of low-density physics in EoS modeling, highlighting the value of integrating both theoretical and observational constraints across densities to robustly describe the EoS. Moreover, the adopted approach can be readily applied to any existing EoS model to provide a solid framework for interpreting upcoming high-precision multimessenger data.

The development of science has been transforming man's view towards nature for centuries. Observing structures and patterns in an effective approach to discover regularities from data is a key step toward theory-building. With increasingly complex data being obtained, revealing regularities systematically has become a challenge. Correlation is a most commonly-used and effective approach to describe regularities in data, yet for complex patterns, spatial inhomogeneity and complexity can often undermine the correlations. We present an algorithm to derive maps representing the type and degree of correlations, by taking the two-fold symmetry of the correlation vector into full account using the Stokes parameter. The method allows for a spatially resolved view of the nature and strength of correlations between physical quantities. In the correlation view, a region can often be separated into different subregions with different types of correlations. Subregions correspond to physical regimes for physical systems, or climate zones for climate maps. The simplicity of the method makes it widely applicable to a variety of data, where the correlation-based approach makes the map particularly useful in revealing regularities in physical systems and alike. As a new and efficient approach to represent data, the method should facilitate the development of new computational approaches to regularity discovery.

Physics has been transforming our view of nature for centuries. While combining physical knowledge with computational approaches has enabled detailed modeling of physical systems' evolution, understanding the emergence of patterns and structures remains limited. Correlations between quantities are the most reliable approach to describe relationships between different variables. However, for complex patterns, directly searching for correlations is often impractical, as complexity and spatial inhomogeneity can obscure correlations. We discovered that the key is to search for correlations in local regions and developed a new method, adjacent correlation analysis, to extract such correlations and represent them in phase space. When multiple observations are available, a useful way to study a system is to analyze distributions in phase space using the Probability Density Function (PDF). Adjacent correlation analysis evaluates vectors representing local correlations, which can be overlaid on the PDF plot to form the adjacent correlation plot. These correlation vectors often exhibit remarkably regular patterns and may lead to the discovery of new laws. The vectors we derive are equivalent to the vector field in dynamical systems on the attracting manifold. By efficiently representing spatial patterns as correlation vectors in phase space, our approach opens avenues for classification, prediction, parameter fitting, and forecasting.

We study the effect of light-bending on the signal of a pulsar in binaries with rotating black hole companions, focusing on stellar mass black holes. We show that the impacts of various parameters on the bending delays visually match with those observed for a non-rotating black holes, because the magnitude of the spin as well as the orientation of the spin axis of the black hole introduce changes in the nanosecond order and other parameters do so in the microsecond order. Consequently, the distortion of the beam and the resulting changes in the pulse shape are minimally influenced by spin-related parameters of the black hole. We also investigate the impact of various parameters on the difference of the delays with and without the spin of the black hole and notice nanosecond scale discontinuities at orbital phases where the path of the light ray changes its direction with respect to the direction of the spin of the black hole. Moreover, as in the Schwarzschild case, the bending delays become irregular (on the microsecond scale) near the superior conjunction. We also explore the effect of bending on the pulse profiles and bending delays if the companion of the pulsar is a rotating super-massive black hole. We find significant enhancement and change in the shape of the profiles at and near the superior conjunction in comparison to stellar mass black holes. Moreover, bending delays are about three orders of magnitude higher than those in case of the stellar mass black holes.

Accurate parameter estimation(PE) of gravitational waves(GW) is essential for GW data analysis. In extreme mass-ratio inspiral binary(EMRI) systems, orbital eccentricity is a critical parameter for PE. However, current software for for PE of GW often neglects the direct estimation of orbital eccentricity. To fill this gap, we have developed the MatBYIB, a MATLAB-based software package for PE of GW with arbitrary eccentricity. The MatBYIB employs the Analytical Kludge (AK) waveform as a computationally efficient signal generator and computes parameter uncertainties via the Fisher Information Matrix (FIM) and the Markov Chain Monte Carlo (MCMC). For Bayesian inference, we implement the Metropolis-Hastings (M-H) algorithm to derive posterior distributions. To guarantee convergence, the Gelman-Rubin convergence criterion (the Potential Scale Reduction Factor R) is used to determine sampling adequacy, with MatBYIB dynamically increasing the sample size until R < 1.05 for all parameters. Our results demonstrate strong agreement between FIM- based predictions and full MCMC sampling. This program is user-friendly and allows for estimation of gravitational wave parameters with arbitrary eccentricity on standard personal computers.

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

Parametric representations of the high-density nuclear equation of state are used in constructing models for interpreting the astrophysical observations of neutron stars. This study explores how accurately equations of state with strong first-order phase transitions can be represented using spectral or piecewise analytic methods that assume no {\it{a priori}} knowledge of the location or the strength of the phase transition. The model equations of state used in this study have phase transitions strong enough to induce a gravitational instability that terminates the sequence of stable neutron stars. These equations of state also admit a second sequence of stable stars with core matter that has undergone this strong first-order phase transition (possibly driven by quark deconfinement). These results indicate that spectral representations generally achieve somewhat higher accuracy than piecewise analytic representations having the same number of parameters. Both types of representation show power-law convergence at approximately the same rate.

In the last decade or so, metasurface optical components have received considerable scientific and industrial interest for a variety of applications. The miniaturization afforded by metasurfaces could benefit astronomy in particular (which is an often-cited potential application area for metasurfaces). However, few developed examples in which metasurface components offer a unique benefit to astronomical instrumentation - substantiated by the production of scientific data - have been shown. Here, we present the Solar Imaging Metasurface Polarimeter (SIMPol), a first-of-its-kind telescope for snapshot imaging polarimetry of the sun around a Sr I line at 460.7 nm enabled by a metasurface polarization-analyzing grating. This high-performance grating exhibits an overall efficiency of nearly 70% and high polarization contrast (diattenuation) across its four observation channels. We demonstrate SIMPol's integration into a major observatory telescope facility with two different imaging modes. In both cases, Zeeman polarization signatures were clearly observed in two adjacent spectral lines of Fe I and Sr I around 460.7 nm. This work demonstrates an early success of metasurface polarization optics in a real application in astronomical instrumentation (here, polarimetric observations of the solar atmosphere), and heralds the application of metasurfaces and emergent nanophotonic technologies in astronomy more broadly.