Papers for Tuesday, Jan 28 2025

This paper is part of a series investigating the observational appearance of planets accreting from their nascent protoplanetary disk (PPD). We evaluate the differences between gas temperature distributions determined in our radiation hydrodynamical (RHD) simulations and those recalculated via post-processing with a Monte Carlo (MC) radiative transport (RT) scheme. Our MCRT simulations were performed for global PPD models, each composed of a local 3D high-resolution RHD model embedded in an axisymmetric global disk simulation. We report the level of agreement between the two approaches and point out several caveats that prevent a perfect match between the temperature distributions with our respective methods of choice. Overall, the level of agreement is high, with a typical discrepancy between the RHD and MCRT temperatures of the high-resolution region of only about 10 percent. The largest differences were found close to the disk photosphere, at the transition layer between optically dense and thin regions, as well as in the far-out regions of the PPD, occasionally exceeding values of 40 percent. We identify several reasons for these discrepancies, which are mostly related to general features of typical radiative transfer solvers used in hydrodynamical simulations (angle- and frequency-averaging and ignored scattering) and MCRT methods (ignored internal energy advection and compression and expansion work). This provides a clear pathway to reduce systematic temperature inaccuracies in future works. Based on MCRT simulations, we finally determined the expected error in flux estimates, both for the entire PPD and for planets accreting gas from their ambient disk, independently of the amount of gas piling up in the Hill sphere and the used model resolution.

The Sun can efficiently capture leptophilic dark matter that scatters with free electrons. If this dark matter subsequently annihilates into leptonic states, it can produce a detectable neutrino flux. Using 10 years of Super-Kamiokande observations, we set constraints on the dark-matter/electron scattering cross-section that exceed terrestrial direct detection searches by more than an order of magnitude for dark matter masses below 100 GeV, and reach cross-sections as low as $\sim$4$\times$10$^{-41}$cm$^{-2}$.

Thermal dark matter models generically include the prompt creation of highly-concentrated dark matter cusps in the early Universe. Recent studies find that these cusps can survive to the present day, as long as they do not fall into extremely dense regions of baryonic structure. In this work, we build models of dark matter annihilation within the prompt cusps that reside in galaxy clusters, showing that they dominate the total $\gamma$-ray annihilation signal. Using 15 years of Fermi-LAT data, we find no evidence for a $\gamma$-ray excess from these sources, and set strong constraints on annihilating dark matter. These constraints generically rule out the thermal annihilation cross-section to the $b\bar{b}$ channel for dark matter masses below $\sim$200~GeV.

We focus on two key diagnostics of stellar physics in red giant branch (RGB) stars: the first dredge-up (FDU) of nuclear processed material and the location of the red giant branch bump (RGBB). We compare asteroseismic and spectroscopic APOKASC-3 data with theoretical MESA models. Our FDU predictions have similar mass and metallicity trends to the data, but the observed magnitude of the change in $[{\rm C}/{\rm N}]$ in data is smaller than theoretical predictions by $0.1615 \pm 0.0760 \,({\rm obs}) \pm 0.0108 \,({\rm sys}) \,{\rm dex}$. These results are insensitive to the input physics, but they are at a level consistent with systematic uncertainties in the abundance measurements. When we include observed trends in birth $[{\rm C}/{\rm Fe}]$ and $[{\rm N}/{\rm Fe}]$ in our models, it modestly stretches the metallicity dependent difference relative to the data. We find a well-defined empirical RGBB locus: $\log g = 2.6604 - 0.1832 (M/{\rm M}_\odot-1) + 0.2824 \,[{\rm Fe}/{\rm H}]$. Our model RGBB loci have mass and composition trends that mirror the data, but we find that the observed RGBB is $0.1509 \pm 0.0017 \,({\rm obs}) \pm 0.0182 \,({\rm sys})$ magnitudes higher than predicted across the board, similar to prior literature results. We find that envelope overshooting, proposed solution to reconcile theory with data, increases ${\rm Li}$ destruction during the FDU at higher metallicities, creating tension with depletion observed in GALAH data. We propose ${\rm Li}$ in the FDU as a sensitive test of the RGBB and FDU, and discuss other potential solutions.

We present a 5 month NICER X-ray monitoring campaign for two low luminosity active galactic nuclei (LLAGNs) -- NGC 4594 and IC 1459 -- with complementary Swift and NuSTAR observations. Utilizing an absorbed power law and thermal source model combined with NICER's SCORPEON background model, we demonstrate the effectiveness of joint source/background modeling for constraining emission from faint, background-dominated targets. Both sources are dominated by nuclear power law emission with photon indices $\Gamma \sim 1.5 - 2$, with NGC 4594 being slightly harder than IC 1459. The thermal contribution in both sources is fainter, but constant, with $kT \sim 0.5$ keV ($\sim 5 \times 10^6$ K). The power law flux and $\Gamma$ are strongly anti-correlated in both sources, as has been seen for other LLAGNs with radiatively inefficient accretion flows. NGC 4594 is the brighter source and exhibits significant aperiodic variability. Its variability timescale with an upper limit of $5 - 7$ days indicates emission originating from $< 100 R_{g}$, at the scale of the inner accretion flow. A spectral break found at $\sim 6$ keV, while tentative, could arise from synchrotron/inverse compton emission. This high-cadence LLAGN X-ray monitoring campaign underlines the importance of multi-wavelength variability studies for a sample of LLAGNs to truly understand their accretion and outflow physics.

Local-type primordial non-Gaussianity (PNG), predicted by many non-minimal models of inflation, creates a scale-dependent contribution to the power spectrum of large-scale structure (LSS) tracers. Its amplitude is characterized by the product $b_\phi f_{\rm NL}^{\rm loc}$, where $b_\phi$ is an astrophysical parameter dependent on the properties of the tracer. However, $b_\phi$ exhibits significant secondary dependence on halo concentration and other astrophysical properties, which may bias and weaken the constraints on $f_{\rm NL}^{\rm loc}$. In this work, we demonstrate that incorporating knowledge of the relation between Lagrangian bias parameters and $b_\phi$ can significantly enhance PNG constraints. We employ the Hybrid Effective Field Theory (HEFT) approach at the field-level and a linear regression model to seek a connection between the bias parameters and $b_{\phi}$ for halo and galaxy samples, constructed using the \textsc{AbacusSummit} simulation suite and mimicking the luminous red galaxies (LRGs) and quasi-stellar objects (QSOs) of the Dark Energy Spectroscopic Instrument (DESI) survey. For the fixed-mass halo samples, our full bias model reduces the uncertainty by more than 70\%, with most of that improvement coming from $b_\nabla$, which we find to be an excellent proxy for concentration. For the galaxy samples, our model reduces the uncertainty on $b_\phi$ by 80\% for all tracers. By adopting Lagrangian-bias informed priors on the parameter $b_\phi$, future analyses can thus constrain $f_{\rm NL}^{\rm loc}$ with less bias and smaller errors.

We present a study of new 7.7-11.3 $\mu$m data obtained with the James Webb Space Telescope Mid-InfraRed Instrument in the starburst galaxy M 82. In particular, we focus on the dependency of the integrated CO(1-0) line intensity on the MIRI-F770W and MIRI-F1130W filter intensities to investigate the correlation between CO content and the 7.7 and 11.3 $\mu$m features from polycyclic aromatic hydrocarbons (PAH) in M 82's outflows. To perform our analysis, we identify CO clouds using archival $^{12}$CO($J$=1-0) NOEMA moment 0 map within 2 kpc from the center of M 82, with sizes ranging between $\sim$21 and 270 pc; then, we compute the CO-to-PAH relations for the 306 validated CO clouds. On average, the power-law slopes for the two relations in M 82 are lower than what is seen in local main-sequence spirals. In addition, there is a moderate correlation between $I_{\rm CO(1-0)}$-$I_{\rm 7.7\mu m} /I_{\rm 11.3\mu m}$ for most of the CO cloud groups analyzed in this this http URL results suggest that the extreme conditions in M 82 translate into CO not tracing the full budget of molecular gas in smaller clouds, perhaps as a consequence of photoionization and/or emission suppression of CO molecules due to hard radiation fields from the central starburst.

The linear stability of waves driven by ion beams produced during solar flare energy release are explored to assess their role in driving abundance enhancements in minority species such as $^3$He and in controlling, through pitch-angle scattering, proton/alpha confinement during energy release. The Arbitrary Linear Plasma Solver (ALPS) is used to solve the linear dispersion relation for a population of energetic, reconnection-accelerated protons streaming into a cold background plasma. We assume equal densities of the two populations, using an anisotropic ($T_\parallel/T_\perp = 10$), one-sided kappa distribution for the energetic streaming population and a cold Maxwellian for the background. We find two unstable modes with parallel propagation: a right-handed wave with a frequency of the order of the proton cyclotron frequency ($\Omega_{cp}$) and a left-handed, lower frequency mode. We also find highly oblique modes with frequencies below $\Omega_{cp}$ that are unstable for higher beam energies. Through resonant interactions, all three modes will contribute to the scattering of the high-energy protons, thereby limiting their transport out of the flare-acceleration region. The higher-frequency oblique mode, which can be characterized as a kinetic Alfvén wave, will preferentially heat $^3$He, making it a good candidate for the driver of the abundance enhancements commonly observed for this species in impulsive events.

In a recent study, Dannen et al. surveyed a large parameter space to study the transition from efficient to inefficient line driving. They found that when the line force significantly weakens due to ionization, the winds are variable, with a characteristic frequency comparable to the Lamb cut-off frequency of a stratified atmosphere, {\omega}c. In this work, we present a set of simulations and perturbation analyses that elucidate the variability source and characteristics. We found that the line force adds wave energy and amplifies perturbations with frequencies near {\omega}c. This selective amplification results from the coupling between the natural tendency of velocity perturbations to grow in a stratified atmosphere and the dependence of the line force on the velocity gradient, per the Castor-Abbott-Klein line-driven wind theory. We also found that the variability stems from self-excitation that occurs in the exponential atmosphere due to the non-linearity introduced by the absolute value of the velocity gradient in the line force prescription. We conclude that self-consistently calculating ionization is insufficient for modeling the dynamics in the subsonic atmosphere. Instead future wind models should relax the Sobolev approximation, or model the radiative transfer to capture the dynamics and instabilities at the base of the wind.

This paper rigorously examines the potential of the $f(T, \mathcal{T})$ theory as a promising framework for understanding the dark sector of the universe, particularly in relation to cosmic acceleration. The $f(T, \mathcal{T})$ theory extends gravitational dynamics by incorporating both the torsion scalar $T$ and the trace of the energy-momentum tensor $\mathcal{T}$. Further, we explore the functional form $f(T, \mathcal{T}) = T + \beta \mathcal{T}$, where $\beta$ is a free parameter that modulates the matter's influence on spacetime evolution. To evaluate this model, we employ an effective EoS parameter dependent on redshift $z$, to solve the field equations and analyze the evolution of the Hubble parameter $H(z)$. Using a joint dataset ($H(z)+Pantheon^+$) and the Markov Chain Monte Carlo (MCMC) method with Bayesian analysis, we obtain the best-fit parameter values: $H_0 = 68.04 \pm 0.64$, $\beta = 0.14 \pm 0.17$, and $\gamma = 0.96^{+0.38}_{-0.69}$, which align well with current observational data. Our findings indicate a deceleration parameter of $q_0 = -0.51$, supporting a present-day accelerated expansion phase, with a transition redshift $z_t = 0.57$ marking the universe's shift from deceleration to acceleration. Moreover, we confirm a positive cosmic fluid energy density, reinforcing stability, and find an EoS parameter value of $\omega_0 = -0.76$, consistent with quintessence-driven acceleration. These results underscore the viability of $f(T, \mathcal{T})$ as a robust framework for addressing the accelerating universe and dark energy dynamics, paving the way for future investigations into its cosmological implications.

This work deals with the presence of the cuscuton term in the otherwise standard dark energy evolution under the usual FLRW background. We disclose a first-order framework similar to the Hamilton-Jacobi formalism, which helps us to solve the equations of motion and find analytical solutions. We explore several possibilities, concentrating mainly on how the cuscuton-like contribution works to modify cosmic evolution. Some results are of current interest since they describe scenarios capable of changing the evolution, adding or excluding possible distinct phases during the Universe's expansion history. Additionally, we present interesting constraints on the cuscuton-like contribution for the dark energy evolution using a set of homogeneous geometrical observational probes. Finally, based on the Akaike Information Criterion (AIC), we perform a statistical comparison of the cuscuton-like model with $\Lambda$CDM, and find strong support for our model.

