Papers for Wednesday, May 07 2025

The question of whether a dynamo can be triggered by gravitational collapse is of great interest, especially for the early Universe. Here, we employ supercomoving coordinates to study the magnetic field amplification from decaying turbulence during gravitational collapse. We perform three-dimensional simulations and show that for large magnetic Reynolds numbers, there can be exponential growth of the comoving magnetic field with conformal time before the decay of turbulence impedes further amplification. The collapse dynamics only affects the nonlinear feedback from the Lorentz force, which diminishes more rapidly for shorter collapse times, allowing nearly kinematic continued growth. We also confirm that helical turbulence is more efficient in driving dynamo action than nonhelical turbulence, but this difference decreases for larger collapse times. We also show that for nearly irrotational flows, dynamo amplification is still possible, but it is always associated with a growth of vorticity -- even if it still remains very small. In nonmagnetic runs, the growth of vorticity is associated with viscosity and grows with the Mach number. In the presence of magnetic fields, vorticity emerges from the curl of the Lorentz force.

We investigate the hierarchical distribution and evolution of young stellar structures in the dwarf starburst galaxy NGC 4449 using data from the GULP survey. By analyzing the spatial distribution of field stars younger than 100 Myr, we identify large-scale stellar complexes and substructures using HDBSCAN -- a density-based clustering algorithm -- and trace their evolution over time. While comparing these stellar structures in different regions of the galaxy, we find that the central bar-like region shows a clear expansion of the structures within the first $\sim$ 60 Myrs, while the arm-like structure in the NE shows no discernible trend, possibly due to external perturbations from tidal interactions with a neighboring galaxy. An age-dependent two point correlation function (TPCF) analysis shows that young stars exhibit a strong hierarchical distribution, with clustering strength decreasing over time. The power-law slope of the TPCF, which starts at $\alpha \sim 0.65$ for stars younger than 5 Myr, shows a slight decline to $\alpha \sim 0.4$ for stars older than 50 Myr, though it does not reach a completely flat (random) distribution. This trend indicates a subtle weakening of structural hierarchy among young ($<$100 Myr) stars, which is primarily driven by internal stellar motions. Future work will extend this analysis to the remaining 26 galaxies in the GULP survey to better constrain the role of the galactic environment in shaping the hierarchical evolution of young stellar populations.

Context. Galaxy clusters assemble over gigayears in a very anisotropic environment, which causes a remarkable diversity in their mass assembly histories (MAH). Aims. In this work, we have aimed to understand how the present-day properties of the dark matter halo and the intracluster medium are related to the whole evolution of these structures. Methods. To this end, we analysed a $\Lambda$CDM hydrodynamical+$N$-Body simulation of a $(100 \, h^{-1} \mathrm{Mpc})^3$ volume, containing over 30 clusters and 300 groups, and looked at the individual and the stacked MAHs (determined from complete merger trees) in relation to properties of the DM haloes and the ICM at fixed cosmic time (indicators of assembly state). Results. The ensemble MAHs are well separated when stacked in bins of these indicators, yielding clear dependencies of evolutionary properties (such as formation redshift) on fix-time halo properties. Additionally, we find that different indicators are informative about distinct epochs of accretion. Finally, by summarising the complex MAH diversity with two parameters, we describe how different indicators bring complementary information in different directions of this biparametric space. Overall, halo spin and a combined indicator appear to be the ones encoding the most information about the MAH. Conclusions. The results shown here add up to the idea that the dynamical state of cosmic structures is a multifaceted concept, and warn that single indicators are incapable of capturing the whole complexity of the process. This work sheds light on the nature of this characterisation by untangling precisely when and how several indicators are informative about. In turn, this can provide clues to better constrain the MAH of observed structures.

We present a spectroscopic search for broad H$\alpha$ emitters at z~3.7-6.5 in the GOODS-N field, utilizing JWST/NIRCam slitless spectroscopy from FRESCO and CONGRESS, complemented by JADES imaging. We identify 19 broad H$\alpha$ emitters with FWHM > 1000 km/s at z~4-5.5, including 9 new sources. The black hole masses and AGN bolometric luminosities, inferred from the broad H$\alpha$ components, indicate that most sources are accreting at ~10% of the Eddington limit. We derive their host stellar masses via SED fitting and find higher $M_{BH}/M_{*}$ ratios relative to the local $M_{BH}$-$M_{*}$ relations, consistent with previous studies. We find that 42% of the sample do not satisfy the widely-used color selection criteria for Little Red Dots (LRDs), with the majority of these sources lacking the characteristic steep red slope. A comparison of the average SEDs between our sample and LRDs selected in the same field reveals that the steep red slopes observed in some LRDs are likely due to line-boosting effects as previously suggested. Furthermore, we find that 68% of color-selected LRDs with H$\alpha$ detections in the NIRCam/Grism spectra do not exhibit broad-line features. While the limited sensitivity of the grism spectra may hinder the detection of broad-line components in faint sources, our findings still highlight the enigmatic nature of the LRD population.

We study the large-scale environments and clustering properties of 28 low-luminosity AGNs at $z=3.9-6$ in the GOODS-N field. Our sample, identified from the JWST NIRCam Imaging and WFSS data in CONGRESS and FRESCO surveys with either broad H$\alpha$ emission lines or V-shape continua, are compared to 782 H$\alpha$ emitters (HAEs) selected from the same data. These AGNs are located in diverse large-scale environments and do not preferentially reside in denser environments compared to HAEs. Their overdensity field, $\delta$, averaged over (15 $h^{-1}$cMpc)$^3$, ranges from $-0.56$ to 10.56, and shows no clear correlation with broad-line luminosity, black hole (BH) masses, or the AGN fraction. It suggests that $> 10$ cMpc structures do not significantly influence BH growth. We measure the two-point cross-correlation function of AGNs with HAEs, finding a comparable amplitude to that of the HAE auto-correlation. This indicates similar bias parameters and host dark matter halo masses for AGNs and HAEs. The correlation length of field AGNs is 4.26 $h^{-1}$cMpc, and 7.66 $h^{-1}$cMpc at $3.9 < z < 5$ and $5 < z < 6$, respectively. We infer a median host dark matter halo mass of $\log (M_h/M_\odot)\approx 11.0-11.2$ and host stellar masses of $\log (M_\star/M_\odot) \approx 8.4-8.6$ by comparing with the UniverseMachine simulation. Our clustering analysis suggests that low-luminosity AGNs at high redshift reside in normal star-forming galaxies with overmassive BHs. They represent an intrinsically distinct population from luminous quasars and could be a common phase in galaxy evolution.

This paper examines the evolution of cosmic filaments across redshifts 1, 0.5, and 0 using the IllustrisTNG100-1 magneto-hydrodynamical simulation. To achieve this, we introduce GrAviPaSt, a simple and efficient filament identification method that leverages gravitational potential, an A*-like path-finding algorithm, and spanning trees. Applying this method to galaxy distributions at different redshifts allows us to analyze various filament properties, including their length, thickness, mass density contrast, and radial profile. Additionally, we investigate dynamic characteristics such as the mean distance of filament galaxies from the skeleton, their weighted mean velocity, and velocity trends normalized by their positions within the filaments. Our findings reveal the evolution of cosmic filaments from redshift 1 to 0, highlighting key differences across classifications. In particular, we examine the mass density contrast radial profile of filaments connecting two galaxy groups and those linking two galaxy clusters, identifying distinct differences in profile shape between these categories. Furthermore, in the context of weighted mean velocity, we analyze cosmic filaments exhibiting either negative or positive weighted mean velocity, demonstrating their differing evolutionary trends in terms of the mean distance of filament galaxies from the skeleton.

Context: Filaments are common features in molecular clouds and they play a key role in star formation (SF). Studying their life cycle is essential to fully understand the SF process. Aims: We aim to characterise the impact of magnetic field ($B$) and stellar feedback on the evolution of filamentary structures in star-forming clouds. Methods: We performed two numerical simulations of a collapsing $10^4\,$M$_\odot$ cloud with different mass-to-flux ratios ($\mu=2$ and $\mu=8$), including early stellar feedback (jets and HII regions). Using DisPerSE, we extracted the three-dimensional filamentary network and analysed its properties as it evolves throughout the SF event. Results: We observed that the filamentary network in the simulations follow two distinct evolutionary pathways. In the strongly magnetised case, the cloud maintains a sparser filamentary network, and the arising filaments are predominantly perpendicular to $B$ lines. With a weak field, the cloud develops a single central hub, with converging filaments favouring a parallel alignment relative to $B$. Furthermore, while always accreting, filaments exhibit faster flows towards the hub relative to the surrounding gas. In the weakly magnetised run, the central hub dominates the dynamics, and filaments exhibit faster flows as they approach the central hub. Finally, once the expanding HII region impacts the filaments, they align to $B$ independently of the initial configuration. Conclusions: Magnetic fields play a critical role in shaping the structure and dynamics of molecular clouds. Stronger magnetic fields slow the cloud's evolution and inhibit the formation of central hubs, promoting a broader filamentary network instead. However, ionising feedback dominates the late-stage evolution, overriding the initial differences and dictating the final filament configuration.

Correlations in the Lyman-$\alpha$ (Ly$\alpha$) forest, both as a function of line of sight separation (1D) and 3D separation, provide a unique window to the distribution of matter at redshifts not accessible by current galaxy surveys. While optimal quadratic estimators have been used to measure 1D correlations, they are computationally expensive and difficult to extend to 3D analyses. On the other hand, estimators based on the Fast Fourier Transform (FFT) are significantly faster, but are affected by missing data in the spectra (masked pixels) and so far have not used pixel weights to reduce the uncertainties in the measurement. In this publication we describe how to compute the window matrix that enables forward-modeling the impact of masked pixels and weights on the FFT-based estimators. We use Gaussian and hydrodynamical simulations with artificially masked pixels to validate the method on the measurement of 1D correlations. Finally, we show that the formalism can be extended to model the impact on 3D correlations, in particular on the cross-spectrum, the correlation of 1D Fourier modes as a function of transverse separation. This work will enable more precise clustering measurements with the Ly$\alpha$ forest dataset recently collected by the Dark Energy Spectroscopic Instrument (DESI).

Over the past 16 years, the Fermi Large Area Telescope (LAT) has significantly advanced our view of the GeV gamma-ray sky, yet several key questions remain - such as the composition of the isotropic gamma-ray background, the origin of the Fermi Bubbles or the potential presence of signatures from exotic physics like dark matter. Addressing these challenges requires sophisticated astrophysical modeling and robust statistical methods capable of handling high-dimensional parameter spaces. In this work, we analyze 14 years of high-latitude ($|b|\geq30^{\circ}$) Fermi-LAT data in the range from 1 to 10 GeV using simulation-based inference (SBI) via neural ratio estimation. This approach allows us to detect individual gamma-ray sources and derive a list of significant gamma-ray emitters containing more than 98\% of all sources listed in the Fermi-LAT Fourth Source Catalog (4FGL) with a flux $S>3\times10^{-10}\;\mathrm{cm}^{-2}\,\mathrm{s}^{-1}$ (about a factor of three larger than the flux above which 4FGL is nearly complete), without any non-4FGL source detected in that flux range. Additionally, we reconstruct the source-count distribution in both parametric and non-parametric forms, achieving large agreement with previous literature results as well as those sources detected by our SBI pipeline. We also quantitatively validate our gamma-ray emission simulator via an anomaly detection technique, demonstrating that the synthetic data closely reproduces the complexity of the real observations.

We apply the framework of Vlasov Perturbation Theory (VPT) to the two-loop matter power spectrum within $\Lambda$CDM cosmologies. The main difference to Standard Perturbation Theory (SPT) arises from taking the velocity dispersion tensor into account, and the resulting screening of the backreaction of UV modes renders loop integrals cutoff-independent. VPT is informed about non-perturbative small scale dynamics via the average value of the dispersion generated by shell-crossing, which impacts the evolution of perturbations on weakly non-linear scales. When using an average dispersion from halo models, the VPT power spectrum agrees with the one from the simulation, up to differences from missing three-loop contributions. Alternatively, treating the average dispersion as free parameter we find a remarkably stable prediction of the matter power spectrum from collisionless dynamics at percent level for a wide range of the dispersion scale. We quantify the impact of truncating the Vlasov hierarchy for the cumulants of the phase-space distribution function, finding that the two-loop matter power spectrum is robust to neglecting third and higher cumulants. Finally, we introduce and validate a simplified fast scheme fVPT that can be easily incorporated into existing codes and is as numerically efficient as SPT.

We present early multi-wavelength photometric and spectroscopic observations of the Type IIb supernova SN 2024uwq, capturing its shock-cooling emission phase and double-peaked light curve evolution. Early spectra reveal broad H-alpha (v ~ 15,500 km s$^{-1}$) and He I P-Cygni profiles of similar strengths. Over time the He I lines increase in strength while the H-alpha decreases, consistent with a hydrogen envelope ($M_{env}$ = 0.7 - 1.35 $M_\odot$ ) overlying helium-rich ejecta. Analytic modeling of early shock cooling emission and bolometric light analysis constrains the progenitor to a partially stripped star with radius R = 10 - 60 $R_\odot$, consistent with a blue/yellow supergiant with an initial ZAMS mass of 12 - 20 $M_\odot$ , likely stripped via binary interaction. SN 2024uwq occupies a transitional position between compact and extended Type IIb supernovae, highlighting the role of binary mass-transfer efficiency in shaping a continuum of stripped-envelope progenitors. Our results underscore the importance of both early UV/optical observations to characterize shock breakout signatures critical to map the diversity in evolutionary pathways of massive stars. Upcoming time domain surveys including Rubin Observatory's LSST and UV missions like ULTRASAT and UVEX will revolutionise our ability to systematically capture these early signatures, probing the full diversity of stripped progenitors and their explosive endpoints.

