Papers for Monday, Jun 12 2023

Diversity, equity and inclusion (DEI) programs can only be implemented successfully if proper work-life balance is possible in Heliophysics (and in STEM field in general). One of the core issues stems from the culture of "work-above-life" associated with mission concepts, development, and implementation but also the expectations that seem to originate from numerous announcements from NASA (and other agencies). The benefits of work-life balance are well documented; however, the entire system surrounding research in Heliophysics hinders or discourages proper work-life balance. For example, there does not seem to be attention paid by NASA Headquarters (HQ) on the timing of their announcements regarding how it will be perceived by researchers, and how the timing may promote a culture where work trumps personal life. The same is true for remarks by NASA HQ program officers during panels or informal discussions, where seemingly innocuous comments may give a perception that work is expected after "normal" work hours. In addition, we are calling for work-life balance plans and implementation to be one of the criteria used for down-selection and confirmation of missions (Key Decision Points: KDP-B, KDP-C).

Understanding the nature of energetic particles in the solar atmosphere is one of the most important outstanding problems in heliophysics. Flare-accelerated particles compose a huge fraction of the flare energy budget; they have large influences on how events develop; they are an important source of high-energy particles found in the heliosphere; and they are the single most important corollary to other areas of high-energy astrophysics. Despite the importance of this area of study, this topic has in the past decade received only a small fraction of the resources necessary for a full investigation. For example, NASA has selected no new Explorer-class instrument in the past two decades that is capable of examining this topic. The advances that are currently being made in understanding flare-accelerated electrons are largely undertaken with data from EOVSA (NSF), STIX (ESA), and NuSTAR (NASA Astrophysics). This is despite the inclusion in the previous Heliophysics decadal survey of the FOXSI concept as part of the SEE2020 mission, and also despite NASA's having invested heavily in readying the technology for such an instrument via four flights of the FOXSI sounding rocket experiment. Due to that investment, the instrumentation stands ready to implement a hard X-ray mission to investigate flare-accelerated electrons. This white paper describes the scientific motivation for why this venture should be undertaken soon.

The circumgalactic medium (CGM) plays a crucial role in galaxy evolution as it fuels star formation, retains metals ejected from the galaxies, and hosts gas flows in and out of galaxies. For Milky Way-type and more massive galaxies, the bulk of the CGM is in hot phases best accessible at X-ray wavelengths. However, our understanding of the CGM remains largely unconstrained due to its tenuous nature. A promising way to probe the CGM is via X-ray absorption studies. Traditional absorption studies utilize bright background quasars, but this method probes the CGM in a pencil beam, and, due to the rarity of bright quasars, the galaxy population available for study is limited. Large-area, high spectral resolution X-ray microcalorimeters offer a new approach to exploring the CGM in emission and absorption. Here, we demonstrate that the cumulative X-ray emission from cosmic X-ray background sources can probe the CGM in absorption. We construct column density maps of major X-ray ions from the Magneticum simulation and build realistic mock images of nine galaxies to explore the detectability of X-ray absorption lines arising from the large-scale CGM. We conclude that the OVII absorption line is detectable around individual massive galaxies at the $3\sigma-6\sigma$ confidence level. For Milky Way-type galaxies, the OVII and OVIII absorption lines are detectable at the $\sim\,6\sigma$ and $\sim\,3\sigma$ levels even beyond the virial radius when co-adding data from multiple galaxies. This approach complements emission studies, does not require additional exposures, and will allow probing of the baryon budget and the CGM at the largest scales.

Local universe measurements of the Hubble constant (H0) using SNe Ia with Cepheids as calibrators yield a value of H0 which is in tension with the value inferred from the CMB and other higher redshift probes. In ref. [1], the authors proposed a rapid transition in the value of the effective Newtonian gravitational constant G in order to alleviate the Hubble tension. The transition point was chosen so as to only affect distance estimates to Hubble flow SNe. However, in this study, the authors made the assumption that SNe Ia peak luminosity $L$ increases with Chandrashekhar mass $M_c$. This hypothesis contradicts a previous semi-analytic study of SN light curves in the presence of G-transition [2] which found that $L\propto M_c^{-0.97}$. Motivated by the results of refs. [1] and [2], we propose a hypothesis of a sudden recent change in the effective G at an epoch which corresponds to a look-back distance between $\sim$ 7 - 80 Mpc. A transition in G at these distances would affect both our estimate of the distances to Cepheids in calibrator galaxies, as well as to the Hubble flow supernovae. Upon fitting the observational data to this hypothesis, we find three interesting results: (i) we find mild evidence for a G-transition at 22.4 Mpc (73 million years ago) which is preferred (using certain estimators) by the calibrator type Ia SNe data over no G-transition, (ii) the H0 parameter inferred under this hypothesis is in good agreement with the value obtained from the CMB for a 4% larger value of G at earlier times, thus potentially resolving the Hubble tension, (iii) we obtain a fit to the scaling relationship between SN peak luminosity $L$ and Chandrasekhar mass $M_c$, as $L\propto M_c^{-1.68 \pm 0.68}$, which is in good agreement with the prediction of the theoretical study of ref. [2]. We also discuss how other probes could be used to verify this transition in the value of G.

As a key driver of the secular evolution of disc galaxies, bar formation is potentially linked to the surrounding tidal field. We systematically investigate the dependence of bars on both the small (${<}2\,\mathrm{Mpc}/h$) and large-scale (${>}5\,\mathrm{Mpc}/h$) tidal fields using galaxies observed between $0.01{<}z{<}0.11$ by the Sloan Digital Sky Survey (SDSS). We characterise bar strength using the ellipticity of the isophote that corresponds to the bar, $e_{\mathrm{bar}}$, derived from the galaxy image after subtracting the 2D disc component. We demonstrate the efficacy of our bar detection method by performing an extensive comparison with the visual identifications from SDSS and the DESI Legacy Surveys. Using the Yang et al. SDSS group catalogue, we confirm the results from a recent study that the average $e_{\mathrm{bar}}$ of galaxies within interacting clusters is higher than that within isolated ones at $0.01{<}z{<}0.06$, but this small-scale tidal enhancement of bars disappears after we increase the cluster sample by a factor of five to $z{=}0.11$. On large scales, we explore the dependence of $e_{\mathrm{bar}}$ on $\alpha_{5}$, the tidal anisotropy of the density field defined over $5\,\mathrm{Mpc}/h$. We do not detect any such dependence for $98\%$ of the galaxies with $\alpha_{5}{<}10$. Intriguingly, among the $2\%$ with $\alpha_{5}{\ge}10$, we detect some hint of a boost in bar strength in the underdense regions and a suppression in the overdense regions. Combining our results on both scales, we conclude that there is little evidence for the tidal dependence of bar formation in the local Universe, except for the extremely anisotropic environments.

The early spectra of the kilonova AT2017gfo have a remarkably smooth blackbody continuum, which reveals information on the thermal properties and radioactive heating within the ejecta. However, the widespread use of a single-temperature blackbody to fit kilonova data is theoretically invalid, because 1) the significant travel-time delays for a rapidly cooling surface result in a broad distribution of temperatures and 2) the relativistic Doppler correction varies across different surface elements. Thus, the observed spectrum will be a modified blackbody with a range of temperatures over the surface. In this paper we quantify the impact of these effects and illustrate the typical wavelength-dependent spectral corrections. We apply the multi-temperature blackbody framework to the first epoch X-shooter AT2017gfo spectrum, to deconvolve the underlying physical temperature at the photosphere from the relativistic Doppler shift. We show that cooling and Doppler effect individually results in a variation of temperatures over the photosphere of up to 30%, but in combination these effects nearly cancel and produce the single-temperature blackbody observed. Furthermore, we show that for more rapidly expanding ejecta, $v \approx 0.4-0.5c$, one may be able to constrain the cooling rate and the expansion velocity purely from the spectral shape.

We study the impact of resonantly scattered X-ray line emission on the observability of the hot circumgalactic medium (CGM) of galaxies. We apply a Monte Carlo radiative transfer post-processing analysis to the high-resolution TNG50 cosmological magnetohydrodynamical galaxy formation simulation. This allows us to model the resonant scattering of OVII(r) X-ray photons within the complex, multi-phase, multi-scale CGM. The resonant transition of the OVII He-like triplet is one of the brightest, and most promising, X-ray emission lines for detecting the hot CGM and measuring its physical properties. We focus on galaxies with stellar masses 10 < log(M*/Msun) < 11 at z ~ 0. After constructing a model for OVII(r) emission from the central galaxy as well as from CGM gas, we forward model these intrinsic photons to derive observable surface brightness maps. We find that scattering significantly boosts the observable OVII(r) surface brightness of the extended and diffuse CGM. This enhancement can be large -- an order of magnitude on average at a distance of 200 projected kpc for high-mass M* = 10^10.7 Msun galaxies. The enhancement is larger for lower mass galaxies, and can even reach a factor of 100, across the extended CGM. Galaxies with higher star formation rates, AGN luminosities, and central OVII(r) luminosities all have larger scattering enhancements, at fixed stellar mass. Our results suggest that next-generation X-ray spectroscopic missions including XRISM, LEM, ATHENA, and HUBS -- which aim to detect the hot CGM in emission -- could specifically target halos with significant enhancements due to resonant scattering.

A standard practice in extragalactic population studies is the fitting of parametric models to galaxy images. From such fits, key structural parameters of galaxies such as total flux and effective radius (size) can be extracted. One of the most popular parametric forms is that of the S\'ersic profile, which is flexible enough to reasonably fit the light distribution of nearly all galaxies. Here we present pysersic, a Bayesian framework created to facilitate the inference of structural parameters from galaxy images. Pysersic is written in pure Python, and is built using the jax framework, allowing for just-in-time compilation, auto-differentiation and seamless execution on CPUs, GPUs or TPUs. Inference is performed with the numpyro package using gradient based methods, e.g., No U-Turn Sampling, for efficient and robust posterior estimation in only a few minutes on a modern laptop. Pysersic is designed to have a user-friendly interface, allowing users to fit single or multiple sources in a few lines of code, while also being flexible enough for integration into current and future analysis pipelines. In addition to sampling, pysersic can produce point estimates of the best model via optimization in several seconds, and approximate the posterior via stochastic variational inference. The use of the numpyro probabilistic language provides future extensibility to arbitrary models beyond the S\'ersic.

