Papers for Tuesday, May 13 2025

We present the ALMA-CRISTAL survey, an ALMA Cycle 8 Large Program designed to investigate the physical properties of star-forming galaxies at $4 \lesssim z \lesssim 6$ through spatially resolved, multi-wavelength observations. This survey targets 19 star-forming main-sequence galaxies selected from the ALPINE survey, using ALMA Band 7 observations to study [CII] 158 $\mu$m line emission and dust continuum, complemented by JWST/NIRCam and HST imaging to map stellar and UV emission. The CRISTAL sample expanded to 39 after including newly detected galaxies in the CRISTAL fields, archival data, and pilot study targets. The resulting dataset provides a detailed view of gas, dust, and stellar structures on kiloparsec scales at the end of the era of reionization. The survey reveals diverse morphologies and kinematics, including rotating disks, merging systems, [CII] emission tails from potential interactions, and clumpy star formation. Notably, the [CII] emission in many cases extends beyond the stellar light seen in HST and JWST imaging. Scientific highlights include CRISTAL-10, exhibiting an extreme [CII] deficit similar to Arp 220; and CRISTAL-13, where feedback from young star-forming clumps likely causes an offset between the stellar clumps and the peaks of [CII] emission. CRISTAL galaxies exhibit global [CII]/FIR ratios that decrease with increasing FIR luminosity, similar to trends seen in local galaxies but shifted to higher luminosities, likely due to their higher molecular gas content. CRISTAL galaxies also span a previously unexplored range of global FIR surface brightness at high-redshift, showing that high-redshift galaxies can have elevated [CII]/FIR ratios. These elevated ratios are likely influenced by factors such as lower metallicity gas, the presence of significant extraplanar gas, and contributions from shock-excited gas.

Non-spherical potentials allow a wide range of trajectories, both regular and chaotic, whose periapse distances can vary orbit to orbit. In particular chaotic trajectories can bring a system arbitrarily close to the central massive black hole leading to a disruption. In this paper, we work with an observationally benchmarked model of the innermost 200 pc of the Milky Way and show that low z-angular momentum trajectories are commonly chaotic. We compute the timescales and properties of close pericenter passages, and compare the implied collisionless disruption rate to the well-studied collisional rate from 2-body scatterings. We find that the relative collisionless rate can dominate by orders of magnitude. Our calculations are relevant for a wide range of disruption phenomena, including the production of hypervelocity stars (HVSs) and tidal disruption events (TDEs). Most of these disruptions involve stars come from the Nuclear Stellar Cluster, with a pericenter distribution that strongly favours shallow encounters, and a preference for high inclination interactions. The latter implies that unbound disrupted material - whether ejected stars or stellar debris - would be preferentially directed towards the galactic poles. Many of our conclusions apply generally to any galaxy with a non-spherical galactic centre potential and central massive black hole.

The formation of compact high-redshift star-forming clumps, the physical processes driving their evolution and their potential connection to present-day Globular Clusters are key open questions in galaxy formation. In this work, we aim to shed light on these aspects using the SImulating the Environment where Globular clusters Emerged (SIEGE) project, a suite of cosmological zoom-in simulations with sub-parsec resolution specifically designed to investigate the physical conditions behind the origin of compact stellar systems in high-redshift environments. The simulation object of this study focuses on a dwarf galaxy with a virial mass of a few $10^9$ $M_\odot$ at $z=6.14$, where the spatial resolution reaches 0.3 pc $h^{-1}$. Individual stars are formed directly by sampling the initial mass function with a 100\% star formation efficiency, a setup designed to explore the impact of a high star formation efficiency under high-redshift conditions. The simulation reveals the emergence of numerous stellar clumps with sizes of 1-3 pc, stellar surface densities up to almost $10^4$ $M_\odot$ pc$^{-2}$, and masses predominantly spanning from $10^3$ $M_\odot$ to several $10^4$ $M_\odot$, with a few reaching $10^5$ $M_\odot$ and up to $10^6$ $M_\odot$. All clumps form during intense, short bursts of star formation lasting less than a Myr, often with negligible dark matter content (dark-to-stellar mass ratios below 1 within three times their effective radii). We measure a clear correlation between mass and size, and a clump mass function described by a power-law with a slope of -2. Star formation conditions in the simulation behave similarly to those of a feedback-free starburst scenario, where dense clumps form due to inefficient stellar feedback over small timescales. Notably, some clumps exhibit properties closely resembling those of present-day globular clusters.

The steady-state gamma-ray emission from the Sun is thought to consist of two emission components due to interactions with Galactic cosmic rays: (1) a hadronic component covering the solar disk, and (2) a leptonic component peaking at the solar edge and extending into the heliosphere. The flux of these components is expected to vary with the 11-year solar cycle, being highest during solar minimum and lowest during solar maximum, because it is correlated with the cosmic-ray flux. No study has yet analyzed the flux variation of the two components separately over solar cycles. In this work, we measure the temporal variations of the flux of each component over 15 years of Fermi Large Area Telescope observations and compare them with the sunspot number and Galactic cosmic-ray flux from AMS-02 near the Earth. We find that the flux variation of the disk anticorrelates with solar activity and correlates with cosmic-ray protons, confirming its emission mechanism. The flux variation of the extended component anticorrelates with solar activity only until mid 2012. After that, we no longer observe any correlation or anticorrelation, even with the CR electron flux. This most likely suggests that cosmic-ray transport and modulation in the inner heliosphere are unexpectedly complex and different for electrons and protons or, alternatively, the presence of an additional, unknown component of gamma rays or cosmic rays. These findings impact space weather research and emphasize the need for close monitoring of Cycle 25 and the ongoing polarity reversal.

JWST observations uncovered a large number of massive quiescent galaxies (MQGs) at z>3, which theoretical models struggle to reproduce. Explaining the number density of such objects requires extremely high conversion efficiency of baryons into stars in early dark matter halos. Using stellar kinematics, we can investigate the processes shaping the mass assembly histories of MQGs. We present high-resolution JWST/NIRSpec integral field spectroscopy of GS-9209, a massive, compact quiescent galaxy at z=4.66 ($\log \left (M_{\ast}/M_{\odot} \right) = 10.52 \pm 0.06 $, $R_{eff} = 220 \pm 20$ pc). Full spectral fitting of the spatially resolved stellar continuum reveals a clear rotational pattern, yielding a spin parameter of $\lambda_{R_{eff}} = 0.65 \pm 0.12$. With its high degree of rotational support, this galaxy challenges the scenario of MQGs growing mainly by dry major mergers. This study suggests that at least a fraction of the earliest quiescent galaxies were fast rotators and that quenching was dynamically gentle process, preserving the stellar disc even in highly compact objects. Using Jeans anisotropic modelling (JAM) and a NFW profile, we measure a dark matter fraction of $f_{\rm DM} \left (< R_{eff} \right ) = 6.3^{+2.8}_{-1.7}%$, which is plausible given that this galaxy is extremely compact. Our findings use kinematics to independently confirm the massive nature of early quiescent galaxies, previously inferred from stellar population modelling. We suggest that GS-9209 has a similar structure to low-redshift 'relic' galaxies. However, unlike relic galaxies which have bottom-heavy initial mass functions (IMF), the dynamically inferred stellar mass-to-light ratio of GS-9209 is consistent with a Milky-Way like IMF. The kinematical properties of GS-9209 are different from those of z<1 early-type galaxies and more similar to those of recently quenched post-starburst galaxies at z>2.

During the epoch of reionization (EoR), the first generation of luminous sources in our Universe emitted ionizing photons that almost completely ionized the gas in the intergalactic medium (IGM). The growth of ionized bubbles and the persistence of neutral islands within the IGM hold vital clues to understanding the morphology and timeline of cosmic reionization. We explore the potential of photometric IGM tomography using deep narrow-band (NB) imaging to observe the Lyman-$\alpha$ forest transmission in background galaxies with the Subaru/Hyper-Suprime Cam (HSC). Based on our simulations, we find that the currently available NB filter is suitable for mapping the IGM at $z\simeq 5.7$, corresponding to the late stages of reionization. Our findings indicate that over $\sim$500 background galaxies are needed to accurately reconstruct the IGM at scales greater than 200 Mpc, achieving more than 40 per cent correlation with the true distribution. This technique can help detect final remaining neutral islands that span more than 20 Mpc lengths. Using the superpixel method built to identify physical patterns in noisy image data, we find that the neutral island size distribution can be recovered with an accuracy of $\sim$0.3 dex. Furthermore, we demonstrate that these reconstructed maps are correlated with the galaxy distribution and anti-correlated with the cosmological 21-cm signal from neutral hydrogen in the IGM. Lastly, we find that these reconstructed maps are anti-correlated with the patchy optical depth to the cosmic microwave background. As such, multiple measurements can be employed for confirmed detection of neutral islands during the end stages.

We investigate the scaling relation between the Compton Y parameter and mass in a gravity-selected sample of galaxy clusters, selected based on their gravitational lensing effects on background galaxies. Unlike ICM-selected samples, gravity-selected clusters do not require corrections for selection biases in intracluster medium (ICM) properties, such as Compton Y, since these properties are not used for selection. Our sample consists of 13 gravity-selected clusters with weak-lensing signal-to-noise ratios greater than 7 and redshifts in the range $0.1 < z < 0.4$, drawn from the peak catalog of shear-selected clusters in Oguri et al. (2021). We determined spectroscopic redshifts using SDSS spectroscopy, derived tangential radial profiles from HSC shear data, and calculated cluster masses by fitting a Navarro, Frenk, \& White (1997) profile. Compton Y measurements were obtained from ACT maps. Our sample reveals a broader scatter in the Compton Y-mass scaling relation compared to ICM-selected samples, even after accounting for selection effects and mass bias of the latter, that we find to be $(1-b) = 0.72 \pm 0.09$. Additionally, we identify a second population of clusters with unusually small Compton Y for their mass, which is absent in ICM-selected samples. This second population comprises $13^{+10}_{-8}$\% of the entire sample, with O19 currently being the only secure representative.

Young (<500 Myr) multi-planet transiting systems are valuable environments for understanding planet evolution by offering an opportunity to make direct comparisons between planets from the same formation conditions. TOI-2076 is known to harbor three, 2.5-4 $R_\oplus$ planets on 10-35 day orbits. All three are JWST cycle 3 targets (for transmission spectroscopy). Here, we present the detection of TOI-2076 e; a smaller (1.35 $R_\oplus$), inner (3.02 day) planet in the system. We update the age of the system by analyzing the rotation periods, Lithium equivalent widths, color-magnitude diagram, and variability of likely co-moving stars, finding that TOI-2076 and co-moving planetary system TOI-1807 are 210 $\pm$ 20 Myr. The discovery of TOI-2076 e is motivation to revisit known transiting systems in search of additional planets that are now detectable with new TESS data and updated search methods.

JWST spectroscopy has built large emission line samples at $z\gtrsim4$, but it has yet to confidently reveal many galaxies with the hard radiation fields commonly associated with AGN photoionization. While this may indicate a weaker UV ionizing spectrum in many $z>4$ AGNs or obscuration from dense neutral gas and dust, the complete picture remains unclear owing to the small number of deep rest-UV spectra. Here we characterize the strength of high ionization lines in $53$ new galaxies observed with NIRSpec $R=2700$ grating spectroscopy. We present new detections of narrow NV$\lambda1240$ in two galaxies. One is a previously-confirmed $z=6.98$ Little Red Dot (LRD) with broad H$\beta$, and the other is a $z=8.72$ galaxy with a narrow line spectrum. Neither source exhibits CIV or HeII emission, indicating large NV/CIV and NV/HeII ratios that may reflect a combination of nitrogen-enhancement and resonant scattering effects. We investigate the incidence of narrow high ionization lines in a large database of $851$ NIRSpec grating spectra, and we separately quantify the fraction of LRDs with narrow high ionization UV emission lines. Our results likely suggest that hard radiation fields are indeed present in a small subset of LRDs ($12.5^{+23.7}_{-10.4}\%$) and UV-selected galaxies ($2.2^{+1.7}_{-1.0}\%$) at $z>4$. The identification of narrow high ionization lines in the population of LRDs with strong Balmer absorption suggests the dense neutral hydrogen gas may not uniformly cover the nucleus. The NV/CIV and NV/HeII ratios suggest that efforts to identify high ionization lines should extend down in wavelength to the NV doublet.

To explore the hypothesis of a common source of variability in two time series, observers may estimate the magnitude-squared coherence (MSC), which is a frequency-domain view of the cross correlation. For time series that do not have uniform observing cadence, MSC can be estimated using Welch's overlapping segment averaging. However, multitaper has superior statistical properties to Welch's method in terms of the tradeoff between bias, variance, and bandwidth. The classical multitaper technique has recently been extended to accommodate time series with underlying uniform observing cadence from which some observations are missing. This situation is common for solar and geomagnetic datasets, which may have gaps due to breaks in satellite coverage, instrument downtime, or poor observing conditions. We demonstrate the scientific use of missing-data multitaper magnitude-squared coherence by detecting known solar mid-term oscillations in simultaneous, missing-data time series of solar Lyman $\alpha$ flux and geomagnetic Disturbance Storm Time index. Due to their superior statistical properties, we recommend that multitaper methods be used for all heliospheric time series with underlying uniform observing cadence.