We investigate potential deviations from the standard adiabatic evolution of the cosmic microwave background (CMB) temperature, $T_{\rm CMB}(z)$, using the latest Sunyaev-Zeldovich (SZ) effect measurements and molecular line excitation data, covering a combined redshift range of $0 < z \lesssim 6$. We follow different approaches. First, we reconstruct the redshift evolution of $T_{\rm CMB}(z)$ in a model-independent way using Gaussian Process regression. The tightest constraints come from SZ measurements at $z < 1$, while molecular line data at $z > 3$ yield broader uncertainties. By combining both datasets, we find good consistency with the standard evolution across the full analysed redshift range, inferring a present-day CMB monopole temperature of $T_0 = 2.744 \pm 0.019$ K. Next, we test for deviations from the standard scaling by adopting the parameterisation $T_{\rm CMB}(z) = T_0(1+z)^{1-\beta}$, where $\beta$ quantifies departures from adiabaticity, with $\beta = 0$ corresponding to the standard scenario. In this framework, we use Gaussian Process reconstruction to test the consistency of $\beta = 0$ across the full redshift range and perform $\chi^2$ minimisation techniques to determine the best-fit values of $T_0$ and $\beta$. In both cases, we find good consistency with the standard temperature-redshift relation. The $\chi^2$-minimisation analysis yields best-fit values of $\beta = -0.0106 \pm 0.0124$ and $T_0 = 2.7276 \pm 0.0095$ K, in excellent agreement with both $\beta = 0$ and independent direct measurements of $T_0$ from FIRAS and ARCADE. We discuss the implications of our findings, which offer strong empirical support for the standard cosmological prediction and place tight constraints on a wide range of alternative scenarios of interest in the context of cosmological tensions and fundamental physics.

We develop an orbit-superposition method for edge-on barred galaxies, and evaluate its capability to recover the bar pattern speed $\rm\Omega_p$. We select three simulated galaxies Au-18, Au-23, and Au-28 with known pattern speeds from the Auriga simulations, and create MUSE-like mock data sets with edge-on views (inclination angles $\theta_{\rm T}\ge85^\circ$) and various bar azimuthal angles $\varphi_{\rm T}$. For mock data sets with side-on bars ($\varphi_{\rm T}\ge50^\circ$), the model-recovered pattern speeds $\rm\Omega_p$ encompass the true pattern speeds $\rm\Omega_T$ within the model uncertainties ($1\sigma$ confidence levels, $68\%$) for 10 of 12 cases. The average model uncertainty within the $1\sigma$ confidence levels is equal to $10\%$. For mock data sets with end-on bars ($\varphi_{\rm T}\le30^\circ$), the model uncertainties of $\rm\Omega_p$ depend significantly on the bar azimuthal angles $\varphi_{\rm T}$, with the uncertainties of cases with $\varphi_{\rm T}=10^\circ$ approaching $\sim30\%$. However, by imposing a stricter constraint on the bar morphology ($p_{\rm bar}\le0.50$), the average uncertainties are reduced to $14\%$, and $\rm\Omega_p$ still encompass $\rm\Omega_T$ within the model uncertainties for three of four cases. For all the models that we create in this paper, the $2\sigma$ ($95\%$) confidence levels of the model-recovered pattern speeds $\rm\Omega_p$ always cover the true values $\rm\Omega_T$.

The origin of carbon in the Universe remains uncertain. At solar metallicity, binary-stripped massive stars -- stars that lost their envelope through stable interaction with a companion -- have been suggested to produce twice as much carbon as their single-star counterparts. However, understanding the chemical evolution of galaxies over cosmic time requires examining stellar yields across a range of metallicities. Using the stellar evolution code MESA, we compute the carbon yields from wind mass loss and supernova explosions of single and binary-stripped stars across a wide range of initial masses ($10$-$46,M_\odot$), metallicities ($Z = 0.0021$, $0.0047$, $0.0142$), and initial orbital periods ($10$-$5000$ days). We find that metallicity is the dominant factor influencing the carbon yields of massive stars, outweighing the effects of binarity and detailed orbital parameters. Since the chemical yields from binary massive stars are highly sensitive to metallicity, we caution that yields predicted at solar metallicity should not be directly extrapolated to lower metallicities. At sub-solar metallicities, particularly below $1/7$ solar, weak stellar winds and inefficient binary stripping result in carbon yields from binary-stripped stars that closely resemble those of single stars. This suggests that binary-stripped massive stars are unlikely to explain the presence of carbon-enhanced metal-poor stars or the carbon enrichment observed in high-redshift galaxies as probed by the James Webb Space Telescope. Our findings only concern the stripped stars in massive binaries. The impact of other outcomes of binary star evolution, in particular stellar mergers and accretors, remains largely unexplored but will be necessary for a full understanding of the role of massive binaries in nucleosynthesis.

The correlations between Supermassive Black Holes (SMBHs) and their host galaxies still defy our understanding from both the observational and theoretical perspectives. Here we perform pairwise residual analysis on the latest sample of local inactive galaxies with a uniform calibration of their photometric properties and with dynamically measured masses of their central SMBHs. The residuals reveal that stellar velocity dispersion $\sigma$ and, possibly host dark matter halo mass $M_{\rm halo}$, appear as the galactic properties most correlated with SMBH mass, with a secondary (weaker) correlation with spheroidal (bulge) mass $M_{\rm sph}$, as also corroborated by additional Machine Learning tests. These findings may favour energetic/kinetic feedback from Active Galactic Nuclei (AGN) as the main driver in shaping SMBH scaling relations. Two state-of-the-art hydrodynamic simulations, inclusive of kinetic AGN feedback, are able to broadly capture the mean trends observed in the residuals, although they tend to either favour $M_{\rm sph}$ as the most fundamental property, or generate too flat residuals. Increasing AGN feedback kinetic output does not improve the comparison with the data. In the Appendix we also show that the galaxies with dynamically measured SMBHs are biased high in $\sigma$ at fixed luminosity with respect to the full sample of local galaxies, proving that this bias is not a byproduct of stellar mass discrepancies. Overall, our results suggest that probing the SMBH-galaxy scaling relations in terms of total stellar mass alone may induce biases, and that either current data sets are incomplete, and/or that more insightful modelling is required to fully reproduce observations.

Context: Globular clusters (GCs) host multiple populations characterised by abundance variations in a number of light elements. In many cases, these populations also show spatial and/or kinematic differences, which vary in strength from cluster to cluster and tend to decrease with the clusters' dynamical ages. Aims: In this work, we aim to study the dynamical mixing of multiple populations and establish a link between the more theoretical aspects of the mixing process and various observational parameters that quantify differences between the populations' spatial concentration and velocity anisotropy. Methods: We follow the dynamical mixing of multiple populations in a set of numerical simulations through their distribution in the energy and angular momentum phase space and quantify the evolution of their degree of dynamical mixing. Results: We present the degree of dynamical mixing traced by the intrinsic differences in the phase space distribution of the populations. We compare the differences in phase-space with three observable quantities that describe the degree of mixing in the structural and kinematic differences of the populations: $A^{+}$, commonly used in the literature for spatial differences; and we introduce two new parameters, $\Delta A_{\beta}$ that traces the difference in velocity anisotropy and $\sigma_{\text{Lz}}$, that traces the angular momentum distribution of stars. Conclusions: Our study provides new insights into the dynamics of phase space mixing of multiple populations in globular clusters. We show that differences between the 1P and the 2P observed in old clusters contain key information on the cluster's dynamics and the 1P and the 2P spatial and kinematic properties set by the formation processes, but caution is necessary in using the strength of the present-day differences to quantitatively constrain those imprinted at the time of formation.

The formation channels of the gravitational-wave (GW) sources detected by LIGO/Virgo/KAGRA (LVK) remain poorly constrained. Active galactic nucleus (AGN) has been proposed as one of the potential hosts but the fraction of GW events originating from AGNs has not been quantified. Here, we constrain the AGN-origin fraction $f_{\rm agn}$ by analyzing the spatial correlation between GW source localizations ($O1\!-\!O4$a) and AGN distributions (SDSS DR16). We find evidence of an excess of low-luminosity ($L_{\rm bol} \le 10^{45}~\!\mathrm{erg~s}^{-1}$) as well as low-Eddington ratio ($\lambda_{\rm Edd} \le 0.05$) AGNs around the LVK events, the explanation of which requires $f_{\rm agn} = 0.39^{+0.41}_{-0.32}$ and $0.29^{+0.40}_{-0.25}$ (90\% confidence level) of the LVK events originating from these respective AGN populations. Monte Carlo simulations confirm that this correlation is unlikely to arise from random coincidence, further supported by anomalous variation of the error of $f_{\rm agn}$ with GW event counts. These results provide the first observational evidence for GW sources coming from either low-luminosity or low-accretion-rate AGNs, offering critical insights into the environmental dependencies of the formation of GW sources.

Virtually all extragalactic use cases of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) require the use of galaxy redshift information, yet the vast majority of its sample of tens of billions of galaxies will lack high-fidelity spectroscopic measurements thereof, instead relying on photometric redshifts (photo-$z$) subject to systematic imprecision and inaccuracy best encapsulated by photo-$z$ probability density functions (PDFs). We present the version 1 release of Redshift Assessment Infrastructure Layers (RAIL), an open source Python library for at-scale probabilistic photo-$z$ estimation, initiated by the LSST Dark Energy Science Collaboration (DESC) with contributions from the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) Frameworks team. RAIL's three subpackages provide modular tools for end-to-end stress-testing, including a forward modeling suite to generate realistically complex photometry, a unified API for estimating per-galaxy and ensemble redshift PDFs by an extensible set of algorithms, and built-in metrics of both photo-$z$ PDFs and point estimates. RAIL serves as a flexible toolkit enabling the derivation and optimization of photo-$z$ data products at scale for a variety of science goals and is not specific to LSST data. We thus describe to the extragalactic science community, including and beyond Rubin the design and functionality of the RAIL software library so that any researcher may have access to its wide array of photo-$z$ characterization and assessment tools.

We explore an extension of the $\Lambda$CDM model in which the pressure $p$ of the dark energy (DE) fluid evolves with the expansion of the Universe, expressed as a function of the scale factor $a$. The corresponding energy density $\rho$ is derived from the continuity equation, resulting in a dynamical equation-of-state parameter $w \equiv p/\rho$ during the late-time expansion of the Universe. The pressure is modeled using a Taylor expansion around the present epoch ($a = 1$), introducing deviations from a cosmological constant within the dynamical dark energy (DDE) framework. At first order, a single new parameter $\Omega_1$ captures linear deviations, while a second-order parameter, $\Omega_2$, accounts for quadratic evolution in the pressure. We constrain the first- and second-order DDE models using multiple observational datasets and compare their performance against $\Lambda$CDM and the CPL parameterization. A joint analysis of Planck CMB, DESI, and DESY5 data yields the strongest evidence for DDE, with a $2.7\sigma$ deviation in the first-order model and over $4\sigma$ in the second-order model, providing strong statistical support for a departure from a cosmological constant. The reconstructed DE evolution in the second-order case reveals a distinctive non-monotonic behavior in both energy density and $w_{\rm DE}(a)$, including clear phantom-crossing phenomena. Notably, the late-time evolution of $w_{\rm DE}(a)$ remains consistent across datasets and shows strong agreement with the CPL parameterization, underscoring the robustness of the pressure-based approach.

Ultraluminous X-ray sources (ULXs) have captivated researchers for decades due to their exceptionally high luminosities and unique spectral characteristics. Some of these sources defy expectations by exhibiting super-Eddington luminosities with respect to stellar mass sources even in their low-hard state. Numerical steady-state calculations suggest that ULXs in this state can be explained as highly magnetized advective accretion sources around stellar-mass black holes. To explore this further, we employ GRMHD simulations using the publicly available code, BHAC (Black Hole Accretion Code), to model the behavior of highly magnetized advective accretion flows around a black hole. Our simulations demonstrate that such systems can indeed produce the intense luminosities observed in ULXs. Additionally, we validate that the magnetic fields required for these high emissions are of the order of $10^7$ Gauss, consistent with previous numerical steady-state findings.

We present extensive ultraviolet to optical photometric and optical to near-infrared (NIR) spectroscopic follow-up observations of the nearby intermediate-luminosity ($M_V = -$16.81$\pm$0.19~mag) Type Iax supernova (SN) 2024pxl in NGC 6384. SN~2024pxl exhibits a faster light curve evolution than the high-luminosity members of this class, and slower than low-luminosity events. The observationally well-constrained rise time of $\sim$10 days and an estimated synthesized $^{56}$Ni mass of 0.03 M$_\odot$, based on analytical modeling of the pseudobolometric light curve, are consistent with models of the weak deflagration of a carbon-oxygen white dwarf. Our optical spectral sequence of SN~2024pxl shows weak \ion{Si}{2} lines and spectral evolution similar to other high-luminosity Type Iax SNe, but also prominent early-time \ion{C}{2} line, like lower-luminosity Type Iax SNe. The late-time optical spectrum of SN~2024pxl closely matches that of SN 2014dt, and its NIR spectral evolution aligns with those of other well-studied, high-luminosity Type Iax SNe. The spectral-line expansion velocities of SN~2024pxl are at the lower end of the Type Iax SN velocity distribution, and the velocity distribution of iron-group elements compared to intermediate-mass elements suggests that the ejecta are mixed on large scales, as expected in pure deflagration models. SN~2024pxl exhibits characteristics intermediate between those of high-luminosity and low-luminosity Type~Iax SNe, further establishing a link across this diverse class.

In spite of accounting for only a small fraction of the mass of the Interstellar Medium (ISM), dust plays a primary role in many physical and chemical processes in the Universe. It is the main driver of extinction of radiation in the UV/optical wavelength range and a primary source of thermal IR emission. Dust grains contain most of the refractory elements of the ISM and they host chemical processes that involve complex molecular compounds. However, observational evidence suggests that grain structure is highly non-trivial and that dust particles are characterized by granularity, asymmetry and stratification, which significantly affect their interaction with radiation fields. Accurate modeling of such interaction is fundamental to properly explain observational results, but it is a computationally demanding task. Here we present the possibility to investigate the effects of radiation/particle interactions in non-spherically symmetric conditions using a novel implementation of the Transition Matrix formalism, designed to run on scalable parallel hardware facilities.

The process of mass-exchange between the components of High Mass X-ray Binary (HMXB) systems with neutron stars undergoing wind-fed accretion is discussed. The X-ray luminosity of these systems allows us to evaluate the mass capture rate by the neutron star from the stellar wind of its massive companion and set limits on the relative velocity between the neutron star and the wind. We find that the upper limit to the wind velocity in the orbital plane during the high state of the X-ray source is in the range 120--1000 km/s, which is by a factor of 2--4 lower than both the terminal wind velocity and the speed of the wind flowing out from the polar regions of massive stars for all the objects under investigation. This finding is valid not only for the systems with Be stars, but also for those system in which the optical components do not exhibit the Be phenomenon. We also show that the lower limit to the radial wind velocity in these systems can unlikely be smaller than a few per cent of the orbital velocity of the neutron star. This provide us with a new constraint on the mass-transfer process in the outflowing disks of Be-type stars.

