Papers for Monday, Jun 23 2025

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

The profound comprehension of the evolution and phenomenology of an Active Galactic Nucleus requires an accurate exploration of the dynamics of the magnetized gaseous disk surrounding the massive black hole in the centre. Many numerical simulations have studied this environment using elaborate grid-based codes, but in recent years, new mesh-less schemes have exhibited excellent conservation properties and good accuracy at a more moderate computational cost. Still, none implement general relativistic magnetic fields, a fundamental ingredient to model an accretion disk around a massive black hole. We present here a general relativistic magnetohydrodynamics (GRMHD) scheme working within the mesh-less framework of the code \texttt{GIZMO}. We implement the hyperbolic divergence cleaning procedure, consistently extended to general relativistic effects, to keep the magnetic field divergence under safe levels. We benchmark the scheme against various relativistic magnetohydrodynamics stress tests, considering different dimensionalities and both a Minkowski or a Schwarzchild/Kerr background. To date, this is the first GRMHD scheme working in a mesh-free environment.

The mass assembly and chemical enrichment of the first galaxies provide key insights into their star-formation histories and the earliest stellar populations at cosmic dawn. Here we compile and utilize new, high-quality spectroscopic JWST/NIRSpec Prism observations from the JWST archive. We extend the wavelength coverage beyond the standard pipeline cutoff up to 5.5$\mu$m, enabling a detailed examination of the rest-frame optical emission-line properties for galaxies at $z\approx 10$. The improved calibration allows us to detect H$\beta$ and the [OIII]$\lambda\lambda 4959,5007$ doublet and resolve the auroral [OIII]$\lambda 4363$ line for the 11 galaxies in our sample ($z=9.3-10.0$) to obtain direct $T_e$-based metallicity measurements. We find that all galaxies show high ionisation fields and electron temperatures, with derived metallicities in the range $12+\log {\rm (O/H)} = 7.1 - 8.3$, consistent with previous strong-line diagnostics at high-z. We derive an empirical relation for $M_{\rm UV}$ and 12+log(O/H) at $z\approx 10$, useful for future higher-z studies, and show that the sample galaxies are `typical' star-forming galaxies though with relatively high specific star-formation rates and with evidence for bursty star formation. Combining the rest-frame optical line analysis and detailed UV to optical SED modelling, we determine the mass-metallicity relation and the fundamental-metallicity relation of the sample, pushing the redshift frontier of these measurements to $z=10$. These results, together with literature measurements, point to a gradually decreasing MZR at higher redshifts, with a break in the FMR at $z\approx 3$, decreasing to metallicities $\approx 3\times$ lower at $z=10$ than observed during the majority of cosmic time at $z=0-3$, likely caused by massive pristine gas inflows diluting the observed metal abundances during early galaxy assembly at cosmic dawn.

We demonstrate that there are no scale-invariant one-loop corrections to the superhorizon tensor perturbations from small-scale (potentially enhanced) scalar perturbations, irrespective of the details of inflationary background time evolution. For this purpose we derive a soft tensor effective field theory at leading order in the gradient expansion by integrating out small-scale scalar fluctuations in a general time-dependent background over the Schwinger-Keldysh contour, i.e., we perform loop calculations in the soft limit of external momentum. The absence of scale-invariant corrections originates from the diffeomorphism invariance of general relativity and is therefore unavoidable.

The existence of $\sim10^{7-8}{\rm M}_\odot$ supermassive black holes at $z\gtrsim 8$ challenges conventional growth channels. One attractive possibility is that light seeds ($M_{\rm BH}\lesssim10^{3}{\rm M}_\odot$) undergo short, super-Eddington episodes when they cross, and are captured by, dense massive gas clumps. We revisit this ``BH-clump-capture'' model using analytic arguments supported by toy-model simulations that follow Bondi-scale inflow, radiative feedback, gas dynamical friction and the recently discovered forward acceleration effect caused by the ionised bubble. For substantial growth the black hole must remain trapped for many dynamical times, which imposes three simultaneous constraints. The clump must be heavier than the black hole (mass doubling condition); its cooling time must exceed the super-Eddington growth time (lifetime condition); and dynamical friction must dominate shell acceleration (BH-trapping condition). These requirements confine viable clumps to a narrow density-temperature region, $n_{\rm H}\simeq10^{7-8}{\rm cm}^{-3}$ and $T\simeq(2-6)\times10^{3}{\rm K}$, for a $10^{3}{\rm M}_\odot$ seed. Even inside this sweet spot a $10^{3}{\rm M}_\odot$ seed grow up at most $4\times10^{3}{\rm M}_\odot$; the maximum growth ratio $M_{\rm BH, fin}/M_{\rm BH}$ falls approximately as $M_{\rm BH}^{-0.4}$ and is negligible once $M_{\rm BH}\gtrsim10^{4}{\rm M}_\odot$. The forward-acceleration effect is essential, expelling the black hole whenever photon trapping fails. We conclude that BH-clump-capture model, and potentially broad super-Eddington models, cannot produce the $>10^{4}{\rm M}_\odot$ seeds required for subsequent Eddington-limited growth, suggesting that alternative pathways, such as heavy seed formation, remain necessary.

Horizon-scale imaging of supermassive black holes has opened a new window onto the studies of strong-field gravity and plasma physics in low-luminosity accretion flows. As future efforts aim to image fainter and smaller angular-size targets, primarily through space-based very long baseline interferometry (VLBI), it is important to identify optimal sources and observing strategies for such studies. In this work, we assess the prospects for imaging black hole shadows in a broad population of nearby supermassive black holes by modeling their accretion flows using a covariant semi-analytic model for the flow and general relativistic ray tracing. We explore the influence of black hole and accretion flow parameters on spectra, image morphology, and the critical frequency at which the flows become optically thin. We identify three general classes of sources: those that become transparent at traditional imaging frequencies; those requiring higher frequencies; and those unlikely to be transparent down to the black hole shadow in the submillimeter band. Our results will inform target selection and wavelength optimization for future VLBI arrays, where both resolution and transparency are essential for resolving black hole shadows.

We explore how a realistic surface brightness detection limit of $\mu_V \approx 32.5$ mag arcsec$^{-2}$ for stars at the edges of ultra-faint galaxies affects our ability to infer their underlying properties. We use a sample of 19 galaxies with stellar masses $\approx 400 - 40,000~{\rm M}_\odot$ simulated with FIRE-2 physics and baryonic mass resolution of $30~M_{\odot}$. The surface brightness cut leads to smaller sizes, lower stellar masses, and lower stellar velocity dispersions than the values inferred without the cut. However, by imposing this realistic limit, our inferred galaxy properties lie closer to observed populations in the mass-size plane, better match observed velocity dispersions as a function of stellar mass, and better reproduce derived circular velocities as a function of half-light radius. For the most massive galaxies in our sample, the surface brightness cut leads to higher mean $\rm [Fe/H]$ values, but the increase is not enough to match the observed MZR. Finally, we demonstrate that the common Wolf et al. (2010) mass estimator is less accurate when the surface brightness cut is applied. For our lowest-mass galaxies, in particular, excluding the low-surface brightness outskirts causes us to overestimate their central dark-matter densities and virial masses. This suggests that attempts to use mass estimates of ultra-faint galaxies to constrain dark-matter physics or to place constraints on the low-mass threshold of galaxy formation must take into account surface brightness limits or risk significant biases.

Synthesizer is a fast, flexible, modular, and extensible Python package that empowers astronomers to turn theoretical galaxy models into realistic synthetic observations - including spectra, photometry, images, and spectral cubes - with a focus on interchangeable modelling assumptions. By offloading computationally intensive tasks to threaded C++ extensions, Synthesizer delivers both simplicity and speed, enabling rapid forward-modelling workflows without requiring users to manage low-level data processing and computational details.

Observations of Europa's leading hemisphere reveal elevated H2O2 in the warmer, low latitude chaos terrains compared to the colder, polar regions. This distribution disagrees with prior laboratory radiolysis studies of pure water ice, which show higher H2O2 yields at colder temperatures. The regions with higher peroxide abundance, Tara and Powys Regiones, also present increased amounts of CO2, possibly sourced from Europa's interior. To investigate whether CO2 influences radiolysis of water ice to boost H2O2 production, we irradiated water ice doped with varying amounts of CO2 with 10 keV electrons at 70 and 100 K. Our results indicate that CO2, even in trace amounts (< 3%), significantly enhances H2O2 yields at temperatures relevant to Europa. We discuss the mechanisms by which CO2 boosts peroxide synthesis and quantify H2O2 creation and destruction cross sections and G-values across different CO2 concentrations. These findings provide a plausible explanation for the perplexing H2O2 distribution on Europa and has implications for understanding peroxide on other icy bodies such as Ganymede and Charon, where it has been detected alongside CO2.

IXPE has enabled the X-ray polarizations of many blazars to be measured. We analyze X-ray spectropolarimetric data for all high synchrotron peaked blazars observations published by the IXPE team. We find that there are no statistically significant correlations between the X-ray polarimetric properties of the blazars and their X-ray spectral properties. This implies that there is no connection between the energy distribution of the X-ray emitting electrons and the uniformity of the magnetic field in the X-ray emitting regions of these blazars. We also find that the X-ray polarizations vary significantly between blazars and observations of the same blazar, implying the level of uniformity of the magnetic field in the X-ray emitting region is highly variable. This result will inform future theoretical work and potentially help narrow down the acceleration process of the synchrotron electrons.

We present a study of the neutral atomic hydrogen (HI) content of spatially resolved, low-redshift galaxies in the SIMBA cosmological simulations. We create synthetic HI data cubes designed to match observations from the Apertif Medium-Deep HI imaging survey, and follow an observational approach to derive the HI size-mass relation. The HI size-mass relation for SIMBA is in broad agreement with the observed relation to within 0.1 dex, but SIMBA galaxies are slightly smaller than expected at fixed HI mass. We quantify the HI spectral ($A_{\mathrm{flux}}$) and morphological ($A_{\mathrm{mod}}$) asymmetries of the galaxies and motivate standardizing the relative spatial resolution when comparing values in a sample that spans several orders of magnitude in HI mass. Galaxies are classified into three categories (isolated, interacted, or merged) based on their dynamical histories over the preceding ~2 Gyr to contextualize disturbances in their HI reservoirs. We determine that the interacted and merged categories have higher mean asymmetries than the isolated category, with a larger separation between the categories' $A_{\mathrm{mod}}$ distributions than between their $A_{\mathrm{flux}}$ distributions. For the interacted and merged categories, we find an inverse correlation between baryonic mass and $A_{\mathrm{mod}}$ that is not observed between baryonic mass and $A_{\mathrm{flux}}$. These results, coupled with the weak correlation found between $A_{\mathrm{flux}}$ and $A_{\mathrm{mod}}$, highlight the limitations of only using $A_{\mathrm{flux}}$ to infer the HI distributions of spatially unresolved HI detections.

Among the various classes of fast optical transients (FOTs), kilonovae (KNe), which can emerge as a result of neutron star mergers, are extremely challenging to observe because of not only the rapid timescale on which they fade (on the order of days), but also due to the relative scarcity of their occurrence. This scarcity is compounded by the large number of other FOTs that may initially resemble the characteristic rise of a KNe. While these objects can be ruled out as candidate KNe by taking spectroscopy, a method of confidently ruling out transients based on photometric analysis alone would be incredibly valuable. We describe the compilation of various ``imposter" transients, including a plurality of IIb SNe, and investigate a number of comparative metrics by which one might be able to remove transients from consideration without the use of spectroscopy. We provide a list of these objects and their classifications as well as a glossary of the transient types included in the sample.

ComPair, a prototype of the All-sky Medium Energy Gamma-ray Observatory (AMEGO), completed a short-duration high-altitude balloon campaign on August 27, 2023 from Fort Sumner, New Mexico, USA. The goal of the balloon flight was the demonstration of ComPair as both a Compton and Pair telescope in flight, rejection of the charged particle background, and measurement of the background $\gamma$-ray spectrum. This analysis compares measurements from the balloon flight with Monte Carlo simulations to benchmark the instrument. The comparison finds good agreement between the measurements and simulations and supports the conclusion that ComPair accomplished its goals for the balloon campaign. Additionally, two charged particle background rejection schemes are discussed: a soft ACD veto that records a higher charged particle event rate but with less risk of event loss, and a hard ACD veto that limits the charged particle event rate on board. There was little difference in the measured spectra from the soft and hard ACD veto schemes, indicating that the hard ACD veto could be used for future flights. The successes of ComPair's engineering flight will inform the development of the next generation of ComPair with upgraded detector technology and larger active area.

Non-thermal velocities, derived from spectral line broadening, can provide crucial insights into plasma dynamics before and during solar flares. To systematically study the pre-flare phase, we constructed a Hinode/Extreme-ultraviolet Imaging Spectrometer (EIS) flare catalog of 1,449 flares from 2011-2024. This enabled a large-scale analysis of flare loop footpoint non-thermal velocity evolution across different flare magnitudes (C, M, X-classes). Analyzing Fe VIII - Fe XXIV emission lines formed at log(T/K) ~5.7-7.3 with piecewise linear fits in the pre-flare period, we find that non-thermal velocities consistently increase 4-25 minutes before GOES soft X-ray start in C and M-class flares. Onset timing patterns vary with flare magnitude: smaller flares show temperature-dependent progression, while larger flares exhibit more compressed, near-simultaneous onsets across temperatures. While our limited X-class sample (N=18) also show onset before GOES, larger statistics are needed to confirm its behavior. M-class flares show a systematic precursor non-thermal velocity onset ~30-60 minutes before GOES peak. In a subset of M-class flares (2011-2018), eruptive events show even earlier non-thermal velocity increases ($>$1 hour before GOES peak) than confined events, suggesting a link between extended pre-flare non-thermal broadening and successful eruptions. This large scale study establishes pre-flare non-thermal velocity increase at footpoints as a common precursor observable before any X-ray signature.

We report on large-scale radio observations of the Chamaeleon star-forming region obtained with the Australia Telescope Compact Array (ATCA) that led to the definite detection of five young stars and the tentative detection of five more. As in other regions surveyed in the radio domain, the majority of detected sources are fairly evolved low-mass T Tauri stars, but we also detect one protostellar object (Ced 110 IRS4) and one Herbig Ae/Be star. With the exception of the protostellar source, the radio emission mechanism is likely of non-thermal origin. The three brightest radio stars identified with ATCA were subsequently observed with the Australian Long Baseline Array (LBA) and one, J11061540-7721567 (Ced 110 IRS2), was detected at three epochs. This confirms the non-thermal nature of the radio emission in that specific case, and enabled accurate radio position measurements. Comparison with predictions from Gaia DR3 strongly suggests that this star is a binary system with an orbital period of order 40 years; additional LBA observations in the next decades would enable accurate determinations of the individual stellar masses in that system.

Joy's law describes the tilt of bipolar active regions on the Sun away from an east-west orientation, where the flux of the polarity concentrated at the prograde side tends to be closer to the equator than the polarity on the retrograde side. Joy's law is attributed to the Coriolis force because of the observed increase in tilt angle at higher latitudes. This tilt plays a crucial role in some solar dynamo models. Our goal is to model the effects of the Coriolis force on a flux tube as it rises through the near-surface convection zone. We use a three-dimensional Cartesian magnetohydrodynamic simulation of an untwisted flux tube ascending from a depth of 11 Mm. We model the Coriolis effect using the f-plane approximation, that only considers and acts on horizontal flows. On the Sun, Joy's law is weak and is only evident as an average over many active regions. To achieve a measurable effect in a single simulation, we consider a rotation rate 110 times faster than the Sun. The simulation shows that the flux tube emerges at the surface with a tilt angle consistent with Joy's law when scaled to the Sun's slower rotation, and the tilt angle does not substantially change after emergence. This shows that the Coriolis force acting on flows horizontal to the surface within the near-surface convection zone is consistent with Joy's law.

We investigate the chemical abundance distributions of the Fornax, Sculptor, Ursa Minor, and Draco dwarf galaxies using Subaru/HSC photometric data. The HSC dataset, which includes broadband g and i filters and the narrowband NB515 filter, offers sensitivity to iron and magnesium abundances as well as surface gravity, enabling the identification of giant stars and foreground dwarfs. For analysis, we selected a total of 6713 giant candidates using a Random Forest regressor trained on medium-resolution (R ~ 6000) Keck/DEIMOS spectroscopic data. Our analysis reveals the extent of radial metallicity gradients in the galaxies. Such trends, not detectable in earlier studies, are now captured owing to the substantially enlarged sample size and areal coverage provided by the HSC data. These results are also consistent with chemical abundance patterns previously observed in the central regions through spectroscopic studies. Furthermore, we infer that Fornax underwent extended star formation, whereas Sculptor formed both metal-poor and metal-rich stars over a shorter time. Ursa Minor and Draco appear to have experienced brief, intense star formation episodes leading to nearly extinguished star formation. This study underscores the critical role of the expanded HSC dataset in revealing chemical gradients that were previously inaccessible. Future work incorporating additional spectra of metal-poor stars and age-sensitive isochrone modeling will enable more accurate maps of chemical abundance distributions.

The observation of neutron star mergers with gravitational waves (GWs) has provided a new method to constrain the dense-matter equation of state (EOS) and to better understand its nuclear physics. However, inferring nuclear microphysics from GW observations necessitates the sampling of EOS model parameters that serve as input for each EOS used during the GW data analysis. The sampling of the EOS parameters requires solving the Tolman-Oppenheimer-Volkoff (TOV) equations a large number of times -- a process that slows down each likelihood evaluation in the analysis on the order of a few seconds. Here, we employ emulators for the TOV equations built using multilayer perceptron neural networks to enable direct inference of nuclear EOS parameters from GW strain data. Our emulators allow us to rapidly solve the TOV equations, taking in EOS parameters and outputting the associated tidal deformability of a neutron star in only a few tens of milliseconds. We implement these emulators in \texttt{PyCBC} to directly infer the EOS parameters using the event GW170817, providing posteriors on these parameters informed solely by GWs. We benchmark these runs against analyses performed using the full TOV solver and find that the emulators achieve speed ups of nearly \emph{two orders of magnitude}, with negligible differences in the recovered posteriors. Additionally, we constrain the slope and curvature of the symmetry energy at the 90\% upper credible interval to be $L_{\rm sym}\lesssim106$ MeV and $K_{\rm sym}\lesssim26$ MeV.

We apply principal component analysis (PCA) to the integrated intensity maps of 13 molecular lines of the nearby type-2 Seyfert galaxy NGC 1068 obtained by Atacama Large Millimeter/sub-millimeter Array (ALMA) to objectively visualize the features of its center, (1) within a radius of about 2 kpc ($\sim$ 27".5; hereafter the "overall region") and (2) the ring shaped starburst region between 750 pc ($\sim$ 10") and 2 kpc ($\sim$ 27".5) of the galaxy (hereafter the "SB ring region"). PCA is a powerful unsupervised machine learning technique that extracts key information through dimensionality reduction. The PCA results for the overall region have a possibility to reconstruct a map representing the approximate H$_2$ column density and difference of volume density and/or chemical composition between the circumnuclear disk (CND) and the starburst ring (SB ring). Additionally, the PCA results for the SB ring region have a possibility to reconstruct a map representing the approximate H$_2$ column density and distinction between starburst dominated region and shock dominated region. Furthermore, the PCA results for the SB ring region indicate a possible interaction between the Active Galactic Nucleus (AGN) outflow and gas in the SB ring. Although further investigation is required, we suggest that the AGN outflow interacts with gas in the SB ring, as this feature is consistent with the direction of the AGN outflow and is contributed by CN, C$_2$H and HCN, which are known to be enhanced by the AGN outflow. These results demonstrate that PCA can effectively extract features even for galaxies with complex structures, such as AGN + SB ring. This study also implies that PCA has the potential to uncover previously unrecognized phenomena by visualizing latent structures in multi-line data.

Recent cosmological parameter analyses combining DESI DR2 Baryon Acoustic Oscillation (BAO) data with external probes, such as Pantheon+ Supernovae (SNe) observations, have reported deviations of the dark energy equation-of-state parameters ($\oo, \oa$) from the standard $\Lambda$CDM model predictions ($\oo=-1, \oa=0$). A notable aspect of these results is the role of $\Omo$ prior information from SNe, which is known to exhibit tension with BAO-only constraints. In this study, we rigorously investigate this effect through a statistical analysis using 100 mock DESI DR2 BAO data realizations. We demonstrate that the strong degeneracy between $\oo$, $\oa$, and $\Omo$ causes significant biases in the estimated dark energy parameters when the $\Omo$ prior mean deviates from its true underlying value. Specifically, applying an $\Omo$ prior mean of 0.33 (consistent with some SNe-only constraints) to mock data, assuming a true $\Lambda$CDM universe ($\Omo=0.30, \oo=-1, \oa=0$), yields biased estimates such as $\oo \approx -0.82 \pm 0.06$ and $\oa \approx -0.82 \pm 0.4$. This systematic shift, driven by the $\Omo$ prior, moves the estimated parameters towards the non-$\Lambda$CDM region, offering a qualitative resemblance to outcomes reported in current combined DESI DR2 BAO + Pantheon+ SNe analyses (e.g., $\oo = -0.888^{+0.055}_{-0.064}$, $\oa = -0.17 \pm 0.46$). Our findings suggest that these observed non-$\Lambda$CDM parameters may largely arise from statistical biases due to $\Omo$ prior tensions between datasets. This study proposes a potential resolution to current cosmological tensions without necessarily invoking new physics.

Optical observations with high time resolution are essential for understanding the origin of sub-millisecond timescale astronomical phenomena, including giant radio pulses from the Crab Pulsar. We have developed a high-speed imaging system called the Imager of MPPC-based Optical photoN counter from Yamagata (IMONY). The system uses a customized Multi-Pixel Photon Counter (MPPC), which independently reads out signals from all $8\times8$ pixels and functions as an imager based on a Geiger-mode avalanche photodiode array. This system assigns timestamps to detected photons with a time resolution of 100 ns. We installed IMONY on the 3.8 m aperture Seimei Telescope in Okayama, Japan. We have successfully detected the 34-ms period of the Crab Pulsar and imaged stars in the sensor's field of view. However, we have also found that a small fraction of the pixels have shown double or multiple pulses that are used for photon arrival timing. This situation is likely due to circuit noise and may unfortunately result in overestimating the number of photons detected. In order to precisely estimate the photon flux of targets or the sky background, calibration of such over-counts is important. We measured the number of detected photons relative to the light intensity of each pixel in a laboratory environment. We estimated the number of spurious hit pulses caused by signal tail fluctuations exceeding the comparator threshold, based on the exponential distribution of time intervals between pulses. These are distinct from typical SiPM afterpulses and originate from electronic effects in our readout system. After applying the calibration to the observed data, we confirmed the linearity between the V-band magnitudes of stars and the number of detected photons.

Extragalactic carbon monoxide (CO) line emission will likely be an important signal in current and future Cosmic Microwave Background (CMB) surveys on small scales. However, great uncertainty surrounds our current understanding of CO emission. We investigate the implications of this modeling uncertainty on CMB surveys. Using a range of star formation rate and luminosity relations, we generate a suite of CO simulations across cosmic time, together with the broadband cosmic infrared background (CIB). From these, we quantify the power spectrum signatures of CO that we would observe in a CMB experiment at 90, 150, and 220 GHz. We find that the resulting range of CO auto-spectra spans up to two orders of magnitude and that while CO on its own is unlikely to be detectable in current CMB experiments, its cross-correlation with the CIB will be a significant CMB foreground in future surveys. We then forecast the bias on CMB foregrounds that would result if CO were neglected in a CMB power spectrum analysis, finding shifts that can be comparable to some of the uncertainties on CMB foreground constraints from recent surveys, particularly for the thermal and kinetic Sunyaev-Zel'dovich effects and radio sources, and many times greater than the expected uncertainties expected for future data. Finally, we assess how the broad range of multifrequency CO$\times$CIB spectra we obtain is captured by a reduced parameter set by performing a principal component analysis, finding that three amplitude parameters suffice for a CMB-S4-like survey. Our results demonstrate the importance for future CMB experiments to account for a wide range of CO modeling, and that high-precision CMB experiments may help constrain extragalactic CO models.

Light unaffiliated mass clumps (LUMCs), i.e. dark matter (DM) components without any stellar counterparts, have been reported in strong lensing mass reconstructions of MACS 0416, both on galaxy and galaxy cluster scales. On galaxy cluster scale, the most recent parametric study based on 303 multiple images features a LUMC, in the south of the cluster. On galaxy scale, the most recent GRALE non-parametric study based on 237 multiple images features two LUMCs, M1 and M2. Given the implications of these findings in the context of structure formation and evolution, we test these features parametrically, using LENSTOOL. First, we show that a mass model where each large scale DM component introduced in the modelling is associated with a stellar counterpart can reproduce the 303 multiple images, removing the need for any cluster scale LUMC in \lens. We then update the GRALE model using the 303 multiple images, finding that one of the two galaxy scale LUMC, M1, is no longer significant, while M2 remains. We test M2 by explicitely including it in our parametric model. We find that the inclusion of this LUMC does not improve the global RMS, but mildly improves locally the RMS for one multiple image located close to M2. Besides, the preferred mass for M2 corresponds to the lowest mass allowed by the adopted prior. M2 is therefore not strongly favoured by a parametric approach but it is not ruled out. We present a comparison of parametric and non-parametric models in the M2 area. Both approaches show very similar surface mass density at this location, with a 5-6 % difference between the mass maps. The difference is that GRALE favors a distinct mass substructure when LENSTOOL favors a more diffuse mass distribution. We have been able to propose a parametric mass model without including any LUMCs, providing further evidences for DM being associated with light in clusters.

Atmospheric observations by JWST raise a growing evidence that atmospheric metallicity exhibits an anti-correlation with masses of giant exoplanets. While such a trend was anticipated by planetesimal-based planet formation models, it remains unclear what kind of atmospheric metallicity trends emerge from pebble-based planet formation. Moreover, while recent studies of solar system Jupiter suggest that uppermost observable atmosphere may not represent the bulk envelope composition, it remains uncertain how the envelope inhomogeneity influences the atmospheric metallicity trend. In this study, we develop disk evolution and planet formation models to investigate the possible atmospheric metallicity trends of giant exoplanets formed via pebble accretion and how they depend on the metallicity inhomogeneity within the envelope. We find that pebble-based planet formation produces two distinct mass-metallicity relations depending on planetary birthplace. Planets formed beyond the H2O snowline exhibit a mass-metallicity anti-correlation similar to that predicted by planetesimal-based models if their atmospheres are fully convective. This anti-correlation disappears if the convective mixing is inefficient. In contrast, planets formed inside the H2O snowline show a shallower mass-metallicity anti-correlation, regardless of the efficiency of atmospheric mixing. Many gas giants observed by JWST observations lie around the mass-metallicity relation predicted for formation at close-in orbits, although some planets with sub-stellar atmospheric metallicity appear to require unmixed envelopes and formation beyond the H2O snowline. We also examine the relationship between bulk and atmospheric metallicity and find a clear correlation that closely follows atmospheric metallicity that is comparable to bulk metallicity.

6000+ exoplanets are now detected using various techniques, with estimates of billions of planets existing in our Galaxy alone. They are called super-Earths, hot Earths, mini-Neptunes, hot Neptunes, sub-Neptunes, Saturns, Jupiters, hot Jupiters, gas giants, ice giants, rocky, terran, subterran, superterran, and so on. This prompted the emergence of many recent works on taxonomy, or classification, of exoplanets. However, there is still no basic, fundamental definition of 'What is a planet?'. IAU has forwarded a definition in 2006 which, however, raised more questions than it solved. The first ambitious task here is to establish if there are limits on the size/mass of planets. The lower mass limit may be assumed as of Mimas - approximately minimum mass required to attain a nearly spherical hydrostatic equilibrium shape. The upper mass limit may be easier - there is a natural lower limit to what constitutes a star: ~0.08 SU. But then there are brown dwarfs: IAU has defined brown dwarfs as objects that exceed the deuterium burning limit (~13 JU), and giant exoplanets generally have masses of ~0.3 to ~60 JU. The resolution requires assembling the basic physical parameters that define planets quantitatively. Mass and radius are the two fundamental properties, and we propose to use a third correlated parameter: the moment of inertia. Based on this, we create the fundamental planet plane where the two parameters are correlated with the third, similar to a fundamental galactic plane. The fundamental planetary plane with turn-off point diagrams is constructed for the first time in this work. We define a planet as 'A celestial spherical object, bound to a star or unbound, that lies on the fundamental planetary plane, within a mass range between 0.02 EU to 13 JU'.

While massive stars are known to shape galactic ecosystems, their formation has long been assumed to require the high-density environments of inner galactic disks. This paradigm is challenged by mounting evidence of young massive stars in extended galaxy outskirts, yet direct confirmation of in situ massive star formation in such extreme low-density environments remains scarce. Here, we present the discovery of LAMOST J0048+4154, a massive yellow supergiant situated at a deprojected galactocentric distance of ~34 kpc in M31, making it the most distant massive star confirmed in this galaxy. Through spectroscopic and photometric analyses, we classify J0048+4154 as an F5-F8I supergiant with an effective temperature of $6357^{+121}_{-118}$ K and a luminosity of $\log L/L_{\odot} = 5.00^{+0.06}_{-0.06}$, corresponding to an ~18 $M_{\odot}$ progenitor and an age of ~10 Myr. FAST H I observations reveal close spatial and kinematic alignment between the star and a faint H I external arm, suggesting in situ formation in a region of low gas density. The presence of other UV-bright, early-type stars in the vicinity further supports low-level recent star formation in M31's very outer disk. These findings challenge the prevailing assumption that massive star formation is confined to inner disks or classical star-forming regions and underscore the need to re-examine the role of spiral galaxy outskirts in fueling and sustaining star formation. J0048+4154 thereby expands our understanding of the extent of M31's young stellar component and exemplifies how outer disks may harbor conditions conducive to forming massive stars, despite low-density environments.

The formation and evolution of ring structures in galaxies are crucial for understanding the nature and distribution of dark matter, galactic interactions, and the internal secular evolution of galaxies. However, the limited number of existing ring galaxy catalogs has constrained deeper exploration in this field. To address this gap, we introduce a two-stage binary classification model based on the Swin Transformer architecture to identify ring galaxies from the DESI Legacy Imaging Surveys. This model first selects potential candidates and then refines them in a second stage to improve classification accuracy. During model training, we investigated the impact of imbalanced datasets on the performance of the two-stage model. We experimented with various model combinations applied to the datasets of the DESI Legacy Imaging Surveys DR9, processing a total of 573,668 images with redshifts ranging from z_spec = 0.01-0.20 and magr <17.5. After applying the two-stage filtering and conducting visual inspections, the overall Precision of the models exceeded 64.87%, successfully identifying a total of 8052 newly discovered ring galaxies. With our catalog, the forthcoming spectroscopic data from DESI will facilitate a more comprehensive investigation into the formation and evolution of ring galaxies.

Fireballs (bolides) are high-energy luminous phenomena produced when meteoroids and small asteroids enter Earth's atmosphere at hypersonic speeds, often resulting in fragmentation or complete disintegration accompanied by significant energy release. The resulting bolide light curves capture temporal brightness variations as these objects traverse increasingly dense atmospheric layers, providing essential information on meteoroid entry dynamics, fragmentation behavior, and atmospheric energy deposition processes. The Center for Near-Earth Object Studies' (CNEOS) continuously expanding fireball database offers a globally comprehensive archive of bolide events, including light curves and associated metadata. Events associated with infrasound detections allow direct correlations between acoustic signatures and light-curve features, therefore enabling detailed analyses of fragmentation dynamics and energy deposition. Here, we introduce BLADE (Bolide Light-curve Analysis and Discrimination Explorer), a robust and high-fidelity framework specifically designed to analyze bolide light curves for objects detected from space. BLADE incorporates a processing pipeline integrating Savitzky-Golay filtering, prominence-based peak detection, and gradient analysis, enabling systematic identification and classification of fragmentation events and their associated energy release characteristics. Preliminary results demonstrate that BLADE reliably distinguishes distinct bolide behaviors, providing an objective, scalable methodology for characterization and analysis of large bolide light curve datasets. This foundational work establishes a novel pathway for advanced bolide research, with promising applications in planetary defense and global atmospheric monitoring.

Stellar parameters and abundances provide crucial insights into stellar and Galactic evolution studies. In this work, we developed a convolutional neural network (CNN) to estimate stellar parameters: effective temperature ($T_{\text{eff}}$), surface gravity (log $g$) and metallicity (both [Fe/H] and [M/H]) as well as six $\alpha$-elements (C, N, O, Mg, Si, Ca) and [$\alpha$/M]. We selected giant stars with \( 3500 \, \text{K} < T_{\text{eff}} < 5500 \, \text{K} \) and \( 0 \, \text{dex} < \log g < 3.6 \, \text{dex} \) from the LAMOST and APOGEE surveys, while requiring (S/N)$_g$ of the LAMOST low-resolution spectra $>$ 10, which leaves 1,100,858 giant stars. The spectral from LAMOST and the labels from APOGEE for 62,511 common stars were used as our training set. The corresponding test set yields scatters 50 K, 0.06 dex and 0.13 dex for $T_{\text{eff}}$, [Fe/H] and log $g$, respectively. For $\alpha$ elements O, Mg, Si and Ca, the scatters are 0.05 dex, 0.04 dex, 0.03 and 0.04 dex, respectively. For C and N elements, the scatters are 0.07 dex and 0.05 dex. For [$\alpha$/M] and [M/H], the scatters are 0.03 dex and 0.06 dex. The mean absolute error (MAE) of most elements are between 0.02 $-$ 0.04 dex. The predicted abundances were cross-matched with previously identified substructures PG1 and PG2, with their origins subsequently analyzed. Finally, the catalog is available at this https URL.

The Evolutionary Map of the Universe (EMU) survey with ASKAP is transforming our understanding of radio galaxies, AGN duty cycles, and cosmic structure. EMUCAT efficiently identifies compact radio sources, yet struggles with extended objects, requiring alternative approaches. The Radio Galaxy Zoo: EMU (RGZ EMU) project proposes a general framework that combines citizen science and machine learning to identify around 4 million extended sources in EMU. This framework is expected to enhance the EMUCAT cataloging on extended sources and can be further empowered with the introduction of cross-matched external data from surveys such as POSSUM and WALLABY.

SMBHs are theorised to undergo significant growth in the early Universe, however, the X-ray Luminosity Function (XLF), used as a principal tracer of the SMBH accretion density, lacks observational constraints at z>6, until now. We present new measurements of the z=4-10 XLF at intermediate luminosities, taking advantage of recent deep near-IR imaging from UltraVISTA that enables us to identify galaxies and AGN at high redshifts within which we identify X-ray sources using Chandra COSMOS data. We first performed a cross-match to a deep Chandra source list, for which the X-ray sensitivity can be accurately quantified, before exploiting available X-ray data further through direct extraction of X-ray counts at the positions of COSMOS2020 galaxies. With the resulting z=4-10 X-ray AGN sample, comprised of 21 blind detections and 11 directly extracted detections, we have measured the early space density of AGN, at moderate-luminosities where the majority of early SMBH growth occurred. These measurements reveal higher space-densities than expected, based on the extrapolation of XLF models from lower redshifts. Whilst our measured space densities at z=4-5 are consistent with model predictions, at z=5-7 we find space densities of the order of 10$\times$ the extrapolated model predictions and could be as high as 220$\times$ the model extrapolations at z=7-10. In addition, we find evidence that a large fraction of the early AGN population are heavily obscured, with an obscured fraction of 0.982$^{+0.007}_{-0.008}$; correcting for this obscuration further increases the measured space densities. Comparing to recent JWST results, these measurements begin to bridge the gap between the bright-end of the quasar luminosity function and the latest JWST observations of very early, low-luminosity AGN, indicating a larger fraction of the first galaxies play host to rapidly growing SMBH than previously thought.

We present high-cadence photometric and low-resolution (R $\sim$ 400--700) optical spectroscopic observations of Type IIP supernova, SN~2018pq, which exploded on the outskirts of the galaxy IC~3896A. The optically thick phase (``plateau'') lasts approximately 97 d, the plateau duration of normal Type IIP supernovae. SN~2018pq has a {\em V}-band absolute magnitude of $-16.42 \pm 0.01$ mag at 50 d, resembles normal-luminous supernova, and the V-band decline rate of 0.42$\pm$0.06 mag 50 d$^{-1}$ during the plateau phase. A steeper decline rate of 11.87$\pm$1.68 mag 100 d$^{-1}$ was observed compared to that of typical Type IIP supernovae during the transition between plateau to nebular phase. We employ detailed radiative transfer spectra modelling, TARDIS, to reveal the photospheric temperature and velocity at two spectral epochs. The well-fitted model spectra indicate SN~2018pq is a spectroscopically normal Type IIP supernova. Semi-analytical light curve modelling suggests the progenitor as a red supergiant star with an ejecta mass of $\sim$11 $M_\odot$ and an initial radius of 424 $R_\odot$. On the contrary, hydrodynamical modelling suggests a higher mass progenitor between 14--16 $M_\odot$.

This paper presents the development and application of SpectrumMate LR, a broadband, low-resolution spectrograph for small telescope use. SpectrumMate LR is designed to offer affordable, accessible spectroscopic capabilities for amateur astronomers and educators, inspired by the need for versatile instruments in amateur and educational settings. Utilising a 300 grooves/mm grating, 80 mm collimator and objective lenses, SpectrumMate LR is optimised to perform analyses across the visible spectrum, enabling users to classify stars by spectral type, measure stellar temperatures, and test filter transmission ranges. Tests demonstrate SpectrumMate LR's ability to capture accurate spectral data, validating its efficacy in observing both celestial and terrestrial light sources. This instrument fills a niche for cost-effective spectroscopy, empowering a broader audience to engage in detailed observational astronomy

We investigate constraints on the least explored, smallest mass scales of primordial black holes (PBHs), which evaporate prior to Big Bang Nucleosynthesis (BBN). Our study examines the impact of Planck-mass relics on the allowed fraction of dark matter composed of PBHs ($f_{PBH}$), as well as on the resulting stochastic gravitational wave background and the formation of primordial binaries. We discuss how these binaries and early mergers lead to longer PBH lifetimes, extending the reach of the stringent BBN constraints to smaller masses than usually expected. We make comprehensive constraint plots on the collapse fraction $\beta$ and $f_{PBH}$ (including relics), focusing on ultra-light PBHs.

Tidal disruption events (TDEs) are excellent tools for probing low mass supermassive black holes (SMBHs) that may otherwise remain undetected. Here, we present an extended SMBH--Bulge mass scaling relationship using these lower mass TDE black holes and their host galaxies. Bulge masses are derived using Prospector fits to UV-MIR spectral energy distributions for the hosts of 40 TDEs that have a detected late-time UV/optical plateau emission, from which a SMBH mass is derived. Overall, we find that TDE plateaus are a successful method for probing BH scaling relations. We combine the observed TDE sample with a higher mass SMBH sample and extend the known relationship, recovering a steeper slope ($m = 1.34 \pm 0.03$) than current literature estimates, which focus on the high mass regime. For the TDE only sample, we measure an equally significant but shallower relationship with a power-law slope of $m = 1.17 \pm 0.10$ and significance $<0.001$. Forward modelling is used to determine whether known selection effects can explain both the comparatively flatter TDE only relation and the overall steepening across the full SMBH mass range. We find that the flattening at TDE masses can be accounted for, however the steepening can not. It appears that if a single slope extends for the whole BH mass range, it must be steeper to include the TDE population.

Recent analyses from the Pierre Auger Collaboration suggest correlations between the arrival directions of Ultra-High-Energy Cosmic Rays (UHECRs) and catalogs of starburst galaxies (SGBs) and jetted active galactic nuclei (AGNs). We revisit these analyses using the same methodology as \auger , but explicitly incorporating UHECR deflections in turbulent extragalactic magnetic fields (EGMFs). We demonstrate that while for SBGs the same sources as for the generic \auger\ analysis dominate the catalog correlations, jetted AGNs are dominated by Centaurus~A when accounting for source distances and deflections. Using our framework, we derive 90\% confidence level upper limits on the local EGMF strength of 4.4~nG~Mpc$^{1/2}$ for SBGs and 6.7~nG~Mpc$^{1/2}$ for jetted AGNs. Assuming instead that the UHECR deflections predominantly arise from the Galactic magnetic field (GMF), we obtain a GMF upper limit of $1.4 \, \mu$G~kpc$^{1/2}$ for a Galactic halo size of 30~kpc.

Coronal Mass Ejections (CMEs) are space weather phenomena capable of causing significant disruptions to both space- and ground-based infrastructure. The timely and accurate detection and prediction of CMEs is a crucial steps towards implementing strategies to minimize the impacts of such events. CMEs are commonly observed using coronagraphs and heliospheric imagers (HIs), with some forecasting methods relying on manually tracking CMEs across successive images in order to provide an estimate of their arrival time and speed. This process is time-consuming, and the growing volume of available data makes manual identification of CMEs increasingly impractical. We investigate the application of machine learning (ML) techniques to the problem of automated CME detection, focusing on data from the HI instruments aboard the STEREO spacecraft. HI data facilitates the tracking of CMEs through interplanetary space, providing valuable information on their evolution along the way. Building on advances in image segmentation, we present the Solar Transient Recognition Using Deep Learning (STRUDL) model. STRUDL is designed to automatically detect and segment CME fronts in HI data. We address the challenges inherent to this task and evaluate the model's performance across a range of solar activity conditions. To complement segmentation, we implement a basic tracking algorithm that links CME detections across successive frames, thus allowing us to automatically generate time-distance profiles for all CMEs under study. Our results demonstrate the feasibility of applying ML-based segmentation techniques to HI data, while highlighting areas for future improvement, particularly regarding the accurate segmentation and tracking of faint and interacting CMEs.

We present the first band-2 (120 - 250 MHz) uGMRT (upgraded Giant Metrewave Radio Telescope) observations of the bimodal galaxy cluster Abell 1644 (z = 0.0471), complemented by Chandra X-ray data. While weak lensing measurements reveal a third substructure in Abell 1644, our radio analysis reveals only two compact sources coinciding with the respective brightest cluster galaxies (BCGs) of the northern (A1644N) and southern (A1644S) substructures, seen in the X-ray observations. Radio analysis yields compact active galactic nuclei (AGN) powered sources with radio power $P_{A1644S} = 1.1\times 10^{23} W/Hz$ and $P_{A1644N} = 7.3\times 10^{23} W/Hz$ at 200MHz. We find no evidence of non-thermal diffuse radio emission, such as halos or relics, within the sensitivity of our band-2 image. We measured the flux density of each radio source and performed spectral analysis. A1644N exhibits a synchrotron power law spectrum while A1644S shows spectral turnover suggestive of synchrotron self-absorption. X-ray analysis reveals two shock fronts at the southern substructure with a Mach number $M= 3.21 \pm 0.51$ and $M= 2.22 \pm 0.07$, indicating an ongoing merger. Our findings reinforce the complex dynamical nature of Abell 1644 and contribute to a deeper understanding of the cluster's thermodynamic state. Future deep radio observations with improved radio frequency interference (RFI) mitigation will be crucial for probing non-thermal phenomena in this system.

HD 209458b is the prototypical hot Jupiter and one of the best targets available for precise atmosphere characterisation. Now that spectra from both Hubble Space Telescope (HST) and James Webb Space Telescope (JWST) are available, we can reveal the atmospheric properties in unprecedented detail. In this study, we perform a new data reduction and analysis of the original HST/WFC3 spectrum, accounting for the wavelength dependence of the instrument systematics that was not considered in previous analyses. This allows us to precisely and robustly measure the much-debated H$_2$O abundance in HD 209458b's atmosphere. We combine the newly reduced spectrum with archival JWST/NIRCam data and run free chemistry atmospheric retrievals over the 1.0 - 5.1 $\mu$m wavelength range, covering possible features of multiple absorbing species, including CO$_2$, CO, CH$_4$, NH$_3$, HCN, Na, SO$_2$, and H$_2$S. We detect H$_2$O and CO$_2$ robustly at above 7 $\sigma$ significance, and find a 3.6 $\sigma$ preference for cloudy models compared to a clear atmosphere. For all other absorbers we tested, only upper limits of abundance can be measured. We use Bayesian model averaging to account for a range of different assumptions about the cloud properties, resulting in a water volume mixing ratio of $0.95^{+0.35}_{-0.17} \:\times$ solar and a carbon dioxide abundance of $0.94^{+0.16}_{-0.09} \:\times$ solar. Both results are consistent with solar values and comparable to predictions from the VULCAN 1D photochemistry model. Combining these values with a prior on the CO abundance from ground-based measurements, we derive an overall atmospheric composition comparable to solar metallicity of $\mathrm{[M/H]} = 0.10^{+0.41}_{-0.40}$ and very low C/O of $0.054^{+0.080}_{-0.034}$ with a 3 $\sigma$ upper limit of 0.454. This indicates a strong enrichment in oxygen and depletion in carbon during HD 209458b's formation.

Observational astronomy relies on visual feature identification to detect critical astrophysical phenomena. While machine learning (ML) increasingly automates this process, models often struggle with generalization in large-scale surveys due to the limited representativeness of labeled datasets -- whether from simulations or human annotation -- a challenge pronounced for rare yet scientifically valuable objects. To address this, we propose a conditional diffusion model to synthesize realistic galaxy images for augmenting ML training data. Leveraging the Galaxy Zoo 2 dataset which contains visual feature -- galaxy image pairs from volunteer annotation, we demonstrate that our model generates diverse, high-fidelity galaxy images closely adhere to the specified morphological feature conditions. Moreover, this model enables generative extrapolation to project well-annotated data into unseen domains and advancing rare object detection. Integrating synthesized images into ML pipelines improves performance in standard morphology classification, boosting completeness and purity by up to 30\% across key metrics. For rare object detection, using early-type galaxies with prominent dust lane features ( $\sim$0.1\% in GZ2 dataset) as a test case, our approach doubled the number of detected instances from 352 to 872, compared to previous studies based on visual inspection. This study highlights the power of generative models to bridge gaps between scarce labeled data and the vast, uncharted parameter space of observational astronomy and sheds insight for future astrophysical foundation model developments. Our project homepage is available at this https URL.

To investigate the interstellar medium (ISM) and Galactic cosmic rays (CRs) in the solar neighborhood, we analyzed ${\gamma}$-ray data by Fermi Large Area Telescope (LAT) for five nearby molecular cloud regions. Our data includes the MBM/Pegasus region (MBM~53, 54, 55 clouds and Pegasus loop), R CrA region (R Coronae Australis clouds), Chamaeleon region (Chamaeleon clouds), Cep/Pol region (Cepheus and Polaris flare), and Orion region (Orion clouds). The ISM templates are constructed by a component decomposition of the 21~cm {\HI} line, the Planck dust emission model, and the carbon monoxide (CO) 2.6~mm line. Through $\gamma$-ray data analysis the ISM gas is successfully decomposed into non-local {\HI}, narrow-line and optically thick {\HI}, broad-line and optically thin {\HI}, CO-bright {\Htwo}, and CO-dark {\Htwo} for all five regions. CR intensities evaluated by the ${\gamma}$-ray emissivity of broad {\HI} agree well with a model based on directly-measured CR spectra at the Earth, with a gradient giving a higher CR intensity toward the inner Galaxy at the 10\% level in ${\sim}$ 500~pc. The ratio of CO-dark {\Htwo} to CO-bright {\Htwo} anti-correlates with the {\Htwo} mass traced by the CO 2.6~mm line, and reaches 5--10 for small systems of ${\sim}$1000 solar mass.

Conventional galaxy generation methods rely on semi-analytical models and hydrodynamic simulations, which are highly dependent on physical assumptions and parameter tuning. In contrast, data-driven generative models do not have explicit physical parameters pre-determined, and instead learn them efficiently from observational data, making them alternative solutions to galaxy generation. Among these, diffusion models outperform Variational Autoencoders (VAEs) and Generative Adversarial Networks (GANs) in quality and diversity. Leveraging physical prior knowledge to these models can further enhance their capabilities. In this work, we present GalCatDiff, the first framework in astronomy to leverage both galaxy image features and astrophysical properties in the network design of diffusion models. GalCatDiff incorporates an enhanced U-Net and a novel block entitled Astro-RAB (Residual Attention Block), which dynamically combines attention mechanisms with convolution operations to ensure global consistency and local feature fidelity. Moreover, GalCatDiff uses category embeddings for class-specific galaxy generation, avoiding the high computational costs of training separate models for each category. Our experimental results demonstrate that GalCatDiff significantly outperforms existing methods in terms of the consistency of sample color and size distributions, and the generated galaxies are both visually realistic and physically consistent. This framework will enhance the reliability of galaxy simulations and can potentially serve as a data augmentor to support future galaxy classification algorithm development.

The Giant Radio Array for Neutrino Detection (GRAND) is an envisioned large-scale radio array designed to detect ultra-high-energy cosmic rays (UHECRs, $E > 100$ PeV) and neutrinos. Employing cost-effective antennas distributed across vast areas, GRAND is optimized to observe the rare flux of ultra-high-energy particles with high precision. The GRANDProto300 (GP300) pathfinder array, currently under deployment, targets the $10^{16.5} - 10^{18}$ eV range and is anticipated to achieve approximately 15\% energy resolution and 20g/cm$^2$ $X_{\mathrm{max}}$ precision. This level of precision enables accurate measurements of the fine structure of the energy spectrum, mean logarithmic mass ($\langle \ln A \rangle$), and proton flux within this range. After five years of data collection, the sensitivity for detecting anisotropy could reach $5 \times 10^{-3}$ for energies below $10^{17.1}$ eV. With its substantially larger effective area, GRAND extends these capabilities to the highest energies ($\sim 10^{20}$ eV), offering enhanced statistics and sensitivity for spectral, composition, and anisotropy measurements within one year for UHECRs.

Studying the radio spectral energy distribution (SED) of distant galaxies is essential for understanding their assembly and evolution over cosmic time. We present rest-frame radio SEDs of a sample of 160 starburst galaxies at redshifts 1.5 to 3.5 in the COSMOS field, as part of the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) project. MeerKAT observations, combined with archival VLA and GMRT data, allow us to determine the integrated mid-radio (1-10 GHz) continuum (MRC) luminosity and magnetic field strength. A Bayesian method is used to model the SEDs and separate free-free and synchrotron emission. We calibrate the star formation rate (SFR) in radio both directly through SED analysis and indirectly via the infrared-radio correlation (IRRC). With a mean synchrotron spectral index of approximately 0.7, we find that the index flattens with redshift and specific SFR, suggesting that cosmic rays are more energetic in the early universe due to higher star formation activity. The magnetic field strength increases with redshift (B is proportional to (1 + z)^0.7) and with star formation rate (B is proportional to SFR^0.3), indicating a small-scale dynamo as the dominant amplification mechanism. Accounting for SED evolution, the IRRC remains redshift-invariant and does not vary with stellar mass at 1.5 < z < 3.5, though the correlation deviates from linearity. Similarly, we show that SFR estimates based on integrated MRC luminosity are also redshift-invariant.

Recent studies have highlighted the potential of intracluster light (ICL) as a dark matter tracer. Moreover, the ICL co-evolves with the brightest cluster galaxy (BCG) and the host cluster, making it a valuable tool for understanding cluster dynamics. In this study, we utilize 426 galaxy clusters (with total mass $M_{\rm tot}>10^{14} M_{\odot}$ at $z=0$) simulated in the cosmological hydrodynamical simulation Illustris TNG300 to compare the spatial distributions of dark matter, member galaxies, gas, and ICL and assess their effectiveness as dark matter tracers in the central regions of clusters at $R_{\rm vir}<0.3$. We apply the Weighted Overlap Coefficient (WOC), a methodology for quantifying the similarity of two-dimensional spatial distributions, to various components of the galaxy clusters at different dynamical stages. Our findings reveal that the spatial distributions of both ICL combined with the BCG and gas closely resemble the dark matter distribution, with higher fidelity observed in more relaxed galaxy clusters with earlier half-mass epochs. These results demonstrate that the BCG+ICL component serves as an effective tracer of dark matter, consistent with previous observational studies linking cluster light to mass. Moreover, the degree of spatial similarity between the BCG+ICL and dark matter distributions appears to reflect the dynamical state of the cluster.

In this study, we present a comprehensive analysis of transit timing variations (TTVs) in the ultra-short-period gas giant WASP-19b, which orbits a G-type main-sequence star. Our analysis is based on a dataset comprising 204 transit light curves obtained from the Transiting Exoplanet Survey Satellite (TESS), the Exoplanet Transit Database (ETD), and the ExoClock project, supplemented by 18 publicly available light curves. Mid-transit times were extracted from these data, and an additional 98 mid-transit times compiled from the literature were incorporated, resulting in a combined dataset spanning approximately 14 years. After excluding light curves significantly impacted by stellar activity, such as starspot anomalies, the final dataset consisted of 252 high-quality mid-transit times. Initial inspection of the transit timing residuals using an apsidal precession model suggested the possible presence of an additional planetary companion. However, subsequent frequency analysis and sinusoidal model fitting indicate that the observed TTVs are more consistently explained by apsidal precession of WASP-19b's orbit. We also considered alternative mechanisms, including the Applegate mechanism and the Shklovskii effect. Our findings suggest that stellar magnetic activity, potentially linked to the Applegate mechanism, may also contribute to the observed timing variations. To further constrain the origin of the TTVs and assess the contributions of these mechanisms, continued high-precision photometric monitoring of the WASP-19 system is strongly recommended.

High-energy SETI pushes astrobiology to its limits, testing the most fundamental needs of life and the most extreme limits of technology. It has lagged behind the rest of the field, but the increased respectability of SETI could spark interest in the coming years. This white paper reviews the case for SETI in X-rays, gamma rays, and neutrinos, including rationales, challenges, and possible technosignatures, and suggests future strategies for observational work.

We study the energy-dependent pulse profile of 4U 1538-52 and its phase-dependent spectral variability, with emphasis on the behavior around the cyclotron resonant scattering feature at around 21 keV. We analyze all available NuSTAR observations of 4U 1538-52. We decompose energy-resolved pulse profiles into Fourier harmonics to study their energy dependence. Specifically, we compute pulsed fraction spectra, cross-correlation and lag spectra, identifying discontinuities and linking them to features in the phase-averaged spectra. We perform phase-averaged and phase-resolved spectral analyses to probe spectral variability and its relation to pulse profile changes. Finally, we interpret our findings via physical modeling of energy- and angle-dependent pulse profile emission, performing radiative transfer in a homogeneous slab-like atmosphere under conditions relevant to 4U 1538-52. The emission is projected onto the observer's sky plane to derive expected observables. In contrast to the dips in pulsed fraction spectra observed in other sources (e.g., Her X-1), we find a broad bump near the cyclotron resonance energy in 4U 1538-52. This increase is driven primarily by phase-dependent spectral variability, especially by strong variations in cyclotron line depth across different phase intervals. We interpret the observed contrast between dips and bumps in various sources as arising from phase-dependent variations of cyclotron line depth relative to the phase-modulated flux. We model the X-ray emission from an accreting neutron star and find that our simulations indicate high values of both the observer's inclination and the magnetic obliquity, along with a 10-15 degrees asymmetry between the locations of the magnetic poles. Assuming this geometry, we satisfactorily reproduce the observed pulse profiles and introduce general trends in the observables resulting from the system's geometry.

Forecasting cosmological constraints from halo-based statistics often suffers from instability in derivative estimates, especially when the number of simulations is limited. This instability reduces the reliability of Fisher forecasts and machine learning based approaches that use derivatives. We introduce a general framework that addresses this challenge by stabilizing the input statistic and then systematically identifying the optimal subset of summary statistics that maximizes cosmological information while simultaneously minimizing the instability of predicted constraints. We demonstrate this framework using the Voronoi volume function (VVF), a summary statistic that captures beyond two-point clustering information. Applying our two-step procedure -- random sub-sampling followed by optimization -- improves the constraining power by up to a factor of 2, while also enhancing the stability of the forecasts across realizations. As surveys like Euclid, DESI, and LSST push toward tighter constraints, the ability to produce stable and accurate theoretical predictions is essential. Our results suggest that new summary statistics such as the VVF, combined with careful data curation and stabilization strategies, can play a key role in next-generation precision cosmology.

Weak lensing surveys are often summarized by constraints on the derived parameter ${S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}}$, obscuring the rich scale and redshift information encoded in the data, and limiting our ability to identify the origin of any tensions with $\Lambda$CDM predictions from the cosmic microwave background. In this work, we introduce a fast and flexible framework to extract the scale-dependent matter power spectrum $P(k, z)$ from cosmic shear and CMB lensing measurements, parameterizing deviations from the Planck $\Lambda$CDM prediction as a free function $\alpha(k)$. Using public data from DES Y3, KiDS-1000, HSC Y3, and ACT DR6, we constrain $\alpha(k)$ with fast Hamiltonian Monte Carlo inference. Our results show a consistent 15-30\% suppression in the matter power spectrum at intermediate scales ($k \sim 0.1-1{\rm Mpc}^{-1}$) in galaxy-lensing data relative to Planck, with combined tensions reaching up to $4\sigma$ at fixed cosmology. In contrast, ACT CMB lensing is consistent with $\Lambda$CDM at ${k\lesssim 0.1 {\rm Mpc}^{-1}}$. We validate our method using mock data, quantify consistency between datasets, and demonstrate how the resulting $\alpha(k)$ likelihoods can be used to test specific models for the power spectrum. All code, data products, and derived likelihoods are publicly released. Our results highlight the importance of reporting lensing constraints on $P(k, z)$ and pave the way for model-agnostic test of growth of structure with upcoming surveys such as LSST, Euclid, and Roman.

Star formation in ultra-faint dwarf galaxies (UFDs, $M_* <10^5M_\odot$) is suppressed by reionization, but may not be completely quenched. The metallicity distribution function (MDF) of stars in ultra-faint dwarf galaxies could show these signatures of reionization. However, past studies of UFD MDFs have been limited, because there are only a few dozen red giant branch (RGB) stars in such low-mass galaxies. We present low-resolution Magellan/IMACS spectroscopy of 167 stars in the UFD Reticulum II ($M_* \approx 3000 M_\odot$), increasing the number of stellar metallicities by 6.5 times and resulting in the most populated spectroscopic metallicity distribution function of any UFD. This is possible because we determined the first spectroscopic metallicities of main sequence turn-off stars in any UFD. The MDF of Reticulum II is clearly a bimodal distribution, displaying two peaks with about $80\%$ of the stars in the metal-poor peak at $\rm[Fe/H]=-3.0$ and $20\%$ of the stars in the more metal-rich peak at $\rm[Fe/H]=-2.1$. Such a large metallicity gap can be explained by Type Ia supernova enrichment during a long quiescent period. This supports the currently-favored two-burst star formation history for Reticulum II and shows that such low-mass galaxies clearly can form stars after reionization.

We report the serendipitous detection of the SO $J_N=6_5-5_4$ (219.949 GHz) rotational transition in archival Atacama Large Millimeter/submillimeter Array (ALMA) observations of the spiral hosting protoplanetary disks around CQ Tau (with $\approx4.9\sigma$ significance) and MWC 758 (with $\approx3.4\sigma$ significance). In the former, the SO emission comes in the shape of a ring, arises from the edge of the continuum cavity, and is qualitatively consistent, at the currently available spectral resolution, with being in Keplerian rotation. In the latter, instead, while arising primarily from inside the continuum cavity, the SO emission also extends to the continuum ring(s), and its morphology and kinematics are less clear. We put these sources in the context of the other protoplanetary disks where SO detections have been previously reported in the literature and discuss the possible origins of SO in terms of (thermal) desorption or formation in the gas phase. We argue that these processes might be fostered by dynamical perturbations caused by unseen embedded massive companions, shadows, or late-time infall, thus suggesting a possible link between perturbed dynamics and SO emission in (these) protoplanetary disks. If confirmed, our interpretation would imply that chemical evolution timescales could be significantly shorter in these systems than is commonly assumed, indicating that dynamical perturbations might influence the composition of newborn (proto-)planets by altering the volatile makeup of their formation environment.

We explore how the correlation between host star metallicity and giant planets shapes hot Jupiter occurrence as a function of Galactic birth radius (\rbirth) and phase-space density in the Milky Way disk. Using the GALAH and APOGEE surveys and a galaxy from the NIHAO simulation suite, we inject hot Jupiters around stars based on metallicity power laws, reflecting the trend that giant planets preferentially form around metal-rich stars. For \rbirth\ $\geq 5$ kpc, hot Jupiter occurrence decreases with \rbirth\ by $\sim -0.1%$ per kpc; this is driven by the Galaxy's chemical evolution, where the inner regions of the disk are more metal-rich. Differences in GALAH occurrence rates versus APOGEE's and the simulation's at \rbirth\ $< 5$ kpc arise from survey selection effects. APOGEE and the NIHAO simulation have more high-$\alpha$ sequence stars than GALAH, resulting in average differences in metallicity (0.2--0.4 dex), $\alpha$-process element enrichment (0.2 dex), and vertical velocities (7--14 km/s) at each \rbirth\ bin. Additionally, we replicate the result of \cite{Winter20}, which showed that over 92% of hot Jupiters are associated with stars in phase-space overdensities, or "clustered environments." However, our findings suggest that this clustering effect is primarily driven by chemical and kinematic differences between low and high-$\alpha$ sequence star properties. Our results support stellar characteristics, particularly metallicity, being the primary drivers of hot Jupiter formation, which serves as the "null hypothesis" for interpreting planet demographics. This underscores the need to disentangle planetary and stellar properties from Galactic-scale effects in future planet demographics studies.

We present a catalog of $\sim 10,000$ resolved triple star systems within 500 pc of the Sun, constructed using Gaia data. The triples include main-sequence, red giant, and white dwarf components spanning separations of 10 to 50,000 au. A well-characterized selection function allows us to constrain intrinsic demographics of the triple star population. We find that (a) all systems are compatible with being hierarchical and dynamically stable; (b) mutual orbital inclinations are isotropic for wide triples but show modest alignment as the systems become more compact; (c) primary masses follow a Kroupa initial mass function weighted by the triple fraction; (d) inner binary orbital periods, eccentricities, and mass ratios mirror those of isolated binaries, including a pronounced twin excess (mass ratios greater than 0.95) out to separations of 1000+ au, suggesting a common formation pathway; (e) tertiary mass ratios follow a power-law distribution with slope -1.4; (f) tertiary orbits are consistent with a log-normal period distribution and thermal eccentricities, subject to dynamical stability. Informed by these observations, we develop a publicly available prescription for generating mock triple star populations. Finally, we estimate the catalog's completeness and infer the intrinsic triple fraction, which rises steadily with primary mass: from $5\%$ at $\lesssim 0.5\,{\rm M_\odot}$ to $35\%$ at $2\,{\rm M_\odot}$. The public catalog provides a robust testbed for models of triple star formation and evolution.

We use global MHD galaxy simulations to investigate the effects of spiral arms on the evolution of magnetic fields and star formation within a self-regulated interstellar medium (ISM). The same galaxy is simulated twice: once with self-consistent stellar spiral arms and once more with the stellar spirals suppressed via a novel numerical approach, using the Ramses AMR code. Spiral arms continually promote star formation, with 2.6 times higher rates in the spiral galaxy. The higher rate is due to high gas columns gathered along the spiral arms, rather than increasing the star formation efficiency at a given gas column. In both cases, the magnetic field is initially amplified via a small-scale dynamo driven by turbulence due to supernova feedback. Only the spiral galaxy exhibits late-time, consistent field growth due to a large-scale dynamo (e-folding time $\sim600$ Myr). This results in volume-averaged field strengths of $\sim 1$ $\mu$G after 1 Gyr of evolution. The mean-fields tend to align themselves with the spiral arms and are coherent up to 10 kpc scales. We demonstrate a novel large-scale dynamo mechanism, whereby spiral-driven radial flows enable the mean-field amplification.

Recent cosmological data and astrophysical observations, such as the Hubble tension and the increasing preference from galaxy surveys for dynamical dark energy, have begun to challenge the standard $\Lambda$-cold dark matter cosmological model. Primordial magnetic fields (PMFs) offer a mechanism to alleviate these tensions within the framework of the standard model. These fields source excess small-scale baryon clumping, which can speed up recombination and shrink the comoving sound horizon at the surface of last scattering. Computing the modified recombination history requires coupling the radiative transport of Lyman-$\alpha$ photons to compressible magnetohydronamic simulations. Since doing so is generically computationally intractable, we have developed a linearized treatment which self-consistently computes the modified recombination history in the presence of PMF induced baryon clumping for fields with red-tilted spectra. The clumping factors we find are too small to alleviate outstanding cosmological tensions, but our general framework can be applied to other PMF spectra, and provides a significant theoretical step towards a complete account of recombination in the presence of small-scale baryon clumping.

Cepheus A HW2 is a well-studied massive young stellar object (MYSO) featuring Class II methanol masers. Recently, certain maser components have been found to flare every five years. This period went undetected in (NEO)WISE photometry due to detector saturation. However, difference imaging revealed a light echo (LE), representing a variability record. Its periodicity is similar to that of the masers. Thereby, the mid-IR light curve since 2007 could be reconstructed. Phase and period maps provide information on the circumstellar dust distribution and viewing geometry. This is the first time an LE has been used for reverberation mapping of a YSO environment and expanding its light curve into the past.

Results from pulsar timing arrays (PTAs) show evidence of a gravitational wave background (GWB) consistent with a population of unresolved supermassive black hole binaries (BHBs). The observed spectrum shows a flattening at lower frequencies that can be explained by a population of eccentric BHBs. This study aims to determine the dynamical evolution and merger timescales of the most massive BHBs, which are potential sources of the GWB. We select successive galactic major mergers from the IllustrisTNG100-1 cosmological simulation and re-simulate them at high resolution with the N-body code Griffin, down to binary separations of the order of a parsec. Coalescence timescales are estimated using a semi-analytical model that incorporates gravitational wave emission and stellar hardening. Throughout our investigation, we consider the impact of prior mergers on the remnant galaxy in the form of core scouring and anisotropy, which can influence the subsequent formation and evolution of BHBs. We find that all the binaries in our sample enter the PTA band with an eccentricity e>0.85: such a large eccentricity can impact the shape of the PTA observed GWB spectrum, and it highlights the importance of including the eccentricity of binaries when interpreting the PTA signal. Furthermore, we find that: (i) starting from initial separations of a few tens of kpcs, the dynamical friction phase lasts for a few hundred Myrs; (ii) the binary formation time is not resolution dependent; (iii) the scatter on the eccentricity at binary formation decreases with increasing resolution; (iv) triple systems form whenever a third galaxy interacts with a binary which hasn't yet reached coalescence.

We perform a systematic study of the preconditioning of the interplanetary (IP) medium due to isolated interplanetary coronal mass ejections (ICMEs). Preconditioning is highly relevant when ICMEs, ejected in close succession and direction, modify the IP medium to allow subsequent ICMEs to propagate more freely, decelerate less, and retain higher kinetic energy at larger distances. We base our study on a sample of carefully selected events. The IP medium is analyzed during time intervals of 48 hours before and after the ICMEs in order to statistically quantify their impact on the properties of the solar wind (SW) and interplanetary magnetic field (IMF). We find that the SW behind ICMEs on average exhibits reduced density (-41%) and dynamic pressure (-29%), and increased total velocity (+15%), while the trailing IMF is more intense (+14%) and more radially aligned (13 degrees). The results suggest that even relatively low speed ICMEs can significantly precondition the IP medium. The results are relevant for better understanding of CME propagation and SW interaction, and hold implications for heliospheric models and applied research of space weather.

Supermassive black hole (SMBH) binary systems are unavoidable outcomes of galaxy mergers. Their dynamics encode information about their formation and growth, the composition of their host galactic nuclei, the evolution of galaxies, and the nature of gravity. Many SMBH binaries with separations pc-kpc have been found, but closer (sub-parsec) binaries remain to be confirmed. Identifying these systems may elucidate how binaries evolve past the ``final parsec'' until gravitational radiation drives them to coalescence. Here we show that SMBH binaries in non-active galactic nuclei can be identified and characterized by the gravitational lensing of individual bright stars, located behind them in the host galaxy. The rotation of `caustics' -- regions where sources are hugely magnified due to the SMBH binary's orbit and inspiral -- leads to Quasi-Periodic Lensing of Starlight (QPLS). The extreme lensing magnification of individual bright stars produces a significant variation in the host galaxies' luminosity; their lightcurve traces the orbit of the SMBH binary and its evolution. QPLS probes the population of sources observable by pulsar timing arrays and space detectors (LISA, TianQin), offering advance warning triggers for merging SMBHs for coincident or follow-up GW detections. SMBH population models predict $1-50\; [190-5,000] \left({n_\star}/{\rm pc}^{-3}\right)$ QPLS binaries with period less than $10\; [40]$ yr with comparable masses and $z<0.3$, where $n_\star$ is the stellar number density. Additionally, stellar and orbital motion will lead to frequent instances of single/double flares caused by SMBHBs with longer periods. This novel signature can be searched for in a wealth of existing and upcoming time-domain photometric data: identifying quasi-periodic variability in galactic lightcurves will reveal an ensemble of binary systems and illuminate outstanding questions around them.

