Papers for Friday, Oct 18 2024

The formation details of globular clusters (GCs) are still poorly understood due to their old ages and the lack of detailed observations of their formation. A large variety of models for the formation and evolution of GCs have been created to improve our understanding of their origins, based on GC properties observed at z=0. We present the first side-by-side comparison of six current GC formation models with respect to their predictions for the GC ages and formation redshifts in Milky Way (MW)-like galaxies. We find that all the models are capable of forming most of the surviving GCs at more than 10 Gyr ago, in general agreement with the observation that most GCs are old. However, the measured MW GC ages are still systematically older than those predicted in the galaxies of four of the models. Investigating the variation of modelled GC age distributions for general MW-mass galaxies, we find that some of the models predict that a significant fraction of MW-mass galaxies would entirely lack a GC population older than 10 Gyr, whereas others predict that all MW-mass galaxies have a significant fraction of old GCs. This will have to be further tested in upcoming surveys, as systems without old GCs in that mass range are currently not known. Finally, we show that the models predict different formation redshifts for the oldest surviving GCs, highlighting that models currently disagree about whether the recently observed young star clusters at high redshifts could be the progenitors of today's GCs.

Massive black hole (MBH) binaries in galactic nuclei are one of the leading sources of $\sim$ mHz gravitational waves (GWs) for future missions such as $\rm{\textit{LISA}}$. However, the poor sky localization of GW interferometers will make it challenging to identify the host galaxy of MBH mergers absent an electromagnetic counterpart. One such counterpart is the tidal disruption of a star that has been captured into mean motion resonance with the inspiraling binary. Here we investigate the production of tidal disruption events (TDEs) through capture into, and subsequent evolution in, orbital resonance. We examine the full nonlinear evolution of planar autoresonance for stars that lock in to autoresonance with a shrinking MBH binary. Capture into the 2:1 resonance is guaranteed for any realistic astrophysical parameters (given a relatively small MBH binary mass ratio), and the captured star eventually attains an eccentricity $e\approx 1$, leading to a TDE. Stellar disks can be produced around MBHs following an active galactic nucleus episode, and we estimate the TDE rates from resonant capture produced when a secondary MBH begins inspiralling through such a disk. In some cases, the last resonant TDE can occur within a decade of the eventual $\rm{\textit{LISA}}$ signal, helping to localize the GW event.

Theory and observations reveal that the circumgalactic medium (CGM) and the cosmic web at high redshifts are multiphase, with small clouds of cold gas embedded in a hot, diffuse medium. A proposed mechanism is `shattering' of large, thermally unstable, pressure-confined clouds into tiny cloudlets of size $\ell_{\rm shatter}\sim {\rm min}(c_{\rm s}t_{\rm cool})$. These cloudlets can then disperse throughout the medium like a `fog', or recoagulate to form larger clouds. We study these processes using idealized numerical simulations of thermally unstable gas clouds. We expand upon previous works by exploring the effects of cloud geometry (spheres, streams, and sheets), metallicity, and the inclusion of an ionizing UV background. We find that `shattering' is triggered by clouds losing sonic contact and rapidly imploding, leading to a reflected shock which causes the cloud to re-expand and induces Richtmyer-Meshkov instabilities at its interface. In all cases, the expansion velocity of the cold gas is of order the cold gas sound speed, $c_{\rm s,c}$, slightly smaller in sheets than in streams and spheres due to geometrical effects. After fragmentation the cloudlets experience a drag force from the surrounding hot gas, dominated by condensation rather than ram pressure due to their low relative velocity, leading to recoagulation into larger clouds. We apply our results to the case of cold streams feeding massive ($M_{\rm v}\gtrsim 10^{12}M_\odot$) high-$z$ ($z\gtrsim 2$) galaxies from the cosmic web, finding that streams are likely to shatter upon entering the CGM through the virial shock. This offers a possible explanation for the large clumping factors and covering fractions of cold gas in the CGM around such galaxies, and may be related to galaxy quenching by providing a mechanism to prevent cold streams from reaching the central galaxy.

X-ray binaries (XRBs) exhibit spectral hysteresis for luminosities in the range $10^{-2}\lesssim L/L_\mathrm{Edd}\lesssim 0.3$, with a hard X-ray spectral state that persists from quiescent luminosities up to $\gtrsim 0.3L_\mathrm{Edd}$, transitioning to a soft spectral state that survives with decreasing luminosities down to $\sim 10^{-2}L_\mathrm{Edd}$. We present a possible approach to explain this behavior based on the thermal properties of a magnetically arrested disk simulation. By post-processing the simulation to include radiative effects, we solve for all the thermal equilibrium solutions as the accretion rate, $\dot{M}$, varies along the XRB outburst. For an assumed scaling of the disk scale height and accretion speed with temperature, we find that there exists two solutions in the range of $ 10^{-3}\lesssim\dot{M}/\dot{M}_{\rm Eddington} \lesssim 0.1$ at $r=8\:r_g$ ($ 4\times10^{-2}\lesssim\dot{M}/\dot{M}_{\rm Eddington} \lesssim 0.5$ at $r=3\:r_g$) : a cold, optically thick one and a hot, optically thin one. This opens the possibility of a natural thermal hysteresis in the right range of luminosities for XRBs. We stress that our scenario for the hysteresis does not require to invoke the strong-ADAF principle nor does it require for the magnetization of the disk to change along the XRB outburst. In fact, our scenario requires a highly magnetized disk in the cold, soft state to reproduce the soft-to-hard state transition at the right luminosities. Hence, a prediction of our scenario is that there should be a jet, although possibly very weakly dissipative, in the soft state of XRBs. We also predict that if active galactic nuclei (AGN) have similar hysteresis cycles and are strongly magnetized, they should undergo a soft-to-hard state transition at much lower $L/L_\mathrm{Edd}$ than XRBs.

We investigate a novel probe of spatial geometry of the Universe through the observation of gravitational waves (GWs) induced by first order curvature perturbations. The existence of spatial curvature leaves imprints on the gravitational wave spectrum and formation of primordial black holes. Given the peaked scalar spectrum, the induced spectrum deviates from the flat space power spectrum and the deviation is dependent on the spatial curvature K and reheating temperature $T_{rh}$. For prolonged reheating and negative spatial curvature the spectrum is amplified enough and exhibits an additional peak solely due to K indicating a possible detection by future gravitational wave experiments including LISA and DECIGO. We also observe that the presence of negative spatial curvature improves the constraints on PBH formation, increasing the mass of black holes which are viable dark matter candidates.

Unresolved binaries have a strong influence on the observed parameters of stellar clusters (SCs). We quantify this influence and compute the resulting mass underestimates and stellar mass function (MF). N-body simulations of realistic SCs were used to investigate the evolution of the binary population in a SC and its tidal tails. Together with an empirically gauged stellar mass-luminosity relation, the results were then used to determine how the presence of binaries changes the photometric mass and MF of the SC and its tails as deduced from observations. Tail 1 (T1), which is the tidal tail caused by gas expulsion, contains a larger fraction of binaries than both the SC and tail 2 (T2), which forms after gas expulsion. Additionally, T1 has a larger velocity dispersion. Using the luminosity of an unresolved binary, an observer would underestimate its mass. This bias sensitively depends on the companion masses due to the structure of the stellar mass-luminosity relation. Combining the effect of all binaries in the simulation, the total photometric mass of the SC is underestimated by 15%. Dark objects (black holes/neutron stars) increase the difference between the real and observed mass of the SC further. For both the SC and the tails, the observed power-law index of the MF between a stellar mass of 0.3 and 0.7 $M_\odot$ is smaller by up to 0.2 than the real one, the real initial mass function (IMF) being steeper by this amount. This difference is larger for stars with a larger velocity dispersion or binary fraction. Since the stars formed in SCs are the progenitors of the Galactic field stars, this work suggests that the binary fractions of different populations of stars in the Galactic disc will differ as a function of the velocity dispersion. The direction of this correlation is currently unclear and a complete population synthesis will be needed to investigate this effect.