The Hubble constant ($H_0$) tension is one of the biggest challenges in modern cosmology. This consists of the discrepancy, at around $5\sigma$, between the local value of $H_0$ measured through Supernovae Ia (SNe Ia) constrained with the Cepheids and the value inferred from the observations of Cosmic Microwave Background (CMB) by Planck data. According to the most appealing cosmological models, such as the flat $\Lambda$CDM, the $H_0$ should not vary according to the measurement method or the redshift $z$ of the probe used for estimating it. Thus, many ideas have been proposed in the literature to face this tension. In the current work, we summarize the results obtained with the binned analysis of SNe Ia, showing a decreasing trend for $H_0$ with $z$ with an evolutionary coefficient $\eta \sim 0.01$, and we further discuss the impact of high-$z$ probes such as Gamma-ray Bursts (GRBs) and quasars (QSOs) that allow reaching constraints on the cosmological parameters that will extend the Hubble diagram to high-$z$ values.

The cool temperatures of M dwarf atmospheres enable complex molecular chemistry, making robust characterization of M dwarf compositions a long-standing challenge. Recent modifications to spectral synthesis pipelines have enabled more accurate modeling of M dwarf atmospheres, but these methods are too slow for characterizing more than a handful of stars at a time. Data-driven methods such as The Cannon are viable alternatives, and can harness the information content of many M dwarfs from large spectroscopic surveys. Here, we train The Cannon on M dwarfs with FGK binary companions from the Sloan Digital Sky Survey-V/Milky Way Mapper (SDSS-V/MWM), with spectra from the Apache Point Observatory Galactic Evolution Experiment (APOGEE). The FGK-M pairs are assumed to be chemically homogeneous and span $-$0.56 $<$ [Fe/H] $<$ 0.31 dex. The resulting model is capable of inferring M dwarf $T_{\textrm{eff}}$ and elemental abundances for Fe, Mg, Al, Si, C, N, O, Ca, Ti, Cr, and Ni with median uncertainties of 13 K and 0.018$-$0.029 dex, respectively. We test the model by verifying that it reproduces reported abundance values of M dwarfs in open clusters and benchmark M dwarf datasets, as well as expected metallicity trends from stellar evolution. We apply the model to 16,590 M dwarfs in SDSS-V/MWM and provide their detailed abundances in our accompanying catalog.

The Stability of Planetary Orbital Configurations Klassifier (SPOCK) package collects machine learning models for predicting the stability and collisional evolution of compact planetary systems. In this paper we explore improvements to SPOCK's binary stability classifier (FeatureClassifier), which predicts orbital stability by collecting data over a short N-body integration of a system. We find that by using a system-specific timescale (rather than a fixed $10^4$ orbits) for the integration, and by using this timescale as an additional feature, we modestly improve the model's AUC metric from 0.943 to 0.950 (AUC=1 for a perfect model). We additionally discovered that $\approx 10\%$ of N-body integrations in SPOCK's original training dataset were duplicated by accident, and that $<1\%$ were misclassified as stable when they in fact led to ejections. We provide a cleaned dataset of 100,000+ unique integrations, release a newly trained stability classification model, and make minor updates to the API.

The Wide Field Survey Telescope (WFST) is a dedicated photometric surveying facility built jointly by the University of Science and Technology of China (USTC) and the Purple Mountain Observatory (PMO). Since many of its scientific objectives rely on near-real-time data for effective analysis, prompt processing of WFST images is of great significance. To meet this need, we adapted the Rubin Observatory Legacy Survey of Space and Time (LSST) science pipelines to handle the data collected by WFST. This paper presents the complete data processing workflow, from ingestion of raw images to the distribution of alerts, and details the primary data products generated by our pipeline. Researchers using data processed by this pipeline can refer to this document to fully understand the data processing procedures.

We constrain the spin of the black hole (BH) candidate MAXI J1727-203 using Insight-HXMT data. Due to limited HXMT observations covering only part of the outburst, NICER data were used to analyze the full outburst's state transitions, we identified two of three HXMT observations in the high soft state and applied the continuum-fitting method to measure the spin. Based on previous estimates and continuum spectral fittings, we explored the parameter space and found that the best-fitting values were $(D, i, M) \approx (6\ \text{kpc}, 30^\circ, 12 M_{\odot})$. We also tested the variation of these parameters using Monte Carlo simulations, sampling over 3000 sets within the parameter ranges: $5.9 \text{kpc}< D<7 \text{kpc}$, $24^\circ<i< 35^\circ$, and $10 M_{\odot}<M<14 M_{\odot}$, yielding a spin measurement of $a=0.34_{-0.19}^{+0.15}$ (1$\sigma$). In addition, we analyzed NuSTAR data in low hard state and found a good fit with the {\tt tbabs*(diskbb+powerlaw)} model, with no significant iron line features observed in the residuals, then the previous reflection model results suggesting an extremely high spin would over-estimate the BH spin.

Studying the interstellar medium in nearby starbursts is essential for understanding the physical mechanisms driving these objects, thought to resemble young star-forming galaxies. This study aims to analyze the physical properties of the first spatially-resolved multi-wavelength SED of an extragalactic source, spanning six decades in frequency (from near-UV to cm wavelengths) at an angular resolution of 3$^{\prime\prime}$ (51 pc at the distance of NGC,253). We focus on the central molecular zone (CMZ) of NGC,253, which contains giant molecular clouds (GMCs) responsible for half of the galaxy's star formation. We use archival data, spanning optical to centimeter wavelengths, to compute SEDs with the GalaPy and CIGALE codes for validation, and analyze stellar optical spectra with the \textsc{starlight} code. Our results show significant differences between central and external GMCs in terms of stellar and dust masses, star formation rates (SFRs), and bolometric luminosities. We identify the best SFR tracers as radio continuum bands at 33 GHz, radio recombination lines, and the total infrared luminosity (L$_{\rm IR}$; 8-1000$\mu$m), as well as 60$\mu$m IR emission. BPT and WHAN diagrams indicate shock signatures in NGC~253's nuclear region, associating it with AGN/star-forming hybrids, though the AGN fraction is negligible ($\leq$7.5%). Our findings show significant heterogeneity in the CMZ, with central GMCs exhibiting higher densities, SFRs, and dust masses compared to external GMCs. We confirm that certain centimeter photometric bands can reliably estimate global SFR at GMC scales.

This study presents the findings of using the Square Kilometre Array (SKA) telescope to measure redshift drift via the HI 21cm signal, employing semi-annual observational interval within redshift around z $\sim$ 1 with main goal is to directly gauge the universe's expansion acceleration rate with millimeter-per-second (mm/s) precision. The SKA can detect over a billion HI 21cm emissions from individual galaxies to the redshift z $\sim$ 2 and thousands of absorption lines from Damped Lyman-alpha (DLA) systems against bright quasars to the redshift z $\sim$ 13, with the sensitivity limit of 100 mJy. By utilizing SKA's high spectral resolution settings (0.001, 0.002, 0.005, 0.01 Hz) to detect redshift drift, particularly focusing on the 0.001 and 0.002 Hz configuration, one aims to achieve the necessary mm/s in precision measurement by the 0.5-year observation period. The velocity drift rate, crucially determined by the two operational regimes within 0.01 to 0.21 mm/s and 0.031 to 0.17 mm/s, respectively, exceeds the theoretical accuracy limit of 1.28 mm/s. The analysis thoroughly restricts cosmological parameters related to dark energy using the Sandage-Loeb (SL) signal from the HI 21cm emission and absorption lines. It estimates $\rm H_0$ of about 70 km/s/Mpc, $\rm \Omega_m$ near 0.3, with w close to -1, $\rm w_0$ around -1, and $\rm w_a$ approaching -0.1. These results strongly endorse the SL effect as an effective method for confirming cosmic acceleration and exploring the dark sector in real-time cosmology with the SKA.

We demonstrate that generative deep learning can translate galaxy observations across ultraviolet, visible, and infrared photometric bands. Leveraging mock observations from the Illustris simulations, we develop and validate a supervised image-to-image model capable of performing both band interpolation and extrapolation. The resulting trained models exhibit high fidelity in generating outputs, as verified by both general image comparison metrics (MAE, SSIM, PSNR) and specialized astronomical metrics (GINI coefficient, M20). Moreover, we show that our model can be used to predict real-world observations, using data from the DECaLS survey as a case study. These findings highlight the potential of generative learning to augment astronomical datasets, enabling efficient exploration of multi-band information in regions where observations are incomplete. This work opens new pathways for optimizing mission planning, guiding high-resolution follow-ups, and enhancing our understanding of galaxy morphology and evolution.

PG1553+113 has drawn significant attention for its quasi-periodic oscillation (QPO) in gamma-ray variability, though the origin of its variability remains uncertain. In this study, we propose a physical mechanism to explain the observed gamma-ray variability within the framework of a supermassive black hole binary (SMBHB) system, supported by a newly identified component hidden in the light curve. A detailed analysis for its about 16-year light curve obtained from Fermi-LAT observations is performed by Gaussian process (GP). As anticipated, the QPO of 2.1 years is effectively captured by the stochastically-driven damped simple harmonic oscillator (SHO) kernel within the under-damped regime, and the overall stochastic nature of the variability is described by the damped random walk (DRW) kernel albeit with an unconstrained damping timescale. Additionally, our results reveal a previously unrecognized component in active galactic nuclei variability, characterized by the Matérn-3/2 kernel, which is typically associated with systems undergoing abrupt energy release. These findings can be consistently interpreted within the SMBHB framework. The QPO of about 2.1 years is likely attributed to the orbital motion in a SMBHB system. The Matérn-3/2 component is interpreted as resulting from magnetic reconnection events triggered by gravitational perturbations of the magnetic field within the jet, occurring as one black hole approaches the other. Meanwhile, in this case, the damping timescale of the common DRW kernel remains unconstrained due to the influence of new perturbations within the system.

This work introduces a systematic method for identifying analytical and semi-analytical solutions of force-free magnetic fields with plane-parallel and axial symmetry. The method of separation of variables is used, allowing the transformation of the non-linear partial differential equation, corresponding to force-free magnetic fields, to a system of decoupled ordinary differential equations, which nevertheless, are in general non-linear. It is then shown that such solutions are feasible for configurations where the electric current has a logarithmic dependence to the magnetic field flux. The properties of the magnetic fields are studied for a variety of physical parameters, through solution of the systems of the ordinary differential equations for various values of the parameters. It is demonstrated that this new logarithmic family of solutions has properties that are highly distinct from the known linear and non-linear equations, as it allows for bounded solutions of magnetic fields, for periodic solutions and for solutions that extend to infinity. Possible applications to astrophysical fields and plasmas are discussed as well as their use in numerical studies, and the overall enrichment of our understanding of force-free configurations.

The cosmic distance duality relates the angular-diameter and luminosity distances and its possible violation may puzzle the standard cosmological model. This appears particularly interesting in view of the recent results found by the DESI Collaboration, suggesting that a dynamical dark energy scenario seems to be favored than a genuine cosmological constant. Accordingly, we take into account possible violations by considering four different parameterizations, namely: a Taylor expansion around $z\simeq 0$, a slightly-departing logarithmic correction, a (1;2) Padé rational series to heal the convergence problem and a Chebyshev polynomial expansion, reducing \emph{de facto} the systematic errors associated with the analysis. We test each of them in a model-independent (-dependent) way, by working out Monte-Carlo Markov chain analyses, employing the Bézier interpolation of the Hubble rate $H(z)$ for the model-independent approach while assuming the flat (non-flat) $\Lambda$CDM and $\omega_0\omega_1$CDM models, motivating the latter paradigm in view of the DESI findings. Subsequently, we explore two analyses, employing observational Hubble data, galaxy clusters from the Sunyaev-Zeldovich effect and type Ia supernovae, investigating the impact of the DESI data catalog, first including then excluding the entire data set. Afterwards, we adopt statistical model selection criteria to assess the statistically favored cosmological model. Our results suggest \emph{no violation} of the cosmic distance duality. While a slight spatial curvature cannot be entirely excluded, the preferred cosmological model remains the flat $\Lambda$CDM background, even when incorporating DESI data. Finally, concerning the Hubble tension, our findings match the Riess estimates, as BAO data points are excluded.

We present a sample of 12 quasar candidates with highly variable soft X-ray emission from the 4th XMM-newton Serendipitous Source Catalog (4XMM-DR13) using random forest. We obtained optical to mid-IR photometric data for the 4XMM-DR13 sources by correlating the sample with the SDSS DR18 photometric database and the AllWISE database. By cross-matching this sample with known spectral catalogs from the SDSS and LAMOST surveys, we obtained a training data set containing stars, galaxies, and quasars. The random forest algorithm was trained to classify the XMM-WISE-SDSS sample. We further filtered the classified quasar candidates with $\it{Gaia}$ proper motion to remove stellar contaminants. Finally, 53,992 quasar candidates have been classified, with 10,210 known quasars matched in SIMBAD. The quasar candidates have systematically lower X-ray fluxes than quasars in the training set, which indicates the classifier is helpful to single out fainter quasars. We constructed a sample of 12 sources from these quasars candidates which changed their soft X-ray fluxes by a factor of 10 over $\sim$ 20 years in the 4XMM-newton survey. Our selected highly variable quasar candidates extend the quasar sample, characterized by extreme soft X-ray variability, to the optically faint end with magnitudes around $r \sim 22$. None of the 12 sources were detected in ROSAT observations. Given the flux limit of ROSAT, the result suggests that quasars exhibiting variations of more than two orders of magnitudes are extremely rare.