Globular clusters (GCs) provide statistically significant coeval populations of stars spanning various evolutionary stages, allowing robust constraints on stellar evolution model parameters and ages. We analyze eight old Milky Way GCs with metallicities between [Fe/H] $=-2.31$ and $-0.77$ by comparing theoretical isochrone sets from the Dartmouth Stellar Evolution Program to HST observations. The theoretical isochrones include uncertainties introduced by $21$ stellar evolution parameters such as convective mixing, opacity, diffusion, and nuclear reactions, capturing much of the quantifiable physics used in our code. For each isochrone, we construct synthetic color-magnitude diagrams (CMD) near the main-sequence turn-off region and apply two full-CMD-fitting methods to fit HST ACS data across a range of distance and reddening and measure the absolute age of each GC from the resulting posterior distribution, which accounts for uncertainties in the stellar models, observations, and fitting method. The resulting best-fitting absolute ages range from $\approx 11.5$ to $13.5$ Gyr, with a typical error of $0.5-0.75$ Gyr; the data show a clear trend toward older ages at lower metallicities. Notably, distance and reddening account for over $50\%$ of the uncertainty in age determination in each case, with metallicity, $\alpha$ abundance, mixing length, and helium diffusion being the most important stellar physics parameters for the error budget. We also provide an absolute age-metallicity relation for Milky Way GCs.

The presence of infrared excesses around stars directly correlates to spatially-resolved imaging detections of circumstellar disks at both mm and optical/near-infrared wavelengths. High contrast imagers have resolved dozens of circumstellar disks with scattered light polarimetric imaging. Many of these detections are members of the Scorpius-Centaurus OB association, demonstrating it to be a rich sample for investigating planetary system architectures and planet-disk interactions. With the goal of expanding the sample of directly imaged debris disks in Scorpius-Centaurus, we conducted the Disks In Scorpius-Centaurus Survey (DISCS) leveraging knowledge of high-IR excesses and the power of high contrast polarimetric differential imaging. In combination with the GPIES polarimetric disk survey, we observe seven new Scorpius-Centaurus targets to achieve a 60% complete survey of debris disks with IR excesses exceeding $2.5\times10^{-4}$, resolving four new debris disks. HD 98363, HD 109832, and HD 146181 are resolved for the first time, and HD 112810 is resolved for the first time in polarized intensity. We identify morphological structures that may be indications of planet-disk interactions. We place the systems in the greater context of resolved debris disks, identifying factors of ten differences in scattered light contrast for a given IR excess and implying gaps in our understanding of the smallest and largest dust grains of a system. We conclude that while thermal emission measurements are correlated with scattered light detection, they poorly predict the magnitude of scattered light brightness. We also establish Scorpius-Centaurus debris disks as critical benchmarks in understanding the properties of disks in the scattering regime.

We present the results from the spectral analysis of Cygnus X-3 using simultaneous data from AstroSat and Insight-HXMT during its soft state. A pure reflection spectrum, including emission lines of iron, silica, and sulfur, provides a good fit to the spectra. Orbital phase-resolved analysis shows no significant spectral parameter variations, except for the normalization. Leveraging IXPE polarization results, we model the funnel-shaped geometry and estimate scattered flux and observed polarization for various funnel parameters and observer inclinations. We consider two scenarios: reflection from the funnel walls and scattering by gas within the funnel. Our results reconfirm previous findings, showing that reflection can produce a high polarization degree (PD) of 23$\%$, but not a low PD of 12$\%$. Conversely, scattering can produce a PD of 10-12$\%$, but not as high as 23$\%$ for a fixed observer inclination of $30^\circ$. Scattering results align with previous findings without absorption, but with absorption, PD drops significantly with increasing funnel opening angle. Thus, we can identify common funnel parameters that can produce the different observed PDs in the soft and hard states. The intrinsic luminosity of the source was estimated by comparing the results from a plane disk and the funnel model, to be $\sim7$ $\times$ 10$^{40}$ erg/s for 12$\%$ PD (scattering) and $\sim5$ $\times$ 10$^{41}$ erg/s for 23$\%$ PD (reflection). However, for the reflection model, the luminosity may decrease to $\sim$ 10$^{40}$ erg/s when the 23$\%$ PD observed is taken as a lower limit.

The ASIACO project is an atlas of cometary imaging and spectroscopy obtained from the Asiago Astrophysical Observatory between 2012 and 2024, still ongoing. Up to now, we monitored 41 comets, acquiring broad-band images using the INAF-OAPd 1.82m Copernico and the 0.67/0.92m Schmidt telescopes, as well as low-to-medium resolution spectra with the DFA-UniPd 1.22m Galileo Telescope. Our dataset offers significant scientific potential, enabling studies on the true extent of cometary tails, orientation of nuclei spin axes, retrieval of compositional data, and determination of production rates for key molecular and atomic species. The Atlas will become publicly available for scientific research, and the project remains active, continuously adding monitored comets.

We study the spectroscopic appearance of very massive stars and their effect on the integrated light of starbursts at low metallicity (Z). We adopt two frameworks for the mass loss rates of VMS: in one case we assume no Z dependence, in the other case we assume a linear scaling with Z. We compute evolutionary models for masses 150, 200, 250 and 300 Msun at Z=0.2, 0.1 and 0.01 Zsun. We compute the associated synthetic spectra at selected points along the evolutionary tracks. Finally we build population synthesis models including VMS. We find that the evolution of VMS critically depends on the assumptions regarding mass loss rates. In case of no Z dependence VMS remain hot for all their lifetime. Conversely when mass loss rates are reduced because of lower Z VMS follow a classical evolution towards the red part of the HR diagram. VMS display HeII 1640 emission in most phases of their evolution, except when they become too cool. This line is present in the integrated light of population synthesis models down to 0.1 Zsun whatever the star formation history, and is also sometimes seen at Z=0.01 Zsun. HeII 1640 is weaker in models that include a Z scaling of the mass loss rates. The optical spectra of starbursts, especially the Wolf-Rayet bumps, sometimes display VMS signatures when these stars are present. At low Z, adding VMS to population synthesis models produces more ionising photons down to 45 eV. At higher energy the ionising flux depends on age, Z, assumption regarding VMS mass loss rates, and on the very short phases at the end of VMS evolution. HeII ionising fluxes large enough to produce some amount of nebular HeII 4686 emission can be produced under specific circumstances. Our models are able to reproduce qualitatively and sometimes also quantitatively the UV spectra of star-forming regions. However we are not able to clearly identify which mass loss framework is favoured.

Recent in situ observations from Parker Solar Probe (PSP) near perihelia reveal ion beams, temperature anisotropies, and kinetic wave activity. These features are likely linked to solar wind heating and acceleration. During PSP Encounter 17 (at $11.4R_s$) on Sep-26-2023, the PSP/FIELDS instrument detected enhanced ion-scale wave activity associated with deviations from local thermodynamic equilibrium in ion velocity distribution functions (VDFs) observed by the PSP/Solar Probe Analyzers-Ion (SPAN-I). Dense beams (secondary populations) were present in the proton VDFs during this wave activity. Using bi-Maxwellian fits to the proton VDFs, we found that the density of the proton beam population increased during the wave activity and, unexpectedly, surpassed the core population at certain intervals. Interestingly, the wave power was reduced during the intervals when the beam population density exceeded the core density. The drift velocity of the beams decreases from 0.9 to 0.7 of the Alfvén speed and the proton core shows a higher temperature anisotropy ($T_\perp/T_\parallel>2.5$) during these intervals. We conclude that the observations during these intervals are consistent with a reconnection event during a heliospheric current sheet crossing. During this event, $\alpha$ particle parameters (density, velocity, and temperature anisotropy) remained nearly constant. Using linear analysis, we examined how the proton beam drives instability or wave dissipation. Furthermore, We investigated the nonlinear evolution of ion kinetic instabilities using hybrid kinetic simulations. This study provides direct clues about energy transfer between particles and waves in the young solar wind.

The exact location of the $\gamma$-ray emitting region in blazar jets has long been a matter of debate. However, the location has important implications about the emission processes, geometric and physical parameters of the jet, as well as the nature of interaction of the jet with the interstellar and intergalactic medium. Diverse conclusions have been drawn by various authors based on a variety of methods applied to different data sets of many blazars, e.g., the location is less than 0.1 pc from the central engine within the broad line region (BLR) or a few or tens of pc downstream beyond the dusty torus or at some intermediate distance. Here we use a method, established in a previous work, in which the location of the GeV/optical emission is determined using the ratio of energy dissipated during contemporaneous outbursts at those wave bands. We apply it to a total of 47 multi-wavelength outbursts in 10 blazars. We find that the location of the GeV/optical emission is beyond the BLR for all cases. This result is consistent with other studies, in which the location has been determined for a large sample of blazars. We compare the location determined by our method for several GeV outbursts of multiple blazars to that obtained by other authors using different methods. We find that our results are consistent in such one-to-one comparison in most cases, for which the required data were available.

Ellerman bombs (EBs) are observable signatures of photospheric small-scale magnetic reconnection events. The reliable automatic detection of EBs would enable the study of the impact of magnetic reconnection on the Sun's dynamics. We aim to develop a method to automatically detect EBs in H$\alpha$ observations from the Swedish 1-m Solar Telescope (SST) and in SDO/AIA observations using the 1600Å, 1700Å, 171Å and 304Å passbands. We trained models based on neural networks (NNs) to perform automatic detection of EBs. Additionally, we used different types of NNs to study how different properties contribute to the detection of EBs. We find that for SST observations, the NN-based models are proficient at detecting EBs. With sufficiently high spectral resolution, the spatial context is not required to detect EBs. However, as we degrade the spectral and spatial resolution, the spatial information becomes more important. Models that include both dimensions perform best. For SDO/AIA, the models struggle to reliably distinguish between EBs and bright patches of different origin. Permutation feature importance revealed that the H$\alpha$ line wings (around $\pm$1 Å from line center) are the most informative features for EB detection. For the SDO/AIA case, the 1600Å channel is the most relevant one when used in combination with 171Å and 304Å. The combination of the four different SDO/AIA passbands is not informative enough to accurately classify EBs. From our analysis of a few sample SDO/AIA 1600Å and 1700Å light curves, we conclude that inclusion of the temporal variation may be a significant step towards establishing an effective EB detection method that can be applied to the extensive SDO/AIA database of observations. Abstract modified for ArXiv purposes.

Examining the spin-orbit coupling effects for highly elliptical orbits is relevant to the mission design and operation of cislunar space assets, such as the Lunar Gateway. In high-eccentricity orbits, the gravity-gradient moment is here modelled as an instantaneous excitation at each periapsis passage. By approximating the gravity-gradient moment through Dirac pulses, we derive a recursive discrete map describing the rotational state of the satellite at the periapsis passage. Thanks to the recursive map, we are able to find the initial attitude corresponding to an unbounded growth of angular velocity, and to identify initial conditions whose evolution is such that the pulses have the same sign (in-phase condition) or the alternate sign (counterphase condition) at successive periapsis passages. In the recursive map, we perform the numerical analysis up to ten periapsis passages. In order to justify the introduction of the discrete map, we compare the results of the discrete map with those found in the spin-orbit problem. Because of numerical errors due to the high eccentricity, we restrict the investigation up to three periapsis passages in the spin-orbit problem. Moreover, we apply the Fast Lyapunov Indicators method to draw the phase portrait and detect the initial conditions fulfilling the counterphase condition.

We argue that $^{23}$Na is a potentially significant source of galactic axions. For temperatures $\gtrsim 7 \times 10^8$K -- characteristic of carbon burning in the massive progenitors of supernovae and ONeMg white dwarfs -- the 440 keV first excited state of $^{23}$Na is thermally populated, with its repeated decays pumping stellar energy into escaping axions. Odd-A nuclear abundances are typically very low in high-temperature stellar environments (or absent entirely due to burn-up). $^{23}$Na is an exception: $\approx 0.1 M_\odot$ of the isotope is synthesized during carbon burning then maintained at $\approx 10^9$K for times ranging up to $6 \times 10^4$y. Using MESA simulations, a galactic model, and sampling over progenitor masses, locations, and evolutionary stages, we find a continuous flux at earth of $\langle \phi_a \rangle \approx 22$/cm$^2$s for $g^\mathrm{eff}_{aNN} = 10^{-9}$. Some fraction of these axions convert to photons as they propagate through the galactic magnetic field, producing a distinctive 440 keV line $\gamma$ ray detectable by all-sky detectors like the Compton Spectrometer and Imager (COSI). Assuming a 1$\mu$G galactic magnetic field and a sufficiently light axion mass, we find that COSI will be able to probe $| g_{aNN}^\mathrm{eff} g_{a \gamma \gamma} | \gtrsim1.8 \times 10^{-22}$ GeV$^{-1}$ at $3\sigma$ after two years of surveying.

How effective are methods for estimating bolide energies from infrasound signal period-yield relationships? A single global period-energy relation can obscure significant variability introduced by parameters such as the atmospheric Doppler wind profile and the bolide's energy deposition profile as a function of altitude. Bolide speed, entry angle, burst altitude, and multi-episode fragmentation all may play a role in defining the detected period of the shockwave. By leveraging bolide light curve data from the Center for Near Earth Object Studies (CNEOS), we re-examined the period-energy relation as a function of these parameters. Through a bootstrap approach, we show that various event subsets can deviate from widely cited period-energy models and we identify which specific conditions most strongly reshape the period-energy scaling. The results define both the fidelity and reliability of period-energy relations when no additional data beyond the infrasound record is available and improve the outcome when supporting data from bolide trajectories and light curves are included. Ultimately, these findings expand the scope of earlier models, providing a nuanced and robust framework for infrasound-only yield estimation under a range of bolide scenarios.