Context. The duration of star formation (SF) in globular clusters (GCs) is an essential aspect for understanding their formation. Contrary to previous presumptions that all stars above 8 M explode as core-collapse supernovae (CCSNe), recent evidence suggests a more complex scenario. Aims. We analyse iron spread observations from 55 GCs to estimate the number of CCSNe explosions before SF termination, thereby determining the SF duration. This work for the first time takes the possibility of failed CCSNe into account, when estimating the SF duration. Methods. Two scenarios are considered: one where all stars explode as CCSNe and another where only stars below 20 M lead to CCSNe, as most CCSN models predict that no failed CCSNe happen below 20 M . Results. This establishes a lower ($\approx$ 3.5 Myr) and an upper ($\approx$ 10.5 Myr) limit for the duration of SF. Extending the findings of our previous paper, this study indicates a significant difference in SF duration based on CCSN outcomes, with failed CCSNe extending SF by up to a factor of three. Additionally, a new code is introduced to compute the SF duration for a given CCSN model. Conclusions. The extended SF has important implications on GC formation, including enhanced pollution from stellar winds and increased binary star encounters. These results underscore the need for a refined understanding of CCSNe in estimating SF durations and the formation of multiple stellar populations in GCs.

Inverse Compton (IC) emission associated with the non-thermal component of the intracluster medium (ICM) has been a long sought phenomenon in cluster physics. Traditional spectral fitting often suffers from the degeneracy between the two-temperature thermal spectrum (2T) and the one-temperature plus IC power-law spectrum (1T+IC). We present a semi-supervised deep learning approach to search for IC emission in galaxy clusters. We employ a conditional autoencoder (CAE), which is based on an autoencoder with latent representations trained to constrain the thermal parameters of the ICM. The algorithm is trained and tested using synthetic NuSTAR X-ray spectra with instrumental and astrophysical backgrounds included. The training data set only contains 2T spectra, which is more common than 1T+IC spectra. Anomaly detection is performed on the validation and test datasets, consisting of 2T spectra as the normal set and 1T+IC spectra as anomalies. With a threshold anomaly score, chosen based on cross-validation, our algorithm is able to identify spectra that contain an IC component in the test dataset, with a balanced accuracy (BAcc) of 0.64, which outperforms traditional spectral fitting (BAcc = 0.55) and ordinary autoencoder (BAcc = 0.55). Traditional spectral fitting is better at identifying IC cases among true IC spectra (a better recall), while IC predictions made by CAE have a higher chance of being true IC cases (a better precision), demonstrating their mutual complement to each other.

This introductory guide aims to provide insight to new researchers in the field of cosmic microwave background (CMB) map analysis on best practices for several common procedures. I will discuss common map-modifying procedures such as masking, downgrading resolution, the effect of the beam and the pixel window function, and adding white noise. I will explore how these modifications affect the final power spectrum measured from a map. This guide aims to describe the best way to perform each of these procedures, when the different steps and measures should be carried out, and the effects of incorrectly performing or applying any of them.

We present the results of Monte Carlo simulations aimed at exploring the evolution towards energy equipartition of first- (1G) and second-generation (2G) stars in multiple-population globular clusters and how this evolution is affected by the initial differences between the spatial distributions of the two populations. Our results show that these initial differences have fundamental implications for the evolution towards energy equipartition of the two populations. We find that 2G stars, which are assumed to be initially more centrally concentrated than 1G stars, are generally characterized by a more rapid evolution towards energy equipartition. The evolution towards energy equipartition depends on the velocity dispersion component and is more rapid for the tangential velocity dispersion. The extent of the present-day differences between the degree of energy equipartition of 2G and 1G stars depends on the cluster's dynamical age and may be more significant in the tangential velocity dispersion and at intermediate distances from the cluster's center around the half-mass radius.

Astronomical solutions provide calculated orbital and rotational parameters of solar system bodies based on the dynamics and physics of the solar system. Application of astronomical solutions in the Earth sciences has revolutionized our understanding in at least two areas of active research. (i) The Astronomical (or Milankovic) forcing of climate on time scales > ~10 kyr and (ii) the dating of geologic archives. The latter has permitted the development of the astronomical time scale, widely used today to reconstruct highly accurate geological dates and chronologies. The tasks of computing vs. applying astronomical solutions are usually performed by investigators from different backgrounds, which has led to confusion and recent inaccurate results on the side of the applications. Here we review astronomical solutions and Milankovic forcing in the Earth sciences, primarily aiming at clarifying the astronomical basis, applicability, and limitations of the solutions. We provide a summary of current up-to-date and outdated astronomical solutions and their valid time span. We discuss the fundamental limits imposed by dynamical solar system chaos on astronomical calculations and geological/astrochronological applications. We illustrate basic features of chaotic behavior using a simple mechanical system, i.e., the driven pendulum. Regarding so-called astronomical "metronomes", we point out that the current evidence does not support the notion of generally stable and prominent metronomes for universal use in astrochronology and cyclostratigraphy. We also describe amplitude and frequency modulation of astronomical forcing signals and the relation to their expression in cyclostratigraphic sequences. Furthermore, the various quantities and terminology associated with Earth's axial precession are discussed in detail. Finally, we provide some suggestions regarding practical considerations.

We analyze two galaxy-scale strong gravitational lenses, SPT0418-47 and SPT2147-50, using JWST NIRCam imaging across multiple filters. To account for angular complexity in the lens mass distribution, we introduce multipole perturbations with orders $m=1, 3, 4$. Our results show strong evidence for angular mass complexity in SPT2147, with multipole strengths of 0.3-1.7 $\%$ for $m=3, 4$ and 2.4-9.5 $\%$ for $m=1$, while SPT0418 shows no such preference. We also test lens models that include a dark matter substructure, finding a strong preference for a substructure in SPT2147-50 with a Bayes factor (log-evidence change) of $\sim 60$ when multipoles are not included. Including multipoles reduces the Bayes factor to $\sim 11$, still corresponding to a $5\sigma$ detection of a subhalo with an NFW mass of $\log_{10}(M_{200}/M_{\odot}) = 10.87\substack{+0.53\\ -0.71}$. While SPT2147-50 may represent the fourth detection of a dark matter substructure in a strong lens, further analysis is needed to confirm that the signal is not due to systematics associated with the lens mass model.

The model explaining the spectral "zebra" pattern of the high-frequency interpulse (HFIP) of the Crab pulsar radio emission is proposed. The observed emission bands are diffraction fringes in the spectral domain. The pulsar's own plasma-filled magnetosphere plays a role of a frequency-dependent "diffraction screen". The observed features such as the proportional band spacing, high polarization, constant position angle, and others are explained. The model allows one to perform "tomography" of the pulsar magnetosphere. Indeed, we have obtained the plasma density profile directly from observations, without assuming a particular magnetosphere. Our model is testable and several predictions are made. The two "high-frequency components" observed at the same frequencies as the HFIP are proposed to be related to HFIP.