The horizontal branch (HB) phase of stellar evolution plays a critical role in understanding the life cycle of stars, particularly in the context of low-mass stars. However, generating theoretical evolutionary tracks for HB stars is computationally intensive, especially when attempting to match a wide range of observational data. This work presents an extensive grid of HB evolutionary tracks computed using the PGPUC code, covering a broad range of chemical compositions, progenitor masses, and alpha-element distributions. The aim is to provide a robust tool for interpreting HB stellar populations and advancing our understanding of their diverse properties. The evolutionary tracks are made publicly available through the PGPUC Online database for easy access and interpolation. We computed over 19,000 HB evolutionary tracks encompassing a wide range in terms of mass, metallicity, helium abundance, and alpha-element enhancement, along with different progenitor masses. Using these tracks, we calculated zero-age, middle-age, and terminal-age loci. The PGPUC database now includes a fine grid of HB evolutionary tracks, allowing for precise interpolation. Key findings include the dependence of HB morphology on progenitor mass, metallicity, and helium abundance.

We present high-resolution ($\sim$0.05"; 8 au) dust continuum and molecular line observations toward the Class I protostellar system IRAS 04169+2702 in the Taurus B213 region, as part of the ALMA Large Program Early Planet Formation in Embedded Disks (eDisk). The 1.3-mm dust continuum emission traces a circumstellar disk with a central depression toward the protostar. Our VLA observations of the same target reveal a single central peak dominated by the free-free emission, which coincides with the depression of the thermal dust emission. The mean spectral index of the thermal dust emission from 1.3 mm to 1.4 cm is approximately 2.8, suggestive of the presence of grains grown to millimeter or centimeter sizes in the disk. Velocity gradients along the disk major axis are seen in emission from $^{12}$CO (2-1), $^{13}$CO (2-1), and C$^{18}$O (2-1) molecular lines. The position-velocity diagrams of these lines unveil a Keplerian-rotating disk with a radius of $\sim$21 au around a 1.3 $M_{\odot}$ protostar, as well as an infalling and rotating envelope with the angular momentum conserved. In addition to the compact disk, large-scale infalling spiral structures extending up to approximately 1400 au, streamers, are discovered in C$^{18}$O (2-1), SO (6$_5$-5$_4$), and H$_2$CO (3$_{0, 3}$-2$_{0, 2}$) as well as in the 1.3-mm continuum emission. Notably, in the region closer to the protostar, the spatial coincidence of C$^{18}$O and SO may indicate the presence of a shock related to accretion through the spiral arms.

The inner circumgalactic medium (CGM) of galaxies, where disk and halo processes intersect, remains poorly characterized despite its critical role in regulating galaxy evolution. We present results from Project AMIGA Insider, mapping Andromeda's (M31) inner CGM within 0.25 R_vir (~75 kpc) using 11 QSO sightlines, bringing our total sample to 54 sightlines from the disk to 2 R_vir. We detect a clear transition between M31's thick disk and CGM at R < 30 kpc, where low/intermediate ions show thick-disk corotating components with higher column densities than the CGM ones, while high ions exhibit similar column densities in both the CGM and thick disk. Beyond this region, all ion column densities decrease with impact parameter, with steeper gradients for low ions than high ions. The inner CGM (R < 100 kpc) shows more complex gas phases and multi-component absorption compared to the predominantly single-component outer CGM. We find no significant azimuthal dependence for any observed ions, suggesting M31's CGM is shaped by radial processes (e.g., cooling flows, precipitation) rather than disk-aligned outflows. We estimate the total metal mass in M31's cool (SiII, SiIII, SiIV) CGM within R_vir to be (1.9+/-0.3_stat+/-0.7_sys)x10^7 M_sun, leading to a cool gas mass of approximately 6x10^9 (Z/0.3 Z_sun)^-1 M_sun. The warmer OVI gas may contain at least 10 times more metal and gas mass. Compared to the COS-Halos L* galaxies, M31's cool CGM shows lower Si column densities at R < 0.4 R_200 and lower cool CGM masses, possibly resulting from M31's higher halo mass or different environments.

We use FIRE-2 zoom cosmological simulations of Milky Way size galaxy halos to characterize the shape of the J-factor emission signal on the sky. We find that, at a fixed dJ/d{\Omega} contour,the shape is well-fit by an ellipse, with semi-major axis Rmajor and semi-minor axis Rminor measured in degrees on the sky. We use least squares fitting to fit ellipses to the J-factor emission, viewed from the solar location. The ratio of minor to major axes (Rminor/Rmajor) allows us to characterize the shape of each contour, with ratio 1.0 corresponding to circular/spherical emission. These results provide new expectations for the shape of dark matter annihilation emission signals we might expect to see if dark matter is annihilating with its own anti-particle. We find that both the shape and angular orientation of the emission signal is different than results predict from dark matter only simulations. In terms of shape, we find that the ratios of semi-minor to semi-major axis is always consistently at 0.8, revealing a consistent circular emission shape, whereas with dark matter only emission signals the range of ratios from halo to halo is broader, and often closer to 0.5, showing a much more elliptical shape in general. In terms of angular orientation, we find that the major axis of the J-factor emission signal maps for FIRE halos are consistently aligned with the galactic plane within a few degrees, meaning that excess emission out of the plane would be hard to explain with a dark matter annihilation signal. However, we also find that the expected emission signal would be consistent with Fermi-LAT measurements showing a galactic center excess in full hydrodynamical simulations with stars and gas included, while dark-matter-only simulations do not produce the expected signal.

The classification of exoplanets has been a longstanding challenge in astronomy, requiring significant computational and observational resources. Traditional methods demand substantial effort, time, and cost, highlighting the need for advanced machine learning techniques to enhance classification efficiency. In this study, we propose a methodology that transforms raw light curve data from NASA's Kepler mission into Gramian Angular Fields (GAFs) and Recurrence Plots (RPs) using the Gramian Angular Difference Field and recurrence plot techniques. These transformed images serve as inputs to the Vision Transformer (ViT) model, leveraging its ability to capture intricate temporal dependencies. We assess the performance of the model through recall, precision, and F1 score metrics, using a 5-fold cross-validation approach to obtain a robust estimate of the model's performance and reduce evaluation bias. Our comparative analysis reveals that RPs outperform GAFs, with the ViT model achieving an 89.46$\%$ recall and an 85.09$\%$ precision rate, demonstrating its significant capability in accurately identifying exoplanetary transits. Despite using under-sampling techniques to address class imbalance, dataset size reduction remains a limitation. This study underscores the importance of further research into optimizing model architectures to enhance automation, performance, and generalization of the model.

In this work, we have developed a 6-dimensional joint probability density function for the 3-dimensional position and 3-dimensional velocity vectors of space objects in the Low Earth Orbit (LEO) based on the Principle of Maximum Entropy (MaxEnt), adhering to the principle of energy conservation. For the problem under consideration, maximizing entropy subject to energy conservation ensures that the derived probability density function (PDF) is the best representation of the uncertainty of a space object while the sampled position and velocity vectors from the PDF adhere to the orbital dynamics. We approach the entropy maximization by constructing a Lagrangian functional incorporating the energy conservation constraint and the normalization constraint of the PDF using Lagrange multipliers, setting the functional derivative of the Lagrangian to zero. This PDF can be used to generate position and velocity samples for space objects without any prior assumption and can further be utilized for orbital uncertainty propagation either using the Monte-Carlo method or by direct propagation of the PDF through the Fokker-Planck Equation.

Galaxy model fitting is widely employed to estimate properties such as galaxy shape, size, and color. Understanding how the outputs of galaxy model fitting respond to weak-lensing shear distortions is crucial for accurate shear estimation and mitigating shear-related systematics in weak lensing image analyses. In this paper, we investigate how the fitted parameters - specifically flux, size, and shape - respond to weak-lensing shear distortions within the AnaCal framework. To achieve this, we introduce quintuple numbers, a novel algebraic system inspired by dual numbers from automatic differentiation. Quintuple numbers enable the propagation of shear response information throughout the entire model-fitting process by linking analytical pixel shear responses to those of the fitted parameters. We integrate quintuple numbers into the AnaCal framework to derive the shear responses of shapes estimated with model fitting and validate the pipeline using image simulations that include realistic blending. Our results demonstrate that the multiplicative bias remains below 0.003 for ground-based, oversampled images.

We present high-resolution Keck Cosmic Web Imager (KCWI) and MUSE IFU spectroscopy of VV 114, a local infrared-luminous merger undergoing a vigorous starburst and showing evidence of galactic-scale feedback. The high-resolution data allow for spectral deblending of the optical emission lines and reveal a broad emission line component ($\sigma_{\rm{broad}} \sim$~100--300 km s$^{-1}$) with line ratios and kinematics consistent with a mixture of ionization by stars and radiative shocks. The shock fraction (percent ionization due to shocks) in the high velocity gas is anticorrelated with projected surface number density of resolved star clusters, and we find radial density profiles around clusters are well fit by models of adiabatically expanding cluster winds driven by massive stellar winds and supernovae (SNe). The total kinetic power estimated from the cluster wind models matches the wind+SNe mechanical energy deposition rate estimated from the soft band X-ray luminosity, indicating that at least 70\% of the shock luminosity in the galaxy is driven by the star clusters. \emph{Hubble Space Telescope} narrow band near-infrared imaging reveals embedded shocks in the dust-buried infrared nucleus of VV 114E. Most of the shocked gas is blueshifted with respect to the quiescent medium, and there is a close spatial correspondence between the shock map and the \emph{Chandra} soft band X-ray image, implying the presence of a galactic superwind. The energy budget of the superwind is in close agreement with the total kinetic power of the cluster winds, confirming the superwind is driven by the starburst.

Type-I X-ray bursts observed from neutron stars originate from intermittent unstable thermonuclear burning of accreted matter on these stars. Such bursts, particularly those reaching the Eddington luminosity and having a temporary photospheric radius-expansion due to radiation pressure, provide a testbed to study nuclear fusion processes in intense radiation, gravity, and magnetic fields. Here, we investigate time-resolved spectroscopic properties of a type-I burst from the accretion-powered millisecond X-ray pulsar IGR J17591-2342. Our basic spectral model includes an absorbed blackbody to describe the burst emission and an absorbed power law to represent the non-burst emission. The blackbody normalisation shows two consecutive humps aligned with blackbody temperature dips during the burst. Such an unusual behaviour could imply two consecutive photospheric radius-expansion events during the same burst or a systematic metallicity evolution in the neutron star atmosphere. However, our spectral analysis suggests the latter option is less likely to be happening for IGR J17591-2342. The novel former option implies that sufficient fuel survived after the first photospheric radius-expansion event to power a second similar event a few seconds later, challenging the current theoretical understanding. If confirmed, the double photospheric radius-expansion event observed in IGR J17591-2342 suggests the possibility of avoiding photospheric expansion at luminosities exceeding Eddington. Mechanisms such as temporary enhancement of the magnetic field by convection and confinement of the plasma could be invoked to explain the peculiar behaviour of the source.

Context. Intracluster Light (ICL) is a faint stellar component of galaxy groups and clusters bound to the cluster potential, and making up a significant fraction of the cluster mass. ICL formation and evolution is strongly linked to the Brightest Cluster Galaxies of clusters. Aims. To compare the properties and progenitor galaxies of the Intracluster Light (ICL) and Brightest Cluster Galaxies (BCGs) of clusters and groups at redshift z = 0, and determine how they coevolve. Methods. We select 127 clusters and groups in the hydrodynamic Illustris-TNG100 simulation above a mass of $10^{13} M_\odot$. We divide the ICL from the BCG by applying a surface brightness cut at the Holmberg radius of 26.5 mag/arcsec$^{-2}$, where star particles within this radius are defined as being attached to the BCG, and outside, the ICL. We then study the properties and formation history of the ICL and BCG. Results. We find the ICL is generally composed of material from stripped or merged intermediate mass galaxies, with a smaller in-situ component, while the BCG is composed of more massive merged galaxies and has a larger in-situ fraction. The ICL mass fraction increases weakly with cluster mass, declines with concentration and increases with time since the BCGs most recent major merger. The ICL is bluer and more metal-poor than the BCG, but there is no significant difference in the age of the material. Universally, BCG+ICL systems have negative colour and metallicity gradients. The ICL and BCG share a high fraction of progenitor galaxies, but the most significant progenitor is frequently not shared. Conclusions. ICL properties and formation are tied to the formation histories of the host cluster and BCG, and thus their properties are individual to each system. Although the ICL and BCG coevolve, they have distinct formation histories and properties.

We decomposed the G333 complex and the G331 giant molecular cloud into multi-scale hub-filament systems (HFs) using the high-resolution $^{13}$CO (3$-$2) data from LAsMA observations. We employed the filfinder algorithm to identify and characterize filaments within HFs. Compared with non-HFs, HFs have significantly higher density contrast, larger masses and lower virial ratios. Velocity gradient measurements around intensity peaks provide evidence of gas inflow within these structures. There may be an evolutionary sequence from non-HFs to HFs. Currently, non-HFs lack a distinct gravitational focusing process that would result in significant density contrast. The density contrast can effectively measure the extent of gravitational collapse and the strength of the gravitational center of the structure that definitively shape the hub-filament morphology. Combined with the kinematic evidence in our previous studies, we suggest that molecular clouds are network structures formed by the gravitational coupling of multi-scale hub-filament structures. The knots in the networks are the hubs, they are the local gravitational centers and the main star-forming sites. Actually, clumps in molecular clouds are equivalent to the hubs. The network structure of molecular clouds can naturally explain that feedback from protoclusters does not significantly change the kinematic properties of the surrounding embedded dense gas structures, as concluded in our previous studies.

We introduce a geometric embedding method that directly constrains $H_0$ by exploiting the exact kinematic relation $\dot{z} = H_0(1+z) - H(z)$ in a unified observational framework. This approach reconstructs expansion history without requiring specific cosmological models, incorporating both existing Cosmic Chronometer data and forecasted Sandage-Loeb measurements. Through mathematical embedding of observational constraints, we achieve $\sim 1.6\%$ precision in $H_0$ determination, matching current Cepheid distance ladder accuracy while maintaining model independence. The method demonstrates how fundamental kinematic constraints enable precision cosmology, providing a promising pathway to address the Hubble tension as future redshift drift observations become available.

We present results of more than 13 years of Fermi-LAT data analysis for the Vela pulsar from 60 MeV to 100 GeV and its pulsar wind nebula (PWN), Vela-X, for E > 1 GeV in the off-pulse phases. We find the Vela-X PWN can be best characterized using two extended components: a large radial Gaussian accompanied by an off-set, compact radial disk, both with a similar spectral index, \Gamma \sim 2.3. The common spectral properties support a common PWN origin, but a supernova remnant component is plausible for the compact radial disk. With an updated Vela-X model, the phase resolved spectral properties of the Vela pulsar are explored through a phase-resolved analysis. The phase-resolved spectral properties of the pulsar are presented, such as the SED peak energy E$_p$, the width of the SED at its peak, d$_p$, and the asymptotic (low-energy) spectral index, $\Gamma_0$, are presented. The best-fit spectral models for each LAT pulse peak (Peak 1 and Peak 2) are extrapolated to UV energies and compared to archival, phase-resolved spectra at UV, X-ray, soft \gamma-ray and TeV energies. We also discuss the physical implications of our modeling and the data comparisons.

We analyze the expansion signatures of 35 HII regions mapped in [CII] 158 micron emission by the Stratospheric Observatory for Infrared Astronomy (SOFIA). The [CII] emission primarily traces photodissociation regions (PDRs) at the transition between ionized and neutral gas. The brightness and narrow linewidth of [C II] allow us to measure PDR expansion. Bubble-shaped regions often exhibit expansion, while irregular-shaped ones are less likely to. Of the 35 HII regions, 12 (~34%) exhibit clear expansion in position-velocity (PV) diagrams, making them expansion candidates (ECs), with an average expansion velocity of ~12.2 km/s. The remaining 23 regions show no clear expansion signatures, though they may still be expanding below detection limits. Blueshifted expansion is more common (eight ECs solely blueshifted; one redshifted; three both), with mean velocities of ~10.9 km/s (blueshifted) and ~13.2 km/s (redshifted). A comparison of our observations to spherical expansion models supports expansion in eight of 12 ECs. Estimated dynamical ages are 10 to 100 times shorter than the ionizing star lifetimes, in agreement with the results of previous studies. Of the 35 regions, 14 (~40%) appear as [CII] bubbles; nine of the 12 ECs are bubble-shaped. Thermal pressure likely drives expansion in M43, while stellar winds dominate in M17, M42, RCW 120, and RCW 79. For other ECs, available data do not allow a definitive conclusion. Larger samples and more information about ionizing sources are needed to refine our understanding of HII region feedback and evolution.

It is shown that a simple quasi-equilibrium analysis of a multi-component plasma can be harnessed to explain catastrophic energy transformations in astrophysical objects. We limit ourselves to the particular class of binary systems for which the typical plasma consists of one classical ion component, and two relativistic electron components - the bulk degenerate electron gas with a small contamination of hot electrons. We derive, analytically, the conditions conducive to such a catastrophic change. The pathway to such sudden changes is created by the slow changes in the initial parameters so that the governing equilibrium state can no longer be sustained and the system must find a new equilibrium that could have vastly different energy mix - of thermal, flow-kinetic and magnetic energies. In one such scenario, macro-scale flow kinetic, and magnetic energies abound in the final state. For the given multi-component plasma, we show that the flow (strongly Super-Alfvénic) kinetic energy is mostly carried by the small hot electron component. Under specific conditions, it is possible to generate strong macro-scale magnetic (velocity) field when all of the flow (magnetic) field energy is converted to the magnetic (velocity) field energy at the catastrophe. The analysis is applied to explain various observed characteristics of white dwarf (WD) systems, in particular, of the magnetic and dense/degenerate type.

The Indian Pulsar Timing Array (InPTA) employs unique features of the upgraded Giant Metrewave Radio Telescope (uGMRT) to monitor dozens of the International Pulsar Timing Array (IPTA) millisecond pulsars (MSPs), simultaneously in the 300-500 MHz and the 1260-1460 MHz bands. This dual-band approach ensures that any frequency-dependent delays are accurately characterized, significantly improving the timing precision for pulsar observations, which is crucial for pulsar timing arrays. We present details of InPTA's second data release that involves 7 yrs of data on 27 IPTA MSPs. This includes sub-banded Times of Arrival (ToAs), Dispersion Measures (DM), and initial timing ephemerides for our MSPs. A part of this dataset, originally released in InPTA's first data release, is being incorporated into IPTA's third data release which is expected to detect and characterize nanohertz gravitational waves in the coming years. The entire dataset is reprocessed in this second data release providing some of the highest precision DM estimates so far and interesting solar wind related DM variations in some pulsars. This is likely to characterize the noise introduced by the dynamic inter-stellar ionised medium much better than the previous release thereby increasing sensitivity to any future gravitational wave search.

The analysis of biorelevant molecules in returned mission samples such as from the carbonaceous asteroid (162173) Ryugu is key to unravelling the role of extraterrestrial organics in the evolution of life. Coordinated analyses using minimally destructive techniques at the finest length-scales on pristine samples are particularly important. Here, we identify the chemical signature of unique globular and nitrogen-containing diffuse organic matter in asteroid Ryugu and map the distribution of these biorelevant molecules with unprecedented detail. Using a novel electron-microscopy-based combination of vibrational and core-level spectroscopy, we disentangle the chemistry and nanoscale petrography of these organics. We show that some of these organics contain soluble and highly aliphatic components as well as NHx functional groups, that must have formed in outer solar nebula environments before parent body incorporation. These novel coordinated analyses will open up new avenues of research on these types of precious and rare asteroidal dust samples.

Solar spectropolarimetric inversion -- inferring atmospheric conditions from the Stokes vector -- is a key diagnostic tool for understanding solar magnetism, but traditional inversion methods are computationally expensive and sensitive to local minima. Advances in artificial intelligence (AI) offer faster solutions, but are often restricted to shallow models or a few spectral lines. We present a proof-of-concept study using a transformer machine learning (ML) model for multi-line, full-Stokes inversion, to infer stratified parameters from synthetic spectra produced from 3D magnetohydrodynamic simulations. We synthesise a large set of Stokes vectors using forward modelling across 15 spectral lines spanning the deep photosphere towards the chromosphere. The model maps full-Stokes input to temperature, magnetic field strength, inclination, azimuth (encoded as $\sin2\phi$, $\cos2\phi$), and line-of-sight velocity as a function of optical depth. The transformer incorporates an attention mechanism that allows the model to focus on the most informative regions of the spectrum for each inferred parameter, and uses positional embedding to encode wavelength and depth order. We benchmark it against a multilayer perceptron (MLP), test robustness to noise, and assess generalisation. The transformer outperforms the MLP, especially in the higher layers and for magnetic parameters, yielding higher correlations and more regularised stratifications. The model retains strong performance across a range of noise levels typical for real observations, with magnetic parameter inference degrading predictably while temperature and velocity remain stable. We establish transformer architectures as a powerful tool for spectropolarimetric regression. This approach paves the way for analysis of large datasets from large solar telescopes.

During Type III solar radio bursts, electromagnetic waves are radiated at plasma frequency $\omega_p$ and its harmonics by electrostatic wave turbulence generated by electron beams ejected by Sun in randomly inhomogeneous solar wind and coronal plasmas. These emissions, detected since decades by spacecraft and radiotelescopes, are split by the plasma magnetic field into three modes $\mathcal{X}$, $\mathcal{O}$ and $\mathcal{Z}$ of different dispersion, polarization and radiation properties. This work demonstrates, using three independent and converging approaches, that only a small fraction of electromagnetic energy radiated at $\omega_p$ ($\lesssim10\%$) is escaping from beam-generated radio sources, mainly as $\mathcal{O}$-mode waves and, depending on plasma conditions, as $\mathcal{X}$-mode waves. Most energy is radiated in $\mathcal{Z}$-mode and can therefore be only observed close to sources. Results have major implications for solar radio emission and provide strong support for interpretation of observations performed up to close distances to Sun by Parker Solar Probe and Solar Orbiter spacecraft.

The masses of the white dwarfs in a binary carry information about previous mass-transfer phases. The core mass -- radius relation of low-mass giants gives the size of the progenitor of a helium white dwarf at the moment it last filled its Roche lobe. Previously, we used this information for a few observed systems to propose a new mass-transfer type, based on an angular momentum balance. Our aim is to investigate if stable mass transfer instead of the angular momentum prescription is consistent with the observed double helium white dwarf masses. We reconstruct the progenitor evolution of observed double helium white dwarfs using the core mass -- radius relation and evaluate if the periods at the start of the second phases of mass transfer are consistent with the outcome of stable mass transfer. More generally, we calculate the mass distribution of double helium white dwarfs for three different progenitors scenarios: double common envelope (with parameter $\alpha \lambda$), angular momentum prescription (with parameter $\gamma$) and stable mass transfer. We find that the observed systems are generally not consistent with stable mass transfer. Stable mass transfer leads to a tight correlation between the two white dwarf masses in a binary that is not consistent with the observed mass distribution. Double common envelope evolution is a particularly poor fit to the observations. The angular momentum prescription can populate the observed mass distribution, but not perfectly. We conclude that the first phase of mass transfer initiated on the red giant branch in low-mass systems does not generally proceed as stable mass transfer nor as common envelope, and thus is poorly understood. This may be related to the fact that for many observed binaries that have finished the first phase of mass transfer the orbit is eccentric, which is an unexpected outcome of mass transfer.

Modelling star-forming galaxies is crucial for upcoming observations of large-scale matter and galaxy distributions with galaxy redshift surveys and line intensity mapping (LIM). We introduce CosmoGLINT (Cosmological Generative model for Line INtensity mapping with Transformer), a Transformer-based generative framework designed to create realistic galaxy populations from dark matter (DM)-only simulations. CosmoGLINT auto-regressively generates sequences of galaxy properties -- including star formation rate (SFR), distance to the halo centre, and radial and tangential velocities relative to the halo -- conditioned on halo mass. Trained on the IllustrisTNG hydrodynamic simulation, the model reproduces key statistical properties of the original data, including the voxel intensity distribution and the power spectrum both in real and redshift space. It can efficiently generate a number of different realisations of the designated galaxy populations, enabling the creation of mock LIM/redshift survey catalogues from large halo catalogues produced by fast DM-only simulations. We show that our model trained at multiple redshifts can be applied to DM halo lightcone data to generate a realistic mock galaxy lightcone with incorporating the redshift evolution of the galaxy population. The mock catalogues can be readily used to derive statistical quantities and to develop data analysis pipelines for ongoing and future wide-field surveys.

We aim to obtain the equations which govern the evolution of the mass density profile of a dense irrotational molecular cloud (MC). Our study is based on the notion of "ensemble of MCs", introduced in our previous work. The MCs are modeled by use of the "ensemble abstract representative member": a spherically symmetric, isotropic and isothermal MC which accretes radially matter from its surroundings. Applying the equations of hydrodynamics to a self-gravitating isothermal spherical gas cloud, we obtain a system of 2 first-order non-linear partial differential equations, which govern the evolution of two unknown fields: the exponent of density profile and the accretion velocity. Assuming a steady-state flow, we get approximate solutions using the method of leading-order terms. Far from the cloud centre the obtained density profile is $\varrho=\ell^{-2}$ and the accretion velocity is constant, while near to the centre $\varrho=\ell^{-3/2}$ and $v_{\rm a}\propto\ell^{-1/2}$. Through our dynamical equations, the obtained solutions coincide completely with those found using the equation of energy conservation of a fluid element (in our previous work). Also, combining the equations of energy balance for a fluid element, we arrive at the conclusion that the cloud layers far from the centre are in a stable dynamical state if the accretion velocity flow is sub- or transsonic; otherwise they are marginally stable (moderately supersonic flow) or unstable (supersonic flow). Both solutions are consistent only if the accretion flow is subsonic and hence the outer layers are stable. Finally, under the assumption that both the accretion velocity and density scale with $\ell$ and their power-law exponents are position-independent, we show that the density scaling exponent far away from the centre is $p=2$ and this value is an attractor. Hence this value should be observable in dense MCs.

The Magellanic Clouds are two nearby dwarf irregular galaxies whose study can help us in understanding galaxy and stellar evolution. In particular, the Large Magellanic Cloud, the larger one, contains approximately 30 billion stars at various evolutionary stages. In this work, we present an SDSS $gri$-bands photometric analysis based on multiple images acquired with DECam, the Dark Energy Camera, installed on the Blanco telescope at Cerro Tololo Inter-American Observatory (Chile). We performed a full image analysis and photometric calibration, resulting in a photometric catalog named COSMIC-L, consisting of 57,997,665 stars, of which 18,676,294 contain estimates for all three $gri$ magnitudes, resulting in a completeness magnitude of $\simeq 21$ and a limiting magnitude of $\simeq 22$ in all three bands.

Context. The classification of active galactic nuclei (AGNs) is a challenge in astrophysics. Variability features extracted from light curves offer a promising avenue for distinguishing AGNs and their subclasses. This approach would be very valuable in sight of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). Aims. Our goal is to utilize self-organizing maps (SOMs) to classify AGNs based on variability features and investigate how the use of different subsets of features impacts the purity and completeness of the resulting classifications. Methods. We derived a set of variability features from light curves, similar to those employed in previous studies, and applied SOMs to explore the distribution of AGNs subclasses. We conducted a comparative analysis of the classifications obtained with different subsets of features, focusing on the ability to identify different AGNs types. Results. Our analysis demonstrates that using SOMs with variability features yields a relatively pure AGNs sample, though completeness remains a challenge. In particular, Type 2 AGNs are the hardest to identify, as can be expected. These results represent a promising step toward the development of tools that may support AGNs selection in future large-scale surveys such as LSST.