We perform three-dimensional supernova simulations with a phenomenological treatment of neutrino flavor conversions. We show that the explosion energy can increase to as high as ~10^51 erg depending on the critical density for the onset of flavor conversions, due to a significant enhancement of the mean energy of electron antineutrinos. Our results confirm previous studies showing such energetic explosions, but for the first time in three-dimensional configurations. In addition, we predict neutrino and gravitational wave (GW) signals from a nearby supernova explosion aided by flavor conversions. We find that the neutrino event number decreases because of the reduced flux of heavy-lepton neutrinos. In order to detect GWs, next-generation GW telescopes such as Cosmic Explorer and Einstein Telescope are needed even if the supernova event is located at the Galactic center. These findings show that the neutrino flavor conversions can significantly change supernova dynamics and highlight the importance of further studies on the quantum kinetic equations to determine the conditions of the conversions and their asymptotic states.

We discuss the implications of the DESI 2024 BAO data on scalar-tensor models of gravity. We consider four representative models: induced gravity (IG, equivalent to Jordan-Brans-Dicke), where we either fix today's value of the effective gravitational constant on cosmological scales to the Newton's constant or allow them to differ, Jordan-Brans-Dicke supplemented with a Galileon term (BDG), and early modified gravity (EMG) with a conformal coupling. In this way it is possible to investigate how different modified gravity models compare with each other when confronted with DESI 2024 BAO data. Compared to previous analyses, for all of these models, the combination of Planck and DESI data favors a larger value of the key parameter of the theory, such as the nonminimal coupling to gravity or the Galileon term, leading also to a larger value of $H_0$, due to the known degeneracy between these parameters. These new results are mainly driven by the first two redshift bins of DESI. In BDG, in which we find the largest value for $H_0$ among the models considered, the combination of Planck and DESI is consistent with CCHP results and reduces the $H_0$ tension with the SH0ES measurement to $1.2\sigma$ (compared to $4.5\sigma$ of $\Lambda$CDM in our Planck + DESI analysis).

With a systematic timing investigation of the persistent X-ray binary pulsar Her X-1 based on a large number of Insight-HXMT observations between 2017 to 2019, we confirm the presence of X-ray millihertz quasi-periodic oscillations (mHz QPOs) at $\sim 0.01$ Hz. By applying wavelet analysis in our data analysis procedures, we firstly identified $\sim 0.005-0.009$ Hz QPOs coexisting with the $\sim 0.01$ Hz QPOs. Wavelet analysis suggests that these QPO features show transient behaviors, frequencies of mHz QPOs evolved in short time scales. There exists a positive relation between QPO centroid frequency (from $\sim 0.005-0.009$ Hz) and the X-ray luminosity, while the 10 mHz QPO frequencies keep nearly constant for different luminosities, which suggests different physical mechanisms for two types of mHz QPOs. The 10 mHz QPOs in both X-ray and UV bands would have the same origin related to the beat frequency where the Alfv$\acute{e}$n radius is close to the corotation radius, and the 5 mHz QPOs may originate from magnetic disk precession.

Understanding the nature of the accretion disk, its interplay with the X-ray corona, and assessing black hole spin demographics are some open challenges in astrophysics. In this work, we examine the predictions of the standard $\alpha$-disk model, origin of the soft X-ray excess, and measure the black hole spin parameter by applying the updated high-density disk reflection model to the XMM-Newton/NuSTAR broadband (0.3$-$78 keV) X-ray spectra of a sample of Type-1 AGN. Our Bayesian analysis confirms that the high-density relativistic reflection model with a broken power-law emissivity profile can simultaneously fit the soft X-ray excess, broad iron K line, and Compton hump for $\sim$70% of the sample, while an additional warm Comptonization model is still required to describe the observed soft X-ray excess for the remaining sources. Our first-ever calculation of the disk-to-corona power transfer fraction reveals that the fraction of power released from the accretion disk into the hot corona can have diverse values, the sample median of which is $0.7_{-0.4}^{+0.2}$. We find that the transferred power from the accretion disk can potentially soften the X-ray spectrum of the hot corona. The median values of the hot coronal temperature and optical depth for the sample are estimated to be $63_{-11}^{+23}$ keV and $0.85_{-0.27}^{+0.12}$, respectively. Finally, through joint XMM-Newton+NuSTAR relativistic reflection spectroscopy, we systematically constrain the black hole spin parameter across the broad range of black hole masses, $\log(M_{\rm BH}/M_{\odot}) \sim 5.5-9.0$, and increase the available spin measurements in the AGN population by $\sim$20%.

The Doppler factor ($\delta$) is a fundamental quality for the relativistic jets from active galactic nuclei (AGNs). It is also fundamental to assessing the number of the entire population of the jetted AGN at high redshift ($z$), and therefore is important for tracing the growth of supermassive black holes (SMBHs) across cosmic time. Here we present the identification of the positive cosmic evolution of the Doppler factor in {\it Fermi}-detected bright $\gamma-$ray blazars. The redshift dependence of the Doppler factor, $\delta\propto(1+z)^{0.8}$, is measured from the observed characteristic energies in the gamma-ray spectra of 141 blazars. Moreover, the analysis of the characteristic timescales derived from modeling the long-term optical light curves of 89 blazars with Gaussian process regression gives $\delta\propto(1+z)^{1.1}$, but with a larger scatter. \textbf{Note that each index is derived from the entire sample, representing an average evolution. Interestingly, the index itself also appears to evolve, with low-luminosity sources showing stronger evolution, as indicated by a larger index.} This detection straightly suggests that relativistic jets from AGNs are much more common at high redshifts than what is previously estimated.

Celestial observations often exhibit inexplicable planetary dependencies when the timing of an observable is projected onto planetary heliocentric positions. This is possible only for incident, non-relativistic streams. Notably, the celebrated dark matter (DM) in the Universe can form streams in our vicinity with speeds of about 240 km/s. Since gravitational impact scales with $1/(\text{velocity})^2$, all solar system objects, including the Sun and the Moon, act as strong gravitational lenses, with their focal planes located within the solar system. Even the Moon can focus penetrating particles toward the Earth at speeds of up to approximately 400 km/s, covering a large portion of the phase space of DM constituents. Consequently, the unexpected planetary dependencies of solar system observables may provide an alternative to Zwicky's tension regarding the overestimated visible cosmic mass. In this work, an overlooked but unexpected planetary dependency of any local observable serves as an analogue to Zwicky's cosmic measurements, particularly if a similar mysterious behavior has been previously noted. Thus, a persistent, unexpected planetary dependency represents a new tension between observation and expectation. The primary argument supporting DM in line with Zwicky's paradigm is this planetary dependency, which, on a local scale, constitutes the novel tension between observation and expectation. In particular, the recurrent planetary dependency of diverse observables mirrors Zwicky's cosmic tension with the overestimated visible mass. No other approach accounts for so many otherwise striking and mysterious observations in physics and medicine.

Millisecond pulsars (MSPs) are recycled pulsars which have been spun-up due to mass accretion during the phase of mass exchange in binaries. Although the interactions with companion stars play important roles on the spin-up process, the global properties of pulsars determined by the equation of state (EoS), such as mass, radius and the moment of inertia, should also play a role. We investigate the spin-up of MSPs in neutron star (NS) and strangeon star (SS) models, both of which have passed the tests by the existence of high-mass pulsars and the tidal deformability of GW170817. Combining the spin-up condition and the transferred angular momentum, we can constrain the accreted mass and the magnetic field strength. The results indicate that it is easier for an SS to form a fully recycled low-mass MPS, whose spin-period is below 10 ms and mass is below about $1.5 M_\odot$, than that for an NS. Finding more MSPs with low mass and and short spin-period could help to put more strict constraints on the EoS of pulsars.

Dust growth is often indirectly inferred observationally in star-forming environments, theoretically predicted to produce mm-sized particles in circumstellar discs, and also presumably witnessed by the predecessors of the terrestrial meteoritic record. For those reasons it is believed that young gas giants under formation in protoplanetary discs with putative circumplanetary discs (CPDs) surrounding them, such as PDS 70c, should be containing mm-sized particles. We model the spectra of a set of CPDs, which we obtained from radiation hydrodynamic simulations at varying Rosseland opacities kappa_R. The kappa_R from the hydrodynamic simulations are matched with consistent opacity sets of ISM-like composition, but grown to larger sizes. Our high kappa_R hydro data nominally corresponds to 10 mum-sized particles, and our low kappa_R-cases correspond to mm-sized particles. We investigate the resulting broad spectral features at first while keeping the overall optical depth in the planetary envelope constant. Dust growth to size distributions dominated by millimeter particles generally results in broad, featureless spectra with black-body like slopes in the far-infrared, while size distributions dominated by small dust develop steeper slopes in the far-infrared and maintain some features stemming from individual minerals. We find that significant dust growth from microns to millimeters can explain the broad features of the PDS 70c data, when upscaling the dust masses from our simulations by x100. Furthermore our results indicate that the spectral range of 30-500 mum is an ideal hunting ground for broadband features arising from the CPD, but that longer wavelengths observed with ALMA can also be used for massive circumplanetary discs

Fast Radio Bursts (FRBs) have emerged as one of the most dynamic areas of research in astronomy and cosmology. Despite increasing number of FRBs have been reported, the exact origin of FRBs remains elusive. Investigating the intrinsic distributions of FRBs could provide valuable insights into their possible origins and enhance the power of FRBs a cosmological probe. In this paper, we propose a hierarchical Bayesian inference approach combining with several viable models to investigate the population information of FRBs released in the CHIME catalog 1. By utilizing this method, we aim to uncover the underlying patterns and characteristics of the FRB population. Taking into account the uncertainties and complex relationships within the data. We find that the distribution of FRBs does not trace the history of star formation, and there is evidence that the FRB population has time delay with respect to the history of star formation.

In November 2020, A0535+26 underwent one of its brightest outbursts, reaching nearly 12 Crab in X-ray flux. Observed by \textit{Insight-HXMT}, \textit{NuSTAR}, \textit{NICER}, and \textit{AstroSat}, this event provided valuable insights into Be/X-ray binaries. The pulse profiles evolved significantly with luminosity, transitioning from pencil-beam to fan-beam geometries. A0535+26, known for its fundamental cyclotron line at $\sim$44 keV, became only the second source to exhibit a negative correlation between cyclotron line energy and flux at high luminosities, with a plateau phase preceding the transition from positive to negative correlation. We report the discovery of a phase-transient low-energy cyclotron line, detected in a narrow phase range ($\sim$16\%) across all seven \textit{NuSTAR} observations during the rising, peak, and declining phases of the outburst. The new line exhibited dramatic variations with pulse phase and luminosity. We explain this behavior using an accretion geometry where the accretion column sweeps across the line of sight.

The origin of the variability in accretion disks of active galactic nuclei (AGN) is still unknown, but its behavior can be characterized by modeling the time series of optical wavelength fluxes coming from the accretion disks with damped random walk (DRW) being the most popular model for this purpose. The DRW is modeled by a characteristic fluctuation amplitude and damping timescale {\tau}, with the latter being potentially related to the mass and accretion rate onto the massive black hole. The estimation of {\tau} is challenging, with commonly used methods such as the maximum likelihood (ML) and the least square error (LSE) resulting in biased estimators. This problem arises most commonly for three reasons: i) the light curve has been observed with additive noise; ii) some cadence scheme; iii) when the autocorrelation parameter is close to one. The latter is called the unit root problem. In order to improve these parameter estimation procedures for estimating {\tau}, we developed a simulation-extrapolation (SIMEX) methodology in the context of time series analysis. We consider both the standard class of autoregressive processes observed at regular intervals and a recently developed class of irregularly autoregressive processes (iAR). The performance of the SIMEX estimation method was evaluated under conditions of near-unit-root behavior and additive noise through extensive Monte Carlo simulations. Monte Carlo experiments confirm that the SIMEX approach outperforms MLE and LSE methods, offering a reduction in estimation bias ranging from 30% to 90%. Real-data applications further validate the methodology, yielding better model fits and lower mean squared errors (MSE). The more accurate estimation of the damping timescale parameter can help to better understand the currently tentative links between DRW and physical parameters in AGN.

Context. Matter that falls onto a protoplanetary disk (PPD) from a protostellar envelope is heated before it cools again. This induces sublimation and subsequent re-adsorption of ices that accumulated during the prestellar phase. Aims. We explore the fate of ices on multiple-sized dust grains in a parcel of infalling matter. Methods. A comprehensive kinetic chemical model using five grain-size bins with different temperatures was applied for an infalling parcel. The parcel was heated to 150 K and then cooled over a total timescale of 20 kyr. Effects on ice loss and re-accumulation by the changed gas density, the maximum temperature, the irradiation intensity, the size-dependent grain temperature trend, and the distribution of the ice mass among the grain-size bins were investigated. Results. A massive selective redistribution of ices exclusively onto the surface of the coldest grain-size bin occurs in all models. The redistribution starts already during the heating stage, where ices that are sublimated from warmer grains re-adsorb onto colder grains before complete sublimation. During the cooling stage, the sublimated molecules re-freeze again onto the coldest grains. In the case of full sublimation, this re-adsorption is delayed and occurs at lower temperatures because a bare grain surface has lower molecular desorption energies in our model. Conclusions. Most protostellar envelope grains enter the PPD ice poor (bare). Ices are carried by a single coldest grain-size bin, here representing 12 % of the total grain surface area. This bare ice-grain dualism can affect the rate of the grain coagulation. The ice components are stratified on the grains according to their sublimation temperatures.