A new Cherenkov telescope of the SPHERE type is under development. Its main goal is to promote the solution of the problem of the primary cosmic ray mass composition at ultra high energies (1--100 PeV) using a newly developed technique of the primary mass assignment to EAS event on event-by-event basis. The telescope will carry out measurements of both the Cherenkov light reflected from the snow surface as well as the direct one. Sensitivity of the direct Cherenkov images' shapes to the primary mass is demonstrated.

We numerically investigate a crucial parameter for understanding particle acceleration theory via turbulence-induced magnetic reconnection: the particle acceleration time. We examine particles accelerated either during the jet's dynamic evolution or in a post-processing, nearly stationary regime. We derive the particle acceleration time and compare it with theoretical predictions for both the Fermi and drift regimes identified in the simulations. In the Fermi regime, the acceleration time is expected to be independent of the particles' energy, for constant reconnection velocity, as energy increases exponentially with time. Conversely, we expect the reconnection acceleration time to depend on the current sheet's thickness and the reconnection velocity, a dependence recently revisited by xu and lazarian 2023. They identified three conditions for \(t_{acc}\). We tested these relations using statistical distributions of the current sheets' thickness and reconnection velocities in the turbulent jet over time. The resulting average value of \(t_{acc}\) was found to be nearly constant with particle energy. We compared this acceleration time with the average acceleration time derived directly from 50,000 particles accelerated in situ in the same relativistic jet. When considering a longer time interval for particle acceleration in a nearly stationary snapshot of the turbulent jet, we find that the acceleration time during the Fermi regime remains nearly independent of particle energy and aligns with the acceleration time theoretical relations up to the threshold energy, attained when the particles Larmor radius becomes as large as the thickness of the largest current sheets. Beyond this threshold, the acceleration regime shifts to the slower drift regime, showing strong energy dependence, as predicted. The results also indicate a clear dominance of the Fermi regime of acceleration.

Lenticular galaxies (S0s) are formed mainly from the gas stripping of spirals in the cluster. But how S0s form and evolve in the field is still untangled. Based on spatially resolved observations from the optical Hispanic Astronomical Center in Andalusia 3.5-m telescope with the PPAK Integral Field Spectroscopy instrument and NOrthern Extended Millimeter Array, we study a dwarf (M*<10^9 Msun) S0, PGC 44685, with triple star-forming regions in the central region, namely A, B, and C, respectively. In northwest region C, we clearly detect the spectral features of Wolf-Rayet (WR) stars and quantify the WR population by stacking spectra with high WR significance. Most of the molecular gas is concentrated in the region C(WR), and there is diffuse gas around regions A and B. The WR region possesses the strongest intensities of Ha, CO(1-0), and 3mm continuum, indicating its ongoing violent star formation (gas depletion timescale $\lesssim$25 Myr) with tentative hundreds (<500) km/s stellar winds accompanied by the WR phase. Most (~96%) of three star-forming regions show relatively low metallicity distributions, suggesting possible (minor) accretions of metal-poor gas that trigger the subsequent complex star formation in a field S0 galaxy. We speculate that PGC 44685 will become quiescent in less than 30 Myr if there is no new molecular gas to provide raw materials for star formation. The existence of this dwarf star-forming S0 presents an example of star formation in the low-mass/metallicity S0 galaxy.

Crucial for star formation is the interplay between gravity and turbulence. The observed cloud virial parameter, $\alpha_{\mathrm{vir}}$, which is the ratio of twice the turbulent kinetic energy to the gravitational energy, is found to vary significantly in different environments, where the scatter among individual star-forming clouds can exceed an order of magnitude. Therefore, a strong dependence of the initial mass function (IMF) on $\alpha_{\mathrm{vir}}$ may challenge the notion of a universal IMF. To determine the role of $\alpha_{\mathrm{vir}}$ on the IMF, we compare the star-particle mass functions obtained in high-resolution magnetohydrodynamical simulations including jet and heating feedback, with $\alpha_{\mathrm{vir}}=0.0625$, $0.125$, and $0.5$. We find that varying $\alpha_{\mathrm{vir}}$ from $\alpha_{\mathrm{vir}}\sim0.5$ to $\alpha_{\mathrm{vir}}<0.1$ shifts the peak of the IMF to lower masses by a factor of $\sim2$ and increases the star formation rate by a similar factor. The dependence of the IMF and star formation rate on $\alpha_{\mathrm{vir}}$ is non-linear, with the dependence subsiding at $\alpha_{\mathrm{vir}}<0.1$. Our study shows a systematic dependence of the IMF on $\alpha_{\mathrm{vir}}$. Yet, it may not be measurable easily in observations, considering the uncertainties, and the relatively weak dependence found in this study.

Solar energetic particles (SEPs) are an important space radiation source, especially for the space weather environment in the inner heliosphere. The energy spectrum of SEP events is crucial both for evaluating their radiation effects and for understanding their acceleration process at the source region and their propagation mechanism. In this work, we investigate the properties of the SEP peak flux spectra and the fluence spectra and their potential formation mechanisms using statistical methods. We aim to advance our understanding of both SEPs' acceleration and propagation mechanisms. Employing the dataset of European Space Agency's Solar Energetic Particle Environment Modelling (SEPEM) program, we have obtained and fitted the peak-flux and fluence proton spectra of more than a hundred SEP events from 1974 to 2018. We analyzed the relationship among the solar activity, X-ray peak intensity of solar flares and the SEP spectral parameters. Based on the assumption that the initial spectrum of accelerated SEPs generally has a power-law distribution and also the diffusion coefficient has a power-law dependence on particle energy, we can assess both the source and propagation properties using the observed SEP event peak flux and fluence energy spectra. We confirm that SEPs' spectral properties are influenced by the solar source and the interplanetary conditions and their transportation process can be influenced by different phases of solar cycle. This study provides an observational perspective on the double power-law spectral characteristics of the SEP energy spectra, showing their correlation with the adiabatic cooling and diffusion processes during the particle propagation from the Sun to the observer. This contributes to a deeper understanding of the acceleration and propagation of SEP events, in particular the possible origins of the double-power law.

Coupling black hole (BH) feeding and feedback involves interactions across vast spatial and temporal scales that is computationally challenging. Tracking gas inflows and outflows from kilo-parsec scales to the event horizon for non-spinning BHs in the presence of strong magnetic fields, Cho et al. (2023, 2024) report strong suppression of accretion on horizon scales and low (2%) feedback efficiency. In this letter, we explore the impact of these findings for the supermassive BHs M87* and Sgr A*, using high-resolution, non-cosmological, magnetohydrodynamic (MHD) simulations with the Feedback In Realistic Environments (FIRE-2) model. With no feedback, we find rapid BH growth due to "cooling flows," and for 2% efficiency feedback, while accretion is suppressed, the rates still remain higher than constraints from Event Horizon Telescope (EHT) data (Event Horizon Telescope Collaboration et al. 2021, 2022) for M87* and Sgr A*. To match EHT observations of M87*, a feedback efficiency greater than 15% is required, suggesting the need to include enhanced feedback from BH spin. Similarly, a feedback efficiency of $>15\%$ is needed for Sgr A* to match the estimated observed star formation rate of $\lesssim 2 {\rm M_\odot}$ yr$^{-1}$. However, even with 100% feedback efficiency, the accretion rate onto Sgr A* matches with EHT data only on rare occasions in the simulations, suggesting that Sgr A* is likely in a temporary quiescent phase currently. Bridging accretion and feedback across scales, we conclude that higher feedback efficiency, possibly due to non-zero BH spin, is necessary to suppress "cooling flows" and match observed accretion and star formation rates in M87* and Sgr A*.