Recent discoveries of Ultra Faint dwarf galaxies (UFG's) infalling onto the Milky Way, namely Leo K \& M at $r \simeq 450$kpc, considerably strengthens the case that UFG's constitute a distinct galaxy class that is inherently smaller, fainter and metal poor compared to the classical dwarf spheroidals (dSph). This distinction is at odds with the inherent continuity of galaxy halo masses formed under scale-free gravity for any standard dark matter model. Here we show that distinct galaxy classes do evolve in cosmological simulations of multiple light bosons representing the 'Axiverse' proposal of String Theory, where a discrete mass spectrum of axions is generically predicted to span many decades in mass. In this context, the observed UFG class we show corresponds to a relatively heavy boson of $3\times 10^{-21}$eV, including Leo K \& M, whereas a lighter axion of $10^{-22}$eV comprises the bulk of dark matter in all larger galaxies including the dSph's. Although Leo M is larger in size than Leo K, we predict its velocity dispersion to be smaller $\simeq 1.7$km/s, compared to $\simeq 4.5$km/s for Leo K, since soliton cores are required by the Uncertainty Principle to be wider at lower momentum. This scenario can be definitively tested using millisecond pulsars close to the Galactic center, where the Compton frequencies of the heavy and light bosons imprint monotone timing residuals that may be detected by SKA on timescales of approximately 1 week and 4 months, respectively.

The inner regions of protoplanetary discs, which encompass the putative habitable zone, are dynamically complex, featuring a well ionised, turbulent active inner region and a poorly ionised dead outer region. In this first paper, we investigate a base level model of the magnetohydrodynamic processes around this dead-active zone interface, using five three-dimensional global magnetohydrodynamic simulations in the zero-net flux regime. We employ physically motivated profiles for Ohmic resistivity and ambipolar diffusion, alongside a simplified thermodynamic model comprising a hot disc and cool corona. Our results show that, first, large-scale coherent poloidal magnetic field loops form in the magnetorotational instability active region. These loops lead to the accumulation of tightly wound magnetic flux at the disc-corona temperature transition, driving strong, localised accretion flows in the surface layers of the active region. Second, an axisymmetric pressure maximum, that extends across multiple disc scale heights, develops as a result of outward mass transport from the active region. This triggers the Rossby wave instability and leads to the development of anticyclonic vortices, though their coherence is weakened due to the turbulent dynamics at the interface. Third, the dead zone develops magnetic field with a distinct morphology, likely resulting from the outward diffusion of large-scale poloidal loops in the active zone. This self-consistently generated field exhibits a vertical structure consistent with that required to drive accretion in the inner dead zone under the weak magnetic wind paradigm. In the second paper in the series, we extend this work to the vertical-net flux regime, where global magnetic flux transport and magnetically driven outflows become dynamically significant.

Interferometric observations of protoplanetary disks by VLTI and ALMA have greatly improved our understanding of the detailed structure of these planetary birthplaces. These observations have revealed a variety of large-scale disk substructures, including rings, gaps, and spirals, spanning tens to hundreds of au, supporting the predictions of planet formation models. Recent instruments, such as MATISSE at the VLTI, allow one to resolve and investigate the inner few au of protoplanetary disks in nearby star formation regions, shedding light on the traces of planet formation and evolution at these small scales. The aim of this work is to assess the feasibility of interferometric observations of small-scale planet-induced substructures in protoplanetary disks in nearby star-forming regions. We aim to characterize these substructures in multi-wavelength and multi-epoch observations and subsequently differentiate between simulation parameters. On the basis of 3D hydrodynamic simulations of embedded planetary companions and subsequent 3D Monte Carlo radiative transfer simulations, we calculated and analyzed interferometric observables, assuming observations with the VLTI in the K, L, M, and N bands. The hydrodynamic simulations exhibit mass-dependent planet-induced density waves that create observable substructures, most notably for the considered case of a 300 $M_{\oplus}$ planet. These substructures share similarities with observed large-scale structures and feature a prominent accretion region around the embedded planet. The visibilities show a detectable variability for multi-epoch VLTI/GRAVITY and VLTI/MATISSE observations, caused by the orbital motion of the planet, that are distinguishable from other sources of variability due to their unique combination of timescale and amplitude.

The evolution of magnetic braking and dynamo processes in subgiant stars is essential for understanding how these stars lose angular momentum. We investigate the magnetic braking and dynamo evolution of $\beta$ Hydri, a G-type subgiant, to test the hypothesis of weakened magnetic braking and the potential rejuvenation of large-scale magnetic fields. We analyze spectropolarimetric observations from HARPS (HARPSpol; polarimetric mode of High Accuracy Radial velocity Planet Searcher), and combine them with archival X-ray data and asteroseismic properties from TESS (Transiting Exoplanet Survey Satellite) to estimate the current wind braking torque of $\beta$ Hydri. Despite experiencing weakened magnetic braking during the second half of its main-sequence lifetime, our results indicate that $\beta$ Hydri has regained significant magnetic activity and a large-scale magnetic field. This observation aligns with the "born-again" dynamo hypothesis. Furthermore, our estimated wind braking torque is considerably stronger than what would be expected for a star in the weakened magnetic braking regime. This suggests that subgiants with extended convective zones can temporarily re-establish large-scale dynamo action. These results provide critical constraints on stellar rotation models and improve our understanding of the interplay between magnetic field structure, stellar activity cycles, and angular momentum evolution in old solar-type stars.

We present the first IXPE spectro-polarimetric observation of the black hole candidate MAXI J1744$-$294, a transient X-ray source discovered during a bright 2025 outburst in the Galactic center region. During the $\sim$150 ks observation, the source was found in the soft state, and its spectrum was well described by an absorbed multicolor disk with a minor high-energy tail. No significant polarization was detected, and we derived a 3$\sigma$ upper limit on the polarization degree of $1.3\%$ in the 2--8 keV energy band. This result is consistent with previous findings for soft-state black hole binaries observed at low to intermediate inclination angles. By comparing the polarization degree upper limit with theoretical predictions for standard accretion disk emission, we constrain the disk inclination to $i \lesssim 38^\circ$--$71^\circ$, depending on the black hole spin and the disk atmosphere albedo.

Around our Sun, terrestrial planets did not grow beyond Earth in mass, while super-Earths are found to orbit approximately every other solar-like star. It remains unclear what divides these super-Earth systems from those that form terrestrial planets, and what role wide-orbit gas giants play in this process. Here, we show that the key uncertainty is the degree of viscous heating in the inner disc, which regulates the pebble accretion efficiency. In this parameter study, we assume pebble sizes limited by fragmentation and radial drift. The initial seed planetesimals for embryo growth are taken from the top of the streaming instability mass distribution. We then evaluate the important role of the pebble scale height and the assumed pebble fragmentation velocity. In systems with maximally efficient viscous heating, where all the accretion heating is deposited in the disc midplane, pebble accretion in the terrestrial region is suppressed. More realistic levels of viscous heating, at higher elevations, allow terrestrial embryo formation at Earth-like orbits. We also find that the role of the water iceline is minor, unless it is paired with extreme volatile loss and a change in the pebble fragmentation velocity. Furthermore, we show that in systems with gas-giant formation, the role of mutual pebble filtering by outer pebble-accreting embryos is limited, unless some mechanism of delaying inner disc growth, such as viscous heating or the presence of an iceline, is simultaneously employed. This latter point appears to be consistent with the fact that no strong suppression is seen in the occurrence rate of super-Earths in systems with known gas giants in wider orbits. We conclude that the diversity in inner-disc systems may largely be driven by complex, and as of yet poorly understood, disc accretion physics inside the water iceline.

Detecting orbital eccentricity in a stellar-mass black-hole merger would point to a non-isolated formation channel. Eccentric binaries can form in dense stellar environments such as globular clusters or active galactic nuclei, or from triple stellar systems in the Galactic field. However, confidently measuring eccentricity is challenging -- short signals from high-mass eccentric mergers can mimic spin-induced precession, making the two effects hard to disentangle. This degeneracy weakens considerably for longer-duration signals. Here, GW200208_222617 provides a rare opportunity. Originating from a relatively low-mass binary with source-frame chirp mass $\sim20$ M$_\odot$, its gravitational-wave signal spanned $\sim14$ orbital cycles in band, with no indication of data quality issues. Previous analyses for quasi-circular binaries found no evidence for spin precession and, multiple subsequent studies found the data to favour an eccentric merger despite notable technical differences. All in all, we believe GW200208_222617 is the black-hole merger event from GWTC-3 with the least ambiguous detection of eccentricity. We present a critical discussion of properties and astrophysical interpretation of GW200208_222617 as an eccentric black-hole merger using models of field triples, globular clusters, and active galactic nuclei. We find that if GW200208_222617 was indeed eccentric, its origin is consistent with a field triple or globular cluster. Formation in the inner regions of an active galactic nucleus is disfavoured. The outer regions of such a disk remain a viable origin for GW200208_222617; we demonstrate how future detections of eccentric mergers formed in such environments could be powerful tools for constraining the disk geometry.

Accurate atomic physics reference data are a crucial requirement for analysis and interpretation of observed spectra, even more so for observations with high spectral resolution. This document provides a curated list of atomic physics references frequently used for plasma diagnostics in X-ray spectroscopy, outside of comprehensive plasma models that typically come with their own underlying atomic databases. The list includes references to physical constants, laboratory benchmarks, transition energies, position and line shapes of neutral fluorescence lines, radiative branching ratios, and commonly used notation for prominent transitions. Quick-look tables for transition energies in H-, He-, and Li-like ions and line positions and shapes for fluorescence lines in neutrals. The main focus is on K-shell transitions. For the H- and He-like tables, we cite state-of-the art calculations that we consider currently the best available reference energies, which are considered high accuracy and thus typically used for energy scale calibration in laboratory measurements. Omissions in these tables are due to the lack of availability in the chosen references, and are not a statement about the relevance of these lines. Due to their complex and highly source-dependent line shape, the atomic data for neutrals is of lower accuracy than that for the highly charged ions, and the best reference data for these line shapes typically consist of empirical models derived from very high-resolution laboratory measurements. The table for neutrals provided here is consistent with the reference used for the energy gain scale calibration of XRISM/Resolve. This document is meant to serve as a resource to help find relevant references and conveniently formatted overview tables. When making use of the information found in these papers, credit should be given to their original authors by citing the appropriate references.

A real-time broadband RFI mitigation system has been developed for the Upgraded GMRT (uGMRT). Here, we report the engineering, imaging, and beam/pulsar results from the test observations in all the bands of the uGMRT during the monsoon of 2018 and 2019. Broadband power-line RFI is almost always present in uGMRT bands 2 (120 - 240 MHz), 3 (300 - 500 MHz) and 4 (550 - 850 MHz), but the monsoon period is chosen to understand the performance in the worst RFI situations arising due to an increase in the wind-speed and moisture.

Quasi-periodic eruptions (QPEs) are repeating soft X-ray bursts from the nuclei of galaxies, tantalizingly proposed to be extreme mass ratio inspirals. Here, we report the discovery of a new galaxy showing X-ray QPEs, the fifth found through a dedicated blind search in the \emph{SRG}/eROSITA all-sky survey data, hereafter named eRO-QPE5. Its QPE duration ($t_{\rm dur}\sim0.6$\,d), recurrence time ($t_{\rm recur}\sim3.7\,$d), integrated energy per eruption ($\sim3.4 \times 10^{47}\,$erg), and black hole mass ($M_{\rm BH}\sim 2.8 \times 10^7\,M_{\odot}$) sit at the high end of the known population. Like other eROSITA or X-ray-discovered QPEs, no previous or concurrent optical-IR transient is found in archival photometric datasets, and the optical spectrum looks almost featureless. With a spectroscopic redshift of $0.1155$, eRO-QPE5 is the most distant QPE source discovered to date. Given the number of recent discoveries, we test for possible correlations and confirm a connection between $t_{\rm dur}$ and $t_{\rm recur}$, while we do not find any significant correlation involving either $M_{\rm BH}$ or the QPE temperature. The slope of the $t_{\rm dur}-t_{\rm recur}$ relation ($1.14\pm0.16$) is roughly consistent with predictions from star-disk collision models, with a preference for those that suggest that QPEs are powered by stellar debris streams around the orbiter. Considering this and previous discoveries, eROSITA has proved extremely successful in finding many QPE candidates given its grasp, namely its sensitivity and large field of view, and scanning capabilities over the full sky. We advocate the need of sensitive wide-area and time-domain oriented surveys from future-generation soft X-ray missions.

Context. AGNs are strong X-ray emitters shaped by disk-corona interactions. The soft excess (0.5-2.0 keV) reveals key information about the "warm corona" bridging the disk and hot corona. Yet, how this feature evolves with accretion properties remains poorly constrained, especially in large samples using spectral stacking. Aims. The eROSITA All-Sky Survey (eRASS:5) provides an unprecedented sample. We investigate how the average AGN X-ray spectra evolve with accretion parameters, and explore disk-corona connection by further combining stacked UV data. Methods. We developed Xstack, a novel tool that stacks rest-frame X-ray spectra and responses while preserving spectral shape through optimized weighting. We stack 17929 AGNs ("spec-z" sample, 23 Ms) with similar X-ray loudness alpha_ox, UV luminosity L_UV, and 4159 AGNs ("BH-mass" sample, 3 Ms) with similar Eddington ratio lambda_Edd and black hole mass M_BH. The resulting stacked X-ray spectra are analyzed with a phenomenological model. We further fit the stacked optical-UV-Xray SED with AGNSED model. Results. Soft excess strengthens strongly with alpha_ox and lambda_Edd (~5), while the hard X-ray spectral shape remains largely unchanged, supporting that soft excess is dominated by warm corona rather than reflection. AGNSED modeling reveals that warm corona radius (R_g units) generally increases with lambda_Edd and decreases with M_BH, or equivalently the disk-to-warm-corona transition consistently occurs near 1e4 K. The hot corona contracts with lambda_Edd and is unaffected by M_BH, aligning with disk evaporation predictions. Conclusions. The soft excess likely originates from a warm corona, with the disk to warm corona transition tied to hydrogen ionization near 1e4 K - supporting earlier eFEDS-HSC stacking results (Hagen et al. 2024). This study shows the strength of spectral stacking in probing AGN disk-corona physics.

Detecting anisotropic screening of the cosmic microwave background (CMB) holds the promise of revealing the distribution of gas in the Universe, characterizing the complex processes of galaxy formation and feedback, and studying the epoch of reionization. Estimators for inhomogeneous screening, including some recently proposed small-scale (stacked) estimators, are quadratic or higher order in the CMB temperature or polarization fields and are therefore subject to contamination from CMB lensing. We review the origin of this lensing bias and show that, when stacking on unWISE galaxies, the expected lensing bias dominates the signal if left unmitigated. Hardening techniques that null the lensing bias have been proposed for standard quadratic estimators, whereas only approximate methods have been proposed for stacked estimators. We review these techniques and apply the former to stacked estimators, presenting several strategies (including the optimal strategy) to null lensing contamination when stacking on any large-scale structure (LSS) tracer.

The hot circumgalactic medium (CGM), probed by X-ray observations, plays a central role in understanding gas flows that drive a galaxy's evolution. While CGM properties have been widely studied, the influence of a galaxy's large-scale cosmic environment on the hot gas content remains less explored. We investigate how the large-scale cosmic web affects the X-ray surface brightness (XSB) profiles of galaxies in the context of cosmological simulations. We use our novel IllustrisTNG-based lightcone, spanning $0.03 \leq z \leq 0.3$, first developed in our previous work, and generate self-consistent mock X-ray observations, using intrinsic gas cell information. We apply the filament finder DisPerSE on the galaxy distributions to identify the cosmic filaments within the lightcone. We classify central galaxies into five distinct large-scale environment (LSE) categories: clusters and massive groups, cluster outskirts, filaments, filament-void transition regions, and voids/walls. We find that the X-ray surface brightness profiles (XSB) of central galaxies of dark matter halos in filaments with $M_{\rm 200m} >10^{12}\ M_\odot$ are X-ray brighter than those in voids and walls, with $20-45%$ deviations in the radial range of $\sim (0.3-0.5) R_{\rm 200m}$. We investigate the source of this enhancement and find that the filament galaxies show higher average gas densities, temperatures, and metallicities compared to voids/walls galaxies. Our results demonstrate that the impact of the large-scale cosmic environment is imprinted on the hot CGM's X-ray emission. Future theoretical work on studying the effect of assembly history, connectivity, and gas accretion on galaxies in filaments and voids would help to further our understanding of the impact of the environment on X-ray observations.

We study modifications of gravitational wave observables, such as the wave amplitude and frequency, which follow from the quantum equivalence principle, and are expressed in terms of the inertial, gravitational and rest masses of the LIGO/Virgo mirrors. We provide bounds on the violations of the quantum equivalence principle by comparing the results with the most resolved gravitational wave events observed by the LIGO/Virgo collaboration. The formalism is equally applicable to other future ground and space-based gravitational wave detectors.

We analyze DESI DR2 data with a model-independent method and find that: (a) the expansion of the universe may speed up with a confidence level more than 2.3 $\sigma$ at redshift $z_{51}\in (0.51, 0.955)$; (b) the expansion of the universe may speed down with a confidence level greater than 1.7 $\sigma$ at redshift $z_{75}\in (0.955, 1.484)$; (c) $w_{\rm{x}}\leq w_{\rm{t}}<-1$ with confidence level exceeding 1.6 $\sigma$ at redshift $z_{53}\in (0.922, 0.955)$.

Smooth hybrid inflation is a hybrid inflation model which is free from topological defects and predicts the density perturbation with the spectral index of about $0.97$. We show that the prediction on the spectral index is robust regardless the power of nonrenormalizable terms and meets with the latest results reported by Atacama Cosmology Telescope. Viable scenarios for baryogenesis and dark matter are also discussed.

Systematic errors in detector calibration can bias signal analyses and potentially lead to incorrect interpretations suggesting violations of general relativity. In this study, we investigate how calibration systematics affect black hole (BH) spectroscopy, a technique that uses the quasinormal modes (QNMs) emitted during the ringdown phase of gravitational waves (GWs) to study remnant BHs formed in compact binary coalescences. We simulate a series of physically motivated, tunable calibration errors and use them to intentionally miscalibrate numerical relativity waveforms. We then apply a QNM extraction method -- the rational QNM filter -- to quantify the impact of these calibration errors. We find that current calibration standards (errors within $10\%$ in magnitude and $10^\circ$ in phase across the most sensitive frequency range of 20--2000 Hz) are adequate for BH ringdown analyses with existing observations, but insufficient for the accuracy goals of future upgraded and next-generation observatories. Specifically, we show that for events with a high ringdown signal-to-noise ratio of $\sim 120$, calibration errors must remain $\lesssim 6\%$ in magnitude and $\lesssim 4^\circ$ in phase to avoid introducing biases. While this analysis addresses a particular aspect of BH spectroscopy, the results offer broadly applicable quantitative benchmarks for calibration standards crucial to fully realizing the potential of precision GW astrophysics in the next-generation detector era.

FLARE is an open source data workflow orchestration tool designed for the FCC Analysis software and Key4HEP stack. Powered by b2luigi, FLARE automates and orchestrates the fccanalysis stages from start to finish. Furthermore, FLARE is capable of managing the Monte Carlo (MC) data workflow using generators inside the Key4HEP stack such as Whizard, MadGraph5 aMC@NLO, Pythia8 and Delphes. In this paper the FLARE v0.1.4 package will be explored along with its extensible capabilities and a feature rich work environment. Examples of FLARE will be discussed in a variety of use-cases, all of which can be found at this https URL. The open source repository of FLARE can be found at this https URL

We test an analytic cosmological solution within the framework of Weyl Integrable Spacetime using current observational data. In this model, dark energy is described by a pressureless fluid, while a scalar field arises naturally through the definition of the connection. This gravitational theory reveals a Chameleon Mechanism leading to a nonzero interaction between the scalar field and the matter sector. This model extends the standard $\Lambda $CDM cosmology by introducing one additional degree of freedom, allowing for deviations from $\Lambda$CDM dynamics. For the observational constraints we consider the Supernova data of Pantheon+ collaboration, the Cosmic Chronometers, the Baryonic Acoustic Oscillators of DESI DR2 collaboration and the gamma-ray bursts. We find that Weyl Integrable Spacetime fits the data in a better way than the $\Lambda$CDM. When all datasets are considered, the statistical comparison indicates a weak to moderate preference for Weyl Integrable Spacetime according to the Akaike Information Criterion and Jeffrey's scale for the Bayesian evidence.

Machine learning is often viewed as a black box when it comes to understanding its output, be it a decision or a score. Automatic anomaly detection is no exception to this rule, and quite often the astronomer is left to independently analyze the data in order to understand why a given event is tagged as an anomaly. We introduce here idea of anomaly signature, whose aim is to help the interpretability of anomalies by highlighting which features contributed to the decision.

If the discovered Higgs boson with $m_H$=125 GeV is interpreted as a $t\bar{t}$-boson where the $t$-quarks are bound by Higgs-exchange with binding energy 220 GeV, then the $2t$-baryons should have approximately the same mass as the Higgs-boson. As $m_H<m_t$, the life time of $2t$-baryons must be much bigger then the life time of $1t$-baryons. If in the primordial Universe the number of $2\bar{t}$-antibaryons was bigger than the number of $2t$-baryons, then the excess should be compensated by nucleons. The relatively long living heavy $2\bar{t}$-antibaryons could in primordial Universe fast evolve to antimatter black halls and disappear in the world of matter under the Schwarzschild spheres.

Observations of $\gamma$-ray from blazars suggest the presence of magnetic fields in the intergalactic medium, which may require a primordial origin. Intense enough primordial magnetic fields can arise from theories of dynamical electroweak symmetry breaking during the big bang, where supercooling is ended by a strongly first order phase transition. We consider theories involving new scalars and possibly vectors, including thermal particle dark matter candidates. Intense enough magnetic fields can arise if the reheating temperature after the phase transition is below a few TeV. The same dynamics also leaves testable primordial gravitational waves and possibly primordial black holes.

The tidal deformability is a key observable to test the nature of compact objects in a binary coalescence. Within vacuum General Relativity, the tidal Love numbers of a four-dimensional black hole are strictly zero, while they are non-zero and model-dependent for material objects, matter distributions around black holes, or in alternative theories of gravity. Here, we develop a model-agnostic framework based on the membrane paradigm, where the tidal properties of a spherically symmetric object are encoded in the response of a fictitious membrane in terms of viscosity coefficients. We show that both neutron stars and exotic compact objects (in particular, we provide an explicit example for thin-shell gravastars) are included in this framework, provided that the bulk and shear viscosity coefficients of the membrane have a nontrivial frequency dependence. We derive a general result for the electric and magnetic tidal Love numbers for non-rotating and spherically symmetric compact objects in terms of the viscosity coefficients, discussing their behavior with the compactness of the object, and identifying the conditions for the logarithmic behavior found in previous literature for various ultracompact objects. Finally, we show that the membrane shear viscosity coefficients associated with neutron star models feature novel quasi-universal relations, which depend only mildly on the neutron-star equation of state.

Quantum computers represent a new computational paradigm with steadily improving hardware capabilities. In this article, we present the first study exploring how current quantum computers can be used to classify different neutrino event types observed in neutrino telescopes. We investigate two quantum machine learning approaches, Neural Projected Quantum Kernels (NPQKs) and Quantum Convolutional Neural Networks (QCNNs), and find that both achieve classification performance comparable to classical machine learning methods across a wide energy range. By introducing a moment-of-inertia-based encoding scheme and a novel preprocessing approach, we enable efficient and scalable learning with large neutrino astronomy datasets. Tested on both simulators and the IBM Strasbourg quantum processor, the NPQK achieves a testing accuracy near 80 percent, with robust results above 1 TeV and close agreement between simulation and hardware performance. A simulated QCNN achieves approximately a 70 percent accuracy over the same energy range. These results underscore the promise of quantum machine learning for neutrino astronomy, paving the way for future advances as quantum hardware matures.

The recent discovery of a central compact object (CCO) within the supernova remnant HESS J1731-347, with mass $0.77^{+0.20}_{-0.17} \ M_\odot $ and radius $10.4^{+0.86}_{-0.78}$ km is the lightest and smallest compact object ever observed. We identify it as an ultra-light Neutron star (NS) and constrain the chiral invariant mass of nucleon $m_0$ from the observational data of NS using an extended parity doublet model with including the isovector scalar meson $a_0(980)$. We study the higher order asymmertic matter properties such as the symmetry incompressibility $K_{sym}$ and the symmetry skewness $Q_{sym}$ in the presence of $a_0$ meson. We find that $K_{sym}$ and $Q_{sym}$ is sensitive to the chiral invariant mass of nucleon $m_0$ in the presence of $a_0$ meson. We show that the equation of state in the present model satisfies all observational constraints within $2\sigma$ credible region including the HESS J1731-347 observation, as well as the constraint from $K_{sym}$ when $740 \,\text{ MeV} \lesssim m_0 \lesssim 860 \,\text{ MeV}$ for $L_0 = $ 57.7 MeV. Yet, the $1\sigma$ constraint from neutron stars appears to be not fully compatible with the constraint from $K_{sym}$ from the present model.

Significant uncertainties persist in describing the equation of state and internal structure of hyperon stars due to the limited understanding of the mechanisms underlying hyperon interactions. Constraining the interaction parameter space through a combination of the latest astronomical observations and hypernuclear physics experiments is therefore essential. In this study, we incorporate experimental constraints from $\Xi$ hypernuclear physics on top of $\Lambda$ hyperons considered in \citet{Sun2023APJ942.55}. Specifically, based on updated measurements of hyperon separation energies from $\Xi$ hypernuclear experiments, sets of $\Xi N$ effective interactions are constructed and a linear correlation between their scalar ($\sigma$) and vector ($\omega$) coupling strength ratios is proposed as a constraint derived from $\Xi$ hypernuclear physics. Together with experimental correlations and astronomical observational data, four types of analyses are performed to constrain hyperon-nucleon interactions and the properties of hyperon stars. Compared to the vector $\omega$ meson-hyperon coupling, the introduction of linear correlations in hypernuclear physics imposes a more substantial constraint on the scalar $\sigma$ meson-hyperon coupling, significantly enhancing its coupling strength and thereby ensuring the stiffness of the equation of state, highlighting the crucial role of hypernuclear studies in solving the hyperon puzzle problem. Consequently, a maximum mass of around $2M_{\odot}$ can be achieved with all five interactions considered in this study under the combined constraints from astronomical observations and nuclear physics. With more reliably estimated hyperon-nucleon contributions, the uncertainties in both the fractions and the threshold densities at which hyperons appear inside neutron stars are notably reduced, along with those in the mass-radius predictions.

The extreme conditions within the supernova core, a high-temperature and high-density environment, create an ideal laboratory for the search for new physics beyond the Standard Model. Of particular interest are low-energy supernovae, characterized by their low explosion energies, which place strong constraints on the new-physics energy transfer from the core to the mantle. We compute low-energy supernova constraints on lepton-flavor-violating axions and axion-like particles that couple to both electrons and muons. For axion mass above the muon mass, the electron-muon coalescence and the axion decay are dominant production and reabsorption processes, respectively. We find that the low-energy supernovae provide the most stringent constraints on the axions in the mass range of $\sim (110,550)$ MeV, probing the coupling constant down to $g_{ae\mu} \simeq {\cal O}(10^{-11})$.

General Relativity (GR) was created in November 1915 and since its creation and up to now this theory has undergone many tests. The first realistic cosmological models were proposed in the works of Friedman, written in the 1920s. For a long time Friedman's cosmological works were actually banned in Soviet Union due to philosophical reasons, since the models where the birth and evolution of the Universe occurs were considered ideologically unacceptable. Due to great achievements in relativity and cosmology and due to increasing interest to these branches of of science in last decades we recall a development of relativistic astrophysics and contribution of Russian researchers in these studies. Since one of the world leaders in physical cosmology A. A. Friedman passed away in September 1925, it is reasonable to outline the main achievements of physical cosmology over the past 100 years. We discuss also observational and theoretical achievements in confirmations of relativistic observational predictions for black holes, including the closest supermassive black hole in our Galactic Center. We outline an evolution of black hole shadow from the purely theoretical concept to observable quantities for supermassive black holes in Sgr A* and M87*.

In the present work, quasilocal Brown-York charges are derived that coincide in the large sphere limit with the conserved supertranslation hair and superrotation charges introduced by Hawking, Perry and Strominger in [45, 46]. Given these charges, a general scenario is outlined in which a non-rotating black hole completely evaporates after its collapse due to particle creation effects, whereby a genuine one-way traversable event horizon is never formed, but merely a two-way traversable dynamical (resp. future trapping) horizon. The formation of such a dynamical horizon has the consequence, as is demonstrated, that quasilocal energy transported by the considered charges, and thus information, can continuously escape through the black hole horizon to infinity; a mechanism which, as is argued, could possibly prevent information loss once the black hole formation and evaporation process comes to an end.

The beta decay of $^{77}$Ge and $^{77\mathrm{m}}$Ge, both produced by neutron capture on $^{76}$Ge, is a potential background for Germanium based neutrinoless double-beta decay search experiments such as GERDA or the LEGEND experiment. In this work we present a search for $^{77}$Ge decays in the full GERDA Phase II data set. A delayed coincidence method was employed to identify the decay of $^{77}$Ge via the isomeric state of $^{77}$As (9/2$^+$, 475 keV, ${T_{1/2} = 114}\,\mu $s, $^{77\mathrm{m}}$As). New digital signal processing methods were employed to select and analyze pile-up signals. No signal was observed, and an upper limit on the production rate of was set at $<0.216$ nuc/(kg$\cdot$yr) (90% CL). This corresponds to a total production rate of $^{77}$Ge and $^{77\mathrm{m}}$Ge of $<0.38$ nuc/(kg$\cdot$ yr) (90% CL), assuming equal production rates. A previous Monte Carlo study predicted a value for in-situ $^{77}$Ge and $^{77\mathrm{m}}$Ge production of (0.21$\pm$0.07) nuc/(kg$\cdot$yr), a prediction that is now further corroborated by our experimental limit. Moreover, tagging the isomeric state of $^{77\mathrm{m}}$As can be utilised to further suppress the $^{77}$Ge background. Considering the similar experimental configurations of LEGEND-1000 and GERDA, the cosmogenic background in LEGEND-1000 at LNGS is estimated to remain at a sub-dominant level.

In this paper, we study flat FLRW cosmology for a Poincaré gauge theory containing cubic invariants that is free from ghosts in arbitrary backgrounds in the axial and vector sectors of the torsion tensor. The new degrees of freedom can be related to hypermomentum but continue to be dynamical even in vacuum. These extra degrees of freedom open a more natural way in which to construct potential gravitational models that provide possible ways to modify astrophysical and cosmological physics. In this framework, we study two particular branches of the theory where preliminary routes of exploring these new variables are exposed. The first is the branch where the hypermomentum vanishes, while the second branch involves the setting where the perfect fluid and hypermomentum parts of the sources are independently conserved. In both settings, we find generically faster expanding cosmologies with similar estimates of the cosmic matter content as in the standard model of cosmology. Cubic Poincaré Gauge gravity offers an interesting theoretical basis on which to study cosmology, and indicates some preliminary positive constraints when compared with observational constraints.