The ionising radiation of young and massive stars is a crucial form of stellar feedback. Most ionising (Lyman-continuum; LyC, $\lambda < 912A$) photons are absorbed close to the stars that produce them, forming compact HII regions, but some escape into the wider galaxy. Quantifying the fraction of LyC photons that escape is an open problem. In this work, we present a semi-novel method to estimate the escape fraction by combining broadband photometry of star clusters from the Legacy ExtraGalactic UV Survey (LEGUS) with HII regions observed by the Star formation, Ionized gas, and Nebular Abundances Legacy Survey (SIGNALS) in the nearby spiral galaxy NGC 628. We first assess the completeness of the combined catalogue, and find that 49\% of HII regions lack corresponding star clusters as a result of a difference in the sensitivities of the LEGUS and SIGNALS surveys. For HII regions that do have matching clusters, we infer the escape fraction from the difference between the ionising power required to produce the observed HII luminosity and the predicted ionising photon output of their host star clusters; the latter is computed using a combination of LEGUS photometric observations and a stochastic stellar population synthesis code SLUG (Stochastically Lighting Up Galaxies). Overall, we find an escape fraction of $f_{esc} = 0.09^{+0.06}_{-0.06}$ across our sample of 42 HII regions; in particular, we find HII regions with high $f_{esc}$ are predominantly regions with low H$\alpha$-luminosity. We also report possible correlation between $f_{esc}$ and the emission lines [O ii]/[N ii] and [O ii]/H$\beta$.

We use the rotation curve from Gaia data release (DR) 3 to estimate the mass of the Milky Way. We consider an Einasto density profile to model the dark matter component. We extrapolate and obtain a dynamical mass $M=2.75^{+3.11}_{-0.48}\times 10^{11} M_\odot$ at $112$ kpc. This lower-mass Milky Way is consistent with the significant declining rotation curve, and can provide new insights into our Galaxy and halo inhabitants.

We investigate the potential of low-redshift Lyman alpha (Ly$\alpha$) forest for constraining active galactic nuclei (AGN) feedback models by analyzing the Illustris and IllustrisTNG simulation at z=0.1. These simulations are ideal for studying the impact of AGN feedback on the intergalactic medium (IGM) as they share initial conditions with significant differences in the feedback prescriptions. Both simulations reveal that the IGM is significantly impacted by AGN feedback. Specifically, feedback is stronger in Illustris and results in reducing cool baryon fraction to 23% relative to 39% in IllustrisTNG. However, when comparing various statistics of Ly$\alpha$ forest such as 2D and marginalized distributions of Doppler widths and H I column density, line density, and flux power spectrum with real data, we find that most of these statistics are largely insensitive to the differences in feedback models. This lack of sensitivity arises because of the fundamental degeneracy between the fraction of cool baryons and the H I photoionization rate ($\Gamma_{\rm HI}$) as their product determines the optical depth of the Ly$\alpha$ forest. Since the $\Gamma_{\rm HI}$ cannot be precisely predicted from first principles, it needs to be treated as a nuisance parameter adjusted to match the observed Ly$\alpha$ line density. After adjusting $\Gamma_{\rm HI}$, the distinctions in the considered statistics essentially fade away. Only the Ly$\alpha$ flux power spectrum at small spatial scales exhibits potentially observable differences, although this may be specific to the relatively extreme feedback model employed in Illustris. Without independent constraints on either $\Gamma_{\rm HI}$ or cool baryon fraction, constraining AGN feedback with low-redshift Ly$\alpha$ forest will be very challenging.

Accretion disk winds launched close to supermassive black holes (SMBHs) are a viable mechanism to provide feedback between the SMBH and the host galaxy. We aim to characterize the X-ray properties of the inner accretion disk wind of the nearby active galactic nucleus (AGN) PG 1126-041, and to study its connection with the ultraviolet (UV)-absorbing wind. We perform spectroscopic analysis of eight XMM-Newton observations of PG 1126-041 taken between 2004 and 2015, using both phenomenological models and the most advanced accretion disk wind models available. For half of the dataset, we can compare the X-ray analysis results with the results of quasi-simultaneous, high-resolution spectroscopic UV observations taken with the Cosmic Origins Spectrograph (COS) on board the Hubble Space Telescope. The X-ray spectra of PG 1126-041 are complex and absorbed by ionized material which is highly variable on multiple time scales, sometimes as short as 11 days. Accretion disk wind models can account for most of the X-ray spectral complexity of PG 1126-041, with the addition of massive clumps, represented by a partially covering absorber. Variations in column density ($N_H \sim 5-20 \times 10^{22}$ cm$^{-2}$) of the partially covering absorber drive the observed X-ray spectral variability of PG 1126-041. The absorption from the X-ray partially covering gas and from the blueshifted C IV troughs appear to vary in a coordinated way. The line of sight toward PG 1126-041 offers a privileged view through a highly dynamic nuclear wind originating on inner accretion disk scales, making the source a very promising candidate for future detailed studies of the physics of accretion disk winds around SMBHs.

New mid-infrared HgCdTe (MCT) detector arrays developed in collaboration with Teledyne Imaging Sensors (TIS) have paved the way for improved 10-$\mu$m sensors for space- and ground-based observatories. Building on the successful development of longwave HAWAII-2RGs for space missions such as NEO Surveyor, we characterize the first 13-$\mu$m GeoSnap detector manufactured to overcome the challenges of high background rates inherent in ground-based mid-IR astronomy. This test device merges the longwave HgCdTe photosensitive material with Teledyne's 2048x2048 GeoSnap-18 (18-$\mu$m pixel) focal plane module, which is equipped with a capacitive transimpedance amplifier (CTIA) readout circuit paired with an onboard 14-bit analog-to-digital converter (ADC). The final assembly yields a mid-IR detector with high QE, fast readout (>85 Hz), large well depth (>1.2 million electrons), and linear readout. Longwave GeoSnap arrays would ideally be deployed on existing ground-based telescopes as well as the next generation of extremely large telescopes. While employing advanced adaptive optics (AO) along with state-of-the-art diffraction suppression techniques, instruments utilizing these detectors could attain background- and diffraction-limited imaging at inner working angles <10 $\lambda/D$, providing improved contrast-limited performance compared to JWST MIRI while operating at comparable wavelengths. We describe the performance characteristics of the 13-$\mu$m GeoSnap array operating between 38-45K, including quantum efficiency, well depth, linearity, gain, dark current, and frequency-dependent (1/f) noise profile.

We present the PRIYA suite of cosmological simulations, based on the code and hydrodynamic model of the ASTRID simulation, and designed for cosmological analyses of the Lyman-$\alpha$ forest. Our simulation suite spans a $9$-dimensional parameter space, including $4$ cosmological parameters and $5$ astrophysical/thermal parameters. We have run $48$ low fidelity simulations with $1536^3$ particles in a $120$ Mpc/h box and $3$ high fidelity simulations with $3072^3$ particles in a $120$ Mpc/h box. All our simulations include a full physics model for galaxy formation, including supernova and AGN feedback, and thus also contain a realistic population of DLAs. We advance on earlier simulations suites by larger particle loads, by incorporating new physical models for patchy hydrogen and helium reionization, and by self-consistently incorporating a model for AGN feedback. We show that patchy helium reionization imprints an excess in the 1D flux power spectrum on large scales, which may allow future measurements of helium reionization bubble sizes. Simulation parameters are chosen based on a Latin hypercube design and a Gaussian process is used to interpolate to arbitrary parameter combinations. We build a multi-fidelity emulator for the 1D flux power spectrum and the mean IGM temperature. We show that our final interpolation error is $< 1\%$ and that our simulations produce a flux power spectrum converged at the percent level for $z=5.4$ - $2.2$. Our simulation suite will be used to interpret Lyman-$\alpha$ forest 1D flux power spectra from SDSS and future DESI data releases.

Tidal disruption events (TDEs) take place when a star ventures too close to a supermassive black hole (SMBH) and becomes ruptured. One of the leading proposed physical mechanisms often invoked in the literature involves weak two-body interactions experienced by the population of stars within the host SMBH's sphere of influence, commonly referred to as two-body relaxation. This process can alter the angular momentum of stars at large distances and place them into nearly radial orbits, thus driving them to disruption. On the other hand, gravitational perturbations from an SMBH companion via the eccentric Kozai-Lidov (EKL) mechanism have also been proposed as a promising stellar disruption channel. Here we demonstrate that the combination of EKL and two-body relaxation in SMBH binaries is imperative for building a comprehensive picture of the rates of TDEs. Here we explore how the density profile of the surrounding stellar distribution and the binary orbital parameters of an SMBH companion influence the rate of TDEs. We show that this combined channel naturally produces disruptions at a rate that is consistent with observations and also naturally forms repeated TDEs, where a bound star is partially disrupted on multiple orbits. Recent observations show stars being disrupted in short-period orbits, which is challenging to explain when these mechanisms are considered independently. However, the diffusive effect of two-body relaxation, combined with the secular nature of the eccentricity excitations from EKL, is found to drive stars on short eccentric orbits at a much higher rate.

Primordial black holes (PBHs) can answer a plethora of cosmic conundra, among which the origin of the cosmic magnetic fields. In particular, supermassive PBHs with masses $M_\mathrm{PBH}>10^{10} M_\odot$ and furnished with a plasma-disk moving around them can generate through the Biermann battery mechanism a seed primordial magnetic field which can later be amplified so as to provide the magnetic field threading the intergalactic medium. In this Letter, we derive the gravitational wave (GW) signal induced by the magnetic anisotropic stress of such a population of magnetised PBHs. Interestingly enough, by using GW constraints from Big Bang Nucleosynthesis (BBN) and an effective model for the galactic/turbulent dynamo amplification of the magnetic field, we set a conservative upper bound constraint on the abundances of supermassive PBHs at formation time, $\Omega_\mathrm{PBH,f}$ as a function of the their masses, namely that $\Omega_\mathrm{PBH,f}\leq 2.5\times 10^{-10}\left(\frac{M}{10^{10}M_\odot}\right)^{45/22}$. Remarkably, these constraints are comparable, and, in some mass ranges, even tighter compared to the constraints on $\Omega_\mathrm{PBH,f}$ from large-scale structure (LSS) probes; hence promoting the portal of magnetically induced GWs as a new probe to explore the enigmatic nature of supermassive PBHs.

We extend the multi-tracer (MT) formalism of the effective field theory of large-scale structure to redshift space, comparing the results of MT to a single-tracer analysis when extracting cosmological parameters from simulations. We used a sub-halo abundance matching method to obtain more realistic multi-tracer galaxy catalogs constructed from N-body simulations. Considering different values for the sample shot noise and volume, we show that the MT error bars on $A_s$, $\omega_{\rm cdm}$, and $h$ in a full-shape analysis are approximately $50\%$ smaller relative to ST. We find that cosmological and bias coefficients from MT are less degenerate, indicating that the MT parameter basis is more orthogonal. We conclude that using MT combined with perturbation theory is a robust and competitive way to accommodate the information present in the mildly non-linear scales.