Origin of energetic upstream ions propagating towards the Sun from the Earth's bow shock is not understood clearly. In this letter, relationship between solar wind suprathermal and upstream ions has been investigated by analyzing fluxes of H, 4He, and CNO obtained from multidirectional in-situ measurements at the first Lagrange point of the Sun-Earth system during 2012-2014. 49 upstream events have been selected based on flux enhancements of the upstream ions in comparison with the solar wind suprathermal ions. An energy cut-off at less than 300 keV is observed for the upstream events. This is attributed to the efficacy of the particle acceleration process near the bow shock. Interestingly, spectra of upstream ions soften systematically as compared to the spectra of their solar wind counterpart with decreasing mass of elements. The degree of spectral softening increases with decreasing mass-to-charge ratio of the species. Since during most of the events the interplanetary magnetic field was radial, we argue that cross-field diffusion of upstream ions gives rise to the modulation (spectral softening) of upstream ions, which is dependent on the mass-to-charge ratio of species. Our work indicates towards a systematic change in solar wind suprathermal ions after interaction with the bow shock.

We present results from a high-resolution interstellar turbulence simulation and show that it closely reproduces recent $Planck$ measurements. Our model captures the scaling of $EE$ and $BB$ spectra, and the $EE/BB$ ratio in the inertial range. The PDF of the dust polarization fraction is also consistent with observations. The $TE$ cross-correlation is in broad agreement with observations. This simulation provides new insights into the physical origins of the observed $E/B$ asymmetry and positive $TE$ signal, facilitating the development of advanced Galactic dust emission models for current and future CMB experiments.

The origin of Lyman Continuum (LyC) photons responsible for reionizing the universe remains a mystery, with the fraction of escaping LyC photons from galaxies at z$\sim$ 6 to 12 being highly uncertain. While direct detection of LyC photons from this epoch is hindered by absorption from the intergalactic medium, lower redshift analogs offer a promising avenue to study LyC leakage. We present Hubble Space Telescope Cosmic Origins Spectrograph (HST COS) observations of five low redshift (z$\sim$ 0.3) massive starburst galaxies, selected for their high stellar mass and weak [SII] nebular emission - an indirect tracer of LyC escape. Three of the five galaxies show LyC leakage, highlighting the reliability of weak [SII] as a tracer, especially in light of recent JWST discoveries of z $>$ 5 galaxies with similarly weak [SII] emission. The dust corrected LyC escape fractions, which represent the LyC photons that would escape in the absence of dust, range from 33% to 84%. However, the absolute escape fractions, which show the LyC photons escaping after passing through both neutral hydrogen absorption and dust attenuation, are significantly lower, ranging between 1% and 3%. This suggests that while the galaxies are nearly optically thin to HI, their high dust content significantly suppresses LyC photon escape. These [SII] weak, massive leakers are distinct from typical low-redshift LyC emitters, showing higher metallicity, lower ionization states, more dust extinction and higher star formation surface densities. This suggests that these galaxies constitute a distinct population, likely governed by a different mechanism facilitating LyC photon escape. We propose that the feedback-driven winds in these compact starbursts create ionized channels through which LyC photons escape, aligning with a picket-fence model.

Metallicities of both gas and stars decline toward large radii in spiral galaxies, a trend known as the radial metallicity gradient. We quantify the evolution of the metallicity gradient in the Milky Way as traced by APOGEE red giants with age estimates from machine learning algorithms. Stars up to ages of $\sim$9 Gyr follow a similar relation between metallicity and Galactocentric radius. This constancy challenges current models of Galactic chemical evolution, which typically predict lower metallicities for older stellar populations. Our results favor an equilibrium scenario, in which the gas-phase gradient reaches a nearly constant normalization early in the disk lifetime. Using a fiducial choice of parameters, we demonstrate that one possible origin of this behavior is an outflow that more readily ejects gas from the interstellar medium with increasing Galactocentric radius. A direct effect of the outflow is that baryons do not remain in the interstellar medium for long, which causes the ratio of star formation to accretion, $\dot{\Sigma}_\star / \dot{\Sigma}_\text{in}$, to quickly become constant. This ratio is closely related to the local equilibrium metallicity, since its numerator and denominator set the rates of metal production by stars and hydrogen gained through accretion, respectively. Building in a merger event results in a perturbation that evolves back toward the equilibrium state on $\sim$Gyr timescales. Under the equilibrium scenario, the radial metallicity gradient is not a consequence of the inside-out growth of the disk but instead reflects a trend of declining $\dot{\Sigma}_\star / \dot{\Sigma}_\text{in}$ with increasing Galactocentric radius.

The cloud cover and meteorological parameters serve as fundamental criteria for the qualification of an astronomical observatory working in optical and infrared wavelengths. In this paper, we present a systematic assessment of key meteorological parameters at the Lenghu site. The datasets adopted in this study includes the meteorological parameters collected at the local weather stations at the site and in the Lenghu Town, the sky brightness at the local zenith acquired by the Sky Quality Meters and night sky all-sky images from a digital camera, the ERA5 reanalysis database and global climate monitoring data. From 2019 to 2023, the fractional observable time of photometric condition is 69.70%, 74.97%, 70.26%, 74.27% and 65.12%, respectively. The fractional observing time is inversely correlated with surface air temperature, relative humidity, precipitable water vapor, and dew temperature, demonstrating that the observing conditions are influenced by these meteorological parameters. Large-scale air-sea interactions affect the climate at Lenghu site, which in fact delivers a clue to understand the irregularity of 2023. Specifically, precipitable water vapor at Lenghu site is correlated to both the westerly wind index and the summer North Atlantic Oscillation index, the yearly average temperature of Lenghu site is observed to increase significantly during the occurrence of a strong El Niño event and the relative humidity anomaly at Lenghu site is correlated to the Pacific Decadal Oscillation index. The decrease of fractional observing time in 2023 was due to the ongoing strong El Niño event and relevant global climate change. We underscore the substantial role of global climate change in regulating astronomical observing conditions and the necessity for long-term continuous monitoring of the astronomical meteorological parameters at Lenghu site.

Magnetic fields (B-fields) are ubiquitous in the interstellar medium (ISM), and they play an essential role in the formation of molecular clouds and subsequent star formation. However, B-fields in interstellar environments remain challenging to measure, and their properties typically need to be inferred from dust polarization observations over multiple physical scales. In this work, we seek to use a recently proposed approach called the Velocity Gradient Technique (VGT) to study B-fields in star-forming regions and compare the results with dust polarization observations in different wavelengths. The VGT is based on the anisotropic properties of eddies in magnetized turbulence to derive B-field properties in the ISM. We investigate that this technique is synergistic with dust polarimetry when applied to a turbulent diffused medium for the purpose of measuring its magnetization. Specifically, we use the VGT on molecular line data toward the NGC~1333 star-forming region ($\rm ^{12}CO$, $\rm ^{13}CO$, $\rm C^{18}O$, and $\rm N_{2}H^{+}$), and we compare the derived B-field properties with those inferred from 214 and 850~$\mu$m dust polarization observations of the region using SOFIA/HAWC+ and JCMT/POL-2, respectively. We estimate both the inclination angle and the 3D Alfvénic Mach Number $M_A$ from the molecular line gradients. Crucially, testing this technique on gravitationally bound, dynamic, and turbulent regions, and comparing the results with those obtained from polarization observations at different wavelength, such as the plane-of-the-sky field orientation, is an important test on the applicability of the VGT in various density regimes of the ISM.

