Papers for Wednesday, Feb 26 2025

Star formation through the dynamical magnetized collapse remains an active area of astrophysical research. We carry out a comprehensive exploration on the magnetized gravitational collapse of a non-rotating self-gravitating initially spherically symmetric prestellar cloud core using two-dimensional nonideal magnetohydrodynamic simulations incorporating ambipolar diffusion and Ohmic dissipation. Our study encompasses a broader range of equations of state (EOSs) in the form of $P(\rho) \propto \rho^{\Gamma}$, with the aim of constraining the choice of EOSs for allowing star formation. Our results reveal that the collapse with a $\Gamma$ no stiffer than $4/3$, complemented by magnetized virial theorem, allows the dynamical contraction of the prestellar core to happen continuously where a central point mass forms and steadily builds up its mass from the infalling envelope, with a mass accretion rate of a scale of the order of $c_{\rm s}^3/G$. The choice of an isothermal EOS most naturally facilitates the collapse as a magnetic analog of the inside-out collapse. In addition to that, our study exhibits that the nonisothermal magnetized collapse models with a $\Gamma$ no stiffer than 4/3 qualitatively demonstrate similar infall features to those of an isothermal EOS. Furthermore, the collapse models with a $\Gamma$ stiffer than $4/3$ fail to ensure the sufficient cooling to allow the direct mass growth of the central point mass, thus delaying the infall. Our work can offer deeper insights in understanding the significance of EOSs on the magnetized gravitational collapse, enabling star formation.

Context. Massive contact binaries are both stellar merger and gravitational wave progenitors, but their evolution is still uncertain. An open problem in the population synthesis of massive contact binaries is the predicted mass ratio distribution. Current simulations evolve quickly to mass ratios close to unity, which is not supported by the sample of observed systems. It has been shown that modifying the near core mixing properties of massive stars can alter the evolution of contact binaries, but it has not been tested on a whole population. Aims. We implement a prescription of the convective core overshooting based on the molecular gradient. The goal of the implementation is to limit the rejuvenation efficiency of the accretor. We aim to investigate the effects of the reduced rejuvenation on the mass ratio distribution of massive contact binaries. Methods. We calculate a grid of 4896 models with the binary-evolution code MESA with our convective overshoot formulation, and we compare the simulations with the known observations of massive contact binaries. Results. We find that by limiting the core rejuvenation through the convective core overshoot, the predicted mass-ratio distribution shifts significantly to values away from unity. This improves the theoretical predictions of the mass ratios of massive contact binaries. Conclusions. The core rejuvenation of the components in massive contact binaries is a key parameter for their predicted mass ratio distribution. Establishing the rejuvenation efficiency with the contact binary population should therefore be possible. The sample size and uncertainties associated to the characterization of contact binaries however prevents us from doing so, and other methods like asteroseismology can place constraints on the rejuvenation process.

The observed star formation and wind outflow rates in galaxies suggest cold gas must be continually replenished via infalling clouds or streams. Previous studies have highlighted the importance of cooling-induced condensation on such gas, which enables survival, mass growth, and a drag force which typically exceeds hydrodynamic drag. However, the combined effects of magnetic fields, cooling, and infall remains unexplored. We conduct 3D magnetohydrodynamic (MHD) simulations of radiatively cooling infalling clouds and streams in uniform and stratified backgrounds. For infalling clouds, magnetic fields aligned with gravity do not impact cloud growth or dynamics significantly, although we see enhanced survival for stronger fields. By contrast, even weak transverse magnetic fields significantly slow cloud infall via magnetic drag, due to the development of strong draped fields which develop at peak infall velocity, before the cloud decelerates. Besides enhancing survival, long, slow infall increases total cloud mass growth compared to the hydrodynamic case, even if reduced turbulent mixing lowers the rate of mass growth. Streams often result in qualitatively different behavior. Mass growth and hence accretion drag are generally much lower in hydrodynamic streams. Unlike in clouds, aligned magnetic fields suppress mixing and thus both mass growth or loss. Transverse fields do apply magnetic drag and allow streams to grow, when the streams have a well-defined 'head' pushing through the surrounding medium. Overall, regardless of the efficacy of drag forces, streams are surprisingly robust in realistic potentials, as the destruction time when falling supersonically exceeds the infall time. We develop analytic models which reproduce cloud/stream trajectories.

Subhalo abundance matching (SHAM) is widely used for connecting galaxies to dark matter haloes. In SHAM, galaxies and (sub-)haloes are sorted according to their mass (or mass proxy) and matched by their rank order. In this work, we show that SHAM is the solution of the optimal transport (OT) problem on empirical distributions (samples or catalogues) for any metric transport cost function. In the limit of large number of samples, it converges to the solution of the OT problem between continuous distributions. We propose SHAM-OT: a formulation of abundance matching where the halo-galaxy relation is obtained as the optimal transport plan between galaxy and halo mass functions. By working directly on these (discretized) functions, SHAM-OT eliminates the need for sampling or sorting and is solved using efficient OT algorithms at negligible compute and memory cost. Scatter in the galaxy-halo relation can be naturally incorporated through regularization of the transport plan. SHAM-OT can easily be generalized to multiple marginal distributions. We validate our method using analytical tests with varying cosmology and luminosity function parameters, and on simulated halo catalogues. The efficiency of SHAM-OT makes it particularly advantageous for Bayesian inference that requires marginalization over stellar mass or luminosity function uncertainties.

We present optical-to-near infrared (NIR) photometry and spectroscopy of the Type Ia supernova (SN Ia) 2024epr, including NIR spectra observed within two days of first light. The early-time optical spectra show strong, high-velocity Ca and Si features near rarely-observed velocities at $\sim$0.1$c$, and the NIR spectra show a \CI\ "knee." Despite these high-velocity features at early times, SN~2024epr evolves into a normal SN Ia, albeit with stronger peak-light Ca absorption than other SNe Ia with the same light curve shape. Although we infer a normal decline rate, $\Delta m_{15}(B)=1.09\pm0.12$ mag, from the light-curve rise, SN 2024epr is a Branch "cool" object and has red early-time colors ($g-r\approx0.15$ mag at $-10$ days). The high velocities point to a density enhancement in the outer layers of the explosion, but thick-shell He-detonation models do not match the smoothly rising light curve or lack of He in our early-time NIR spectra. No current models (e.g., delayed detonation or thin He shell double detonation) appear to reproduce all of the observed properties. Such constraints are only possible for SN 2024epr from the earliest optical and NIR observations, highlighting their importance for constraining SN Ia models. Finally, we find several other SNe Ia with intermediate mass elements at $\sim$30\,000 km s$^{-1}$ within days after the explosion that evolve into otherwise normal SNe Ia at peak light, suggesting the early-time spectra of SNe Ia may hide a broad diversity of observational characteristics.

Relativistic collimated outflows, or jets, provide a crucial mode of active galactic nucleus feedback. Although the jets extract their energy from the black hole (BH) rotation, their effect on the BH spin is poorly understood. Because the spin controls radiative and mechanical BH feedback, lack of first-principles models for BH spin evolution limits our ability to interpret observations, including the recent LIGO-Virgo-KAGRA spin constraints. Particularly important are luminous disks, which rapidly grow and strongly torque their BHs. Jet-less and weakly magnetized standard luminous disks spin-up their BHs to near-maximum dimensionless spin, $a_{eq,NT}=0.998$. However, sufficient large-scale vertical magnetic flux can cause the inner disk to enter a magnetically arrested disk (MAD) state, whose jets can efficiently extract BH rotational energy and significantly spin-down the BH. Indeed, Lowell et al. 2024 found that non-radiative MADs spin-down their BHs to very low $a_\text{eq,MAD}^\text{nr}=0.07$. Moreover, their analytic model predicted that luminous MADs also spin-down their BHs to low $a_{eq,MAD}^{lum}\sim0.3-0.5$. To test this prediction, we perform 3D general relativistic (radiation) magnetohydrodynamic (GR(R)MHD) simulations of MADs across a wide range of BH spin ($-0.9\le{}a\le0.99$) and disk dimensionless thickness ($0.03\le{}h/r\le0.3$, which corresponds to Eddington ratio, $0.35\le{}\dot{m}/\dot{m}_\text{Edd}\le\infty$). We find that luminous MADs spin-down their BHs to a low universal equilibrium spin value, $a_{eq,MAD}^{lum}\simeq0.3$ for $0.03\le{}h/r\le0.1$. Moreover, we find evidence for quadratic convergence, $a_{eq,MAD}\simeq0.31-2.71(h/r)^2\to0.31$ as $h/r\to0$. We attribute this to disk thermodynamics becoming irrelevant as the cooling becomes more aggressive and magnetic forces start to dominate. We finish by discussing the astrophysical implications.

The UltraViolet EXplorer (UVEX) is a wide-field ultraviolet space telescope selected as a NASA Medium-Class Explorer (MIDEX) mission for launch in 2030. UVEX will undertake deep, cadenced surveys of the entire sky to probe low mass galaxies and explore the ultraviolet (UV) time-domain sky, and it will carry the first rapidly deployable UV spectroscopic capability for a broad range of science applications. One of UVEX's prime objectives is to follow up gravitational wave (GW) binary neutron star mergers as targets of opportunity (ToOs), rapidly scanning across their localization regions to search for their kilonova (KN) counterparts. Early-time multiband ultraviolet light curves of KNe are key to explaining the interplay between jet and ejecta in binary neutron star mergers. Owing to high Galactic extinction in the ultraviolet and the variation of GW distance estimates over the sky, the sensitivity to kilonovae can vary significantly across the GW localization and even across the footprint of a single image given UVEX's large field of view. Good ToO observing strategies to trade off between area and depth are neither simple nor obvious. We present an optimal strategy for GW follow-up with UVEX in which exposure time is adjusted dynamically for each field individually to maximize the overall probability of detection. We model the scheduling problem using the expressive and powerful mathematical framework of mixed integer linear programming (MILP), and employ a state-of-the-art MILP solver to automatically generate observing plan timelines that achieve high probabilities of kilonova detection. We have implemented this strategy in an open-source astronomical scheduling software package called the Multi-Mission Multi-Messenger Observation Planning Toolkit (M4OPT), on GitHub at this https URL.

The James Webb Space Telescope (JWST) has opened a window on many new puzzles in the early Universe, including a population of high-redshift star clusters with extremely high stellar surface density, suggesting unique star formation conditions in the Universe's early evolution. We study the formation and evolution of these first star clusters and galaxies using an AREPO cosmological simulation box designed to resolve the intricate environments of the smallest halos hosting Population III star clusters at $z \geq 12$. Our approach, which prioritizes baryonic structure identification through a friends-of-friends algorithm, provides new insights into early star cluster formation and delivers predictions directly relevant to observations. We investigate the dynamical properties of these first star clusters and use numerical and analytical methods to understand the populations of virialized and non-virialized systems. Our findings indicate that high-$z$ star clusters in a feedback-free regime can achieve extreme surface densities, consistent with the systems detected by JWST. These results imply that JWST may have the opportunity to uncover stellar systems at high redshift whose dynamical state preserves evidence of the hierarchical structure formation process.

Stars and their associated planets originate from the same cloud of gas and dust, making a star's elemental composition a valuable indicator for indirectly studying planetary compositions. While the connection between a star's iron (Fe) abundance and the presence of giant exoplanets is established (e.g. Gonzalez 1997; Fischer & Valenti 2005), the relationship with small planets remains unclear. The elements Mg, Si, and Fe are important in forming small planets. Employing machine learning algorithms like XGBoost, trained on the abundances (e.g., the Hypatia Catalog, Hinkel et al. 2014) of known exoplanet-hosting stars (NASA Exoplanet Archive), allows us to determine significant "features" (abundances or molar ratios) that may indicate the presence of small planets. We test on three groups of exoplanets: (a) all small, R$_{P}$ $<$ 3.5 $R_{\oplus}$, (b) sub-Neptunes, 2.0 $R_{\oplus}$ $<$ R$_{P}$ $<$ 3.5 $R_{\oplus}$, and (c) super-Earths, 1.0 $R_{\oplus}$ $<$ R$_{P}$ $<$ 2.0 $R_{\oplus}$ -- each subdivided into 7 ensembles to test different combinations of features. We created a list of stars with $\geq90\%$ probability of hosting small planets across all ensembles and experiments ("overlap stars"). We found abundance trends for stars hosting small planets, possibly indicating star-planet chemical interplay during formation. We also found that Na and V are key features regardless of planetary radii. We expect our results to underscore the importance of elements in exoplanet formation and machine learning's role in target selection for future NASA missions: e.g., the James Webb Space Telescope (JWST), Nancy Grace Roman Space Telescope (NGRST), Habitable Worlds Observatory (HWO) -- all of which are aimed at small planet detection.