In the next decade, there is an opportunity for very high return on investment of relatively small budgets by elevating the priority of smallsat funding in heliophysics. We've learned in the past decade that these missions perform exceptionally well by traditional metrics, e.g., papers/year/\$M (Spence et al. 2022 -- arXiv:2206.02968). It is also well established that there is a "leaky pipeline" resulting in too little diversity in leadership positions (see the National Academies Report at https://www.nationalacademies.org/our-work/increasing-diversity-in-the-leadership-of-competed-space-missions). Prioritizing smallsat funding would significantly increase the number of opportunities for new leaders to learn -- a crucial patch for the pipeline and an essential phase of career development. At present, however, there are far more proposers than the available funding can support, leading to selection ratios that can be as low as 6% -- in the bottom 0.5th percentile of selection ratios across the history of ROSES. Prioritizing SmallSat funding and substantially increasing that selection ratio are the fundamental recommendations being made by this white paper.

Collisionless shock waves are one of the main forms of energy conversion in space plasmas. They can directly or indirectly drive other universal plasma processes such as magnetic reconnection, turbulence, particle acceleration and wave phenomena. Collisionless shocks employ a myriad of kinetic plasma mechanisms to convert the kinetic energy of supersonic flows in space to other forms of energy (e.g., thermal plasma, energetic particles, or Poynting flux) in order for the flow to pass an immovable obstacle. The partitioning of energy downstream of collisionless shocks is not well understood, nor are the processes which perform energy conversion. While we, as the heliophysics community, have collected an abundance of observations of the terrestrial bow shock, instrument and mission-level limitations have made it impossible to quantify this partition, to establish the physics within the shock layer responsible for it, and to understand its dependence on upstream conditions. This paper stresses the need for the first ever spacecraft mission specifically designed and dedicated to the observation of both the terrestrial bow shock as well as Interplanetary shocks in the solar wind.

To fully take advantage of the data provided by large-scale structure surveys, we need to quantify the potential impact of baryonic effects, such as feedback from active galactic nuclei (AGN) and star formation, on cosmological observables. In simulations, feedback processes originate on scales that remain unresolved. Therefore, they need to be sourced via subgrid models that contain free parameters. We use machine learning to calibrate the AGN and stellar feedback models for the FLAMINGO cosmological hydrodynamical simulations. Using Gaussian process emulators trained on Latin hypercubes of 32 smaller-volume simulations, we model how the galaxy stellar mass function and cluster gas fractions change as a function of the subgrid parameters. The emulators are then fit to observational data, allowing for the inclusion of potential observational biases. We apply our method to the three different FLAMINGO resolutions, spanning a factor of 64 in particle mass, recovering the observed relations within the respective resolved mass ranges. We also use the emulators, which link changes in subgrid parameters to changes in observables, to find models that skirt or exceed the observationally allowed range for cluster gas fractions and the stellar mass function. Our method enables us to define model variations in terms of the data that they are calibrated to rather than the values of specific subgrid parameters. This approach is useful, because subgrid parameters are typically not directly linked to particular observables, and predictions for a specific observable are influenced by multiple subgrid parameters.

Collisionless shocks are fundamental processes that are ubiquitous in space plasma physics throughout the Heliosphere and most astrophysical environments. Earth's bow shock and interplanetary shocks at 1 AU offer the most readily accessible opportunities to advance our understanding of the nature of collisionless shocks via fully-instrumented, in situ observations. One major outstanding question pertains to the energy budget of collisionless shocks, particularly how exactly collisionless shocks convert incident kinetic bulk flow energy into thermalization (heating), suprathermal particle acceleration, and a variety of plasma waves, including nonlinear structures. Furthermore, it remains unknown how those energy conversion processes change for different shock orientations (e.g., quasi-parallel vs. quasi-perpendicular) and driving conditions (upstream Alfv\'enic and fast Mach numbers, plasma beta, etc.). Required to address these questions are multipoint observations enabling direct measurement of the necessary plasmas, energetic particles, and electric and magnetic fields and waves, all simultaneously from upstream, downstream, and at the shock transition layer with observatory separations at ion to magnetohydrodynamic (MHD) scales. Such a configuration of spacecraft with specifically-designed instruments has never been available, and this white paper describes a conceptual mission design -- MAKOS -- to address these outstanding questions and advance our knowledge of the nature of collisionless shocks.

The effects of metallicity on the evolution of protoplanetary disks may be studied in the outer Galaxy where the metallicity is lower than in the solar neighbourhood. We present the VLT/KMOS integral field spectroscopy in the near-infrared of $\sim$120 candidate young stellar objects (YSOs) in the CMa-$\ell$224 star-forming region located at a Galactocentric distance of 9.1 kpc. We characterise the YSO accretion luminosities and accretion rates using the hydrogen Br$\gamma$ emission and find the median accretion luminosity of $\log{(L_{\rm acc})} = -0.82^{+0.80}_{-0.82} L_\odot$. Based on the measured accretion luminosities, we investigate the hypothesis of star formation history in the CMa-$\ell$224. Their median values suggest that Cluster C, where most of YSO candidates have been identified, might be the most evolved part of the region. The accretion luminosities are similar to those observed toward low-mass YSOs in the Perseus and Orion molecular clouds, and do not reveal the impact of lower metallicity. Similar studies in other outer Galaxy clouds covering a wide range of metallicities are critical to gain a complete picture of star formation in the Galaxy.

The formation scenario of the Magellanic Bridge during an encounter between the Large and Small Magellanic Clouds $\sim200\,$Myr ago, as proposed by $N$-body models, would be imprinted in the chemical enrichment and kinematics of its stars, and sites of ongoing star formation along its extension. We present an analysis of 33 Bridge star clusters using photometry obtained with the SOAR 4-m telescope equipped with adaptive optics for the VISCACHA survey. We performed a membership selection and derived self-consistent ages, metallicities, distances and reddening values via statistical isochrone fitting, as well as tidal radii and integrated masses from structure analysis. Two groups are clearly detected: 13 well-studied clusters older than the Bridge, with $0.5-6.8\,$Gyr and $\rm{[Fe/H]}<-0.6\,$dex; and 15 clusters with $< 200\,$Myr and $\rm{[Fe/H]}>-0.5\,$dex, probably formed in-situ. The old clusters follow the overall age and metallicity gradients of the SMC, whereas the younger ones are uniformly distributed along the Bridge. The main results are as follows: $(i)$ we derive ages and metallicities for the first time for 9 and 18 clusters, respectively; $(ii)$ we detect two metallicity dips in the age-metallicity relation of the Bridge at $\sim 200\,$Myr and $1.5\,$Gyr ago for the first time, possibly chemical signatures of the formation of the Bridge and Magellanic Stream; $(iii)$ we estimate a minimum stellar mass for the Bridge of $3-5 \times 10^5\,M_\odot$; $(iv)$ we confirm that all the young Bridge clusters at $\rm{RA} < 3^h$ are metal-rich $\rm{[Fe/H]} \sim -0.4\,$dex.

Given the anisotropic emission from quasar accretion discs, their viewing angle affects estimates of the quasar luminosity, black-hole mass and Eddington ratio. Discs appear overluminous when viewed pole-on and underluminous when viewed at high inclination. In radio-quiet quasars, the viewing angle is usually unknown, although spectroscopic indicators have been proposed. Here, we use a recently discovered universality in the variability structure function (SF) of quasar light curves (LCs), where all quasars show the same SF when clocks run in units of orbital timescale. As an offset from the mean relation can be caused by incorrect orbital timescales and thus incorrect luminosities, we correlate these offsets with suggested inclination indicators. We derive SFs from NASA/ATLAS LCs spanning $\sim 6$ years of observation, using a sample of 183 luminous quasars with measured H$\beta$ lines as well as 753 quasars with CIV and MgII lines. Starting from the proposed orientation indicators, we expect quasars with narrower H$\beta$ lines and with more blueshifted CIV lines to be viewed more pole-on and thus appear overluminous. In contrast, our SF analysis finds that presumed pole-on discs appear underluminous, consistently for both line indicators. We discuss possible explanations for the behaviour of quasars with highly blueshifted CIV lines irrespective of inclination angle, including dusty outflows that might render the accretion disc underluminous and flatter disc temperature profiles with longer orbital timescales than in thin-disc models but reach no satisfying conclusion.

Dynamical perturbations from supermassive black hole (SMBH) binaries can increase the rates of tidal disruption events (TDEs). However, most previous work focuses on TDEs from the heavy black hole in the SMBH binary (SMBHB) system. In this work, we focus on the lighter black holes in SMBHB systems and show that they can experience a similarly dramatic increase in their TDE rate due to perturbations from a more massive companion. While the increase in TDEs around the more massive black hole is mostly due to chaotic orbital perturbations, we find that, around the smaller black hole, the eccentric Kozai-Lidov (EKL) mechanism is dominant and capable of producing a comparably large number of TDEs. In this instance, the mass derived from the light curve and spectra of TDEs triggered by the lighter SMBH companion are expected to be significant smaller than the SMBH mass estimated from galaxy scaling relations, which are dominated by the more massive companion. This apparent inconsistency can help find SMBHB candidates that are not currently accreting as active galactic nuclei (AGN) and that are at separations too small to be resolved as two distinct sources. TDEs thus provide us with an exciting opportunity to study SMBHB, in particular when they might be detected from galaxies hosting SMBHs that are too massive to disrupt sun-like stars.