The instantaneous emission from a relativistic surface endowed with a Lorentz factor $\Gamma$ that decreases away from the outflow symmetry axis can naturally explain the three phases observed by Swift/XRT in GRBs and their afterglows (GRB tail, afterglow plateau and post-plateau). We expand the analytical formalism of the "Larger-Angle Emission" model previously developed for "Power-Law" outflows to "n-Exponential" outflows (e.g. exponential with $n=1$ and Gaussian with $n=2$) and compare their abilities to account for the X-ray emission of XRT afterglows. We assume power-law $\Gamma$-dependences of two spectral characteristics (peak-energy and peak intensity) and find that, unlike Power-Law outflows, n-Exponential outflows cannot account for plateaus with a temporal dynamical range larger than 100. To include all information existing in the Swift/XRT measurements of X-ray aferglows (0.3-10 keV unabsorbed flux and effective spectral slope), we calculate 0.3 keV and 10 keV light-curves using a broken power-law emission spectrum of peak-energy and low-and high-energy slopes that are derived from the effective slope measured by XRT. This economical peak-energy determination is found to be consistent with more expensive spectral fits. The angular distributions of the Lorentz factor, comoving frame peak-energy, and peak-intensity ($\Gamma (\theta), E'_p (\theta), i'_p(\theta)$) constrain the (yet-to-be determined) convolution of various features of the production of relativistic jets by solar-mass black-holes and of their propagation through the progenitor/circumburst medium, while the $E'_p (\Gamma)$ and $i'_p (\Gamma)$ dependences may constrain the GRB dissipation mechanism and the GRB emission process.

We present the public release of EXP, a basis function expansion C++ library and Python package for running N-body galactic simulations and dynamical discovery. EXP grew out of the need for methodology that seamlessly connects theoretical descriptions of dynamics, N-body simulations, and compact descriptions of observed data. EXP decomposes a galaxy into multiple bases for a variety of scales and geometries and is thus able to represent arbitrarily complex simulations with many components (e.g., disk, bulge, dark matter halo, satellites). pyEXP provides a full Python interface to the EXP C++ libraries. Example workflows based on previously published work are available and distributed as accompanying examples and tutorials. The examples and tutorials flatten the learning curve for adopting basis function expansion tools to generate and analyze the significance of coefficients and discover dynamical relationships using time series analysis. The current release uses a powerful, non-parametric technique for time series analysis that decomposes basis function expansion coefficients into interpretable components without assuming a specific model structure. This enables the automated discovery of dynamical signals from simulations. We provide a full online manual hosted by Read the Docs.

We study the apparent optical brightness of satellite streaks in ground-based observatory cameras as the orbit height is varied. This is done via simulations. Several factors are in play: the range to the satellite, satellite size and geometry, the diameter of the telescope aperture, and the angular velocity of the satellite image as it streaks across the camera's focal plane. For a large telescope aperture, a satellite in a lower orbit becomes out-of-focus in a camera focused at infinity. In the case of Rubin Observatory's LSSTCam, we find that these factors nearly cancel as a large satellite is moved from an orbit at 550 km down to 350 km.

The relation between Faraday Rotation Measure (RM) and differential Faraday rotation by unresolved structure of a turbulent plasma is investigated for extragalactic radio sources. The RM scatter of a sample of sources affected by turbulent Faraday screens with identical power spectra of Faraday depth structure, is referred to as RM jitter. For fixed amplitude and slope of the power spectrum, the range of possible RMs depends on the wavelength coverage of the survey. RM jitter is independent of Faraday depth resolution as it results from the true Faraday depth dispersion and effects of wavelength-dependent depolarization. RM jitter for a flux density limited sample is sensitive to the power law index gamma of the power spectrum of Faraday depth structure. Assuming depolarization by a turbulent Faraday screen for all sources, a simulated flux-density-limited sample can reproduce the high RM scatter found by Vanderwoude et al. (2024) for sources that are less than 3% polarized. RM jitter of sources that are more than 3% polarized, is found to be smaller than the observed scatter, indicating that plasma other than the near-source environment dominates the RM scatter for the more polarized sources. The significance of RM jitter for applications of the RM grid is discussed.

AR Scorpii, the so called white dwarf pulsar, contains a rapidly rotating magnetic white dwarf (WD; Pspin = 117.0564 s) interacting with a cool, red dwarf (RD) companion in a 3.56 hour orbit. It is a strong radio source with an inverted spectral index between 1-200 GHz that is indicative of synchrotron emission. This paper presents the first submillimeter observations of AR Scorpii using the Submillimeter Array, helping to fill a critical gap in the spectral energy distribution between 10-600 GHz. The average flux densities at 220 and 345 GHz are 124 and 86 mJy, respectively. The lower than expected flux density at 345 GHz suggests a break in the synchrotron emission at about 200 GHz. A periodogram analysis of the 220 GHz observations shows a modulation with an amplitude of ~6% at a period of 58.26 s or at twice the spin frequency of the white dwarf. This modulation is the first direct detection of the WD spin period at radio frequencies and implies that the synchrotron emission arises near the WD and and not from an interaction with the photosphere of the RD. A fit to the spectral energy distribution shows that the synchrotron emission arises from a small, low density region with a magnetic field of 43 MG at a distance of 0.6 orbital radii from the WD. This result implies that AR Scorpii contains a weakly magnetic WD (~15 MG) and not a strongly magnetic WD (~500 MG) as previously asserted.

Object LAMOST J020623.21+494127.9 (program star) in the thin disk of the Milky Way (MW) is reported as a highly r-process-enhanced (RPE) r-II star with[Eu/Fe]= +1.32 and [Fe/H]= -0.54. The chemical profile of the star reflects the intrinsic composition of the gas cloud present at its birth. Using an abundance decomposition method, we fit 25 elements from the abundance dataset, including 10 heavy neutron capture elements. We explore the astrophysical origin of the elements in this star through its abundance ratios and component ratios. We find that the contributions from the massive stars played a significant role in the production of light elements in the program star. Our analysis reveals that the heavy neutron-capture elements are produced purely by the main r-process. However, the adopted main r-process model does not adequately fit the observed data, suggesting another main r-process pattern may exist.

The properties of X-ray flashes (XRFs) and X-ray rich gamma-ray bursts (XRRs) as compared with classical gamma-ray bursts (C-GRBs) have been widely discussed during the \emph{Swift} era. It has been proposed that XRFs and XRRs are low-energy extensions of the GRB population so that they should follow similar correlations. To further examine this idea, we collect a sample of $303$ GRBs detected by \emph{Swift} over the past two decades, all of which have reliable redshifts and spectral parameters. The bursts are classified into XRFs, XRRs, and C-GRBs based on their softness ratio (SR), which is calculated by dividing the $25-50$ keV fluence with the $50-100$ keV fluence. A strong correlation is found among the isotropic energy $E_{\mathrm{iso}}$, peak luminosity $L_{\mathrm{p}}$, and rest frame burst duration $T_{90, \mathrm{rest}}$, i.e., $E_{\mathrm{iso}} \propto L_{\mathrm{p}}^{0.88\pm0.02} T_{90, \mathrm{rest}}^{0.58\pm0.02}$. Additionally, two tight three-parameter correlations involving SR and the rest-frame peak energy $E_{\mathrm{p}}$ are also derived, i.e. $E_{\mathrm{p}} \propto E_{\mathrm{iso}}^{0.20\pm0.02} \mathrm{SR}^{-2.27\pm0.15}$ and $E_{\mathrm{p}} \propto L_{\mathrm{p}}^{0.17\pm0.02} \mathrm{SR}^{-2.33\pm0.14}$. It is interesting to note that XRFs, XRRs, and C-GRBs all follow the above correlations. The implications of these correlations and their potential application in cosmology are discussed.

We report a detailed spectroscopic study of the gas dynamics and hydrostatic mass bias of the galaxy cluster Abell 2029, utilizing high-resolution observations from XRISM Resolve. Abell 2029, known for its cool core and relaxed X-ray morphology, provides an excellent opportunity to investigate the influence of gas motions beyond the central region. Expanding upon prior studies that revealed low turbulence and bulk motions within the core, our analysis covers regions out to the scale radius $R_{2500}$ (670~kpc) based on three radial pointings extending from the cluster center toward the northern side. We obtain accurate measurements of bulk and turbulent velocities along the line of sight. The results indicate that non-thermal pressure accounts for no more than 2% of the total pressure at all radii, with a gradual decrease outward. The observed radial trend differs from many numerical simulations, which often predict an increase in non-thermal pressure fraction at larger radii. These findings suggest that deviations from hydrostatic equilibrium are small, leading to a hydrostatic mass bias of around 2% across the observed area.

We investigate the relativistic jet of the powerful radio-emitting blazar J1429+5406 at redshift z=3.015. Our understanding of jet kinematics in z>3 quasars is still rather limited, based on a sample of less than about 50 objects. The blazar J1429+5406 was observed at a high angular resolution using the method of very long baseline interferometry over more than two decades, between 1994 and 2018. These observations were conducted at five radio frequencies, covering a wide range from 1.7 to 15 GHz. The outer jet components at ~20-40 milliarcsecond (mas) separations from the core do not show discernible apparent motion. On the other hand, three jet components within the central 10 mas region exhibit significant proper motion in the range of (0.045-0.16)mas/year, including one that is among the fastest-moving jet components at z>3 known to date. Based on the proper motion of the innermost jet component and the measured brightness temperature of the core, we estimated the Doppler factor, the bulk Lorentz factor, and the inclination angle of the jet with respect to the line of sight. The core brightness temperature is at least 3.6 x 10^{11} K, well exceeding the equipartition limit, indicating Doppler-boosted radio emission. The low jet inclination (<5.4 deg) firmly places J1429+5406 into the blazar category.

The bulk silicate Earth (BSE) is depleted in moderately volatile elements, indicating Earth formed from a mixture of volatile-rich and -poor materials. To better constrain the origin and nature of Earth's volatile-rich building blocks, we determined the mass-dependent isotope compositions of Ge in carbonaceous (CC) and enstatite chondrites. We find that, similar to other moderately volatile elements, the Ge isotope variations among the chondrites reflect mixing between volatile-rich, isotopically heavy matrix and volatile-poor, isotopically light chondrules. The Ge isotope composition of the BSE is within the chondritic range and can be accounted for as a ~2:1 mixture of CI and enstatite chondrite-derived Ge. This mixing ratio appears to be distinct from the ~1:2 ratio inferred for Zn, reflecting the different geochemical behavior of Ge (siderophile) and Zn (lithophile), and suggesting the late-stage addition of volatile-rich CC materials to Earth. On dynamical grounds it has been argued that Earth accreted CC material through a few Moon-sized embryos, in which case the Ge isotope results imply that these objects were volatile-rich, presumably because they were either undifferentiated or accreted volatile-rich objects themselves before being accreted by Earth.

The Gaia-Sausage-Enceladus was the last major merger and central turning point in the Milky Way's story. This event, comparable in mass to the Large Magellanic Cloud today, left behind significant debris that provides valuable insights into the assembly history of our Galaxy and the chemical evolution of dwarf galaxies. By examining the aftermath of the GSE merger, we can delve deeper into understanding how the Milky Way's formation unfolded and how dwarf galaxies evolved chemically. Specifically, the distinct patterns of neutron capture elements such as Eu and Ba, along with Mg, offer clues about the star formation history. Through a comprehensive analysis of data compiled in the SAGA database, we investigated the Gaia Sausage-Enceladus' star formation history. Elemental abundance ratios ([Eu/Mg], [Ba/Mg], and [Eu/Ba]) derived from this study, when compared with those of surviving Milky Way satellites, indicate that the GSE experienced a prolonged period of slow star formation, lasting over 2 Gyr, until it was eventually quenched by merging with the Milky Way. Consequently, these elemental signatures serve as a unique window into the complex history of both surviving and accreted satellites orbiting our Galaxy.

We present a systematic analysis of radio supernovae (SNe) to investigate the statistical tendencies of SN progenitors' mass-loss rates and shock acceleration efficiencies. We conduct parameter estimation through Markov chain Monte Carlo (MCMC) analysis for 32 radio SN samples with a clear peak observed in their light curves, and successfully fit 27 objects with the widely-used radio SN model. We find the inferred mass-loss rates of stripped-envelope SN progenitors are by an order of magnitude greater ($\sim 10^{-3}\,M_\odot{\rm yr}^{-1}$) than those of SN II progenitors ($<10^{-4}\,M_\odot{\rm yr}^{-1}$). The efficiencies of electron acceleration and magnetic field amplification are found to be less than $10^{-2}$, and the possibility of their energy equipartition is not ruled out. On the other hand, we find the following two properties that might be related to limitations of the standard model for radio SNe; one is the extremely high magnetic field amplification efficiency, and the other is the shallower density gradient of the outer ejecta. We suggest the new interpretation that these peculiar results are misleading due to the setup that is not included in our model, and we identify the missing setup as a dense CSM in the vicinity of the progenitor. This means that a large fraction of radio SN progenitors might possess dense CSM in the vicinity of the progenitor, which is not smoothly connected with the outer CSM.