Type Ia Supernovae (SNe Ia) have become the most precise distance indicators in astrophysics due to their incredible observational homogeneity. Increasing discovery rates, however, have revealed multiple sub-populations with spectroscopic properties that are both diverse and difficult to interpret using existing physical models. These peculiar events are hard to identify from sparsely sampled observations and can introduce systematics in cosmological analyses if not flagged early; they are also of broader importance for building a cohesive understanding of thermonuclear explosions. In this work, we introduce DiTSNe-Ia, a variational diffusion-based generative model conditioned on light curve observations and trained to reproduce the observed spectral diversity of SNe Ia. In experiments with realistic light curves and spectra from radiative transfer simulations, DiTSNe-Ia achieves significantly more accurate reconstructions than the widely used SALT3 templates across a broad range of observation phases (from 10 days before peak light to 30 days after it). DiTSNe-Ia yields a mean squared error of 0.108 across all phases-five times lower than SALT3's 0.508-and an after-peak error of just 0.0191, an order of magnitude smaller than SALT3's 0.305. Additionally, our model produces well-calibrated credible intervals with near-nominal coverage, particularly at post-peak phases. DiTSNe-Ia is a powerful tool for rapidly inferring the spectral properties of SNe Ia and other transient astrophysical phenomena for which a physical description does not yet exist.

The Milky Way is a dynamic and evolving system shaped by numerous merger events throughout its history. These mergers bring stars with kinematic and dynamic properties differing from the main stellar population. However, it remains uncertain whether any of the Galactic supernova remnants can be attributed to such a merger this http URL this work, we compare the progenitor of Kepler's supernova to its nearby stars, ``alien'' stars, and in-situ Milky Way stellar populations. We uncover the abnormal kinematics and dynamics of Kepler's supernova and propose that its progenitor may have an extragalactic origin. We call the SNe Ia produced by stars accreted into the Milky Way through merger events ``Alien-Type-Ia supernovae'' (Alien-SNe-Ia) since they are cosmic immigrants. We estimate the rate of Alien-SNe-Ia exploded recently, using two methods: through galactic chemical evolution (using GalCEM), and through a method without considering exact star formation history, introduced for the first time in this paper. We consider the past accretion of a few major satellite galaxies -- Kraken, Gaia-Enceladus/Sausage, Helmi streams, Sequoia, Sagittarius, Wukong/LMS-1, and Cetus -- assuming these were dry mergers. The first method yields $1.5\times 10^{-5} - 5.0\times10^{-5}\rm\,yr^{-1}$, while the second method yields a comparable ${3.1}^{+1.8}_{-{1.1}}\times10^{-5}\rm\,yr^{-1}$ as the rate estimates for recent Alien-SNe-Ia. These estimates represent lower bounds because we assumed no post-merger star formation.

We present the morphological properties of 17 spectroscopically confirmed massive quiescent galaxies ($10.2 < \log(M_{\ast}/M_{\odot}) < 11.2$) at $3.0 < z < 4.3$, observed with JWST/NIRSpec and NIRCam. Using Sérsic profile fits to the F277W and F444W imaging, we derive the size-mass relation and find typical sizes of $\sim$0.6-0.8 kpc at $M_{\ast} = 5 \times 10^{10}~M_{\odot}$, consistent with a factor of $\sim$7 growth from $z \sim 4$ to the present day, including a factor of $\sim$2 by $z \sim 2$. Galaxies with earlier formation times ($z_{50} \gtrsim 5$) tend to be larger (up to 0.3 dex above the relation), while late-forming galaxies are more compact. Using a random forest regressor, we identify local environmental density, quantified by $\log(1+\delta^{\prime}_{3})$ from the distances to the three nearest neighbors, as the strongest predictor of bulge-to-total ratio ($B/T$), which ranges from 0.25 to 1. Applying the same method to the IllustrisTNG simulation, the ex-situ stellar mass fraction ($f_{\ast,\mathrm{ex\text{-}situ}}$) -- a proxy for mergers -- is the dominant predictor of $B/T$. Galaxies with high $B/T$ in dense environments show bursty star formation and short quenching timescales ($\lesssim0.4$ Gyr), consistent with bulge growth via merger-induced starbursts; in simulations, such galaxies exhibit elevated ex-situ fractions ($\sim20\%$-$30\%$). In contrast, some high-$B/T$ galaxies in intermediate-density environments have low ex-situ fractions, suggesting that additional mechanisms -- such as violent disk instabilities -- may contribute. These findings point to multiple bulge growth pathways at high redshift, unified by rapid gas accretion, central starbursts, and AGN feedback, as predicted by cosmological simulations.

For time-domain astronomy, it is crucial to frequently image celestial objects at specific depths within a predetermined cadence. To fulfill these scientific demands, scientists globally have started or planned the development of non-interferometric telescope arrays in recent years. Due to the numerous parameters involved in configuring these arrays, there is a need for an automated optimization framework that selects parameter sets to satisfy scientific needs while minimizing costs. In this paper, we introduce such a framework, which integrates optical design software, an exposure time calculator, and an optimization algorithm, to balance the observation capabilities and the cost of optical telescope arrays. Neural networks are utilized to speed up results retrieval of the system with different configurations. We use the SiTian project as a case study to demonstrate the framework's effectiveness, showing that this approach can aid scientists in selecting optimal parameter sets. The code for this framework is published in the China Virtual Observatory PaperData Repository, enabling users to optimize parameters for various non-interferometric telescope array projects.

We present a catalog of quasar pairs (QPs) constructed from the DESI DR1 quasar sample, which includes approximately 1.6 million spectroscopically confirmed quasars. Using a redshift-dependent self-matching procedure and applying physical constraints on projected separation (up to 110 kpc) and line-of-sight velocity difference (up to 2000 km/s), we identified 1,842 candidate quasar pairs. Each pair is spectroscopically confirmed, providing reliable redshift and velocity information. We visually classified these systems using DESI Legacy Imaging and SPARCL spectral data into four categories: QP (quasar pairs), QPC (quasar pair candidates), LQC (lensed quasar candidates), and Uncertain. The redshift distribution peaks at around z = 1 to 2.5, and 64.3 percent of QPs have line-of-sight velocity differences below 600 km/s, suggesting that many systems may be physically associated. This catalog provides a statistically meaningful sample for future studies of dual AGNs, gravitational lensing, and quasar clustering.

We present the multi-band photometric and low-resolution spectroscopic analysis of four W Ursae Majoris eclipsing binaries (EWs) - J080510.1+141528 (hereinafter as J0805a), J080516.3+143138 (hereinafter as J0805b), J143358.7+053953 (hereinafter as J1433), and J143458.4+054143 (hereinafter as J1434). The multi-band ground based photometric data are collected using the 1.3-m Devasthal Fast Optical Telescope (DFOT) while we also make use of TESS photometric observations. The spectroscopic analysis is based on the low-resolution observations by 4-m Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). The ephemeris of these systems are updated using the photometric data from TESS and other photometric surveys. The system J0805b shows a secular change in the orbital period with a period change rate of 4.2 ($\pm$ 0.1) $\times10^{-7}$ days per year while the orbital period change rate for J1433 is calculated as -1.1 ($\pm$ 0.1) $\times10^{-6}$ days per year. The mass-transfer rate for J0805b is found to be dM$_{1}$/dt = -1.56($\pm$0.07) $\times 10^{-6}$ M$_{\odot}$/year and dM$_{1}$/dt = -7.95($\pm$0.87) $\times 10^{-7}$ M$_{\odot}$/year for J1433. All the systems have inclination > 79$^{\circ}$ except J1433 which has inclination of 72.8$^{\circ}$. The mass-ratios (less massive to more massive component) for these targets are < 0.5. All the system except J0805b are A-subtype contact binaries. The absolute parameters of the systems are determined using GAIA DR3 parallax and reddening information. The LAMOST spectra are analyzed using spectral subtraction technique. A small excess emission is detected for J1433 in H$_{\alpha}$ and H$_{\beta}$ region. The systems are plotted on Hertzsprung-Russell (HR) diagram and compared with previously studied systems. The mass-ratio vs radius-ratio relation is also investigated for these systems.

With current and upcoming experiments on the horizon, the global 21-cm signal can open up new avenues for probing dark matter (DM) physics at redshifts that are otherwise inaccessible to other observables. This work investigates the effects of elastic scattering between DM and baryons on the global 21-cm signal in two distinct interacting DM (IDM) models: Coulomb-like and velocity-independent interactions. Our analysis incorporates key astrophysical parameters essential for accurately modeling the global signal, including star formation efficiency, escape fraction of ionizing photons, normalization of the X-ray luminosity, the number of Lyman-Werner photons emitted per stellar baryon, the minimum virial temperature of star-forming halos, as well as the IDM particle mass and cross section. We perform a Fisher analysis to forecast the sensitivity of four global 21-cm signal experimental scenarios as probes of DM-baryon scattering. We find that global signal experiments, even at the sensitivity of the current facilities such as EDGES and SARAS3, could improve existing cosmological and astrophysical constraints on DM-baryon scattering. Our results also highlight the degeneracies among the DM-baryon interaction cross section and astrophysical quantities. In particular, degeneracies between the IDM cross section and two astrophysical parameters, the minimum virial temperature, and Lyman-Werner photon production, can significantly impact the DM interaction inference. Conversely, the velocity-independent cross section is found to be insensitive to uncertainties in the X-ray luminosity. These findings underscore the necessity of accurately characterizing the uncertainties in astrophysical parameters to leverage the full potential of the 21-cm global signal experiments in probing IDM physics.

We present the results of our recent HI observations conducted on the edge-on galaxy NGC2683 using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). In comparison to previous observations made by the VLA, FAST has detected a more extensive distribution of HI. Particularly noteworthy is that the detections made by FAST extend approximately four times farther than those of the VLA in the vertical direction from the galactic plane. The total HI flux measured for NGC2683 amounts to $F_{\rm HI} = 112.1\,\rm{Jy\,km\,s^{-1}}$ (equivalent to a total HI mass of $M_{\rm HI} = 2.32 \times 10^9\,{\rm M_\odot}$), which is slightly higher than that detected by VLA. FAST has also identified three dwarf galaxies in close proximity to NGC2683, namely KK69, NGC2683dw1 (hereafter dw1), and NGC2683dw3$?$ (hereafter dw3$?$). dw3$?$ is situated within the extended HI distribution of NGC2683 in projection and lies near the tail of KK69 extending towards NGC2683. These observations suggest that dw3$?$ is likely a result of the accretion process from NGC2683 to KK69. Furthermore, FAST has detected three high-velocity clouds (HVCs), with complex B potentially undergoing accretion with NGC2683. Based on the model from Vollmer et al. 2016 and incorporating the HI halo component, we found that the model with the added HI halo aligns more closely with our FAST observations in NGC2683. The estimated mass of this HI halo is $3 \times 10^8\,{\rm M_\odot}$, constituting approximately 13% of the total HI mass of the galaxy. We suggest that the origination of this HI halo is more likely attributed to external gas accretion.

Magnetic fields play a significant role in stellar evolution. In the last few years, asteroseismology has enabled the measurement of strong magnetic fields $10^4$--$10^6\,\mathrm{G}$ in the cores of dozens of red giants, and is the only known way to directly measure internal stellar magnetic fields. However, current data are still interpreted assuming that these fields are too weak or too axisymmetric to affect the orientation of the pulsations (i.e., make the pulsations ``oblique''), rendering stronger field strengths beyond the reach of existing asteroseismic searches. We show that, even when an oblique pulsator is also stochastic (such as in a red giant with a strong non-axisymmetric magnetic field), geometric effects will cause the signal to contain frequency components which remain in perfect relative phase with each other. This perfect phase relationship persists even over timescales in which stochasticity erases absolute phase information. This perfect relative coherence is a distinctive observational signature of oblique pulsation that does not require a model for mode frequencies to search for. However, due to its dependence on phase, this effect will not be evident in the power spectral density alone, and phase information should be retained in order to detect it. Coherence-based searches for oblique pulsations may pave the way to measurements of magnetic fields of currently inaccessible strengths in red giants, as well as some main-sequence and compact pulsators.

A significant fraction of little red dots (LRDs) exhibit nearby extended emission of unknown origin. If physically associated with the LRD, this component may trace stellar emission from an off-centered host galaxy, neighboring companions, or nebular gas illuminated by the active nucleus. We investigate the detailed spectral energy distribution of the extended emission near four LRDs in the JWST UNCOVER and MegaScience surveys. We accurately decompose the extended emission from the dominant point source by simultaneously fitting the images in eight broad-band and nine medium-band filters. After considering both the results from photometric redshift fitting and the probability of galaxies at different redshift overlapping, we confirm that the off-centered blobs in three sources are physically associated with the LRDs, with two of them showing strong [\ion{O}{3}] $\lambda\lambda 4959,\,5007$ emission captured by the medium-band filters. While the spectral energy distributions of all three blobs can be modeled assuming star-forming galaxies with stellar mass $\sim 10^8\,M_{\odot}$, the exceptionally strong [\ion{O}{3}] emission of two sources is best interpreted as pure nebular emission from low-density ($n<10\, {\rm cm}^{-3}$), low-metallicity ($Z\approx 0.05\,Z_{\odot}$) gas photoionized by the ultraviolet radiation from the nearby LRD. Adopting LRD halo masses constrained by clustering measurements and theoretical considerations, we estimate a typical baryonic halo mass accretion rate of $\sim 2-9\, M_{\odot}\,{\rm yr}^{-1}$. If the halo accretion rate is sustained to $z=4$ and stars form with an efficiency of 10\%, the accreted gas would form a galaxy with stellar mass $\sim 10^9\,M_{\odot}$, potentially rendering them spatially resolved at lower redshift.

High-precision infrared spectroscopic measurements now enable detailed characterization of sub-Neptune atmospheres, potentially providing constraints on their interiors. Motivated by this, atmospheric models have been developed to explore chemical interactions between hydrogen-dominated atmospheres and possibly underlying magma oceans with various redox states. Recent models have predicted monosilane (SiH$_4$) as a potential atmospheric species derived from magma oceans in sub-Neptunes, but suggested that it is highly depleted in the observable atmospheric layers. Here, we propose that SiH$_4$ can persist throughout the atmospheres of sub-Neptunes with FeO-free reduced magma oceans by considering the dissolution of H$_2$O into the magma oceans, a factor not accounted for in previous models. We construct a one-dimensional atmospheric model to simulate the chemical equilibrium composition of hydrogen-dominated atmospheres overlying FeO-free magma oceans, incorporating H-O-Si chemistry. Our results show that the dissolution of H$_2$O enhances the SiH$_4$ molar fraction to levels of 0.1--10~\%, preventing it from reverting to silicates in the upper atmospheric layers. We find that SiH$_4$-rich atmospheres can exist across a broad parameter space with ground temperatures of 2000--6000~K and hydrogen pressures of 10$^2$--10$^5$~bar. We discuss that SiH$_4$-rich atmospheres could contain the other silanes but lack C-/N-/O-bearing species. The detection of SiH$_4$ in future observations of sub-Neptunes would provide compelling evidence for the presence of a rocky core with a reduced magma ocean. However, the accuracy of our model is limited by the lack of data on the non-ideal behavior and radiative properties of SiH$_4$, highlighting the need for further numerical and laboratory investigations.