We describe an algorithm for application of the classic `drizzle' technique to produce 3d spectral cubes using data obtained from the slicer-type integral field unit (IFU) spectrometers on board the James Webb Space Telescope. This algorithm relies upon the computation of overlapping volume elements (composed of two spatial dimensions and one spectral dimension) between the 2d detector pixels and the 3d data cube voxels, and is greatly simplified by treating the spatial and spectral overlaps separately at the cost of just 0.03% in spectrophotometric fidelity. We provide a matrix-based formalism for the computation of spectral radiance, variance, and covariance from arbitrarily dithered data and comment on the performance of this algorithm for the Mid-Infrared Instrument's Medium Resolution IFU Spectrometer (MIRI MRS). We derive a series of simplified scaling relations to account for covariance between cube spaxels in spectra extracted from such cubes, finding multiplicative factors ranging from 1.5 to 3 depending on the wavelength range and kind of data cubes produced. Finally, we discuss how undersampling produces periodic amplitude modulations in the extracted spectra in addition to those naturally produced by fringing within the instrument; reducing such undersampling artifacts below 1% requires a 4-point dithering strategy and spectral extraction radii of 1.5 times the PSF FWHM or greater.

Reliable tools to extract patterns from high-dimensionality spaces are becoming more necessary as astronomical datasets increase both in volume and complexity. Contrastive Learning is a self-supervised machine learning algorithm that extracts informative measurements from multi-dimensional datasets, which has become increasingly popular in the computer vision and Machine Learning communities in recent years. To do so, it maximizes the agreement between the information extracted from augmented versions of the same input data, making the final representation invariant to the applied transformations. Contrastive Learning is particularly useful in astronomy for removing known instrumental effects and for performing supervised classifications and regressions with a limited amount of available labels, showing a promising avenue towards \emph{Foundation Models}. This short review paper briefly summarizes the main concepts behind contrastive learning and reviews the first promising applications to astronomy. We include some practical recommendations on which applications are particularly attractive for contrastive learning.

A gravitational-wave background can be detected in pulsar-timing-array data as Hellings--Downs correlations among the timing residuals measured for different pulsars. The optimal statistic implements this concept as a classical null-hypothesis statistical test: a null model with no correlations can be rejected if the observed value of the statistic is very unlikely under that model. To address the dependence of the statistic on the uncertain pulsar noise parameters, the pulsar-timing-array community has adopted a hybrid classical--Bayesian scheme (Vigeland et al. 2018) in which the posterior distribution of the noise parameters induces a posterior distribution for the statistic. In this article we propose a rigorous interpretation of the hybrid scheme as an instance of posterior predictive checking, and we introduce a new summary statistic (the Bayesian signal-to-noise ratio) that should be used to accurately quantify the statistical significance of an observation instead of the mean posterior signal-to-noise ratio, which does not support such a direct interpretation. In addition to falsifying the no-correlation hypothesis, the Bayesian signal-to-noise ratio can also provide evidence supporting the presence of Hellings--Downs correlations. We demonstrate our proposal with simulated datasets based on NANOGrav's 12.5-yr data release. We also establish a relation between the posterior distribution of the statistic and the Bayes factor in favor of correlations, thus calibrating the Bayes factor in terms of hypothesis-testing significance.

The detection of nanoHertz gravitational waves through pulsar timing arrays hinges on identifying a common stochastic process affecting all pulsars in a correlated way across the sky. In the presence of other deterministic and stochastic processes affecting the time-of-arrival of pulses, a detection claim must be accompanied by a detailed assessment of the various physical or phenomenological models used to describe the data. In this study, we propose posterior predictive checks as a model-checking tool that relies on the predictive performance of the models with regards to new data. We derive and study predictive checks based on different components of the models, namely the Fourier coefficients of the stochastic process, the correlation pattern, and the timing residuals. We assess the ability of our checks to identify model misspecification in simulated datasets. We find that they can accurately flag a stochastic process spectral shape that deviates from the common power-law model as well as a stochastic process that does not display the expected angular correlation pattern. Posterior predictive likelihoods derived under different assumptions about the correlation pattern can further be used to establish detection significance. In the era of nanoHertz gravitational wave detection from different pulsar-timing datasets, such tests represent an essential tool in assessing data consistency and supporting astrophysical inference.

We present the results of our X-ray, ultraviolet, and optical follow-up campaigns of 1RXS J165424.6-433758, an X-ray source detected with the \textit{Swift} Deep Galactic Plane Survey (DGPS). The source X-ray spectrum (\textit{Swift} and \textit{NuSTAR}) is described by thermal bremsstrahlung radiation with a temperature of $kT=10.1\pm1.2$ keV, yielding an X-ray ($0.3-10$ keV) luminosity $L_X=(6.5\pm0.8)\times10^{31}$ erg s$^{-1}$ at a \textit{Gaia} distance of 460 pc. Spectroscopy with the Southern African Large Telescope (SALT) revealed a flat continuum dominated by emission features, demonstrating an inverse Balmer decrement, the $\lambda4640$ Bowen blend, almost a dozen HeI lines, and HeII $\lambda4541$, $\lambda4686$ and $\lambda 5411$. Our high-speed photometry demonstrates a preponderance of flickering and flaring episodes, and revealed the orbital period of the system, $P_\textrm{orb}=2.87$ hr, which fell well within the cataclysmic variable (CV) period gap between $2-3$ hr. These features classify 1RXS J165424.6-433758 as a nearby polar magnetic CV.

This document is an abbreviated version of the law review, led by Alexander Q. Gilbert, entitled: "Major Federal Actions Significantly Affecting the Quality of the Space Environment: Applying NEPA to Federal and Federally Authorized Outer Space Activities." Here, we discuss the future of the space environment, and how it is increasingly becoming a human environment with regard to continued robotic and human presence in orbit, planned and proposed robotic and human presence on bodies such as the Moon and Mars, planned space mining projects, the increase use of low-Earth orbit for communications satellites, and other human uses of space. As such, we must evaluate and protect these environments just as we do on Earth. In order to prioritize mitigating threat of contamination, avoiding conflict, and promoting sustainability in space, all to ensure that actors maintain equal and safe access to space, we propose applying the National Environmental Policy Act, or NEPA, to space missions. We put forward three examples of environmental best practices for those involved in space missions to consider: adopting precautionary and communicative structure to before, during, and after missions taking place off-world, environmental impact statements, and transparency in tools that may impact the environment (including radioisotope power sources, plans in case of vehicle loss or loss of trajectory, and others). For additional discussion related to potential space applications of NEPA, NEPA's statutory text, and NEPA's relation to space law and judicial precedent for space, we recommend reading the full law review.

The recent pre-print by Saio et al. 2023 argues that the supergiant Betelgeuse is already undergoing carbon burning, based on the assumption that all of its light variations are caused by radial pulsations. However, the angular diameter measurements of the star are in conflict with the stellar radius required by their models, as we show in this note. We discuss the feasibility that the Great Dimming was caused by constructive mode interference using long-term brightness measurements and comment on differences in modeling frameworks adopted in Saio et al. 2023 vs Joyce et al. 2020.

The phenomenon of sub-pulse drifting is an important single-pulse phenomenon that can potentially provide important insights into the elusive radio emission mechanism in pulsars. We analyze the frequency behaviour of the single pulses of B0818-41, observed from 300 to 500 MHz (Band 3 of the uGMRT), and compare it to the evolution of the average profile to place constraints on the geometry of the pulsar's emission beam. We show that a circular carousel of discrete beamlets, where each beamlet has radial symmetry, is not consistent with the observed behaviour, and describe an alternative, consistent range of possible elliptical carousel geometries. We also combine the uGMRT data with some archival MWA observations and several other published profiles to characterize the profile evolution across a frequency range spanning ~170 MHz to ~1.4 GHz.

Search for Doppler-boosting effect in flux evolution of superluminal components in blazars has been an important subject, which can help clarify their kinematic and emission properties. The kinematics and flux evolution observed at 15GHz for the three superluminal components (knot-c, -i and -k) in QSO B1308+326 (z=0.997) were investigated in detail. It is shown that the precessing jet nozzle model previously proposed by Qian et al. (1991, 2014, 2017, 2022a, 2022b) can be used to fully simulate their kinematics on pc-scales with a nozzle precession period of 16.9 yr. With the acceleration/deceleration in their motion found in the model-simulation of their kinematics we can derive their bulk Lorentz factor and Doppler factor as function of time and predict their Doppler-boosting effect. Interestingly, the flux evolution of the three superluminal components can be well interpreted in terms of their Doppler-boosting effect. The full explanation of both their kinematic behavior and flux evolution validates our precessing nozzle model and confirms that superluminal components are physical entities moving relativistically toward us at small viewing angles.

In many exoplanetary systems with `hot Jupiters', it is observed that the spin axes of host stars are highly misaligned to planetary orbital axes. In this study, a possible channel is investigated for producing such a misalignment under a hierarchical three-body system where the evolution of stellar spin is subjected to the gravitational torque induced from the planet inside Kozai--Lidov (KL) resonance. In particular, two special configurations are explored in detail. The first one corresponds to the configuration with planets at KL fixed points, and the second one corresponds to the configurations with planets moving on KL librating cycles. When the planet is located at the KL fixed point, the corresponding Hamiltonian model is of one degree of freedom and there are three branches of libration centres for stellar spin. When the planet is moving on KL cycles, the technique of Poincar\'e section is taken to reveal global structures of stellar spin in phase space. To understand the complex structures, perturbative treatments are adopted to study rotational dynamics. It shows that analytical structures in phase portraits under the resonant model can agree well with numerical structures arising in Poincar\'e sections, showing that the complicated dynamics of stellar spin are governed by the primary resonance under the unperturbed Hamiltonian model in combination with the 2:1 (high-order and/or secondary) spin-orbit resonances.

A significant excess of the stellar mass density at high redshift has been discovered from the early data release of James Webb Space Telescope ($\it{JWST}$), and it may require a high star formation efficiency. However, this will lead to large number density of ionizing photons in the epoch of reionization (EoR), so that the reionization history will be changed, which can arise tension with the current EoR observations. Warm dark matter (WDM), via the free streaming effect, can suppress the formation of small-scale structure as well as low-mass galaxies. This provides an effective way to decrease the ionizing photons when considering a large star formation efficiency in high-$z$ massive galaxies without altering the cosmic reionization history. On the other hand, the constraints on the properties of warm dark matter can be derived from the $\it JWST$ observations. In this work, we study WDM as a possible solution to reconcile the $\it JWST$ stellar mass density of high-$z$ massive galaxies and reionization history. We find that the WDM particle mass with $m_{\text{W}}\sim0.6$ keV can match both the $\it JWST$ high-$z$ comoving cumulative stellar mass density and the reionization history, for a star formation efficiency parameter range $f_{*}^{0} = 0.1 - 1.0$.

Here we investigate how small amounts of hydrogen (much smaller than the mass of the exoplanet) above a magma ocean on a rocky exoplanet may modify the atmospheric chemistry and atmospheric escape.We use a chemical model of a magma ocean coupled to a gas equilibrium code. An energy-limited model is used to compute atmospheric escape. The composition of the vapor above a magma ocean is drastically modified by hydrogen, even for very modest amounts of H ($\ll 10^{-6}$ planetary mass). Hydrogen consumes much of the O$_2$(g), which, in turn, promotes the evaporation of metals and metal oxides (SiO, Mg, Na, K, Fe) from the magma ocean. Vast amounts of H$_2$O are produced by the same process. At high hydrogen pressures, new hydrogenated species such as SiH$_4$ form in the atmosphere. In all cases, H, H$_2$, and H$_2$O are the dominant nonmetal-bearing volatile species. Sodium is the dominant atmospheric metal-bearing species at T$<$ 2000K and low H content, whereas Fe is dominant at high H content and low temperature, while SiO predominates at T>3000 K. We find that the atmospheric Mg/Fe, Mg/Si, and Na/Si ratios deviate from those in the underlying planet and from the stellar composition. As such, their determination may constrain the planet's mantle composition and H content. As the presence of hydrogen promotes the evaporation of silicate mantles, it is conceivable that some high-density, irradiated exoplanets may have started life as hydrogen-bearing planets and that part of their silicate mantle evaporated (up to a few $10 \%$ of Si, O, and Fe) and was subsequently lost owing to the reducing role of H. Even very small amounts of H can alter the atmospheric composition and promote the evaporation to space of heavy species derived from the molten silicate mantle of rocky planets.