Bars are remarkable stellar structures that can transport gas toward centers and drive the secular evolution of galaxies. In this context, it is important to locate dynamical resonances associated with bars. For this study, we used ${Spitzer}$ near-infrared images as a proxy for the stellar gravitational potential and the ALMA CO(J=2-1) gas distribution from the PHANGS survey to determine the position of the main dynamical resonances associated with the bars in the PHANGS sample of 74 nearby star-forming galaxies. We used the gravitational torque method to estimate the location of the bar corotation radius ($R_{\rm CR}$), where stars and gas rotate at the same angular velocity as the bar. Of the 46 barred galaxies in PHANGS, we have successfully determined the corotation (CR) for 38 of them. The mean ratio of the $R_{\rm CR}$ to the bar radius ($R_{\rm bar}$) is $\mathcal{R} = R_{\rm CR}/R_{\rm bar} = 1.12$, with a standard deviation of $0.39$. This is consistent with the average value expected from theory and suggests that bars are predominantly fast. We also compared our results with other bar CR measurements from the literature, which employ different methods, and find good agreement ($\rho = 0.64$). Finally, using rotation curves, we have estimated other relevant resonances such as the inner Lindblad resonance (ILR) and the outer Lindblad resonance (OLR), which are often associated with rings. This work provides a useful catalog of resonances for a large sample of nearby galaxies and emphasizes the clear connection between bar dynamics and morphology.

With data observed by the Hard X-ray Modulation Telescope (\textit{Insight}-HXMT) and the Neutron star Interior Composition Explorer (\textit {NICER}), we study low-frequency quasi-periodic oscillations (LFQPOs) of the black hole candidate MAXI J1803$-$298 during the 2021 outburst. Based on hardness intensity diagram and difference of the QPOs properties, Type-C and Type-B QPOs are found in the low-hard state and soft intermediate state, respectively. After searching for and classifying QPOs in Fourier domains, we extract the QPO component and study it with wavelet analysis. The QPO and no-QPO time intervals are separated by the confidence level, so that the S-factor, which is defined as the ratio of the QPO time interval to the total length of good time interval, is calculated. We found S-factors decrease with QPOs frequency for Type-C QPOs but stay stable around zero for Type-B QPOs. The relation of S-factor of Type-C QPOs and photon energy, the correlation of S-factor and counts are also studied. Different correlation of S-factor and counts for different energy bands indicates different origins of QPOs in high energy and low energy bands, which may be due to a dual-corona scenario.

The brightest ever gamma-ray burst (GRB) 221009A displays a significant emission line component around 10 MeV. As the GRB central engine is neutron-rich, we propose that the emission line could be originally due to the 2.223 MeV gamma-rays following neutron capture with protons. The measured line profile can be adequately fitted with a neutron capture model that involves thermal broadening and a bulk Doppler shift. The spectral modeling reveals a Doppler factor varying from 5.1 to 2.1 for the neutron-rich component, along with a temperature increase from 300 keV to about 900 keV, during the time interval of 280--360 s since the trigger, with about $10^{-2}$ $M_\odot$ deuteriums produced in the process. We argue that neutron capture can take place in the outer shell of a structure jet. Disk winds could be another possible site.

The Gamma-Ray Integrated Detectors (GRID) are a space science mission that employs compact gamma-ray detectors mounted on NanoSats in low Earth orbit (LEO) to monitor the transient gamma-ray sky. Owing to the unpredictability of the time and location of gamma-ray bursts (GRBs), obtaining the photon responses of gamma-ray detectors at various incident angles is important for the scientific analysis of GRB data captured by GRID detectors. For this purpose, a dedicated Monte Carlo simulation framework has been developed for GRID detectors. By simulating each GRID detector and the NanoSat carrying it, the spectral energy response, detection efficiency, and other angular responses of each detector for photons with different incident angles and energies can be obtained within this framework. The accuracy of these simulations has been corroborated through on-ground calibration, and the derived angular responses have been successfully applied to the data analysis of recorded GRBs.

Can a naked singularity (NkS) be distinguished from a black hole (BH)? We have investigated it with cutting-edge general relativistic magneto-hydrodynamic (GRMHD) simulations, followed by general relativistic radiation transfer (GRRT) calculation for magnetized accretion flow around NkS and BH. Based on our simulations, the accreting matter close enough to the singularity repels due to effective potential. This prevents matter from reaching a NkS and forms a quasi-spherical symmetric density distribution around it, unlike the accretion flows around a BH. We observe an order of magnitude higher mass flux through the jet and much stronger wind from a NkS than a BH. We found that the jet launching mechanism in a NkS differs significantly from that in a BH. In the horizon-scale images, a NKs shows a photon arc instead of a photon ring that is shown around a BH. In summary, the flow dynamics and radiative properties around an NkS are distinctly different from a BH. These properties would be useful to either confirm or rule out such exotic compact objects through future observations.

To investigate how the radio-identified active galactic nuclei (AGN) fraction varies with cluster-centric radius, we present the projected and de-projected distributions of a large sample of LOFAR-identified radio AGN out to $30R_{500}$ around galaxy clusters. The AGN fraction experiences a $\sim 25\%$ increase above the field fraction in the cluster outskirts at around $10R_{500}$, a $\sim 20\%$ decrease around $\sim 0.5R_{500}$, and an increase of over three times the field fraction value in the very cluster core. We label these three radial windows the outer, intermediate and inner regions respectively, and investigate how these radial trends might arise due to intrinsic properties of the AGN population. The only difference seen in host galaxy stellar mass is in the inner region, where there is a much higher fraction of massive host galaxies. Analysing AGN radio luminosity, regions with a higher AGN fraction tend to have more radio luminous AGN, and vice versa. We discuss the physical mechanisms that might be responsible for these results with reference to the literature.

The increased brightness temperature of young rocky protoplanets during their magma ocean epoch makes them potentially amenable to atmospheric characterization to distances from the solar system far greater than thermally equilibrated terrestrial exoplanets, offering observational opportunities for unique insights into the origin of secondary atmospheres and the near surface conditions of prebiotic environments. The Large Interferometer For Exoplanets (LIFE) mission will employ a space-based mid-infrared nulling interferometer to directly measure the thermal emission of terrestrial exoplanets. Here, we seek to assess the capabilities of various instrumental design choices of the LIFE mission concept for the detection of cooling protoplanets with transient high-temperature magma ocean atmospheres, in young stellar associations in particular. Using the LIFE mission instrument simulator (LIFEsim) we assess how specific instrumental parameters and design choices, such as wavelength coverage, aperture diameter, and photon throughput, facilitate or disadvantage the detection of protoplanets. We focus on the observational sensitivities of distance to the observed planetary system, protoplanet brightness temperature using a blackbody assumption, and orbital distance of the potential protoplanets around both G- and M-dwarf stars. Our simulations suggest that LIFE will be able to detect (S/N $\geq$ 7) hot protoplanets in young stellar associations up to distances of $\approx$100 pc from the solar system for reasonable integration times (up to $\sim$hours). Detection of an Earth-sized protoplanet orbiting a solar-sized host star at 1 AU requires less than 30 minutes of integration time. M-dwarfs generally need shorter integration times. The contribution from wavelength regions $<$6 $\mu$m is important for decreasing the detection threshold and discriminating emission temperatures.