It has been demonstrated that at least 10 percent of the brightest blazars in the fourth Fermi-LAT catalog of $\gamma$-ray sources exhibit repeating patterns of $\gamma$-ray flares. These events may be associated with the presence of a non-uniform sheath surrounding a fast jet spine in some blazars. Theoretical models suggest that such a sheath could facilitate neutrino production in these structured jets. We aim to test the previously reported marginal statistical evidence for a connection between repeating patterns of $\gamma$-ray flares in blazars and high-energy neutrino events that are positionally consistent with these sources. We identified a repeating pattern of flares in the $\gamma$-ray light curve of the blazar PKS 1502+106, which lies within the 50\% uncertainty region of the IC190730A neutrino candidate event. This occurrence is combined with two other high-energy ($\ge 200$ TeV) neutrino events from ICECAT-1 that arrived in both positional and temporal coincidence with two blazars exhibiting ongoing repeating flare patterns. A Monte-Carlo simulation was conducted to evaluate the likelihood of accidental coincidences between the repeating flare patterns and neutrino events, accounting for potential unrecognized systematic uncertainties in the arrival directions of the ICECAT-1 events. Our findings indicate that the probability of a random coincidence in both time and arrival direction for three high-energy neutrino candidates and three blazars with ongoing recurring patterns of $\gamma$-ray flares is $1.56\times 10^{-3}$ ($3.2\sigma$).

The near-Earth asteroid (98943) Torifune, previously designated 2001 CC$_{21}$, is the flyby target of the Hayabusa2 extended mission, nicknamed Hayabusa2$\#$ (SHARP: Small Hazardous Asteroid Reconnaissance Probe). The ground-based telescope observations offer a key science input for the mission's scientific investigation. During 2022 - 2024 this asteroid was at visible apparent magnitudes brighter than 18.5, allowing for a detailed characterization using ground-based telescope observations. We determined its rotation period $P~=~5.021516\pm0.000106$ h and its absolute magnitude H = 18.78 $\pm$ 0.14 and. The large number of lightcurves allows to estimate its axes ratio, its convex shape and its pole orientation $\lambda = 301^{\circ} \pm 35^{\circ}$, $\beta = {89^{+1}_{-6}}^{\circ}$ and $\epsilon = 5^{\circ} \pm 3^{\circ}$ which indicate a prograde rotation. We report the semi-axis of the equivalent ellipsoid, $a$ = 0.42$^{+0.08}_{0.06}$ km, $b$ = 0.16$^{+0.05}_{0.04}$ km, and $c$ = $0.17\pm0.03$ km. Consequently, the volume equivalent diameter is $D_{eq}$ = $0.44 \pm 0.06$ km . Using observations conducted simultaneously with four broadband filters, we determined $(g-r) = 0.663 \pm 0.022$ mag, $(r-i) = 0.177 \pm 0.012$ mag, and $(i-z_s) = -0.061 \pm 0.032$ mag. Additionally, we found that Torifune exhibits no detectable large-scale heterogeneity. We classified the object using a high signal-to-noise ratio spectrum (over the visible and near-infrared region) as Sq-type in the Bus-DeMeo taxonomy. We estimate a mineralogy similar to LL/L ordinary chondrites, with an ol/(ol+px) = 0.60, a Fa content of 28.5 mol$\%$, and a Fs content of 23.4 mol$\%$. The spectral data indicate a surface affected by moderate space weathering effects.

Extreme Ultraviolet (EUV) waves, frequently produced by eruptions, propagate through the non-uniform magnetic field of the solar corona and interact with distant prominences, inducing their global oscillations. However, the generation, propagation, and interaction of these waves with distant prominences remain poorly understood. We aim to study the influence of an eruptive flux rope (EFR) on a distant prominence by means of extreme-resolution numerical simulations. We cover a domain of horizontal extent of 1100 Mm while capturing details down to 130 km using automated grid refinement. We performed a 2.5D numerical experiment using the open-source MPI-AMRVAC 3.1 code, modeling an eruption as a 2.5D catastrophe scenario augmented with a distant dipole magnetic field to form the flux rope prominence. Our findings reveal that the EFR becomes unstable and generates a quasicircular front. The primary front produces a slow secondary front when crossing equipartition lines where the Alfvén speed is close to the sound speed. The resulting fast and slow EUV waves show different behaviors, with the fast EUV wave slightly decelerating as it propagates through the corona, while the slow EUV wave forms a stationary front. The fast EUV wave interacts with the remote prominence, driving both transverse and longitudinal oscillations. Additionally, magnetic reconnection at a null point below the prominence flux rope is triggered by the fast EUV wave, affecting the flux rope magnetic field and the prominence oscillations. Our study unifies important results of the dynamics of eruptive events and their interactions with distant prominences, including details of (oscillatory) reconnection and chaotic plasmoid dynamics. We demonstrate for the first time the full consequences of remote eruptions on prominence dynamics and clarify the damping mechanisms of prominence oscillations.

Astrophysical transients can be powered by a broad range of energy sources including shock-heating (internal and external shocks), decay of radioactive isotopes, and long-lived central engines (magnetar and fallback). The dominant energy source for astrophysical transients depends on the nature of the explosive engine and its progenitor. To model all transients, light-curve codes must include all of these energy sources. Here we present a supernova light-curve code implementing analytic source models to compare the role of different energy sources in these transients. To demonstrate the utility of this code, we conduct an extensive study of type Ic broad-line supernovae. A diverse set of energy sources have been linked to Ic broad-line supernovae making them an excellent candidate for this light-curve code. In this paper, we explore which features of the explosion (mass, velocity, etc.) affect the type Ic supernovae light-curves, focusing on shock-interaction and radioactive decay energy sources. Although the explosion properties under both energy sources can be tuned to match the peak emission, matching the light-curve evolution in many Ic broad-line supernovae requires fine-tuned conditions. We find that shock interactions in the stellar wind are likely to be the dominant energy source at peak for these supernovae.

Recent measurements of the Hubble constant using type Ia supernovae explicitly correct for their estimated peculiar velocities using the 2M++ reconstruction of the local density field. The amount of uncertainty from this reconstruction procedure has thus far been unquantified. To rectify this, we use mock universe realisations of the 2M++ catalogue and generate predicted peculiar velocities using the same method as the predictions that are used to correct for the Pantheon+ catalogue. We find that the method yields uncertainties of 0.3 km/s/Mpc and hence subdominant to the total uncertainty in H0.

We present a composite machine learning framework to estimate posterior probability distributions of bulge-to-total light ratio, half-light radius, and flux for Active Galactic Nucleus (AGN) host galaxies within $z<1.4$ and $m<23$ in the Hyper Supreme-Cam Wide survey. We divide the data into five redshift bins: low ($0<z<0.25$), mid ($0.25<z<0.5$), high ($0.5<z<0.9$), extra ($0.9<z<1.1$) and extreme ($1.1<z<1.4$), and train our models independently in each bin. We use PSFGAN to decompose the AGN point source light from its host galaxy, and invoke the Galaxy Morphology Posterior Estimation Network (GaMPEN) to estimate morphological parameters of the recovered host galaxy. We first trained our models on simulated data, and then fine-tuned our algorithm via transfer learning using labeled real data. To create training labels for transfer learning, we used GALFIT to fit $\sim 20,000$ real HSC galaxies in each redshift bin. We comprehensively examined that the predicted values from our final models agree well with the GALFIT values for the vast majority of cases. Our PSFGAN + GaMPEN framework runs at least three orders of magnitude faster than traditional light-profile fitting methods, and can be easily retrained for other morphological parameters or on other datasets with diverse ranges of resolutions, seeing conditions, and signal-to-noise ratios, making it an ideal tool for analyzing AGN host galaxies from large surveys coming soon from the Rubin-LSST, Euclid, and Roman telescopes.

The size of a protoplanetary disk is a fundamental property, yet most remain unresolved, even in nearby star-forming regions (d $\sim$ 140-200 pc). We present the complete continuum size distribution for the $105$ brightest protoplanetary disks (M$_{\text{dust}}$ $\gtrsim$ 2 M$_{\oplus}$) in the Ophiuchus cloud, obtained from ALMA Band 8 (410 GHz) observations at 0.05$^{\prime\prime}$ (7 au) to 0.15$^{\prime\prime}$ (21 au) resolution. This sample includes 54 Class II and 51 Class I and Flat Spectrum sources, providing a comprehensive distribution across evolutionary stages. We measure the Half Width at Half Maximum (HWHM) and the radius encircling $68\%$ of the flux ($R_{68\%}$) for most non-binary disks, yielding the largest flux-limited sample of resolved disks in any star-forming region. The distribution is log-normal with a median value of $\sim$14 au and a logarithmic standard deviation $\sigma_{\log} = 0.46$ (factor of 2.9 in linear scale). Disks in close binary systems ($<$ 200 au separation) have smaller radii, with median value of $\sim$5 au, indicating efficient radial drift as predicted by dust evolution models. The size distribution for young embedded objects (SED Class I and Flat Spectrum, age $\lesssim$ 1 Myr) is similar to that of Class II objects (age $\sim$ a few Myr), implying that pressure bumps must be common at early disk stages to prevent mm-sized particle migration at au scales.

We present a detailed analysis of acetylene (C$_2$H$_2$) and its isotopologues in the Orion IRc2 region, focusing on the determination of $^{12}$C/$^{13}$C isotopic ratios using high-resolution infrared spectra from SOFIA. By employing a robust $\chi^2$ fitting method, we simultaneously determined temperature and column density, achieving a $^{12}$C/$^{13}$C ratio of $18.72^{+1.54}_{-1.46}$ for the blue clump and $15.07^{+1.61}_{-1.60}$ for the red clump. These results revealed significant discrepancies with the traditional rotational diagram method, which overestimated the ratios by 12.1% and 23.9%, respectively. Our $\chi^2$ approach also reduced uncertainties by up to 80% , providing more precise and reliable isotopic ratios. Additionally, we extended the analysis to isotopologues not covered in HITRAN, calculating vibrational and rotational constants through quantum chemical calculations. This allowed us to model subtle isotopic shifts induced by $^{13}$C and deuterium substitution, enabling accurate isotopologue detection in astrophysical environments. The Python package developed in this study facilitates efficient $\chi^2$ fitting and isotopic ratio analysis, making it a valuable tool for future high-resolution observations. This work highlights the critical role of advanced spectral models and fitting techniques in understanding isotopic fractionation and the chemical evolution of interstellar matter.

Accretion and migration usually proceeds concurrently for giant planet formation in the natal protoplanetary disks. Recent works indicate that the concurrent accretion onto a giant planet imposes significant impact on the planetary migration dynamics in the isothermal regime. In this work, we carry out a series of 2D global hydrodynamical simulations with Athena++ to explore the effect of thermodynamics on the concurrent accretion and migration process of the planets in a self-consistent manner. The thermodynamics effect is modeled with a thermal relaxation timescale using a $\beta$-cooling prescription. Our results indicate that radiative cooling has a substantial effect on the accretion and migration processes of the planet. As cooling timescales increase, we observe a slight decrease in the planetary accretion rate, and a transition from the outward migrating into inward migration. This transition occurs approximately when the cooling timescale is comparable to the local dynamical timescale ($\beta\sim1$), which is closely linked to the asymmetric structures from the circumplanetary disk (CPD) region. The asymmetric structures in the CPD region which appear with an efficient cooling provide a strong positive torque driving the planet migrate outward. However, such a positive torque is strongly suppressed, when the CPD structures tend to disappear with a relatively long cooling timescale ($\beta\gtrsim10$). Our findings may also be relevant to the dynamical evolution of accreting stellar-mass objects embedded in disks around active galactic nuclei.

PSR J2021+4026 is a gamma-ray pulsar having variations in its spin-down rate and gamma-ray flux. Its variations in timing and emission are correlated, e.g., a larger spin-down rate for a low gamma-ray flux. We show that the mode change in PSR J2021+4026 can be understood in the precession scenario. In the precession model, the inclination angle is modulated due to precession. At the same time, the wobble angle may decay with time. This results in damping of the precession. Combined with magnetospheric torque model and the outer gap model, the damped precession can explain: (1) when the inclination angle is larger, the spin-down rate will be larger, accompanied by a lower gamma-ray flux. (2) The variation amplitude of the gamma-ray flux and spin-down rate is smaller than previous results due to the damping of the precession. The modulation period is becoming shorter due to a smaller wobble angle. In the end, we propose that there are two kinds of modulations in pulsars. Long-term modulations in pulsars may be due to precession. Short-term modulations may be of magnetospheric origin.