We report the detection of three compact ($< 0.001$ pc) radio sources (CRSs) at K$_{a}$-band (0.9 cm) in the \uchiirs G040.54+2.59 (two CRSs) and G034.13+0.47 (one CRS). These CRSs have weak flux densities and are located at the center of their respective \uchii regions. We found no clear association between massive ionizing stars and CRSs but some radiative influence on the latter, as suggested by their large emission measures (> $10^7 \mathrm{cm}^{-6}\mathrm{pc}$), typical of photo evaporating neutral objects close to or associated with massive stars. Our modelling of G40.54+2.59 shows that their CRSs supply enough ionized material to shape its morphology while significantly extending its observable lifetime. On the other hand, despite the possible relation of the CRS with the large-scale outflow signatures observed in G034.13+0.47, the influence of this CRS on the evolution of the \uchii region is unlikely. Our results show that the presence of CRSs can alleviate the so-called lifetime problem of UCHII regions. Still, to address their dynamical evolution adequately, the scenario must include additional mechanisms like ambient confinement, or the role of the kinematics of their associated stellar objects.

Our recent multi-band photometric study of the colour width of the lower main sequence of the open cluster M37 has revealed the presence of a sizeable initial chemical composition spread in the cluster. If initial chemical composition spreads are common amongst open clusters, this would have major implications for cluster formation models and the foundation of the chemical tagging technique. Here we present a study of the unevolved main sequence of the open cluster M38, employing Gaia DR3 photometry and astrometry, together with newly acquired Sloan photometry. We have analysed the distribution of the cluster's lower main sequence stars with a differential colour-colour diagram made of combinations of Gaia and Sloan magnitudes, like in the study of M37. We employed synthetic stellar populations to reproduce the observed trend of M38 stars in this diagram, and found that the observed colour spreads can be explained simply by the combined effect of differential reddening across the face of the cluster and the presence of unresolved binaries. There is no need to include in the synthetic sample a spread of initial chemical composition as instead necessary to explain the main sequence of M37. Further photometric investigations like ours, as well as accurate differential spectroscopic analyses on large samples of open clusters, are necessary to understand whether chemical abundance spreads are common among the open cluster population.

The shapes of Stokes profiles contain much information about the atmospheric conditions that produced them. However, a variety of different atmospheric structures can produce very similar profiles. Thus, it is important for proper interpretation of observations to have a good understanding of how the shapes of Stokes profiles depend on the underlying atmosphere. An excellent tool in this regard is forward modeling, i.e. computing and studying synthetic spectra from realistic simulations of the solar atmosphere. Modern simulations routinely produce several hundred thousand spectral profiles per snapshot. With such numbers, it becomes necessary to use automated procedures in order to organize the profiles according to their shape. Here we illustrate the use of two complementary methods, k-means and k-Shape, to cluster similarly shaped profiles, and demonstrate how the resulting clusters can be combined with knowledge of the simulation's atmosphere to interpret spectral shapes. We generate synthetic Stokes profiles for the Ca II 854.2 nm line using the Multi3D code from a Bifrost simulation snapshot. We then apply the k-means and k-Shape clustering techniques to group the profiles together according to their shape. We show and compare the classes of profile shapes we retrieve from applying both k-means and k-Shape to our synthetic intensity spectra. We then show the structure of the underlying atmosphere for two particular classes of profile shapes retrieved by the clustering, and demonstrate how this leads to an interpretation for the formation of those profile shapes. Furthermore, we apply both methods to the subset of our profiles containing the strongest Stokes V signals, and demonstrate how k-Shape can be qualitatively better than k-means at retrieving complex profile shapes when using a small number of clusters.

We use two cosmological simulations to study the impact of spurious heating of stellar motions within simulated galaxies by dark matter (DM) particles. The simulations share the same numerical and subgrid parameters, but one used a factor of 7 more DM particles. We find that many galaxy properties are unaffected by spurious heating, including their masses, star formation histories, and the spatial distribution of their gaseous baryons. The distribution and kinematics of their stellar and DM particles, however, are affected. Below a resolution-dependent virial mass, $M_{200}^{\rm spur}$, galaxies have higher characteristic velocities, larger sizes, and more angular momentum in the simulation with lower DM mass resolution; haloes have higher central densities and lower velocity dispersions. Above $M_{200}^{\rm spur}$, galaxies and haloes have similar properties in both runs. The differences arise due to spurious heating, which transfers energy from DM to stellar particles, causing galaxies to heat up and haloes to cool down. The value of $M_{200}^{\rm spur}$ can be derived from an empirical disc heating model, and coincides with the mass below which the predicted $spurious$ velocity dispersion exceeds the $measured$ velocity dispersion of the simulated galaxies. We predict that galaxies in the $100^3\, {\rm Mpc}^3$ EAGLE run and IllustrisTNG-100 are robust to spurious collisional effects at their half-mass radii provided their halo mass exceeds $M_{200}^{\rm spur}\approx 10^{11.7}{\rm M_\odot}$; for the $25^3\, {\rm Mpc}^3$ EAGLE run and IllustrisTNG-50, we predict $M_{200}^{\rm spur}\approx 10^{11}{\rm M_\odot}$. Suppressing spurious heating at smaller/larger radii, or for older/younger stellar populations, requires haloes to be resolved with more/fewer DM particles.

The mechanism by which galaxies stop forming stars and get rid of their interstellar medium (ISM) remains elusive. Here, we study a sample of more than two thousand elliptical galaxies in which dust emission has been detected. This is the largest sample of such galaxies ever analysed. We infer the timescale for removal of dust in these galaxies and investigate its dependency on physical and environmental properties. We obtain a dust removal timescale in elliptical galaxies of $\tau$ = 2.26 $\pm$ 0.18 Gyr, corresponding to a half-life time of 1.57 $\pm$ 0.12 Gyr. This timescale does not depend on environment, stellar mass or redshift. We observe a departure of dusty elliptical galaxies from the star formation rate vs. dust mass relation. This is caused by the star-formation rates declining faster than the dust masses and indicates that there exists an internal mechanism, which affects star formation, but leaves the ISM intact. Morphological quenching together with ionisation or outflows caused by older stellar populations (supernova type Ia or planetary nebulae) are consistent with these observations.

In order to obtain a good understanding of astrochemistry, it is crucial to better understand the key parameters that govern grain-surface chemistry. For many chemical networks, these crucial parameters are the binding energies of the species. However, there exists much disagreement regarding these values in the literature. In this work, a Bayesian inference approach is taken to estimate these values. It is found that this is difficult to do in the absence of enough data. The Massive Optimised Parameter Estimation and Data (MOPED) compression algorithm is then used to help determine which species should be prioritised for future detections in order to better constrain the values of binding energies. Finally, an interpretable machine learning approach is taken in order to better understand the non-linear relationship between binding energies and the final abundances of specific species of interest.

A great portion of the cataclysmic variable population, between 40% and 70%, is predicted to be made up of period-bouncers, systems with degenerate donors that have evolved past the period minimum. However, due to their intrinsic faintness, only a few of these systems have been observed and confidently identified so far. We have searched for X-ray emission as a proof of accretion in order to confirm period-bounce cataclysmic variables. In a dedicated XMM-Newton observation of the period-bounce candidate SDSS J151415.65+074446.5 we discovered X-ray modulation at the binary orbital period confirming it as an accreting system. The X-ray light curve and the X-ray spectrum display characteristics of magnetic Polar-type systems allowing for the first time the determination of the X-ray luminosity and mass accretion rate for this system. Catalog data from eROSITA on the SRG satellite for V379 Vir and SDSS J125044.42+154957.4 allowed a first look into the X-ray behavior of period-bounce candidates with this new all-sky instrument. From the eROSITA measurements the X-ray luminosity and mass accretion rate were determined for the first time for SDSS J125044.42+154957.4, and the earlier result for V379 Vir from XMM-Newton was confirmed. All three cataclysmic variables with a magnetic white dwarf and very low-mass donor studied in this work present evidence for X-ray emission at a similar level of $L_{\rm x}\,{\rm [erg/s]} \approx 10^{29}$, which, together with the detection of X-ray orbital modulation in two of them, V379 Vir and SDSS J151415.65+074446.5, unambiguously proves the presence of accretion in these systems. The detection of these period-bouncers at faint X-ray luminosity levels with the all-sky X-ray survey eROSITA offers new prospects for the identification of additional period-bouncers, providing impetus for theoretical studies of binary evolution.

Bayesian Estimation Applied to Multiple Species (BEAMS) has been implemented in the BEAMS with Bias Corrections (BBC) framework to produce a redshift-binned Hubble diagram (HD) for Type Ia Supernovae (SN Ia). The BBC method corrects for selection effects and non-SNIa contamination, and systematic uncertainties are described by a covariance matrix with dimension matching the number of BBC redshift bins. For spectroscopically confirmed SNIa samples, a recent "Binning is Sinning" article (BHS21, arxiv:2012.05900) showed that an unbinned HD and covariance matrix reduces the systematic uncertainty by a factor of 1.5 compared to the binned approach. Here we extend their analysis to obtain an unbinned HD for a photometrically identified sample processed with BBC. To test this new method, we simulate and analyze 50 samples corresponding to the Dark Energy Survey (DES) with a low-redshift anchor; the simulation includes SNe Ia, and contaminants from core collapse SNe and peculiar SNe Ia. The analysis includes systematic uncertainties for calibration, and measures the dark energy equation of state parameter (w). Compared to a redshift-binned HD, the unbinned HD with nearly 2,000 events results in a smaller systematic uncertainty, in qualitative agreement with BHS21, and averaging results among the 50 samples we find no evidence for bias in measured cosmological parameters. To reduce computation time for fitting an unbinned HD with large samples, we propose an HD-rebinning method that defines the HD in bins of redshift, color, and stretch; the rebinned HD results in similar uncertainty as the unbinned case, and shows no evidence for biased cosmology parameters.

We infer initial masses of the synthesized radioactive nickel-56 in a sample of recent Type Ia supernovae applying a new formalism introduced recently by Khatami & Kasen (2019). It is shown that the nickel masses we derive do not differ significantly from previous estimates based on the traditional Arnett-model. We derive the $\beta$ parameter for our sample SNe and show that these are consistent with the fiducial value of $\sim 1.6$ given by Khatami & Kasen (2019) from SN Ia hydrodynamical simulations.

