Papers for Tuesday, Mar 04 2025

Gravitational time delays offer unique, independent measurements of the Hubble constant, $H_0$. Precise measurements of $H_0$ stand as one of the most pressing challenges in modern cosmology, and to do so with time delays requires precise lens models. While much work has focused on streamlining the modeling process to keep pace with the erumpent discovery of strongly-lensed systems, a critical step toward reducing uncertainty in $H_0$ comes from increasing the precision of individual lens models themselves. In this work, we demonstrate that the unprecedented imaging capabilities of JWST make this goal attainable. We present the first lens model for time-delay cosmography derived from JWST data, applied to the quadruply-imaged quasar WFI2033--4723. While the primary source of systematic uncertainty in time-delay cosmography is currently the mass-sheet degeneracy (MSD), the sensitivity of models to this MSD varies on how the point spread function (PSF) errors are mitigated. As the PSF is also the primary source of uncertainty in lens modeling, we focus on a comparison of different PSF modeling methods. Within the context of power-law models, we recover results in agreement with previous Hubble Space Telescope (HST)-based models, but with better precision of key lensing parameters through implementation of new PSF modeling techniques. Despite the record-holding precision of this system's HST modeling, we achieve an additional 22% increase in precision of the Fermat potential difference, directly reducing uncertainties of cosmological inference. These results would produce a 3% ($1\sigma$ of the modeling error) shift of $H_0$ towards a higher value for this lens, keeping all else constant. This work substantiates the groundbreaking potential of JWST for time-delay cosmography and lays the groundwork for modeling systems previously too faint to provide meaningful constraints on $H_0$.

We develop a sample of 234 dusty star-forming galaxies (DSFGs) from the A2744 and GOODS-S fields using JWST/NIRCam-selected galaxies as priors for SCUBA-2 measurements. This provides a large number of galaxies both above an 850 $\mu$m flux of 2 mJy (47 bright DSFGs) and below (187 faint DSFGs). It represents the largest sample of individually identified (i.e., not stacked) faint DSFGs to date. We identify a tight negative correlation between redshift and both F444W and F150W flux, suggesting that the observed NIR flux may be an effective way of selecting high-redshift DSFGs. We study the physical properties of the DSFGs through spectral energy distribution fitting and morphological analysis. Other than the lower star formation rates (SFRs) and total infrared luminosities in the faint DSFGs, the two populations have similar properties. The stellar masses do not appear to be strongly dependent on either the SFRs or the submillimeter flux. These results suggest that the faint DSFGs are drawn from the same population of galaxies as the bright DSFGs. We find a lower merger fraction (~21%) relative to previous HST-based studies.

Characterizing the evolution of the star-forming main sequence (SFMS) at high redshift is crucial to contextualize the observed extreme properties of galaxies in the early Universe. We present an analysis of the SFMS and its scatter in the THESAN-ZOOM simulations, where we find a redshift evolution of the SFMS normalization scaling as $\propto (1+z)^{2.64\pm0.03}$, significantly stronger than is typically inferred from observations. We can reproduce the flatter observed evolution by filtering out weakly star-forming galaxies, implying that current observational fits are biased due to a missing population of lulling galaxies or overestimated star-formation rates. We also explore star-formation variability using the scatter of galaxies around the SFMS ($\sigma_{\mathrm{MS}}$). At the population level, the scatter around the SFMS increases with cosmic time, driven by the increased importance of long-term environmental effects in regulating star formation at later times. To study short-term star-formation variability, or ''burstiness'', we isolate the scatter on timescales shorter than 50 Myr. The short-term scatter is larger at higher redshift, indicating that star formation is indeed more bursty in the early Universe. We identify two starburst modes: (i) externally driven, where rapid large-scale inflows trigger and fuel prolonged, extreme star formation episodes, and (ii) internally driven, where cyclical ejection and re-accretion of the interstellar medium in low-mass galaxies drive bursts, even under relatively steady large-scale inflow. Both modes occur at all redshifts, but the increased burstiness of galaxies at higher redshift is due to the increasing prevalence of the more extreme external mode of star formation.

The simplest inflationary models predict the primordial power spectrum (PPS) of curvature perturbations to be nearly scale-invariant. However, various other models of inflation predict deviations from this behaviour, motivating a data-driven approach to reconstruct the PPS and constrain its shape. In this work, we present a novel method that employs a fully differentiable pipeline to reconstruct the PPS using Gaussian Processes and uses neural network emulators for fast and differentiable theoretical predictions. By leveraging gradient-based sampling techniques, such as Hamiltonian Monte Carlo, our approach efficiently samples the high-dimensional parameter space of cosmological parameters and the free-form PPS, enabling joint constraints on both. Applying this framework to Planck 2018 Cosmic Microwave Background (CMB) temperature anisotropy data we find our reconstructed PPS to be consistent with near scale-invariance on small scales, while exhibiting large uncertainties at large scales, driven mostly by cosmic variance. Our results show an overestimation of the PPS amplitude compared to $\Lambda$CDM predictions from the Planck 2018 analysis, which we attribute to our choice of a conservative prior on the optical depth $\tau$ based on Planck 2015 measurements. Adopting a prior consistent with Planck 2018 measurements brings our results into full agreement with previous work. To ensure robustness of our results, we validate our differentiable pipeline against a non-differentiable framework, and also demonstrate that our results are insensitive to the choice of Gaussian process hyperparameters. These promising results and the flexibility of our pipeline make it ideally suited for application to additional data sets such as CMB polarisation as well as Large-Scale Structure probes, thus moving towards multi-probe primordial power spectrum reconstruction.

We present a machine learning search for local, low-mass galaxies ($z < 0.02$ and $10^6 M_\odot < M_* < 10^9 M_\odot$) using the combined photometric data from the DESI Imaging Legacy Surveys and the WISE survey. We introduce the spectrally confirmed training sample, discuss evaluation metrics, investigate the features, compare different machine learning algorithms, and find that a 7-class neural network classification model is highly effective in separating the signal (local, low-mass galaxies) from various contaminants, reaching a precision of $95\%$ and a recall of $76\%$. The principal contaminants are nearby sub-$L^*$ galaxies at $0.02 < z < 0.05$ and nearby massive galaxies at $0.05 < z < 0.2$. We find that the features encoding surface brightness information are essential to achieving a correct classification. Our final catalog, which we make available, consists of 112,859 local, low-mass galaxy candidates, where 36,408 have high probability ($p_{\rm signal} > 0.95$), covering the entire Legacy Surveys DR9 footprint. Using DESI-EDR public spectra and data from the SAGA and ELVES surveys, we find that our model has a precision of $\sim 100\%$, $96\%$, and $97\%$, respectively, and a recall of $\sim 51\%$, $68\%$ and $53\%$, respectively. The results of those independent spectral verification demonstrate the effectiveness and efficiency of our machine learning classification model.

Based on their cosmological origin, the stars of a galaxy can be divided into two categories: those that enter through merger events (ex situ) and those born in the main progenitor (in situ). We used the TNG50 cosmological magnetohydrodynamical simulation and its Lagrangian tracer particles to explore and quantify the origin of gas that ultimately forms the in situ stars of galaxies. We tracked back the baryonic mass contributing to the $z=0$ in situ stellar populations of galaxies, studying trends in mass from dwarfs to group-scale halos. We find that more massive halos acquire this matter earlier than lower-mass halos, reflecting an overall earlier assembly of their in situ stellar mass. Defining the Lagrangian half-mass radius R$_{\rm L, 1/2}$ of a galaxy as the distance containing half of the mass that will form its in situ stars by $z=0$, we find that R$_{\rm L, 1/2}$ is larger for more massive halos at early times, reflecting larger "in situ Lagrangian regions." However, the dependence of this radius on halo mass becomes flat at $z \simeq 3$ and then inverts toward $z=0$. In addition, R$_{\rm L, 1/2}$ increases rapidly with redshift, surpassing the virial radii of halos at $z \sim 2$. This marks the cosmic epoch at which most of the gas that eventually forms the in situ stars of galaxies leaves the intergalactic medium (IGM) and enters halos, a transition that occurs earlier for more massive halos. The formation redshift of the in situ stellar component increases with halo mass, while the formation redshift of the dark matter halo decreases, indicative of a differential assembly history between these two components. Finally, we decomposed the $z=0$ in situ stellar mass into its distinct modes of accretion. Smooth accretion from the IGM is the most important for low-mass galaxies, while mergers and satellite-stripped gas become relevant and even dominant only for high-mass galaxies.

We present the first follow-up with JWST of radio-selected NIRfaint galaxies as part of the COSMOS-Web survey. By selecting galaxies detected at radio frequencies ($S_{\rm 3 GHz}>11.5$ $\mu$Jy; i.e. S/N$>5$) and with faint counterparts at NIR wavelengths (F150W$>26.1$ mag), we collect a sample of 127 likely dusty star-forming galaxies (DSFGs). We estimate their physical properties through SED fitting, compute the first radio luminosity function for these types of sources, and their contribution to the total cosmic star formation rate density. Our analysis confirms that these sources represent a population of highly dust-obscured ($\langle A_{\rm v} \rangle \sim3.5$ mag), massive ($\langle M_\star \rangle \sim10^{10.8}$ M$_\odot$) and star-forming galaxies ($\langle {\rm SFR} \rangle\sim300$ M$_\odot$ yr$^{-1}$) located at $\langle z \rangle\sim3.6$, representing the high-redshift tail of the full distribution of radio sources. Our results also indicate that these galaxies could dominate the bright end of the radio luminosity function and reach a total contribution to the cosmic star formation rate density equal to that estimated only considering NIR-bright sources at $z\sim4.5$. Finally, our analysis further confirms that the radio selection can be employed to collect statistically significant samples of DSFGs, representing a complementary alternative to the other selections based on JWST colors or detection at FIR/(sub)mm wavelengths.

The integral-field unit mode of the Near-Infrared Spectrograph (NIRSpec+IFU) mounted on the James Webb Space Telescope has now enabled kinematic studies of smaller and less massive compact stellar systems in which to search for central massive black holes (BHs) than ever before. We present here the first such detection using NIRSpec+IFU in its highest resolution (R~2700) mode. We report a $3\sigma$ detection of a central black hole with mass ${\cal M}_{BH}=2.2\pm1.1\times10^6\,M_\odot$ in UCD736 orbiting within the Virgo galaxy cluster based on Schwarzschild's modeling of the 1D kinematic profile. The presence of such a massive BH strongly argues against a globular cluster origin of this UCD, and rather suggests a tidally stripped formation route from a former $\gtrsim10^9\,M_\odot$ dwarf galaxy host. Two other methods produce results consistent with Schwarzschild's modelling, but can only provide upper-limits on ${\cal M}_{BH}$. This represents the detection of a BH in the most compact ($r_h\approx15\,{\rm pc}$) stellar system to date, with a ${\cal M}_{BH}$ corresponding to ~9 percent of the system's stellar mass, roughly in line with previously reported UCD BH detections and comparable to the BH detected in the compact elliptical galaxy NGC4486B.

One of the main challenges in modeling massive stars to the onset of core collapse is the computational bottleneck of nucleosynthesis during advanced burning stages. The number of isotopes formed requires solving a large set of fully-coupled stiff ordinary differential equations (ODEs), making the simulations computationally intensive and prone to numerical instability. To overcome this barrier, we design a nuclear neural network (NNN) framework with multiple hidden layers to emulate nucleosynthesis calculations and conduct a proof-of-concept to evaluate its performance. The NNN takes the temperature, density and composition of a burning region as input and predicts the resulting isotopic abundances along with the energy generation and loss rates. We generate training sets for initial conditions corresponding to oxygen core depletion and beyond using large nuclear reaction networks, and compare the predictions of the NNNs to results from a commonly used small net. We find that the NNNs improve the accuracy of the electron fraction by $280-660\:\%$ and the nuclear energy generation by $250-750\:\%$, consistently outperforming the small network across all timesteps. They also achieve significantly better predictions of neutrino losses on relatively short timescales, with improvements ranging from $100-10^{6}\:\%$. While further work is needed to enhance their accuracy and applicability to different stellar conditions, integrating NNN trained models into stellar evolution codes is promising for facilitating large-scale generation of core-collapse supernova (CCSN) progenitors with higher physical fidelity.

The Einstein Telescope (ET), a proposed next-generation gravitational wave (GW) observatory, will expand the reach of GW astronomy of stellar-mass compact object binaries to unprecedented distances, enhancing opportunities for multi-messenger observations. Here we investigate multi-messenger emission properties of binary neutron star (NSNS) and black hole-neutron star (BHNS) mergers detectable by ET, providing projections to optimize observational strategies and maximize scientific insights from these sources. Using a synthetic population of compact binary mergers, we model each source's GW signal-to-noise ratio, sky localization uncertainty, kilonova (KN) light curves in optical and near-infrared bands, fluence of the relativistic jet gamma-ray burst (GRB) prompt emission and afterglow light curves across radio, optical, X-ray and very high energy wavelengths. We analyze multi-messenger detectability prospects for ET as a standalone observatory with two different configurations and within a network of next-generation GW detectors. ET will detect over $10^4$ NSNS mergers annually, enabling potential observation of tens to hundreds of electromagnetic (EM) counterparts. BHNS mergers have more limited multi-messenger prospects, but joint GW-EM rates will increase by an order of magnitude compared to current-generation instruments. We quantify uncertainties due to the NS equation of state (EoS) and mass distribution of NSNSs, as well as the NS EoS and BH spin for BHNSs. While a single ET will achieve an impressive GW detection rate, the fraction of well-localized events is orders of magnitude lower than in a network with additional detectors. This significantly limits efficient EM follow-up and science cases requiring well-characterized counterparts or early observations. The challenge is even greater for BHNS mergers due to their low EM rate.

We present the COSMOS Spectroscopic Redshift Compilation encompassing ~ 20 years of spectroscopic redshifts within the 2 deg$^2$ COSMOS legacy field. This compilation contains 165,312 redshifts of 97,929 unique objects from 108 individual observing programs up to $z \sim 8$ with median stellar mass $\sim 10^{9}$ to $10^{10}$ M$_\odot$ (redshift dependent). Rest-frame $NUVrJ$ colors and SFR -- stellar mass correlations show the compilation primarily contains low- to intermediate-mass star-forming and massive, quiescent galaxies at $z < 1.25$ and mostly low-mass bursty star-forming galaxies at $z > 2$. Sources in the compilation cover a diverse range of environments, including protoclusters such as ``Hyperion''. The full compilation is 50\% spectroscopically complete by $i \sim 23.2$ and $K_s \sim 21.3$ mag; however, this is redshift dependent. Spatially, the compilation is $>50$\% complete within the CANDELS area, while the outer regions of COSMOS are $>10$\% complete limited to $i < 24$ mag and $K_S < 22.5$ mag, separately. We demonstrate how the compilation can be used to validate photometric redshifts and investigate calibration metrics. By training self-organizing maps on COSMOS2020/Classic and projecting the compilation onto it, we find key galaxy subpopulations that currently lack spectroscopic coverage including $z < 1$ intermediate-mass quiescent galaxies and low-/intermediate-mass bursty star-forming galaxies, $z \sim 2$ massive quiescent galaxies, and $z > 3$ massive star-forming galaxies. This approach highlights how combining self-organizing maps with our compilation can provide guidance for future spectroscopic observations to get a complete spectroscopic view of galaxy populations. Lastly, the compilation will undergo periodic data releases that incorporate new spectroscopic redshift measurements, providing a lasting legacy resource for the community.