We present a pilot method to estimate the high-mass initial mass function (IMF) across the arm, interarm, and spur regions in galaxies and apply it to NGC 628. We extracted star-forming complexes (SFCs) from H$\alpha$ VLT/MUSE and UVIT (FUV and NUV) observations of NGC 628 and used ALMA observations to define the molecular gas distribution. We find that the extinction-corrected H$\alpha$ and FUV luminosities correlate well. Using the fact that O stars have a shorter lifetime (10$^7$ yr) compared to B stars (10$^8$ yr), we estimated the approximate number of O stars from H$\alpha$ emission, and the number of B0 ($M_{*} > 10 M_{\odot}$), and B1 ($10 M_{\odot} \geq M_{*} \geq 3 M_{\odot}$) stars using FUV, NUV observations. We derived the IMF index ($\alpha$) for different regions using O to B0 ($\alpha_{1}$) and B0 to B1 ($\alpha_{2}$) stellar ratios. Our findings indicate that if we assume H$\alpha$ arises only from O8-type stars, the resulting $\alpha_{1}$ value is consistent with the canonical IMF index. It steepens when we assume O stars with masses up to 100 $M_{\odot}$ with mean $\alpha_{1}= 3.16 \pm 0.62$. However, the $\alpha_{2}$ does not change for large variations in the O-star population, and the mean $\alpha= 2.64 \pm 0.14$. When we include only blue SFCs ($ FUV-NUV\leq0.3$), mean $\alpha_{2}$ is $2.43 \pm 0.06$. The IMF variation for SFCs in arms and spurs is insignificant. We also find that $\alpha_{2}$ correlates with different properties of the SFCs, the most prominent being the extinction-corrected UV color (FUV-NUV).

Compact Symmetric Objects (CSOs) are a distinct category of jetted active galactic nuclei (AGN) whose optical variability characteristics have not been well investigated. We present here the results of our investigation on the optical flux and colour variability properties of a bona fide sample of 38 CSOs. We used the g-, r- and i-bands data from the Zwicky Transient Facility survey that spans a duration of about 5 years. We also considered a comparison sub-sample of blazars that includes 5 flat spectrum radio quasars and 12 BL Lac objects with redshifts and g-band magnitudes similar to the limited sub-sample of 9 CSOs. These two sub-samples of AGN, chosen for this comparative study of their long-term optical variability, represent different orientations of their relativistic jets with respect to the observer. We found that both CSOs and blazars exhibit optical flux variations, although variability of CSOs is lower than that of blazars. The observed variability in both CSOs and blazars is attributed to the relativistic jets and the increased optical variations in blazars relative to CSOs are likely due to beaming effects. CSOs and blazars exhibit similar colour variations, with both of them showing a bluer when brighter trend. Such a colour variability pattern is expected due to processes associated with their relativistic jets.

We used ALMA to perform a line survey of the high-mass star forming region Sgr B2(N), called ReMoCA. We modeled under the assumption of LTE the spectra obtained toward the sources embedded in the secondary hot core Sgr B2(N2). We compared the chemical composition of these sources to that of sources from the literature and to predictions of the chemical kinetics model MAGICKAL. We detected up to 58 molecules toward Sgr B2(N2)'s hot cores, including up to 24 COMs, as well as many less abundant isotopologs. The compositions of some pairs of sources are well correlated, but differences also exist in particular for HNCO and NH2CHO. The abundances of series of homologous molecules drop by about one order of magnitude at each further step in complexity. The nondetection of radicals yields stringent constraints on the models. The comparison to the chemical models confirms previous evidence of a high cosmic-ray ionization rate in Sgr B2(N). The comparison to sources from the literature gives new insight into chemical differentiation. The composition of most hot cores of Sgr B2(N2) is tightly correlated to that of the hot core G31.41+0.31 and the hot corino IRAS 16293-2422B after normalizing the abundances by classes of molecules (O-, N-, O+N-, and S-bearing). There is no overall correlation between Sgr B2(N2) and the shocked region G+0.693-0.027 also located in Sgr B2, and even less with the cold starless core TMC-1. The class of N-bearing species reveals the largest variance among the four classes of molecules. The S-bearing class shows in contrast the smallest variance. These results imply that the class of N-bearing molecules reacts more sensitively to shocks, low-temperature gas phase chemistry after non-thermal desorption, or density. The abundance shifts observed between the N- and O-bearing molecules may indicate how violently and completely the ice mantles are desorbed. [abridged]

A comprehensive understanding of solar coronal heating and charged particle acceleration remains one of the most critical challenges in space and astrophysical plasma physics. In this study, we explore the contribution of Alfvén waves, both in their kinetic (KAWs) and inertial (IAWs) regimes, to particle acceleration processes that ultimately lead to coronal heating. Using a kinetic plasma framework based on the generalized Vlasov-Maxwell model, we analyze the dynamics of these waves with a focus on the perpendicular components of the Poynting flux vectors and the net resonance speed of the particles. Our results show that both the magnitude and dissipation rate of the Poynting flux for KAWs and IAWs decrease with increasing electron-to-ion temperature ratio (T_e/T_i) and normalized perpendicular electron inertial length (c k_x / omega_pe). We evaluate the associated electric potentials and find that KAWs are significantly influenced in the high wavenumber (k_x rho_i) regime. IAWs, on the other hand, show a decrease in electric potential along the magnetic field and an increase across it when the perpendicular electric field (E_x) is enhanced. We also determine the net resonant speeds of particles in the perpendicular direction and show that these wave-particle interactions can efficiently heat the solar corona over large distances (R_Sun). Finally, we quantify the power transported by KAWs and IAWs through solar flux loop tubes, finding that both wave types deliver greater energy with increasing T_e/T_i and c k_x / omega_pe. These findings offer deeper insights into wave-driven heating and are relevant to solar wind and magnetospheric physics.

Blue horizontal-branch (BHB) stars are ideal tracers for mapping the structure of Galactic stellar halo. Traditionally, BHB sample stars are built from large-scale spectroscopic surveys utilizing their spectral features, however, the resulting sample sizes have been quite limited. In this paper, we construct a catalog of BHB stars based on synthetic colors $(u-v)_{0}$ and $(g-i)_{0}$ in SkyMapper photometric systems, which are convolved from Gaia XP spectra. A total of 49,733 BHB stars are selected from nearly the entire sky (excluding regions of low Galactic latitudes $|b| \le 8^{\circ}$ with heavy reddening), with a completeness and purity exceeding 90\%. Using member stars of globular clusters with precise distance determinations, we carefully calibrate the relationship between the $g$-band absolute magnitude and $(g-i)_{0}$, achieving a precision of 0.11\,mag, which corresponds to a 5\% uncertainty in distance. This relation is applied to derive distances for all BHB stars in the constructed sample. Given current capabilities of Gaia XP observations, the constructed BHB sample is primarily located within 20 kpc, enabling detailed mapping of the inner stellar halo. To extend this depth to the outer halo or even the edge of our Galaxy, we explore the potential of the Chinese Space Station Telescope (CSST) and its broad-band photometry for detecting BHB stars. Using mock data from synthetic spectra, we find that it is feasible to distinguish BHB stars from blue stragglers (BS) stars using CSST near-ultraviolet bands ($NUV, u$) photometry. Thanks to the deep limiting magnitude of CSST, its data will provide a groundbreaking perspective on our Galaxy, particularly regarding the outer halo, in an unprecedented volume.

The study of cyclic hydrocarbons is of utmost relevance in current astrochemical research, as they are considered to be among the most significant reservoirs of carbon in the interstellar medium. However, while unsaturated cyclic hydrocarbons with three, five, and six carbon atoms have been widely investigated, the highly strained antiaromatic cyclobutadiene, c-C4H4, still remains uncharted. Here, we employed high-level CCSD(T)-F12/cc-pVTZ-F12//B2PLYPD3/aug-cc-pVTZ theoretical calculations to analyze whether the cyano (CN), ethynyl (CCH), and hydroxy (OH) derivatives of c-C4H4 and its structural isomers butatriene (H2CCCCH2) and vinylacetylene (H2CCHCCH) can readily form via the gas-phase reaction: C4H4 + X -> C4H3X + H (where X = CN, CCH, and OH). For each system, we thoroughly explored the corresponding potential energy surfaces, identifying their critical points to enable a detailed analysis of the thermochemistry. Hence, we found various exothermic pathways for the formation of CN and CCH derivatives of butatriene and vinylacetylene, with no net activation barriers, while the formation of the OH derivatives is in general less favorable. Prior to the mechanistic study, we also analyzed the complete conformational panorama and stability of all the derivatives at the CCSD(T)-F12/cc-pVTZ-F12 level. Overall, c-C4H3CN and c-C4H3CCH emerge as particularly promising candidates for interstellar detection, provided that the parental c-C4H4 is present in the gas phase. These findings highlight the potential for detecting polar derivatives of c-C4H4 as indirect evidence of its presence in the ISM, as it appears to be "invisible" to radioastronomical observations. Also, this study underscores the need for future laboratory and theoretical efforts to characterize the spectroscopic properties of the proposed derivatives, paving the way for their eventual identification in space.

Satellites around large asteroids are preferentially found among those with the most rapid rotation and elongated shape. The taxonomic statistics are similarly skewed; in total, 13 asteroids larger than 100 km are known to have satellites, but none have been discovered among S-type asteroids. Previous modeling suggests that satellites could be generated by impacts, but spin and shape have never been tracked in models to relate collisional circumstances with those two observed properties concerning the primary. Here we show, by combining simulations of impacts into porous low-density asteroids, their subsequent disruption, reaccumulation and long-term satellite stability, a direct pathway for the formation of satellites. The immediate distortion and elongation of a rotating target body provides a launching point for some debris distinct from simple ballistic ejecta trajectories. The debris that are found to originate from the distorted long-axis is sourced primarily from 10-20 km below the surface and can be placed directly onto eccentric orbits with sufficiently large pericenter distances that avoid rapid re-impact. The specific energy and resultant total mass loss in satellite-forming collisions are not constraining, which explains the observed lack of correlation between asteroids with satellites and those that are part of large asteroid families.

Aims: The enigma of the missing baryons poses a prominent and unresolved problem in astronomy. Dispersion measures (DM) serve as a distinctive observable of fast radio bursts (FRBs). They quantify the electron column density along each line of sight and reveal the missing baryons that are described in the Macquart (DM-z) relation. The scatter of this relation is anticipated to be caused by the variation in the cosmic structure. This is not yet statistically confirmed, however. We present statistical evidence that the cosmological baryons fluctuate. Methods: We measured the foreground galaxy number densities around 14 and 13 localized FRBs with the WISE-PS1-STRM and WISE x SCOS photometric redshift galaxy catalog, respectively. The foreground galaxy number densities were determined through a comparison with measured random apertures with a radius of 1 Mpc. Results: We found a positive correlation between the excess of DM that is contributed by the medium outside galaxies (DM_cosmic) and the foreground galaxy number density. The correlation is strong and statistically significant, with median Pearson coefficients of 0.6 and 0.6 and median p-values of 0.012 and 0.032 for the galaxy catalogs, respectively, as calculated with Monte Carlo simulations. Conclusions: Our findings indicate that the baryonic matter density outside galaxies exceeds its cosmic average along the line of sight to regions with an excess galaxy density, but there are fewer baryons along the line of sight to low-density regions. This is statistical evidence that the ionized baryons fluctuate cosmologically on a characteristic scale of $\lesssim$6 Mpc.

Solar flares are among the most powerful and dynamic events in the solar system, resulting from the sudden release of magnetic energy stored in the Sun's atmosphere. These energetic bursts of electromagnetic radiation can release up to 10^32 erg of energy, impacting space weather and posing risks to technological infrastructure and therefore require accurate forecasting of solar flare occurrences and intensities. This study evaluates the predictive performance of three machine learning algorithms: Random Forest, k-Nearest Neighbors (KNN), and Extreme Gradient Boosting (XGBoost) for classifying solar flares into 4 categories (B, C, M, X). Using the dataset of 13 SHARP parameters, the effectiveness of the models was evaluated in binary and multiclass classification tasks. The analysis utilized 8 principal components (PC), capturing 95% of data variance, and 100 PCs, capturing 97.5% of variance. Our approach uniquely combines binary and multiclass classification with different levels of dimensionality reduction, an innovative methodology not previously explored in the context of solar flare prediction. Employing a 10-fold stratified cross-validation and grid search for hyperparameter tuning ensured robust model evaluation. Our findings indicate that Random Forest and XGBoost consistently demonstrate strong performance across all metrics, benefiting significantly from increased dimensionality. The insights of this study enhance future research by optimizing dimensionality reduction techniques and informing model selection for astrophysical tasks. By integrating this newly acquired knowledge into future research, more accurate space weather forecasting systems can be developed, along with a deeper understanding of solar physics.

This study presents a prototype D-band waveguide orthomode transducer (OMT) fabricated using chemically etched brass platelets. This method offers a fast, cost-effective, and scalable approach for producing waveguide OMTs above 100 GHz, making it well-suited for current and future Cosmic Microwave Background polarization experiments, where large focal planes with thousands of receivers are required to detect the faint primordial \textit{B}-modes. Chemical etching has already demonstrated its effectiveness in manufacturing corrugated feedhorn arrays with state-of-the-art performance up to 150 GHz. Here, we evaluate its applicability to more complex structures, such as OMTs. We designed a single OMT prototype operating in the 130-170 GHz range, fabricated by chemically etching 0.15 mm-thick brass plates, which were then stacked, aligned, and mechanically clamped. Simulations based on metrological measurements of the OMT profile predict return losses below $-$10 dB, isolation better than $-$30 dB, and transmission around $-$0.5 dB. The measured transmission and isolation, however, is around $-$1.5/$-$2 dB and $<-$20 dB, respectively. Further simulations show that the degradation in the transmission is related to defects and roughness along the etched profile ($\mathrm{RMS}\simeq$3 $\mu$m), which is a typical and unavoidable effect of chemical etching. The discrepancy in isolation, instead, could arise from a slight rotation ($\sim$3$^{\circ}$) of the polarization angle within the measurement chain. Our results show that chemical etching is a fast, low-cost, and scalable technique for producing waveguide OMTs with state-of-the-art performance in terms of return loss and isolation. However, at frequencies above 100 GHz the transmission coefficient may degrade due to the mechanical precision limitations of chemical etching.