Due to the technical time delay of the XRT instrument on board the Neil Gehrels Swift Observatory satellite, we cannot observe the X-ray emission occurring less than $\sim 40$~s after a gamma-ray burst (GRB) trigger time. We here indicate a new strategy of using the cosmological time dilatation in high redshift GRBs to observe the earliest X-ray emission by Swift/XRT in the GRB cosmological rest-frame. We illustrate this procedure using $354$ GRBs with a well-defined cosmological redshift selected from the Swift GRB catalog. We compare and contrast the time delay between the trigger of the source and the first observation by Swift/XRT as measured in the observer frame (OTD) and the corresponding delay measured in GRBs' cosmological rest-frame (RTD). We consider as specific prototypes three binary-driven hypernovae of type I (BdHNe I): GRB 090423 at $z=8.2$ with an RTD of $8.2$~s, GRB 090429B at $z\sim 9.4$ with an RTD of $10.1$~s, as well as the GRB 220101A at $z=4.6$ with an RTD of $14.2$~s. This opens a new possibility for probing Episode (1) of BdHNe, linked to the origin and early appearance of the newborn neutron star ($\nu$NS) and its transition from a Jacobi triaxial ellipsoid (JTE) to a Maclaurin spheroid configuration that originates the GRB afterglow onset. We also present the methodology to compute the sweeping frequencies and the energetics of the associated conspicuous gravitational wave emission.

SRG/eROSITA is situated in a halo orbit around L2 where the highly variable solar wind charge exchange (SWCX) emission from Earth's magnetosheath is expected to be negligible. The soft X-ray foreground emissions from the local hot bubble (LHB) and the remaining heliospheric SWCX emissions could be studied in unprecedented detail with eROSITA All-Sky Survey (eRASS) data in a 6-month cadence and better spectral resolution than ROSAT. We aim to use eRASS data of the sight lines towards three giant molecular clouds away from the Galactic plane to isolate and study the soft X-ray diffuse foreground emission. These X-ray shadows will serve as calibration baselines for the future three-dimensional structural study of the LHB. We conducted spectral analysis on the diffuse X-ray spectra of these clouds from the first four eRASSs to estimate and separate the heliospheric SWCX contribution from the LHB emission. We find the density of the LHB to be independent of the sight line with $n_e \sim 4 \times 10^{-3}\,$cm$^{-3}$, but not the temperature. We report a lower temperature of $kT_{\mathrm{LHB}}=0.084\pm0.004\,$keV towards Chamaeleon$~$II & III (Cha$~$II & III) than Ophiuchus (Oph) and Corona Australis (CrA), in which we measured $0.102\pm0.006$ and $0.112\pm0.009\,$keV, respectively. We measured the emission measure of the LHB to be $\sim 2\times10^{-3}\,$cm$^{-6}\,$pc at medium Galactic latitudes ($|b| \sim 20^{\circ}$). A monotonic increase in the SWCX contribution has been observed since the start of 2020, coincidental with the beginning of solar cycle 25. For Oph, SWCX has dominated the LHB in the $0.3$-$0.7\,$keV band intensity since eRASS2. We observed lower SWCX contributions in Cha$~$II & III and CrA, consistent with the expected decreasing solar wind ion density at high heliographic latitudes.