Precise astrometric measurement with Gaia satellite resulted in the discovery of tens of wide binary systems consisting of a Sun-like star and an invisible component. The latter can be a white dwarf, a neutron star, or a black hole. In this paper, we model magneto-rotational evolution of neutron stars in wide low-mass binaries accounting for the orbital eccentricity. We aim to calculate when neutron stars in such systems can start to accrete matter from the stellar wind of the companion. We show that the transition from the ejector to the propeller stage occurs earlier in more eccentric systems, thus increasing the time that neutron stars can spend accreting matter. Our calculations show that in the case of efficient spin-down at the propeller stage, a neutron star in an eccentric orbit with $e\gtrsim0.6$ and a standard magnetic field $B=10^{12}$ G can start accreting within a few Gyr. For neutron stars with $B=10^{13}$ G the onset of accretion occurs earlier regardless of the orbital eccentricity. Otherwise, with a lower spin-down rate, such a neutron star will remain at the propeller stage for most of its life.

Space weather predictions are necessary to avoid damage caused by intense geomagnetic storms. Such strong storms are usually caused by a co-rotating interaction region (CIR) passing at Earth or by the arrival of strong coronal mass ejections (CMEs). To mitigate the damage, the effect of propagating CMEs in the solar wind must be estimated accurately at Earth and other locations. Modelling solar wind accurately is crucial for space weather predictions, as it is the medium for CME propagation. The Icarus heliospheric modelling tool is upgraded to handle dynamic inner heliospheric driving instead of using steady boundary conditions. The ideal magnetohydrodynamic (MHD) solver and the automated grid-adaptivity are adjusted to the latest MPI-AMRVAC version. The inner boundary conditions, prescribed at 0.1 AU for the heliospheric model, are updated time-dependently throughout the simulation. The coronal model is computed repeatedly for selected magnetograms. The particle sampling within MPI-AMRVAC is extended to handle stretched spherical grid information. The solar wind obtained in the simulation is dynamic and shows significant variations throughout the evolution. When comparing the results with the observations, the dynamic solar wind results are more accurate than previous results obtained with purely steady boundary driving. The CMEs propagated through the dynamic solar wind background produce more similar signatures in the time-series data than in the steady solar wind. Dynamic boundary driving in Icarus results in a more self-consistent solar wind evolution in the inner heliosphere. The obtained space weather modelling tool for dynamic solar wind and CME simulations is better suited for space weather forecasting than a steady solar wind model.

Protoplanetary disks have been found around free-floating objects with masses comparable to those of giant planets. The frequency and properties of these disks around planetary-mass objects are still debated. Here we present ultradeep mid-infrared images for the young cluster IC348, obtained through stacking of time series images from Spitzer. We measure fluxes at 3.6 and 4.5 microns for known free-floating planetary-mass objects (FFPMOs, spectral type M9 or later) in this cluster. By comparing the observed infrared spectral energy distributions with photospheric templates, we identify six planetary-mass objects with disks, plus three which may or may not have a disk. This corresponds to a disk fraction of 46% (34-59%). The disk fraction among planetary-mass objects is comparable to more massive brown dwarfs. We show the disk fraction among free-floating planetary-mass objects as a function of age, demonstrating that these objects retain disks for several million years, similar to low-mass stars and brown dwarfs.

We report an extensive rotational spectroscopic analysis of singly deuterated methyl mercaptan (CH$_{2}$DSH) using both millimeter and far-infrared synchrotron spectra to achieve a global torsional analysis of the three lowest torsional substates (e$_{0}$, e$_{1}$, and o$_{1}$) of this non-rigid species. A fit including 3419 millimeter wave transitions along with 43 infrared torsional subband centers was performed with root mean square deviations of 0.233 MHz and 0.270 cm$^{-1}$, respectively, resulting in 68 fit parameters. A spectroscopic catalogue built from this analysis for a temperature of 125 K has led to the first interstellar detection of CH$_{2}$DSH towards the Solar-like protostar IRAS 16293-2422 B. We report the identification of 46 transitions, including eight relatively unblended lines, resulting in a derived column density of (3.0$\pm$0.3)$\times$10$^{14}$ cm$^{-2}$. The column density ratio for HDCS/CH$_{2}$DSH compared to HDCO/CH$_{2}$DOH suggests a difference in the interstellar chemistry between the sulphur and oxygen complex organics, in particular a different link between H$_{2}$CO and CH$_{3}$OH and between H$_{2}$CS and CH$_{3}$SH. This is the first interstellar detection of a deuterated sulphur-bearing COM and therefore an important step into understanding the chemical origin of sulphur-based prebiotics.

During its most recent return, comet 12P/Pons-Brooks experienced 14 well-documented outbursts, observed between June 13, 2023, and April 2024, at heliocentric distances ranging from $4.26\,$au to $0.85\,$au. After perihelion, the comet experienced two additional outbursts between June and August 31, 2024, at heliocentric distances of $1.20\,$au and $2.26\,$au, respectively. Using observational data, we developed a numerical model to estimate the mass ejected during these events, focusing on the sublimation of ice through the porous cometary nucleus. The key factors affecting ejected mass estimates are the outburst amplitude and the active surface area during both quiet sublimation and the outburst phases. Pogson's law was used to express outburst magnitude, incorporating scattering cross-sections of cometary agglomerates. The model iteratively determined the mass ejected in observed outbursts, considering various ice types (H$_{2}$O and CO$_{2}$) controlling sublimation activity. Our results indicate that the mass ejected during these outbursts ranges from 10$^{10}$ to 10$^{13}$ kg. Our findings highlight the significant role of surface morphology and thermodynamic conditions in cometary outbursts, providing insights into the mechanisms driving these phenomena and their implications for cometary evolution and dust trail formation. Based on the analysis of observational data, we propose a six-level classification system for cometary outbursts.

PS16dtm is one of the earliest reported candidate tidal disruption events (TDEs) in active galactic nuclei (AGNs) and displays a remarkably bright and long-lived infrared (IR) echo revealed by multi-epoch photometry from the Wide-field Infrared Survey Explorer (WISE). After a rapid rise in the first year, the echo remains persistently at a high state from July 2017 to July 2024, the latest epoch, and keeps an almost constant color. We have fitted the extraordinary IR emission with a refined dust echo model by taking into account the dust sublimation process. The fitting suggests that an extremely giant dust structure with a new inner radius of $\sim1.6$ pc and an ultra-high peak bolometric luminosity, i.e., $\sim6\times10^{46} \rm erg~s^{-1}$ for typical 0.1$\mu$m-sized silicate grain, is required to account for the IR echo. This work highlights the distinctive value of IR echoes in measuring the accurate intrinsic bolometric luminosity, and thus the total radiated energy of TDEs, which could be severely underestimated by traditional methods, i.e. probably by more than an order of magnitude in PS16dtm. Such large energetic output compared to normal TDEs could be boosted by the pre-existing accretion disk and gas clouds around the black hole. Our model can be validated in the near future by IR time-domain surveys such as Near-Earth Object (NEO) Surveyor, given the recent retirement of WISE. In addition, the potential for spatially resolving a receding dusty torus after a TDE could also be an exciting subject in the era of advanced IR interferometry.

Some white dwarfs undergo significant changes in atmospheric composition owing to the diffusion and mixing of residual hydrogen in a helium-rich envelope. Of particular interest are a few objects exhibiting hydrogen and helium line variations modulated by rotation, revealing surface composition inhomogeneities. Recently, the hot ultramassive white dwarf ZTF J203349.80+322901.1 emerged as the most extreme such specimen, with hydrogen and helium lines successively appearing and vanishing in anti-phase, suggesting a peculiar double-faced configuration. However, standard atmosphere models fail to reproduce the observed spectrum at all rotation phases, hampering further interpretation. Here, we perform a new analysis of ZTF J203349.80+322901.1 using stratified atmosphere models, where hydrogen floats above helium, and obtain excellent fits to the phase-resolved spectra. Our results imply that an extremely thin hydrogen layer covers the entire surface but varies from optically thick to optically thin across the surface, thus producing the observed spectral variations. We present new envelope models indicating that the hydrogen layer arises from a delicate interplay between diffusion and convection. We discuss possible explanations for the surface layer asymmetry, including an asymmetric magnetic field and a non-uniform internal hydrogen distribution. Finally, we highlight implications for expanding and understanding the emerging class of inhomogeneous white dwarfs.

The cicumgalactic medium (CGM) is a reservoir of metals and star-forming fuel. Most baryons in the universe are in the circumgalactic medium (CGM) or intergalactic medium (IGM). The baryon cycle -- how mass and metals reach the CGM from the inner regions of the galaxy and how gas from the CGM replenishes star-forming activity in the inner regions -- is an essential question in galaxy evolution. In this paper, we study the flow of mass and metals in a stacked sample of 2770 isolated halos from the IllustrisTNG cosmological hydrodynamic simulation. The mean gas flow as a function of radius and angle is similar across a large galactic mass range when accounting for different feedback modes. Although both star formation and black holes cause powerful outflows, the flows from star formation are more angularly restricted. Black hole feedback dominates massflow throughout the halo, while star-formation feedback mainly affects the inner region. When scaling by virial radius ($R_v$), large dynamical changes occur at $0.2R_v$ for most halos, suggesting a characteristic size for the inner galaxy. Despite radio mode feedback from black holes being the primary quenching mechanism in IllustrisTNG, a small population of high mass radio mode disks are able to form stars.

Stratospheric balloons offer cost-effective access to space and grant the opportunity for fast scientific innovation cycles and higher-risk explorations. In addition to science pathfinders, they serve as platforms for technology advancement, and offer a unique opportunity to train future instrument scientists and PIs. However, the increase in complexity of some projects (sub-arcsecond pointing, numerous degrees of freedom, advanced cooling systems, real-time communication and data transfer, low readiness level technologies, etc.) elevates the scale of challenges. This paper discusses the challenges brought by the increase of instrument complexity in the constrained area of ballooning projects. We use the example of the multi-object slit UV spectrograph FIREBall-2 but also discuss other ambitious payloads to expose some lessons learned. We will then propose strategies for future balloon projects and potential interesting adaptations for funding agencies to accommodate these emerging complexities. The objective being to foster the added value of more ambitious balloon-borne projects, and to initiate discussion about how to generate more robust strategies for handling high-complexity undertakings. Although this paper is enriched by answers from a survey sent to balloon-borne payload PIs, it remains influenced by my personal experience on FIREBall-2. As such, it does not necessarily represent the perspectives of other members of the project, let alone the broader balloon community.

Gamma-ray bursts (GRBs) are among the most luminous electromagnetic transients in the universe, providing unique insights into extreme astrophysical processes and serving as promising probes for cosmology. Unlike Type Ia supernovae, which have a unified explosion mechanism, GRBs cannot directly act as standard candles for tracing cosmic evolution at high redshifts due to significant uncertainties in their underlying physical origins. Empirical correlations derived from statistical analyses involving various GRB parameters provide valuable information regarding their intrinsic properties. In this brief review, we describe various correlations among GRB parameters involving the prompt and afterglow phases, discussing possible theoretical interpretations behind them. The scarcity of low-redshift GRBs poses a major obstacle to the application of GRB empirical correlations in cosmology, referred to as the circularity problem. We present various efforts aiming at calibrating GRBs to address this challenge and leveraging established empirical correlations to constrain cosmological parameters. The pivotal role of GRB sample quality in advancing cosmological research is underscored. Some correlations that could potentially be utilized as redshift indicators are also introduced.

The polarisation of light induced by aligned interstellar dust serves as a significant tool in investigating cosmic magnetic fields, dust properties, and poses a challenge in characterising the polarisation of the cosmic microwave background and other sources. To establish dust polarisation as a reliable tool, the physics of the grain alignment process needs to be studied thoroughly. The Magnetically enhanced Radiative Torque (MRAT) alignment is the only mechanism that can induce highly efficient alignment of grains with magnetic fields required by polarisation observations of the diffuse interstellar medium. Here, we aim to test the MRAT mechanism in starless cores using the multiwavelength polarisation from optical/NIR to far-IR/submm. Our numerical modelling of dust polarisation using the MRAT theory demonstrated that the alignment efficiency of starlight polarisation ($p_{\rm ext}/A_{\rm V}$) and the degree of thermal dust polarisation ($p_{\rm em}$) first decrease slowly with increasing visual extinction ($A_{\rm V}$) and then falls steeply as $\propto A^{-1}_{\rm V}$ at large $A_{\rm V}$ due to the loss of grain alignment, which explains the phenomenon known as polarisation holes. Visual extinction at the transition from shallow to steep slope ($A^{\rm loss}_{\rm V}$) increases with the maximum grain size. By applying physical profiles suitable for a starless core Pipe-109, our model successfully reproduces the existing observations of starlight polarisation at R-band ($0.65\,\mu$m) and H-band ($1.65\,\mu$m), as well as emission polarisation at submillimeter ($850\,\mu$m). Successful modelling of observational data requires perfect alignment of large grains as evidence of the MRAT mechanism, and larger maximum size with higher elongation at higher $A_{\rm V}$. The latter reveals the first evidence for the new model of anisotropic grain growth induced by magnetic grain alignment.