Using observations from GMRT and FAST, we conducted multi-wavelength studies on PSR J1948+3540 and analyzed its intensity modulation characteristics in detail. We found that the intensity modulation of this pulsar exhibits broad low-frequency modulation features. The modulation frequency/period is time-dependent, but the dominant modulation component varies with the observing frequency. Specifically, at low frequencies, the modulation is dominated by the first half of the middle component, while at high frequencies, it is dominated by the second half of the middle component. Spectral analysis revealed that the intensities of the leading and trailing components vary with the observing frequency, but the middle component does not change significantly. Besides, the polarization analyses reveal that the peak of the radiation intensity is located in the latter half of the middle component, whereas the linear polarization is dominant in the former half. However, due to the low degree of linear polarization, the change of the dominant modulation component with the observed frequency is not caused by the variation in linear polarization. The phenomenon of the dominant modulation component varying with observing frequency has not been reported before and remains difficult to understand within the current theoretical framework.

Compact symmetric objects (CSOs) are intrinsically compact extragalactic radio sources that are thought to be the progenitors of classical radio galaxies. To date, evolutionary models have mainly focused on the formation and growth of high-power radio sources, leaving unanswered many questions related to low-power objects, whose relativistic jets are likely more prone to instabilities. We present a new sample of candidate low-power CSOs selected from the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. The main selection criteria are (i) a parsec-scale double radio morphology from archival Very Long Baseline Array (VLBA) images and (ii) a VLBA total flux density consistent with that from the FIRST survey, which rules out the presence of significant radio emission extending beyond the parsec scale. The final sample consists of 60 sources with radio luminosities between 10$^{24}$ and 10$^{27}$ W Hz$^{-1}$ at 1.4 GHz and projected linear sizes between 45 and 430 pc, which fill a region in the radio power-size plane that is currently underpopulated. We carried out VLBA observations at 4.98 GHz of a sub-sample of 20 sources among the brightest candidate CSOs with the aim of confirming their classification. We classify 12 sources as CSOs on the basis of their radio structure and spectral index distribution. In two out of the four CSOs with core identification, the asymmetry in the flux density of the outer components is in agreement with light travel time effects, and there is no evidence of jet-cloud interaction. If we assume a simplistic parametric model, most of the sources in the total sample have a jet power of $\sim 10^{44} - 10^{45}$ erg s$^{-1}$, making their evolutionary paths sensitive to the individual conditions of the jet and its environment.

We present the result of a proof-of-concept simulation designed to classify LEO satellites based on their occultations of background stars. We generate satellite shapes drawn from two broad shape classes, 'boxwing' and 'square'. We then simulate the resulting occultation photometry that would be caused by these satellites orbiting in LEO and intersecting with background stars. The resulting data is then inverted to attempt to recover the input shape and classify the satellite correctly. We find that the technique is theoretically sound, but ambitious with current telescope capabilities. We construct an equation for the required success rate of the method, as a function of exposure time and density of background stars. We find that successful classification requires short exposure times and high background stellar densities. For success rates in excess of 75%, we find a required exposure time of $\sim2.0\times10^{-3}$s, and $\sim500$ stars along the satellites path. Results are presented in terms of these two key parameters, and are discussed in the context of current observational capabilities and alternative satellite characterisation methods.

State-of-the-art adaptive optics (AO) systems perform non-linear Fourier-type wavefront sensing for real-time corrections of dynamic wavefront aberrations. This general class of sensors uses a filtering mask in the focal plane that converts phase fluctuations of the incoming light into intensity variations in the subsequent pupil plane. Due to their high sensitivity, Fourier-type wavefront sensors (WFSs) are the sensors of choice for many current and upcoming AO systems in ophthalmic imaging, free-space optical communications (FSOC) and astronomical ground-based telescopes such as the forthcoming generation of extremely large telescopes (ELTs). Conventionally, linear methods, like a matrix-vector-multiplication (MVM), are used for the inversion of Fourier-type WFSs. However, their non-linear behavior gives rise to severe performance degradations when significant channel perturbations are observed. They are expected to occur during strong atmospheric turbulence conditions, which are commonplace for non-rural sites and daytime observations. This study presents a novel type of iterative reconstructor to overcome non-linear wavefront sensing regimes. The underlying method is the non-linear Landweber iteration with Nesterov acceleration, well known in the field of inverse problems. A significant advantage of the new approach is its direct applicability to any Fourier-type WFS. This is implemented by adapting the filtering mask of the specific Fourier-type WFSs in the model-based algorithm. Several Fourier-type wavefront sensors are considered for ELT-scale instruments and their performance with the new algorithm is compared. The study goes on to concentrate on the pyramid wavefront sensor (PWFS), one of the most well-known Fourier-type WFSs. We demonstrate in end-to-end simulations that this novel approach outperforms linear methods in non-linear sensing regimes.

Deuterated molecules can be used to study the physical conditions and the astro-chemical evolution of molecular clouds. large-sample surveys for deuterated molecules are needed to understand the enhancement of deuterated molecules from diffuse molecular gas to cold cores. A single-pointing survey toward the 559 Planck cold clumps of the Early Cold Core Catalogue (ECC) has been conducted using the Arizona Radio Observatory 12-meter telescope, focusing on the $J$=1-0 transitions of DCO$^+$ and DCN. The survey included observations of 309 cores for DCO$^+$ and DCN 1-0 simultaneously, followed by 71 of these cores where DCO$^+$ 1-0 was detected for H$^{13}$CO$^+$ and H$^{13}$CN 1-0 simultaneously, aiming to determine the deuterated fraction ($D_{\rm frac}$). Additionally, 250 cores were observed for DCO$^+$, DCN, H$^{13}$CO$^+$ and H$^{13}$CN 1-0 simultaneously. Among the 309 sources, DCO$^+$ and DCN 1-0 were detected in 79 and 11 sources, with a detection rates of 25.6 % and 3.6 % respectively. In the 250 sources observed for all four species, DCO$^+$, DCN, H$^{13}$CO$^+$ and H$^{13}$CN 1-0 were detected in 58, 9, 57 and 13 sources, with a detection rate of 23.2 %, 3.6 %, 22.8 % and 5.2 % respectively. The $D_{\rm frac}$(HCO$^+$) values in 112 sources range from 0.89 % to 7.4 % with a median value of 3.1 %, while $D_{\rm frac}$(HCN) values in 11 sources range from 1.5 % to 5.5 % with a median value of 2.3 %. The line widths of DCO$^+$ and H$^{13}$CO$^+$ 1-0 detections are mostly within 1 km s$^{-1}$. The similarity in $D_{\rm frac}$ values between HCO$^+$ and HCN indicates that the higher detection rate of DCO$^+$ 1-0 compared with DCN 1-0 is due to the lower critical density of DCO$^+$ 1-0. We suggest that the enhancement of DCO$^+$ and DCN likely begins in the early diffuse stage of the molecular cloud, rather than during the cold core formation stage.

Astronomical archives contain vast quantities of unexplored data that potentially harbour rare and scientifically valuable cosmological phenomena. We have systematically searched approximately 100 million image cutouts from the entire Hubble Legacy Archive using the recently developed AnomalyMatch method, which combines semi-supervised and active learning techniques for the efficient detection of astrophysical anomalies. This comprehensive search rapidly uncovered a multitude of astrophysical anomalies presented here that significantly expand the inventory of known rare objects. Among our discoveries are 138 new candidate gravitational lenses, 18 jellyfish galaxies, and 417 mergers or interacting galaxes. The efficiency and accuracy of our iterative detection strategy allows us to trawl the complete archive within just 2-3 days, highlighting its potential for large-scale astronomical surveys. We present a detailed overview of these newly identified objects, discuss their astrophysical significance, and demonstrate the considerable potential of AnomalyMatch to efficiently explore extensive astronomical datasets, e.g., upcoming Euclid data releases.

Distance measurements to extragalactic systems that are both accurate and precise are cornerstones of modern astrophysics, underpinning the calibration of standard candles and the determination of the Hubble constant. Dwarf galaxies, such as Wolf-Lundmark-Melotte (WLM) and Sextans A, provide valuable laboratories for testing distance scales across different stellar populations. In this work, we utilize the high sensitivity and spatial resolution of the James Webb Space Telescope (JWST) to measure the distances to WLM and Sextans A using the tip of the red giant branch (TRGB) method. Adopting the TRGB absolute magnitude calibrated by NGC 4258, we determine distance moduli of $\mu_{\mathrm{0,WLM}} = 24.977 \pm 0.018 (\mathrm{stat}) \pm 0.056 (\mathrm{sys})$ mag for WLM and $\mu_{\mathrm{0,SexA}} = 25.740 \pm 0.011 (\mathrm{stat}) \pm 0.057 (\mathrm{sys})$ mag for Sextans A. Our results are consistent within a 3% distance uncertainty with previous measurements based on TRGB, Cepheids, and J-Region Asymptotic Giant Branch (JAGB) methods. With improved distance measurements in the future, these two galaxies have the potential to serve as additional anchor points for TRGB calibration, aiming to reduce the TRGB-based distance uncertainty to below 2%.

Flux emergence is crucial for the formation of solar active regions and triggering of various eruptions. However, the detailed mechanisms by which flux emergence drives these eruptions remain unclear and require numerical investigation. Using 2.5-dimensional magnetohydrodynamic simulations, we investigate the interaction between emerging flux and background magnetic fields and dynamics of the induced eruptions. By systematically varying the strength and angle of the emerging magnetic field relative to the background field, we investigate its impact on the initiation and evolution of solar eruptions. The simulations show that magnetic reconnection between the emerging flux and background field drives the formation of current sheets, magnetic islands and multithermal jets. Stronger magnetic fields result in earlier eruptions, more energetic jets, and enhanced heating. The formation and ejection of magnetic islands affect the structure and dynamics of the jet. When the hot and cool components of jets reach the other footpoint of magnetic loops, they will generate spicules near the transition region. Varying the angle between the emerging flux and the background field, we find that larger angles delay filament ascent and eruption timing but facilitate filament formation. Filaments form a hot shell and oscillate with a period of 10 minutes driven by periodic plasma ejections. Repetitive reconnection events inject cold plasma into the self-consistently formed filament channel, introducing a new prominence formation mechanism by flux-emergence-fed injection. Our analysis highlights the dynamic interplay between magnetic reconnection, plasma cooling and heating, and filament dynamics. These findings provide insights into solar eruptions and their observational signatures, emphasizing the role of multi-thermal structures in the corona.

This study employs stellar parameters from spectroscopic surveys and Zhang et al. based on Gaia XP spectra, along with photometric data from GALEX, UVOT, and Gaia, to obtain extensive ultraviolet color excess information for the relevant bands of GALEX and UVOT. By considering the impact of stellar parameters (\teff, \feh, and \logg) on intrinsic color indices, and utilizing the blue-edge method combined with a random forest algorithm, an empirical relationship between stellar parameters and intrinsic ultraviolet color indices is established. By combining observed color indices, the study derives color excesses for 11,624,802 and 65,531 stars in the GALEX/near-UV and far-UV bands, and 336,633, 137,739, and 253,271 stars in the UVOT/\uvwa, \uvmb, and \uvwb\ bands, constructing corresponding ultraviolet extinction maps. Notably, the color excess data for the GALEX/near-UV band shows a tenfold increase from previous results, with the extinction map covering approximately two-thirds of the sky at a resolution of 0.4$^{\circ}$. The typical uncertainties in the ultraviolet color excesses are 0.21 mag, 0.30 mag, 0.19 mag, 0.24 mag, and 0.21 mag for \enb, \efb, \ewab, \embb, and \ewbb, respectively. By comparing the spatial distributions of \rv\ derived from ultraviolet and Gaia optical band measurements with those obtained from infrared and optical data in previous works, it is evident that the \rv\ distributions based on the ultraviolet data show noticeable differences, with some regions even exhibiting opposite trends. This suggests that a single-parameter \rv\ extinction law may not be sufficient to simultaneously characterize the extinction behavior across infrared, optical, and ultraviolet bands.

Aiming at a general description of four basic solutions describing the accretion processes, i.e., the Shakura-Sunyaev thin disc (SSD), the Shapiro-Lightman-Eardley solution (SLE), the slim disc, and the advection-dominated accretion flow (ADAF), we present generalized axisymmetric height-averaged equations, where the entropy advection, the radiation pressure, and photon trapping effect are all included self-consistently. Our generalized solution can reproduce the ADAF, SLE, SSD, and slim disc branches in a wide range of accretion rates from sub- to super-Eddington accretion. An S-curve in the $\dot{m}-\Sigma$ plane is also reproduced, representing the SSD branch, the radiation-pressure-dominant branch, and the slim disc branch. The solution gives a natural transition between SSD and slim disc when photon trapping occurs in the accretion flow, producing a radially hybrid SSD-slim disc structure in a certain range of accretion rates. The coexistence of ADAF and SSD below a critical accretion rate is clearly shown with distinct advection fraction of accretion energy. We also present the luminosity, radiation efficiency, and spectrum from the generalized solutions for a large range of accretion rates in stellar mass black holes.

We study the Faraday rotation measure (RM) variability of flat spectrum radio quasar 1633+382 on 5 epochs spanning from 2004 to 2008. We used 4 to 43~GHz VLBI polarization data from VLBA. Core RM across 4 to 15 GHz scales with a power index $a\sim2$, for $a$ in RM$\propto \nu^a$. We detected sign changes across epochs, both in the core and in the jet region. RM time variability in the core and jet region are not correlated, hence limiting the size of a possible Faraday screen. We relate the core RM variability to a new component emerging from the core region. For the jet region, we consider the jet-medium interaction to be a less likely cause of the RM variability because of the uniform spectral index distribution. The observed RM value variation requires a huge fluctuation in electron density or magnetic field, hence a foreground Faraday screen is less favoured. We further discuss other possibilities of the RM variability based on jet kinematics.