High-velocity particles ( $v>v_\mathrm{esc}$) in the Milky Way are rare but nonetheless important to characterize due to their impact on dark matter (DM) direct detection experiments. We select halos similar in mass to the Milky Way in a large-volume dark matter simulation and measure the incidence of high-velocity particles, finding that an average fraction $\sim 1.3\times 10^{-5}$ of the DM particles have $\Delta v > 600$ km/s within 5-11 kpc of the halo centers. However, some systems have dramatically higher fractions. Milky Way-like systems with high-speed satellites can have high-velocity DM fractions of order $\sim 100\times$ higher than average. The environment also affects high-velocity DM fractions; massive nearby halos ( $>10^{13} M_{\odot}$) can boost high-velocity DM density by $\sim 10\times$, although there is little effect for nearby Andromeda-like systems. We confirm previous predictions from zoom-in simulations that the high-velocity particles in the Milky Way with heliocentric speeds $>700$ km/s primarily originate from the Large Magellanic Cloud, and provide a table of expected high-velocity DM densities at different heliocentric velocity thresholds.

We systematically perform long-term (millions of Schwarzschild time) axisymmetric viscous hydrodynamics simulations for tori around black holes in general relativity supposing the super Eddington accretion flow. The initial condition for the tori is modeled simply by the Fishbone-Moncrief torus with a constant specific angular momentum $j$ but with a wide variety of $j$. We find that for a given density profile, the fraction of the mass infall onto the black hole is approximately proportional to $j^{-1}$, indicating that only a minor fraction of the matter in the torus formed far from the black hole falls into the black hole while the majority is ejected with the typical average velocity of a few percent of the speed of light. We also find that the mass ejection is driven only outside $\approx 2\,r_\mathrm{ISCO}$ where $r_\mathrm{ISCO}$ is the areal radius of the innermost stable circular orbit around black holes, which depends strongly on the black hole spin. We derive an approximate fitting formula for the spin-dependence on the mass infall fraction as $\propto r_\mathrm{ISCO}^{0.7}$, which suggests that the rapid growth of supermassive black holes proceeded primarily by the accretion of the matter with the angular momentum counter-rotating with the black hole spin.

It has been a decade since the IceCube collaboration began detecting high-energy (HE) neutrinos originating from cosmic sources. Despite a few well-known individual associations and numerous phenomenological, observational, and statistical multiwavelength studies, the origin of astrophysical HE neutrinos largely remains a mystery. To date, the most convincing associations link HE neutrinos with active galactic nuclei (AGN). Consequently, many studies have attempted population-based correlation tests between HE neutrinos and specific AGN subpopulations (such as blazars). While some of associations are suggestive, no definitive population-based correlation has been established. This could result from either a lack of a population-based correlation or insufficient detection power, given the substantial atmospheric neutrino background. By leveraging blazar variability, we performed spatio-temporal blazar-neutrino correlation tests aimed at enhancing detection power by reliably incorporating temporal information into the statistical analysis. We used simulations to evaluate the detection power of our method under various test strategies. We find that: (1) with sufficiently large source samples, if 20% of astrophysical HE neutrinos originate from blazars, we should robustly observe $\sim$4$\sigma$ associations; (2) a counting-based test-statistic combined with a top-hat weighting scheme (rather than a Gaussian one) provides the greatest detection power; (3) applying neutrino sample cuts reduces detection power when a weighting scheme is used; (4) in top-hat-like weighting schemes, small p-values do not occur arbitrarily with an increase in HE neutrino error region size (any such occurrence is indicative of an underlying blazar-neutrino correlation).

Neutrinos are excellent probes of the inner structures of supernovas. However, an understanding of their dynamics remains incomplete, which is crucial for properly interpreting the detector data. The huge matter density inside a core-collapse supernova affects the self-coupling of the neutrinos through a geometrical four-fermion interaction induced by spacetime torsion generated by the fermions themselves. For normal mass hierarchy, the effect is negligible; for inverted mass hierarchy however, we find that the gravitational coupling can significantly alter the flavor dynamics. The possibility of constraining the said interactions through the relative abundance of different flavors of the neutrinos is discussed.

We show that the blue 18.3-minute variable object discovered in the Galactic disk by the OGLE-III survey and named OGLE-GD-WD-0001 is a pulsating pre-white dwarf of PG 1159 spectral type. With an effective temperature of about 160,000 K it is among the hottest known pulsators being located close to the blue edge of the GW Virginis instability strip. The long-term OGLE observations indicate that the object has a positive period change rate of about $5 \times 10^{-10}$ $s~s^{-1}$ and thus already contracts. There are no traces of a planetary nebula around this star.

We analyze a new set of 275 n-body calculations designed to place limits on the masses of the small circumbinary satellites in the Pluto-Charon system. Together with calculations reported in previous papers, we repeat that a robust upper limit on the total mass of the four satellites is roughly $9.5 \times 10^{19}$ g. For satellite volumes derived from New Horizons, this mass limit implies a robust upper limit on the bulk densities of Nix and Hydra, $\lesssim 1.7$ g cm$^{-3}$, that are comparable to the bulk density of Charon. Additional calculations demonstrate that satellite systems with mass $\lesssim 8.25 \times 10^{19}$ g are robustly stable over the current age of the Sun. The bulk densities of Nix and Hydra in these lower mass systems are clearly smaller than the bulk density of Charon. These new n-body results enable accurate measurements of eccentricity and inclination for Nix, Kerberos, and Hydra that agree well with orbital elements derived from numerical calculations with new HST and New Horizons state vectors. With these new state vectors, Styx has a 37% larger eccentricity and an 85% smaller inclination, which makes it more prone to gravitational perturbations from Nix.

New multi-wavelength Karl G. Jansky VLA observations of CKR02A, the compact radio source in the center of the compact HII region NGC 6334A, are presented. The observations were carried out in five epochs and included the frequency ranges 8.0 - 12.0 GHz (X-band), 18.0 - 26.0 GHz (K-band), and 29.0 - 37.0 GHz (Ka-band). The source is detected and resolved in all the observed epochs and in all bands. The source shows a clear arc-shaped structure consistent with a bow shock. The analysis of the spectral index maps indicates that its spectral index is $\alpha=-0.68\pm0.17$, suggesting that the emission is non-thermal. Two astronomical objects can explain the emission nature and morphology of the source: a colliding wind region of two massive stars or the bow shock of a massive runaway star. However, no massive stars are reported so far in the center of NGC 6334A, though its presence is also suggested by the free-free radio emission of the C-HII region itself. Using ancillary VLA data, we measured a preliminary proper motion of $19\pm6$ mas yr$^{-1}$, equivalent to a velocity of $120\pm40$ km s$^{-1}$. A detailed discussion of the implications of both scenarios is provided. Finally, a list of compact radio sources in the vicinity of NGC 6334A is given and briefly discussed.

Active galactic nuclei (AGN) are known to be variable sources across the entire electromagnetic spectrum, in particular at optical/ultraviolet and X-ray energies. Over the past decades, a growing number of AGN have displayed type transitions: from type 1 to type 2 or viceversa within a few years or even several months. These galaxies have been commonly referred to as changing-look AGN (CLAGN). Here we report on a new CLAGN, NGC 4614, which transitioned from a type 1.9 to a type 2 state. NGC 4614 is a nearly face-on barred galaxy at redshift $z = 0.016$, classified as a low-luminosity AGN. Its central black hole has a mass of about $1.6\times 10^7 M_\odot$ and an Eddington ratio around 1 percent. We recently acquired optical spectra of NGC 4614 at the Telescopio Nazionale Galileo and the data clearly suggest that the broad H$\alpha$ component has strongly dimmed, if not disappeared. A very recent Swift observation confirmed our current optical data, with the AGN weakened by almost a factor of 10 with respect to previous X-ray observations. Indeed, NGC 4614 had been also observed by Swift/XRT 6 times in 2011, when the source was clearly detected in all observations. By fitting the stack of the 2011 Swift observations we obtain a photon index of $\Gamma=1.3\pm0.3$ and an equivalent hydrogen column density of $N_{\rm H}$=$1.2\pm0.3$ $\times$10$^{22}$ cm$^{-2}$, indicating that NGC 4614 can be moderately absorbed in the X-rays. Although a significant change in the foreground gas absorption that may have obscured the broad line region cannot be entirely ruled out, the most likely explanation for our optical and X-ray data is that NGC 4614 is experiencing a change in the accretion state that reduces the radiative efficiency of the X-ray corona.

This paper describes a method for parametric radial profile modelling of radio interferometric visibility data. Image-based parametric modelling is common in the field of circumstellar debris disks, and high resolution ALMA observations make this method computationally expensive because many high resolution images must be generated and Fourier transformed. Most circumstellar disks are axisymmetric, so here a method that avoids generating images and instead models radial profiles is outlined. The two main elements are i) $u,v$ space visibility averaging, and ii) using the Discrete Hankel Transform to convert arbitrary parametric radial surface brightness profiles to visibilities. Vertical scale height can be included, but is fixed with radius; in principle this is a limitation but in practise radially-dependent vertical structure information is beyond even ALMA's capability for all but a few debris disks. Combined, these ingredients yield a method that is one to two orders of magnitude faster than image based methods and can be run on a laptop, is sufficient for most ALMA observations of debris disks, and can be used with gradient-based fitting methods such as Hamiltonian Monte-Carlo. A small software package is provided, which is dubbed vis-r.