This study investigates the thermal stability of thin accretion disks in high-energy astrophysical systems, incorporating the effects of magnetic fields, winds, and coronae. We analyze how these factors influence disk stability, focusing on conditions under which magnetic fields enhance stability and on scenarios where winds and coronae can either stabilize or destabilize the disk. Our results reveal that increasing corona parameters raises disk thickness and reduces temperature, thereby affecting gas, radiation, and magnetic pressures. These interactions underscore the complex dependencies that shape accretion disk dynamics, offering insights into their structural and thermal behavior under varying physical conditions. The findings contribute to advancing theoretical models and numerical simulations of accretion processes in environments such as active galactic nuclei (AGN) and X-ray binaries, where disk stability plays a critical role in observed emissions and variability patterns.

Context. Star clusters are often invoked as contributors to the flux of Galactic cosmic rays and as sources potentially able to accelerate particles to $\sim$PeV energies. The gamma radiation with $E\gtrsim$ TeV recently observed from selected star clusters has profound implications for the origin of Galactic cosmic rays. Aims. We show that if the gamma rays observed from the Cygnus cocoon and Westerlund 1 are of hadronic origin, then the cosmic rays escaping the cluster at energies $\gtrsim$ 10 TeV must cross a grammage inside the cluster that exceeds the Galactic grammage. At lower energies, depending on the model adopted to describe the production of gamma rays, such grammage may exceed or be comparable with the grammage inferred from propagation on Galactic scales. Methods. The flux of gamma rays is analytically computed for a few models of injection of cosmic rays in star clusters, and compared with the flux measured from selected clusters. Results. In all models considered here, comparing the inferred and observed gamma ray fluxes at $E\gtrsim$ TeV, we conclude that CRs must traverse a large grammage inside or around the cluster before escaping. Clearly these implications would not apply to a scenario in which gamma rays are produced due to radiative losses of leptons in the cluster. Leptonic models typically require weaker magnetic fields, which in turn result in maximum energies of accelerated particles much below $\sim$ PeV. Conclusions. We conclude that if gamma ray emission in SCs is a generic phenomenon and if this radiation is due to hadronic interactions, either star clusters cannot contribute but a small fraction of the total cosmic ray flux at the Earth, or their contribution to the grammage cannot be neglected and the paradigm of Galactic transport should be profoundly revisited.

Several observational lines of evidence imply that a fraction of the dark matter in the Galaxy may be comprised of small cold clouds of molecular hydrogen. Such objects are difficult to detect because of their small size and low temperature, but they can reveal themselves with gamma radiation arising in interactions between such clouds and cosmic rays or as dark shadows cast on the optical, UV and X-ray sky background. In our work we use the data of Fermi LAT 4FGL-DR4 catalogue of gamma-ray sources together with the data of GALEX UV All-Sky Survey to search for small dark clouds of molecular hydrogen in the Solar neighbourhood. This approach allows us to put an upper limit on the local concentration of such objects: $n < 2.2 \times 10^{-11} {\mathrm{AU}^{-3}}$. Constraints (upper limits) on the total amount of matter in this form bound to the Sun strongly depend on the radial profile of the clouds' distribution and reside in $0.05-30~M_{\odot}$ mass range.

The cosmic microwave background (CMB) experiments have reached an era of unprecedented precision and complexity. Aiming to detect the primordial B-mode polarization signal, these experiments will soon be equipped with $10^{4}$ to $10^{5}$ detectors. Consequently, future CMB missions will face the substantial challenge of efficiently processing vast amounts of raw data to produce the initial scientific outputs - the sky maps - within a reasonable time frame and with available computational resources. To address this, we introduce \texttt{BrahMap}, a new map-making framework that will be scalable across both CPU and GPU platforms. Implemented in C++ with a user-friendly Python interface for handling sparse linear systems, \texttt{BrahMap} employs advanced numerical analysis and high-performance computing techniques to maximize the use of super-computing infrastructure. This work features an overview of the \texttt{BrahMap}'s capabilities and preliminary performance scaling results, with application to a generic CMB polarization experiment.

This paper presents the first effort to Extend the Investigation of Stellar Populations In RElics (E-INSPIRE). We present a catalogue of 430 spectroscopically-confirmed ultra-compact massive galaxies (UCMGs) from the Sloan Digital Sky Survey at redshifts $0.01<z<0.3$. This increases the original INSPIRE sample eightfold, bridging the gap with the local Universe. For each object, we compute integrated stellar velocity dispersion, age, metallicity, and [Mg/Fe] through spectroscopic stellar population analysis. We infer star formation histories (SFHs), metallicity evolution histories (MEHs) and compute the Degree of Relicness (DoR) of each object. The UCMGs, covering a wide range of DoR from 0.05 to 0.88, can be divided into three groups, according to how extreme their SFH was. The first group consists of 81 extreme relics ($\text{DoR}\gtrsim0.6$) that have formed the totality of their stellar mass by $z\sim2$ and have super-solar metallicities at all cosmic epochs. The second group ($0.3\lesssim\text{DoR}\lesssim0.6$) contains 293 objects also characterised by peaked SFHs but with a small percentage of later-formed stars and with a variety of MEHs. The third group ($\text{DoR}\lesssim0.3$), has 56 objects that cannot be considered relics since they have extended SFHs and formed a non-negligible fraction ($>25\%$) of their stellar mass at $z<2$. We confirm that an efficient method of finding relics is to select UCMGs with large velocity dispersion values but we believe that the most efficient way is to select high velocity dispersion objects that also have super-solar metallicities and high [Mg/Fe].

Nanodust particles produced near the Sun by collisional breakup of larger grains are accelerated in the magnetised solar wind and reach high speeds outwards of 1 AU. Vaporisation and ionisation of fast dust grains impacting a spacecraft produce voltage pulses on wave instruments that enable them to act as dust detectors. Wave instruments on STEREO and on Cassini during its cruise phase detected a highly variable flux of fast nanodust. Both detections took place when the orientation of the solar magnetic dipole produced an interplanetary electric field that focused nanoparticles towards the heliospheric current sheet (HCS) - a geometry that is recurring because of the periodicity of solar activity.

The stellar-to-halo mass relation (SHMR) embodies the joint evolution of galaxies and their host dark matter halos. However, the relation is poorly constrained at sub-galactic masses, because the stellar emission from such objects is so faint. However, it is possible to directly detect the mass of halos along the line of sight to a strong gravitational lens, when they perturb one of its multiple images. Space telescopes including Euclid, CSST, and Roman will soon discover millions of galaxy-galaxy strong lensing systems. We simulate Euclid-like imaging of a typical lens galaxy, and find that a lensing reconstruction is sensitive to $3\times10^{10}$ subhalos with various positions and concentrations, at statistical signficance $>$$3.6\sigma$. The subhalo mass can be measured without bias, provided the model simultaneously fits light from both the main lens and the subhalo. A future sample of 48 subhalos with $\geqslant$$5\sigma$ detection significance would constrain the SHMR at this mass range with $1\sigma$ uncertainty of 0.045 dex: distinguishing between different theoretical predictions at the sub-galactic scale. Follow-up spectroscopy is needed to measure lens and source redshifts; follow-up imaging at greater spatial resolution and depth would substantially improve the measurement, and eliminate false-positives at even lower halo masses.

The spectral energy distribution (SED) of observed stars in wide-field images is crucial for chromatic point spread function (PSF) modelling methods, which use unresolved stars as integrated spectral samples of the PSF across the field of view. This is particularly important for weak gravitational lensing studies, where precise PSF modelling is essential to get accurate shear measurements. Previous research has demonstrated that the SED of stars can be inferred from low-resolution observations using machine-learning classification algorithms. However, a degeneracy exists between the PSF size, which can vary significantly across the field of view, and the spectral type of stars, leading to strong limitations of such methods. We propose a new SED classification method that incorporates stellar spectral information by using a preliminary PSF model, thereby breaking this degeneracy and enhancing the classification accuracy. Our method involves calculating a set of similarity features between an observed star and a preliminary PSF model at different wavelengths and applying a support vector machine to these similarity features to classify the observed star into a specific stellar class. The proposed approach achieves a 91\% top-two accuracy, surpassing machine-learning methods that do not consider the spectral variation of the PSF. Additionally, we examined the impact of PSF modelling errors on the spectral classification accuracy.

We study the deflections of ultra-high-energy cosmic rays in several widely used models of the coherent Galactic magnetic field (GMF), including PT11 (Pshirkov et al. [1]), JF12 (Jansson and Farrar [2]), UF23 (Unger and Farrar [3]) and KST24 (Korochkin, Semikoz, and Tinyakov [4]). We propagate particles with rigidities of 5, 10, and 20 EV and analyze the differences in deflection predictions across these GMF models. We identify the GMF components responsible for deflections in various regions of the sky and discuss the uncertainties in modeling these components, as well as potential future improvements.

Although dust constitutes only about 1% of the mass of a protoplanetary disk, recent studies demonstrate that it can exert a significant torque on low- and intermediate-mass planetary cores. We compute and quantify for the first time the influence of the dust torque on the evolution of growing planetary embryos as they move in a protoplanetary disk while growing via gas and pebble accretion. Our global model evolves the gaseous disk via viscous accretion and X-ray photoevaporation, while accounting for dust growth and evolution including coagulation, drift, and fragmentation. Our research indicates that dust torque significantly influences planetary migration, particularly driving substantial outward migration for planets forming within the water ice-line. This effect occurs due to an increased dust-to-gas mass ratio in the inner disk, resulting from inward pebble drift from outer regions. In contrast, for planets initially located beyond the water ice-line, the dust torque mitigates inward migration but does not significantly alter their paths, as the dust-to-gas ratio diminishes rapidly due to rapid pebble drift and the brief timescales of planet formation in these areas. These findings underscore the pivotal role of dust torque in shaping the migration patterns of low- and intermediate-mass planets, especially when enhanced dust concentrations in the inner disk amplify its effects

Context: This paper focuses on a class of galaxies characterised by an extremely low surface brightness: the ultra-diffuse galaxies (UDGs). We used new integral-field spectroscopic data from the ESO Large Programme Looking into the faintEst WIth MUSE (LEWIS) project. Aims: Our main goals are addressing the formation channels and investigating possible correlations of their observational properties. In particular, we derive their stellar kinematics and dynamical properties. Methods: We extract the 1D stacked spectrum inside the effective radius to obtain an unbiased measure of $\sigma_{\rm eff}$. To derive the spatially-resolved stellar kinematics, we first apply the Voronoi tessellation algorithm to bin the spaxels in the datacube and then follow the same prescription adopted for the 1D case. In addition, we extract the velocity profiles along the galaxy's major and minor axes. Results: We find that 7 out of 18 UDGs in LEWIS show a mild rotation, 5 do not have evidence of any rotation, and the remaining 6 UDGs are unconstrained cases. This is the first large census of velocity profiles for UDGs. On average, UDGs in LEWIS are characterised by low values of $\sigma_{\rm eff}$, comparable with available values from the literature. In the Faber-Jackson relation plane, we found a group of UDGs consistent with the relation within the errorbars, whereas outliers are objects with non-negligible rotation components. UDGs and LSBs in LEWIS have larger dark matter content than dwarf galaxies with similar total luminosity. We do not find clear correlations between the derived properties and the local environment. Conclusions: Based on the stellar kinematics, two classes of UDGs are found in the Hydra I cluster: the rotating and non-rotating systems. This result, combined with other structural properties, can help to discriminate between the several formation scenarios proposed for UDGs.

Crystallization of the lunar magma ocean yielded a chemically unique liquid residuum named KREEP. This component is expressed as a large patch on the near side of the Moon, and a possible smaller patch in the northwest portion of the Moon's South Pole-Aitken basin on the far side. Thermal models estimate that the crystallization of the lunar magma ocean (LMO) could have spanned from 10 and 200 Myr, while studies of radioactive decay systems have yielded inconsistent ages for the completion of LMO crystallization covering over 160 Myr. Here, we show that the Moon achieved over 99 percent crystallization at 4429+/-76 Myr, indicating a lunar formation age of 4450 Myr or possibly older. Using the 176Lu-176Hf decay system (t1/2=37 Gyr), we found that the initial 176Hf/177Hf ratios of lunar zircons with varied U-Pb ages are consistent with their crystallization from a KREEP-rich reservoir with a consistently low 176Lu/177Hf ratio of 0.0167 that emerged ~140 Myr after solar system formation. The previously proposed younger model age of 4.33 Ga for the source of mare basalts (240 Myr after solar system formation) might reflect the timing of a large impact. Our results demonstrate that lunar magma ocean crystallization took place while the Moon was still battered by planetary embryos and planetesimals leftover from the main stage of planetary accretion. Study of Lu-Hf model ages for samples brought back from the South Pole-Aitken basin will help to assess the lateral continuity of KREEP and further understand its significance in the early history of the Moon.