Rocky exoplanets accessible to characterisation often lie on close-in orbits where tidal heating within their interiors is significant, with the L 98-59 planetary system being a prime example. As a long-term energy source for ongoing mantle melting and outgassing, tidal heating has been considered as a way to replenish lost atmospheres on rocky planets around active M-dwarfs. We simulate the early evolution of L 98-59 b, c and d using a time-evolved interior-atmosphere modelling framework, with a self-consistent implementation of tidal heating and redox-controlled outgassing. Emerging from our calculations is a novel self-limiting mechanism between radiative cooling, tidal heating, and mantle rheology, which we term the 'radiation-tide-rheology feedback'. Our coupled modelling yields self-limiting tidal heating estimates that are up to two orders of magnitude lower than previous calculations, and yet are still large enough to enable the extension of primordial magma oceans to Gyr timescales. Comparisons with a semi-analytic model demonstrate that this negative feedback is a robust mechanism which can probe a given planet's initial conditions, atmospheric composition, and interior structure. The orbit and instellation of the sub-Venus L 98-59 b likely place it in a regime where tidal heating has kept the planet molten up to the present day, even if it were to have lost its atmosphere. For c and d, a long-lived magma ocean can be induced by tides only with additional atmospheric regulation of energy transport.

The nuclear regions of active galactic nuclei (AGNs) likely host clumpy clouds that occasionally obscure the central X-ray source, causing eclipse events. These events offer a unique opportunity to study the properties and origins of such clouds. However, these transient events are rarely reported due to the need for extensive, long-term X-ray monitoring for years. Here, we conduct a systematic search for eclipse events in 40 AGNs well-monitored by the Swift X-ray Telescope (XRT) over the past 20 years, comprising a total of $\sim$11,000 observations. Our selection criteria rely on significant variations in X-ray flux and spectral shape. We identify 3 high-confidence events in 3 AGNs and 8 candidates in 6 AGNs, all in type I AGNs. The observed clouds have column densities of $N_{\rm H}~\sim$ (0.2 $-$ 31.2) $\times$ 10$^{22}$ cm$^{-2}$ and ionization degrees of log $\xi$ $\sim$ (-1.3 $-$ 2.2). For the 5 events with well-constrained duration, their distances from the central black hole range from (2.4 $-$ 179) $\times$ 10$^{4}$ $R_{\rm g}$, with 2 clouds near the dust sublimation zone, 2 farther out. Interestingly, we find tentative correlations between the cloud properties (i.e. ionization state and column density) and the black hole mass and mass accretion rate, implying their strong connection to the accretion process, potentially via outflows. Our study also demonstrates the potential of the new X-ray all-sky monitor, Einstein Probe, in providing more detection and physical constraints for such events.

Effective Field Theory (EFT) modeling is expected to be a useful tool in the era of future higher-redshift galaxy surveys such as DESI-II and Spec-S5 due to its robust description of various large-scale structure tracers. However, large values of EFT bias parameters of higher-redshift galaxies could jeopardize the convergence of the perturbative expansion. In this paper we measure the bias parameters and other EFT coefficients from samples of two types of star-forming galaxies in the state-of-the-art MilleniumTNG and Astrid hydrodynamical simulations. Our measurements are based on the field-level EFT forward model that allows for precision EFT parameter measurements by virtue of cosmic variance cancellation. Specifically, we consider approximately representative samples of Lyman-break galaxies (LBGs) and Lyman-alpha emitters (LAEs) that are consistent with the observed (angular) clustering and number density of these galaxies at $z=3$. Reproducing the linear biases and number densities observed from existing LAE and LBG data, we find quadratic bias parameters that are roughly consistent with those predicted from the halo model coupled with a simple halo occupation distribution model. We also find non-perturbative velocity contributions (Fingers of God) of a similar size for LBGs to the familiar case of Luminous Red Galaxies. However, these contributions are quite small for LAEs despite their large satellite fraction values of up to $\sim 30\%$. Our results indicate that the effective momentum reach $k_{\rm{Max}}$ at $z=3$ for LAEs (LBGs) will be in the range $0.3-0.6 ~h\rm{Mpc}^{-1}$ ($0.2-0.8~h\rm{Mpc}^{-1}$), suggesting that EFT will perform well for high redshift galaxy clustering. This work provides the first step toward obtaining realistic simulation-based priors on EFT parameters for LAEs and LBGs.

We measure the true obliquity of TOI-2364, a K dwarf with a sub-Saturn-mass ($M_p = 0.18\,M_J$) transiting planet on the upper edge of the hot Neptune desert. We used new Rossiter-McLaughlin observations gathered with the Keck Planet Finder to measure the sky-projected obliquity $\lambda = 7$$^\circ$$^{+10^\circ}_{-11^\circ}$. Combined with a stellar rotation period of $23.47\pm0.29$ days measured with photometry from the Tierras Observatory, this yields a stellar inclination of $90^{\circ} \pm 13^{\circ}$ and a true obliquity $\psi = 15.6$$^\circ$$^{+7.7^\circ}_{-7.3^\circ}$, indicating that the planet's orbit is well aligned with the rotation axis of its host star. The determination of $\psi$ is important for investigating a potential bimodality in the orbits of short-period sub-Saturns around cool stars, which tend to be either aligned with or perpendicular to their host stars' spin axes.

High-amplitude $\delta$ Scuti stars (HADS) that pulsate in their first three radial modes are rare in current samples. Here, we analyse five such triple-mode HADS observed by the Transiting Exoplanet Survey Satellite (TESS) and report that the previously considered second overtone mode ($f_2$) is actually the non-radial component of a resonating integration mode (RI mode, resulting from resonance between a radial p-mode and a non-radial p-g mixed mode), which shows significant amplitude and frequency variations over short timescales (approximately 20 days). This RI mode appears to be widespread among these stars. Notably, all five stars are in the post-main sequence evolutionary phase, actively crossing the Hertzsprung gap. These stars offer valuable insights into stellar evolution during the Hertzsprung gap, which is one of the most rapid evolutionary stages in a star's life.

The extreme-ultraviolet (EUV) brightenings identified by Solar Orbiter, commonly known as campfires, are the smallest transient brightenings detected to date outside active regions in the solar corona. We investigate the spatio-temporal distribution of a large ensemble of the finest-scale EUV transient brightenings observed by the Extreme Ultraviolet Imager (EUI) aboard Solar Orbiter. We perform a statistical analysis of the EUV brightenings by using quiet-Sun observations at the highest possible spatial resolution ever obtained by EUI. We use observations in the 17.4 nm passband of the High Resolution EUV Imager (HRIEUV) of EUI acquired during the closest perihelia of Solar Orbiter in 2022 and 2023. Solar Orbiter being at a distance 0.293 AU from the Sun, these observations have an exceptionally high image scale of 105 km, recorded at a fast cadence of 3 seconds. We use a wavelet-based automatic detection algorithm to detect and characterise the events of interest, and study their morphological and photometrical properties. We report the detection of smallest and shortest lived EUV brightenings to date in the quiet-Sun. The size and lifetime of the detected EUV brightenings appear power-law distributed down to a size of 0.01 Mm$^{2}$ and a lifetime of 3 seconds. In general their sizes lie in the range of 0.01 Mm$^{2}$ to 50 Mm$^{2}$ , and their lifetimes vary between 3 seconds and 40 minutes. We find an increasingly high number of EUV brightenings at smaller spatial and temporal scales. We estimate that about 3600 EUV brightenings appear per second on the whole Sun. The HRIEUV brightenings thus represent the most prevalent, localised and finest-scale transient EUV brightenings in the quiet regions of the solar corona.

The TRAPPIST-1 planetary system is observationally favorable for studying if planets orbiting M stars can retain atmospheres and host habitable conditions. Recent JWST secondary eclipse observations of TRAPPIST-1 c rule out a thick \ch{CO2} atmosphere but do not rule out atmospheric water vapor or its photochemical product, oxygen. Given the high expected escape rate, maintenance of atmospheric water vapor would require a present-day water source, such as volcanic outgassing. Here, we simulate water outgassing on the TRAPPIST-1 planets over a broad phase space based on solar system terrestrial bodies. We then apply two filters based on observation and geochemistry that narrow this phase space and constrain the plausible outgassing scenarios. For all seven TRAPPIST-1 planets, we find that the water outgassing rate is most likely $\sim$0.03x Earth's but has upper limits of $\sim$8x Earth's. The allowed range also implies low, Mars-like magma emplacement rates and relatively dry, Earth-like mantles, although mantle water mass fractions up to 1 wt\% are possible. We also present scenarios with magma emplacement rates similar to Mars, Earth, and Io, resulting in different preferred mantle water content and outgassing rates. We find that water outgassing rates are potentially high enough to balance water escape rates, providing a theoretical pathway for the TRAPPIST-1 planets to maintain surface water or water-vapor-containing atmospheres over long timescales. The bounds on outgassing rates and interior properties can be used in atmospheric chemistry and escape models to contextualize future observations of the TRAPPIST-1 planets, and may be applicable to other terrestrial exoplanets.

Planet migration within inner protoplanetary disks significantly influences exoplanet architectures. We investigate various migration mechanisms for young planets close to young stars. To quantify the stochastic migration driven by turbulent disks, we incorporate planets into existing 3-D MHD disk simulations of magnetospheric accretion. Besides the stochastic torque, we identify periodic torques from slowly evolving disk substructures farther out. We quantify these turbulent torques analytically using a modified Gaussian process. Then, using the disk structure in our simulation, we calculate migration timescales of various processes, including the smooth Type I/II migration, planet-star tidal interaction, magnetic dipole-dipole interaction, unipolar induction, and aerodynamical drag with the magnetosphere. Since our inner MHD turbulent disk reveals a very low surface density ($\sim 0.01$ g/cm$^{2}$), the resulting disk migration is significantly slower than previously estimated. Earth-mass planets have the migration timescale in the inner MHD turbulent disk exceeding the Hubble time, effectively stalling at the deadzone inner boundary ($R_{\mathrm{DZIB}}$). Only giant planets could migrate inward within the turbulent disk, and may stall at the magnetospheric truncation radius ($R_T$). A simplified planet population synthesis demonstrates that, at the end of the disk phase, all planets around solar-mass stars typically stall at $\lesssim$0.1 au since $R_T\sim R_{\mathrm{DZIB}}$. However, around 2 $M_{\odot}$ stars, higher-mass planets stall significantly closer to the star compared to low-mass planets, due to $R_T\ll R_{\mathrm{DZIB}}$. These results are consistent with recent observations on exoplanet demographics around different types of stars. Finally, turbulence in the low-density disk is unable to break the resonant planets, and thus young planets in resonances may be abundant.

We analyse the kinematical properties of the CHEX-MATE (Cluster HEritage project with XMM-Newton - Mass Assembly and Thermodynamics at the Endpoint of structure formation) galaxy cluster sample. [...] We derive cluster mass profiles for 75 clusters using the \textsc{MG-MAMPOSSt} procedure, which recovers the gravitational potential and the anisotropy profiles from line-of-sight velocities and projected positions of galaxy members. The standard NFW and the Burkert models with flatter cores than NFW both adequately fit the kinematic data, with only marginal statistical preference for one model over the other. An estimation of the mass bias $(1-B_1) = M^{SZ}_{500}/M^{M}_{500} $ is performed from the comparison with SZ-X-ray-calibrated mass estimates, resulting in a value of $ 0.54 \pm 0.11$ when four evidently disturbed clusters are removed from the sample. We assess the dynamical state of the clusters by inferring the Anderson-Darling coefficient $(A^2)$ and the fraction of galaxies in substructures ($f_\text{sub}$). Except for a few cases, we found relatively low values for $A^2$, suggesting that CHEX-MATE clusters are not too far from relaxation. Moreover, no significant trends emerge among $A^2,\,f_\text{sub}$ and the difference between the log-masses estimated by \textsc{MG-MAMPOSSt} and by SZ-X-ray. We study the concentration-mass relation for the sample; despite the large scatter, we observe signs of an increasing trend for large-mass clusters, in agreement with recent theoretical expectations. Finally, the analysis of radial anisotropy profiles of member galaxies - stacked in five bins of mass and redshift - reveals that orbits tend to be isotropic at the center and more radial towards the edge, as already found in previous studies. A slight trend of increasing radial orbits at $r_{200}$ is observed in clusters with larger velocity dispersion

We present results from a large spectroscopic survey of the nearest young association to the Sun, Upper Scorpius, conducted using 2dF/HERMES on the Anglo-Australian Telescope. We use spectroscopic youth criteria such as Li-equivalent widths to identify >1000 pre-main sequence (PMS) members across the region and measure radial velocities, combining these with Gaia EDR3 5-parameter astrometry to obtain 6D kinematic information. We separate confirmed PMS association members into distinct kinematic groups and measure expansion and rotation trends in each. We find strong evidence for asymmetric expansion in several groups and derive expansion timescales from the greatest rates of expansion in each group. We also trace the past motion of these groups using an epicycle approximation and estimate the time since their most compact configuration. These kinematic properties are compared to literature ages and the star formation history of Upper Scorpius is discussed. We find evidence that a scenario in which star formation in the subgroups of Upper Scorpius proceeded independently, either by self-instability or external feedback from Upper Centaurus-Lupus, is more likely than a recently proposed "cluster chain" scenario in which these subgroups have triggered each other.

Space weather events in planetary environments sourced from transient host star emissions, including stellar flares, coronal mass ejections, and stellar proton events, can substantially influence an exoplanet's climate and atmospheric evolution history. These time-dependent events may also affect our ability to measure and interpret its properties by modulating reservoirs of key chemical compounds and changing the atmosphere's brightness temperature. The majority of previous work focusing on photochemical effects used single-column models. Here, using three-dimensional (3D) general circulation models with interactive photochemistry, we simulate the climate and chemical impacts of stellar flare-sourced energetic particle precipitation. We examine their effects on synchronously rotating TRAPPIST-1e-like planets on a range of spatiotemporal scales. We find that sudden thermospheric cooling is associated with radiative cooling of NO and CO$_2$, and middle-to-lower atmospheric warming is associated with transient increases in infrared absorbers such as N$_2$O and H$_2$O. In certain regimes, these temperature changes are strongly affected by O$_3$ variability depending on the flare's spectral shape and energy. However, the role of O$_3$ in temperature change is reduced in the most extreme flaring scenarios explored in our simulations. In addition to effects on planetary climate and atmospheric chemistry, we also find that intense flares can energize the middle atmosphere, causing enhancements in wind velocities up to 40 m s$^{-1}$ in sub-stellar night-sides between 30 and 50 km in altitude. Our results suggest that successive, energetic eruptive events from younger stars may be pivotal factors in determining the atmosphere dynamics of their planets.