We report the first detections of [Ne V] 14.3 {\mu}m and [Ne VI] 7.7 {\mu}m at high confidence (S/N>=6) in the nuclear region of the nearby spiral galaxy M83. Emission line maps of these high ionization lines show several compact structures. Specifically, the [Ne VI] emission is located at 140 pc from the optical nucleus and appears as a point source of size ~<18 pc (FWHM =<0.8"). We investigate the possible source of this extreme emission through comparison with photoionization models and ancillary data. We find that photoionization models of fast radiative shocks are able to reproduce the observed high excitation emission line fluxes only for the lowest preshock density available in the library, n =0.01 cm^-3. Additionally, tailored active galactic nuclei (AGN) photoionization models assuming a two-zone structure are compatible with the observed high ionization fluxes. Our simple AGN model shows that the emission at the location of the [Ne VI] source can be the result of a cloud being ionized by the radiation cone of an AGN. We stress, however, that to definitively confirm an AGN as the main source of the observed emission, more complex modeling accounting for different geometries is required. Previously known as a purely starburst system, these new findings of the nuclear region of M83 will require a reassessment of its nature and of objects similar to it, particularly now that we have access to the unparalleled infrared sensitivity and spatial resolution of the James Webb Space Telescope.

Discovered over fifty years ago, neutron stars exhibit a remarkable variety of behaviors depending on their age, magnetic field strength, rotational dynamics, emission mechanisms, and surrounding environments. This diversity in their observational manifestations has led astronomers to classify neutron stars into numerous categories, much like wandering through a zoo and admiring various species. This chapter focuses on isolated (i.e., non-accreting) neutron stars. We review extensive observational results of rotation-powered pulsars, magnetars, X-ray dim isolated neutron stars and central compact objects, highlighting the unique properties and recent discoveries of these neutron star classes. We briefly touch on theoretical models that have significantly advanced our understanding of these classes, emphasizing the valuable insights they provide into the underlying physics of neutron stars.

The spectral energy distribution (SED) of a galaxy represents the distribution of electromagnetic radiation emitted across all wavelengths, from radio waves to gamma rays. The galaxy SED is akin to its fingerprint, and serves as a fundamental tool in modern astrophysics. It enables researchers to determine crucial properties of galaxies, including their star-formation rates, stellar populations, dust content, and evolutionary state. By analyzing galactic SEDs, astronomers can reconstruct the physical processes occurring within galaxies and trace their evolutionary histories. This article explores our current understanding of the components that contribute to galactic SEDs, the observational techniques used to measure them, and their applications in understanding galaxy formation and evolution.

Over a hundred gravitational-wave (GW) detections and candidates have been reported from the first three observing runs of the Advanced LIGO-Virgo-KAGRA (LVK) detectors. Among these, the most intriguing events are binary black hole mergers that result in a 'lite' intermediate-mass black hole (IMBH) of ${\sim}10^2~\mathrm{M}_\odot$, such as GW170502 and GW190521. In this study, we investigate 11 GW candidates from LVK's Third Observing Run (April 2019-March 2020) that have a total detector-frame masses in the lite IMBH range. Using the Bayesian inference algorithm \texttt{RIFT}, we systematically analyze these candidates with three state-of-the-art waveform models that incorporate higher harmonics, which are crucial for resolving lite IMBHs in LVK data. For each candidate, we infer the pre-merger and post-merger black hole masses in the source frame, along with black hole spin projections across all three models. Under the assumption that these are binary black hole mergers, our analysis finds that 5 of them have a post-merger lite IMBH with masses ranging from $110\sim 350~\mathrm{M}_\odot$ with over 90\% confidence interval. Additionally, we note that one of their pre-merger black holes is within the pair-instability supernova mass gap ($60-120~\mathrm{M}_\odot$) with more than 90\% confidence interval, and additional two pre-merger black holes above the mass-gap. Furthermore, we report discrepancies among the three waveform models in their mass and spin inferences of lite IMBHs, with at least three GW candidates showing deviations beyond accepted statistical limits. While the astrophysical certainty of these candidates cannot be established, our study provides a foundation to probe the lite IMBH population that emerge within the low-frequency noise spectrum of LVK detectors.

3C 454.3 is a flat spectrum radio quasar (FSRQ) known for its high variability across the electromagnetic spectrum, showing structural and flux variability in its pc-scale jet, and correlated variability among frequency bands. This study aims to identify the structure, dynamics, and radiative processes common to the innermost regions of the blazar 3C 454.3. We investigate whether any jet component can be associated with $\gamma-$ray emission and variability. We analyze the relationship between the variable $\gamma-$ray emission and pc-scale jet properties in 3C 454.3 by combining $\gamma-$ray data spanning twelve years with contemporaneous VLBA multi-epoch images at 15 and 43 GHz. Spearman rank correlation tests are conducted to determine if the flux variability of any jet component is associated with $\gamma-$ray variability. Core emission at 43 and 15 GHz strongly correlates with $\gamma-$ray emission. The 43 GHz core (Q0) contributes around 37$\%$ of the observed $\gamma-$ray variability, while the 15 GHz core (K0) accounts for 30$\%$. A quasi-stationary component at 43 GHz, at a projected distance of 4.6 pc, correlates with the $\gamma-$ray flux, accounting for 20$\%$ of its emission between 2016 and 2021. We found a mobile component (Q3 between 2010.18 and 2011.16) at 43 GHz with a projected distance between 0.8 and 2.3 pc and apparent velocity of $\beta_{app} = 9.9 \pm 1.1$ c, accounting for approximately 28% of the $\gamma-$ray emission. The observed simultaneous variability in emission regions beyond the central parsec strongly suggests synchrotron self-Compton (SSC) as the primary mechanism for $\gamma-$ray production in these regions. Our findings demonstrate the existence of multiple $\gamma-$ray emission regions within the blazar jet but also suggest that some of these regions are non-stationary over time.

We combined Gaia DR3 and TESS photometric light curves to estimate the internal physical properties of 2,497 gravity-mode pulsators. We relied on asteroseismic properties of Kepler $\gamma\,$Dor and SPB stars to derive the near-core rotation frequency, $f_{\rm rot}$, of the Gaia-discovered pulsators from their dominant prograde dipole gravito-inertial pulsation mode. We offer a recipe based on linear regression to deduce $f_{\rm rot}$ from the dominant gravito-inertial mode frequency. It is applicable to prograde dipole modes with an amplitude above 4mmag and occurring in the sub-inertial regime. By applying it to the 2,497 pulsators, we have increased the sample of intermediate-mass dwarfs with such an asteroseismic observable by a factor of 4. We used the estimate of $f_{\rm rot}$ to deduce spin parameters between 2 and 6, while the sample's near-core rotation rates range from 0.7% to 25% of the critical Keplerian rate. We used $f_{\rm rot}$, along with the Gaia effective temperature and luminosity to deduce the (convective core) mass, radius, and evolutionary stage from grid modelling based on rotating stellar models. We derived a decline of $f_{\rm rot}$ with a factor of 2 during the main-sequence evolution for this population of field stars, which covers a mass range from 1.3M$_\odot$ to 7M$_\odot$. We found observational evidence for an increase in the radial order of excited gravity modes as the stars evolve. For 307 pulsators, we derived an upper limit of the radial differential rotation between the convective core boundary and the surface from Gaia's vbroad measurement and found values up to 5.4. Our recipe for the near-core rotation frequency from the dominant gravito-inertial mode detected in the independent Gaia and TESS light curves is easy to use, facilitates applications to large samples, and allows to map their angular momentum and evolutionary stage in the Milky Way.

Gaia DR3, released in June 2022, included low-resolution XP spectra that have been used for the classification of various types of emission-line objects through machine-learning techniques. The Gaia Extended Stellar Parametrizer for Emission-Line Stars (ESP-ELS) algorithm identified 565 sources as potential Wolf-Rayet (WR) stars. Over half of them were already known as WR stars in the Milky Way and Magellanic Clouds. This study aimed to utilize Gaia DR3 data to identify new Galactic WR stars. We extracted all sources classified as WC or WN type stars by the ESP-ELS algorithm from the Gaia catalog. By applying judicious 2MASS color selection criteria, leveraging Gaia H$\alpha$ measurements, and filtering out objects already cataloged in various databases, we selected 37 bright candidates ($G \leq $ 16 mag) and 22 faint candidates ($G > $ 16 mag). Spectroscopic follow-up observations of these candidates were conducted using the 2SPOT facilities in Chile and France, as well as 1-m C2PU's Epsilon telescope at the Calern Observatory. This paper focuses on the brighter sample. Among the 37 targets, we confirmed 17 and 16 new Galactic WC and WN type WR stars, respectively. Three of them were recently reported as new WR stars in an independent study. The Gaia mission provides a valuable resource for uncovering WR stars missed in earlier surveys. While this work concentrated on a relatively small starting sample provided by the ESP-ELS algorithm, our findings highlight the potential for refining selection criteria to identify additional candidates not included in the outputs of the algorithm. Furthermore, the observation program underscores the utility of small telescopes in acquiring initial spectral data for sources with magnitudes up to $G \sim 16$ mag.

The light curves of old G-dwarfs obtained in the visible and near-infrared wavelength ranges are highly irregular. This significantly complicates the detectability of the rotation periods of stars similar to the Sun in large photometric surveys, such as Kepler and TESS. In this study, we show that light curves collected in the ultraviolet wavelength range are much more suitable for measuring rotation periods. Motivated by the observation that the Sun's rotational period is clearly discernible in the UV part of the spectrum, we study the wavelength dependence of the rotational period detectability. We employ the Spectral and Total Solar Irradiance Reconstructions model, SATIRE-S, to characterize the detectability of the solar rotation period across various wavelengths using the autocorrelation technique. We find that at wavelengths above 400 nm, the probability of detecting the rotation period of the Sun observed at a random phase of its activity cycle is approximately 20\%. The probability increases to 80\% at wavelengths shorter than 400 nm. These findings underscore the importance of ultraviolet stellar photometry.

With the growing number of gravitational wave detections, achieving a competitive measurement of $H_0$ with dark sirens is becoming increasingly feasible. The expansion of the Ligo-Virgo-KAGRA Collaboration into a four detector network will reduce both the localisation area and the luminosity distance uncertainty associated with each gravitational wave event. It is therefore essential to identify and mitigate other major sources of error that could increase the uncertainty in $H_0$. In this work, we explore three scenarios relevant to the dark siren method in future observing runs. First, we demonstrate that there is a precision gain offered by a catalogue of spectroscopic-like redshifts compared to photometric-like redshifts, with the greatest improvements observed in smaller localisation areas. Second, we show that redshift outliers (as occur in realistic photometric redshift catalogues), do not introduce bias into the measurement of $H_0$. Finally, we find that uniformly sub-sampling spectroscopic-like redshift catalogues increases the uncertainty in $H_0$ as the completeness fraction is decreased; at a completeness of 50% the benefit of spectroscopic redshift precision is outweighed by the degradation from incompleteness. In all three scenarios, we obtain unbiased estimates of $H_0$. We conclude that a competitive measurement of $H_0$ using the dark siren method will require a hybrid catalogue of both photometric and spectroscopic redshifts, at least until highly complete spectroscopic catalogues become available.

Context. Omega Cen is the largest known globular cluster in the Milky Way. It is also quite a complex object with a large metallicity spread and multiple stellar populations. Despite a number of studies over the past several decades, the series of events that led to the formation of this cluster is still poorly understood. One of its peculiarities is the presence of a metal-rich population which does not show the phenomenon of light-element anti-correlations (C-N, Na-O, Mg-Al), a trait that is considered as characteristic of Galactic Globular Clusters, and is in fact present among more metal-poor Omega Cen stars, leading to speculations that such anomalous population was accreted by the cluster. In this paper, we aim at investigating the kinematics of Omega Cen populations to gain insight on the formation scenario of the cluster. Using the newly released Gaia FPR and DR3 catalogue, we conducted a detailed kinematical analysis of cluster members within Omega Cen. The cluster members were divided into four metallicity populations and their mean proper motion in radial and tangential components were compared with each other. We also performed Gaussian Mixture Model fitting on the metallicity distribution to estimate the number of populations within our sample and an independent analysis of the HST catalogue as confirmation. The mean proper motions (mu_r and mu_t) of the metallicity populations do not show any significant differences. It is also not dependent on the approach chosen to determine the number of metallicity populations. We do find clear signature of rotation in all of the populations (including the metal-rich) with similar velocities.

An important step in understanding the formation and evolution of the Nuclear Star Cluster (NSC) is to investigate its chemistry and chemical evolution. Additionally, exploring the relationship of the NSC to the other structures in the Galactic Center and the Milky Way disks is of great interest. Extreme optical extinction has previously prevented optical studies, but near-IR high-resolution spectroscopy is now possible. Here, we present a detailed chemical abundance analysis of 19 elements - more than four times as many as previously published - for 9 stars in the NSC of the Milky Way, observed with the IGRINS spectrometer on the Gemini South telescope. This study provides new, crucial observational evidence to shed light on the origin of the NSC. We demonstrate that it is possible to probe a variety of nucleosynthetic channels, reflecting different chemical evolution timescales. Our findings reveal that the NSC trends for the elements F, Mg, Al, Si, S, K, Ca, Ti, Cr, Mn, Co, Ni, Cu, and Zn, as well as the s-process elements Ba, Ce, Nd, and Yb, generally follow the inner bulge trends within uncertainties. This suggests a likely shared evolutionary history and our results indicate that the NSC population is consistent with the chemical sequence observed in the inner Galaxy (the inner-disk sequence). However, we identify a significant and unexplained difference in the form of higher Na abundances in the NSC compared to the inner-bulge. This is also observed in few Galactic globular clusters, and may suggest a common enrichment process at work in all these systems.

We present Atacama Large Millimeter/submillimeter Array Band 3 observations of N$_2$H$^+$ (1-0) and CH$_3$CN (5-4), as well as Band 7 observations of the H$_2$CO molecular line emissions from the protostellar system GGD 27-MM2(E). Through position-velocity diagrams along and across the outflow axis, we study the kinematics and structure of the outflow. We also fit extracted spectra of the CH$_3$CN emission to obtain the physical conditions of the gas. We use the results to discuss the impact of the outflow on its surroundings. We find that N$_2$H$^+$ emission traces a dense molecular cloud surrounding GGD 27-MM2(E). We estimate that the mass of this cloud is $\sim$13.3-26.5 M$_\odot$. The molecular cloud contains an internal cavity aligned with the H$_2$CO-traced molecular outflow. The outflow, also traced by $\mathrm{CH_3 CN}$, shows evidence of a collision with a molecular core (MC), as indicated by the distinctive increases in the distinct physical properties of the gas such as excitation temperature, column density, line width, and velocity. This collision results in an X-shape structure in the northern part of the outflow around the position of the MC, which produces spray-shocked material downstream in the north of MC as observed in position-velocity diagrams both along and across of the outflow axis. The outflow has a mass of 1.7-2.1 M$_\odot$, a momentum of 7.8-10.1 M$_\odot$ km s$^{-1}$, a kinetic energy of 5.0-6.6$\times 10^{44}$ erg, and a mass loss rate of 4.9--6.0$\times10^{-4}$ M$_\odot$ yr$^{-1}$. The molecular outflow from GGD 27-MM2(E) significantly perturbs and erodes its parent cloud, compressing the gas of sources such as MC and ALMA 12. The feedback from this powerful protostellar outflow contributes to maintain the turbulence in the surrounding area.

Radio interferometers must grapple with apparent fields of view that distort the true radio sky. These so-called 'A-term' distortions may be direction-, time- and baseline-dependent, and include effects like the primary beam and the ionosphere. Traditionally, properly handling these effects has been computationally expensive and, instead, less accurate, ad-hoc methods have been employed. Image domain gridding (IDG; van der Tol et al., 2018) is a recently developed algorithm that promises to account for these A-terms both accurately and efficiently. Here we describe a new implementation of IDG known as the Parallel Interferometric GPU Imager (Pigi). Pigi is capable of imaging at rates of almost half a billion visibilities per second on modest hardware, making it well suited for the projected data rates of the Square Kilometre Array, and is compatible with both NVIDIA and AMD GPU hardware. Its accuracy is principally limited only by the degree to which A-terms are spatially sampled. Using simulated data from the Murchison Widefield Array, we demonstrate the extraordinary effectiveness of Pigi in correcting for extreme ionospheric effects, and point to future work that would enable these results on real-world data.

We constrain two vacuum decay models ($\Lambda(t)$CDM, proposed by the authors of~\cite{Brito:2024bhh}) utilizing the baryon acoustic oscillations (BAO) data released by the Dark Energy Spectroscopic Instrument (DESI), distance prior from the cosmic microwave background (CMB) observed by the Planck satellite, Hubble rate data obtained via the cosmic chronometers (CC) method and type Ia supernova (SNIa) data. The interaction terms between dark matter and dark energy are defined as $Q=3\varepsilon H\rho_{\Lambda}$ for model I and $3\varepsilon aH\rho_{\Lambda}$ for model II. We find that the decay parameter is constrained to be $\varepsilon=0.0094^{+0.0037}_{-0.0033}$ for model I and $\varepsilon=0.0119\pm{0.0045}$ for model II, respectively, indicating a potential interaction between dark matter and dark energy at the $2\sigma$ confidence level. The current Hubble parameter values are estimated to be $H_{0}=70.30\pm{0.67}$ for model I and $H_{0}=70.28\pm{0.64}$ for model II. These values of $H_0$ fall between those derived from the Planck and SH0ES data, suggesting that these two vacuum decay models could provide a potential solution to alleviate the Hubble tension problem.

The objective of this study is to develop a fully compressible magnetohydrodynamic solver for fast simulations of the global dynamo of the Sun using unstructured grids and GPUs. Accurate modeling of the Sun's convective layers is vital to predicting the Sun's behavior, including the solar dynamo and sunspot cycles. Currently, there are many efficient codes capable of conducting these large simulations; however, many assume an anealastic density distribution. The anelastic assumption is capable of producing accurate results for low mach numbers; however, it fails in regions with a higher mach number and a fully compressible flow must be considered. To avoid these issues, Wang et al. [1] created a Compressible High-ORder Unstructured Spectral difference (CHORUS) code for simulating fluid dynamics inside stars and planets. CHORUS++ augmented the CHORUS code to adopt a higher degree of polynomials by using cubed-sphere meshing and transfinite mapping to perform simulations on unstructured grids [2]. Recently, CHORUS++ was further developed for parallel magnetohydrodynamic (MHD) solutions on GPUs at Clarkson University. In this study the solar benchmark problems presented by Chen et al. [2] are extended to unsteady solar dynamo problems, with two different density scale heights. The CHORUS-MHD code is further accelerated by multiple GPUs and used to successfully solve these solar dynamo benchmark problems. [1] Wang, J., Liang, C., and Miesch, M. S., "A Compressible High-Order Unstructured Spectral Difference Code for Stratified Convection in Rotating Spherical Shells," Journal of Computational Physics, Vol. 290, 2015, pp. 90-111. [2] Chen, K., Liang, C., and Wan, M., "Arbitrarily high-order accurate simulations of compressible rotationally constrained convection using a transfinite mapping on cubed-sphere grids," Physics of Fluids, Vol. 35, 2023, p. 086120.

Lightsails are a highly promising spacecraft concept that has attracted interest in recent years due to its potential to travel at near-relativistic speeds. Such speeds, which current conventional crafts cannot reach, offer tantalizing opportunities to probe nearby stellar systems within a human lifetime. Recent advancements in photonics and metamaterials have created novel solutions to challenges in propulsion and stability facing lightsail missions. This review introduces the physical principles underpinning lightsail spacecrafts and discusses how photonics coupled with inverse design substantially enhance lightsail performance compared to plain reflectors. These developments pave the way through a previously inaccessible frontier of space exploration.

The aromatic molecule benzene is considered the essential building block for larger polycyclic aromatic hydrocarbons (PAHs) in space. Despite benzene's importance in the formation of PAHs, the formation mechanisms of interstellar benzene are not well understood. A single ion-molecule reaction sequence is considered when modeling the formation of benzene in the interstellar medium, beginning with the protonation of acetylene. Although this process has been used to model the initial steps for formation of PAHs, it has not been experimentally measured. To explore this reaction mechanism, we have carried out the first experimental study of sequential ion-molecule reactions beginning with protonation of acetylene at single-collision conditions. Surprisingly, we find that the reaction sequence does not result in benzene and instead terminates at C6H5+, which is unreactive toward both acetylene and hydrogen. This result disproves the previously proposed mechanism for interstellar benzene formation, critically altering our understanding of interstellar PAH formation.

Recurrent and propagating intensity perturbations are frequently observed in extreme ultraviolet (EUV) channels along coronal fan loops above sunspots, and these perturbations are suggested to be slow magnetoacoustic waves. Numerous studies have been conducted to investigate their propagation speeds, damping, and excitation sources; however, there have been limited observational analyses on whether these waves are dispersive despite some theoretical studies. In this study, we apply cross-correlation analysis in the Fourier domain on slow magnetoacoustic waves using three different datasets: EUV intensity observed by SDO/AIA, differential emission measure (DEM) temperature maps, and Doppler velocities from Hinode/EIS spectrometer observations. The apparent phase velocities of the waves, which are the plane-of-sky component of the waves' phase velocities, are derived as functions of frequency for all the three datasets. It is found that the phase velocities show clear frequency dependency, with a general trend of increase with frequency, ranging from approximately 30 km/s around 3 mHz to about 80 km/s around 10 mHz. The frequency dependency of the phase velocities demonstrates that the slow magnetoacoustic waves in the coronal loops are dispersive. The dispersiveness of these waves can provide a useful tool for the diagnosis of physical conditions inside the coronal loops along which these waves travel.

This study presents an energy-dependent analysis of seasonal variations in the atmospheric muon neutrino spectrum, using 11.3 years of data from the IceCube Neutrino Observatory. By leveraging a novel spectral unfolding method, we explore the energy range from 125 GeV to 10 TeV for zenith angles between 90° to 110°, corresponding to the Antarctic atmosphere. Our findings reveal that the seasonal variation amplitude decreases with energy reaching ($-4.6 \pm 1.1$)\% during Austral winter and increases ($+3.9 \pm 1.2$)\% during Austral summer relative to the annual average at 10TeV. While the unfolded flux exceeds the model predictions by up to 30\%, the differential measurement of seasonal variations remains unaffected. The measured seasonal variations of the muon neutrino spectrum are consistent with theoretical predictions using the MCEq code and the NRLMSISE-00 atmospheric model.

(Reduced) Using the NewHorizon simulation, we have studied ten gas-versus-star counter-rotating galaxies in field environments with a stellar mass of M*~[1-5]x10^10 Msun. For all of them, the retrograde accretion of gas either from gas stripping from a nearby companion or from the circumgalactic medium is the starting point of the formation process. Then follows the co-existence of two distinct disks of gas rotating in opposite directions, the pre-existing one in the inner parts and the accreted gas in the outer parts of the galaxy. The latter progressively replaces the former leading to the final gas-star kinetic misalignment configuration. During the process, the star formation is first enhanced and then progressively decreases. We roughly estimate that a higher fraction of the pre-existing gas is converted into stars rather than being expelled. We also found that the black hole activity (BH) tends to be enhanced during the removal of the pre-existing gas. Furthermore, our analysis suggests that the formation of a counter-rotating gas component is always accompanied with the formation of counter-rotating stellar disks. These stellar disks can have diverse properties but host in general a younger and more metal rich population of stars with respect to the main disc, depending on the star formation history and BH activity. The central part of counter-rotating disks tend also to be characterized by a younger population, an enhanced star formation rate and a higher metallicity than their outer parts. The high metallicity comes the progressive metal enrichment of the accreted gas by mixing with the pre-existing gas and by supernovae activity as it sinks toward the center of the galaxy. In case of major mergers, a large amount of accreted stars from the companion would be distributed at large distances from the remnant center due to conservation of the initial orbital angular momentum.

Millisecond pulsars, representing the older neutron star population, are believed to have undergone a prolonged period of dark matter accumulation, resulting in a higher dark matter content. Their extreme rotation makes them unique laboratories for studying rapidly rotating neutron stars admixed with dark matter. In this work, we model uniformly rotating neutron stars with a dark matter component that rotates independently from the baryon matter, allowing for the investigation of both co-rotating and counter-rotating scenarios. We examine the impact of dark matter rotation on the macroscopic properties of neutron stars, including the mass-radius relation, the mass-shedding Keplerian limit, and moments of inertia, for various dark matter particle masses and total fractions, considering both core and halo distributions. Our findings provide a more comprehensive understanding of how dark matter influences the equilibrium properties of rotating neutron stars, offering new insights into the astrophysical implications of self-interacting dark matter.

HERMES Pathfinder is an in-orbit demonstration consisting of a constellation of six 3U cubesats hosting simple but innovative X-ray/gamma-ray detectors for the monitoring of cosmic high-energy transients. HERMES-PF, funded by ASI and by the EC Horizon 2020 grant, is scheduled for launch in Q1 2025. An identical X-ray/gamma-ray detector is hosted by the Australian 6U cubesat SpIRIT, launched on December 1st 2023. The main objective of HERMES-PF/SpIRIT is to demonstrate that high energy cosmic transients can be detected efficiently by miniatured hardware and localized using triangulation techniques. The HERMES-PF X-ray/gamma-ray detector is made by 60 GAGG:Ce scintillator crystals and 12 2x5 silicon drift detector (SDD) mosaics, used to detect both the cosmic X-rays directly and the optical photons produced by gamma-ray interactions with the scintillator crystals. This design provides a unique broad band spectral coverage from a few keV to a few MeV. Furthermore, the use of fast GAGG:Ce crystals and small SDD cells allows us to reach an exquisite time resolution better than a microsecond. We present a progress report on the missions focusing the discussion on the scientific innovation of the project and on the main lessons learned during the project development including: the importance and the challenges of using distributed architectures to achieve ambitious scientific objectives; the importance of developing critical technologies under science agreements for the realization of high-performing but low-cost payloads; best use of COTS technologies in scientific missions. We finally discuss the prospects of applying these concepts for the creation of an all-sky, all-time monitor to search for the high-energy counterparts of gravitational wave events that Advanced LIGO/Virgo/Kagra will find at the end of this decade and the Einstein Telescope during the 2030s.

Solar active regions (ARs) are the places hosting the majority of solar eruptions. Studying the evolution and morphological features of ARs is not only of great significance to the understanding of the physical mechanisms of solar eruptions, but also beneficial for the hazardous space weather forecast. An automated DBSCAN-based Solar Active Regions Detection (DSARD) method for solar ARs observed in magnetograms is developed in this work, which is based on an unsupervised machine learning algorithm called Density-Based Spatial Clustering of Applications with Noise (DBSCAN). The method is then employed to identify ARs on the magnetograms observed by the Helioseismic and Magnetic Imager (HMI) onboard Solar Dynamics Observatory (SDO) during solar cycle 24 and the rising phase of solar cycle 25. The distributions of the number, area, magnetic flux, and the tilt angle of bipolar of ARs in latitudes and time intervals during solar cycle 24, as well as the butterfly diagram and drift velocities are obtained. Most of these statistical results based on our method are in agreement with previous studies, which also guarantees the validity of the method. In particular, the dipole tilt angles in ARs in solar cycle 24 and the rising phase of solar cycle 25 are analyzed which reveal that 13% and 16% of ARs, respectively, violate Hale's law.

In previous works, we have shown that stars in the Orion and the Lagoon Nebula Clusters, and simulations of collapsing clouds, exhibit constant velocity dispersion as a function of mass, a result described by Lynden-Bell 50 years ago as an effect of a violent relaxation mechanism. In contrast, numerical simulations of turbulent clouds show that newborn massive stars experience stronger dynamical heating than low-mass stars. We analyzed turbulent numerical simulations and found that this effect arises from the fact that, in clouds that are globally turbulence-supported against collapse, massive stars are formed within more massive and denser clumps and in more crowded environments compared to low-mass stars. This allows them to accrete more mass and interact with other stars simultaneously. As they become more massive, their orbits tighten, increasing their velocity dispersion. In contrast, low-mass stars are formed in the periphery of such cores, more separated, and at lower densities. Thus, their velocity dispersion remains lower because they do not accrete as vigorously as massive stars and tend to be more isolated. We call this mechanism "accretion-induced orbital tightening." Our results and previous findings about violent relaxation provide a key observational diagnostic of how to distinguish the dynamic state of star-forming molecular clouds through the kinematics of their newborn stars.

We report observations of stripe-like features in Enceladus' plumes captured simultaneously by Cassini's VIMS-IR and ISS NAC instruments during flyby E17, with similar patterns seen in VIMS-IR data from flyby E13 and E19. These parallel stripes, inclined at approximately 16$^{\circ}$ to the ecliptic and 43$^{\circ}$ to Saturn's ring plane, appear continuous across images when projected in the J2000 frame. A bright stripe, most visible at wavelengths around 5 $\mu$m, acts as the zeroth-order diffraction peak of a reflection grating with an estimated groove spacing of 0.12$-$2.60 mm, while adjacent stripes are attributed to higher-order diffraction peaks. We suggest that this light-dispersing phenomenon originates from an inclined periodic structure within Saturn's E ring. This structure, constrained between Saturn's G-ring and Rhea's orbit, likely consists of fresh ice particles supplied by Enceladus' plumes.

In the era of the James Webb Space Telescope (JWST), the dramatic improvement in the spectra of exoplanetary atmospheres demands a corresponding leap forward in our ability to analyze them: atmospheric retrievals need to be performed on thousands of spectra, applying to each large ensembles of models (that explore atmospheric chemistry, thermal profiles and cloud models) to identify the best one(s). In this limit, traditional Bayesian inference methods such as nested sampling become prohibitively expensive. We introduce FASTER (Fast Amortized Simulation-based Transiting Exoplanet Retrieval), a neural-network based method for performing atmospheric retrieval and Bayesian model comparison at a fraction of the computational cost of classical techniques. We demonstrate that the marginal posterior distributions of all parameters within a model as well as the posterior probabilities of the models we consider match those computed using nested sampling both on mock spectra, and for the real NIRSpec PRISM spectrum of WASP-39b. The true power of the FASTER framework comes from its amortized nature, which allows the trained networks to perform practically instantaneous Bayesian inference and model comparison over ensembles of spectra -- real or simulated -- at minimal additional computational cost. This offers valuable insight into the expected results of model comparison (e.g., distinguishing cloudy from cloud-free and isothermal from non-isothermal models), as well as their dependence on the underlying parameters, which is computationally unfeasible with nested sampling. This approach will constitute as large a leap in spectral analysis as the original retrieval methods based on Markov Chain Monte Carlo have proven to be.

Transitional millisecond pulsars (tMSPs) represent a crucial link between the rotation-powered and accretion-powered states of binary pulsars. During their active X-ray state, tMSPs are the only low-mass X-ray binary systems detected up to GeV energies by the Fermi Large Area Telescope (LAT). CXOU J110926.4-650224 is a newly discovered tMSP candidate in an active X-ray state, potentially spatially compatible with a faint gamma-ray source listed in the latest Fermi-LAT point-source catalogue as 4FGL J1110.3-6501. Confirming the association between CXOU J110926.4-650224 and the Fermi source is a key step toward validating its classification as a tMSP. In this study, we analyse Fermi-LAT data collected from August 2008 to June 2023 to achieve a more accurate localisation of the gamma-ray source, characterise its spectral properties, and investigate potential time variability. By thoroughly reconstructing the gamma-ray background around the source using a weighted likelihood model, we obtain a new localisation that aligns with the position of the X-ray source at the 95% confidence level, with a Test Statistic value of $\sim 42$. This establishes a spatial association between the gamma-ray source and CXOU J110926.4-650224. The gamma-ray emission is adequately described by a power-law model with a photon index of $\Gamma = 2.5 \pm 0.1$ and a corresponding flux of $(3.7\pm0.9) \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ in the 0.1-300 GeV range.

We present a novel approach for identifying members of open star clusters using Gaia DR3 data by combining Minimum Spanning Tree (MST) and Gaussian Mixture Model (GMM) techniques. Our method employs a three-step process: initial filtering based on astrometric parameters, MST analysis for spatial distribution filtering, and GMM for final membership probability determination. We tested this methodology on 12+1 open clusters of varying ages, distances, and richness. The method demonstrates superior performance in distinguishing cluster members from field stars, particularly in regions with overlapping populations, as evidenced by its application to clusters like NGC 7790. By effectively reducing the number of probable field stars through MST analysis before applying GMM, our approach enhances both computational efficiency and membership determination accuracy. The results show strong agreement with previous studies while offering improved precision in member identification. This method provides a robust framework for analyzing the extensive datasets provided by Gaia DR3, addressing the challenges of processing large-scale astronomical data while maintaining high accuracy in cluster membership determination.

The formation of the first stars and the subsequent population of X-ray binaries represents a fundamental transition in the state of the Universe as it evolves from near homogeneity to being abundant in collapsed structures such as galaxies. Due to a lack of direct observations, the properties of these stars remain highly uncertain. By considering the impact of the first stars and their remnant X-ray binaries on the cosmological 21-cm signal, we demonstrate that upcoming observations have the potential to significantly improve our understanding of these objects. We find a 25 mK sensitivity measurement of the 21-cm global signal by a wide-beam radiometer, such as REACH, or 3,000 hours of foreground avoidance observations of the 21-cm power spectrum by SKA-Low, could provide three-sigma constraints on the mass distribution of the first stars. Such measurements will fill a critical gap in our understanding of the early Universe and aid in interpreting high-redshift galaxy observations.

Pulsating variables are prevalent in the classical {\delta} Scuti instability strip of intermediate-age open star clusters. The cluster membership of these stars facilitates a comparative analysis of their evolution in analogous environments. In this study, we integrate ground-based observations, TESS Full Frame Images (FFIs), and Gaia DR3 data to investigate variable stars in the intermediate-age open star cluster NGC 2126. We performed ground-based time-series observations of NGC 2126 to identify variable stars within its vicinity. Next, we determined the membership of these stars using parallax and the proper motions from Gaia DR3 archive. Then, we searched the TESS Full Frame Images (FFIs) for counterparts to the variables identified above and performed their frequency analysis and classification. Finally, we modeled the light curves (LCs) of detected eclipsing binaries (EBs), including V551 Aur, which has a pulsating component. We found 25 members and 85 field variable stars. In TESS FFIs, we found LCs for 11 known variables and a new rotational variable. We determined that the pulsating EB V551 Aur is a member of the cluster. The low- and medium-resolution spectra revealed the line profile variation and the basic parameters for the star, respectively. Simultaneous modeling of the eclipses and the embedded pulsations resulted in improved orbital parameters for the binary system. We also report the determination of orbital parameters for the previously uncharacterized EB system UCAC4 700-043174.

The review aims to give an overview of atmospheric escape processes from exoplanets. I briefly discuss the physics of various escape processes responsible for atmospheric escape across different types of exoplanets. Transmission spectroscopy is one of the major workhorses to observe the escaping atmosphere from exoplanets. I discuss recent observations that established the fact that atmospheric escape is very common in exoplanets, especially during the early phase of their evolution when stellar high-energy radiation (X-ray and extreme ultraviolet, hence XUV) is strong. There are many theoretical efforts/models to understand atmospheric escape processes. Stellar radiation is one of the major drivers of atmospheric escape, but other stellar environments (e.g., stellar flares, stellar winds, stellar coronal mass ejections, and stellar magnetic field) also have control over how the escape process will be affected for a given property of exoplanet, as the planetary properties (e.g., gravity, thermal energy, magnetic field) plays an important role for atmospheric escape. I discuss all governing factors for the atmospheric escape process and corresponding theoretical models in detail. I also discuss how atmospheric escape plays a crucial role in the overall atmospheric evolution of exoplanets and can lead us to understand some features in recently observed exoplanet demographics.

Late-type stars are known to host numerous exoplanets, and their photometric variability, primarily caused by rotational modulation, provides a unique opportunity to study starspots. As exoplanets transit in front of their host stars, they may occult darker, spotted regions on the stellar surfaces. The monitoring of starspots from planetary transits, known as transit mapping, offers a possibility to detect small dark regions on magnetically active, late-type stars. These spots may be so small that they would be undetectable to other methods used to reconstruct stellar magnetic activity. We describe a Bayesian analysis framework on the transit light curves of planets orbiting K- and M-type main-sequence stars in search for spot occultation event candidates. We present a systematic analysis of high-precision, high-cadence light curves from Kepler and TESS to detect and characterise starspots during exoplanetary transits. According to our tests, the set of criteria applied in the analysis is robust and not prone to false positives. Our sample comprises K- and M-dwarfs hosting transiting exoplanets observed by the Kepler or TESS space telescopes at a high cadence, totalling 99 planets meeting our selection criteria. After analysing 3273 transit light curves from 99 planets, we find 105 candidates for starspot occultation events by six planets. We report new spot occultation candidates for the K-dwarfs HD 189733 and TOI-1268. The identified dark regions have a lower limit for radii between 1.6 degrees and 29.5 degrees and contrasts up to 0.69. We estimate a spot detection frequency of 3.7% and 4.2% for K- and M-dwarfs by TESS, and 37.5% for K-dwarfs by Kepler.

Parker Solar Probe (PSP) observations revealed that most of the solar wind acceleration occurs very close to the Sun. This acceleration is partly due to the global electric potential originating from the mass disparity between electrons and protons, coupled with the constraints of charge quasi-neutrality and zero-current conditions in the solar wind plasma. However, the exact mechanism that accounts for the remaining acceleration is still not identified. We aim to provide a framework incorporating the electric field-driven component of the acceleration while also introducing an additional acceleration mechanism via velocity-space diffusion of the particles. This will help us determine the extent of extra acceleration, beyond the electric field-driven component, required to fully reproduce the acceleration of the solar wind in theoretical models. We modified an existing kinetic exospheric model to account for the unexplained solar wind acceleration by including velocity-space diffusion. We compared the electric field derived from the sunward deficit of velocity distribution functions observed by PSP between 13.3 and 50 solar radii ($R_s$) with the electric field found self-consistently by the kinetic exospheric model. The effect of velocity-space diffusion is found to reduce the temperature anisotropy and impact the solar wind acceleration while maintaining the electric potential unchanged. Our approach allows for adjusting the diffusion to effectively reduce or increase the solar wind acceleration, providing a way to account for the unexplained acceleration within the kinetic exospheric framework. It is found that, even without diffusion, the model can reproduce the anticorrelation between the electric potential and the solar wind terminal velocity found by PSP. This suggests that the electric potential might still be of major importance in explaining the solar wind acceleration.

We identified a sample of 88 bona fide Population II Cepheids (henceforth referred to as Cepheids) and 44 candidates in Galactic globular clusters (GCs). Seventy-eight of the Cepheids in the sample align within $2\,\sigma$ of the period-luminosity relation for Milky Way Type II Cepheids (T2CEPs). Nine align with the period-luminosity relation for fundamental-mode anomalous Cepheids (ACEPs), and only one (BL Boötis) follows the relation for first-overtone ACEPs, as determined from observations of ACEPs in the Large Magellanic Cloud. For sources in common between our catalog and the OGLE catalog, the classification agrees in 94% of cases. In comparison, for sources shared between the Gaia Specific Object Study (SOS) and OGLE, the agreement is 74%. In the dense environments of GCs, our analysis shows that the completeness of the Gaia catalogs for Cepheids is 64% for the SOS and 74% for the classifier of variable stars. We determined the red and blue edges of the instability strip for T2CEPs using linear MESA-RSP models. We find that the best-fit models, with $M = 0.6\,{\rm M}_\odot$ and Z = 0.0003, are able to fit $90%$ of the stars in our sample. This percentage is the same for helium abundances $Y = 0.220$ and $0.245$. Higher values of $Y$ lower this percentage, and the same effect is observed with lower values of $Z$. In the future, combining the sample of T2CEPs with the precise parallaxes obtained from GCs will strengthen the geometric calibration of a distance ladder based on Population II stars. This will be useful for determining distances within the Milky Way and for cross-checking distances to Local Group galaxies determined through other methods.

We present the results of deep Hubble Space Telescope photometry of the dwarf galaxy DDO 68-C, proposed as possibly associated with the isolated peculiar dwarf DDO 68. The new data resolve for the first time the stars of DDO 68-C down to well below the tip of the Red Giant Branch (RGB), revealing a low mass (M$_{*}$ $\sim$ 1.5 $\times$ 10$^7$ M$_{\odot}$) star forming galaxy with a backbone of old stars. By means of a fully homogeneous analysis and using the RGB tip as a standard candle, we find that DDO 68 and DDO 68-C lie at the same distance from us, within the uncertainties (D = 12.6 $\pm$ 0.3 Mpc and D = 12.7 $\pm$ 0.4 Mpc, respectively), thus confirming that the two dwarfs are physically associated. While paired dwarf galaxies with mutual projected distance similar to DDO 68 and DDO 68-C are not exceptional in the Lynx-Cancer Void where they live, DDO 68 remains a unicum as, in addition to the newly confirmed companion, it records the evidence of at least two other satellites.

This study proposes a 3D-printed, quad-ridged, dual-pol, flared horn antenna feeder to meet the demands of a modern radio-telescope. The research work presented in this paper involves both the theoretical and the customized / adapted mechanical design of the proposed horn antenna, as well as a description of the 3D-printing and coating process. A novel slotted waveguide approach is followed for perfectly homogeneous coating. Additional mechanical adaptions (i.e. optimized feed section, radome) are employed to provide realistic simulation results. The proposed antenna offers a dual polarization capability and an attractive radiation beam to properly illuminate a 6m mesh dish reflector. A satisfying impedance matching from 1-3 GHz is achieved in both polarizations while maintaining an excellent isolation between the two N-type ports. Reflection coefficient and radiation pattern simulation results are presented, derived from the theoretical and the adapted mechanical design, exhibiting a close agreement.

Nest200047 is a clear example of multiple radio bubbles from an Active Galactic Nucleus (AGN) in a galaxy group, featuring non-thermal filaments likely shaped by buoyancy, gas motions, and stabilized by magnetic fields. This study presents high-quality data obtained from uGMRT, MeerKAT, and VLA, alongside existing LOFAR data, to analyze the system's morphology and spectrum over a broad frequency range (53-1518 MHz). Our findings reveal new filamentary emission in the inner 60 kpc, surrounding and extending from the inner bubbles and jets, suggesting complex dynamical evolution of the non-thermal plasma in the group core. The filaments have widths of a few kpc and lengths from tens to hundreds of kpc, with a steep and curved radio spectrum ($\rm \alpha=1\sim2$). They exhibit a constant spectral index profile along their length, implying particles are either (re-)accelerated together or move at super-Alfvenic speeds. Spectral aging analysis yields jet active times between 50 and 100 Myr with short inactive phases, suggesting continuous energy injection typical of AGN feedback in galaxy groups. This study highlights the potential of combining high-quality radio data to understand recurrent jet activity and feedback, with implications for future research with the SKA observatory.

The LHAASO collaboration has recently released the spectrum and the angular distribution of the $\gamma$-ray Galactic diffuse emission from 1 TeV to 1 PeV measured with the Kilometer-2 Array (KM2A) and Water Cherenkov Detector Array (WCDA). We show that these data are in remarkably good agreement with a set of models that assume the emission to be produced by the Galactic population of cosmic rays if its spectral shape traces that measured by CALET and DAMPE as well as KASCADE at higher energies. No extra-components besides the CR sea is needed to explain LHAASO results. Accounting for unresolved sources, we consistently reproduce Tibet AS$\gamma$ as well as a wide set of $\gamma$-ray data at lower energy. To do this, we consider two different transport setups: a conventional one and a $\gamma$-optimized spatial-dependent one (a development of the widely adopted KRA$_\gamma$ model). We demonstrate that both setups are compatible with LHAASO results. However, the latter is preferred if one takes into account Fermi-LAT gamma-ray data and neutrino measurements. In fact, we also compute the associated Galactic neutrino diffuse emission finding that the contribution from sources cannot be dominant and showing that spatial-dependent propagation models closely match the ANTARES and IceCube best fits for the Galactic Center Ridge and the Galactic Plane emissions. We argue that our $\gamma$-optimized model should be used as a template for future analyses of upcoming data from the Global Neutrino Network.

Solar flares and coronal mass ejections (CMEs) are manifestations of energy release in the solar atmosphere, which can be accompanied by dynamic mass motions and waves in the surrounding atmosphere. Here, we present observations of plasma moving in a helical trajectory along a set of coronal loops formed following the eruption of a CME on 2024 May 14. This helical motion was observed in extreme ultraviolet (EUV) images from the Solar Dynamic Observatory (SDO), which provides new insights into plasma properties in a set of post-eruption coronal loops. We utilize images from the SDO Atmospheric Imaging Assembly (AIA) instrument to track the helical motion of plasma and to characterize its speed, acceleration, and physical properties. Additionally, we explore the evolution of the plasma density and temperature along the helical structure using the differential emission measure technique. The helical structure was visible in AIA for approximately 22 minutes, having a diameter of 22 Mm, and a total trajectory of nearly 184 Mm. Analysis of the AIA observations reveals that the plasma flow along this helical coronal loop exhibits speeds of 77-384 km s$^{-1}$ and temperatures ranging from 3.46 to 10.2 MK. Additionally, the densities were estimated to be between 4.3.106 and 1.55.107 cm-3, with an estimated magnetic field strength of 0.05-0.3 G. Following the launch of a CME, we find clear evidence for impulsive heating and expansion of plasma that travels a helical trajectory along a set of post-eruption loops. These observations provide an insight into impulsive plasma flows along coronal loops and indeed the topology of coronal loops.

We present a study the initial stages of ice growth on pristine and oxygen-functionalized highly oriented pyrolytic graphite (O-HOPG), combining low-temperature scanning tunneling microscopy (LT-STM) and machine-learning structural searches. LT-STM images show that oxygen atoms act as nucleation sites for ice growth, and that the size, structure and porosity of the nanometer-sized ice clusters depend strongly on the growth temperature. Machine learning-assisted structural searches and first-principles energy calculations confirm that clusters of water molecules are likely to bind to chemisorbed oxygen atoms through hydrogen bonding. During the early stages of the cluster growth clusters of water molecules are likely to be immobilized by binding to more than one chemisorbed oxygen atom through hydrogen bonding. However, the energy gain by hydrogen bond formation of a molecule, upon incorporation into smaller clusters only bound to a single oxygen atom, is large enough to induce cluster diffusion and favor the growth of larger ice clusters. Our results demonstrate that the mobility of water molecules is significantly lowered in the presence of defects on the surface. The observed lower mobility on defected carbon presented here offers an enhanced understanding of macroscopic anti-icing properties observed for functionalized HOPG under ambient conditions and provides insight into the early stages of ice growth on dust grain surfaces in interstellar space.

We identify galaxies hosting ionised-gas high-velocity outflows from the complete sample of medium resolution (R1000) JWST/NIRSpec MSA spectroscopy taken as part of the JWST Advanced Deep Extragalactic Survey (JADES). From a total sample of 1087 [OIII]5007 emitters we identify 40 galaxies with a blue-side broadened [OIII]5007 line at z=2.5-9. Of these, 34 are strong outflow candidates whilst 6 sources have broadening potentially driven by rotating clumps. Our outflow candidate sample is mainly composed of star-forming galaxies, including ~65 starbursts, which span the stellar mass range log10(M*/Msun)=7.5-11.0. It also includes two candidate type-2 active galactic nuclei (AGN) and a 'little red dot' (LRD). We report a median outflow velocity of 531(+146)(-159) km/s and an overall incidence rate of 3.4%. These values are significantly higher and lower respectively than recent similar works, which we accredit to the limiting resolution of the R1000 spectroscopy and a stricter outflow selection criterion. We find no correlation between the outflow velocity and the galaxy stellar mass or star-formation rate. The median ratio between outflow velocity and escape velocity is 0.77(+0.36)(-0.32), indicating that most outflows cannot escape the galaxy gravitational potentials. We do find an anti-correlation between mass loading factor and stellar mass up to M* ~(10^10)Msun, with most lowest stellar-mass M*<(10^9)Msun galaxies reaching values well above unity, as is the case for local starburst galaxies.

Cosmological observables rely heavily on summary statistics such as two-point correlation functions. In many practical cases (e.g. the weak-lensing cosmic shear), those correlation functions are estimated from a finite, discrete sample of measurements that are randomly distributed in the sky. The result then inevitably depends on the sample at hand, regardless of any experimental noise. This sample dependence is a source of uncertainty for cosmological observables which I call sparsity covariance. This article proposes a mathematical definition and a generic method to compute sparsity covariance. It is then applied to the concrete case of cosmic shear, showing that sparsity covariance mostly enhances shape noise, whose amplitude is determined by the apparent ellipticity of galaxies rather than their intrinsic ellipticity. In general, sparsity covariance is non-negligible when the signal-to-noise ratio of individual measurements in the sample is comparable to, or larger than, unity.

While dust is a key parameter of Mars climate, its behaviour from one year to the next can appear erratic. This variability is notably related to Global Dust Storms (GDS) which occur only certain years with different onset, duration and intensity. The interannual variabilities of the dust cycle may notably explain some characteristics of Recurring Slope Lineae (RSL), slope flows once thought to be caused by liquid water. Long-term monitoring of dust dynamics is thus required to better understand surface-atmosphere dust exchanges on Mars. Here we present a new method to detect atmospheric dust as a function of space and time in the OMEGA Near-InfraRed (NIR) dataset. This dataset covers more than three Martian years; it includes the 2007 GDS which seasonality differs from the preceding (2001) and later (2018) GDS. The method is based on the decrease of the atmospheric optical path caused by dust, that can be measured by OMEGA with the 2 $\mu$m CO$_2$ absorption band. This measure is converted to a 0.9 $\mu$m NIR dust optical depth using notably comparisons with Mars Exploration Rovers measurements. We derive dust optical depth maps and comment on the variability of the dust seasonal cycle before, during and after the 2007 GDS. We also compare OMEGA NIR optical depths to Thermal InfraRed (TIR) ones derived by other studies. We found a NIR/TIR dust extinction optical depth ratio of 1.8 on average, with some variations notably related to dust particle size. Finally, we show in the northern hemisphere that atmospheric dust and RSL activity is correlated. This may indicate that dust lifting or transport mechanisms working at regional scale also participate to local RSL activity.

Analysing extended emission in photometric observations of star-forming regions requires maps free from compact foreground, embedded, and background sources, which can interfere with various techniques used to characterise the interstellar medium. Within the framework of the NEMESIS project, we apply machine learning techniques to improve our understanding of the star formation timescales, which involves the unbiased analysis of the extended emission in these regions. We present a deep learning-based method for separating the signals of compact sources and extended emission in photometric observations made by the Herschel Space Observatory, facilitating the analysis of extended emission and improving the photometry of compact sources. Central to our approach is a modified U-Net architecture with partial convolutional layers. This method enables effective source removal and background estimation across various flux densities, using a series of partial convolutional layers, batch normalisation, and ReLU activation layers within blocks. Our training process utilised simulated sources injected into Herschel images, with controlled flux densities against known backgrounds. A dynamic, signal-to-noise ratio-based adaptive masking system was implemented to assess how prominently a compact source stands out from the surrounding background. The results demonstrate that our method can significantly improve the photometric accuracy in the presence of highly fluctuating backgrounds. Moreover, the approach can preserve all characteristics of the images, including the noise properties. The presented approach allows users to analyse extended emission without the interference of disturbing point sources or perform more precise photometry of sources located in complex environments. We also provide a Python tool with tutorials and examples to help the community effectively utilise this method.

Understanding the evolution of the complex magnetic fields found in solar active regions is an active area of research. There are numerous models for such fields which range in their complexity due to the number of known physical effects included in them, the one common factor being they all extrapolate the field up from the photosphere. In this study we focus on the fact that, above the photosphere, and below the corona, lies the relatively cool and dense chromosphere -- which is often neglected in coronal models due to it being comparatively thin and difficult hard to model. We isolate and examine the effect including this boundary layer has on a 2.5D class of driven MHD models of an active region eruption. We find that it can result in significant changes to the dynamics of an erupting field far higher in the atmosphere than the chromosphere itself, generally delaying eruption and increasing the magnetic energy released in each eruption. We also test whether these effects can be approximated using a variation of the more computationally efficient magnetofrictional model, finding a number of simple adaptations of the standard magnetofrictional model capture the effect the chromospheric stratification well.

OGLE-2002-BLG-360 is an example of a Galactic red nova, the aftermath of a non-compact stellar merger. The dusty nature of the circumstellar environment makes observations of this particular source difficult, meaning the properties of the central star and its surrounding environment are poorly understood. We aim to establish the parameters of the merger remnant, as well as of the dusty environment and its structure. We attempt to establish similarities with other Galactic red novae and try to understand how such an environment has been formed. We use optical, infrared, and sub-millimetre observations to construct the spectral energy distribution (SED) between 2 um and 1.27 mm for an epoch 15-21 years after the red nova eruption. We use various radiative transfer codes to model the SED and retrieve the physical properties of both the central star and the surrounding dust. We show that the SED is best replicated by a spherically symmetric model consisting of an M-type supergiant surrounded by 0.012 Msun of dust concentrated within two spherical shells. The dust extends out to a maximum distance of 9500 AU from the central source, whilst the inner shell extends out to 1000 AU. The dust composition is dominated by iron grains (58%), but also contains olivine silicates (25%) and alumina dust (17%), which are both required to reproduce the profile of the observed 10 um absorption feature. The inner shell likely originates from merger and post-merger ejecta, whilst the outer shell consists of material lost much earlier, before the merger event occurred. Evolution of the SED indicates continued dust formation within the expanding inner shell, which may be analogous to winds of red supergiants. The object is extremely dusty compared to other red nova remnants.

Coronal mass ejections (CMEs) are the eruptions of magnetised plasma from the Sun and are considered the main driver of adverse space weather events. Hence, undrstanding its formation process, particularly the magnetic topology, is critical for accurate space weather prediction. Here, based on imaging observations and three-dimensional (3D) data-constrained thermodynamic magnetohydrodynamical (MHD) simulation in spherical coordinates, we exhibit the birth of a CME with intricate magnetic structure from multiple active regions (ARs) due to 3D magnetic reconnection. It is observed as a coronal jet between active regions, accompanied by the back-flowing of filament materials along the jet spine after the passage of the eruptive filament. This jet connects two dimming regions within different active regions. This is an observational proxy of 3D magnetic reconnection between the CME flux rope and the null-point magnetic field lines crossing active regions. Hereafter, the thermodynamic data-constrained MHD simulation successfully reproduces the observed jet and the reconnection process that flux ropes partake in, leading to a CME flux rope with a complex magnetic structure distinct from its progenitor. The generality of this scenario is then validated by data-inspired MHD simulations in a simple multipolar magnetic configuration. This work demonstrates the role of multiple active regions in forming CMEs with intricate magnetic structures. On the one hand, a non-coherent flux rope where not all twisted magnetic field lines wind around one common axis is naturally formed. On the other hand, our findings suggest that the topology of a real CME flux rope may not be solely determined by a single active region, particularly during periods of solar maximum.

Using results of numerical simulations and astrophysical observations (mainly in the WMAP and Planck frequency bands) it is shown that Galactic foreground emission becomes more sensitive to the mean magnetic field with the frequency, that results in the appearance of two levels of its randomization due to chaotic/turbulent dynamics of magnetized interstellar medium dominated by the magnetic helicity. The galactic foreground emission is more randomized at higher frequencies. The Galactic synchrotron and polarized dust emissions have been studied in detail. It is shown that the magnetic field imposes its level of randomization on the synchrotron and dust emission. The main method for the theoretical consideration used in this study is the Kolmogorov-Iroshnikov phenomenology in the frames of distributed chaos notion. Despite the vast differences in the values of physical parameters and spatio-temporal scales between the numerical simulations and the astrophysical observations, there is a quantitative agreement between the results of the astrophysical observations and the numerical simulations in the frames of the distributed chaos notion.

At least some short gamma-ray bursts (SGRBs) originate from neutron star mergers, systems that release both a relativistic collimated jet and slower, wider ejecta. These jets evolve through a dense, anisotropic medium shaped by the merger process, resulting in interactions that affect their morphology and observable signatures. We investigate the propagation of SGRB jets through funnel-like structures that can be static or expanding with mildly relativistic speed using 2D axisymmetric relativistic hydrodynamic simulations. Our models are based on radial and angular distributions of density and pressure from general-relativistic magnetohydrodynamic simulations of binary neutron star mergers. We explore different values of the funnel opening angle and density contrast, observing that narrow funnels constrain jet expansion, leading to larger collimation. The polar structure of the ejecta mainly affects the jet evolution in the early stages, whereas the effect of expanding ejecta dominates in later stages. Jets propagating through a low-dense polar funnel move initially faster, while the presence of mildly relativistic ejecta maintains the outflow more collimated after the breakout. Our results highlight the importance of the structure and dynamical properties of the ejecta in shaping SGRB jets and their electromagnetic emission.

We present abundance ratio estimates of individual elements, namely C, N, Na, and the so-called alpha elements, Mg, O, Si, Ca, and Ti, for the bulge of M31. The analysis is based on long-slit, high-quality spectroscopy of the bulge, taken with the OSIRIS spectrograph at the Gran Telescopio CANARIAS (GTC). Abundance ratios, [X/Fe]s, are inferred by comparing radially binned spectra of M31 with different state-of-the-art stellar population models, averaging out results from various methods, namely full-spectral, full-index, and line-strength fitting, respectively. For the bulk of the bulge, we find that O, N, and Na are significantly enhanced compared to Fe, with abundances of about 0.3dex, followed by C, Mg, and Si, with [X/Fe] about 0.2dex, and lastly, Ti and Ca, mostly tracking Fe ([X/Fe]<0.1dex), within the error bars. Performing the same analysis on SDSS stacked spectra of early-type galaxies with different velocity dispersion, we find that the abundance pattern of the M31 bulge is very similar to that of most massive galaxies, supporting a scenario where most of the bulge formed in a fast and intense episode of star-formation.

We investigate deviations from $\Lambda$CDM by independently parameterizing modifications in the background evolution and the growth of structures. The background is characterized by two parameters, $A$ and $B$, which reduce to $\Omega_{m0}$ and $2/3$ in the $\Lambda$CDM limit, while deviations in the growth of structures are captured through a new fitting function for $f\sigma_8$ involving the growth index $\gamma$. Using recent observational datasets, we find significant evidence for departures from $\Lambda$CDM, with data suggesting a preference for modifications in the background evolution, likely due to evolving dark energy, rather than modifications to gravity. Additionally, our framework alleviates the $H_0$ tension and addresses the $\sigma_8$ and $S_8$ tensions. We also demonstrate that future high-precision RSD data could unveil correlations between background and perturbation parameters that are currently suppressed by observational uncertainties, offering deeper insights into cosmic evolution.

Studying the redshifted 21-cm signal from the the neutral hydrogen during the Epoch of Reionization and Cosmic Dawn is fundamental for understanding the physics of the early universe. One of the challenges that 21-cm experiments face is the contamination by bright foreground sources, such as Cygnus A, for which accurate spatial and spectral models are needed to minimise the residual contamination after their removal. In this work, we develop a new, high-resolution model of Cygnus A using Low Frequency Array (LOFAR) observations in the $110{-}250$ MHz range, improving upon previous models by incorporating physical spectral information through the forced-spectrum method during multi-frequency deconvolution. This approach addresses the limitations of earlier models by providing a more accurate representation of the complex structure and spectral behaviour of Cygnus A, including the spectral turnover in its brightest hotspots. The impact of this new model on the LOFAR 21-cm signal power spectrum is assessed by comparing it with both simulated and observed North Celestial Pole data sets. Significant improvements are observed in the cylindrical power spectrum along the Cygnus A direction, highlighting the importance of having spectrally accurate models of the brightest foreground sources. However, this improvement is washed out in the spherical power spectrum, where we measure differences of a few hundred mK at $k<0.63\,h\,\text{cMpc}^{-1}$, but not statistically significant. The results suggest that other systematic effects must be mitigated before a substantial impact on 21-cm power spectrum can be achieved.

Utilizing deep NIRCam/WFSS data from JWST's FRESCO program, we spectroscopically survey the 3.3 $\mu m$ aromatic and 3.4 $\mu m$ aliphatic C--H stretching emission bands of polycyclic aromatic hydrocarbon (PAH) molecules in galaxies at redshifts $z$$\sim$0.2--0.5. Unlike pre-JWST studies, largely limited to infrared (IR)-bright galaxies ($L_{\rm IR}\gtrsim10^{11}~L_\odot$) at $z\lesssim0.1$, we probe 200 galaxies down to $L_{\rm IR}$$\sim$$10^{8.5}$--$10^{10}~L_\odot$ well beyond the local Universe. The 3.3 $\mu m$ emission is detected at $\geq$3-$\sigma$ in 88 out of 187 galaxies, correlating tightly with galaxy IR luminosity and star formation rate (SFR) and confirming the 3.3 $\mu m$ PAH as a viable SFR tracer. Despite a large scatter, the 3.3 $\mu m$-to-IR luminosity ratio ($L_{3.3}/L_{\rm IR}$) exhibits a strong metallicity dependence with a drop of $L_{3.3}/L_{\rm IR}$ by a factor of $\gtrsim10$ at 12+log(O/H)$\sim$8.4--8.5 towards lower metallicities. The 3.4 $\mu m$ emission is detected in 37 out of 159 galaxies, with the 3.4 $\mu m$-to-3.3 $\mu m$ luminosity ratio ($L_{3.4}/L_{3.3}$) spanning from $\sim$0.05 to $\sim$0.58 (median $\sim$0.19), corresponding to PAH aliphatic fractions of $\sim$0.78%--8.3% (median $\sim$2.9%) in terms of fractional carbon atoms in aliphatic units. While $L_{3.4}/L_{3.3}$ does not depend significantly on redshift, stellar mass, metallicity, or galaxy morphology, it does decrease with various SFR tracers, suggesting that ultraviolet photons in active star-forming regions may strip aliphatic sidegroups from PAH molecules. Our study showcases the unique power of JWST's NIRCam/WFSS to systematically map PAH aromatic and aliphatic content in statistically significant, less-biased galaxy samples, providing critical insights into PAH chemistry and its connection to galaxy properties.

This review explores modified theories of gravity, particularly $f(R)$ gravity, as extensions to General Relativity (GR) that offer alternatives to dark energy for explaining cosmic acceleration. These models generalize the Einstein-Hilbert action to include functions of the Ricci scalar, providing new insights into cosmology and astrophysics. The detection of gravitational waves (GWs) has enabled rigorous tests of $f(R)$ gravity, as deviations in GW propagation, speed, and polarization can signal modifications to GR. Constraints on $f(R)$ models arise from LIGO-Virgo observations of binary mergers, the stochastic gravitational wave background (SGWB), and complementary tests in cosmology and weak-field regimes. Future GW detectors, such as LISA and the Einstein Telescope, will enhance sensitivity to smaller deviations from GR, necessitating advancements in theoretical modeling. Among competing theories -- including scalar-tensor, massive gravity, and Horndeski models -- $f(R)$ gravity remains a pivotal framework for understanding fundamental gravitational physics and cosmology. This review highlights key developments, challenges, and future directions in the field.

The discovery of solar neutrinos confirmed that the inner workings of the Sun generally match our theoretical understanding of the fusion process. Solar neutrinos have also played a role in discovering that neutrinos have mass and that they oscillate. We combine the latest solar neutrino data along with other oscillation data from reactors to determine the Sun's density profile. We derive constraints given the current data and show the anticipated improvements with more reactor neutrino data from JUNO constraining the true oscillation parameters and more solar neutrino data from DUNE which should provide a crucial measurement of $hep$ neutrinos.

We present the first extension of the special-relativistic Lattice-Boltzmann Method for radiative transport developed by Weih et al. (2020), to solve the radiative-transfer equation in curved spacetimes. The novel approach is based on the streaming of carefully selected photons along null geodesics and interpolating their final positions, velocities, and frequency shifts to all photons in a given velocity stencil. Furthermore, by transforming between the laboratory frame, the Eulerian frame, and the fluid frame, we are able to perform the collision step in the fluid frame, thus retaining the collision operator of the special-relativistic case with only minor modifications. As a result, with the new method we can model the evolution of the frequency-independent (grey) radiation field as it interacts with a background fluid via absorption, emission, and scattering in a curved background spacetime. Finally, by introducing a refined adaptive stencil, which is suitably distorted in the direction of propagation of the photon bundle, we can reduce the computational costs of the method while improving its performance in the optically-thin regime. A number of standard and novel tests are presented to validate the approach and exhibit its robustness and accuracy.

Ultra-heavy dark matter is a class of candidates for which direct detection experiments are ineffective due to the suppressed dark matter flux. We explore the potential of large underwater acoustic arrays, developed for ultra-high energy neutrino detection, to detect ultra-heavy dark matter. As ultra-heavy dark matter traverses seawater, it deposits energy through nuclear scattering, generating thermo-acoustic waves detectable by hydrophones. We derive the dark matter-induced acoustic pressure wave from first principles and characterise attenuation effects, including frequency-dependent modifications due to viscous and chemical relaxation effects in seawater, providing an improved framework for signal modelling. Our sensitivity analysis for a hypothetical 100 cubic kilometre hydrophone array in the Mediterranean Sea shows that such an array could probe unexplored regions of parameter space for ultra-heavy dark matter, with sensitivity to both spin-independent and spin-dependent interactions. Our results establish acoustic detection as a promising dark matter search method, paving the way for analysing existing hydrophone data and guiding future detector designs.

Using the parity doublet model (PDM) for hadronic matter and a modified Nambu-Jona-Lasinio (NJL) model for quark matter, we investigate the potential existence of two- and three-flavor quark matter in neutron star cores. Both models respect chiral symmetry, and a sharp first-order phase transition is implemented via Maxwell construction. We find stable neutron stars with quark cores within a specific parameter space that satisfies current astronomical observations. Typical neutron stars with masses around $1.4 \ M_\odot$ may possess deconfined quark matter in their centers. The hybrid star scenario with a two-flavor quark core offers enough parameter space to allow the neutron stars with large quark cores exceeding $\sim 1\ M_\odot$, and allow the early deconfinement position before $2\ \rho_0$, where $\rho_0$ is the nuclear saturation density. The observations of gravitational wave event GW170817 suggest a relatively large chiral invariant mass $m_0=600\ \rm MeV$ in the PDM for scenarios involving three-flavor quark matter cores. The maximum mass of the hybrid star with a quark core is found to be approximately $2.2\ M_\odot$ for both two- or three-flavor quark matter in their centers.

We present a novel approach for reconstructing the $f(Q)$ gravitational theory using parameterizations of the deceleration parameter or alternative options. This method enables the development of modified gravity scenarios that align with cosmological observations. We analyze two deceleration parameter models and one effective equation of state model from the literature. By leveraging the matter density parameter's asymptotic behaviour, we further refine the viable parameter space for these models. One tested model demonstrates that even slight modifications can yield viable cosmic scenarios, where the dark energy (DE) sector becomes dynamic and is entirely described instead of using a cosmological constant, we use modified gravity. These scenarios are qualitatively distinct yet quantitatively consistent with $\Lambda$CDM.

Incorporating microscopic nuclear-structure information into the discussion of bulk properties of astronomical objects such as neutron stars has always been a challenging issue in interdisciplinary nuclear astrophysics. Using the $rp$-process nucleosynthesis as an example, we studied the effective stellar $\beta^+$ and electron capture (EC) rates of eight waiting-point (WP) nuclei with realistic stellar conditions and presence of strong magnetic fields. The relevant nuclear transition strengths are provided by the projected shell model. We have found that, on average, due to the magnetic field effect, the $\beta^+$ and EC rates can increase by more than an order of magnitude for all combinations of density and temperature, as well as in each of WPs studied. We relate the onset field strength, at which the weak rates begin to increase, to nuclear structure quantities, $Q$ value or electronic chemical potential $\mu_e$. The enhanced weak rates may change considerably the lifetime of WPs, thereby modifying the current understanding of the $rp$-process.

Gravitational wave (GW) signals of astrophysical origin are typically weak. This is because gravity is a weak force, the weakest among the four forces we know of. In order to detect GW signals, one must make differential measurements of effective lengths less than a thousandth of the size of a proton. In spite of the detectors achieving extraordinary sensitivity, the detector noise typically overwhelms the signal, so that GW signals are deeply buried in the data. The challenge to the data analyst is of extracting the GW signal from the noise, that is, first deciding whether a signal is present or not then if present, measuring its parameters. However, in the search for coalescing compact binary (CBC) signals, short-duration non-Gaussian noise transients (glitches) in the detector data significantly affect the search sensitivity. $\chi^2$ discriminators are therefore employed to mitigate their effect. We show that the underlying mathematical structure of any $\chi^2$ is a vector bundle over the signal manifold $\mathcal{P}$, that is, the manifold traced out by the signal waveforms in the Hilbert space of data segments $\mathcal{D}$. The $\chi^2$ is then defined as the square of the $L_2$ norm of the data vector projected onto a finite-dimensional subspace $\mathcal{S}$ (fibre) of $\mathcal{D}$ chosen orthogonal to the triggered template waveform. Any such fibre leads to a $\chi^2$ discriminator and the full vector bundle comprising the subspaces $\mathcal{S}$ and the base manifold $\mathcal{P}$ contitute the discriminator. We show that this structure paves the way for constructing effective discriminators against different morphologies of glitches. Here we specifically demonstrate our method on blip glitches, which can be modelled as sine-Gausians, which then generates an optimal $\chi^2$ statistic for blip glitches.

The recent discovery of the stochastic gravitational-wave background via pulsar timing arrays will likely be followed by the detection of individual black hole binaries that stand out above the background. However, to confidently claim the detection of an individual binary, we need not only more and better data, but also more sophisticated analysis techniques. In this paper, we develop two new approaches that can help us more robustly ascertain if a candidate found by a search algorithm is indeed an individual supermassive black hole binary. One of these is a coherence test that directly compares the full signal model to an incoherent version of that. The other is a model scrambling approach that builds null distributions of our detection statistic and compares that with the measured value to quantify our confidence in signal coherence. Both of these rely on finding the coherence between pulsars characteristic to gravitational waves from a binary system. We test these methods on simple simulated datasets and find that they work well in correctly identifying both true gravitational waves and false positives. However, as expected for such a flexible and simple signal model, confidently identifying signal coherence is significantly harder than simply finding a candidate in most scenarios. Our analyses also indicate that the confidence with which we can identify a true signal depends not only on the signal-to-noise ratio, but also on the number of contributing pulsars and the amount of frequency evolution shown by the signal.

In a recent paper, we computed the bispectrum of primordial density perturbations in CMB to second order in the slow-roll parameters of single field inflation, and found logarithmic infrared contributions that diverge in both large physical distances and squeezed limit where one momentum vanishes. In this work, we provide an independent test of the result by checking its conservation and the validity of the consistency relation between the squeezed limit of the bispectrum and the square of the power spectrum. Despite the violation of the main assumption for its general proofs which is the finiteness of the relevant observables in these limits, we find that the identity continues to hold in the vicinity of the squeezed limit and large time.

Lorentz invariance is a fundamental symmetry of spacetime and serves as the cornerstone of modern physics, supporting the constancy of the speed of light. A crucial implication of this principle is that no particle can propagate faster than this universal speed limit. In this study, we present a stringent neutrino-based test of Lorentz invariance, utilizing the highest-energy neutrino ever detected, known as event KM3-230213A. The detection of this neutrino, with measured energy of approximately 220 PeV, allows us to establish a lower bound on the scale of second-order Lorentz invariance violation, quantified as \(\Lambda_2>5.0\times 10^{19}\) GeV at 90 \% confidence level.

We revisit constraints on sub-GeV dark matter (DM) annihilation via $s$-wave, $p$-wave, and resonance processes using current and future CMB data from Planck, FIRAS, and upcoming experiments as LiteBird, CMB-S4, PRISTINE, and PIXIE. For $s$-wave annihilation, we provide updated limits for both $e^{+}e^{-}$ and $\pi\pi$ channels, with the profile likelihood method yielding stronger constraints than marginal posterior method. In the $p$-wave case, we comprehensively conduct a model-independent inequality for the 95\% upper limits from FIRAS, PRISTINE, and PIXIE, with future experiments expected to surpass current BBN limits. For resonance annihilation, we report -- for the first time -- the $95\%$ upper limits on the decay branching ratio of the mediator particle, when resonances peak during the recombination epoch. Overall, our study highlights the complementary strengths of $\mu$-distortion and CMB anisotropies in probing sub-GeV DM annihilation.

We obtain universal relations for fluid spheres without rotation made of dark energy assuming the Extended-Chaplygin Gas equation-of-state. After integrating the relevant differential equations, we make a fit to obtain the unknown coefficients of the functions a) normalized moment of inertia versus dimensionless deformability and b) normalized moment of inertia versus factor of compactness. We find that the form of the functions does not depend on the details of the underlying equation-of-state.