In this work we present source-tailored WISE mid-infrared photometry (at 3.4$\mu$m, 4.6$\mu$m, 12$\mu$m, and 23$\mu$m) of 2812 galaxies in the extended Spitzer Survey of Stellar Structure in Galaxies (S$^4$G) sample, and characterise the mid-infrared colors and dust properties of this legacy nearby galaxy data set. Informed by the relative emission between W3 (12$\mu$ m) and W4 (23$\mu$ m), we re-derive star formation rate (SFR) scaling relations calibrated to L$_{\rm TIR}$, which results in improved agreement between the two tracers. By inverse-variance weighting the W3 and W4-derived SFRs, we generate a combined mid-infrared SFR that is a broadly robust measure of star formation activity in dusty, star-forming galaxies in the nearby Universe. In addition, we investigate the use of a W3-derived dust density metric, $\Sigma_{\rm 12\mu m}$ (L$_\odot$/kpc$^2$), to estimate the SFR deficit of low mass, low dust galaxies. This is achieved by combining WISE with existing GALEX ultraviolet (UV) photometry, which we further use to explore the relationship between dust and UV emission as a function of morphology. Finally, we use our derived SFR prescriptions to examine the location of galaxies in the log SFR - log M$_\textrm{stellar}$ plane, as a function of morphological type, which underscores the complexity of dust-derived properties seen in galaxies of progressively earlier type.

Supernova remnants (SNRs) are considered as the most promising source class to account for the bulk of the Galactic cosmic-ray flux. Yet amongst the population of ultra-high energy (UHE) sources that has recently emerged, due to high-altitude particle detector experiments such as LHAASO and HAWC, remarkably few are associated with known SNRs. These observations might well indicate that the highest energy particles would escape the remnant early during the shock evolution as a result of its reduced confinement capabilities. This flux of escaping particles may then encounter dense targets (gas and dust) for hadronic interactions in the form of both atomic and molecular material such as interstellar clouds, thereby generating a UHE gamma-ray flux. We explore such a scenario here, considering known SNRs in a physically driven model for particle escape, and as coupled to molecular clouds in the Galaxy. Our analysis allows the investigation of SNR-illuminated clouds in coincidence with sources detected in the first LHAASO catalogue. Indeed, the illuminated interstellar clouds may contribute to the total gamma-ray flux from several unidentified sources, as we discuss here. Yet we nevertheless find that further detailed studies will be necessary to verify or refute this scenario of passive UHE gamma-ray sources in future.

Ongoing large stellar spectroscopic surveys of the Milky Way seek to reconstruct the major events in the assembly history of the Galaxy. Chemical and kinematic observations can be used to separate the contributions of different progenitor galaxies to the present-day stellar halo. Here we compute the number of progenitors that contribute to the accreted stellar halos of simulated Milky Way-like galaxies as a function of radius (the radial diversity) in three suites of models: Bullock & Johnston, Aquarius and Auriga. We show that there are significant differences between the predictions of these three models, beyond the halo-to-halo scatter expected in $\Lambda$CDM. Predictions of diversity from numerical simulations are sensitive to model-dependent assumptions regarding the efficiency of star formation in dwarf galaxies. We compare, at face value, to current constraints on the radial diversity of the Milky Way's accreted halo. These constraints imply that the halo of our Galaxy is dominated by $\sim2$ progenitors in the range $8-45\,\mathrm{kpc}$, in contrast to averages of $7$ progenitors in the Bullock & Johnston models, $3.5$ in Aquarius and $4.2$ in Auriga over the same region. We additionally find that the models with radial diversity most similar to that of the Milky Way are predominantly those with ongoing merger events. The Milky Way therefore appears unusual in having an accreted stellar halo dominated by a small number of progenitors accreted at very early times.

We report the discovery of a bright pulse having a dispersion measure (DM) equal to 134.4 \pm 2 pc cm^{-3}, a peak flux density (S_p) equal to 20 \pm 4 Jy and a half-width (W_e) equal to 211 \pm 6 ms. The excessive DM of the pulse, after taking into account the Milky Way contribution, is 114 pc cm^{-3} that indicates its extragalactic origin. Such value of DM corresponds to the luminosity distance 713 Mpc. The above parameters make the pulse to be a reliable candidate to the fast radio burst (FRB) event, and then it is the second FRB detected at such a large \lambda \sim 2.7 m wavelength and the first one among non-repeating FRBs. The normalized luminosity L_\nu of the event, which we have designated as FRB 20190203, estimated under assumption that the whole excessive DM is determined by the intergalactic environment toward the host galaxy, is equal to \simeq 10^{34} erg s^{-1} Hz{-1}. In addition to the study of radio data we analyzed data from the quasi-simultaneous observations of the sky in the high energy (\ge 80 keV) band by the omnidirectional detector SPI/ACS aboard the INTEGRAL orbital observatory (in order to look for a possible gamma-ray counterpart of FRB 20190203). We did not detect any transient events exceeding the background at a statistically significant level. In the INTEGRAL archive, the FRB 20190203 localization region has been observed many times with a total exposure of \sim 73.2 days. We have analyzed the data but were unable to find any reliable short gamma-ray bursts from the FRB 20190203 position. Finally we note that the observed properties of FRB 20190203 can be reproduced well in the framework of a maser synchrotron model operating in the far reverse shock (at a distance of \sim 10^{15} cm) of a magnetar. However, triggering the burst requires a high conversion efficiency (at the level of 1%) of the shock wave energy into the radio emission.

The Auger Engineering Radio Array (AERA), part of the Pierre Auger Observatory, is a facility designed to detect radio emissions from extensive air showers at high energies. Consisting of 153 autonomous radio-detector stations spread over 17 km$^2$, it detects radio waves in the frequency range of 30 to 80 MHz. Accurate characterization of the detector response is essential for correct data interpretation, previously achieved through laboratory measurements of the analog chain and measurements of the antenna's directional response. In this study, we perform an absolute calibration using the continuously monitored sidereal modulation of the diffuse Galactic radio emission. Calibration is done by comparing the average spectra recorded by the stations with seven different models of the full radio sky propagated through the system response, including antennas, filters, and amplifiers. The Galactic calibration is in good agreement with the original laboratory calibration. In addition, we analyze the time-dependence of the calibration constants over a period of seven years. No aging effects are observed in AERA stations over a timescale of a decade, which shows that radio detectors could help monitor possible aging effects of other detector systems during long-term operations and highlight their importance in determining an absolute cosmic-ray energy scale.

Non-axisymmetric structures, such as bars and spiral arms, are known to concentrate molecular gas and star formation in galaxy centers, actively building up the pseudo-bulges. However, a direct link between the neutral (i.e., molecular and atomic) gas distribution and the exerted torque forces over a broader radial range of galactic disks still remains to be explored. In the present work, we investigate this link by carefully evaluating the torque force field using the $3.6\, \mathrm{\mu m}$ images for 17 The H I Nearby Galaxy Survey (THINGS) galaxies, and measuring neutral gas distribution on resolved atomic and molecular line maps. We find that galaxies with stronger torque forces show a more concentrated neutral gas distribution over the disk-scale, defined as half the isophotal radius at $25.5\, \mathrm{mag\, arcsec^{-2}}$. The correlation holds regardless of whether the neutral gas fraction, or the effective stellar mass surface density is controlled for. In addition, $\mathrm{kpc}$-scale neutral gas over-densities tend to be located close to the local maxima of torque forces. Most of these correlations involving the torque forces are comparatively stronger than those using the traditional Fourier amplitudes to quantify the non-axisymmetric structures. These results are consistent with the scenario that non-axisymmetric structures exert torque forces, and trigger dissipative processes to transport gas inward, not only to build the pseudo-bulges, but also fuel the inner disk growth. In this regard, non-axisymmetric structures inducing stronger torque forces appear to be more efficient in these processes.