The difference between the total neutrino mass estimates ($\sum m_\nu$) derived from cosmological data within the standard $\Lambda$CDM model and those obtained from terrestrial particle physics experiments underscores the need to explore alternative scenarios. Recent analyses have shown that a dynamic dark energy modeled by the CPL parameterization of the dark energy equation of state (EoS) can ease the constraints on $\sum m_\nu$, thus alleviating this tension. This study investigates the robustness of this discrepancy by assessing the extent to which the CPL assumption influences the results. We examine how other EoS parameterizations, such as the Barboza-Alcaniz (BA) and Jassal-Bagla-Padmanabhan (JBP) parameterizations, affect $\sum m_\nu$ estimates. We perform a Markov Chain Monte Carlo (MCMC) analysis combining the latest baryon acoustic oscillation data from DESI with the Planck 2018 cosmic microwave background data, which includes information on temperature, polarization, and lensing, as well as the Pantheon+ Type Ia supernovae observations. While both the BA and JBP parameterizations can also resolve the tension, our results show a correlation between the dark energy EoS and the constraints on neutrino mass.

Recent observations by the James Webb Space Telescope have revealed massive galaxies at very high redshift ($z\simeq 7-15$). The question of whether the existence of such galaxies is expected in the corresponding JWST surveys has received a lot of attention, though the answer straddles areas of cosmology and complex astrophysical details of high-redshift galaxy formation. The growth rate of density fluctuations determines the amplitude of overdensities that collapse to form galaxies. Late-time modifications of growth, combined with measurements at both $z\sim 1$ from large-scale structure and $z\sim 1000$ from the cosmic microwave background, affect the predictions for the abundance of first galaxies in the universe. In this paper, we point out that the late-time growth rate of structure affects the statistical significance of high-redshift, high-mass objects very weakly. Consequently, if the existence and abundance of these objects are confirmed to be unexpected, the variations in the late-time growth history are unlikely to explain these anomalies.

Neural posterior estimation (NPE), a type of amortized variational inference, is a computationally efficient means of constructing probabilistic catalogs of light sources from astronomical images. To date, NPE has not been used to perform inference in models with spatially varying covariates. However, ground-based astronomical images have spatially varying sky backgrounds and point spread functions (PSFs), and accounting for this variation is essential for constructing accurate catalogs of imaged light sources. In this work, we introduce a method of performing NPE with spatially varying backgrounds and PSFs. In this method, we generate synthetic catalogs and semi-synthetic images for these catalogs using randomly sampled PSF and background estimates from existing surveys. Using this data, we train a neural network, which takes an astronomical image and representations of its background and PSF as input, to output a probabilistic catalog. Our experiments with Sloan Digital Sky Survey data demonstrate the effectiveness of NPE in the presence of spatially varying backgrounds and PSFs for light source detection, star/galaxy separation, and flux measurement.

The Tully-Fisher relation (TFR) is one of the most important and widely used empirical correlations in extragalactic astronomy. Apart from its importance as a secondary distance indicator, the TFR relation serves as a test for galaxy evolution models, because it connects the baryonic and dark matter components of galaxies. We aim at simulating the multi-wavelength TFR relation from UV to mid-infrared wavelengths for the TNG50 cosmological simulation at $z = 0$, and at comparing the results with observational TFR studies. We want to compare the wavelength dependence of the slope and scatter with the observed values, and search for secondary parameters that reduce the scatter in the TFR. We select a large sample of simulated late-type, disc-dominated galaxies from the TNG50 simulation. For each galaxy, we use the SKIRT radiative transfer code to generate realistic synthetic global fluxes in 12 UV to mid-infrared broadbands and synthetic integrated HI line profiles. We use bivariate linear regression to determine the TFR in each band, and we search for a second TFR parameter by correlating the residuals with different physical parameters. Our TNG50 TFR reproduces the characteristic behaviour of the observed TFR in many studies: the TFR becomes steeper and tighter as we move from UV/optical to infrared wavelengths. The slope changes from $-7.46 \pm 0.14~{\text{mag}}~{\text{dex}}^{-1}$ in the NUV band to $-9.66 \pm 0.09~{\text{mag}}~{\text{dex}}^{-1}$ in the IRAC [4.5] band. Quantitatively, our slopes are well within the spread of different observational results. The ${\textit{u}} - {\textit{r}}$ colour or the sSFR can significantly reduce the scatter in the UV and optical bands. Using ${\textit{u}} - {\textit{r}}$ colour as second parameter, the modified TFR has a roughly constant intrinsic tightness of over the entire UV to MIR range. The combination of the TNG50...

We present a catalog of uniformly processed 3.6-$\mu$m and 4.5-$\mu$m band exoplanet thermal phase curves based on Infrared Array Camera observations obtained from the Spitzer Heritage Archive. The catalog includes phase curve measurements for 34 planets, 16 of which contain full orbit coverage and have detectable secondary eclipses in both channels. The data are processed in the EXCALIBUR pipeline using a uniform analysis consisting of aperture photometry and modeling of instrument effects along with the exoplanet signal. Nearest-neighbors regression with a Gaussian kernel is used to correct for instrumental systematics correlated to the star's centroid position and shape in conjunction with a novel test to avoid overfitting. These methods may have utility in addressing sub-pixel gain variations present in modern infrared detectors. We analyze the 3.6-$\mu$m and 4.5-$\mu$m phase curve properties and find a strong wavelength-dependent difference in how the properties correlate with physical parameters as well as evidence that the phase curve properties are determined by multiple physical parameters. We suggest that differences between the 3.6-$\mu$m and 4.5-$\mu$m phase curve properties are due to 3.6~$\mu$m observations probing regions of the atmosphere which could include a cloud layer. Taken together, the observed phase curve behavior suggests that different physical processes are responsible for establishing the thermal phase curve at different pressures, which are probed by different wavelengths, and that further 3D GCM modeling is required to investigate the reason for this complex dependence on planetary properties.

Constraining the mass-sheet degeneracy (MSD) is crucial for improving the precision and accuracy of time-delay cosmography. Joint analyses of lensing and stellar kinematics have been widely adopted to break the MSD. A 3D mass and stellar tracer population is required to accurately interpret the kinematics data. We aim at forward-modeling the projection effects in strong lensing and kinematics observables, and finding the best model assumption for the stellar kinematics analysis which leads to unbiased interpretation of the MSD. We numerically simulate the projection and selection effects for both a triaxial ETG sample from the IllustrisTNG simulation and an axisymmetric sample which matches the properties of slow-rotator galaxies representative of the strong lens galaxy population. Using the axisymmetric sample, we generate mock kinematics observables with spherically-aligned axisymmetric Jeans Anisotropic Modeling (JAM) and assess kinematic recovery under different model assumptions. Using the triaxial sample, we quantify the random uncertainty introduced by modeling triaxial galaxies with axisymmetric JAM. We show that a spherical JAM analysis of spatially unresolved kinematic data introduces a bias of up to 2%-4% (depending on the intrinsic shape of the lens) in the inferred MSD. Our model largely corrects this bias, resulting in a residual random uncertainty in the range of 0-2.1% in the stellar velocity dispersion (0-4.2% in $H_0$) depending on the projected ellipticity and the anisotropy of the stellar orbits. This residual uncertainty can be further mitigated using spatially resolved kinematic data which constrain the intrinsic shape. We also show that the random uncertainty in the velocity dispersion recovery using axisymmetric JAM for axisymmetric galaxies is at the level of < 0.17%, and the uncertainty using axisymmetric JAM for triaxial galaxies is at the level of < 0.25%.

The stable rotation of young pulsars is often interrupted by two non-deterministic phenomena: glitches and red timing noise. Timing noise provides insights into plasma and nuclear physics under extreme conditions. The framework leverages rotational symmetry in pulsar spin-down models and temporal symmetry in noise processes to achieve computational efficiency, aligning with the journal's focus on symmetry principles in physical systems. In this paper, we apply a novel frequentist framework developed within the PINT software package (v0.9.8) to analyze single-pulsar noise processes. Using 17.5 years of pulse time-of-arrival (TOA) data for the young pulsar PSR J1741-3016, observed with the Nanshan 26 m radio telescope, we investigate its timing properties. In this study, we employed the Downhill Weighted Least-Squares Fitter to estimate the pulsar's spin parameters and position. The Akaike Information Criterion (AIC) was used for model parameter selection. The results obtained with PINT were compared to those from ENTERPRISE and TEMPONEST, two Bayesian-based frameworks. We demonstrate that PINT achieves comparable results with significantly reduced computational costs. Additionally, the adequacy of the noise model can be readily verified through visual inspection tools. Future research will utilize this framework to analyze timing noise across a large sample of young pulsars.

The conservation of stellar actions is a fundamental assumption in orbit reconstruction studies in the Milky Way. However, the disc is highly dynamic, with time-dependent, non-axisymmetric features like transient spiral arms and giant molecular clouds (GMCs) driving local fluctuations in the gravitational potential on top of the near-axisymmetric background. Using high-resolution magnetohydrodynamic simulations that incorporate gas dynamics and star formation, we quantify the rate at which these effects drive non-conservation of the actions of young stars from Myr to Gyr timescales. We find that action evolution is well described as a logarithmic random walk, with vertical action evolving more rapidly than radial action; the diffusion rate associated with this random walk is weakly dependent on the stellar birth environment and scales approximately linearly with the galactic orbital frequency at a star's position. The diffusion rates we measure imply a fundamental limit of $\sim 100$ Myr as the timescale over which stellar orbits can be reliably reconstructed using methods that assume action conservation. By comparing diffusion rates for younger stars to those measured for an older and more vertically-extended control population, we conclude that radial action evolution is driven primarily by transient spiral arms, while vertical action evolution is driven by gravitational scattering off gaseous structures. Our results have significant implications for galactic archaeology and disc dynamics studies, necessitating a closer look at the timescales over which actions are assumed to be conserved in the disc.

The physical origins of the diverse emission-line asymmetries observed in the spectra of active galactic nuclei (AGNs) remain incompletely understood. Monitoring the temporal variations of line profiles offers a promising approach to investigating the underlying physics. In this study, we present an analysis of the broad H$\beta$ emission line profiles of eight AGNs observed from the end of 2016 to May 2023 as part of the reverberation mapping campaign titled "Monitoring AGNs with H$\beta$ Asymmetry" (MAHA), utilizing data obtained from the Wyoming Infrared Observatory (WIRO) 2.3-meter telescope. We measure the temporal variations of line asymmetry, width, and central velocity shift for the eight objects. Our findings reveal that the variation in asymmetry is positively correlated with H$\beta$ flux in five of the eight objects, while the remaining objects exhibit negative or complex correlations. Furthermore, we observe anti-correlations between line width and H$\beta$ flux for most objects, indicating the presence of the "breathing" phenomenon in their H$\beta$ emission lines. In contrast, two objects demonstrate an "anti-breathing" phenomenon or complex behavior. We discuss the physical origins of the temporal variations in line profiles and propose the possibility of decomposing the variations in H$\beta$ asymmetry and width into components: one that corresponds to short-term variations in H$\beta$ flux and another that reflects long-term variations in continuum light curves, perhaps driven by radiation pressure.

This review talk explores the diverse outcomes of white dwarf mergers, emphasising that not all double degenerate mergers result in supernovae. Possible outcomes include the formation of a more massive WD, partial explosions that leave behind unusual remnants, and merger-induced collapse leading to the formation of rapidly spinning neutron stars or magnetars. Additionally, merging WDs have been proposed as potential progenitors for long gamma-ray bursts and Fast Radio Bursts.