Spectropolarimetric observations show that many low-mass stars possess large-scale poloidal magnetic fields with considerable dipole component, which in some cases exhibit temporal dynamics - cycles or reversals. Although it is widely accepted that their magnetic fields are generated by the dynamo process, it is hard to reproduce coherent oscillations of large-scale magnetic fields with a dipolar symmetry as observed for the Sun when turbulent and compressible regimes are explored. Most previous 3D numerical studies partially avoided this problem by considering a numerical domain with low density stratification, which may correspond to neglecting surface effects where density drops considerably. To address this question, we perform direct numerical simulations of convective dynamos in extreme parameter regimes of both strong turbulence and strong density stratification, using software MagIC. Our simulations exhibit rotationally-influenced large-scale convective motions surrounded by a turbulent compressible surface layer. We find complex time variations of the magnetic field in flow regimes of predominantly dipolar configuration with respect to the few large-scale harmonics. In such regimes, turbulent surface layer induces global magnetic pumping mechanism, transporting magnetic energy into the deep interiors of our dynamo model. Dipole magnetic fields are found in regimes of transition between solar- and anti-solar differential rotation, and interact dynamically with it. The spatial distribution and temporal behavior of the large-scale fields is consistent with observations of low-mass stars, which suggest magnetic pumping could promote time-dependent magnetic fields with a dipolar symmetry as observed for the Sun and other solar-like stars. Our results suggest a parameter path in which dynamo models with a complex multiscale dynamics should be explored.

This study investigates the variability of the theoretical correction factor, $f_{\Delta \nu}$, used in red giant branch (RGB) scaling relations, arising from different assumptions in stellar model computations. Adopting a commonly used framework, we focused on a 1.0 $M_{\odot}$ star and systematically varied seven input parameters: the reference solar mixture, the initial helium abundance, the inclusion of microscopic diffusion and mass loss, the method for calculating atmospheric opacity, the mixing-length parameter, and the boundary conditions. Each parameter was tested using two distinct but physically plausible values to mimic possible choices of different evolutionary codes. For each resulting stellar model, we computed the oscillation frequencies along the RGB and derived the large frequency spacing, $\Delta \nu_0$. The correction factor $f_{\Delta \nu}$ was then calculated by comparing the derived $\Delta \nu_0$ with that predicted by the uncorrected scaling relations. We found substantial variability in $f_{\Delta \nu}$ across the different models. The variation ranged from approximately 1.3% in the lower RGB to about 3% at $\log g = 1.4$. This level of variability is significant, as it corresponds to roughly half the values typically quoted in the literature and leads to a systematic change in derived masses from 5% to more than 10%. The most significant contribution to this variability came from the choice of atmospheric opacity calculation (approximately 1.2%), with a smaller contribution from the inclusion of microscopic diffusion (approximately 0.4%). These results indicate that the choice of the reference stellar model has a non-negligible impact on the calculation of correction factors applied to RGB star scaling relations.

Narrow eccentric planetary ringlets have sharp edges, sizable eccentricity gradients, and a confinement mechanism that prevents radial spreading due to ring viscosity. Most proposed ringlet confinement mechanisms presume that there are one or more shepherd satellites whose gravitational perturbations keeps the ringlet confined radially, but the absence of such shepherds in Cassini observations of Saturn's rings casts doubt upon those ringlet confinement mechanisms. The following uses a suite of N-body simulations to explore an alternate scenario, whereby ringlet self-gravity drives a narrow eccentric ringlet into a self-confining state. These simulations show that, under a wide variety of initial conditions, an eccentric ringlet's secular perturbations of itself causes the eccentricity of its outer edge to grow at the expense of its inner edge. This causes the ringlet's nonlinearity parameter $q$ to grow over time until it exceeds the $q\simeq\sqrt{3}/2$ threshold where the ringlet's orbit-averaged angular momentum flux due to viscosity + self-gravity is zero. The absence of any net radial angular momentum transfer through the ringlet means that the ringlet has settled into a self-confining state, i.e. it does not spread radially due to its viscosity, and simulations also show that such ringlets have sharp edges. Nonetheless, viscosity still circularizes the ringlet in time $\tau_e\sim10^6$ orbits $\sim1000$ years, which will cause the ringlet's nonlinearity parameter to shrink below the $q\simeq\sqrt{3}/2$ threshold and allows radial spreading to resume. Either sharp-edged narrow eccentric ringlets are transient phenomena, or exterior perturbations are also sustaining the ringlet's eccentricity. We then speculate about how such ringlets might come to be.

We analyze the velocity field of peripheral members of the Local Group. The Hubble flow at distances from 400 to 1400~kpc, formed by 7 of 11 nearby galaxies, is characterized by an extremely small line-of-sight velocity dispersion of 15 km/s, which differs significantly from the predictions of cosmological simulations of about 70 km/s. This fact allows us to determine the total mass of the Local Group as $M_{LG} = (2.47 \pm 0.15) \times 10^{12}$ $M_\odot$ using an analytical model of the Hubble flow around a spherical overdensity in the standard flat \LCDM{} universe. The practical equality of this mass to the sum of the masses of our Galaxy and the Andromeda Galaxy, as well as the absence of mass growth in the range of distances under consideration, gives grounds to conclude that the entire mass of the Local Group is confined within the virial radii around its two main galaxies. The barycenter, found from the minimal scatter of mass estimates, corresponds to the mass ratio of the Milky Way and the Andromeda Galaxy equal to $M_{MW}/M_{M31} = 0.74\pm0.10$. The velocity of our Galaxy to the barycenter turns out to be $62.6\pm2.6$ km/s. This allows us to determine the apex of the Sun relative to the barycenter of the Local Group to be $(l,b,V) = ( +94.0^\circ \pm 0.7^\circ, -2.7^\circ \pm 0.3^\circ, 301 \pm 3$ km/s in the Galactic coordinates.

A framework is presented within an eleven-dimensional M-theory scenario wherein dynamical geometric moduli, arising from a G2-holonomy compactification, yield an evolving cosmological term Lambda(z). This "geometric vacuum energy" is distinct from conventional dark energy in two main aspects: (i) its origin lies in extra-dimensional fluxes and instanton-like corrections embedded within the moduli potential, and (ii) it exhibits a moderate peak at intermediate redshifts. This feature provides a mechanism to partially mitigate the Hubble tension, elevating the inferred H0 value from ~67 to ~69.5 km s^-1 Mpc^-1 without disrupting the overall concordance of the cosmological model. Furthermore, the inclusion of mild spatial openness (Omega_k0 ~ 0.097) and a slightly reduced matter fraction (Omega_m0 ~ 0.25) relative to standard LambdaCDM allows the geometry-driven Lambda(z) to naturally maintain a cosmic age near 13.8 Gyr. Initial numerical checks indicate that the calibrated model achieves chi^2/nu ~ 1 for H(z) data (chi^2 ~ 26 for 24 points) and yields a structure growth amplitude S8 ~ 0.74, consistent with current observational bounds. While not a definitive solution, these results illustrate a potential pathway for reconciling certain late-time cosmological puzzles through a slight deviation from a pure cosmological constant, rooted in G2-compactified M-theory. The findings suggest a UV-complete geometric origin for dark energy. Further developments, including detailed statistical analyses, N-body simulations, and explicit G2 constructions, are anticipated to refine and test the parameter space of this model.

The role of mergers in the evolution of massive black holes is still unclear, and their dynamical evolution, from the formation of pairs to binaries and the final coalescence, carries large physical uncertainties. The identification of the elusive population of close massive binary black holes (MBBHs) is crucial to understand the importance of mergers in the formation and evolution of supermassive black holes. It has been proposed that MBBHs may display periodic optical/ultra-violet variability. Optical surveys provide photometric measurements of a large variety of objects, over decades and searching for periodicities coming from galaxies in their long-term optical/UV lightcurves may help identify new MBBH candidates. Using the Catalina Real-Time Transient Survey (CRTS) and Zwicky Transient Facility (ZTF) data, we studied the long-term periodicity of variable sources in the centre of galaxies identified using the galaxy catalogue Glade+. We report 36 MBBHs candidates, with sinusoidal variability with amplitudes between 0.1 and 0.8 magnitudes over 3-5 cycles, through fitting 15 years of data. The periodicities are also detected when adding a red noise contribution to the sine model. Moreover, the periodicities are corroborated through Generalized Lomb Scargle (GLS) periodograms analysis, providing supplementary evidence for the observed modulation. We also indicate 58 objects, that were previously proposed to be MBBH candidates from analysis of CRTS data only. Adding ZTF data clearly shows that the previously claimed modulation is due to red noise. We also created a catalogue of 221 weaker candidates which require further observations over the coming years to help validate their nature. Based on our 36 MBBHs candidates, we expect ~20 MBBHs at z<1, which is commensurate with simulations. Further observations will help confirm these results.

A new SPHERE seires complex extensive air showers detector is under development. The main goal of its mission is to study the mass composition of cosmic ray nuclei in the 1-100 PeV energy range at a new level. The already well-established telescope of Cherenkov light reflected from the snow-covered ice surface of Lake Baikal from an altitude of 500-1000 m will be supported by a detector of direct light pointed upward. Since the two detectors will study the same shower at different stages of its development, it could be called a 3D detection, which is completely new for the EAS method. The development is based on an extensive MC modeling of the shower and the detection process using the Supercomputer Complex of the Lomonosov Moscow State University.

High velocity stars move through the interstellar medium with V > 30 km/s. When the star has powerful winds, under the appropriate conditions, the interaction of the wind with the interstellar material produces a system of shocks. The outer shock, called the bowshock, perturbs the ambient medium, heating and compressing the gas. The dust in the compressed bowshock cools, producing infrared radiation. This emission appears as extended coma-shape structures. The discovery of radio nonthermal emission from two stellar bowshock nebulae indicates that these sources might be accelerating electrons up to relativistic energies. The produced nonthermal radio emission is most probably synchrotron which has a high degree of polarization. In this work we model the synchrotron emission of runaway massive star bowshocks aiming to produce synthetic radio emission and polarization maps for two frequencies: 1.40 and 4.86 GHz. We model the interacting plasmas in a steady-state regime by means of magnetohydrodynamics simulations and we compute the injection and transport of the relativistic electrons in the diffusion approximation. We include in the model the most important depolarization effects. Our main conclusions are i) the effects of Faraday rotation within the source are important at the lowest frequency considered, ii) inferring the local magnetic field direction from polarization measurements only can be misleading, iii) thermal radio emission produced by ionized plasma within the bowshock structure and surroundings can surpass the polarized one for the considered frequencies, and iv) the contribution from the background electrons is minor.

Galaxies reside within dark matter halos, but their properties are influenced not only by their halo properties but also by the surrounding environment. We construct an interpretable neural network framework to characterize the surrounding environment of galaxies and investigate the extent to which their properties are affected by neighboring galaxies in IllustrisTNG300 data ($z=0$). Our models predict galaxy properties (stellar mass and star formation rate) given dark matter subhalo properties of both host subhalo and of surrounding galaxies, which serve as an explainable, flexible galaxy-halo connection model. We find that prediction accuracy peaks when incorporating only the nearest neighboring galaxy for stellar mass prediction, while star formation rate prediction benefits from information from up to the third-nearest neighbor. We determine that environmental influence follows a clear hierarchical pattern, with the nearest neighbor providing the dominant contribution that diminishes rapidly with additional neighbors. We confirm that central and satellite galaxies, as well as different galaxy categories based on mass and star-forming activity, exhibit distinct environmental dependencies. Environmental dependence for low-mass galaxies ($\log(M_*/M_\odot) < 10$) shows 35-50% environmental contribution compared to just 8-30% for massive centrals, while satellite galaxies experience consistently stronger environmental effects than centrals across all populations. Furthermore, we find that the most significant attribute from neighboring subhalos for predicting target galaxy properties is its distance to the nearest neighboring galaxy. These quantitative results offer guidance for constructing more sophisticated empirical and semi-analytic models of galaxy formation that explicitly include environmental dependence as a function of galaxy type and mass.

The Star-Forming Main Sequence (SFMS) serves as a critical framework for understanding galaxy evolution, highlighting the relationship between star formation rates (SFR) and stellar masses M_* across cosmic time. Despite its significance, the origin of the 0.3-0.4 dex dispersion in the SFMS remains a key unresolved question. Uncovering the origin of dispersion is crucial for understanding the evolution of galaxies. Using a large sample of approximately 500,000 galaxies, we reveal an unprecedented symmetry in the distribution of key structural properties-effective radius (R_{\rm e}), stellar surface density (M_*/R_{\rm e}^2), and morphology on the SFMS. This symmetry implies that galaxies with high (above SFMS) and low (below SFMS) SFRs share similar fundamental parameters. Moreover, galaxies with smaller R_{\rm e} or higher M_*/R_{\rm e}^2 exhibit greater dispersion in SFR. This dispersion reflects the response to fluctuations in cosmic accretion flows, while the SFR itself represents the time-averaged effect over the gas consumption timescale. Shorter gas consumption timescales, associated with higher M_*/R_{\rm e}^2, lead to greater SFR dispersion. Our results reveal that the variation of SFR originates from the oscillation of accretion flow and is regulated by the stellar surface density.