In this paper, we review recent and ongoing work by our group on numerical simulations of relativistic jets. Relativistic outflows in Astrophysics are related to dilute, high energy plasmas, with physical conditions out of the reach of current laboratory capabilities. Simulations are thus imperative for the study of these objects. We present a number of such scenarios that have been studied by our group at the Universitat de Val\`encia. In particular, we have focused on the evolution of extragalactic outflows through galactic and intergalactic environments, deceleration by interaction with stars or clouds, or the propagation of jets in X-ray binaries and interaction with stellar winds from massive companions. All also share their role as particle acceleration sites and production of non-thermal radiation throughout the electromagnetic spectrum. Therefore, our work is not only aimed to understand the impact of outflows on their environments and thus their role in galaxy and cluster evolution, but also the nature and capabilities of these sites as generators of high and very-high energy radiation and cosmic rays.

We analyse the relation between the activity cycle length and the Rossby number and collected a sample of 44 main sequence stars with well-known activity cycle periods and rotation periods. We find a linear behaviour in the double-logarithmic relation between the Rossby number and cycle period. The bifurcation into a long and a short period branch is clearly real but it depends, empirically, on the colour index B-V, indicating a physical dependence on effective temperature and position on the main sequence. Furthermore, there is also a correlation between cycle length and convective turnover time with the relative depth of the convection zone. Based on this, we derive empirical relations between cycle period and Rossby number, and for the short period cycle branch relations, we estimate a scatter of the relative deviation between 14% and 28% on the long-period cycle branch. With these relations, we obtain a good match with the 10.3 yr period for the well known 11-year solar Schwabe cycle and a long-period branch value of 104 yr for the Gleissberg cycle of the Sun. Finally, we suggest that the cycles on the short-period branch appear to be generated in the deeper layers of the convective zone, while long-period branch cycles seem to be related to fewer deep layers in that zone. We show that for a broader B-V range, the Rossby number is a more suitable parameter for universal relation with cycle-rotation than just the rotation period alone. As proof, we demonstrate that our empirical stellar relations are consistent with the 11-year solar Schwabe cycle, in contrast to earlier studies using just the rotation period in their relations. Previous studies have tried to explain the cycle position of the Sun in the cycle-rotation presentation via other kinds of dynamo, however, in our study, no evidence is found that would suggest another type of dynamo for the Sun and other stars.

IceCube-Gen2, the next generation of the IceCube Neutrino Observatory at the South Pole, will consist of three co-located arrays: a deep Optical Array and a more shallow and larger Radio Array for neutrino detection in the ice, and a Surface Array above the footprint of the Optical Array. The Surface Array will be comprised of hybrid stations featuring elevated radio antennas and scintillation detectors, following the design of a prototype station successfully operating at the South Pole since 2020. Besides providing a veto for neutrino detection, the Surface Array will make IceCube-Gen2 a unique laboratory for cosmic-ray air showers. Compared to the current IceCube detector with its IceTop surface array, the aperture for coincident air-shower measurements detected by both, the deep optical and surface arrays, will increase by about a factor of 30. In addition to particle physics questions, such as the production of PeV muons and neutrinos in prompt decays, these surface-deep coincidences will be used to target astrophysical questions of the most energetic Galactic cosmic rays. The combination of particle and radio measurements at the surface and high-energy muons measured in the ice promises unprecedented accuracy for the mass composition in the energy range of the presumed Galactic-to-extragalactic transition - complementing the multimessenger science case of IceCube-Gen2. This proceeding provides an overview of the IceCube-Gen2 Surface Array and, in particular, its radio component.

Signal-to-noise ratios are a widely used concept for astroparticle radio detectors, such as air-shower radio arrays for cosmic-ray measurements or detectors searching for radio signals induced by neutrino interactions in ice. Nonetheless, no common standards or methods are established for the determination of the signal-to-noise ratio: values cannot be compared between experiments, and for the same signal and noise, various methods differ by large factors on the signal-to-noise ratio. This was the motivation to discuss a community-specific standardization at the ARENA conference 2022. No agreement on a common method to calculate signal-to-noise ratios was reached, however, awareness was raised that signal-to-noise ratios need to be well defined in any publications. This includes providing sufficient information on the procedure used to determine the signal-to-noise ratio, in addition to simply stating the formula. Even when using the same definition of the signal-to-noise ratio, there is still a significant dependence on the procedure of calculation, e.g., the signal-to-noise ratio of waveforms containing only background can vary significantly depending on the size of the time interval used as signal search window. To facilitate the interpretation of any signal-to-noise ratios in a specific study, the recommendation is to also state the mean value of the signal-to-noise ratio that the used method yields when applied to noise used in the study, e.g., the radio background measured by the corresponding experiment.

We report observations of the Didymos-Dimorphos binary asteroid system using the Atacama Large Millimeter/Submillimeter Array (ALMA) and the Atacama Compact Array (ACA) in support of the Double Asteroid Redirection Test (DART) mission. Our observations on UT 2022 September 15 provided a pre-impact baseline and the first measure of Didymos-Dimorphos' spectral emissivity at $\lambda=0.87$ mm, which was consistent with the handful of siliceous and carbonaceous asteroids measured at millimeter wavelengths. Our post-impact observations were conducted using four consecutive executions each of ALMA and the ACA spanning from T+3.52 to T+8.60 hours post-impact, sampling thermal emission from the asteroids and the impact ejecta. We scaled our pre-impact baseline measurement and subtracted it from the post-impact observations to isolate the flux density of mm-sized grains in the ejecta. Ejecta dust masses were calculated for a range of materials that may be representative of Dimorphos' S-type asteroid material. Depending on the material assumed, the average ejecta mass over our observations is consistent with 0.9--5.2$\times10^7$ kg, representing 0.2--1.2% of Dimorphos' total mass and in agreement with ejecta mass estimates based on measurements at optical wavelengths. Our results provide the most sensitive measure of mm-sized material in the ejecta and demonstrate the power of ALMA for providing supporting observations to spaceflight missions.

Bulge-disk (B-D) decomposition is an effective diagnostic to characterize the galaxy morphology and understand its evolution across time. So far, high-quality data have allowed detailed B-D decomposition to redshift below 0.5, with limited excursions over small volumes at higher redshifts. Next-generation large sky space surveys in optical, e.g. from the China Space Station Telescope (CSST), and near-infrared, e.g. from the space EUCLID mission, will produce a gigantic leap in these studies as they will provide deep, high-quality photometric images over more than 15000 deg2 of the sky, including billions of galaxies. Here, we extend the use of the Galaxy Light profile neural Network (GaLNet) to predict 2-S\'ersic model parameters, specifically from CSST data. We simulate point-spread function (PSF) convolved galaxies, with realistic B-D parameter distributions, on CSST mock observations to train the new GaLNet and predict the structural parameters (e.g. magnitude, effective radius, Sersic index, axis ratio, etc.) of both bulge and disk components. We find that the GaLNet can achieve very good accuracy for most of the B-D parameters down to an $r$-band magnitude of 23.5 and redshift $\sim$1. The best accuracy is obtained for magnitudes, implying accurate bulge-to-total (B/T) estimates. To further forecast the CSST performances, we also discuss the results of the 1-S\'ersic GaLNet and show that CSST half-depth data will allow us to derive accurate 1-component models up to $r\sim$24 and redshift z$\sim$1.7.

Super-Earth population, as one of the representatives of exoplanets, plays an important role in constraining the planet formation theories. According to the prediction from core-accretion models, super-Earths should be rare because their masses are in the range of the critical mass above which they would grow to be gas giants by runaway gas accretion. In this work, we investigate the effect of ohmic dissipation on the planetary thermal structure and cooling contraction as planets accrete gas from their surrounding disks. We find that the extra heating energy from Ohmic heating deposited into planetary envelopes can push the planetary radiative-convective boundaries inward and prevent the planets from cooling, and can even halt accretion. We explore parameter space to study the dependence of cooling timescale on the input parameters of the ohmic-dissipation model. Numerical results show that gas accretion can be halted before runaway gas accretion and the envelope mass is only several percent of planetary core mass for some parameter sets. Our results suggest that ohmic dissipation is a potential mechanism to delay the gas accretion and promote the formation of super-Earths. Future observations may help to constrain the importance of ohmic dissipation on the super-Earth formation.

We present a new library of semi-empirical stellar population models that are based on the empirical MILES and semi-empirical sMILES stellar libraries. The models span a large range of age and metallicity, in addition to an [$\alpha$/Fe] coverage from $-$0.2 to $+$0.6 dex, at MILES resolution (FWHM=2.5$ \mathring {\mathrm A}$) and wavelength coverage (3540.5-7409.6$ \mathring {\mathrm A}$). These models are aimed at exploring abundance ratios in the integrated light from stellar populations in star clusters and galaxies. Our approach is to build SSPs from semi-empirical stars at particular [$\alpha$/Fe] values, thus producing new SSPs at a range of [$\alpha$/Fe] values from sub-solar to super-solar. We compare these new SSPs with previously published and well-used models and find similar abundance pattern predictions, but with some differences in age indicators. We illustrate a potential application of our new SSPs, by fitting them to the high signal-to-noise data of stacked SDSS galaxy spectra. Age, metallicity and [$\alpha$/Fe] trends were measured for galaxy stacks with different stellar velocity dispersions and show systematic changes, in agreement with previous analyses of subsets of those data. These new SSPs are made publicly available.

LS I +61$^\circ$303 is a rare representative of the gamma-ray binaries with a compact object known to be a pulsar. We report on the periodicity and spectral analysis of this source performed with more than 14 years of Fermi/LAT data. The periodicity of LS I +61$^\circ$303 is strongly energy dependent. Two periods $P_1 = 26.932\pm 0.004 (stat)\pm 0.008 (syst)$ and $P_2 = 26.485 \pm 0.004 (stat)\pm 0.007 (syst)$ are detected only at $E>1$ GeV and at $E<0.3$ GeV correspondingly. Within $1\sigma$ (stat+syst) the periods are consistent with orbital ($P_2$) and beat orbital/superorbital ($P_1$) periods. We present the orbital light curves of the system in several energy bands and the results of the spectral analysis. We discuss the possible origin of the change in the variability pattern between 0.1 and 1 GeV energy.

Data analysis from upcoming large galaxy redshift surveys, such as Euclid and DESI will significantly improve constraints on cosmological parameters. To optimally extract the information from these galaxy surveys, it is important to control with a high level of confidence the uncertainty and bias arising from the estimation of the covariance that affects the inference of cosmological parameters. In this work, we are addressing two different but closely related issues: (i) the sampling noise present in a covariance matrix estimated from a finite set of simulations and (ii) the impact on cosmological constraints of the non-Gaussian contribution to the covariance matrix of the power spectrum. We focus on the parameter estimation obtained from fitting the matter power spectrum in real space, using the DEMNUni N-body simulations. Regarding the first issue, we adopt two different approaches to reduce the sampling noise in the precision matrix that propagates in the parameter space: on the one hand using an alternative estimator of the covariance matrix based on a non-linear shrinkage, NERCOME; and on the other hand employing a method of fast generation of approximate mock catalogs, COVMOS. We find that NERCOME can significantly reduce the noise induced on the posterior distribution of parameters, but at the cost of a systematic overestimation of the error bars on the cosmological parameters. We show that using a COVMOS covariance matrix estimated from a large number of realisations (10~000) results in unbiased cosmological constraints. Regarding the second issue, we quantify the impact on cosmological constraints of the non-Gaussian part of the power spectrum covariance purely coming from non-linear clustering. We find that when this term is neglected, both the errors and central values of the estimated parameters are affected for a scale cut $\kmax > 0.2\ \invMpc$.

We simulate the hydrodynamics of the wind flow in the B[e] supergiant binary system GG~Carinae and obtain the mass accretion rate onto the secondary and the observed lightcurve. We find an inhomogeneous Bondi-Hoyle-Lyttleton accretion into a curved accretion tail, and confirm that the accretion rate is modulated along the orbit, with a maximum close to periastron. We show that the accretion itself cannot account for the periodical variation in brightness. Instead, we explain the observed variation in the light curve with absorption by the accretion tail. Our results are in general agreement with previously derived stellar masses, orbital parameters, and the system orientation, but imply that the B[e] supergiant is more luminous. We find an effect related to the orbital motion of the two stars, in which the accretion tail is cut by the primary and we term it the Lizard Autotomy Effect. As part of the effect, the primary is self accreting wind that it ejected earlier. The Lizard Autotomy Effect creates an outwardly expanding spiral shell made up of broken segments. We suggest that such a tail exists in other B[e] supergiant systems and can be the source of the circumstellar material observed in such systems. The accretion also forms a disc around the secondary near periastron that later vanishes. We suggest that the formation of such a disc can launch jets that account for the bipolar structure observed around some B[e] supergiants.

Magnetic fields are extremely rare in close, hot binaries, with only 1.5\% of such systems known to contain a magnetic star. The eccentric $\epsilon$ Lupi system stands out in this population as the only close binary in which both stars are known to be magnetic. We report the discovery of strong, variable radio emission from $\epsilon$ Lupi using the upgraded Giant Metrewave Radio Telescope (uGMRT) and the MeerKAT radio telescope.The light curve exhibits striking, unique characteristics including sharp, high-amplitude pulses that repeat with the orbital period, with the brightest enhancement occurring near periastron. The characteristics of the light curve point to variable levels of magnetic reconnection throughout the orbital cycle, making $\epsilon$ Lupi the first known high-mass, main sequence binary embedded in an interacting magnetosphere. We also present a previously unreported enhancement in the X-ray light curve obtained from archival XMM-Newton data. The stability of the components' fossil magnetic fields, the firm characterization of their relatively simple configurations, and the short orbital period of the system make $\epsilon$ Lupi an ideal target to study the physics of magnetospheric interactions. This system may thus help us to illuminate the exotic plasma physics of other magnetically interacting systems such as moon-planet, planet-star, and star-star systems including T Tauri binaries, RS CVn systems, and neutron star binaries.

IC 10 as a starburst dwarf galaxy in the Local Group (LG) has a large population of newly formed stars that are massive and intrinsically very bright in comparison with other LG galaxies. Using the Isaac Newton Telescope (INT) with the Wide Field Camera (WFC) in the i-band and V-band, we performed an optical monitoring survey to identify the most evolved asymptotic giant branch stars (AGBs) and red supergiant stars (RSGs) in this star-forming galaxy, which can be used to determine the star formation history (SFH). The E(B - V) as an effective factor for obtaining the precise magnitude of stars is measured for each star using a 2D dust map (SFD98) to obtain a total extinction for each star in both the i-band and V-band. We obtained the photometric catalog for 53579 stars within the area of 0.07 deg$^{2}$ (13.5 kpc$^{2}$), of which 762 stars are classified as variable candidates after removing the foreground stars and saturated ones from our catalog. To reconstruct the SFH for IC 10, we first identified 424 long-period variable (LPV) candidates within the area of two half-light radii (2r$_{h}$) from the center of the galaxy. We estimated the recent star formation rate (SFR) at $\sim$ 0.32 M$_{\odot}$ yr$^{-1}$ for a constant metallicity Z = 0.0008, showing the galaxy is currently undergoing high levels of star formation. Also, a total stellar mass of 0.44 $\times$ 10$^{8}$ M$_{\odot}$ is obtained within 2r$_{h}$ for that metallicity.

We report the analysis of circumbinary discs formed in a radiation hydrodynamical simulation of star cluster formation. We consider both pure binary stars and pairs within triple and quadruple systems. The protostellar systems are all young (ages < $10^5$ yrs). We find that the systems that host a circumbinary disc have a median separation of $\approx 11$ au, and the median characteristic radius of the discs is $\approx 64$ au. We find that $89$ per cent of pure binaries with semi-major axes $a<1$ au have a circumbinary disc, and the occurrence rate of circumbinary discs is bimodal with log-separation in pure binaries with a second peak at $a \approx 50$ au. Systems with $a>100$ au almost never have a circumbinary disc. The median size of a circumbinary disc is between $\approx 5-6\ a$ depending on the order of the system, with higher order systems having larger discs relative to binary separation. We find the underlying distribution of mutual inclinations between circumbinary discs and binary orbit of both observed and simulated discs to not differ statistically.

Broadband radio waves emitted from pulsars are distorted and delayed as they propagate toward the Earth due to interactions with the free electrons that compose the interstellar medium, with lower radio frequencies being more impacted than higher frequencies. Multipath propagation in the interstellar medium results in both later times of arrival for the lower frequencies and causes the observed pulse to arrive with a broadened tail described via the pulse broadening function. We employ the CLEAN deconvolution technique to recover both the intrinsic pulse shape and pulse broadening function. This work expands upon previous descriptions of CLEAN deconvolution used in pulse broadening analyses by parameterizing the efficacy on simulated data and developing a suite of tests to establish which of a set of figures of merit lead to an automatic and consistent determination of the scattering timescale and its uncertainty. We compare our algorithm to simulations performed on cyclic spectroscopy estimates of the scattering timescale. We test our improved algorithm on the highly scattered millisecond pulsar J1903+0327, showing the scattering timescale to change over years, consistent with estimates of the refractive timescale of the pulsar.

A possible key element for large-scale energy release in the solar corona is an MHD kink instability in a single twisted magnetic flux tube. An initial helical current sheet fragments in a turbulent way into smaller-scale sheets, similarly to a nanoflare storm. As the loop expands in the radial direction during the relaxation process, an unstable loop can disrupt nearby stable loops and trigger an MHD avalanche. Exploratory investigations have been conducted in previous works with relatively simplified loop configurations. Here, we address a more realistic environment that comprehensively accounts for most of the physical effects involved in a stratified atmosphere, typical of an active region. The question is whether the avalanche process will be triggered, with what timescales, and how it will develop, as compared with the original, simpler approach. Three-dimensional MHD simulations describe the interaction of magnetic flux tubes, which have a stratified atmosphere, including chromospheric layers, the thin transition region to the corona, and the related transition from high-beta to low-beta regions. The model also includes the effects of thermal conduction and of optically thin radiation. Our simulations address the case where one flux tube among a few is twisted at the footpoints faster than its neighbours. We show that this flux tube becomes kink unstable first, in conditions in agreement with those predicted by analytical models. It rapidly involves nearby stable tubes, instigating significant magnetic reconnection and dissipation of energy as heat. The heating determines the development of chromospheric evaporation, while the temperature rises up to about 10 MK, close to microflares observations. This work confirms that avalanches are a viable mechanism for the storing and release of magnetic energy in plasma confined in closed coronal loops, as a result of photospheric motions.

The cross correlation between the CMB Doppler mode and the 21 cm line brightness temperature is calculated in the presence of a stochastic primordial magnetic field. Potential detectability is estimated for Planck 2018 bestfit parameters in combination with configuration and survey design parameters of 21 cm line radio telescopes such as LOFAR and the future SKAO. Homogeneous as well as inhomogeneous reionization has been considered. In particular the latter in combination with SKA1-mid shows promising signal-over-noise ratios.

Linear theory is a well developed framework for characterizing instabilities in weakly collisional plasmas, such as the solar wind. In the previous instalment of this series, we analyzed ~1.5M proton and alpha particle Velocity Distribution Functions (VDFs) observed by Helios I and II to determine the statistical properties of the standard instability parameters such as the growth rate, frequency, the direction of wave propagation, and the power emitted or absorbed by each component, as well as to characterize their behavior with respect to the distance from the Sun and collisional processing. In this work, we use this comprehensive set of instability calculations to train a Machine Learning algorithm consisting of three interlaced components that: 1) predict if an interval is unstable from observed VDF parameters; 2) predict the instability properties for a given unstable VDF; and 3) classify the type of the unstable mode. We use these methods to map the properties in multi-dimensional phase space to find that the parallel-propagating, proton-core-induced Ion Cyclotron mode dominates the young solar wind, while the oblique Fast Magnetosonic mode regulates the proton beam drift in the collisionally old plasma.

We present the densely sampled early light curve of the Type II supernova (SN) 2023ixf, first observed within hours of explosion in the nearby Pinwheel Galaxy (Messier 101; 6.7 Mpc). Comparing these data to recently updated models of shock cooling emission, we find that the progenitor likely had a radius of $410 \pm 10\ R_\odot$ (statistical uncertainty only), consistent with a red supergiant. These models provide a good fit to the data starting about 1 day after the explosion, despite the fact that the classification spectrum shows signatures of circumstellar material around SN 2023ixf during that time. Photometry during the first day after the explosion, provided almost entirely by amateur astronomers, does not agree with the shock cooling models or a simple power-law rise fit to data after 1 day. We consider the possible causes of this discrepancy, including precursor activity from the progenitor star, circumstellar interaction, and emission from the shock before or after it breaks out of the stellar surface. The very low luminosity ($-11\mathrm{\ mag} > M > -14\mathrm{\ mag}$) and short duration of the initial excess leads us to prefer a scenario related to prolonged emission from the SN shock traveling through the progenitor system.

Turbulence in protoplanetary disks, when present, plays a critical role in transporting dust particles embedded in the gaseous disk component. When using a field description of dust dynamics, a diffusion approach is traditionally used to model this turbulent dust transport. However, it has been shown that classical turbulent diffusion models are not fully self-consistent. Several shortcomings exist, including the ambiguous nature of the diffused quantity and the nonconservation of angular momentum. Orbital effects are also neglected without an explicit prescription. In response to these inconsistencies, we present a novel Eulerian turbulent dust transport model for isotropic and homogeneous turbulence on the basis of a mean-field theory. Our model is based on density-weighted averaging applied to the pressureless fluid equations and uses appropriate turbulence closures. Our model yields novel dynamic equations for the turbulent dust mass flux and recovers existing turbulent transport models in special limiting cases, thus providing a more general and self-consistent description of turbulent particle transport. Importantly, our model ensures the conservation of global angular and linear momentum unconditionally and implicitly accounts for the effects of orbital dynamics in protoplanetary disks. Furthermore, our model correctly describes the vertical settling-diffusion equilibrium solutions for both small and large particles. Hence, this work presents a generalized Eulerian turbulent dust transport model, establishing a comprehensive framework for more detailed studies of turbulent dust transport in protoplanetary disks.

After a brief review of the different approaches to predict the possible quantum gravity corrections to quantum field theory, we discuss in some detail the formulation based on a Gaussian reference frame fixing. Then, we implement this scenario to the determination of the inflationary spectrum of primordial perturbations. We consider the quantization of an inhomogeneous free massless scalar field on a quasi-classical isotropic Universe, developing a WKB expansion of the dynamics at the next order in the Planckian parameter, with respect to the one at which standard QFT emerges. The quantum gravity corrections to the scale invariant spectrum are discussed in a specific primordial cosmological setting and then in a general minisuperspace formalism, showing that there is no mode-dependent effect and thus the scale invariant inflationary spectrum is preserved. Such result is discussed in connection to the absence of a matter backreaction on the gravitational background in the considered paradigm.

We propose a novel framework where baryon asymmetry of the universe can arise due to forbidden decay of dark matter (DM) enabled by finite temperature effects in the vicinity of a first order phase transition (FOPT). In order to implement this novel cogenesis mechanism, we consider the extension of the standard model by one scalar doublet $\eta$, three right handed neutrinos (RHN), all odd under an unbroken $Z_2$ symmetry, popularly referred to as the scotogenic model of radiative neutrino mass. While the lightest RHN $N_1$ is the DM candidate and stable at zero temperature, there arises a temperature window prior to the nucleation temperature of the FOPT assisted by $\eta$, where $N_1$ can decay into $\eta$ and leptons generating a non-zero lepton asymmetry which gets converted into baryon asymmetry subsequently by sphalerons. The requirement of successful cogenesis forces the nucleation temperature to be lower than 34 GeV. This not only keeps the mass spectrum of new particles in specific ballpark but also leads to observable stochastic gravitational wave spectrum within the reach of planned experiments like LISA and BBO.

The effective number of relativistic neutrino species is a fundamental probe of the early Universe and its measurement represents a key constraint on many scenarios beyond the Standard Model of Particle Physics. In light of this, an accurate prediction of $N_{\rm eff}$ in the Standard Model is of pivotal importance. In this work, we consider the last ingredient needed to accurately calculate $N_{\rm eff}^{\rm SM}$: standard zero and finite temperature QED corrections to $e^+e^- \leftrightarrow \nu\bar{\nu}$ interaction rates during neutrino decoupling at temperatures around $T\sim {\rm MeV}$. We find that this effect leads to a reduction of $-0.0007$ in $N_{\rm eff}^{\rm SM}$. This NLO correction to the interaction rates, together with finite temperature QED corrections to the electromagnetic density of the plasma, and the effect of neutrino oscillations, implies that $N_{\rm eff}^{\rm SM} = 3.043$ with a theoretical uncertainty that is much smaller than any projected observational sensitivity.

Low-background liquid xenon detectors are utilized in the investigation of rare events, including dark matter and neutrinoless double beta decay. For their calibration, gaseous $^{220}$Rn can be used. After being introduced into the xenon, its progeny isotope $^{212}$Pb induces homogeneously distributed, low-energy ($<30$ keV) electronic recoil interactions. We report on the characterization of such a source for use in the XENONnT experiment. It consists of four commercially available $^{228}$Th sources with an activity of 55 kBq. These sources provide a high $^{220}$Rn emanation rate of about 9 kBq. We find no indication for the release of the long-lived $^{228}$Th above 1.7 mBq. Though an unexpected $^{222}$Rn emanation rate of about 3.6 mBq is observed, this source is still in line with the requirements for the XENONnT experiment.

We study with a 3D PIC simulation discontinuities between an electron-positron pair plasma and magnetized electrons and protons. A pair plasma is injected at one simulation boundary with a speed 0.6$c$ along its normal. It expands into an electron-proton plasma and a magnetic field that points orthogonally to the injection direction. Diamagnetic currents expel the magnetic field from within the pair plasma and pile it up in front of it. It pushes electrons, which induces an electric field pulse ahead of the magnetic one. This initial electromagnetic pulse (EMP) confines the pair plasma magnetically and accelerates protons electrically. The fast flow of the injected pair plasma across the protons behind the initial EMP triggers the filamentation instability. Some electrons and positrons cross the injection boundary and build up a second EMP. Electron-cyclotron drift instabilities perturb the plasma ahead of both EMPs seeding a Rayleigh-Taylor-type instability. Despite equally strong perturbations ahead of both EMPs, the second EMP is much more stable than the initial one. We attribute the rapid collapse of the initial EMP to the filamentation instability, which perturbed the plasma behind it. The Rayleigh-Taylor-type instability transforms the planar EMPs into transition layers, in which magnetic flux ropes and electrostatic forces due to uneven numbers of electrons and positrons slow down and compress the pair plasma and accelerate protons. In our simulation, the expansion speed of the pair cloud decreased by about an order of magnitude and its density increased by the same factor. Its small thickness implies that it is capable of separating a relativistic pair outflow from an electron-proton plasma, which is essential for collimating relativistic jets of pair plasma in collisionless astrophysical plasma.

Analysis and synthesis are key steps of the radio-interferometric imaging process, serving as a bridge between visibility and sky domains. They can be expressed as partial Fourier transforms involving a large number of non-uniform frequencies and spherically-constrained spatial coordinates. Due to the data non-uniformity, these partial Fourier transforms are computationally expensive and represent a serious bottleneck in the image reconstruction process. The W-gridding algorithm achieves log-linear complexity for both steps by applying a series of 2D non-uniform FFTs (NUFFT) to the data sliced along the so-called $w$ frequency coordinate. A major drawback of this method however is its restriction to direction-cosine meshes, which are fundamentally ill-suited for large field of views. This paper introduces the HVOX gridder, a novel algorithm for analysis/synthesis based on a 3D-NUFFT. Unlike W-gridding, the latter is compatible with arbitrary spherical meshes such as the popular HEALPix scheme for spherical data processing. The 3D-NUFFT allows one to optimally select the size of the inner FFTs, in particular the number of W-planes. This results in a better performing and auto-tuned algorithm, with controlled accuracy guarantees backed by strong results from approximation theory. To cope with the challenging scale of next-generation radio telescopes, we propose moreover a chunked evaluation strategy: by partitioning the visibility and sky domains, the 3D-NUFFT is decomposed into sub-problems which execute in parallel, while simultaneously cutting memory requirements. Our benchmarking results demonstrate the scalability of HVOX for both SKA and LOFAR, considering state-of-the-art challenging imaging setups. HVOX is moreover computationally competitive with W-gridder, despite the absence of domain-specific optimizations in our implementation.