We present Surveying the Whirlpool at Arcseconds with NOEMA (SWAN), a high-resolution, high-sensitivity survey to map molecular lines in the 3mm band in M51 (the Whirlpool galaxy). SWAN has obtained the largest high-sensitivity map (5x7 kpc2) of N2H+ emission at cloud-scale resolution (3" ~125 pc) in an external galaxy to date. We describe the observations and data reduction of ~214 hours of interferometric data from NOEMA, ~55 hours of tailored new observations with the IRAM-30m telescope and the combination of NOEMA, new and ~14 hours of archival 30m observations. We detect widespread emission from 9 molecular transition lines. The J=1-0 transitions of CO isotopologues 13CO and C18O are detected at high significance across the full observed field-of-view (FoV). HCN, HNC, HCO+, and N2H+(1-0) are detected in the center, molecular ring and spiral arms of the galaxy, while the shock tracer HNCO(4-3), (5-4) and PDR tracer C2H(1-0) are detected in the central ~1 kpc and molecular ring only. For most of the lines that we detect, average line ratios with respect to CO are increased by up to a factor of ~3 in the central 1 kpc, where an AGN and its low-inclination outflow are present, compared to the disk. Across the full SWAN FoV, 13CO, C18O, HCN, HNC, HCO+ and N2H+ are 8\pm2, 29\pm6, 17\pm3,37\pm5, 26\pm5 and 63\pm38 times fainter than 12CO, respectively, in pixels where each line is significantly detected. Although we observe variations in line ratios between larger-scale environments like the center and disk of M51, the scatter within each environment also indicates the influence of smaller-scale processes. The ability to measure these effects is only possible thanks to the high resolution and high sensitivity of the SWAN dataset across multiple environments. This provides the sharpest view of these molecular transitions over the largest physical area ever captured in an external galaxy.

The effect of tidal forces on transport within a relic accretion disk in binary black holes is studied here with a suite of two-dimensional hydrodynamic simulations. As the binary contracts due to the emission of gravitational waves, the accretion disk is truncated, and a two-armed spiral wave is excited, which remains stationary in the rotating reference frame of the coalescing binary. Such spiral waves lead to increased transport of mass and angular momentum. Our findings suggest that even in the case of weakly ionized accretion disks, spiral density waves will drain the disk long before the orbit of the two black holes decays enough for them to merge, thus dimming prospects for a detectable electromagnetic counterpart.

$\gamma$ Persei is a long-period ($P \approx 14.6$ yr) eclipsing binary system. Its period makes it a difficult target to fully understand: so far, only two primary eclipses are known in the literature, from 1990 and from 2019, whereas the 2005 one was missed due to its closeness to the Sun at the time. We aimed to fill in this gap by processing the quasi-continuous photometry collected by the Solar Mass Ejection Imager (SMEI) between 2003 and 2011, which was ideally positioned to observe such a bright targets. In order to do that, we first determined a color-dependent conversion formula from the SMEI measurements into {\it Gaia} \textit{G} magnitudes. We applied various corrections to the photometry and provide the longest continuous light curve of $\gamma$ Persei. We successfully detected the 2005 primary eclipse of the system, with the yearly observations ending during the egress of the companion. We predicted the position of a possible secondary eclipse by forward modeling the binary system with PHOEBE, and successfully recovered the secondary eclipse in the 2006 SMEI observations. The existence of the secondary eclipse puts strong constraints on the orbital configuration, which will be an important constraint for future studies of the system.

On electron kinetic scales, ions and electrons decouple, and electron velocity shear on electron inertial length $\sim d_e$ can trigger electromagnetic (EM) electron Kelvin-Helmholtz instability (EKHI). In this paper, we present an analytic study of EM EKHI in an inviscid collisionless plasma with a step-function electron shear flow. We show that in incompressible collisionless plasma the ideal electron frozen-in condition $\mathbf{E} + \mathbf{v}_e \times \mathbf{B}/c = 0$ must be broken for the EM EKHI to occur. In a step-function electron shear flow, the ideal electron frozen-in condition is replaced by magnetic flux conservation, i.e., $\nabla \times (\mathbf{E} + \mathbf{v}_e\times \mathbf{B}/c) = 0$, resulting in a dispersion relation similar to that of the standard ideal and incompressible magnetohydrodynamics KHI. The magnetic field parallel to the electron streaming suppresses the EM EKHI due to magnetic tension. The threshold for the EM mode of the EKHI is $(\mathbf{k}\cdot\Delta\mathbf{U}_e)^2>\frac{n_{e1}+n_{e2}}{n_{e1} n_{e2}}[n_{e1}(\mathbf{v}_{Ae1}\cdot\mathbf{k})^2+n_{e2}(\mathbf{v}_{Ae2}\cdot\mathbf{k})^2]$, where $\mathbf{v}_{Ae} =\mathbf{B}/(4\pi m_e n_e)^{1/2}$, $\Delta\mathbf{U}_e$ and $n_e$ are the electron streaming velocity shear and densities, respectively. The growth rate of the EM mode is $\gamma_{em} \sim \Omega_{ce}$, the electron gyro-frequency.

A classical nova is a thermonuclear runaway initiated on a white dwarf accreting solar-like material from its stellar companion. Once the white dwarf accretes enough mass, the pressure at the base of the accreted layer reaches a critical point, leading to the ignition of the hydrogen fuel at their interface. This paper presents a set of two-dimensional CO classical nova simulations with an extended buffer zone of a fixed low density and temperature between the top of the accreted layer and the upper boundary, allowing us to capture the thermonuclear outburst in the domain. Our domain reduces the role of the upper-outflow boundary condition that has affected previous simulations and allows us to explore the nucleosynthesis evolution in detail. We also study the effects of the initial temperature perturbation and buffer size to explore their sensitivity in our simulations. Finally, we start our simulations with a lower temperature at the base of the accreted layer ($7\times 10^7\,\mathrm{K}$) than previous work, allowing us to capture mixing earlier in the evolution, reducing the effects of the mixing-length-theory assumptions. This allows for a more realistic description of convective transport in our models.

The study of primordial black holes (PBHs) and the Hawking radiation they produce represents a significant step toward understanding the role of these phenomena in the cosmological evolution of the Universe. PBHs could be a component of dark matter, the seeds of supermassive black holes, and sources of Hawking radiation, which, unlike the radiation from other black holes, might be observable. Over the course of the evolution of the Universe from the Big Bang to the present day, PBHs have lost a substantial portion of their mass through neutrino radiation. This occurs because, for black holes with masses $M\lesssim10^{23}\,$g, the radiation of the lightest massive particles - neutrinos - becomes significant alongside the emission of photons. This neutrino radiation eventually dominates the emission process. By the present time ($t_0=13.8$ Gyr), only black holes with masses $\lesssim 10^{15}\,$g have largely evaporated, meaning that the total emission spectrum of PBHs is now dominated by the neutrino component. In this paper, we present new estimates of the neutrino spectra emitted by PBHs of various masses, with a particular focus on low-energy radiation ($E_{kin} \in [0.01 ÷1]\, $eV) for the first time. Our calculations demonstrate that black holes in the mass range $[10^{9}\div10^{11}]\,$g emit neutrinos with intensities that exceed the background fluxes from known astrophysical sources in the low-energy range. In contrast, in the high-energy range, the emission remains below the background level, consistent with observational constraints. These findings open new avenues for the potential detection of PBH radiation and could stimulate advancements in neutrino detection technologies, particularly in the low-energy regime. The observation of neutrinos in this energy range represents one of the few opportunities to confirm the existence of Hawking radiation.

Classical Cepheids are not only excellent standard candles, but also invaluable tools to test stellar evolution and pulsation theories. Rates of their pulsation period change, quantified usually through $O-C$ diagrams, can be confronted with predictions of stellar evolution theory. On the other hand, period changes on much shorter time scales ($\sim$10$^{2}$-10$^{4}$days), attributed to non-evolutionary effects are often detected and lack detailed explanation. We aim to provide a systematic and quantitative description of irregular or non-linear period changes in Cepheids. We aim to provide a systematic and quantitative description of irregular or non-linear period changes in Cepheids. We analysed part of the OGLE data for classical Cepheids in the Magellanic Clouds (MCs; from both Large Magellanic Cloud, LMC, and the Small Magellanic Cloud, SMC) using the modified Hertzsprung $O-C$ technique. A sample of 3658 stars, with the best quality data and void of additional low-amplitude periodicities (e.g. due to non-radial pulsations), that could impact the results, was selected for analysis. Based on $O-C$ shapes, stars were classified into three categories: no period change (class 1), linear period change (class 2), and irregular change (class 3). In our investigation, $33.5\pm0.7$\% of analysed stars show irregular period changes. Considering the pulsation mode, irregular period changes were detected in $16.5\pm0.7$\% of the analysed fundamental mode stars and in $68.1\pm1.2$\% of the first overtone stars. The amplitude of variability in the $O-C$ diagrams increases with the pulsation period, and at a given pulsation period, it is larger for first overtone stars. While the increase is linear for first overtone stars, for fundamental mode stars it becomes steeper as the pulsation period increases.

MOMOS, the Multi-Object MKID Optical Spectrometer, is a proposed visible wavelength spectrometer that uses MKIDs (Microwave Kinetic Inductance Detectors) targeting an initial resolving power of 3500. With their modest wavelength-resolving abilities, MKIDs take the place of both the cross disperser and detector in the spectrometer. MKIDs lack read noise and dark current enabling noiseless post-observation rebinning and characterization of faint objects, as well as time-resolved photon counting spectroscopy. This work presents an MOMOS simulator customizable for different MOMOS configurations. Treating simulator products as inputs, an algorithm was developed and implemented in the MOMOS data reduction package to calibrate and extract spectra.

A major outstanding challenge in cosmology is the persistent discrepancy between the Hubble constant obtained from early and late universe measurements -- the Hubble tension. Examining cosmological evolution through the lens of information growth within a black hole we show the appearence of two fractal growing processes characterizing the early and late ages. These fractals induce space growth rates of $62.79\pm5.59$ km/s/Mpc and $70.07\pm0.09$ km/s/Mpc; close to the current values of the Hubble constants involved in the tension. These results strongly suggest that the Hubble tension is not given by unexpected large-scale structures or multiple, unrelated errors but by innate properties underlying the universe dynamics.

In a space weather context, the most geoeffective coronal mass ejections (CMEs) are fast CMEs from Earth-facing solar active regions. These CMEs are difficult to characterize in coronagraph data due to their high speed (fewer observations), faintness, Earthward orientation (halo CMEs), and disruptions from associated high-energy particle storms. Any diagnostic aiding in early CME speed identification is valuable. This study investigates whether the 3D speeds of 37 CMEs are correlated with the critical heights of their source regions, to test the hypothesis that if the critical height is located at a higher altitude in the corona, the weaker magnetic field environment will enable a faster CME to be produced. Critical heights near CME onset are calculated by identifying polarity inversion lines (PIL) in magnetogram data using automated and manual methods. 3D speeds are determined by fitting a Graduated Cylindrical Shell (GCS) model to multi-viewpoint coronagraph images. For the automated method, we find a high correlation of 71% +- 8% between CME speed and critical height, dropping to 48% +- 12% when using CME plane-of-sky speeds, on which most previous similar studies are based. An attempt to improve the critical height diagnostic through manual PIL selection yields a lower correlation of 58% +/- 13%. The higher correlation from the automated method suggests that encompassing the full PIL structure is a better measure of the magnetic conditions that influence CME dynamics. Our results highlight the potential for critical height as a continuously computable diagnostic for forecasting the 3D speeds of Earth-directed CMEs.

The first severe (G4) geomagnetic storm of Solar Cycle 25 occurred on 23-24 April 2023, following the arrival of a Coronal Mass Ejection (CME) on 23 April. The characteristics of this CME, measured from coronagraphs (speed and mass), did not indicate that it would trigger such an intense geomagnetic storm. In this work, our aim is to understand why this CME led to such a geoeffective outcome. Our analysis spans from the source active region to the corona and inner heliosphere through 1 au using multiwavelength, multi-viewpoint remote sensing observations and in situ data. We find that rotation and possibly deflection of the CME resulted in an axial magnetic field nearly parallel to the ecliptic plane during the Earth encounter, which might explain the storm's severity. Additionally, we find that imaging away from the Sun-Earth line is crucial in hindcasting the CME Time-of-Arrival at Earth. The position (0.39 au) and detailed images from the SoloHI telescope onboard the Solar Orbiter mission, in combination with SOHO and STEREO images, helped decisively with the three-dimensional (3D) reconstruction of the CME.

Motivated by the morphological measures in assessing the geometrical and topological properties of a generic cosmological stochastic field, we propose an extension of the weighted morphological measures, specifically the $n$th conditional moments of derivative (cmd-$n$). This criterion assigns a distinct weight to each excursion set point based on the associated field. We apply the cmd-$n$ on the Cosmic Microwave Background (CMB) to identify the cosmic string networks (CSs) through their unique Gott-Kaiser-Stebbins effect on the temperature anisotropies. We also formulate the perturbative expansion of cmd-$n$ for the weak non-Gaussian regime up to $\mathcal{O}(\sigma_0^3)$. We propose a comprehensive pipeline designed for analyzing the morphological properties of string-induced CMB maps within the flat sky approximation. To evaluate the robustness of our proposed criteria, we employ string-induced high-resolution flat-sky CMB simulated patches of $7.2$ deg$^2$ size with a resolution of $0.42$ arc-minutes. Our results demonstrate that the minimum detectable value of cosmic string tension is $G\mu\gtrsim 1.9\times 10^{-7}$ when a noise-free map is analyzed with normalized cmd-$n$. Whereas for the ACT, CMB-S4, and Planck-like experiments at 95.45\% confidence level, the normalized cmd-$n$ can distinguish the CSs network for $G\mu\gtrsim2.9 \times 10^{-7}$, $G\mu\gtrsim 2.4\times 10^{-7}$ and $G\mu\gtrsim 5.8\times 10^{-7}$, respectively. The normalized cmd-$n$ exhibits a significantly enhanced capability in the detection of CSs relative to the Minkowski Functionals.

Human living environment is influenced by intense solar activity. The solar activity exhibits periodicity and regularity. Although many deep-learning models are currently used for solar cycle prediction, most of them are based on a multi-step pattern. In this paper a solar cycle prediction method based on a one-step pattern is proposed with the TCN neural network model, in which a number of historical data are input, and only one value is predicted at a time. Through an autoregressive strategy, this predicted value is added to the input sequence to generate the next output. This process is iterated until the prediction of multiple future data. The experiments were performed on the 13-month smoothed monthly total sunspot number data sourced from WDC-SILSO. The results showed that one-step pattern fits the solar cycles from 20-25 well. The average fitting errors are MAE=1.74, RMSE=2.34. Finally, the intensity of Solar Cycle 25 was predicted with one-step pattern. The peak will occur in 2024 October with a magnitude of 135.3 and end in 2030 November. By comparing the prediction results with other methods, our method are more reasonable and better than the most methods. The codes are available on \href{this https URL} {github} and \href{this https URL

Next-generation instruments for ground-based gamma-ray astronomy are marked by a substantial increase in complexity, featuring dozens of telescopes. This leap in scale introduces significant challenges in managing system operations and offline data analysis. Methods, which depend on advanced personnel training and sophisticated software, become increasingly strained as system complexity grows, making it more challenging to effectively support users in such a multifaceted environment. To address these challenges, we propose the development of AI agents based on instruction-finetuned large language models (LLMs). These agents align with specific documentation and codebases, understand the environmental context, operate with external APIs, and communicate with humans in natural language. Leveraging the advanced capabilities of modern LLMs, which can process and retain vast amounts of information, these AI agents offer a transformative approach to system management and data analysis by automating complex tasks and providing intelligent assistance. We present two prototypes that integrate with the Cherenkov Telescope Array Observatory pipelines for operations and offline data analysis. The first prototype automates data model implementation and maintenance for the Configuration Database of the Array Control and Data Acquisition (ACADA). The second prototype is an open-access code generation application tailored for data analysis based on the Gammapy framework.

A unique galaxy at z = 2.2, zC406690, has a striking clumpy large-scale ring structure that persists from rest UV to near-infrared, yet has an ordered rotation and lies on the star-formation main sequence. We combine new JWST/NIRCam and ALMA band 4 observations, together with previous VLT/SINFONI integral field spectroscopy and HST imaging to re-examine its nature. The high-resolution H$\alpha$ kinematics are best fitted if the mass is distributed within a ring with total mass $M_{\rm{ring}} = 2 \times 10^{10} M_\odot$ and radius $R_{ring}$ = 4.6 kpc, together with a central undetected mass component (e.g., a "bulge") with a dynamical mass of $M_{bulge} = 8 \times 10^{10} M_\odot$. We also consider a purely flux emitting ring superposed over a faint exponential disk, or a highly "cuspy" dark matter halo, both disfavored against a massive ring model. The low-resolution CO(4-3) line and 142GHz continuum emission imply a total molecular and dust gas masses of $M_{mol,gas} = 7.1 \times 10^{10}M_\odot$ and $M_{dust} = 3 \times 10^8 M_\odot$ over the entire galaxy, giving a dust-to-mass ratio of 0.7%. We estimate that roughly half the gas and dust mass lie inside the ring, and that $\sim 10\%$ of the total dust is in a foreground screen that attenuates the stellar light of the bulge in the rest-UV to near-infrared. Sensitive high-resolution ALMA observations will be essential to confirm this scenario and study the gas and dust distribution.

We investigated the impact of photometric redshift errors in the ongoing Javalambre Physics of the Accelerating Universe Astrophysical Survey (J-PAS) on void identification and properties using a watershed-based method, aiming to assess the recovery of individual voids and the overall void environment. We created galaxy mock catalogues for redshift z = 0.1 using the IllustrisTNG300-1 simulation, defining two datasets: an $ideal$ sample ($m_r < 21$ mag) and a $perturbed$ sample with the Z-coordinate errors mimicking J-PAS's line-of-sight errors, derived from the precursor miniJPAS survey data. We identified voids using ZOBOV, a watershed algorithm. We found 1065 voids in the $ideal$ sample and 2558 voids in the $perturbed$ sample. The $perturbed$ sample voids have, on average, smaller sizes and denser interiors. We filtered out voids based on density and radius in order to eliminate overdense and small spurious instances. The stacked density profile of filtered voids in the $perturbed$ sample remains close to the average density even at the boundary peak, indicating a strong blurring of structures by the redshift errors. The number of $ideal$ sample voids for which at least $50\%$ of the volume is recovered by a void in the $perturbed$ sample is 53 (29 for the filtered sample). The volume occupied by these voids is less than $10\%$ of the simulation volume. Merging voids in the $perturbed$ sample marginally improves the recovery. The overall volumes defined as voids in the two samples have an overlap of $80\%$, making up $61\%$ of the simulation box volume. While some statistical properties of voids might be recovered sufficiently well, the watershed algorithms may not be optimal for recovering the large-scale structure voids if applied straight to photometric redshift survey data.

The WASP-47 system is notable as the first known system hosting both inner and outer low-mass planetary companions around a hot Jupiter, with an ultra-short-period (USP) planet as the innermost planetary companion. The formation of such an unique configuration poses challenges to the lonely hot Jupiter formation model. Hot Jupiters in multiple planetary systems may have a similar formation process with warm Jupiter systems, which are more commonly found with companions. This implies that the WASP-47 system could bridge our understanding of both hot and warm Jupiter formation. In this work, we propose a possible formation scenario for the WASP-47 system based on its orbital configuration. The mean motion resonance trapping, giant planet perturbations, and tidal effects caused by the central star are key factors in the formation of USP planets in multiple planetary systems with hot Jupiters. Whether a planet can become an USP planet or a short period super-earth (SPSE) planet depends on the competition between eccentricity excitation by nearby giant planet perturbations and the eccentricity damping due to tidal effects. The $Q_p'$ value of the innermost planet is essential for the final planetary configuration. Our results suggest that a $Q_p'$ in the range of [1, 10] is favorable for the formation of the WASP-47 system. Based on the formation scenario, we estimate an occurrence rate of 8.4$\pm$2.4\% for USP planets in systems similar to WASP-47.

The history of the outer solar system is intrinsically related to the Giant Planets migration. A massive disk of material within a radius of 30~au was scattered during the planetary migration, creating different dynamic populations in the Transneptunian region. They were formed in a collisional environment when massive collisions allowed them to grow and form much smaller moons than the primary body. The dynamical group, known as the Cold Classicals, was formed in a sparse disk from 42 to about 47~au and did not suffer much from planet migration. Observations show that many of Cold Classical are binary, consistent with the streaming instability process. The stellar occultation technique, with a spatial resolution of a few kilometres, can be used to search for binaries where other techniques are unable to do so, and to characterise the known satellites of Trans-Neptunian Objects (TNO), constraining their formation scenarios. We review here the first stellar occultations by TNO's satellites (besides Charon), discuss the methods used to detect these events. We also fit new orbital elements and system mass for Vanth (Orcus/1) and Weywot (Quaoar/1), finding reasonable solutions for pure Keplerian orbits. Finally, we discuss the prospects regarding the stellar occultations by TNO binaries and their implications for the study of the history of the Solar System.

In this manuscript, based on SDSS photometry and spectroscopy, a method is proposed to test the hypothesis that the orbital kinematics of kpc-scale dual-core systems can lead to double-peaked narrow emission lines (DPNELs), through analyzing a sample of seven kpc-scale dual-core systems by comparing the upper limits of their orbital velocities (calculated using total stellar mass and projected distance) with the velocity separation of DPNELs (peak separation). To determine accurate total stellar masses, GALFIT is applied to consider the effects of the overlapping components on the photometric images to obtain accurate magnitudes. Then, based on the correlation between absolute Petrosian magnitudes and total stellar masses, the individual masses of the galaxies are determined. Therefore, the maximum orbital velocities can be calculated by combining the projected distance. Meanwhile, the peak separation can be accurately measured after subtracting the pPXF determined host galaxy contributions. Finally, four objects exhibit peak separations almost consistent with their respective maximum orbital velocities under the assumption of a circular orbit, while the remaining three objects display peak separations larger than the maximum orbital velocities. A larger sample will be given later to further test whether DPNELs can arise from kpc-scale dual-core systems.

We present so far the highest resolution ($\sim$0.04") ALMA 1.3 mm continuum observations of three massive star-forming clumps in the Central Molecular Zone, namely 20 km s$^{-1}$ C1, 20 km $^{-1}$ C4, and Sgr C C4, which reveal prevalent compact millimeter emission. We extract the compact emission with $\textit{astrodendro}$ and identify a total of 199 fragments with a typical size of $\sim$370 AU, which represent the first sample of candidates of protostellar envelopes and disks and kernels of prestellar cores in these clumps that are likely forming star clusters. Compared with the protoclusters in the Galactic disk, the three protoclusters display a higher level of hierarchical clustering, likely a result of the stronger turbulence in the CMZ clumps. Compared with the mini-starbursts in the CMZ, Sgr B2 M and N, the three protoclusters also show stronger subclustering in conjunction with a lack of massive fragments. The efficiency of high-mass star formation of the three protoclusters is on average one order of magnitude lower than that of Sgr B2 M and N, despite a similar overall efficiency of converting gas into stars. The lower efficiency of high-mass star formation in the three protoclusters is likely attributed to hierarchical cluster formation.

Radio telescope arrays, such as Square Kilometre Array (SKA) and next-generation Very Large Array (ngVLA), require highly precise synchronization of time-frequency references to ensure high-quality observational data. Fiber-based frequency distribution systems are highly effective. However, their proper functioning can be threatened by risk events. In this paper, we propose a hybrid fiber-based frequency-distribution and vibration detection system tailored for large radio arrays. The system ensures the performance of distributed frequency signals while allowing for the monitoring of potential threats to the optical fiber network. We design and implement a single-to-multiple hybrid system, conducting tests via a 55-km fiber link. Experimental results demonstrate its effectiveness, achieving the relative frequency stability of 3E-14/1 s and 2.7E-17/1E5 s, along with vibration detection and localization capabilities.

We report on the intriguing properties of a variable X-ray source projected at the outskirts of the elliptical galaxy NGC 6099 ($d \approx 139$ Mpc). If truly located near NGC 6099, this is a hyperluminous X-ray source that reached an X-ray luminosity $L_{X} \approx $ a few times $10^{42}$ erg s$^{-1}$ in 2012 February (XMM-Newton data), about 50 to 100 times brighter than in 2009 May (Chandra) and 2023 August (XMM-Newton). The X-ray spectrum was soft at all three epochs, with a thermal component at $kT \approx 0.2$ keV and a power-law photon index $>3$. Such properties make it a strong candidate for an intermediate mass black hole (IMBH). We also discovered a point-like, blue optical counterpart ($m_{g,{Vega}}\approx24.7$~mag, $M_{g,{Vega}}\approx-11.2$~mag), from images taken by the Canada-France-Hawaii Telescope, and later confirmed with Hubble Space Telescope observations. The optical continuum can be modeled as stellar emission from a compact star cluster or an X-ray-irradiated accretion disk, consistent with the IMBH scenario. We discuss alternative explanations for the nature of this system. A possible scenario is tidal stripping of an orbiting star, with repeated X-ray outbursts every few years. An alternative possibility is that the thermal X-ray emission seen in 2009 was from shocked gas in the self-intersecting tidal stream during the rising phase of a tidal disruption event, while the 2012 and 2023 emissions were from the fully-formed accretion disk.

The Largest Cluster Statistics\,(LCS) analysis of the redshifted 21\,cm maps has been demonstrated to be an efficient and robust method for following the time evolution of the largest ionized regions\,(LIRs) during the Epoch of Reionization\,(EoR). The LCS can, in principle, constrain the reionization model and history by quantifying the morphology of neutral hydrogen\,(\HI) distribution during the different stages of the EoR. Specifically, the percolation transition of ionized regions, quantified and constrained via LCS, provides a crucial insight about the underlying reionization model. The previous LCS analysis of EoR 21\,cm maps demonstrates that the convolution of the synthesized beam of the radio interferometric arrays, e.g. SKA1-Low with the target signal, shifts the apparent percolation transition of ionized regions towards the lower redshifts. In this study, we present an optimal thresholding strategy to reduce this bias in the recovered percolation transition. We assess the robustness of LCS analysis of the 21\,cm maps in the presence of antenna-based gain calibration errors and instrumental noise for SKA1-Low. This analysis is performed using synthetic observations simulated by the \textsc{21cmE2E} pipeline, considering SKA1-Low AA4 configuration within a radius of 2\,km from the array centre. Our findings suggest that a minimum of $1500$\,hours of observation (SNR $\gtrapprox 3$) are required for the LCS analysis to credibly suppress the confusion introduced by thermal noise. Further, we also demonstrate that for a maximum antenna-based calibration error tolerance of $\sim 0.05\%$ (post calibration), the reionization history can be recovered in a robust and relatively unbiased manner using the LCS.

HD 45166 was recently reported to be a long-period binary comprising a B7V star and a highly magnetic ($\langle B \rangle = 43.0\pm0.5\,$kG) hot Wolf-Rayet-like component, dubbed as a quasi Wolf-Rayet (qWR) star in literature. While originally proposed to be a short-period binary, long-term spectroscopic monitoring suggested a 22.5 yr orbital period. With a derived dynamical mass of $2.03\pm0.44\,M_\odot$, the qWR component is the most strongly magnetized non-degenerate object ever detected and a potential magnetar progenitor. However, the long period renders the spectroscopic orbital solution and dynamical mass estimates uncertain, casting doubts on whether the qWR component is massive enough to undergo core-collapse. Here, we spatially resolve the HD 45166 binary using newly acquired interferometric data obtained with the GRAVITY instrument of the Very Large Telescope Interferometer. Due to the calibrator star being a binary as well, we implement a new approach for visibility calibration and test it thoroughly using archival GRAVITY data. The newly calibrated HD 45166 data reveal the unmistakable presence of a companion to the qWR component with an angular separation of $10.9\pm0.1$ mas (which translates to a projected physical separation of $10.8\pm0.4$ au), consistent with the long-period orbit. We obtain a model-independent qWR mass $M_{\rm qWR} = 1.96^{+0.74}_{-0.54}\,M_\odot$ using interferometric and spectroscopic data together. This observation robustly confirms that HD 45166 is truly a long-period binary, and provides an anchor point for accurate mass determination of the qWR component with further observations.

JWST has revealed a population of broad-line active galactic nuclei at $z>4$ with remarkably red colors, so-called "Little Red Dots." Ubiquitous Balmer breaks suggest that they harbor old stellar populations in massive, compact host galaxies. We present ALMA observations of three LRDs at $z=3.10$, $4.46$, and $7.04$, targeting molecular and neutral gas via CO(7-6) and [CI](2-1), respectively. We do not detect CO in any target, placing conservative limits on the host molecular gas mass $\lesssim 1$-$5\times10^{10}$ M$_\odot$. We report the tentative ($4.9\sigma$) detection of the [CI](2-1) line in A2744-45924 ($z=4.46$), one of the brightest known LRDs. The [CI] line is narrow (FWHM $\sim 80$ km s$^{-1}$), implying a dynamical mass $\lesssim 10^{10}$ M$_\odot$, adopting conservative limits for the galaxy size. The dynamical mass limit is significantly lower than expected from the local $M_{\rm BH}$-$M_{\rm dyn}$ relation, and is an order of magnitude below the stellar mass derived from SED fitting, potentially supporting a non-stellar origin of the Balmer break. These results, while tentative, paint a picture of LRDs that is markedly different than typical high-$z$ quasars, which live in massive, gas-rich, and actively star-forming host galaxies.

The existence of dynamically young and metal-rich RR Lyrae stars challenges conventional notions of these variable stars. One possible scenario for their formation and evolution is via binary channels involving mass transfer. This study presents the detection of nine fundamental-mode RR Lyrae stars residing in the thin disk of the Milky Way with metallicities higher than [Fe/H] > -1.0 dex and showing proper motion anomalies. Our thin disk classification is based on kinematics and supported by $\alpha$-element abundances, where possible. We searched for indications of the light-travel time effect (LTTE) in the available literature sources and the TESS photometric data of the stars but found no signs of periodic variations induced by companions within the expected period range. This could be because of a lack of observations as well as sparse measurements and large gaps in the data. We propose a continued search for signs of binarity and a subsequent long-term follow-up of nine targets that satisfy all of our search criteria. Beyond these targets, we also report the detection of slow phase changes in the Blazhko star ST Pic, which could be compatible with the LTTE.

Mars' sedimentary rocks record Gyrs of environmental change. New data enable the first global analysis of paleo-environment relevant physical properties of these rocks, including layer thickness and accumulation rate. We find that layer thicknesses of post-3.5 Ga sedimentary rocks across the Martian surface show coherent variations at ~1000 km-scale that are inconsistent with simple volcanic and climatic hypotheses for formation, which are consistent with global compositional homogeneity at orbital scales. These data, in combination with new analyses of outcrop age and total rock volume demonstrate a global decrease in layer thickness that predates the eventual drop off in preserved sedimentary rock volume per Myr. The new constraints confirm a diachronous transition in Mars' global sedimentary rock record while also highlighting a regional dichotomy in young sedimentary rock deposits that has not been quantified before.

Aims. The orientation-based unification scheme of radio-loud active galactic nuclei (AGNs) asserts that radio galaxies and quasars are essentially the same type of object, but viewed from different angles. To test this unification model, we compared the environments of radio galaxies and quasars, which would reveal similar properties when an accurate model is utilized. Methods. Using the second data release of the LOFAR Two-metre Sky Survey (LoTSS DR2), we constructed a sample of 26,577 radio galaxies and 2028 quasars at 0.08 < z < 0.4. For radio galaxies with optical spectra, we further classified them as 3631 low-excitation radio galaxies (LERGs) and 1143 high-excitation radio galaxies (HERGs). We crossmatched these samples with two galaxy cluster catalogs from the Sloan Digital Sky Survey (SDSS). Results. We find that $17.1 \pm 0.2%$ of the radio galaxies and $4.1 \pm 0.4%$ of the quasars are associated with galaxy clusters. Luminous quasars are very rare in clusters, while $18.7 \pm 0.7%$ LERGs and $15.2 \pm 1.1%$ HERGs reside in clusters. We also note that in radio galaxies, both HERGs and LERGs tend to reside in the centers of clusters, while quasars do not show a strong preference for their positions in clusters. Conclusions. This study shows that local quasars and radio galaxies exist in different environments, challenging the orientation-based unification model. This means that factors other than orientation may play an important role in distinguishing radio galaxies from quasars. The future WEAVE-LOFAR survey will offer high-quality spectroscopic data for a large number of radio sources and allow for a more comprehensive exploration of the environments of radio galaxies and quasars.

Low-Luminosity Active Galactic Nuclei (LLAGN) provides a unique view of Comptonization and non-thermal emission from accreting black holes in the low-accretion rate regime. However, to decipher the exact nature of the Comptonizing corona in LLAGN, its geometry and emission mechanism must be understood beyond the limits of spectro-timing techniques. Spectro-polarimetry offers the potential to break the degeneracies between different coronal emission models. Compton-thin LLAGN provide an opportunity for such spectro-polarimetric exploration in the 2-8 keV energy range using IXPE. In this work, we carry out a spectro-polarimetric analysis of the first IXPE observation, in synergy with a contemporaneous NuSTAR observation, of an LLAGN: NGC 2110. Using 554.4 ks of IXPE data from October 2024, we constrain the 99% upper limit on the Polarization Degree (PD) to be less than 8.3% assuming the corresponding Polarization Angle (PA) to be aligned with the radio jet, and less than 3.6% if in the perpendicular direction. In the absence of a significant PD detection, the PA remains formally unconstrained, yet the polarization significance contours appear to be aligned with the radio jet, tentatively supporting models in which the corona is radially extended in the plane of the disk. We also carry out detailed Monte Carlo simulations using MONK and STOKES codes to test different coronal models against our results and compare the polarization properties between NGC 2110 and brighter Seyferts.

We compare Very Large Array observations of GRS 1915+105 made in 1994 and 2023, with nearly three decades of difference. The source has experienced intriguing major changes. The position angle of the bipolar ejecta in the plane of the sky has increased counterclockwise by 24 degrees. The inclination angle of the flow with respect to the line of sight has increased by 17 degrees. Analysis of GRS 1915+105 images over the years suggest that the observed changes took place within a year or less. Our analysis indicates that during 2023 the plane of the accretion disk was aligned with the line of sight, which may explain the deep X-ray obscured state and the high mid-infrared luminosity observed with JWST in that epoch. More recent 2024 observations imply that the position angle of the ejecta has returned to its historic values. We suggest that these abrupt changes could be due to the presence of an undetected tertiary component in the system. Future monitoring of the time evolution of the source may further clarify the cause of these remarkable changes.

The streaming instability is a promising mechanism for planetesimal formation. The instability can rapidly form dense clumps that collapse self-gravitationally, which is efficient for large dust grains with the Stokes number on the order of 0.1. However, dust growth models predict that collisional fragmentation prevents dust grains from growing to such sizes. We perform local simulations of the streaming instability and measure characteristic collision velocities and collision rates of dust grains based on their trajectories in moderate clumping. The collision velocities are on the order of 0.1 percent of the sound speed or lower, implying that dust grains can overcome the fragmentation barrier via the clumping. We also find that the collision rates are appreciably high regardless of the low collision velocities. Corresponding timescales are on the order of ten Keplerian periods or shorter, suggesting that dust grains can overcome the drift barrier as well. This streaming-instability-assisted (SI-assisted) coagulation greatly relaxes the conditions for planetesimal formation as recently implied.

Neptune-size exoplanets are less studied as characterizing their atmospheres presents challenges due to their relatively small radius and atmospheric scale height. As the most common outcome of planet formation, these planets are crucial for understanding planetary formation, migration theories, atmospheric composition, and potential habitability. Their diverse atmospheres, influenced by equilibrium temperature, composition, and cloud presence, offer unique opportunities to study atmospheric dynamics and chemistry. While low-resolution spectroscopy struggles with atmospheric characterization due to clouds, high-resolution observations provide detailed analysis of the atmospheres by detecting molecular lines beyond the cloud deck. This study investigates four subclasses of Neptune atmospheres: HAT-P-11 b (warm Neptune), HD 63433 c (warm sub-Neptune), K2-25 b (temperate Neptune), and TOI-270 d (temperate sub-Neptune), using six ground-based spectrographs: GIANO-B, CARMENES, IGRINS, HISPEC, MODHIS, and ANDES over one and three transits. Our simulation integrates the chemical kinetics model, VULCAN with the 1-D line-by-line radiative transfer model, petitRADTRANS, and estimates detection significance using the ground-based noise simulator, SPECTR. We aim to predict how future Extremely Large Telescopes (ELTs) such as TMT (MODHIS) and E-ELT (ANDES) can utilize their higher resolving powers and larger collecting areas to surpass current observatories in detecting molecular bands. We highlight the importance of photochemistry in these atmospheres and demonstrate how ELTs will help further in constraining nitrogen and sulfur chemistry. Finally, we present a comprehensive picture of cloud presence in the atmospheres and its impact on molecular detectability in Neptune-class atmospheres.

It is still in dispute the existence of intermediate-mass black holes (IMBHs) with a mass of ~10^3-10^5 solar masses (Msun), which are the missing link between stellar-mass black holes (5-50 Msun) and supermassive black holes (10^6-10^10 Msun). The bright flares from tidal disruption events (TDEs) provide a new and direct way to probe IMBHs. 3XMM J215022.4-055108 is a unique off-nuclear X-ray transient which can be best explained as the TDE by an IMBH in a massive star cluster, though its mass is not well determined. Here, we report the discovery of a transient X-ray quasi-periodicity signal from 3XMM J215022.4-055108 with a period of ~85 seconds (at a significance of >3.51sigma) and fractional root-mean-squared amplitude of ~10%. Furthermore, the signal is coherent with a quality factor ~16. The significance drops to >3.13sigma if considering all light curves with sufficient quality for QPO search. Combining with the results from X-ray continuum fittings, the detection of QPO allows for joint constraints on the black hole mass and dimensionless spin in the range [9.9*10^3-1.6*10^4 Msun]$ and [0.26-0.36], respectively. This result supports the presence of an IMBH in an off-nuclear massive star cluster and may open up the possibility of studying IMBHs through X-ray timing of TDEs.

Turbulence is a mysterious phenomenon in physical systems and plays a critical role in the interstellar medium (ISM). Previous simulations and observations have shown that the probability density functions (PDFs) of gas densities in supersonic systems tend to be skewed, exhibiting low-density exponential tails on the $s = \ln{\rho}$ scale. In this work, we argue that these exponential tails originate from the convolution of PDF kernels with skewed tails, which can appear at both ends of the PDF. We introduce two density-fraction strategies -- the mass-fraction and volume-fraction approaches -- to explain the physical origins of the low-$s$ and high-$s$ skewed PDF kernels. The PDF kernels constructed in this work satisfy the relation that its structural entropy equals the PDF variance, as we proposed in Paper I. These two types of PDF kernels define two spaces of skewed PDFs, which are dual to each other and possess highly symmetric mathematical structures. We thus speculate that the high-$s$ skewed PDF kernels are physical and may be related to the power-law tails of the column-density PDFs (on the $\rho$ scale) of molecular clouds, which are shaped by both turbulence and gravity. Inspired by this, we further construct a form of an `isothermal' turbulent system that likely favors the volume-fraction strategy.

Stellar-mass isolated black holes (IsoBHs) wandering in interstellar medium (ISM) are expected to be abundant in our Galaxy. Recently, an IsoBH, OGLE-2011-BLG-0462, was unambiguously discovered using astrometric microlensing. We examine prospects for detecting electromagnetic signatures from an accretion flow surrounding the IsoBH. The accretion rate onto the IsoBH should be highly sub-Eddington, which leads to formation of a hot accretion flow. In this paper, we evaluate the detectability of electromagnetic signals from the hot accretion flows in two accretion states: magnetically arrested disk (MAD) and classical radiatively inefficient accretion flows (RIAFs). For the MAD scenario, we find that the optical, infrared, and X-ray signals can be detectable by the current best facilities, such as HST, JWST, and Chandra, if the IsoBH is in a warm neutral medium. In contrast, for the classical RIAF scenario, the optical and X-ray emissions are weaker than MAD scenario, leading to unobservable signals for a typical parameter set. Future follow-up observations of OGLE-2011-BLG-0462 will provide a good test for theories of accretion processes.

It has always been believed that feedback from active galactic nuclei (AGN) has an important impact on star formation in massive galaxies. Black hole spin is an important physical parameter of AGN. We use a large sample of massive star-forming galaxies to study the effects of AGN on star formation. Our main results are as follows: (i) There are significant correlations between black hole spin and star formation rate, specific star formation rate, and star formation activity parameter for massive star-forming early-type and late-type galaxies, respectively. These results indicate that the spin of supermassive black holes regulates the star formation of massive star-forming early-type and late-type galaxies. (2) The slopes of the relationship between black hole spin and star formation rate, specific star formation rate, and star formation activity parameter for massive star-forming early-type galaxies and late-type galaxies are similar within the error range. These results imply that the mechanism of black hole spin regulating star formation may be similar in massive star-forming early-type and late-type galaxies.

A plateau on the background inflaton potential $V_{\rm b}(\phi)$ can lead cosmic inflation into an ultraslow-roll phase, greatly enhancing the primordial power spectrum on small scales, and resulting in significant scalar-induced gravitational waves (GWs) and abundant primordial black holes (PBHs). In this work, we construct an anti-symmetric perturbation $\delta V(\phi)$ on $V_{\rm b}(\phi)$ with three model parameters, the position, width, and slope of $\delta V(\phi)$, and constrain these parameters from the potential stochastic gravitational wave background (SGWB) in NANOGrav 15-year data set. The GW spectrum from the merger of supermassive black hole binaries (SMBHBs) with two model parameters, the amplitude and spectral index, is also investigated for comparison. We perform Bayesian analysis in three steps with increasing number of model parameters, and obtain the allowed parameter ranges. When the constraints on PBH abundance are taken into account, these ranges become further narrower. We find that the increase of model parameters cannot significantly improve the Bayes factors, and the model with an almost perfect plateau on $V_{\rm b}(\phi)$ is favored. Moreover, the interpretation of the SGWB via only the GWs generated by SMBHBs is not preferred by the data.

Recent JWST observations towards Westerlund 1 revealed extensive nebular emission associated with the cluster. Given the age of the region and proximity of that material to massive stars it cannot be primordial star forming gas and the origin is uncertain. We aim to determine whether the nebular emission in Westerlund 1 could be due to ablation flows from Red Supergiant (RSG) stars embedded in the cluster wind driven by the Wolf-Rayet stars in the cluster core. We also aim to explore the efficiency of mass-loading for the RSG wind in this scenario. We use 3D hydrodynamic simulations with the \textsc{pion} code to study the interaction between the cluster and RSG winds. We compare with the JWST observations by generating synthetic dust-emission maps. We find that the ablation flow morphology is consistent with the observations towards Westerlund 1, with clumps and instabilities. Synthetic observations at 11 $\mu$m predict fluxes in the ablation flow of $\sim1000-6000$ MJy ster$^{-1}$ which is consistent with the unsaturated components of the JWST F1130W observations in the vicinity of the red supergiants in the region. This good agreement is achieved without any consideration of polycyclic aromatic hydrocarbons (PAHs), which have a known 11.3 $\mu$m feature that appears in the F1130W band. This suggests that the ablation flow is PAH depleted. Ablation of RSG winds can explain the observed nebulosity in Westerlund 1, at least in the vicinity of the RSGs. Further observations are encouraged to enable detailed studies of these interactions.

Novae are found to have GeV to TeV gamma-ray emission, which reveals the shock acceleration from the white dwarfs. Recently, V1405 Cas was reported to radiated suspicious gamma-ray by \textit{Fermi}-LAT with low significance ($4.1 \sigma$) after the optical maximum. Radio observations reveal that it is one of the five brightest novae surrounded by low-density ionized gas columns. Here we report continuous search for GeV gamma-ray from \textit{Fermi}-LAT. No gamma-ray were found. For V1405 Cas, the flux level is lower than other well-studied \textit{Fermi} novae, and the gamma-ray maximum appear at $t_{0} + 145$ d. Gamma-ray of V1405 Cas are used to search potential gamma-ray periodicity. No gamma-ray periodicity was found during the time of observation. By comparing multi-wavelength data, the gamma-ray upper limit to optical flux ratio with value at around $10^{-4}$ is obtained to constrain the shock acceleration. Long-term analysis from \textit{Swift}-XRT gets X-ray spectral in the post-shock phase, which indicates that V1405 Cas became a super-soft source. The best-fit black body temperature at the super soft state is 0.11 - 0.19 keV.

Core-collapse supernovae (CCSNe) are the explosive end-points of stellar evolution for $M_{ZAMS} \gtrsim 8$ $M_\odot$ stars. The cores of these stars collapse to neutron stars, a process in which high neutrino luminosity drives off the overlying stellar layers, which get ejected with thousands of kilometers per second. These supernovae enrich their host galaxies with elements made both during the star's life and in the explosion, providing the main cosmic source of elements such as oxygen, neon and silicon. Their high luminosities ($\sim$ $10^{42}$ erg s$^{-1}$ at peak) make SNe beacons to large distances, and their light curves and spectra provide rich information on single and binary stellar evolution, nucleosynthesis, and a diverse set of high-energy physical processes. As the SN ejecta sweep up circumstellar and interstellar matter, it eventually enters a supernova remnant phase, exemplified by nearby, spatially resolved remnants such as Cas A and the Crab Nebula. In this phase, shocks and pulsar winds continue to light up the interior of the exploded stars, giving detailed information about their 3D structure. We review the central concepts of CCSNe, from the late stages of evolution of massive stars, through collapse, explosion, and electromagnetic display, to the final remnant phase. We briefly discuss still open questions, and current and future research avenues.

This article explores the different formation scenarios of the Kronian moons system in the context of a highly dissipative Saturn, with the objective of identifying the most likely of these scenarios. First, we review the diversity of objects - moons and rings - orbiting solar system giant planets, and the diversity of their architectures, which formation scenarios must reproduce. We then identify in this broader context the specific features of the Saturn system, such as the particularly large spectrum of its moon masses, the uniqueness of Titan and the presence of both dense and tenuous rings, before discussing the applicability of the different giant planet moon formation scenarios to the Saturn case. We discuss each of the most relevant scenarios and their respective merits. Finally, we tentatively propose a "favorite" scenario and we identify the key observations to be made by future space missions and/or Earth-based telescopic observations to validate this scenario or possibly alternative ones.

Neutrinos from dense environments are unique laboratories for astrophysics, particle physics and many-body physics. They tell us about the last stages of the gravitational core-collapse and the explosion of massive stars. These elusive particles are also tightly linked to heavy elements synthesis in gravitational core-collapse supernovae and binary neutron star mergers, or play a pivotal role at the MeV epoch during the Universe expansion. We highlight theoretical and observational aspects of this interesting domain, in particular for the future measurement of neutrinos from the next core-collapse supernova, and of the diffuse supernova background, whose discovery might lie in the forthcoming future.

Magnetars have inferred polar field strengths in excess of the Schwinger limit, where nonlinear electromagnetic effects can be significant. Their internal fields may be even stronger, suggesting that Maxwellian characterisations of hydromagnetic structure may require revision. A generalised Grad-Shafranov equation, describing static and axisymmetric fluid stars with mixed poloidal-toroidal fields, is introduced and subsequently solved in a perturbative scheme to calculate quadrupolar deformations. In the Born-Infeld theory, we show that the toroidal field has a maximum strength set by the scale parameter, $b$, implying an upper limit to the stellar prolateness, $|\epsilon_{\rm max}| \sim 10^{-5} \left(b/10^{16}\text{ G}\right)^2$, that is independent of field specifics. Observations of magnetar phenomena that are interpreted as evidence for ellipticity, such as precession, can thus implicitly constrain post-Maxwellian parameters in a way that complements terrestrial experiments. Toroidal ceilings also have implications for dynamo theory and gravitational waves, which we revisit together with field evolution in crusts abiding by beyond-Maxwell physics.

The magnitudes of the velocity kicks that neutron stars (NSs) obtain at their formation have long been a topic of discussion, with the latest studies analysing the velocities of young pulsars and favouring a bimodal kick distribution. In previous work, a novel method was proposed to determine kicks based on the eccentricity of Galactic trajectories. We applied this method to the isolated pulsars with known parallax in order to kinematically constrain the NS natal kick distribution and investigate its proposed bimodality. Since this method is applicable to older pulsars, we effectively increase the sample size with ${\sim}50\%$ compared to the pulsars younger than $10$ Myr. We assumed the velocity vectors of the pulsars to be distributed isotropically in the local standard of rest frame. These velocity vectors were used to trace back the trajectories of the NSs through the Galaxy and estimate their eccentricity. Then, we simulated kicked objects in order to evaluate the relationship between kick magnitude and Galactic eccentricity, which was used to infer the kicks corresponding to the estimated eccentricities. The resulting kick distributions indeed show a bimodal structure for young pulsars. However, for older pulsars the bimodality vanishes and instead we find a log-normal kick distribution peaking at ${\sim}200$ km/s and a median of ${\sim}400$ km/s. We therefore conclude that the natal kicks of these isolated NSs are best described by log-normal distribution with $\mu=6.38$ and $\sigma=1.01$. This analysis reveals no evidence for bimodality in the larger sample, and we suggest that the bimodality found by existing literature may be caused by Poisson noise due to their relatively small sample size.

We study the amplification of primordial curvature perturbations through a sound speed resonance mechanism within an inflationary model featuring a non-minimal derivative coupling. In this scenario, the coupling function is composed of a constant term and a periodic term. Our analysis demonstrates that the periodic behavior of the squared sound speed transforms the curvature perturbation equation into a Mathieu equation on sub-horizon scales. This resonance effect leads to a substantial enhancement of curvature perturbations, which, upon re-entering the Hubble horizon during the radiation-dominated era, can trigger the formation of primordial black holes, potentially contributing to the dark matter density. Additionally, the amplified perturbations generate scalar-induced gravitational waves with amplitudes that may lie within the detectable range of future gravitational wave observatories, providing a testable observational signature of the model.

Recently, the R2D2 paradigm, standing for ''Residual-to-Residual DNN series for high-Dynamic-range imaging'', was introduced for image formation in Radio Interferometry (RI) as a learned version of the traditional algorithm CLEAN. The first incarnations of R2D2 are limited to planar imaging on small fields of view, failing to meet the spherical-imaging requirement of modern telescopes observing wide fields. To address this limitation, we propose the spherical-imaging extension S-R2D2. Firstly, as R2D2, S-R2D2 encapsulates its minor cycles in existing 2D-Euclidean deep neural network (DNN) architectures, but adapts its iterative scheme to incorporate the wide-field measurement model mapping a spherical image to visibility data. We implemented this model as the composition of an efficient Fourier-based interpolator mapping the spherical image onto the equatorial plane, with the standard RI operator mapping the equatorial-plane image to visibility data. Importantly, the interpolation step must inevitably be performed at a lower-than-optimal resolution on the plane, to meet the high-resolution requirement on the sphere of wide-field imaging while preserving scalability. Therefore, secondly, we design S-R2D2's DNN training loss to jointly learn to correct the interpolation approximations and identify residual image structures on the sphere, ensuring consistency with the spherical ground truth using the adjoint plane-to-sphere interpolator. Finally, we demonstrate through simulations S-R2D2's capability to perform fast and accurate reconstructions of spherical monochromatic intensity images, across high-resolution, high-dynamic-range settings.

A major challenge in CGM studies is determining the three-dimensional (3D) properties from the observed projected observations. Here, we decompose the 3D gas density and spatial distribution of cool clouds by fitting a cool CGM model with the absorption observations, including the cool gas density, Ly$\alpha$, and Mg II equivalent widths. The clumpiness in the cool CGM is considered by modeling individual clouds. This model has four major components: the radial profile of the cool gas density; the number density of clouds; the absorption properties within individual clouds; and the velocity dispersion in the CGM. The observed cool gas density exhibits a large dispersion of $\approx2-3$ dex within the virial radius ($r_{vir}$). This dispersion can be reproduced with a combination of the projection effect (i.e., distant low-density clouds projected at small radii) and the intrinsic variation in the gas density. By modeling the probability density functions of gas density at different radii, the cool gas density is modeled as a $\beta$-model with a characteristic gas density of $\log n_{H,0}/{\rm cm^{-3}}=-2.57_{-0.25}^{+0.43}$ at $r_{vir}$ and a slope of $\beta_c=0.63_{-0.20}^{+0.16}$, and the intrinsic dispersion is $\sigma_{n_{H}}\approx 0.56_{-0.20}^{+0.19}$ dex. Assuming a cloud mass of $10^4M_\odot$, we further constrain the number density of cool clouds by jointly reproducing Ly$\alpha$ and Mg II equivalent width samples, resulting into a number density of $\log n_{\mathcal{N}_{cl},0}/ r_{vir}^{-3}=4.76^{+0.27}_{-0.21}$ at $r_{vir}$ and a slope of $\beta_N=0.65^{+0.06}_{-0.07}$. This spatial distribution of the cool CGM leads to a total cool gas mass of $\log M_{cool}/M_\odot=10.01^{+0.06}_{-0.06}$ for $L^*$ galaxies, while varying the cloud mass from $10^3M_\odot$ to $10^6M_\odot$ leads to the total cool CGM mass of $9.62_{-0.07}^{+0.05}$ to $10.46_{-0.05}^{+0.05}$.

On January 15 2025, the Gaia mission completed the collection of the astrometric, photometric, and spectroscopic data for about 2.5 billion celestial sources, from the solar system to the Milky Way to the distant universe. Work is ongoing to produce Gaia DR4 based on the first 5.5 years of data, with the release expected in 2026. The full 10.5 year survey will be turned into Giaa DR5 which will open up scientific possibilities beyond Gaia DR4}. In this contribution I give a brief overview of the Gaia mission, summarize results from the GaiaUnlimited project, provide a glimpse of what is to come in Gaia DR4, and summarize the new science opportunities that Gaia DR5 will bring. I close with a look ahead at the successor to Gaia, the GaiaNIR mission, which will survey the Milky Way in the infrared, thus probing the Galactic ecosystem in the regions hidden to the Gaia mission.

Despite more than half a century of research, the dominant radiation mechanism of gamma-ray burst (GRB) prompt emission remains unsolved. Some progress has been made through the analyses of the observational spectra of Swift/BAT, Konus/Wind, and Fermi/GBM, as well as the spectra of the photosphere or synchrotron models, but it is still insufficient to pin down the answer. Recent and upcoming high-sensitivity polarization observations can provide additional instructive information for model evaluation. In this work, we thoughtfully investigate the polarization samples of POLAR and AstroSAT, combining the light curve, the spectral, and the polarization parameters. The power-law shape of the X-ray afterglows, the $T_{90}\propto (L_{\text{iso}})^{-0.5}$ correlation, and the hard low-energy spectral index $\alpha$ are revealed, thus supporting the photosphere origin. Furthermore, we discover the positive correlation of the $\alpha$ and the polarization degree (PD), which can be consistently explained by the photosphere polarization scenario involving the jet asymmetry from a moderate viewing angle of $\theta _{v}$=0.015.

The occurrence rate of giant planets increases with orbital period and turns over at a location that roughly corresponds to the snow line of solar-type stars. Further, the density distribution of cold Jupiters (CJs) on the semi-major axis - mass diagram shows a relatively steep inner boundary, shaping the desert of warm Jupiters. The eccentricities of CJs show a broad distribution with a decreasing number density towards the larger end. Previous planet formation models fail to reproduce all these features at the same time. We use a planet population synthesis (PPS) model with truncated initial planetesimal distribution and compare the mass and orbital distribution of the simulated planets with the observation. We show that the occurrence of CJs with respect to the orbital period, the slope of the inner boundary of CJs on the semi-major axis - mass diagram, and the eccentricity distribution of CJs agree reasonably well with observation, if CJs form from truncated planetesimal disks of 10 au or wider with suppressed migration. While PPS simulations generally overestimate the fraction of giants with eccentricity below 0.2, $N$-body simulations produce a more consistent eccentricity distribution with observation. While the fraction of high-eccentricity planets can be increased by widening the planetesimal disk or reducing the migration speed, a deficit of giants with eccentricity between 0.2-0.4 exists regardless of the choices of parameters. Our results indicate that CJs are more likely born in truncated disks near the snow line than in classical uniform disks.

SN 2024ggi is a Type II supernova that exploded in the nearby galaxy NGC 3621 at a distance of approximately 7 Mpc, making it one of the closest supernovae of the decade. This SN shows clear signs of interaction with a dense circumstellar material, and several studies have investigated the properties of its possible progenitor star using pre-explosion data. In this work we aim to constrain the progenitor properties of SN 2024ggi by performing hydrodynamical modeling of its bolometric light curve and expansion velocities. We present photometric and spectroscopic observations of SN 2024ggi obtained in the Complejo Astronómico El Leoncito, in Las Campanas Observatory, and in Las Cumbres Observatory Global Telescope Network, spanning from 2 to 106 days after explosion. We constructed its bolometric light curve and we characterize it by calculating its morphological parameters. Then, we computed a grid of one dimensional explosion models for evolved stars with varying masses and estimated the properties of the progenitor star of SN 2024ggi by comparing the models to the observations. The observed bolometric luminosity and expansion velocities are well-matched by a model involving the explosion of a star in the presence of a close circumstellar material (CSM), with a zero-age main sequence mass of $\mathrm{M_{ZAMS}}$ = 15 $M_{\odot}$, a pre-SN mass and radius of 14.1 $M_{\odot}$ and 517 $R_{\odot}$, respectively, an explosion energy of $1.3\times10^{51}$ erg, and a nickel mass below 0.035 $M_{\odot}$. Our analysis suggests that the progenitor suffered a mass-loss rate of $4 \times 10^{-3}$ $M_{\odot}$yr$^{-1}$, confined to a distance of 3000 $R_{\odot}$. The CSM distribution is likely a two-component structure that consists of a compact core and an extended tail. This analysis represents the first hydrodynamical model of SN 2024ggi with a complete coverage of the plateau phase.

We present a sample of 146 high-redshift ($z>7.5$) galaxies from the CANUCS/Technicolor surveys, showcasing photometry in every wide- and medium-band NIRCam filter in addition to ancillary HST data sampling $0.4-5 \mu m$ (22 JWST bands out of 29 bands total). Additionally, 48 ($33\%$) galaxies in our sample meet criteria to be classified as extreme emission line galaxies, 15 ($10\%$) of which are completely missed by typical dropout selections due to faint UV emission. By fitting the SEDs covering the rest-frame UV to optical at $z > 7.5$, we investigate the dust obscuration properties, giving an unbiased view of dust buildup in high-redshift galaxies free from spectroscopic follow-up selection effects. Dust attenuation correlates with stellar mass, but more strongly with star formation rate. We find typical galaxies at $z>7.5$ have $\sim 25 \%$ of their star formation obscured. However, since galaxies with higher star formation rates suffer more attenuation, $\sim 50 \%$ of the total star formation rate density at $7.5<z<9$ is obscured. The obscured fraction drops to $\sim 35 \%$ in our $9<z<12$ bin, possibly due to substantial dust grain growth in the interstellar medium not having time to occur. Extrapolating the decline in dust obscuration of galaxies to higher redshifts, we infer that dust obscuration should approach zero at $z > 15$, implying that epoch as when dust first forms in bright galaxies.

The Flat Spectrum Radio Quasar (FSRQ) B2\,1308+326 was in its highest $\gamma$-ray flaring state during 60260-60310\,MJD. During this period, the source was detected in very high energy (VHE) by the large-sized telescope (LST-1). We conducted a detailed broadband spectral study of this source using the simultaneous data available in optical/UV, X-ray, and $\gamma$-ray bands. For the broadband spectral study, we select two gamma-ray high flux states (59750-59800\,MJD, 60260-60310\,MJD) and one low flux state (59250-59320\,MJD). During the epochs, 59750-59800\,MJD (high flux state) and 59250-59320\,MJD (low flux state), the broadband spectral energy distribution (SED) is well fitted using one zone leptonic emission model involving synchrotron, synchrotron self Compton (SSC) and external Compton (EC) processes. However, the flaring state (60260-60310\,MJD) during which the source showed VHE emission requires an additional component. We show that the inclusion of the photo-meson process can successfully explain this excess $\gamma$-ray emission. Further the estimated parameters, also suggest the source is transparent to VHE gamma-rays against pair production process.

Of the ~25 directly imaged planets to date, all are younger than 500Myr and all but 6 are younger than 100Myr. Eps Ind A (HD209100, HIP108870) is a K5V star of roughly solar age (recently derived as 3.7-5.7Gyr and 3.5$^{+0.8}_{-1.3}$Gyr). A long-term radial velocity trend as well as an astrometric acceleration led to claims of a giant planet orbiting the nearby star (3.6384$\pm$0.0013pc). Here we report JWST coronagraphic images that reveal a giant exoplanet which is consistent with these radial and astrometric measurements, but inconsistent with the previously claimed planet properties. The new planet has temperature ~275K, and is remarkably bright at 10.65um and 15.50um. Non-detections between 3.5-5um indicate an unknown opacity source in the atmosphere, possibly suggesting a high metallicity, high carbon-to-oxygen ratio planet. The best-fit temperature of the planet is consistent with theoretical thermal evolution models, which are previously untested at this temperature range. The data indicates that this is likely the only giant planet in the system and we therefore refer to it as ``b", despite it having significantly different orbital properties than the previously claimed planet ``b".

We present the design and science case for a new array of radio antennas to be located at the Pierre Auger Observatory. Six stations of three SKALA antennas each will be deployed around a single water-Cherenkov surface detector triggering the radio readout. The planned antenna layout will allow for the detection of cosmic rays above a few tens of PeV and reach full efficiency for vertical air showers at several hundred PeV in primary energy. The array will thus be a pathfinder to demonstrate that fully-efficient radio detection in combination with the underground muon detectors already present at the location is possible. This will enable combined studies of the muon component and the depth of shower maximum, relevant for hadronic interaction models studies and more accurate determination of the cosmic-ray mass composition in the energy range of the Galactic-to-extragalactic transition.

The interaction between gas and dust in protoplanetary disks (PPDs) plays a crucial role in setting the stage of planet formation. In particular, the streaming instability (SI) is well recognized as the mechanism for planetesimal formation out of this interaction. The outer region of PPDs is likely subject to the vertical shear instability (VSI), representing a major source of disk turbulence characterized by vertical corrugation that leads to strong dust stirring. In the meantime, the VSI turbulence in 3D generates vortices through the Rossby wave instability (RWI), which can trap dust and thereby promote dust concentration. In this study, we use the multifluid dust module in Athena++ to conduct 2D axisymmetric global simulations of PPDs with mesh refinement and 3D global simulations with modest resolution. In 2D, the VSI corrugation mode is weakened by dust back-reaction, while the SI can still survive regardless of initial conditions. Dust clumping occurs and is seeded by VSI-induced zonal flows. In 3D, dust can settle even more with increased dusty buoyancy, suppressing the VSI corrugation mode. Meanwhile, dust back-reaction enhances dust concentration in RWI vortices, though higher resolution is needed to assess dust clumping.

The study of protoclusters at cosmic noon is essential to understand the impact on galaxies of the environment and of the transformational processes occurring in this epoch. This work tests the predictions of the DIANOGA cosmological hydrodynamical simulations of cluster progenitors at z=2.2, comparing them with observations, and investigates the environmental effects on galaxies by comparing protoclusters with an average volume of the Universe. We analyze 14 protoclusters and a cosmological box of 49 cMpc/h per side. We compare predictions and observations of the galaxy properties, including colors of galaxies obtained with radiative transfer, to analyze UVJ diagrams. We showed that the DIANOGA simulations produce a galaxy stellar mass function in broad agreement with observations, with a higher fraction of high-mass galaxies ($M_{\ast}>10^{10} \ M_{\odot}$) in massive halos in protoclusters, compared to the box. The same signal, with lower significance, is also observed in the wide-field protocluster structures, indicating an accelerated evolution of galaxies before their infall into massive halos. Our simulations underestimate SFRs of galaxies both in protoclusters and in the box, compared to observations, due to low gas reservoirs. We find a weak suppression of SFRs in protocluster galaxies (~0.05 dex), compared to the box, increasing up to ~0.25 dex in massive halos. The quenched galaxy fraction varies significantly across different protocluster halos, consistent with observations. The simulations show a strong dependence of quenched fractions on halo mass and an excess of quenched galaxies in the wide-field protocluster region, compared to the cosmological box. UVJ diagram analysis shows qualitative agreement with observed color distributions of star-forming and quenched galaxies, except for few massive galaxies with steeper reddening vectors than typically assumed in observations.

We present UV and X-ray observations, obtained with the Swift Observatory, of Gaia BH3 a binary system containing a 33 solar masses black hole discovered through Gaia astrometry. The system is well detected in all UV and optical filters (1700-6500 Å). We compare our results with the modeling of the non collapsed component using synthetic stellar libraries, a good agreement is found with our UV observations

Carbon dioxide (CO$_2$) is one of the most important interstellar molecules. While it is considered that it forms on the surface of interstellar dust grains, the exact contribution of different chemical mechanisms is still poorly constrained. Traditionally it is deemed that the CO + OH reaction occurring on top of ices is the main reaction path for its formation. Recent investigations showed that in reality the reaction presents a more complex mechanism, requiring an additional H-abstraction step. Building on our previous works, we carried out a detailed investigation of such H abstraction reactions with the hydrogen atom as a reactant for the abstraction reaction. We found an unconventional chemistry for this reaction, markedly depending on the isomeric form of the HOCO radical prior to reaction. The favored reactions are t-HOCO + H -> CO + H$_2$O, c-HOCO + H -> CO$_2$ + H$_2$ and t/c-HOCO + H -> c/t-HCOOH. We estimate bounds for the rate constants of the less favored reaction channels, t-HOCO + H -> CO$_2$ + H and c-HOCO + H -> CO + H$_2$O, to be approximately 10$^{4-6}$ s$^{-1}$. However, these estimates should be interpreted cautiously due to the significant role of quantum tunneling in these reactions and the complex electronic structure of the involved molecules, which complicates their study. Our findings underscore the need for detailed investigation into the chemistry of interstellar CO$_2$ and pave the way for a reevaluation of its primary formation mechanisms in the interstellar medium.

We report on the discovery of rare emergence (31 nights from 360 nights of observations) of narrow absorption features in hydrogen and helium lines in stationary SS433 spectra with velocities ranging from $-650$ to $-1900$ km/s. The components arise independently of the appearance of P-Cygni line profiles which are frequently observed in the SS433 stationary spectra with terminal velocities ranging from $-200$ to $\sim -2500$ km/s. The characteristic rising time of the transient absorptions is about one day and the decay time is about two days. The phenomenology of the absorptions suggests their origin due to hydrodynamic instabilities of wind outflows from a supercritical accretion disk in SS433.

The properties of host galaxies associated with Fast Radio Bursts (FRBs) provide critical information for inferring the progenitors and radiation mechanisms of these bursts. We report on the host galaxy of the repeating FRB 20190520B, a dwarf galaxy at the spectroscopic redshift $z=0.241$ with a stellar mass of $(6.2 \pm 0.8) \times 10^8 \ M_{\odot}$. The emission line ratios suggest that the ionized gas is powered by star formation. The total H$\alpha$-traced star formation rate (SFR) is $0.70 \pm 0.01 \ {M_{\odot} ~ \rm yr^{-1}}$, and the metallicity is $\rm 12+log_{10} ([O/H]) \geq 7.4 \pm 0.1$. The specific star formation rate (sSFR) is $\rm log \ sSFR/yr^{-1} = -9.0 \pm 0.1$, higher than the upper limit of $-9.4$ observed in nearby dwarf galaxies. The dispersion measure contribution from the host galaxy is estimated to be $\rm DM_{host} \approx 950 \pm 220 \ pc \ cm^{-3}$, based on the H$\alpha$ emission. The FRB and the associated persistent radio source are located at the H$\alpha$ emission peak, offset by $\sim 1.4^{\prime\prime}$ (5.5 kpc) in projection from the stellar continuum. At this position, the lower limit of $\rm \log \ sSFR/yr^{-1}$ is $-8.5 \pm 0.1$, more than three times the galaxy's total sSFR. The H$\alpha$ velocity difference between the stellar continuum and the offset gas is $39.6 \pm 0.4$ km s$^{-1}$, which is sufficient to draw conclusions about the nature of the offset.

Recent studies indicate that circumstellar disks exhibit weak turbulence, with their dynamics and evolution being primarily influenced by magnetic winds. However, most numerical studies have focused on planet-disk interactions in turbulent disk models. We aim to explore how wind-driven accretion affects the orbital and eccentricity evolution of a Jovian planet within a magnetized disk. Conversely, we seek to determine in what extent such a planet can modify the accretion behavior and the wind dynamics. We perform high-resolution 3D global non-ideal magneto-hydrodynamic (MHD) simulations of a massive gap-carving planet interacting with a wind-launching disk, using the accelerated code IDEFIX. We consider the influence of the gap shape on planet migration by restarting a "fixed-planet" simulation at three different times, from which the planet evolves freely in the disk. For a strong initial magnetization and a sufficiently deep planet gap, we find that the planet becomes moderately eccentric, and its migration is slow, unsteady and mostly outward. This migration pattern is due to the gap's radial asymmetry which enhances the inner Lindblad torque while reducing the outer Lindblad torque. We show that eccentricity can grow up to 6-8% and is likely driven by a finite-amplitude instability triggered by first-order external Lindblad resonances. These moderate eccentricity values periodically modulate the gap accretion rate and wind mass loss rate, possibly leading to the formation of discrete structures in CO outflows. Slow outward migration and eccentricity growth appear to be common outcomes of planet-disk-wind interactions, which may contribute significantly to both the long orbital periods and the moderate eccentricities of warm jupiters. Additionally, eccentric massive protoplanets embedded in circumstellar disks could play a role in generating structured outflows.

We extend the survey for organics in the southern hemisphere by observing two cores in the Chamaeleon complex using NASA's Deep Space Network 70-m antenna in Canberra, Australia, over the frequency range of 18 to 25 GHz. We surveyed the class 0 protostar Cha-MMS1 and the prestellar core Cha-C2, which represent two stages in the evolution of dense cores. We detect several molecules including HC$_3$N, HC$_5$N, C$_4$H, CCS, C$_3$S, NH$_3$, and c-C$_3$H$_2$. A longer cyanopolyyne, HC$_7$N, is detected with high confidence via spectral stacking analysis. While molecular column densities in the two Chamaeleon cores are typically an order of magnitude lower compared to the cynaopolyyne peak in TMC-1, the molecular abundance ratios are in general agreement with the TMC-1 values. The two exceptions are c-C$_3$H$_2$, which is enhanced by a factor of \about 25 with respect to cyanopolyynes in the Chamaeleon cores, and ammonia, which is enhanced by a factor of ~ 125. The deuterated species c-C$_3$HD is detected in both cores, with a high D/H ratio of ~0.23 in c-C$_3$H$_2$. A rare isotopologue of ammonia, $^{15}$NH$_3$, is also detected in Cha-MMS1 suggesting a high $^{14}$N/$^{15}$N ratio of ~ 690 in ammonia. However, this ratio may be artificially enhanced due to the high optical depth of the $^{14}$NH$_3$ (1,1) line, which increases the effective source size. We use the detections of ammonia, cyanopolyynes, and far-infrared dust continuum to characterize the density and temperature in the Chamaeleon cores and calculate the molecular column densities and their relative ratios. The ring molecule benzonitrile is not detected in either Chamaeleon core. The $3 \sigma$ upper limits for its column density are a factor of 2 higher than the value derived for TMC-1 and the upper limits for its relative abundance with respect to HC$_5$N are a factor of 3 higher than the TMC-1 value.

The isopoly bi-fluid approach assumes an isothermal evolution of the solar wind near the Sun up to the radial distance riso, followed by a polytropic evolution constrained by the observed polytropic indices. This approach provides a more accurate model of the interplanetary properties of the solar wind (u, n, Tp, Te) and their radial evolution (Dakeyo et al. 2022, 2024). In this article, we present an improvement of the isopoly approach by considering a generalized two thermal regime approach, embedding two distinct polytropic evolutions, the "bipoly" modeling. To demonstrate the capability of the approach, the models are fitted to both interplanetary and coronal observations, all classified by wind speed population in the spirit of Maksimovic et al. (2020). The set of observations used as constraints are coronal temperatures inferred from charge-state ratio observations from Solar Orbiter, and interplanetary measurements from Helios and Parker Solar Probe. The relaxation of the isothermal criteria in the near-Sun region permits to significantly improve the fast wind acceleration for low coronal temperature conditions. In summary, the new model matches closely the observational constraints: the coronal temperature and the radial evolution of the wind properties (u, n, Tp, Te) in the interplanetary medium, and this for all the wind speed populations.

The Simons Observatory (SO) Small Aperture Telescopes (SATs) will observe the Cosmic Microwave Background (CMB) temperature and polarization at six frequency bands. Within these bands, the angular response of the telescope (beam) is convolved with the instrument's spectral response (commonly called bandpass) and the signal from the sky, which leads to the band-averaged telescope beam response, which is sampled and digitized. The spectral properties of the band-averaged beam depend on the natural variation of the beam within the band, referred to as beam chromaticity. In this paper, we quantify the impact of the interplay of beam chromaticity and intrinsic frequency scaling from the various components that dominate the polarized sky emission on the tensor-to-scalar ratio, $r$, and foreground parameters. We do so by employing a parametric power-spectrum-based foreground component separation algorithm, namely BBPower, to which we provide beam-convolved time domain simulations performed with the beamconv software while assuming an idealized version of the SO SAT optics. We find a small, $0.02\sigma$, bias on $r$, due to beam chromaticity, which seems to mostly impact the dust spatial parameters, causing a maximum $0.77 \sigma$ bias on the dust $B$-mode spectra amplitude, $A_{d}$, when employing Gaussian foreground simulations. However, we find all parameter biases to be smaller than $1\sigma$ at all times, independently of the foreground model. This includes the case where we introduce additional uncertainty on the bandpass shape, which accounts for approximately half of the total allowed gain uncertainty, as estimated in previous work for the SO SATs.

The explosion of high-precision astrometric data on main-belt asteroids (MBAs) enables new inferences on gravitational and non-gravitational forces present in this region. We estimate the size of MBA motions caused by mutual gravitational encounters with other MBAs that are either omitted from ephemeris models or have uncertain mass estimates. In other words, what is the typical Brownian motion among MBAs that cannot be predicted from the ephemeris, and therefore serves as noise on inferences from the MBAs? We estimate the RMS azimuthal shift $\sigma_\phi$ of this ``Brownian noise'' by numerical estimation of the distribution of impulse sizes and directions among known MBAs, combined with analytical propagation into future positional uncertainties. At current levels of asteroid-mass knowledge, $\sigma_\phi$ rises to $\approx2$ km or $\approx1$ mas over $T=10$ yr, increasing as $T^{3/2},$ large enough to degrade many inferences from Gaia and LSST MBA data. LSST data will, however, improve MBA mass knowledge enough to lower this Brownian uncertainty by $\ge4\times.$ Radial and vertical Brownian noise at $T=10$ yr are factors $\approx7$ and $\approx45,$ respectively, lower than the azimuthal noise, and grow as $T^{3/2}$ and $T^{1/2}.$ For full exploitation of Gaia and LSST MBA data, ephemeris models should include the $\approx1000$ largest asteroids as active bodies with free masses, even if not all are well constrained. This will correctly propagate the uncertainties from these 1000 sources' deflections into desired inferences. The RMS value of deflections from less-massive MBAs is then just $\sigma_\phi\approx60$ m or 30 $\mu$as, small enough to ignore until occultation-based position data become ubiquitously available for MBAs.

Pulsar wind nebulae (PWNe) are important sources for understanding galactic high-energy processes, but it is controversial until now about how high-energy particles in PWNe are accelerated and transported. Lacking radio counterparts of X-ray PWNe (the proposed acceleration sites) introduce difficulties to better understandings in multi wavelengths. Our recent 3, 6, and 16\,cm high-resolution observations of G11.2$-$0.3 PWN with the Australia Telescope Compact Array (ATCA) uniquely show morphological similarity with its X-ray PWN (a torus/jet feature). Spectral indices of the radio torus and jet are around -0.09 and -0.10, respectively. Meanwhile for the jet region, the spectral break between radio and X-ray spectra implies particle acceleration mechanisms other than a diffusive shock acceleration. Polarization results suggest a helical B-field inside the jet, the equipartition B-field strength of which is below 100\,$\mu$G.

In this review we look into the gauge-dependence of scalar-induced gravitational waves (SIGWs) that are second-order tensors produced by first-order scalar-modes. The method includes deriving the background, first- and second-order Einstein field equations without imposing a gauge. We address the gauge-invariant approach and study the source-term of SIGWs in three different gauges, synchronous, Poisson and uniform curvature gauge. We find that numerically computed kernels in all three gauges behave closely with minimal discrepancy. As expected, when going in sub-horizon modes, $x\gg1$, the discrepancy decreases and the behavior matches, pointing to a gauge-invariant observable.

Thermodynamic analyses of dark energy as a relativistic fluid indicates that this intriguing component of the universe mimics a bulk viscous pressure when the parameter of its barotropic equation of state varies with time. Since in cosmology bulk viscosity and creation or destruction of matter are closely linked processes, we propose in this work a brief thermodynamic study of dark energy considering that particles can be created or destroyed in the fluid. We derive new expressions for quantities such as particle density, entropy density etc. that have been shown to be sensitive to this new ingredient. We also obtain new thermodynamic constraints and compare them with those where the number of particles is conserved. In particular, we found that in the presence of a sink, dark energy tends towards the cosmological constant over time regardless of the sign of its chemical potential and without violating the laws of thermodynamics.

We study the possibility of the highest energy neutrino event with 220 PeV energy, detected recently by the KM3NeT experiment to be originating from heavy dark matter (DM) decay. Considering a heavy right handed neutrino (RHN) DM for illustrative purpose, we show that DM mass of 440 PeV, can explain the observed flux. The required DM lifetime to produce the best-fit value of the neutrino flux saturates the existing gamma-ray bounds. Due to the large uncertainty in the flux, it is possible to explain the KM3NeT event from RHN DM decay at $3\sigma$ confidence level (CL) while being in agreement with gamma-ray bounds and non-observation of similar events at IceCube. While we consider a gauged $B-L$ scenario where DM relic can be generated due to other interactions, we also briefly discuss some alternate DM possibilities where the gamma-ray bounds can be alleviated compared to the minimal RHN DM discussed here.

We assess the universal relations among second-order moments of relativistic stars, namely the moment of inertia, tidal deformability, and spin-induced quadrupole moment, via reformulated perturbation equations. After constructing the spherical background configuration by solving two ordinary differential equations as usual, these three moments are obtained by solving four additional ordinary differential equations. They are solved numerically from the stellar center to the surface, and we do not need to derive homogeneous solutions for obtaining the quadrupole moment. This small number of ordinary differential equations to be solved enables us to identify the primary variable for each second-order moment. Investigating the profile of these variables in the star, we speculate that their nonmonotonic behavior, enhanced typically for soft equations of state and/or high compactnesses, introduces the variety to the relations among these second-order moments unless the black-hole limit is approached.

Leptophilic sub-MeV spin-zero dark matter (DM) decays into photons via one-loop processes, a scenario that has been in part overlooked in current literature. In this work, we provide updated and comprehensive upper limits on scalar, pseudo-scalar, and axion-like DM-electron couplings based on the latest NPIPE cosmic microwave background data from Planck. Our bounds on the couplings are not only competitive with astrophysical and terrestrial experiments, but outperform them in certain regions of parameter space. Notably, we present the most stringent limits to date on scalar DM with masses around a few keV and pseudo-scalar DM with masses between 100 eV and a few keV. Additionally, we explore, for the first time, the impact of implementing a cosmology-consistent treatment of energy deposition into the cosmic medium.

Light rings (LRs) - closed circular orbits of null geodesics - are key features of both black holes and horizonless ultracompact objects. While unstable LRs are relevant for the observation of black hole images, stable LRs have been suspected to trigger instabilities, namely in exotic compact objects that could mimic black holes. The underlying mechanism behind this instability remains poorly understood, but a key missing piece is how the backreaction of a perturbation around the stable LR modifies the surrounding spacetime. In this work, some progress in this direction is provided by examining a conceptually simple, yet instructive, toy model: continuum-shell stars, supported solely by tangential pressures. Using both analytical and numerical methods, we show how perturbations around the stable LR deepen the geodesic potential and shifts its location inward, potentially amplifying any instability associated with the LR. By then extending the analysis to more general stars with nonzero radial pressure, we find that the same phenomenon can be expected to persist under reasonable assumptions.

Understanding the properties, especially the magnetohydrodynamic (MHD) invariants, of coronal mass ejections (CMEs) measured in-situ is key to bridging the CME properties from the Sun to interplanetary space. In order to investigate CMEs from the in-situ measurements that provide a one-dimensional (1-D) cut of the CME parameters over the spacecraft trajectory, various magnetic flux rope (MFR) models have been developed, among which the models with a circular cross-section are the most popular and widely used. CMEs are found to be deformed during their propagation in interplanetary space, in which the cross-section may be flattened in the direction of propagation, i.e., to develop an elliptical or even pancake-like shape. We use numerical MHD simulations in 2.5-D to investigate the influence of the CME deformation on the in-situ fitting using two linear force-free MFR models with a circular cross-section, and we focus on the axial and poloidal magnetic fluxes, which are conserved in the ideal MHD frame and simulations. We quantitatively compare the fitted axial and poloidal fluxes with those in simulations. We find that both models underestimate the axial flux compared to that in simulations, and such underestimation depends on the CME deformation. However, the fitting of the poloidal flux is independent of the deformation. We discuss the reasons for the axial flux underestimation and the implication of the CME deformation for the CME in-situ fitting.

We investigate the classical and quantum dynamics of $f(R)$ gravity's rainbow in the presence of a perfect fluid, employing Schutz's formalism to establish a well-defined notion of time. In the classical regime, we derive and solve the equations of motion, obtaining both analytical and numerical solutions. Through canonical quantisation, we formulate the Schrödinger-Wheeler-DeWitt (SWD) equation for the quantum model. By solving its eigenfunctions, we construct the wave function of the Universe and obtain analytical solutions in scenarios dominated by stiff matter. Our results highlight the impact of rainbow gravity on quantum evolution, particularly in modifying the structure of the wave function and shaping the transition from the quantum to the classical regime.

This study investigates the properties of the thin accretion disk around a black hole in quantum fluctuation modified gravity (QFMGBH) using the Novikov-Thorne model. We restrict the parameter \( \alpha \) characterizing the quantum fluctuation of metric and compute the ISCO radius, examining its effect on the energy flux, the radiation temperature, the luminosity spectrum, the energy efficiency, and the shadow size for various values of the parameter \( \omega \) characterizing the matter around the black hole. Our findings show that the ISCO radius increases with \( \alpha \) for \( \omega = 1/3, 0, -2/3 \), but decreases for \( \omega = -4/3 \). As \( \alpha \) increases, the energy flux and the temperature show distinct trends, with changes in the luminosity spectrum and the efficiency. The shadow size increases for \( \omega = 1/3, 0, -2/3 \), and decreases for \( \omega = -4/3 \). The accretion disk around the QFMGBH is smaller, hotter, and brighter than that surrounded the Schwarzschild BH in GR, suggesting that the observable differences can potentially distinguish quantum fluctuation modified gravity from standard GR.

Collective neutrino flavor oscillations are of primary importance in understanding the dynamic evolution of core-collapse supernovae and subsequent terrestrial detection, but also among the most challenging aspects of numerical simulations. This situation is complicated by the quantum many-body nature of the problem due to neutrino-neutrino interactions which demands a quantum treatment. An additional complication is the presence of three flavors, which often is approximated by the electron flavor and a heavy lepton flavor. In this work, we provide both qubit and qutrit encodings for all three flavors, and develop optimized quantum circuits for the time evolution and analyze the Trotter error. We conclude our study with a hardware experiment of a system of two neutrinos with superconducting hardware: the IBM Torino device for qubits and AQT device for qutrits. We find that error mitigation greatly helps in obtaining a signal consistent with simulations. While hardware results are comparable at this stage, we expect the qutrit setup to be more convenient for large-scale simulations since it does not suffer from probability leakage into nonphsycial qubit space, unlike the qubit setup.

We present an exact black hole solution surrounded by massive vector fields predicted by Kaluza-Klein (KK) gravity. KK gravity in four dimensions (4D) is of particular interest, as it predicts a tower of particle states, including gravitons with spin-0 and spin-1 components, in addition to the massless spin-2 gravitons of general relativity. The extra degrees of freedom in the gravitational sector modify the law of gravity, allowing the theory to explain the effects attributed to dark matter in the universe. In this paper, we construct a black hole solution surrounded by massive spin-1 gravitons within KK theory. In addition to the influence of the massive vector fields, we incorporate an interaction term between the black hole and the massive vector field. The black hole solution is affected by the mass of the spin-1 graviton and an additional parameter that encodes corrections to Newton's constant, as well as the coupling between the massive vector field and the black hole mass. We find that the coupling between the massive vector field and the black hole mimics the effect of an electric charge. To this end, we investigate the accretion disk, quasinormal modes (QNMs), and the stability of the black hole spacetime. Finally, we use Event Horizon Telescope (EHT) observations of Sgr A* to constrain the black hole parameters.

The gravitational wave (GW) spectrum from the first-order phase transition can be characterized by a few phenomenological parameters but with high degeneracies in model/data distinguishments. In this paper, we look into the high-frequency power law of the GW spectrum with preliminary numerical simulations for both quantum and semi-classical pictures of vacuum decay. We first reveal an anti-correlation of the high-frequency power law to a certain power of the ratio between the wall thickness and bubble radius at the onset of bubble collisions, which can be further approximated analytically by some other phenomenological model characteristics to break the model degeneracy.

Superradiance can cause the axion cloud around a rotating black hole to reach extremely high densities, and the decay of these axions can produce a powerful laser. The electric field of these lasers is strong enough that the Schwinger effect may become significant, resulting in the production of an electron-positron plasma. We explore the dynamics between axion lasers and this electron-positron plasma. While there are several mechanisms by which the inclusion of a plasma can impact the laser's behavior, the most significant of these mechanisms is that the electron-positron plasma imparts an effective mass on the photon. As the plasma frequency increases, axion decay becomes energetically unfavorable, up to the point where the axion no longer decays into photons, shutting off the laser. We find that the impact of the electron-positron plasma on the dynamics of the system depend heavily on the parameters, specifically the axion mass $m_\phi$ and the superradiant coupling $\alpha$, and that we may divide parameter space into three regimes: the unenhanced, enhanced, and unstable regimes. In the unenhanced and enhanced regime, the system will eventually settle into an equilibrium state, emitting a laser of constant luminosity while the number of axions remains constant. In the unenhanced regime, this equilibrium state can be calculated while neglecting the effects of Schwinger production; in the enhanced regime, the equilibrium luminosity is slightly larger than what it would be without Schwinger production. In the unstable regime, the electron-positron plasma suppresses axion decay to the point where the system is never able to reach equilibrium; instead, the axions continue to grow superradiantly. In all three cases, the production of superradiant axions will eventually cause the black hole to spin down to the point where superradiance ceases.

Recent studies of QFT in cosmological spacetime indicate that the speeding up of the present universe may not just be associated with a rigid cosmological term but with a running one that evolves with the expansion rate: $\Lambda=\Lambda(H)$. This running is inherited from the cosmic evolution of the vacuum energy density (VED), $\rho_{\rm vac}$, which is sensitive to quantum effects in curved spacetime that ultimately trigger that running. The VED is a function of the Hubble rate and its time derivatives: $\rho_{\rm vac}=\rho_{\rm vac}(H, \dot{H},\ddot{H},...)$. Two nearby points of the cosmic evolution during the FLRW epoch are smoothly related as $\delta\rho_{\rm vac}\sim {\cal O}(H^2)$. In the very early universe, in contrast, the higher powers of the Hubble rate take over and bring about a period of fast inflation. They originate from quantum effects on the effective action of vacuum, which we compute. Herein we focus on the lowest possible power for inflation to occur: $H^4$. During the inflationary phase, $H$ remains approximately constant and very large. Subsequently, the universe enters the usual FLRW radiation epoch. This new mechanism (`RVM-inflation') is not based on any supplementary `inflaton' field, it is fueled by pure QFT effects on the dynamical background and is different from Starobinsky's inflation, in which $H$ is never constant.

We investigate the reheating process in an axion inflation model where the inflaton couples to non-Abelian gauge fields via the Chern-Simons coupling. The Chern-Simons coupling leads to the efficient production of gauge fields via a tachyonic instability during inflation, whose implications have been actively studied in the literatures. Moreover, it has been recently pointed out that the produced gauge fields can be even thermalized during inflation, leading to warm inflation. Apparently, these findings seem to imply that the reheating is completed immediately after inflation because the tachyonic instability or the thermal friction induced by the Chern-Simons coupling cause the inflaton condensate to decay rapidly. Contrary to this naive expectation, however, we show that, in most of the parameter space, either the inflaton condensate, the inflaton particles, or the glueballs once dominate the Universe and their perturbative decay completes the reheating.

Electrical resistivity plays a fundamental role in many astrophysical systems, influencing the evolution of the magnetic field and energy dissipation processes. During the coalescence of two neutron stars, resistive effects can significantly affect the dynamics and gravitational and electromagnetic signatures associated with these events, such as gamma-ray bursts and kilonova emission. Here, we present our new code, named MIR. Developed within the EinsteinToolkit framework, MIR solves the general relativistic magnetohydrodynamic equations in 3D Cartesian coordinates and on a dynamical spacetime using the 3+1 Eulerian formalism, in both the ideal and resistive regimes, filling a crucial gap in the toolkit's capabilities.

We discuss a class of solutions of multidimensional gravity which are formally related to black-hole solutions but can observationally look like compact stars whose surface reflects back all particles or signals getting there. Some particular examples of such solutions are presented and studied, including those with a magnetic field in Maxwell or nonlinear electrodynamics (NED) in five dimensions. For NED as a possible source for magnetic mirror stars, we formulate a methodology of solving the 5D Einstein-NED equations and point out the conditions under which there always exist mirror star solutions. We also note that some of the Einstein-Maxwell solutions under consideration are discussed in the literature and called ``topological stars'' due to the circular topology of the fifth dimension.