We investigate the numerical stability of thermonuclear detonations in 1D accelerated reactive shocks and 2D binary collisions of equal mass, magnetized and unmagnetized white dwarf stars. To achieve high resolution at initiation sites, we devised geometric gridding and mesh velocity strategies specially adapted to the unique requirements of head-on collisional geometries, scenarios in which one expects maximum production of iron-group products. We study effects of grid resolution and the limiting of temperature, energy generation, and reactants for different stellar masses, separations, magnetic fields, initial compositions, detonation mechanisms, and limiter parameters across a range of cell sizes from 1 to 100 km. Our results set bounds on the parameter space of limiter amplitudes for which both temperature and energy limiting procedures yield consistent and monotonically convergent solutions. Within these bounds we find grid resolutions of 5 km or better are necessary for uncertainties in total released energy and iron-group products to drop below 10%. Intermediate mass products (e.g., calcium) exhibit similar convergence trends but with somewhat greater uncertainty. These conclusions apply equally to pure C/O WDs, multi-species compositions (including helium shells), magnetized and unmagnetized cores, and either single or multiple detonation scenarios.

Inhomogeneities along the line of sight in strong gravitational lensing distort the images produced, in an effect called shear. If measurable, this shear may provide independent constraints on cosmological parameters, complementary to traditional cosmic shear. We model 50 strong gravitational lenses from the Sloan Lens ACS (SLACS) catalogue with the aim of measuring the line-of-sight (LOS) shear for the first time. We use the `minimal model' for the LOS shear, which has been shown to be theoretically safe from degeneracies with lens model parameters, a finding which has been confirmed using mock data. We use the dolphin automated modelling pipeline, which uses the lenstronomy software as a modelling engine, to model our selected lenses. We model the main deflector with an elliptical power law profile, the lens light with elliptical Sérsic profiles and the source with a basis set of shapelets and an elliptical Sérsic profile. We successfully obtain a line-of-sight shear measurement from 18 of the 50 lenses. We find that these LOS shear measurements are consistent with external shears measured in recent works using a simpler shear model, which are larger than those expected from weak lensing. Neglecting the post-Born correction to the potential of the main deflector due to foreground shear leads to a propagation of degeneracies to the LOS shear measurement, and the same effect is seen if a prior is used to connect the lens mass and light ellipticities. The inclusion of an octupole moment in the lens mass profile does not lead to shear measurements that are in better agreement with the expectations from weak lensing.

Increasingly realistic simulations of the corona are used to predict synthetic observables for instruments onboard both existing and upcoming heliophysics space missions. Synthetic observables play an important role in constraining coronal heating theories. Choosing the spatial resolution of numerical simulations involves a trade-off between accuracy and computational cost. Since the numerical resolution not only affects the scale of structures that can be resolved, but also thermodynamic quantities such as the average coronal density, it is important to quantify the effect on synthesized observables. Using 3D radiative MHD simulations of coronal loops at three different grid spacings, from 60 km down to 12 km, we find that changes in numerical resolution lead to differences in thermodynamic quantities and stratification as well as dynamic behaviour. Higher grid resolution results in a more complex and dynamic atmosphere. The resolution affects the emission intensity as well as the velocity distribution, thereby affecting synthetic spectra derived from the simulation. The distribution of synthetic coronal loop strand sizes changes as more fine-scale structure is resolved. A number of parameters, however, seem to start to saturate from our chosen medium grid resolution on. Our study shows that while choosing a sufficiently high resolution matters when comparing forward-modelled observables with data from current and future space missions, for most purposes not much is gained by further increasing the resolution beyond a grid spacing of 24 km, which seems to be adequate for reproducing bulk loop properties and forward-modelled emission, representing a good trade-off between accuracy and computational resource.

Different stellar populations may be identified through differences in chemical, kinematic, and chronological properties, suggesting the interplay of various physical mechanisms that led to their origin and subsequent evolution. As such, the identification of stellar populations is key for gaining insight into the evolutionary history of the Milky Way galaxy. This task is complicated by the fact that stellar populations share significant overlap in their chrono-chemo-kinematic properties, hindering efforts to identify and define stellar populations. Our goal is to offer a novel and effective methodology that can provide deeper insight into the nonlinear and nonparametric properties of the multidimensional physical parameters that define stellar populations. For this purpose we explore the ability of manifold learning to differentiate stellar populations with minimal assumptions about their number and nature. Manifold learning is an unsupervised machine learning technique that seeks to intelligently identify and disentangle manifolds hidden within the input data. To test this method, we make use of Gaia DR3-like synthetic stellar samples generated from the FIRE-2 cosmological simulations. These represent red-giant stars constrained by asteroseismic data from TESS. We reduce the 5-dimensional input chrono-chemo-kinematic parameter space into 2-dimensional latent space embeddings generated by manifold learning. We then study these embeddings to assess how accurately they represent the original data and whether they contain meaningful information that can be used to discern stellar populations. We conclude that manifold learning possesses promising abilities to differentiate stellar populations when considering realistic observational constraints.

The Time Domain Extragalactic Survey (TiDES) conducted on the 4-metre Multi-Object Spectroscopic Telescope (4MOST) will perform spectroscopic follow-up of extragalactic transients discovered in the era of the NSF-DOE Vera C. Rubin Observatory. TiDES will conduct a 5-year survey, covering ${>}14\,000\,\mathrm{square\, degrees}$, and use around 250 000 fibre hours to address three main science goals: (i) spectroscopic observations of ${>}30 000$ live transients, (ii) comprehensive follow-up of ${>}200 000$ host galaxies to obtain redshift measurements, and (iii) repeat spectroscopic observations of Active Galactic Nuclei to enable reverberation mapping studies. The live spectra from TiDES will be used to reveal the diversity and astrophysics of both normal and exotic supernovae across the luminosity-timescale plane. The extensive host-galaxy redshift campaign will allow exploitation of the larger sample of supernovae and improve photometric classification, providing the largest-ever sample of spec-confirmed type Ia supernovae, capable of a sub-2 per cent measurement of the equation-of-state of dark energy. Finally, the TiDES reverberation mapping experiment of 700--1\,000 AGN will complement the SN Ia sample and extend the Hubble diagram to $z\sim2.5$

We present a 3D shape analysis of both dark matter (DM) and stellar matter (SM) in simulated dwarf galaxies to determine whether stellar shape traces DM shape. Using 80 central and satellite galaxies from three simulation suites (Marvelous Massive Dwarfs, Marvelous Dwarfs, and DC Justice League) spanning stellar masses of $10^6$--$10^{10}$ M$_\odot$, we measure 3D shapes through the moment of inertia tensor at two times the effective radius to derive axis ratios ($C/A$, $B/A$) and triaxiality. We find that stellar shape does indeed follow DM halo shape for our dwarf galaxies. However, the presence of a stellar disk in more massive dwarfs (M$_* \gtrsim 10^{7.5}$ M$_\odot$) pulls the distribution of stellar $C/A$ ratios to lower values, while in lower mass galaxies the gravitational potential remains predominantly shaped by DM. Similarly, stellar triaxiality generally tracks dark matter halo triaxiality, with this relationship being particularly strong for non-disky galaxies though weaker in disky systems. This correlation is reinforced by strong alignment between SM and DM axes, particularly in disk galaxies. Further, we find no detectable difference in either SM or DM shape comparing two different SNe feedback implementations, demonstrating that shape measurements may be robust to different implementations of baryonic feedback in dwarf galaxies. We also observe that a dwarf galaxy's shape is largely unperturbed by recent mergers (with merger ratios $>4$). This comprehensive study demonstrates that stellar shape measurements can serve as a reliable tool for inferring DM shapes in dwarf galaxies.

The neutron star X-ray binary, EXO 0748--676, was observed regularly by the Rossi X-ray Timing Explorer (RXTE) and XMM-Newton during its first detected outburst (1985 - 2008). These observations captured hundreds of asymmetric, energy-dependent X-ray eclipses, influenced by the ongoing ablation of the companion star and numerous Type I thermonuclear X-ray bursts. Here, we present the light curves of 22 Type I X-ray bursts observed by RXTE that coincide, fully or partially, with an X-ray eclipse. We identify nine instances where the burst occurs entirely within totality, seven bursts split across an egress, and six cases interrupted by an ingress. All in-eclipse bursts and split bursts occurred while the source was in the hard spectral state. We establish that we are not observing direct burst emission during eclipses since the companion star and the ablated outflow entirely obscure our view of the X-ray emitting region. We determine that the reflected flux from the outer accretion disc, even if maximally flared, is insufficient to explain all observations of in-eclipse X-ray bursts and instead explore scenarios whereby the X-ray bursts are scattered, either by a burst-induced rise in $N_{\rm{H}}$ that provides extra material, an accretion disc wind or the ablated outflow into our line of sight. However, the rarity of a burst and eclipse overlap makes it challenging to determine their origin.

Weak-emission-line QSOs (WLQs) are an enigmatic subclass of the QSO population, as their optical/UV spectra are marked by abnormally weak (or absent) emission lines. To obtain much-needed additional clues to the origin of this and other known peculiarities of WLQs, we have determined the 'ensemble' optical variability characteristics for a large, well-defined sample of 76 radio-quiet WLQs and also for a matched control sample comprising 603 normal radio-quiet QSOs. This analysis was done using their light-curves recorded in the $g$ and $r$ bands, under the ZTF survey during 2018-2024, with a typical cadence of 3 days. We find that, compared to normal QSOs, WLQs exhibit systematically milder optical variability on month/year-like time scales (by a factor of $\sim$ 1.76$\pm$0.05 in amplitude). We have independently verified this by carrying out an equivalent analysis of the V-band light-curves acquired under the CRTS during 2007- 2014, with a typical cadence of 10 days. This new observational differentiator between WLQs and normal QSOs may provide clues to understanding the intriguing nature of WLQs. It is proposed that the clumpiness of the torus material flowing into the central engine may play a key role in explaining the observed differences between the WLQs and normal QSOs.

We investigate the gravitational impulse by using the generalized formulation of the quantum Boltzmann equation (QBE), wherein the initial states are taken as wave packets rather than plane waves. The QBE equation operates within an open quantum system framework. Using this approach, we can analyze the two-body gravitational scattering by considering one body as part of the environment and the other as the system. Similarly, we apply this approach to the gravitational impulse of a photon due to a massive object, during which we consider the photon as the system and the massive object as the environment. In this methodology, we use the forward scattering term of the extended QBE to compute the time evolution of the momentum operator associated with the system. While this term vanishes when considering point particles, we demonstrate its persistence when using wave packets to describe both massive particles and the photon field. Through this procedure, we obtain the leading order of the system's impulse. The results reported here are entirely consistent with the previous approaches.

Unimodular theory incorporating the Kaluza-Klein construction in five dimensions leads, after reduction to four dimensions, to a new class of scalar-tensor theory. The vacuum cosmological solutions display a bouncing, non singular behavior. From the four dimensional point of view, the solutions are completely regular. However, the propagation of gravitational waves in this geometry displays the presence of instabilities which reflect some features of the original five dimensional structure. Comparison with a four dimensional quantum model with cosmological constant, which has a similar background behavior, is discussed.

This paper discusses an approach to inertial parameter estimation for the case of cargo carrying spacecraft that is based on causal learning, i.e. learning from the responses of the spacecraft, under actuation. Different spacecraft configurations (inertial parameter sets) are simulated under different actuation profiles, in order to produce an optimised time-series clustering classifier that can be used to distinguish between them. The actuation is comprised of finite sequences of constant inputs that are applied in order, based on typical actuators available. By learning from the system's responses across multiple input sequences, and then applying measures of time-series similarity and F1-score, an optimal actuation sequence can be chosen either for one specific system configuration or for the overall set of possible configurations. This allows for both estimation of the inertial parameter set without any prior knowledge of state, as well as validation of transitions between different configurations after a deployment event. The optimisation of the actuation sequence is handled by a reinforcement learning model that uses the proximal policy optimisation (PPO) algorithm, by repeatedly trying different sequences and evaluating the impact on classifier performance according to a multi-objective metric.

We build a theoretical range for the Milky Way's (MW) inner dark matter (DM) distribution informed by the FIRE-2, Auriga, VINTERGATAN-GM, and TNG50 simulation suites assuming the canonical cold dark matter (CDM) model. The DM density profiles in Auriga, VINTERGATAN-GM, and TNG50 can be approximately modeled using the adiabatic contraction prescription of Gnedin et al. 2004, while FIRE-2 has stronger baryonic feedback, leading to a departure from the adiabatic contraction model. The simulated halos that are adiabatically contracted are close to spherical (axis ratio $q \in [0.75-0.9]$ at $5^\circ$), whereas halos that experience strong baryonic feedback are oblate ($q \in [0.5-0.7]$). Using the adiabatic contraction and strong baryonic feedback models, along with the observed stellar distribution of the MW, the inner logarithmic density slope for CDM in the MW is predicted to range from $ -0.5$ to $-1.3$. The $J$-factor, which determines the DM-annihilation flux, averaged over a solid angle of $5^\circ$ ($10^\circ$) is predicted to span the range $0.8$-$30$ ($0.6$-$10$) $\times 10^{23} \rm{GeV}^2/\rm{cm}^5$. The $D$-factor, which determines the flux due to DM decay, is predicted to be in the range $0.6$-$2$ ($0.5-1$) $\times10^{23} \rm{GeV}/\rm{cm}^2$. GitHub: The results for this work can be found at this https URL.

We investigate the post-inflationary evolution of a non-minimally coupled inflaton field in scalar-tensor theories, framed within the flexible framework of Einstein-Cartan gravity. By focusing on a class of simplified Higgs-like scenarios, we simulate the transition from the end of inflation to the formation of oscillons using fully-fledged 3+1 classical lattice simulations. Once oscillons form, we extract their profiles and perform 1+1 simulations to evolve their radial equations. Our findings reveal that these oscillons, unlike typical cases in the literature, are relatively short-lived, due to the presence of self interactions at small field values. The radiation produced by this novel type of oscillons can quickly lead to a radiation-dominated Universe, even in the absence of additional fields or interactions. Finally, we leverage our results to derive precise predictions for the inflationary observables and the produced spectrum of gravitational waves associated with oscillon formation. Significantly, our study establishes also an upper bound on the duration of the heating phase in Einstein-Cartan Higgs inflation scenarios.

We develop a novel way to probe subgalactic-scale matter distribution with diffractive lensing on gravitational waves. Five-year observations from Einstein Telescope and DECIGO are expected to probe $k= 10^5\sim 10^8 \,{\rm Mpc}^{-1}$ down to $P(k) = 10^{-16} \sim 10^{-14} \,{\rm Mpc}^3$ level. These results can be interpreted in terms of primordial black holes in the range $M_{\rm PBH} \gtrsim 10^{-3}M_\odot$ down to $f_{\rm PBH} = 10^{-6}$ level, or QCD axion minihalos in the range $m_a = 10^{-3} \sim 10^{-12} \,{\rm eV}$. A key result of the paper is the approximate relation between the scale $k$ and the gravitational wave frequency $f$, derived in an ensemble of `multi-lensing' events. This relation enables direct measurement of the power spectrum at specific scales, with sensitivities characterized by model-independent kernels $\delta P(k)$. Additionally, we delineate the statistical properties of `multi-lensing' based on the `Fresnel number' $N_F$. When $N_F \gtrsim {\cal O}(1)$, the statistical significance can be approximately calculated by Variance of lensing effects, which is directly related to the power spectrum among other moments of matter distribution.

We employ the renormalization group optimized perturbation theory (RGOPT) resummation method to evaluate the equation of state (EoS) for strange ($N_f=2+1$) and non-strange ($N_f=2$) cold quark matter at NLO. This allows us to obtain the mass-radius relation for pure quark stars and compare the results with the predictions from perturbative QCD (pQCD) at NNLO. Choosing the renormalization scale to generate maximum star masses of order $M=2 - 2.6 M_\odot$, we show that the RGOPT can produce mass-radius curves compatible with the masses and radii of some recently observed pulsars, regardless of their strangeness content. The scale values required to produce the desired maximum masses are higher in the strange scenario since the EoS is softer in this case. The possible reasons for such behavior are discussed. Our results also show that, as expected, the RGOPT predictions for the relevant observables are less sensitive to scale variations than those furnished by pQCD.

The detection of a gravitational wave background in the nano-Hertz frequency range from Pulsar Timing Array (PTA) observations offers new insights into evolution of the early universe. In this work we analyze gravitational wave data from PPTA, EPTA, and NANOGrav, as arising from a supercooled first-order phase transition within a hidden sector, characterized by a broken $U(1)_X$ gauge symmetry. Several previous works have discussed challenges in producing observable {PTA signal} from supercooled phases transitions. We discuss these challenges and show how they are overcome by inclusion in part of the proper thermal history of the hidden and the visible sectors. The analysis of this work demonstrates that thermal histories of hidden and visible sectors profoundly influence the gravitational wave power spectrum, an aspect not previously explored in the literature. Further, the analysis of this work suggests that supercooled phase transitions not only align with the Pulsar Timing Array observations but also show promise for gravitational wave detection by future gravitational wave detectors. Our analysis shows that the dominant contribution to the gravitational wave power spectrum for PTA signals comes from bubble collision while the sound wave and turbulence contributions are highly suppressed. It is also found that all the PTA events are of detonation type while deflagration and hybrid events are absent. The analysis presented in this work provides a robust framework for further investigations on the origin of gravitational wave power spectrum in the early universe and for their experimental observation in the future.

In this paper, cosmic distance duality relation is probed without considering any background cosmological model. The only \textit{a priori} assumption is that the Universe is described by the Friedmann-Lema$\hat{i}$tre-Robertson-Walker (FLRW) metric The strong gravitational lensing (SGL) data is used to construct the dimensionless co-moving distance function $d(z)$ and latest type Ia supernovae (SNe Ia) Pantheon+ data is used to estimate luminosity distances at the corresponding redshifts $z$. Using the distance sum rule along null geodesics of the FLRW metric, the CDDR violation is probed in both flat and non-flat space-time by considering two parametrizations for $\eta(z)$, the function generally used to probe the possible deviations from CDDR. The results show that, CDDR is compatible with the observations at a very high level of confidence for linear parametrization in flat Universe. In non-flat Universe too, CDDR is valid within $1\sigma$ confidence interval with a mild dependence of $\eta$ on the curvature density parameter $\Omega_{K}$ . The results for non-linear parametrization also show no significant deviation from CDDR.

For several independent reasons, the idea that notorious sources of entropy could exist in the Universe has been recently revived. Taking advantage of a new framework accounting for non-equilibrium processes in cosmology, we explicitly investigate the cosmological dynamics as a function of the entropy production, focusing on the stability of the system. An exhaustive investigation is performed. As the main physical conclusion, we show that for a wide class of entropy source terms, the fluid dynamics converges towards an effective cosmological constant. Constraints on the associated entropic force are also obtained.

We present a physics-informed autoencoder designed to encode the equation of state of neutron stars into an interpretable latent space. In particular the input will be encoded in the mass, radius, and tidal deformability values of a neutron star. Unlike traditional black-box models, our approach incorporates additional loss functions to enforce explainability in the encoded representations. This method enhances the transparency of machine learning models in physics, providing a robust proof-of-concept tool to study compact stars data. Our results demonstrate that the proposed autoencoder not only accurately estimates the EoS parameters and central density/pressure but also offers insights into the physical connection between equation of state and observable physical quantities. This framework conceptualizes the physical differential equations themselves as the ``encoders", allowing interpretability of the latent space.

The Vacuum tunneling rate $\Gamma$ from the effective action is a key to studying the cosmological first-order phase transition(FOPT). One solid way to compute the $\Gamma$ is to start with the derivative expansion of the effective action and solve the bounce equation numerically. In this process, the renormalization factor $Z$ of the tunneling field may play a center rule, which is not considered in existing packages. Therefore, we present a \texttt{Mathematica} package \vt to compute the bounce action with or without the renormalization factor. Applying the \vt package, we find that the presence of $Z$ has a significant impact on the action, as well on the tunneling path. We provide some concrete examples to demonstrate the difference between the solution with and without the renormalization factor, both in the action and tunneling path. This package is based on the modified shooting and path deformation method. We also made some optimizations for the super-cooling phase transition(thick wall scenario), in which other numerical package works poorly. This package works as long as the expressions can give values of the potential and the renormalization at a certain field point. This means the input potential and the renormalization can be merely numerical quantities without analytical expressions. The computation time can be as short as 1 second in single-field tunneling and several seconds in multi-field cases.

In the Arkani-Hamed-Dimopoulos-Dvali (ADD) model with n extra compactified dimensions, the gravitational potential scales as 1/r^{n+1} and becomes significantly stronger at short distances. We investigate the possibility of forming small-sized composite systems of Standard Model particles bound by this potential. Such bound states, composed of quarks, neutrinos, axions, or other particles, exhibit a small cross-section-to-mass ratio, making them viable candidates for dark matter.

We investigate analogue gravity phenomena arising as a result of the linear perturbation of the spherically symmetric accretion flows onto non rotating black holes, where the gravitational field is determined by a set of post Newtonian pseudo Schwarzschild black hole potentials and the infaling matter is described by a relativistic multi-species equation of state. The stationary transonic integral accretion solutions corresponding to the steady state of aforementioned type of accreting systems are constructed and the stability analysis of such solutions are performed through the time dependent linear perturbation of the accretion flow. Such linear stability analysis leads to the formation of a black hole like sonic metric embedded within the infalling matter. The acoustic horizons are then identified by constructing the causal structure, i.e., the Carter-Penrose diagrams. The variation of the analogue surface gravity corresponding to the aforementioned sonic metric has been studied as a function of various parameters governing the accretion flow.

The potential hadron-to-quark phase transition in neutron stars has not been fully understood as the property of cold, dense, and strongly interacting matter cannot be theoretically described by the first-principle perturbative calculations, nor have they been systematically measured through terrestrial low-to-intermediate energy heavy-ion experiments. Given the Tolman--Oppenheimer--Volkoff (TOV) equations, the equation of state (EoS) of the neutron star (NS) matter can be constrained by the observations of NS mass, radius, and tidal deformability. However, large observational uncertainties and the limited number of observations currently make it challenging to strictly reconstruct the EoS, especially to identify interesting features such as a strong first-order phase transition. In this work, we study the dependency of reconstruction quality of the phase transition on the number of NS observations of mass and radius as well as their uncertainty, based on a fiducial EoS. We conquer this challenging problem by constructing a neural network, which allows one to parameterize the EoS with minimum model-dependency, and by devising an algorithm of parameter optimization based on the analytical linear response analysis of the TOV equations. This work may pave the way for the understanding of the phase transition features in NSs using future $X$-ray and gravitational wave measurements.

This article is dedicated to the memory of Alexei Starobinsky. I begin with some recollections of him and then review the generalization of his wonderful stochastic formalism from scalar potential models to theories which interact with fermions and photons, and finally to theories with derivative interactions such as nonlinear sigma models and gravity. This entails effective potentials generated by the usual field-dependent masses, as well as by field-dependent field strengths, and by field-dependent Hubble parameters. I also discuss secular loop corrections which cannot be captured by stochastic techniques.

Axion-like particles (ALPs), hypothetical extensions of the Standard Model, can convert into photons in an external magnetic field. Two recent studies~\cite{Todarello:2023ptf, Ruz:2024gkl} explored the phenomenology of ALP-photon conversion in the magnetic field of the solar atmosphere. Dark matter ALPs convert into radio photons and may be detected with next-generation radio interferometers, while ALPs produced in the solar core convert into X-rays. Thanks to solar observation acquired with the NuSTAR X-ray telescope, Ref.~\cite{Ruz:2024gkl} establishes stringent robust bounds on the ALP-photon coupling over a large portion of parameter space.

Connecting inflation with neutrino physics through non-thermal leptogenesis via direct inflaton-right-handed neutrino (RHN) coupling naturally incorporates neutrino reheating, leaving no ambiguity regarding the early history of the Universe. In ref.~\cite{Zhang:2023oyo}, we demonstrate that non-thermal leptogenesis from inflaton decay expands the viable parameter space compared to thermal leptogenesis and provides a natural link to inflation. In this work, we refine our previous findings by closely examining the dynamics of neutrino reheating. We first calculate the duration of neutrino reheating on a general basis, then analyze inflationary observables consistent with neutrino reheating across four models, establishing a direct connection between baryon asymmetry and the spectral index. This approach places these two important observables on the same plane and yields specific predictions that help break the degeneracy among inflationary models. The well-motivated and economical framework offers a simple, natural, and testable description of the early Universe.

We present PRINCESS, a computational tool designed to predict gravitational wave observations from compact binary coalescences (CBCs) in current and future detector networks. PRINCESS uniquely combines predictions of both individual gravitational wave events and the associated astrophysical background, leveraging user-provided CBC catalogs. With the anticipated improvements in detector sensitivity from second-generation (2G) to third-generation (3G) observatories like the Einstein Telescope and Cosmic Explorer, the tool aims to constrain models of stellar formation and compact object evolution. PRINCESS calculates the signal-to-noise ratio (SNR) for individual events and predicts the stochastic background arising from unresolved sources. We detail the code structure, installation, and usage, providing examples of predictions for different astrophysical models and detector configurations. Results include forecasts for binary black hole detections and background spectra, highlighting the potential of future networks to resolve nearly all CBC events. We also discuss ongoing developments to expand PRINCESS's capabilities, including accounting for additional sources such as extreme mass ratio inspirals (EMRIs) and incorporating more advanced detector models.