We propose a novel and natural mechanism for cosmic acceleration driven by primordial black holes (PBHs) exhibiting repulsive behavior. Using a new ``Swiss Cheese'' cosmological approach, we demonstrate that this cosmic acceleration mechanism is a general phenomenon by examining three regular black hole spacetimes - namely the Hayward, Bardeen and Dymnikova spacetimes - as well as the singular de Sitter-Schwarzschild spacetime. Interestingly, by matching these black hole spacetimes with an isotropic and homogeneous expanding Universe, we obtain a phase of cosmic acceleration that ends at an energy scale characteristic to the black hole parameters or due to black hole evaporation. This cosmic acceleration mechanism can be relevant either to an inflationary phase with a graceful exit and reheating or to an early dark energy type of contribution pertinent to the Hubble tension. Remarkably, we find that ultra-light PBHs with masses $m<5\times 10^8\mathrm{g}$ dominating the energy content of the Univese before Big Bang Nucleosynthesis, can drive a successful inflationary expansion era without the use of an inflaton field. Additionally, PBHs with masses $m \sim 10^{12}\mathrm{g}$ and abundances $0.107 < \Omega^\mathrm{eq}_\mathrm{PBH}< 0.5$, slightly before matter-radiation equality, can produce a substantial amount of early dark energy, helping to alleviate the $H_0$ tension.

We study gravitational waves from supercooled cosmological first-order phase transitions. If such a transition is followed by inefficient reheating, the evolution history of the universe is modified by a period of early matter domination. This leaves an imprint on the predicted gravitational-wave spectra. Using Fisher analysis we show the parameter space in reach of upcoming gravitational wave observatories where reheating can be probed due to its impact on the stochastic background produced by the transition. We use both the simplified geometric parametrisation and the thermodynamical one explicitly including the decay rate of the field undergoing the transition as a parameter determining the spectrum. We show the expansion history following the transition can be probed provided the transition is very strong which is naturally realised in classically scale invariant models generically predicting supercooling. Moreover, in such a scenario the decay rate of the scalar undergoing the phase transition, a parameter most likely inaccessible to accelerators, can be determined through the spectrum analysis.