lightcurver is a photometric pipeline for time series astronomical imaging data, designed for the semi-automatic extraction of precise light curves from small, blended targets. Such targets include, but are not limited to, lensed quasars, supernovae, or Cepheids in crowded fields. lightcurver leverages STARRED (Michalewicz et al., 2023; Millon et al., 2024) to generate state-of-the-art empirical point spread function (PSF) models for each image. It then determines the relative zeropoints between epochs by combining the PSF-photometry fluxes of several stars in the field of view. Subsequently, STARRED is used again to simultaneously model the calibrated pixels of the region of interest across all epochs. This process yields light curves of the point sources and a high-resolution image model of the region of interest, cumulating the signal from all epochs. lightcurver aims to be maintainable, fast, and incremental in its processing approach. As such, it can enable the daily photometric analysis of a large number of blended targets in the context of the upcoming Rubin Observatory Legacy Survey of Space and Time.

We study the performance of the spatially-flat dynamical dark energy (DE) $w_0w_a$CDM parameterization, with redshift-dependent DE fluid equation of state parameter $w(z) = w_0 + w_a z/(1+z)$, with and without a varying CMB lensing consistency parameter $A_L$, against Planck cosmic microwave background (CMB) data (P18 and lensing) and a combination of non-CMB data composed of baryonic acoustic oscillation (BAO) measurements that do not include DESI BAO data, Pantheon+ type Ia supernovae (SNIa) observations, Hubble parameter [$H(z)$] measurements, and growth factor ($f\sigma_8$) data points. From our most restrictive data set, P18+lensing+non-CMB, for the $w_0w_a$CDM+$A_L$ parameterization, we obtain $w_0=-0.879\pm 0.060$, $w_a=-0.39^{+0.26}_{-0.22}$, the asymptotic limit $w(z\to\infty) = w_0+w_a=-1.27^{+0.20}_{-0.17}$, and $A_L=1.078^{+0.036}_{-0.040}$ (all $1\sigma$ errors). This joint analysis of CMB and non-CMB data favors DE dynamics over a cosmological constant at $\sim 1\sigma$ and $A_L>1$ at $\sim 2\sigma$, i.e. more smoothing of the Planck CMB anisotropy data than is predicted by the best-fit model. For the $w_0w_a$CDM parameterization with $A_L=1$ the evidence in favor of DE dynamics is larger, $\sim 2\sigma$, suggesting that at least part of the evidence for DE dynamics comes from the excess smoothing of the Planck CMB anisotropy data. For the $w_0w_a$CDM parameterization with $A_L=1$, there is a difference of $2.8\sigma$ between P18 and non-CMB cosmological parameter constraints and $2.7\sigma$ between P18+lensing and non-CMB constraints. When $A_L$ is allowed to vary these tensions reduced to $1.9\sigma$ and $2.1\sigma$ respectively. Our P18+lensing+non-CMB data compilation positively favors the $w_0w_a$CDM parameterization without and with a varying $A_L$ parameter over the flat $\Lambda$CDM model, and $w_0w_a$CDM+$A_L$ is also positively favored over $w_0w_a$CDM.

We utilize serendipitous observations from the Asteroid Terrestrial-impact Last Alert System (ATLAS) and the Zwicky Transient Facility (ZTF) in addition to targeted follow-up observations from the Las Cumbres Observatory (LCO) and Liverpool Telescope to analyze the first observed instance of cometary activity by the newly-discovered Jupiter-family comet C/2023 RN3 (ATLAS), whose orbital dynamics place it close to residing on a Centaur-like orbit. Across our 7-month baseline, we observe an epoch of cometary activity commencing in August 2023 with an increase in brightness of >5.4 mag. The lightcurve of 2023 RN3 indicates the presence of continuous cometary activity across our observations, suggesting the onset of a new period of sustained activity. We find no evidence of any outbursts on top of the observed brightening, nor do we find any significant color evolution across our observations. 2023 RN3 is visibly extended in LCO and Liverpool Telescope observations, indicating the presence of a spatially-extended coma. Numerical integration of 2023 RN3's orbit reveals the comet to have recently undergone a slight increase in semimajor axis due to a planetary encounter with Jupiter, however whether this orbital change could trigger 2023 RN3's cometary activity is unclear. Our estimate for the maximum dust production metric of Afrho ~400 cm is consistent with previous measurements for the Jupiter-family comet and Centaur populations.

We compute and apply the minimum spanning tree (MST) of the binary millisecond pulsar population, and discuss aspects of the known phenomenology of these systems in this context. We find that the MST effectively separates different classes of spider pulsars, eclipsing radio pulsars in tight binary systems either with a companion with a mass in the range of approximately 0.1-0.8 solar masses (redbacks) or with a companion of less than or approximately 0.06 solar masses (black widows), into distinct branches. The MST also separates black widows located in globular clusters from those found in the field and groups other pulsar classes of interest, including transitional millisecond pulsars. Using the MST and a defined ranking for similarity, we identify possible candidates likely to belong to these pulsar classes. In particular, based on this approach, we propose the black widows' classification of J1300+1240, J1630+3550, J1317-0157, J1221-0633, J1627+3219, J1737-0314A, and J1701-3006F, discuss that of J1908+2105, and analyze J1723-2837, J1431-4715, and J1902-5105 as possible transitional systems. We introduce an algorithm that quickly locates where new pulsars fall within the MST and use this to examine the positions of the transitional millisecond pulsar IGR J18245-2452 (PSR J1824-2452I), the transitional millisecond pulsar candidate 3FGL J1544.6-1125, and the accreting millisecond X-ray pulsar SAX J1808.4-3658. Assessing the positions of these sources in the MST assuming a range for their unknown variables (e.g., the spin period derivative of PSR J1824-2452I) we can effectively narrow down the parameter space necessary for searching and determining key pulsar parameters through targeted observations.

We simulate the distribution of very rare, large excursions in the primordial density field produced in models of inflation in the very early universe which include a strong enhancement of the power spectrum. The stochastic $\delta \mathcal{N}$ formalism is used to identify the probability distribution for the primordial curvature perturbation with the first-passage-time distribution, $P(\delta \mathcal{N})$, and we compare our stochastic results with those obtained in the classical $\delta \mathcal{N}$ approach. We extend the PyFPT numerical code to simulate the full 2D phase space, and apply importance sampling which allows very rare fluctuations to be simulated in $\mathcal{O}(10)$ minutes on a single CPU, where previous direct simulations required supercomputers. We demonstrate that the stochastic noise due to quantum fluctuations after a sudden transition to ultra-slow roll can be accurately modelled using an analytical Bessel-function ansatz to identify the homogeneous growing mode. The stochastic noise found in this way is a function of the field value only. This enables us to coarse grain the inflation field at the Hubble scale and include non-linear, stochastic evolution on all super-Hubble length scales.

In recent years, nanosatellites have revolutionized the space sector due to their significant economic and time-saving advantages. As a result, they have fostered the testing of advanced instruments intended for larger space science missions. The case of MELISA is presented in this work. MELISA is a magnetic measurement instrument which aims at demonstrating the in-orbit performance of AMR sensors featuring dedicated noise reduction techniques at sub-millihertz frequencies. Such low frequency ranges are relevant for future space-borne gravitational wave detectors, where the local magnetic environment of the satellite might yield a significant contribution to the overall noise budget of the observatory. The demanding magnetic noise levels required for this bandwidth, down to 0.1 mHz, make measurements arduous. To explore sensing solutions within the H2020 European Commission Programme with the involvement of ESA, the functional performance of MELISA-III will be validated in-orbit. During operations, there is the possibility to measure the low-frequency magnetic contribution stemming from orbiting the Earth's magnetic field, impeding the characterization of the intrinsic performance of the sensor. With the objective of minimizing excess noise during the in-flight operations, the present research aims to simulate the environmental magnetic conditions in LEO to identify and subtract undesired contributions to the measurements. The in-orbit long-term magnetic fluctuations are replicated using a Helmholtz coil system. A fluxgate magnetometer allows the correlation of the generated field with the payload measurements, leading to the subsequent subtraction. Proving the effect of this approach will facilitate the noise characterization of magnetic sensors in LEO, paving the way for the in-orbit validation of MELISA-III for use in magnetically demanding missions with long integration times.