Weak gravitational lensing of the CMB has been established as a robust and powerful observable for precision cosmology. However, the impact of Galactic foregrounds, which has been studied less extensively than many other potential systematics, could in principle pose a problem for CMB lensing measurements. These foregrounds are inherently non-Gaussian and hence might mimic the characteristic signal that lensing estimators are designed to measure. We present an analysis that quantifies the level of contamination from Galactic dust in lensing measurements, focusing particularly on measurements with the Atacama Cosmology Telescope and the Simons Observatory. We employ a whole suite of foreground models and study the contamination of lensing measurements with both individual frequency channels and multifrequency combinations. We test the sensitivity of different estimators to the level of foreground non-Gaussianity, and the dependence on sky fraction and multipole range used. We find that Galactic foregrounds do not present a problem for the Atacama Cosmology Telescope experiment (the bias in the inferred CMB lensing power spectrum amplitude remains below $0.3\sigma$). For Simons Observatory, not all foreground models remain below this threshold. Although our results are conservative upper limits, they suggest that further work on characterizing dust biases and determining the impact of mitigation methods is well motivated, especially for the largest sky fractions.

We evolve for the first time in full general relativity a small, collisional N-body black hole cluster of arbitrary total mass M. The bound cluster is initially compact (radius R/M~10), stable, and consists of 25 equal-mass, nonspinning black holes. The dynamical interactions of compact objects in N-body clusters is of great interest for the formation of black holes in the upper mass gap as well as intermediate and supermassive black holes. These are potential sources of gravitational waves that may be detected by both current and future observatories. Unlike previous N-body Newtonian and post-Newtonian simulations, no "subgrid physics" is required to handle collisions and mergers. We can therefore confirm in full general relativity several predictions from these simulations and analytic estimates: the runaway growth of a large black hole via repeated mergers; spindown of the central black hole with increasing captures; the ejection of a black hole with a large asymptotic velocity due to a several-body interaction; and a regime where mergers occur primarily via direct collisions on highly eccentric orbits instead of quasicircular inspirals. We extract the gravitational wave signal and find it has several distinct features associated with the compact cluster regime. Our results suggest the signal is sufficiently loud that next generation observatories would likely be able to detect similar events across most of the observable universe. This work is a preliminary proof-of-principle study that we hope will open up a new arena for numerical relativity and the study of N-body compact systems.

Rotating black holes can generate boson clouds via superradiance when the boson's Compton wavelength is comparable to the black hole's size. In binary systems, these clouds can produce distinctive observational imprints. Recent studies accounting for nonlinearities induced by orbital backreaction suggest that if the binary forms at a large separation, resonance transitions can significantly deplete the cloud, minimizing later observational consequences except for very specific orbital inclinations. In this paper, we extend this framework to supermassive black hole binaries (SMBHBs), considering the influence of their astrophysical evolutionary histories. We find that, before entering the gravitational wave (GW) radiation stage, the additional energy loss channels can accelerate orbital evolution. This acceleration makes hyperfine resonant transitions inefficient, allowing a sufficient portion of the cloud to remain for later direct observations. We further discuss the ionization effects and cloud depletion occurring at this stage. Based on these theoretical insights, we explore how multi-messenger observations for SMBHBs can be utilized to detect the ionization effects of boson clouds by examining changes in the orbital period decay rate via electromagnetic measurements and variations in GW strain over a wide frequency band. Our findings reveal a complex dependence on the binary's total mass, mass ratio, and boson mass, emphasizing the significant role of astrophysical evolution histories in detecting boson clouds within binaries.

Particle dark matter and primordial black holes (PBH) might coexist with appreciable cosmic abundances, with both contributing to the observed dark matter density $\Omega_{\rm DM}$. Large populations of PBH (with $\Omega_{\rm PBH}\sim \Omega_{\rm DM}$) are tightly constrained for PBH heavier than $10^{-11} M_\odot$. However, large fractional abundances with $ f_{\rm PBH}\simeq \Omega_{\rm PBH}/\Omega_{\rm DM}\sim0.01$ are consistent with the limits on PBH for a wide range of PBH masses. Scenarios with significant populations of both particle dark matter and PBH are intriguing. Notably, if the particle dark matter has interactions with the Standard Model, new constraints arise due to pair-annihilations that are enhanced by the PBHs, resulting in dark matter indirect detection constraints on $f_{\rm PBH}$. Here we derive the bounds on mixed scenarios in which PBHs coexist with particle dark matter whose relic abundance is set via freeze-in (``FIMPs''). We show that while the restrictions on $f_{\rm PBH}$ are less constraining for FIMPs than WIMPs, modest bounds still arise for large classes of models. We examine both IR and UV freeze-in scenarios, including the case of ``superheavy'' particle dark matter with PeV scale mass.

Recent analyses that combine the latest data from the Atacama Cosmology Telescope (ACT) with cosmic microwave background observations by BICEP/Keck and Planck, together with the DESI baryon acoustic-baryonic-oscillation (BAO) measurements, have tightened the limits on inflationary scenarios. The joint data set yields a spectral index of primordial scalar perturbations $n_s = 0.9743 \pm 0.0034$ and an upper bound on the tensor-to-scalar ratio of $r < 0.038$. This slight upward shift in $n_s$ puts the previously favored Starobinsky model, and the conventional metric Higgs(-like) inflation--based on a quartic potential with a large non-minimal coupling in the Jordan frame--under tension with observations. In metric Higgs-like inflation the attractor behavior makes the predictions remarkably stable against higher-order operators, so modifying $n_s$ through such terms is difficult. In this paper, I show that adding a quadratic mass term to the Jordan-frame potential can raise $n_s$ and restore compatibility with the new data. I also discuss how this mass term can naturally arise from threshold effects in Higgs inflation.

The effects of axions on quark matter equation of state (EOS) is studied within the three flavor Nambu--Jona-Lasinio model and its effects on on the non-radial oscillations of neutron stars is investigated. Using such an EOS for quark matter with axions and a EOS for hadronic matter within the relativistic mean field (RMF) theory, we discuss the hadron-quark phase transition (HQPT) using the Gibbs construction. The EOS so obtained is used to investigate the structure of hybrid neutron star (NS)s. It is found that the presence of axions in the core of compact stars stabilizes hybrid NSs in agreement with modern astrophysical constraints. It is further observed that the quadrupolar fundamental modes ($f$-modes) for such hybrid NSs get substantial enhancements both due to a larger quark core in the presence of axions and from the hyperons as compared to a canonical nucleonic neutron stars.

Beamforming is a well-known technique to combine signals from multiple sensors. It has a wide range of application domains. This paper introduces the Tensor-Core Beamformer: a generic, optimized beamformer library that harnesses the computational power of GPU tensor cores to accelerate beamforming computations. The library hides the complexity of tensor cores from the user, and supports 16-bit and 1-bit precision. An extensive performance evaluation on NVIDIA and AMD GPUs shows that the library outperforms traditional beamforming on regular GPU cores by a wide margin, at much higher energy efficiency. In the 16-bit mode, it achieves over 600 TeraOps/s on an AMD MI300X GPU, while approaching 1 TeraOp/J. In the 1-bit mode, it breaks the 3 PetaOps/s barrier and achieves over 10 TeraOps/J on an NVIDIA A100 GPU. The beamforming library can be easily integrated into existing pipelines. We demonstrate its use for medical ultrasound and radio-astronomical instruments.

Universal cosmologies are exact solutions of 10d type IIA supergravity containing a 4d Friedmann-Lemaître-Robertson-Walker factor, that can also be repackaged as solutions of 4d models, i.e. as 4d consistent truncations. We extend the dynamical system analysis of universal cosmologies, beyond the case of a single exponential potential. For an open universe (negative 3d spatial curvature), these models generally possess many desirable features: parametric control of e-folds, late-time acceleration from potentials with steep exponentials (i.e. in accordance with swampland bounds), small string-loop and $\alpha'$-corrections, scale separation and/or absence of decompactification.

Improving the sensitivity of gravitational wave detectors is necessary to enrich scientific outcome of gravitational wave astronomy. Birefringence in test mass mirrors of gravitational wave detectors can become an important factor for both current and next-generation gravitational wave detectors to achieve improved performance. In-situ birefringence characterization can become an essential diagnostic tool for detector performance, and needs to be established. We report a possible in-situ birefringence characterization method and its experimental results with a tabletop setup. The scheme proposed and demonstrated in this paper can be used as a diagnostic tool in large-scale gravitational wave detectors. We also discuss possible technological developments toward implementation in future gravitational wave detectors.

Anomaly detection in large datasets is essential in fields such as astronomy and computer vision; however, supervised methods typically require extensive anomaly labelling, which is often impractical. We present AnomalyMatch, an anomaly detection framework combining the semi-supervised FixMatch algorithm using EfficientNet classifiers with active learning. By treating anomaly detection as a semi-supervised binary classification problem, we efficiently utilise limited labelled and abundant unlabelled images. We allow iterative model refinement in a user interface for expert verification of high-confidence anomalies and correction of false positives. Built for astronomical data, AnomalyMatch generalises readily to other domains facing similar data challenges. Evaluations on the GalaxyMNIST astronomical dataset and the miniImageNet natural-image benchmark under severe class imbalance (1% anomalies for miniImageNet) display strong performance: starting from five to ten labelled anomalies and after three active learning cycles, we achieve an average AUROC of 0.95 (miniImageNet) and 0.86 (GalaxyMNIST), with respective AUPRC of 0.77 and 0.71. After active learning cycles, anomalies are ranked with 71% (miniImageNet) to 93% precision in the 1% of the highest-ranked images. AnomalyMatch is tailored for large-scale applications, efficiently processing predictions for 100 million images within three days on a single GPU. Integrated into ESAs Datalabs platform, AnomalyMatch facilitates targeted discovery of scientifically valuable anomalies in vast astronomical datasets. Our results underscore the exceptional utility and scalability of this approach for anomaly discovery, highlighting the value of specialised approaches for domains characterised by severe label scarcity.

Aether Scalar-Tensor theory is a relativistic alternative gravity model that behaves like cold dark matter on cosmological scales while predicting the MOND force-law in astrophysical systems. The theory correctly predicts the cosmic microwave background and linear matter power spectra, and the mass discrepancies observed across the Universe. We derive and solve the equations governing neutron stars in Aether Scalar Tensor theory at first-order in slow rotation, finding that the theory predicts approximate universal relations between the moment of inertia and the compactness ($I$--$C$ relations) that differ from their general relativity counterparts. These relations may enable tests of Aether Scalar-Tensor theory using X-ray observations of pulsars and gravitational wave observations of binary neutron star mergers.

The chiral quark condensate of QCD, which spontaneously breaks the anomalous axial symmetry, gives rise to axionic type global string-wall systems. If a Peccei-Quinn type axion exists in the theory, the axionic strings are in general accompanied by winding of the QCD quark condensate. Depending on the axion model the winding can proceed either in the $\eta'$ or in the pion direction. This determines the structure of fermionic zero modes and the anomaly inflow which has important astrophysical consequences. We point out that $\eta'$ and pion string-wall systems exist in pure QCD, independently of the hidden axion. Strikingly, even if a hidden axion exists, the early cosmology can be entirely dominated by string-wall systems formed by the QCD quark condensate. We also discuss their role in the QCD phase transition and in heavy-ion physics.

We introduce a modification of the Starobinsky model in the form of an additional cubic Ricci scalar curvature term $\sim \alpha R^3$, scaled by a dimensionless parameter $\alpha$, with the resulting inflaton potential being the standard Starobinsky potential modified to first parametric order by an additive term. The resulting potential is identical to the potential obtained by a modification of the superpotential employed in the construction of the Starobinsky model in the framework of no-scale supergravity, thus, extending the correspondence between a class of no-scale supergravity models and modifications of the Starobinsky model. We analyze the inflationary predictions of the model and find that for $-4.2 \times 10^{-5} \lesssim \alpha \lesssim -1.9 \times 10^{-5}$, the modified Starobinsky model is in full agreement with the recent observational data from the Atacama Cosmology Telescope, for a range of $e$-folds, $N_\star = 50-60$.

We investigate the impact of a non-minimal Yukawa-like coupling between curvature and inflaton field within the \emph{spontaneous baryogenesis} background. We demonstrate that this coupling leads to a significant enhancement in particle production, even for small values of the coupling constant $\xi$. Assuming a perfectly homogeneous and isotropic universe during the reheating phase, we study the inflaton decay into fermion-antifermion pairs by means of a semiclassical approach, treating fermions as quantized fields and considering the inflaton and the Ricci scalar as classical quantities. We adopt the simplest approach in which the inflaton is minimally coupled to baryons, and non-minimally with gravity. In particular, we solve the equations of motion for the inflaton to first order in perturbation theory, with $\xi$ serving as perturbative parameter. Afterwards, we compute the difference in the number densities of baryons and antibaryons produced through the inflaton decay into fermion-antifermion pairs. We show that the non-minimal coupling term \emph{de facto} increases inflaton mass, letting fermion-antifermion decays be more probable, and thus enhancing the overall baryogenesis process. As a further outcome, we find that the non-minimal Yukawa coupling also leads to a renormalization of the inflaton mass and weakly influences the bounds over the gravitational constant. Finally, since the fermionic fields appear not to be mass eigenstates, we specialize the mass-mixing between them only. To this end, we thus include the effects of mass-mixing and cosmic expansion into our calculations. Physical consequences of baryon production are therefore explored.

We investigate black hole (BH) solutions embedded in a dark matter (DM) halo, modeled as extensions of the Schwarzschild metric. The DM density profile is constrained by Hubble Space Telescope data, stellar dynamics, and globular cluster (GC) measurements of the elliptical galaxy NGC 4649 (M60). Using this profile, we construct two distinct spacetime solutions characterized by the black hole mass ($M$), critical velocity ($V_c$), and core radius ($a$), all reducing to the Schwarzschild case in the limit $V_c=0$ and $a=0$. Our results show that the DM halo modifies essential BH features, such as the event horizon radius and spacetime curvature, as reflected in the Kretschmann scalar. We also derive an approximate analytical expression for the BH shadow radius, which increases slightly due to the halo's influence. Comparisons with two observational datasets further validate the analysis. Thermodynamic properties are examined across the two models. In the first, a generalized Smarr formula is obtained via two additional variables, $\gamma$ and $a$. The second incorporates halo effects through $V_c=0$ and $a=0$. These results underscore the role of DM in shaping both geometric and thermodynamic aspects of BHs.