Using Gaussian process methods, we analyzed the light curves of three extreme solar X-ray flares observed by the RHESSI satellite. Their variability characteristics were then compared with those of HXMT-HE X-ray burst (XRB; in SGR 1935+2154) associated with fast radio burst (FRB) 200428 and blazar $\gamma$-ray giant flares, to investigate the origins of these extreme flaring events. The variability patterns of the solar X-ray flares follow the stochastically driven damped simple harmonic oscillator (SHO) model. The derived timescales $t_{\rm B\underline{} steep}$ and $t_{\rm B\underline{~} flat}$ (corresponding to PSD breaks) are in the range of 4-7 s and 16-53 s, respectively. The FRB-associated HXMT-HE burst has a $Q$ value near 0.3, matching those of the solar flares occurred on 23 July 2002 (flare 1) and 3 November 2003 (flare 2). By contrast, blazar $\gamma$-ray giant flares show $Q >$ 0.3, similar to the solar flare occurred on 25 February 2014 (flare 3). We proposed that the critically damped state of the system may be the condition triggering the association between the XRB in SGR 1935+2154 and the FRB. In this scenario, the critical damping $Q$ value of the system is around 0.3, not the theoretical 0.5. The similarity in $Q$ values might imply that the FRB-associated HXMT-HE XRB and solar X-ray flares 1 $\&$ 2 share comparable dynamic behavior, while blazar $\gamma$-ray flares and solar X-ray flare 3 exhibit another distinct but similar dynamic behavior. Like solar X-ray flares, these extreme flares may all be related to the magnetic reconnection process.

Recent observations by the James Webb Space Telescope have uncovered a population of compact, red object ($z\sim 4\text{--}7$) known as little red dots (LRDs). The presence of broad Balmer emission lines indicates active galactic nuclei powered by supermassive black holes (BHs), while LRDs exhibit unusually weak X-ray and radio emission and low variability, suggesting super-Eddington accretion that obscures the central engine. We suggest that such an extreme accretion disc inevitably drives strong outflows, which would disrupt the LRDs themselves unless confined within the nuclear region -- posing a general feedback problem for overmassive BHs. To resolve this, we propose that the BH is embedded in a massive, optically thick envelope that gravitationally confines the outflow, making any outflow a no-go. This envelope, powered by accretion on to the BH, radiates at nearly the Eddington limit, and is sustained by an infall of the interstellar medium at a rate on the order of $\sim 1 M_{\odot}~{\rm yr}^{-1}$. A photosphere emerges either within the envelope or in the infalling medium, with a characteristic temperature of $5000$ - $7000 {\rm K}$, near the Hayashi limit. The resulting blackbody emission naturally explains the red optical continuum of the distinct V-shaped spectrum observed in most LRDs. Furthermore, the dynamical time-scale at the photosphere, $\sim 0.01~{\rm pc}$, is consistent with the observed year-scale variabilities. The nuclear structure and spectral features of LRDs are shaped by this envelope, which not only regulates feedback but also acts as a gas reservoir that sustains rapid BH growth in the early universe.

The fundamental metallicity relation (FMR) - the three-way trend between galaxy stellar masses, star-formation rates (SFRs) and gaseous metallicities - remains amongst the most studied extragalactic relations. Furthermore, metallicity correlates particularly tightly with gravitational potential. Simulations support a shared origin for these relations relating to long-term gas inflow history variations; however, differences between simulated and observed galaxy samples make it unclear whether this holds for real galaxies. We use MaNGA integral field observations to probe these relations in star-forming galaxies at one effective radius. We confirm the FMR and equivalent relations for stellar metallicity (FMR*) and gaseous N/O (fundamental nitrogen relation, FNR). We find that all relations persist when considering gravitational potential in place of stellar mass and/or considering stellar ages in place of SFR, with the gaseous relations strengthened significantly by considering potential. The gaseous FMR disappears at high masses/potentials, while the FNR persists and the FMR* strengthens. Our results suggest a unified interpretation of galaxies' gaseous and stellar metallicities and their N/O abundances in terms of their formation histories. Deeper gravitational potentials correspond to earlier star-formation histories (SFHs) and faster gas consumption, producing tight potential-abundance relations for stars and gas. In weak potentials, galaxy SFR variations primarily result from recent gas inflows, mostly affecting gas abundances. In deeper potentials, SFR variations instead correspond to broad differences in SFH shapes resulting from differences in long-term gas consumption histories, which is most visible in stellar abundances. This unified interpretation could be confirmed with upcoming higher redshift spectroscopic surveys.

This is a brief review of the recent progress in understanding the evolution of the accretion disks in tidal disruption events (TDEs). Special attention is paid to (1) thermal-viscous instability that causes the disk to transition from a thick state to a thin one, and back and forth, (2) interactions between the fallback material and existing disk. Challenges to the current model from late-time X-ray observations are highlighted and possible solutions are discussed.

Since the start of operations in 2011, the VLT Survey Telescope (VST) has been one of the most efficient wide-field imagers in the optical bands. However, in the next years the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will be a game-changer in this field. Hence, the timing is appropriate for specializing the VST with additions that can make it unique in well-defined scientific cases. VSTPOL is a project that aims to provide the addition of wide-field polarimetric capabilities to the VST telescope, making it the first large survey telescope for linear optical polarimetry. Actually, while there are quite a number of optical telescopes, the telescopes providing polarimetric instrumentation are just a few. The number of relatively large mirror polarimetric telescopes is small, although they would be specifically needed e.g. to support many science cases of the Cherenkov Telescope Array (CTA) that, in the southern hemisphere, is co-located with the VST. The VST telescope is equipped with a single instrument, the OmegaCAM wide-field imaging camera operating in the visible bands with a field of view of 1$^\circ\times1^\circ$. The polarimetric mode will be implemented through the insertion of a large rotatable polarizer installed on the field-corrector optics, which will be exchangeable with the non-polarimetric corrector optics. The limiting polarimetric systematic errors due to variable atmospheric conditions and instrumental polarization can be corrected down to a level of $\sim10^{-3}$ by leveraging the large amount of unpolarized stars within each field-of-view. By the user point of view, VSTPOL will be an additional mode for the VST wide-field imaging camera.

Stellar-mass binary black holes (BBH) may form, evolve, and merge within the dense environments of active galactic nuclei (AGN) disks, thereby contributing to the BBH population detected by gravitational wave (GW) observatories. Mergers occurring in AGN disks may be gravitationally lensed by the supermassive black hole at the AGN centre. The probability of such lensing events has been approximately estimated by using the Einstein criterion in previous work. However, a more reasonable approach to calculating the lensing probability should be based on whether the detector can distinguish the lensed GW waveform from the unlensed one. In this work, we calculate the lensing probability of LIGO sources embedded in AGN disk by relating threshold mismatch to the signal-to-noise ratio of the observed events. For the sensitivity of LIGO-Virgo-KAGRA O3 observation runs, our results indicate that the lensing probability is several times higher than previous estimates. If AGNs are indeed the primary formation channel for BBHs, we could quantify the probability of detecting the lensed GW events in such a scenario. The non-detections, on the other hand, will place stricter constraints on the fraction of AGN disk BBHs and even the birthplaces of BBH mergers.

Abundances of M dwarfs, being the most numerous stellar type in the Galaxy, can enhance our understanding of planet formation processes. They can also be used to study the chemical evolution of the Galaxy, where in particular alpha-capture elements play an important role. We aim to obtain abundances for Fe, Ti, and Ca for a small sample of well-known M dwarfs for which interferometric measurements are available. These stars and their abundances are intended to serve as a benchmark for future large-scale spectroscopic studies. We analysed spectra obtained with the GIANO-B spectrograph. Turbospectrum and the wrapper TSFitPy were used with MARCS atmospheric models in order to fit synthetic spectra to the observed spectra. We performed a differential abundance analysis in which we also analysed a solar spectrum with the same method and then subtracted the derived abundances line-by-line. The median was taken as the final abundance for each element and each star. Our abundances of Fe, Ti, and Ca agree mostly within uncertainties when comparing with other values from the literature. However, there are few studies to compare with.

The mean apparent magnitude of the Qianfan satellites is 5.76 +/- 0.04, while the mean of magnitudes adjusted to a distance of 1,000 km is 5.24 +/- 0.04, based on 1,161 observations. Light curves of several spacecraft display rapid periodic fluctuations which indicate that they are tumbling. Nearly all of the non-tumbling satellite observations can be modeled with diffusely reflecting, Earth-facing surfaces. The Qianfan constellation will impact astronomical research and aesthetic appreciation of the night sky unless their brightness is mitigated.

Magnetic AB stars are known to produce periodic radio pulses by the electron cyclotron maser emission (ECME) mechanism. Only 19 such stars, known as 'Main-sequence Radio Pulse emitters' (MRPs) are currently known. The majority of MRPs have been discovered through targeted observation campaigns that involve carefully selecting a sample of stars that are likely to produce ECME, and which can be detected by a given telescope within reasonable amount of time. These selection criteria inadvertently introduce bias in the resulting sample of MRPs, which affects subsequent investigation of the relation between ECME properties and stellar magnetospheric parameters. The alternative is to use all-sky surveys. Until now, MRP candidates obtained from surveys were identified based on their high circular polarisation ($\gtrsim 30\%$). In this paper, we introduce a complementary strategy, which does not require polarisation information. Using multi-epoch data from the Australian SKA Pathfinder (ASKAP) telescope, we identify four MRP candidates based on the variability in the total intensity light curves. Follow-up observations with the Australia Telescope Compact Array (ATCA) confirm three of them to be MRPs, thereby demonstrating the effectiveness of our strategy. With the expanded sample, we find that ECME is affected by temperature and the magnetic field strength, consistent with past results. There is, however, a degeneracy regarding how the two parameters govern the ECME luminosity for magnetic A and late-B stars (effective temperature $\lesssim 16$ kK). The current sample is also inadequate to investigate the role of stellar rotation, which has been shown to play a key role in driving incoherent radio emission.