We present a compendium of HI 21-cm line observations of circumstellar envelopes (CSEs) of 290 evolved stars, mostly (~84%) on the asymptotic giant branch (AGB), made with the 100m-class, single-dish Nancay Radio Telescope. The observational and data reduction procedures were optimised for separating genuine CSE HI emission from surrounding Galactic line features. For most targets (254) the results have not been previously published. Clear detections were made of 34 objects, for 33 of which the total HI flux and the size of the CSE could be determined. Possible detections were made of 21 objects, and upper limits could be determined for 95 undetected targets, while for 140 objects confusion from Galactic HI emission along the line-of-sight precluded meaningful upper limits. The collective results of this survey can provide guidance on detectability of circumstellar HI gas for future mapping and imaging studies.

We report that the gas components in the N59 bubble suffered from sequential multiple cloud-cloud collision (CCC) processes. The molecular gas in the N59 bubble can be decomposed into four velocity components, namely Cloud A [95, 108] km s$^{-1}$, Cloud B [86, 95] km s$^{-1}$, Cloud C [79, 86] km s$^{-1}$ and Cloud D [65, 79] km s$^{-1}$. Four CCC processes occurred among these four velocity components, i.e., Cloud A vs. Cloud B, Cloud A vs. Cloud C, Cloud C vs. Cloud D, and Cloud A vs. Cloud D. Using Spitzer MIR and UKIDSS NIR photometric point source catalogs, we identified 514 YSO candidates clustered in 13 YSO groups, and most of them (~60$\%$) were located at the colliding interfaces, indicating that they were mainly triggered by these four CCC processes. We also found that these four collisions occurred in a time sequential order: the earliest and most violent collision occurred between Cloud A and Cloud D about 2 Myr ago, then Cloud B collided with Cloud A about 1 Myr ago, and finally, Cloud C collided with Clouds A and D simultaneously about 0.4 Myr ago.

Thermal evolution of the central region of a $0.9 \, M_\odot$ C/O white dwarf at the initial stage of the ion mixture crystallization is studied by numerically solving the heat equation on a fine spatial and temporal grid and by including a detailed treatment of the latent heat release. Formation of two spherical shells is observed. The outer one surrounds a region where crystallization has begun. The inner one bounds a fully solidified core which has exhausted its latent heat. The region between the shells is partially liquid and partially solid. It gradually emits the latent heat of crystallization and also it releases light elements (carbon) in the process of element redistribution, accompanying the mixture solidification. Assuming that all released light elements cross the outer shell, we have estimated their flux induced by the mixture crystallization. The resulting flux is not divergent and is much smaller than an estimate derived from the growth rate of the fully crystallized core.

We apply the probabilistic hierarchical SN Ia SED model BayeSN to analyse SALT-based simulations of SNe Ia to probe consistency between the two models. This paper is the first cross-comparison of dust inference methods using SALT and BayeSN, of great importance given the history of conflicting conclusions regarding the distributions of host galaxy dust properties between the two. Overall we find that BayeSN is able to accurately recover our simulated SALT inputs, establishing excellent consistency between the two models. When inferring dust parameters with simulated samples including non-Ia contamination, we find that our choice of photometric classifier causes a bias in the inferred dust distribution; this arises because SNe Ia heavily impacted by dust are misclassified as contaminants and excluded. We then apply BayeSN to a sample of SNe from DES5YR to jointly infer host galaxy dust distributions and intrinsic differences on either side of a `mass step' at $10^{10}$ M$\odot$. We find evidence in favour of an intrinsic contribution to the mass step and a considerably smaller difference in $R_V$ distributions than most SALT-based analyses, at most $\Delta\mu_{R_V}=0.72\pm0.26$. We also build on recent results in favour of an environmental-dependence on the secondary maximum of SNe Ia in $i$-band. Twenty days post-peak, we find a offset in intrinsic $i$-band light curve between each mass bin at a significance in excess of $3\sigma$.

Stars presently identified in the bulge spheroid are probably very old, and their abundances can be interpreted as due to the fast chemical enrichment of the early Galactic bulge. The abundances of the iron-peak elements are important tracers of nucleosynthesis processes, in particular oxygen burning, silicon burning, the weak s-process, and alpha-rich freeze-out. Aims. The aim of this work is to derive the abundances of V, Cr, Mn, Co, Ni, and Cu in 58 bulge spheroid stars and to compare them with the results of a previous analysis of data from APOGEE. We selected the best lines for V, Cr, Mn, Co, Ni, and Cu located within the H-band of the spectrum, identifying the most suitable ones for abundance determination, and discarding severe blends. Using the stellar physical parameters available for our sample from the DR17 release of the APOGEE project, we derived the individual abundances through spectrum synthesis. We then complemented these measurements with similar results from different bulge field and globular cluster stars, in order to define the trends of the individual elements and compare with the results of chemical-evolution models. We verify that the H-band has useful lines for the derivation of the elements V, Cr, Mn, Co, Ni, and Cu in moderately metal-poor stars. The resulting abundances indicate that: V, Cr, and Ni vary in lockstep with Fe; Co tends to vary in lockstep with Fe, but could be showing a slight decrease with decreasing metallicity; and Mn and Cu decrease with decreasing metallicity. These behaviours are well reproduced by chemical-evolution models except for Cu, which appears to drop faster than the models predict for moderate metallicities. Finally, abundance indicators combined with kinematical and dynamical criteria appear to show that our 58 sample stars are likely to have originated in situ.

Rigid (Uniform) rotation is usually assumed when investigating the properties of mature neutron stars (NSs). Although it simplifies their description, it is an assumption because we cannot observe the NS's innermost parts. Here, we analyze the structure of NSs in the simple case of ''almost rigidity,'' where the innermost and outermost parts rotate with different angular velocities. This is motivated by the possibility of NSs having superfluid interiors, phase transitions, and angular momentum transfer during accretion processes. We show that, in general relativity, the relative difference in angular velocity between different parts of an NS induces a change in the moment of inertia compared to that of rigid rotation. The relative change depends nonlinearly on where the angular velocity jump occurs inside the NS. For the same observed angular velocity in both configurations, if the jump location is close to the star's surface-which is possible in central compact objects (CCOs) and accreting stars-the relative change in the moment of inertia is close to that of the angular velocity (which is expected due to total angular momentum aspects). If the jump occurs deep within the NS, for instance, due to phase transitions or superfluidity, smaller relative changes in the moment of inertia are observed; we found that if it is at a radial distance smaller than approximately $40\%$ of the star's radius, the relative changes are negligible. Additionally, we outline the relevance of systematic uncertainties that nonrigidity could have on some NS observables, such as radius, ellipticity, and the rotational energy budget of pulsars, which could explain the X-ray luminosity of some sources. Finally, we also show that non-rigidity weakens the universal $I$-Love-$Q$ relations.