The study of the redshifted spectral lines can provide a measure of the fundamental constants over large look-back times. Current grand unified theories predict an evolution in these and astronomical observations offer the only experimental measure of the values of the constants over large time scales. Of particular interest are the dimensionless constants; the fine structure constant ($\alpha$), the proton-electron mass ratio ($\mu$) and the proton g-factor ($g_p$), since these do not require a `standard metre-stick'. Here we present a re-analysis of the 18-cm hydroxyl (OH) lines at $z=0.89$, recently detected with the MeerKAT telescope, towards the radio source PKS\,1830-211. Utilizing the previous constraint of $\Delta\mu/\mu=(-1.8\pm1.2)\times10^{-7}$, we obtain $\Delta (\alpha g_p^{0.27})/(\alpha g_p^{0.27})\lesssim5.7\times10^{-5}$, $\Delta \alpha/\alpha\lesssim2.3\times10^{-3}$ and $\Delta g_p/g_p\lesssim7.9\times10^{-3}$. These new constraints are consistent with no evolution over a look-back time of 7.3 Gyr and provide another valuable data point in the putative evolution in the constants.

The Phocaea asteroid family, one of the large ancient families located in the inner main belt, may be the sources of near-Earth asteroids (NEAs) due to the nearby 3:1 mean motion resonance with Jupiter, the v6 secular resonance, and the Yarkovsky and YORP effects. Thus, understanding the influence of the Yarkovsky and YORP effects on the Phocaea family is one of the keys to figuring out the source of NEAs. However, the physical properties of most of the Phocaea family members are unknown at present. We perform a photometric analysis for 44 asteroids in the Phocaea region using photometric data obtained by ground-based and space-based telescopes (i.e., TESS and Gaia). Based on the derived physical properties, we find significant footprints of the Yarkovsky and YORP effects on the Phocaea family members. Selecting five asteroids nearby the inside boundary of the V-shape in the absolute-magnitude semimajor-axis (H, a) space, we estimate their densities considering their migration in semimajor-axis arises from the Yarkovsky effect. The bulk density of (852) Wladilena ({3.54 g/cm3) suggests a link to the H chondrite meteorites. Incorporating the grain density of the H chondrites, we estimate the macroporosities of the asteroids (290) Bruna, (1164) Kobolda, and (587) Hypsipyle, respectively 41%, 47%, and 65%, implying rubble pile structures. Considering the H chondrites link to asteroid (25) Phocaea, we suggest the parent body of the Phocaea family has been composed of H chondrite like material and the Phocaea family may be one of the sources of H chondrite meteorites.

To describe a general bound binary black hole system, we need to consider orbital eccentricity and the misalignment of black holes' spin vectors with respect to the orbital angular momentum. While binary black holes produced through many formation channels have negligible eccentricity close to merger, they often have a non-negligible eccentricity at formation, and dynamical interactions could produce binaries with non-negligible eccentricity in the bands of current and proposed gravitational-wave (GW) detectors. Another quantity that carries information about the formation channel is the angle between each black hole's spin vector and the binary's orbital angular momentum (referred to as the spin tilt) at formation. The spin tilts inferred in GW astronomy are usually those when the binary is in the band of a GW detector, but these can differ significantly from those at formation. Therefore, it is necessary to evolve the binary back in time to compute the tilts at formation. For many formation scenarios, the tilts in the formal limit of infinite orbital angular momentum, also known as tilts at infinity, are a good approximation to those at formation. We thus generalize the publicly available \texttt{tilts\_at\_infinity} code to compute the tilts at infinity for eccentric, spin-precessing binaries. This code employs hybrid post-Newtonian evolution, starting with orbit-averaged evolution for higher frequencies and then transitioning to precession-averaged evolution to compute the tilts at infinity. We find that the transition frequency used in the quasicircular case still gives acceptably small errors in the eccentric case, and show that eccentricity and hybrid evolution both have a significant effect on the tilts at infinity for many binaries. Finally, we give examples of cases where the tilts at infinity are and are not a good approximation to the tilts at formation in the eccentric case.

We report the discovery of an actinide-boost, very metal-poor ($\left[\mathrm{Fe/H} \right]=-2.38$), $r$-process-enhanced ($\left[\mathrm{Eu/Fe} \right]=0.80$) star, LAMOST J0804+5740, within the Gaia-Sausage-Enceladus (GSE). Based on the high-resolution ($R\sim36,000\; and \;60,000$) and high signal-to-noise ratio spectra obtained with the High Dispersion Spectrograph on the Subaru Telescope, the abundances of 48 species are determined. Its $\log\epsilon\rm(\mathrm{Th}/\mathrm{Eu}) = -0.22 $ establishes it as the first confirmed actinide-boost star within the GSE. Comparative analysis of its abundance pattern with theoretical $r$-process models reveals that the magnetorotationally driven jet supernova $r$-process model with $\hat{L}v$ = 0.2 provides the best fit and successfully reproduces the actinide-boost signature. Kinematic analysis of actinide-boost stars reveals that approximately two-thirds of them are classified as $\textit{ex-situ}$ stars, suggesting that actinide-boost stars are more likely to originate from accreted dwarf galaxies. As the first actinide-boost star identified within the GSE, J0804+5740 will provide valuable insights into $r$-process nucleosynthesis in accreted dwarf galaxies like the GSE, especially on the production of the heaviest elements.

Flows around compact objects are necessarily transonic. Due to their dissipative nature, finding of sonic points is not trivial. Becker and Le in 2003 (BL03) proposed a novel methodology to obtain global transonic solutions, using iterative relaxation technique and exploiting the inner boundary conditions of the central object. In the current work, we propose a generic methodology -- IRM-SP and IRM-SHOCK to obtain any class of global accretion and wind solutions, given a set of constants of motion. We have considered viscosity in the system, which transports angular momentum outwards. In addition, it heats the system. Radiative processes like bremsstrahlung which cools the system is also incorporated. An interplay between heating and cooling process, along with gravity and centrifugal forces gives rise to multiple sonic points and hence shocks. The proposed methodology successfully generates any class of accretion as well as wind solutions, allowing us to unify them. Additionally, we report here rigorously the mathematical as well as the computational algorithm needed, to find sonic point(s) and thus obtain global transonic flows around compact objects.

The radio quasar NVSS~J171822+423948 (J1718+4239) was proposed as the counterpart of the IceCube neutrino event IC-201221A. To reveal the nature of the source, we conducted new very long baseline interferometry (VLBI) observations of this blazar candidate with the Very Long Baseline Array (VLBA). The observations were carried out in dual-band mode between $4$ and $7$~GHz. Archival radio data from the literature were also collected for comparison. Our analysis revealed highly variable, Doppler-boosted radio emission of the source, with compact structure at both kpc and pc scales, a slightly inverted spectrum, and a maximum jet inclination angle of $\theta\le5\degr$. These results confirm J1718+4239 as a blazar-type object.

3D Galactic magnetic fields are critical for understanding the interstellar medium, Galactic foreground polarization, and the propagation of ultra-high-energy cosmic rays. Leveraging recent theoretical insights into anisotropic magnetohydrodynamic (MHD) turbulence, we introduce a deep learning framework to predict the full 3D magnetic field structure-including the plane-of-sky (POS) position angle, line-of-sight (LOS) inclination, magnetic field strength, sonic Mach number ($M_s$), and Alfvén Mach number ($M_A$)-from spectroscopic H~I observations. The deep learning model is trained on synthetic H~I emission data generated from multiphase 3D MHD simulations. We then apply the trained model to observational data from the Commensal Radio Astronomy FAST Survey, presenting maps of 3D magnetic field orientation, magnetic field strength, $M_s$, and $M_A$ for two H~I clouds, a low-velocity cloud (LVC) and an intermediate-velocity cloud (IVC), which overlap in the POS yet reside at different LOS distances. The deep-learning-predicted POS magnetic field position angles align closely with those determined using the velocity gradient technique, whose integrated results are consistent with independent measurements from Planck 353~GHz polarization data. This study demonstrates the potential of deep learning approaches as powerful tools for modeling the 3D distributions of 3D Galactic magnetic fields and turbulence properties throughout the Galaxy.

The interstellar medium (ISM) consists of a multiphase gas, including the warm neutral medium (WNM), the unstable neutral medium (UNM), and the cold neutral medium (CNM). While significant attention has been devoted to understanding the WNM and CNM, the formation of a substantial fraction of the UNM, with temperatures ranging from a few hundred to a few thousand Kelvin, remains less well understood. In this study, we present three-dimensional hydrodynamical and magnetohydrodynamical simulations of turbulent multiphase ISM to investigate the roles of turbulence and magnetic fields in regulating the multiphase ISM. Our results confirm that turbulence is crucial in redistributing energy and producing the UNM. The turbulent mixing effect smooths the phase diagram, flattens the pressure-density relationship, and increases the fraction of gas in the UNM. We find that magnetic fields not only contribute to sustaining the UNM but also influence the dynamics and distribution of gas across all phases. Magnetic fields introduce anisotropy to the turbulent cascade, reducing the efficiency of turbulent mixing in the direction perpendicular to the magnetic field. We find the anisotropy results in a less flat phase diagram compared to hydrodynamical cases. Furthermore, the inclusion of magnetic fields shallowens the second-order velocity structure functions across multiple ISM phases, suggesting that more small-scale fluctuations are driven. These fluctuations contribute to the formation of the UNM by altering the energy cascade and thermodynamic properties of the gas. Our findings suggest that the combined effects of turbulence and magnetic fields are important in regulating the multiphase ISM.

Long uninterrupted observations of the X-ray binary system Her X-1 were performed with the Mikhail Pavlinsky ART-XC telescope of the Spectrum-Röntgen-Gamma (SRG) X-ray Observatory in the 4--25 keV energy range with a total exposure of about two days around the main turn-on of the X-ray source. We present the results of timing and spectral analysis of these observations. The opening of the X-ray source is determined to occur at the orbital phase $\phi_{b}\approx 0.25$. The analysis of the X-ray light curve reveals a first direct observational evidence of the nutation of a tilted precessing accretion disk with a period of $\simeq0.87$ days. The appearance of X-ray pulsations near the orbital phase $\phi_{b}\simeq 0.77$ prior to the main turn-on at the maximum of the nutation variability has been also detected. During the X-ray eclipse, a non-zero X-ray flux is measured, which is presumably associated with scattering of an X-ray emission in a hot corona around the optical star illuminated by the X-rays from the central neutron star. An increase in the X-ray flux after the main turn-on can be described by the passage of the radiation from the central source through a scattering corona above the precessing accretion disk.

To assess the significance and scale dependence of anomalous large scale modes in the CatWISE quasar data, we generate smoothed number density fields on the sphere and study their extreme values -- maximum, minimum, maximum antipodal difference. By comparing these summary statistics to those obtained from random isotropic realisations of the data, we determine the statistical significance of large scale modes as a function of smoothing scale. We perform our analysis using five different versions of the data -- the original quasar map, the maps after separately subtracting the ecliptic bias and the CMB dipole, the map obtained after subtracting both, and the map after subtracting the ecliptic bias and anomalous dipole inferred in \cite{Secrest2021}. We find that the ecliptic-corrected, CMB dipole-removed map exhibits large scale modes that are in tension with random realisations of the data (p-values $p \sim 10^{-4}$), over a wide range of smoothing scales $\pi/8 \leq \delta \leq \pi/2$. The most prominent feature in the data is an underdensity in the southern galactic plane at $(b,\ell) = (-31^\circ,78^\circ)$, which reaches its highest statistical significance when smoothed on scales $\delta = \pi/6$ ($p \ll 10^{-5}$). Notably, the minima statistics align with the maximum antipodal difference statistics, whereas the maxima do not. This suggests that the observed dipole-like behavior in the data is primarily driven by the underdensity in the southern sky. The ecliptic corrected, anomalous dipole subtracted map reduces the significance of any residual anisotropic features, but an underdensity in the south sky persists with p-value $p =0.0018$.

Bosonic stars,hypothetical astrophysical entities, are generally categorized into two primary classes based on the nature of their constituent particles: Einstein Klein Gordon stars, made up of massive scalar bosons, and Proca stars, their vector ''cousins''. Depending on the boson masses and field frequencies, these objects may exhibit properties of diffuse, massive structures, with sizes comparable to or even exceeding those of galaxies. This concept has inspired the bosonic dark matter halo hypothesis, providing a theoretical framework to effectively model the dark matter content of galactic halos. In this paper we build on our previous work to explore the possibility of using vector and scalar bosons to model the components of galactic dark matter halos and subhalos in order to reproduce the observed rotation curves of galaxies. By employing diverse combinations of those bosonic dark matter models in conjunction with observable data for a sample of galaxies, we show that our two component dark matter approach notably improves the agreement between observations and theoretical predictions with respect to our previous investigation. Our framework may shed new light on the enduring mystery surrounding the apparent matter deficit observed in dwarf and spiral galaxies.

Solar Orbiter conducted a series of flare-optimised observing campaigns in 2024 utilising the Major Flare Solar Orbiter Observing Plan (SOOP). Dedicated observations were performed during two distinct perihelia intervals in March/April and October, during which over 22 flares were observed, ranging from B- to M-class. These campaigns leveraged high-resolution and high-cadence observations from the mission's remote-sensing suite, including the High-Resolution EUV Imager (EUI/HRI_EUV), the Spectrometer/Telescope for Imaging X-rays (STIX), the Spectral Imaging of the Coronal Environment (SPICE) spectrometer, and the High Resolution Telescope of the Polarimetric and Helioseismic Imager (PHI/HRT), as well as coordinated ground-based and Earth-orbiting observations. EUI/HRI_EUV operating in short-exposure modes, provided two-second-cadence, non-saturated EUV images, revealing structures and dynamics on scales not previously observed. Simultaneously, STIX captured hard X-ray imaging and spectroscopy of accelerated electrons, while SPICE acquired EUV slit spectroscopy to probe chromospheric and coronal responses. Together, these observations offer an unprecedented view of magnetic reconnection, energy release, particle acceleration, and plasma heating across a broad range of temperatures and spatial scales. These campaigns have generated a rich dataset that will be the subject of numerous future studies addressing Solar Orbiter's top-level science goal: "How do solar eruptions produce energetic particle radiation that fills the heliosphere?". This paper presents the scientific motivations, operational planning, and observational strategies behind the 2024 flare campaigns, along with initial insights into the observed flares. We also discuss lessons learned for optimizing future Solar Orbiter Major Flare campaigns and provide a resource for researchers aiming to utilize these unique observations.

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN~1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM).The 1.7--10 keV Resolve spectra are accurately represented by a single component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58} \times 10^{11}$ s cm$^{-3}$.The Resolve spectra are also well reproduced by the 3-D magneto-hydrodynamic simulation presented by Orlando et al. (2020) suggesting substantial contribution from the this http URL metal abundances obtained with Resolve align with the LMC value, indicating that the X-rays in 2024 originate from non-metal-rich shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of this http URL widths of the atomic lines from Si, S, and Fe correspond to velocities of 1,500--1,700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a 90\% upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2 -- 10 keV range, aligning with NuSTAR findings by Greco et al. (2022). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4} M_{\odot}$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015; Grebenev et al. 2012; Leising 2006).

Context. We present simulations of a massive young star cluster using \textsc{Nbody6++GPU} and \textsc{MOCCA}. The cluster is initially more compact than previously published models, with one million stars, a total mass of $5.86 \times 10^5~\mathrm{M}_{\odot}$, and a half-mass radius of $0.1~\mathrm{pc}$. Aims. We analyse the formation and growth of a very massive star (VMS) through successive stellar collisions and investigate the subsequent formation of an intermediate-mass black hole (IMBH) in the core of a dense star cluster. Methods. We use both direct \textit{N}-body and Monte Carlo simulations, incorporating updated stellar evolution prescriptions (SSE/BSE) tailored to massive stars and VMSs. These include revised treatments of stellar radii, rejuvenation, and mass loss during collisions. While the prescriptions represent reasonable extrapolations into the VMS regime, the internal structure and thermal state of VMSs formed through stellar collisions remain uncertain, and future work may require further refinement. Results. We find that runaway stellar collisions in the cluster core produce a VMS exceeding $5 \times 10^4~\mathrm{M}_{\odot}$ within 5 Myr, which subsequently collapses into an IMBH. Conclusions. Our model suggests that dense stellar environments may enable the formation of very massive stars and massive black hole seeds through runaway stellar collisions. These results provide a potential pathway for early black hole growth in star clusters and offer theoretical context for interpreting recent JWST observations of young, compact clusters at high redshift.

Small (400 to 4000 km) and short lived (10 to 200 km) extreme ultraviolet (EUV) brightenings, detected by the High Resolution Imager EUV (HRIEUV), have been found to be ubiquitous in the Quiet Sun (QS). Their contribution to coronal heating as well as their physical origin are currently being investigated. We wish to determine whether models of short loops and impulsive heating are compatible with the results from observations. In particular, we used two models of loops with distinct thermal properties: cool (T below 1E5 K) and hot loops (T above 1E5 K). We simulated the evolution of impulsively heated short loops, using the 1D hydrodynamics (HD) code HYDRAD. We computed the synthetic light curves of HRIEUV, four EUV channels of the Atmospheric Imaging Assembly (AIA), and five emission lines measured by the SPectral Imaging of the Coronal Environment (SPICE). We then compared the results from the synthetic light curves with observations. The aim was to reproduce the short delays observed between the intensity peaks of the light curves. Cool loops subjected to impulsive heating are good candidates to explain the physical origin of the EUV brightenings. On the other hand, hot loops are not consistent with observations, except when they are subjected to especially strong impulsive heating.

Plages are small concentrations of strong, nearly vertical magnetic fields in the solar photosphere that expand with height. A high spatial and spectral resolution that can resolve their fine structure is required to characterize them, and spectropolarimetric capabilities are needed to infer their magnetic fields. We constrain the 3D fine structure of the magnetic field in the photosphere of a solar plage from a unique spectropolarimetric dataset with a very high spatial and spectral resolution and a fast temporal cadence. We analyzed spectropolarimetric observations of a solar plage in the two magnetically sensitive spectral lines of neutral iron around 630 nm. The observations were obtained with MiHI, which is an integral field unit attached to the Swedish Solar Telescope. MiHI obtained diffraction-limited, high-cadence observations with high spectral fidelity. These observations were interpreted using the spectropolarimetric inversion with magnetohydrostatic constraints, which allowed us to recover the magnetic and thermodynamic structure of the plage on a geometrical scale. The inversion results reveal that the magnetic field can reach up to 2 kG and that it expands significantly from the deep to the mid-photosphere. Weaker (200 G), and very small (subarcsecond) vertical magnetic loops lie beneath this canopy, rooted in the photosphere. This novel picture of a solar plage, in which weak opposite-polarity field patches surround the main polarity, provides new insight into convection in strongly magnetized plasma.

Mid-infrared spectroscopy of protoplanetary disks provides a chemical inventory of gas within a few au, where planets are readily detected around older stars. With the JWST Disk Infrared Spectral Chemistry Survey (JDISCS), we explore demographic trends among 31 disks observed with MIRI (MRS) and with previous ALMA millimeter continuum imaging at high angular resolution (5-10 au). With these S/N $\sim$200-450 spectra, we report emission from H$_2$O, OH, CO, C$_2$H$_2$, HCN, CO$_2$, [Ne II], [Ne III], and [Ar II]. Emission from H$_2$O, OH and CO is nearly ubiquitous for low-mass stars, and detection rates of all molecules are higher than for similar disks observed with Spitzer-IRS. Slab model fits to the molecular emission lines demonstrate that emission from C$_2$H$_2$, HCN, and possibly CO$_2$ is optically thin; thus since column densities and emitting radii are degenerate, observations are actually sensitive to the total molecular mass. C$_2$H$_2$ and HCN emission also typically originate in a hotter region ($920^{+70}_{-130}$, $820^{+70}_{-130}$ K, respectively) than CO$_2$ ($600^{+200}_{-160}$ K). The HCN to cold H$_2$O luminosity ratios are generally smaller in smooth disks, consistent with more efficient water delivery via icy pebbles in the absence of large dust substructures. The molecular emission line luminosities are also correlated with mass accretion rates and infrared spectral indices, similar to trends reported from Spitzer-IRS surveys. This work demonstrates the power of combining multi-wavelength observations to explore inner disk chemistry as a function of outer disk and stellar properties, which will continue to grow as the sample of observed Class II systems expands in the coming JWST observation cycles.

During the HOSTS survey by the LBTI, an excess emission from the main sequence star $\theta$ Boo (F7V spectral type, 14.5pc distance) was observed. This excess indicates the presence of exozodiacal dust near the habitable zone (HZ) of the star. Previous observations from Spitzer and Herschel showed no sign of outer cold dust within their respective detection limits. Additional nulling and high-contrast AO observations were taken to spatially constrain the dust distribution, search for variability, and directly image potential companions in the system. This study presents the results of these observations and provides an interpretation of the inner system's architecture. The star was observed using the LBTI's N'-band nulling mode during three epochs in 2017, 2018, and 2023. The dust distribution is modeled and constrained for each epoch using the standard LBTI nulling pipeline, assuming a vertically thin disk with a face-on inclination. In addition, high-contrast AO observations are performed in the L'-band and H-band to constrain the presence of substellar companions around the star. Several solutions are found for the dust distribution, and for each epoch. However, the LBTI nulling observations are not able to discriminate between them. Using the upper limits from previous observations, we constrain the representative size of the dust grains around 3-5$\mu$m. A tentative increase in dust brightness is also measured at the Earth-equivalent insolation distance between 2017 and 2023. Several options are considered to explain the origin of the observed dust and its variability, but no clear sources could be identified from the current observations. Partly because our high-contrast AO observations could only constrain the presence of companions down to $11M_\text{Jup}$ at 1.3" separation.

The discovery of extra-terrestrial life is one of the ultimate goals for future exoplanet-seeking missions, with one major challenge being the presence of 'exozodiacal' dust near target stars or within their habitable zone. Therefore, it is critical to identify which stars possess exozodiacal dust and quantify their emission levels. In this study, we conducted a search for exozodi candidates within 10 parsecs using the Reyl'e sample. We performed proper motion calculations and cross-matched the sample with the WISE and 2MASS database, resulting in 339 preliminary target samples. We further analysed the infrared radiation characteristics of these targets, using spectral energy distribution (SED) fitting to predict photometric flux levels in the infrared and searching for 3sigma excesses in the WISE W3 band. During further selection processes, we applied various analysis methods to perform rigorous validation. We identified five exozodi candidates all of which are brown dwarfs (BDs). Given the clustering in candidate spectral types, we expect that these are not true exozodi candidates, rather the apparent excess arises from the inability of the BD photosphere models to accurately represent the SEDs of objects at the L-T transition. Indeed, for the object DENIS J025503.3-470049, excess is likely due to silicate clouds in the BD atmosphere. We suggest that a more stringent 5sigma excess is required to infer excess for this spectral type. The detection rate (0/339) in our sample shows that less than 1% M stars have exozodi above 21% excess levels. This is consistent with the rate of exozodi at similar level towards FGK stars in the Kennedy & Wyatt sample (25/24,174). We provide upper limits on the 12 micron exozodi emission for the sample, which is typically at 21% relative to the star. For most stars, in particular the low mass M stars, this is the first such upper limit in the literature.

We present the discovery of a peculiar X-ray transient, EP241021a, by the Einstein Probe (EP) mission, and the results from multiwavelength follow-up observations. The transient was first detected with the Wide-field X-ray Telescope as an intense flare lasting for ~100 s, reaching a luminosity of L_(0.5-4 keV)~10^48 erg/s at z=0.748. Further observations with EP's Follow-up X-ray Telescope reveal a huge drop in the X-ray flux by a factor of >1000 within 1.5 days. After maintaining a nearly plateau phase for ~7 days, the X-ray flux declines as t^-1.2 over a period of ~30 days, followed by a sudden decrease to an undetectable level by EP and XMM-Newton, making it the longest afterglow emission detected among known fast X-ray transients. A bright counterpart at optical and radio wavelengths was also detected, with high peak luminosities in excess of 10^44 erg/s and 10^41 erg/s, respectively. In addition, EP241021a exhibits a non-thermal X-ray spectrum, red optical color, X-ray and optical rebrightenings in the light curves, and fast radio spectral evolution, suggesting that relativistic jets may have been launched. We discuss possible origins of EP241021a, including a choked jet with supernova shock breakout, a merger-triggered magnetar, a highly structured jet, and a repeating partial tidal disruption event involving an intermediate-mass black hole, but none can perfectly explain the multiwavelength properties. EP241021a may represent a new type of X-ray transients with months-duration evolution timescales, and future EP detections and follow-up observations of similar systems will provide statistical samples to understand the underlying mechanisms at work.

Classical T Tauri stars are newly formed, low mass stars which may display both periodic and random variations in their brightness. The interaction between the star and its circumstellar disk is time-dependent, leading to short or long-term changes in the environment, and hence variability of the system. By compiling a large dataset with high-cadence photometric (Kepler, TESS), and high-resolution spectroscopic observations (CFHT/ESPaDOnS) of the highly variable T Tauri star DR Tau, we aim to examine the short- and long-term variability of the system, and identify the underlying physical mechanisms. Our results reveal that DR Tau exhibits stochastic photometric variability not only on daily, but also on hourly timescale with peak-to-peak amplitude of 1.4 mag probably originating from accretion related variations. Our ground-based multifilter photometry shows that the amplitude of the variability decreases with increasing wavelength. This trend towards the infrared wavelengths suggests that part of the disk may be optically thick and invariable. The spectroscopic analysis showed that the H$\alpha$ line presents the most complex line profile with several components but the significance of the components changes over time. This suggests the presence and variation of both accretion flow and wind. Broad and narrow components can be clearly distinguished in the He I and the Ca II lines, suggesting contribution from both the accretion flow and the post-shock region. DR Tau exhibits high level of photometric and spectroscopic variability on both short- and long-timescales, which is caused by the combination of accretion, wind, stellar activity, and obscuration by circumstellar matter; and the significance of the physical mechanisms causing the observed variability changes over time.

We review recent advancements in cosmology with galaxy clusters. Galaxy clusters are the most massive objects in the Universe. Consequently the cluster number density as a function of cluster mass, or cluster abundance, is sensitive to cosmological parameters, particularly the matter density of the Universe $\Omega_{\rm m}$ and the amplitude of matter density fluctuation $\sigma_8$. In this review, we describe the methods used to detect galaxy clusters through optical near-infrared (O-NIR), X-ray, and CMB observations, outlining the advantages and disadvantages of cluster detection through different wavelengths. We describe methods for measuring cluster mass, with a particular focus on calibration by WL measurements. We then discuss how the connection between observables in different wavelengths and cluster abundance can be modeled through a cluster selection function and MOR, and quantify the impact of marginalization of nuisance parameters on cosmological constraints. Finally, we also walk through the recent results of cosmological constraints by cluster abundance with the O-NIR, X-ray, and CMB observations.

Centimeter-sized dust grains-pebbles-are necessary for planetesimal formation via the streaming instability, they play an important role in forming protoplanetary cores and giant planets, as well as enriching their atmospheres with chemical elements. This work investigates the effect of luminosity outbursts on the abundance of pebbles and their ice mantles in protoplanetary disks. We perform global simulations of formation and evolution of a self-gravitating viscous protoplanetary disk using the 2D hydrodynamic thin-disk FEOSAD code, which self-consistently reproduces luminosity outbursts. The model includes thermal balance, dust evolution and its interaction with gas, development of magnetorotational instability, adsorption and desorption of four volatile compounds (H$_2$O, CO$_2$, CH$_4$ and CO), and the feedback of ice mantles on dust fragmentation properties. We show that luminosity outbursts have a stronger effect on the positions of CO$_2$, CH$_4$ and CO snowlines compared to the water snowline. This is because the H$_2$O snowline falls within the viscous heating dominated region during early disk evolution stages, while snowlines of other molecules are located in regions dominated by stellar irradiation heating and are thus more sensitive to temperature changes during outbursts. Nevertheless, luminosity outbursts reduce the total amount of pebbles in the disk by half due to destruction of dust aggregates into monomers following the loss of water ice that binds the aggregates together. Pebble recovery occurs over several thousand years after the outburst ends due to collisional coagulation, with recovery timescales significantly exceeding water freeze-out times. Ice mantle desorption occurs in a complex non-axisymmetric 2D region of the disk, associated with spiral substructure formation during early evolution of gravitationally unstable disks.

UV radiation from OB stars can drive ``external'' photoevaporative winds from discs in clusters, that have been shown to be important for disc evolution and planet formation. However, cluster dynamics can complicate the interpretation of this process. A significant fraction of OB stars are runaways, propagating at high velocity which might dominate over the wider cluster dynamics in setting the time variation of the UV field in part of the cluster. We explore the impact of a runaway OB star on discs and the observational impact that may have. We find that discs exposed to even short periods of strong irradiation are significantly truncated, and only rebound slightly following the ``flyby'' of the UV source. This is predicted to leave an observable imprint on a disc population, with those downstream of the OB star vector being more massive and extended than those upstream. Because external photoevaporation acts quickly, this imprint is less susceptible to being washed out by cluster dynamics for faster runaway OB stars. The Gaia proper motion vector of the B star 42 Ori in NGC 1977 is transverse to the low mass stellar population and so may make a good region to search for this signature in resolved disc observations.

Cosmological theories suggest that the angular momentum of galaxies should be closely linked to the structure of the cosmic web. The Local Supercluster is the closest and most studied structure where the orientation of galaxy spins can be studied. As noted by Navarro et al. (2004), the use of edge-on galaxies greatly simplifies this task by reducing it to an analysis of the distribution of position angles in the Supergalactic coordinates. We reexamine this correlation using modern catalogs that allow us to perform a more robust statistical analysis. We test the dependence on redshift, spatial position, luminosity, and color. We find that the spins of galaxies with stellar mass $M_*<10^{8.7}$ $M_\odot$ show a weak tendency to be aligned perpendicular to the plane of the Local Supercluster at the 2-sigma level. Other subsamples do not show statistically significant correlations.

Dwarf AGN serve as the ideal systems for identifying intermediate mass black holes (IMBHs) down to the most elusive regimes ($\sim 10^3 - 10^4 M_{\odot}$). However, the ubiquitously metal-poor nature of dwarf galaxies gives rise to ultraluminous X-ray sources (ULXs) that can mimic the spectral signatures of IMBH excitation. We present a novel photoionization model suite that simultaneously incorporates IMBHs and ULXs in a metal-poor, highly star-forming environment. We account for changes in $M_{BH}$ according to formation seeding channels and metallicity, and changes in ULX populations with post-starburst age and metallicity. We find that broadband X-rays and UV emission lines are insensitive to $M_{BH}$ and largely unable to distinguish between ULXs and IMBHs. Many optical diagnostic diagrams cannot correctly identify dwarf AGN. The notable exceptions include He~II~$\lambda$4686 and [O~I]~$\lambda$6300, for which we redefine typical demarcations to account for ULX contributions. Emission lines in the mid-IR show the most promise in separating stellar, ULX, IMBH, and shock excitation while presenting sensitivity to $M_{BH}$ and $f_{\text{AGN}}$. Overall, our results expose the potential biases in identifying and characterizing dwarf AGN purely on strong line ratios and diagnostic diagrams rather than holistically evaluating the entire spectrum. As a proof of concept, we argue that recently discovered over-massive BHs in high-$z$ JWST AGN might not represent the overall BH population, with many galaxies in these samples potentially being falsely classified as purely star-forming.

Clouds are often considered a highly uncertain barrier for detecting biosignatures on exoplanets, especially given intuition gained from transit surveys. However, for direct imaging reflected light observations, clouds could increase the observational signal by increasing reflected light. Here we constrain the impact of clouds on the detection of O2 and O3 by a direct imaging telescope such as the Habitable Worlds Observatory (HWO) using observations simulated with the Planetary Spectrum Generator (PSG). We first perform sensitivity tests to show that low clouds enhance O2 and O3 detectability while high clouds diminish it, and the effect is greater when cloud particles are smaller. We next apply clouds produced by the cloud microphysics model CARMA with varied planetary parameters and clouds drawn from observations of different types of clouds on Earth to PSG. We find that clouds are likely to increase the SNR of O2 and O3 for terrestrial exoplanets under a wide range of scenarios. This work provides important constraints on the impact of clouds on observations by telescopes including HWO.

Stars form from molecular gas under complex conditions influenced by multiple competing physical mechanisms, such as gravity, turbulence, and magnetic fields. However, accurately identifying the fraction of gas actively involved in star formation remains challenging. Using dust continuum observations from the Herschel Space Observatory, we derived column density maps and their associated probability distribution functions (N-PDFs). Assuming the power-law component in the N-PDFs corresponds to gravitationally bound (and thus star-forming) gas, we analyzed a diverse sample of molecular clouds spanning a wide range of mass and turbulence conditions. This sample included 21 molecular clouds from the solar neighborhood ($d<$500 pc) and 16 high-mass star-forming molecular clouds. For these two groups, we employed the counts of young stellar objects (YSOs) and mid-/far-infrared luminosities as proxies for star formation rates (SFR), respectively. Both groups revealed a tight linear correlation between the mass of gravitationally bound gas and the SFR, suggesting a universally constant star formation efficiency in the gravitationally bound gas phase. The star-forming gas mass derived from threshold column densities ($N_{\mbox {threshold}}$) varies from cloud to cloud and is widely distributed over the range of $\sim$1--17$\times$10$^{21}$ cm$^{-2}$ based on N-PDF analysis. But in solar neighborhood clouds, it is in rough consistency with the traditional approach using $A_{\rm V}$ $\ge$ 8 mag. In contrast, in high turbulent regions (e.g., the Central Molecular Zone) where the classical approach fails, the gravitationally bound gas mass and SFR still follow the same correlation as other high-mass star-forming regions in the Milky Way. Our findings also strongly support the interpretation that gas in the power-law component of the N-PDF is undergoing self-gravitational collapse to form stars.

One of the most profound empirical laws of star formation is the Gao-Solomon relation, a linear correlation between the star formation rate (SFR) and the dense molecular gas mass. It is puzzling how the complicated physics in star-formation results in this surprisingly simple proportionality. Using archival Herschel and Atacama Large Millimeter/submillimeter Array Observations, we derived the masses of the most massive cores ($M^{\rm max}_{\rm core}$) and masses of the gravitationally bound gas ($ M_{\rm gas}^{\rm bound}$) in the parent molecular clouds for a sample of low-mass and high-mass star-forming regions. We discovered a significant correlation $\log(M^{\rm max}_{\rm core}/M_{\odot}) = 0.506 \log(M_{\rm gas}^{\rm bound}/M_{\odot})-0.32$. Our discovered $M^{\rm max}_{\rm core}$-$M_{\rm gas}^{\rm bound}$ correlation can be approximately converted to the Gao-Solomon relation if there is (1) a constant 30% efficiency of converting $M^{\rm max}_{\rm core}$ to the mass of the most massive star ($m^{\rm max}_{\rm star}$), and (2) if SFR and $m^{\rm max}_{\rm star}$ are tightly related through $\log({\rm SFR}/(M_{\odot} {\rm yr}^{-1})) = 2.04 \log(m^{\rm max}_{\rm star}/M_{\odot})-5.80$. Intriguingly, both requirements have been suggested by previous theoretical studies (c.f. Yan et al. 2017). Based on this result, we hypothesize that the Gao-Solomon relation is a consequence of combining the following three non-trivial relations (i) SFR vs. $m^{\rm max}_{\rm star}$, (ii) $m^{\rm max}_{\rm star}$ vs. $M^{\rm max}_{\rm core}$, and (iii) $M^{\rm max}_{\rm core}$ vs. $M_{\rm gas}^{\rm bound}$. This finding may open a new possibility to understand the Gao-Solomon relation in an analytic sense.

Analysis of Callisto's moments of inertia, derived from Galileo's gravity data, suggests that its structure is not fully differentiated. This possibly undifferentiated state contrasts sharply with the globally molten state inferred in its counterpart, Ganymede, and poses unique challenges to theories of the formation and evolution of the Galilean moons. During their formation, both moons experienced multiple heating mechanisms, including tidal heating, radiogenic heating from short-lived radionuclides, accretional heating from impacts, and heat from the surrounding circumplanetary disk. Our study investigates the optimal conditions required to account for Callisto's partially differentiated state in contrast to Ganymede's complete differentiation. We investigate crucial accretion parameters, such as the timing of accretion onset, the duration of accretion, and the impactor size distribution. We find that the observed dichotomy between Ganymede and Callisto can be attributed to similar formation conditions, assuming an identical impactor size distribution and composition in the Jovian circumplanetary disk. The key differences in the formation of Ganymede and Callisto are the disk temperature at their respective formation locations and their final radii. Our results indicate that both moons accreted gradually over more than 2 Myr, concluding at least 5.5 Myr after the formation of calcium-aluminum-rich inclusions in the protosolar nebula. Our model demonstrates that Callisto can remain undifferentiated despite accreting a substantial influx of kilometer-sized impactors, potentially contributing up to 30% of the total mass inflow, while still allowing for the complete differentiation of Ganymede.

The generation of large curvature perturbations associated with the production of primordial black holes is studied in the context of a Higgs inflaton. To enable this amplification, we consider an inflationary model in which the tree-level action for gravity and the Standard Model Higgs is modified by quantum corrections, described by a series of higher-dimension operators. Finally within a minimal EFT framework, we present two viable models in which the spectrum of curvature perturbations generated by the Higgs field is consistent with CMB observations and can lead to the formation of primordial black holes in the asteroid mass range, potentially accounting for the entirety of dark matter.

In recent years, numerous arguments have emerged suggesting that the LCDM (Lambda Cold Dark Matter) model may be inconsistent with observational data, requiring more or less radical revisions. Notable examples include the Hubble tension, the discrepancy between early and late-universe measurements of the Hubble constant, as well as tensions in measurements of cosmic structure growth. These issues have led some to question the validity of the LCDM framework and consider possible modifications or alternative models. However, upon closer inspection, many of these critiques stem from methodological or interpretive disagreements rather than from clear falsifications in the strict Popperian sense. Karl Popper proposed that scientific theories must be testable and falsifiable; in other words, a theory should be rejected if it fails a specific, reproducible test. Yet, many of the alleged inconsistencies within LCDM, while statistically significant, are not definitive falsifications but rather indicators of areas needing refinement or more complex modeling within the same framework. Thus, I review the recent claims about LCDM's limitations and analyze why they often reflect individual biases or philosophical preferences, rather than rigorous scientific falsifications. For example, alternative cosmological models such as MOND (Modified Newtonian Dynamics) or models incorporating new physics like quintessence or modified gravity are sometimes advocated based on theoretical appeal rather than direct evidence from critical tests. In many cases, these arguments for falsifying LCDM reveal more about subjective interpretations of data than about concrete observational contradictions.

Polarity reversals are a key feature of Earth's magnetic field, yet the processes governing them are still poorly understood. Dipole reversals have been found in many numerical dynamo simulations and often occur close to the transition between dipolar and multipolar regimes. Simulated conditions are far from those in Earth's liquid iron core because of the long runtimes needed to capture polarity transitions. We develop a unidimensional path theory in an attempt to simplify the search for the dipole-multipole transition at increasingly realistic physical conditions. We build 3 paths, all based on a constant magnetic Reynolds number $Rm$; one aiming for Magnetic, Coriolis, and Archimedean (MAC), and 2 aiming for inertia-MAC force balance. We add inertia due to its role in simulated reversals. Results show reasonable agreement with predictions within the accessible parameter space, but deviate from predicted behaviour for certain quantities, e.g. magnetic field strength and magnetic/kinetic energy ratio. Further, simulations move into the dipolar non-reversing regime as they are advanced along the path. By increasing the buoyancy driving (via higher Rayleigh number) above the values predicted by the path theory, we are able to access the dipole-multipole transition down to an Ekman number $E\sim 10^{-6}$, comparable to the most extreme conditions reported to date. Results demonstrate that our approach is an efficient method for seeking the dipole-multipole transition at low $E$. However, the conditions under which we access the dipole-multipole transition become increasingly hard to access numerically and also increasingly unrealistic because $Rm$ rises beyond plausible bounds inferred from geophysical observations. Future work combining path theory with variations in the core buoyancy distribution, appears a promising approach to accessing the transition at extreme physical conditions.

Virtual meetings have long been the outcast of scientific interaction. For many of us, the COVID-19 pandemic has only strengthened that sentiment as countless Zoom meetings have left us bored and exhausted. But remote conferences do not have to be negative experiences. If well designed, they have some distinct advantages over conventional in-person meetings, including universal access, longevity of content, as well as minimal costs and carbon footprint. This article details our experiences as organizers of a successful fully virtual scientific conference, the KITP program "Fundamentals of Gaseous Halos" hosted over 8 weeks in winter 2021. Herein, we provide detailed recommendations on planning and optimization of remote meetings, with application to traditional in-person events as well. We hope these suggestions will assist organizers of future virtual conferences and workshops.

We extend the holographic V-QCD model by introducing a charged scalar field sector to represent the condensation of paired quark matter in the deconfined phase. By incorporating this new sector into the previously established framework for nuclear and quark matter, we obtain a phase diagram that, in addition to the first-order deconfinement transition and its critical end-point, also features a second-order transition between paired and unpaired quark matter. The critical temperature for quark pairing exhibits only a mild dependence on the chemical potential and can reach values as high as $T_\mathrm{crit} \approx 30~\rm MeV$. Comparison of the growth rate for the formation of homogeneous paired phases to the growth rate of previously discovered modulated phases suggests that the former is subdominant to the latter.

When combining apparently inconsistent experimental results, one often implements errors on errors. The Particle Data Group's phenomenological prescription offers a practical solution but lacks a firm theoretical foundation. To address this, D'Agostini and Cowan have proposed Bayesian and frequentist approaches, respectively, both introducing gamma-distributed auxiliary variables to model uncertainty in quoted errors. In this Letter, we show that these two formulations admit a parameter-by-parameter correspondence, and are structurally equivalent. This identification clarifies how Bayesian prior choices can be interpreted in terms of frequentist sampling assumptions, providing a unified probabilistic framework for modeling uncertainty in quoted variances.

The IceCube Upgrade will provide unprecedented sensitivity to dark matter particles annihilating in the core of the Sun. For dark matter candidates with spin-dependent couplings to nuclei and that annihilate significantly to tau leptons or neutrinos, we find that the IceCube Upgrade will be capable of testing parameter space that is beyond the reach of existing direct detection experiments. After calculating the sensitivity of the IceCube Upgrade to dark matter annihilation in the Sun, we explore dark matter models that could be tested by this experiment, identifying two classes of scenarios as promising targets for such searches.

Axions have attracted attention promising candidates for dark matter (DM). Although axions have been intensively searched for, they have not been observed yet. Recently, novel experiments to search for axion DM have been proposed that use optical cavities to amplify polarization rotation of laser light induced by the axion-photon interaction. One such experiment employs a ring cavity composed of four mirrors. However, its sensitivity to the axion-photon coupling $\gag$ in the low axion mass region is limited due to the reflection phase difference between s- and p-polarizations. In this paper, we propose a new method to improve the sensitivity using zero-phase shift mirrors and a wavelength tunable laser. Moreover, the laser makes it easier to scan the high axion mass region by tuning the reflection phase difference between s- and p-polarizations. We experimentally confirmed that the phase difference generated upon reflection on a zero phase shift mirror satisfies the requirement of $8.6 \times 10^{-3}~\si{deg}$, which corresponds to the half width at half maximum (HWHM) of the p-polarization with the mirror fixed on a folded cavity and a wavelength tunable laser.

In this paper, we report on the implementation of the EOM spin-flip (SF), ionization-potential (IP) and electron-affinity (EA) coupled cluster singles doubles (CCSD) methods for atoms and molecules in strong magnetic fields for energies as well as one-electron properties. Moreover, non-perturbative triples corrections using the EOM-CCSD(T)(a)* scheme were implemented in the finite-field framework for the EE, SF, IP, and EA variants. These developments allow the access to a large variety of electronic states as well as the investigation of IPs and EAs in a strong magnetic field. The latter two indicate the relative stability of the different oxidation states of elements. The increased flexibility to target challenging electronic states and the access to the electronic states of the anion and cation are important for the assignment of spectra from strongly magnetic White Dwarf (WD) stars. Here, we investigate the development of the IPs and EAs in the presence of a magnetic field for the elements of the first and second row of the periodic table. In addition, we study the development of the electronic structure of Na, Mg, and Ca that aided in the assignment of metal lines in a magnetic WD. Lastly, we investigate the electronic excitations of CH in different magnetic-field orientations and strengths, a molecule that has been found in the atmospheres of WD stars.

Gas hydrates are considered fundamental building blocks of giant icy planets like Neptune and similar exoplanets. The existence of these materials in the interiors of giant icy planets, which are subject to high pressures and temperatures, depends on their stability relative to their constituent components. In this study, we reexamine the structural stability and hydrogen content of hydrogen hydrates, (H2O)(H2)n, up to 104 GPa, focusing on hydrogen-rich materials. Using synchrotron single-crystal X-ray diffraction, Raman spectroscopy, and first-principles theoretical calculations, we find that the C2-filled ice phase undergoes a transformation to C3-filled ice phase over a broad pressure range of 47 - 104 GPa at room temperature. The C3 phase contains twice as much molecular H2 as the C2 phase. Heating the C2-filled ice above approximately 1500 K induces the transition to the C3 phase at pressures as low as 47 GPa. Upon decompression, this phase remains metastable down to 40 GPa. These findings establish new stability limits for hydrates, with implications for hydrogen storage and the interiors of planetary bodies.

In this paper we study gravitationally bound compact objects sourced by a string theory inspired Born-Infeld scalar field. Unlike many of their canonical scalar field counterparts, these ``boson stars'' do not have to extend out to infinity and may generate compact bodies. We analyze in detail both the junction conditions at the surface as well as the boundary conditions at the center which are required in order to have a smooth structure throughout the object and into the exterior vacuum region. These junction conditions, although involved, turn out to be relatively easy to satisfy. Analysis reveals that these compact objects have a richer structure than the canonical boson stars and some of these properties turn out to be physically peculiar: There are several branches of solutions depending on how the junction conditions are realized. Further analysis illustrates that in practice the junction conditions tend to require interior geometries reminiscent of ``bag of gold'' spacetimes, and also hide the star behind an event horizon in its exterior. The surface compactness of such objects, defined here as the ratio $2M/r$, can be made arbitrarily close to unity indicating the absence of a Buchdahl bound. Some comments on the stability of these objects is provided to find possible stable and unstable regimes. However, we argue that even in the possibly stable regime the event horizon in the vacuum region shielding the object is potentially unstable, and would cut off the star from the rest of the universe.

We revisit non-minimally coupled scalar field cosmologies in the Einstein frame and present a comprehensive analysis that spans background dynamics, linear perturbations, thermodynamics, quantum gravity constraints and baryogenesis. Using a dynamical systems approach, we classify all analytical critical points for a representative set of scalar field potentials: axions (ALPs), cyclic ekpyrotic, exponential (ekpyrotic) with $\Lambda$, quintessence, and scalar field dark matter. We show that a chameleon-like coupling $f(\phi)=e^{\beta\phi}$ modifies both the expansion history and the growth of structure in a way that remains compatible with current fifth-force searches. Analytical transfer matrices for primordial tensor modes are derived, revealing a scale-dependent break in the gravitational-wave spectrum whose position and amplitude are fixed by the coupling parameter $\beta$. A causal Israel--Stewart treatment demonstrates that the same coupling supplies an effective bulk pressure that drives a smooth bounce while respecting quantum energy inequalities. Penrose's Weyl-curvature hypothesis is recovered dynamically: during an ekpyrotic phase with $\omega\gg1$ the Weyl invariant decays as $a^{-6(1+\omega)}$, resetting gravitational entropy without violating the generalized second law. A time-varying $\phi$ simultaneously generates an effective chemical potential $\mu_{B}\simeq\beta\dot{\phi}/M^{2}$, allowing for spontaneous high-scale baryogenesis whose back reaction on $\phi$ is negligible for $M\gtrsim10^{15}\sqrt{\beta}$ GeV. Swampland distance and de Sitter criteria are automatically satisfied because the ekpyrotic phase confines the field excursion to $\Delta\phi\lesssim\mathcal{O}(1)M_{\rm Pl}$. These results establish new links between dark-sector interactions, entropy production, and late-time acceleration while identifying observational windows.

One of the most remarkable discoveries by the IBEX is the ribbon - a narrow band of enhanced ENA fluxes observed in the sky. The prevailing explanation attributes the IBEX ribbon to the secondary ENA mechanism. In this process, primary H ENAs, produced via charge exchange between solar wind (SW) protons and interstellar H atoms within the heliosphere, travel beyond the heliopause (HP) and undergo further charge exchange with protons of the LISM, generating pickup protons. Some of these pickup protons subsequently experience charge exchange with interstellar H atoms, forming secondary ENAs, some of which travel back toward the Sun and are detected by the IBEX. This paper presents a kinetic model developed to simulate secondary ENA fluxes. Ribbon simulations are performed using global distributions of plasma and H atoms in the heliosphere derived from a kinetic-MHD model of the SW/LISM interaction. The model accounts for all relevant primary ENA populations, including neutralized thermal SW protons, neutralized pickup protons, and ENAs originating in the inner heliosheath (IHS). The transport of pickup protons beyond the HP is described by the focused transport equation for a gyrotropic velocity distribution in the scatter-free limit, assuming no pitch-angle scattering or energy diffusion. Our simulations qualitatively reproduce IBEX-Hi ribbon observations and exhibit good quantitative agreement at low heliolatitudes. However, the model underestimates fluxes at high heliolatitudes, likely due to the omission of non-stationary SW behavior in the stationary framework used in this work. The study highlights the importance of ENAs from the IHS, a population considered for the ribbon production in the frame of the kinetic model of pickup proton transport in the heliosphere for the first time, for accurately reproducing ribbon fluxes observed by IBEX-Hi at the highest energy steps.

Recently, the Atacama Cosmology Telescope (ACT) collaboration has reported a scalar spectral index $n_s~=~0.9743~\pm~0.0034$. This is substantially larger than the classical prediction of Higgs Inflation, $n_s \approx 0.965$, which is 2.74$\sigma$ below the ACT mean value. We show that when an otherwise metastable Standard Model Higgs Inflation potential is stabilised by the addition of vector-like quark pairs and the potential is renormalised in the Jordan frame, the value of $n_s$ is generally larger than 0.965 and can explain the ACT observation. As an example, assuming the 2022 PDG direct measurement central value for the t quark mass, $m_{t} = 172.69$ GeV, and central values for the SM inputs to the renormalisation group equations, we obtain $n_s = 0.9792 - 0.9844$ for the case of three isosinglet vector-like B quarks with mass $m_{Q}$ in the range 1-3 TeV, with the lowest value of the $n_s$ range being 1.44$\sigma$ above the ACT mean value. The model predicts primordial gravitational wave with tensor-to-scalar ratio $r = 7.87 \times 10^{-3} - 1.21 \times 10^{-2}$ for $m_{Q} =$ 1-3 TeV, which will be easily observable in forthcoming CMB experiments. Observation of vector-like quarks of mass close to 1 TeV mass combined with a large tensor-to-scalar ratio $r \sim 0.1$ would support the model.

The SUPerconduction AXion search experiment (Supax) is a haloscope designed to probe axion-like particles (ALPs) as candidates for dark matter and solutions to the strong CP problem. ALPs are predicted to couple to photons, allowing their detection through resonant conversion in electromagnetic cavities placed within strong magnetic fields. \Supax employs a 12 T magnetic field and tunable superconducting cavities with resonance frequencies ranging from 2\,GHz to 7.2\,GHz, enabling the exploration of axion masses between 8\,$\mu$eV and 30\,$\mu$eV. The tuning mechanism, based on piezo motors and gas-pressure regulation, allows for simultaneous scanning of up to three frequencies, significantly improving search efficiency. This paper presents the technical design of the Supax experiment, preliminary R\&D efforts, and results from prototype experiments. In particular, we exclude dark photons with masses around $35\,\mu$eV with a kinetic mixing parameter $\chi > 5\cdot 10^{-14}$, i.e. a region of parameter space which has not been previously explored.

The astrophysical direct nuclear capture reaction $^{12}{\rm C}(p, \gamma)^{13}{\rm N}$ is studied within the framework of a potential model. Parameters of the nuclear $p-^{12}$C interaction potentials of the Woods-Saxon form are adjusted to reproduce experimental $p-^{12}$C scattering phase shifts, as well as the binding energies and empirical values of the asymptotic normalization coefficient (ANC) for the $^{13}$N(1/2$^-$) ground state from the literature. The reaction rates are found to be very sensitive to the description of the value of the ANC of the $^{13}$N($1/2^{-}$) ground state and width of the $^{13}$N($1/2^+$) resonance at the $E_x=2.365$ MeV excitation energy. The potential model, which yields the ANC value of 1.63 fm$^{-1/2}$ for the $^{13}$N($1/2^{-}$) ground state and a value $\Gamma$=39 keV for the $^{13}$N($1/2^+$) resonance width, is able to reproduce the astrophysical $S$ factor in the energy interval up to 2 MeV, the empirical values of the reaction rates in the temperature region up to $T=10^{10}$ K of the LUNA Collaboration and the results of the R-matrix fit. The astrophysical factor $S(0)=1.35$ keV b is found using the asymptotic expansion method of D. Baye. The obtained value is in a good agreement with the Solar Fusion II result. At the same time, the calculated value of 1.44 keV b of the astrophysical $S$ factor at the Solar Gamow energy is consistent with the result of the R-matrix fit of $S(25~\rm{keV})=1.48 \pm 0.09$ keV b by Kettner {\it et al.}, but slightly less than the result of $S(25~\rm{keV})=1.53 \pm 0.06$ keV b the LUNA Collaboration.

Phenomenological calculations of the properties of dense matter, such as relativistic mean-field theories, represent a pathway to predicting the microscopic and macroscopic properties of neutron stars. However, such theories do not generically have well-controlled uncertainties and may break down within neutron stars. To faithfully represent the uncertainty in this breakdown scale, we develop a hybrid representation of the dense-matter equation of state, which assumes the form of a relativistic mean-field theory at low densities, while remaining agnostic to any nuclear theory at high densities. To achieve this, we use a nonparametric equation of state model to incorporate the correlations of the underlying relativistic mean-field theory equation of state at low pressures and transition to more flexible correlations above some chosen pressure scale. We perform astrophysical inference under various choices of the transition pressure between the theory-informed and theory-agnostic models. We further study whether the chosen relativistic mean-field theory breaks down above some particular pressure and find no such evidence. Using simulated data for future astrophysical observations at about two-to-three times the precision of current constraints, we show that our method can identify the breakdown pressure associated with a potential strong phase transition.

In particle physics and cosmology, distinguishing subtle new physics signals from established backgrounds is always a fundamental challenge for individual phenomenologists. This paper presents a simple and robust statistical framework to evaluate the compatibility of highly motivated (HM) theoretical models with the residual of the experimental results, focusing on scenarios where data appears consistent with background predictions. We develop a likelihood ratio test procedure that compares null and alternative hypotheses, emphasizing cases where new physics introduces small deviations from the background. We demonstrate the approach through two concrete examples, a localized excess and a modulation over an oscillatory background. We derive explicit conditions under which the effect of the background on the residual of the data must be accounted for. The framework's practicality is highlighted and in addition to the limitation, strategies to simplify complex background modeling are discussed.