Papers for Monday, Jan 16 2023

Binary neutron star mergers are one of the ultimate events of massive binary star evolution, and our understanding of their parent system is still in its infancy. Upcoming gravitational wave detections, coupled with multi-wavelength follow-up observations, will allow us to study an increasing number of these events by characterising their neighbouring stellar populations and searching for their progenitors. Stellar evolution simulations are essential to this work but they are also based on numerous assumptions. Additionally, the models used to study the host galaxies differ from those used to characterise the progenitors and are typically based on single star populations. Here we introduce a framework to perform an end-to-end analysis and deploy it to the first binary neutron star merger - GW170817. With the Binary Population And Spectral Synthesis (BPASS) codes we are able to retrieve the physical properties of the host galaxy NGC 4993 as well as infer progenitor candidates. In our simulations there is a >98% chance that GW170817 originated from a stellar population with Z=0.010 born between 5 and 12.5 Gyrs ago. By carefully weighing the stellar genealogies we find that GW170817 most likely came from a binary system born with a 13-24 Msol primary and 10-12 Msol secondary which underwent two or three common envelope events over their lifetime.

Gravitational wave (GW) standard sirens may resolve the Hubble tension, provided that standard siren inference of $H_0$ is free from systematic biases. However, standard sirens from binary neutron star (BNS) mergers suffer from two sources of systematic bias, one arising from the anisotropy of GW emission, and the other from the anisotropy of electromagnetic (EM) emission from the kilonova. For an observed sample of BNS mergers, the traditional Bayesian approach to debiasing involves the direct computation of the detection likelihood. This is infeasible for large samples of detected BNS merger due to the high dimensionality of the parameter space governing merger detection. In this study, we bypass this computation by fitting the Hubble constant to forward simulations of the observed GW and EM data under a simulation-based inference (SBI) framework using marginal neural ratio estimation. A key innovation of our method is the inclusion of BNS mergers which were only detected in GW, which allows for estimation of the bias introduced by EM anisotropy. Our method corrects for $\sim$90$\%$ of the bias in the inferred value of $H_0$ when telescope follow-up observations of BNS mergers have extensive tiling of the merger localization region, using known telescope sensitivities and assuming a model of kilonova emission. Our SBI-based method thus enables a debiased inference of the Hubble constant of BNS mergers, including both mergers with detected EM counterparts and those without.

It is routinely assumed that galaxy rotation curves are equal to their circular velocity curves (modulo some corrections) such that they are good dynamical mass tracers. We take a visualisation-driven approach to exploring the limits of the validity of this assumption for a sample of $33$ low-mass galaxies ($60<v_\mathrm{max}/\mathrm{km}\,\mathrm{s}^{-1}<120$) from the APOSTLE suite of cosmological hydrodynamical simulations. Only $3$ of these have rotation curves nearly equal to their circular velocity curves at $z=0$, the rest are undergoing a wide variety of dynamical perturbations. We use our visualisations to guide an assessment of how many galaxies are likely to be strongly perturbed by processes in several categories: mergers/interactions (affecting $6$/$33$ galaxies), bulk radial gas inflows ($19$/$33$), vertical gas outflows ($15$/$33$), distortions driven by a non-spherical DM halo ($17$/$33$), warps ($8$/$33$), and winds due to motion through the IGM ($5$/$33$). Most galaxies fall into more than one of these categories; only $5$/$33$ are not in any of them. The sum of these effects leads to an underestimation of the low-velocity slope of the baryonic Tully-Fisher relation ($\alpha\sim 3.1$ instead of $\alpha\sim 3.9$, where $M_\mathrm{bar}\propto v^\alpha$) that is difficult to avoid, and could plausibly be the source of a significant portion of the observed diversity in low-mass galaxy rotation curve shapes.

The distribution of mass in galaxy-scale strong gravitational lenses is often modelled as an elliptical power law plus `external shear', which notionally accounts for neighbouring galaxies and cosmic shear. We show that it does not. Except in a handful of rare systems, the best-fit values of external shear do not correlate with independent measurements of shear: from weak lensing in 45 Hubble Space Telescope images, or in 50 mock images of lenses with complex distributions of mass. Instead, the best-fit shear is aligned with the major or minor axis of 88% of lens galaxies; and the amplitude of the external shear increases if that galaxy is disky. We conclude that `external shear' attached to a power law model is not physically meaningful, but a fudge to compensate for lack of model complexity. Since it biases other model parameters that are interpreted as physically meaningful in several science analyses (e.g. measuring galaxy evolution, dark matter physics or cosmological parameters), we recommend that future studies of galaxy-scale strong lensing should employ more flexible mass models.

Simulations of molecular clouds often begin from highly idealised initial conditions, such as a uniform-density sphere with an artificially imposed turbulent velocity field. While the resulting structures may appear qualitatively similar to those detected in continuum and line observations, it is unclear whether they are genuinely representative of real molecular clouds. Recent observational work has discovered a tight, often close-to-linear relationship between the integrated intensity of molecular lines and the total column density of the cloud material. We combine magnetohydrodynamical simulations, time-dependent chemistry, and radiative transfer to produce synthetic molecular line observations of model clouds. We find similarly tight correlations between line intensity and column density to those observed, although the linear behaviour is only seen in isolated (as opposed to colliding) model clouds. This linear relationship is not due to optically thin emission; all lines investigated have high optical depths, and the increase in integrated intensity with column density is due to higher velocity dispersion along the line of sight. Overall, the idealised models commonly used in the literature appear to be reasonably accurate representations of real molecular clouds.

Nearby dwarf irregular galaxies are ideal laboratories for studying the interstellar medium (ISM) at low metallicity, which is expected to be common for galaxies at very high redshift that will be observed by the James Webb Space Telescope. We present the first high-resolution (~0.2 pc) hydrodynamical simulations of an isolated low-metallicity ($0.1~Z_\odot$) dwarf galaxy coupled with a time-dependent chemistry network and a dust evolution model where dust is locally produced and destroyed by various processes. To accurately model carbon monoxide (CO), we post-process the simulations with a detailed chemistry network including the time-dependent effect of molecular hydrogen (H$_2$). Our model successfully reproduces the observed star formation rate and CO(1-0) luminosity ($L_{\rm CO}$). We find that dust growth in dense gas is required to reproduce the observed $L_{\rm CO}$ as otherwise CO would be completely photodissociated. In contrast, the H$_2$ abundance is extremely small and is insensitive to dust growth, leading to a CO-to-H$_2$ conversion factor similar to the Milky Way value despite the low metallicity. Observationally inferred dust-to-gas ratio is thus underestimated if adopting the metallicity-dependent CO-to-H$_2$ conversion factor. The newly-produced dust in dense gas mixes with the ISM through supernova feedback without being completely destroyed by sputtering, which leads to galactic outflows 20% - 50% dustier than the ISM, providing a possible source for intergalactic dust.

The spectral line energy distribution of carbon monoxide contains information about the physical conditions of the star forming molecular hydrogen gas; however, the relation to local radiation field properties is poorly constrained. Using ~ 1-2 kpc scale ALMA observations of CO(3-2) and CO(4-3), we characterize the CO(4-3)/CO(3-2) line ratios of local analogues of main sequence galaxies at z ~ 1-2, drawn from the DYNAMO sample. We measure CO(4-3)/CO(3-2) across the disk of each galaxy and find a median line ratio of $R_{43} = 0.54^{+0.16}_{-0.15}$ for the sample. This is higher than literature estimates of local star-forming galaxies and is consistent with multiple lines of evidence that indicate DYNAMO galaxies, despite residing in the local Universe, resemble main-sequence galaxies at z ~ 1-2. Comparing to existing lower resolution CO(1-0) observations, we find $R_{41}$ and $R_{31}$ values in the range $\sim 0.2-0.3$ and $\sim 0.4-0.8$ respectively. We combine our kpc-scale resolved line ratio measurements with HST observations of H$\alpha$ to investigate the relation to star formation rate surface density and compare this relation to expectations from models. We find increasing CO(4-3)/CO(3-2) with increasing star formation rate surface density; however, models over-predict the line ratios across the range of star formation rate surface densities we probe, particularly at the lower range. Finally, SOFIA observations with HAWC+ and FIFI-LS reveal low dust temperatures and no deficit of [CII] emission with respect to the total infrared luminosity.

The Cassini spacecraft discovered that Saturn's moon Enceladus possesses a series of jets erupting from its South Polar Terrain. Previous studies of in situ data collected by Cassini's Ion and Neutral Mass Spectrometer (INMS) have identified H$_2$O, CO$_2$, CH$_4$, H$_2$, and NH$_3$ within the plume of ejected material. Identification of minor species in the plume remains an ongoing challenge, owing to the large number of possible combinations that can be used to fit the INMS data. Here, we present the discovery of several new compounds of strong importance to the habitability of Enceladus, including HCN, CH$_2$O, C$_2$H$_2$, and C$_3$H$_6$. Our analyses of the low velocity INMS data coupled with our detailed statistical framework enable discriminating between previously ambiguous species in the plume by alleviating the effects of high-dimensional model fitting. Together with plausible mineralogical catalysts and redox gradients derived from surface radiolysis, these compounds could potentially support extant microbial communities or drive complex organic synthesis leading to the origin of life.

Debris discs consist of belts of bodies ranging in size from dust grains to planetesimals; these belts are visible markers of planetary systems around other stars that can reveal the influence of extrasolar planets through their shape and structure. Two key stirring mechanisms -- self-stirring by planetesimals and secular perturbation by an external giant planet -- have been identified to explain the dynamics of planetesimal belts; their relative importance has been studied independently, but are yet to be considered in combination. In this work we perform a suite of 286 N-body simulations exploring the evolution of debris discs over 1~Gyr, combining the gravitational perturbations of both dwarf planets embedded in the discs, and an interior giant planet. Our systems were somewhat modeled after the architecture of the outer Solar system: a Solar mass star, a single massive giant planet at 30~au ($M_{\rm GP} =$ 10 to 316~$\mathrm{M}_{\oplus}$), and a debris disc formed by 100 massive dwarf planets and 1000 massless particles ($M_{\rm DD} =$ 3.16 to 31.6~$\mathrm{M}_{\oplus}$). We present the evolution of both the disc and the giant planet after 1~Gyr. The time evolution of the average eccentricity and inclination of the disc is strongly dependent on the giant planet mass as well as on the remaining disc mass. We also found that efficient stirring is achieved even with small disc masses. In general, we find that a mixed mechanism is more efficient in the stirring of cold debris discs than either mechanism acting in isolation.

Earlier, we suggested the "reload" concept of the polarimetric reverberation mapping of active galactic nuclei (AGN), proposed for the first time more than 10 years ago. We have successfully tested this approach of reverberation mapping of the broad emission line on the galaxy Mrk 6. It was shown that such an idea allows one to look at the AGN central parsec structure literally in a new light. However, the method originally assumed the use of spectropolarimetric observations, expensive in terms of telescope time, and implemented on rare large telescopes. Currently, we propose an adaptation of the polarimetric reverberation mapping of broad lines in medium-band filters following the idea of the photometric reverberation mapping, when filters are selected so that their bandwidth is oriented to the broad line and the surrounding continuum near. In this paper, we present the progress status of such monitoring conducted jointly at the Special astrophysical observatory and Asiago Cima Ekar observatory (OAPd/INAF) with support from Rozhen National Astronomical Observatory (NAO), some first results for the most frequently observed AGNs Mrk 335, Mrk 509, and Mrk 817, and the discussion of the future perspectives of the campaign.

NGC 4839 is the brightest galaxy (cD) of the NGC 4839 group at $R\approx 1$ Mpc in the south-west of the Coma cluster, which is known to be falling into Coma. However, it has been controversial whether it is in the first phase of infall or in the second phase of infall after passing the Coma center. We present a wide field study of globular clusters (GCs) in NGC 4839 and its environment based on Hyper Suprime-Cam $gr$ images in the Subaru archive. We compare the GC system of NGC 4839 with that of NGC 4816, which is the brightest member (S0) of the nearby group and lies at a similar distance in the west from the Coma center. Interestingly the spatial distribution of the GCs in NGC 4839 is significantly more compact than that of the GCs in NGC 4816. In addition, the radial number density profile of the GCs in NGC 4839 shows an abrupt drop at $R_{N4839}\approx 80$ kpc, while that of the GCs in NGC 4816 shows a continuous slow decline even in the outer region at $80<R_{N4816}<500$ kpc. The effective radius of the NGC 4839 GC system is about three times smaller than that of the NGC 4816 GC system. This striking difference can be explained if NGC 4839 lost a significant fraction of the GCs in its outskirt when it passed through Coma. This supports strongly the second infall scenario where the NGC 4839 passed the Coma center about 1.6 Gyr ago, and began the second infall after reaching the apocenter in the south-west recently.

We use the Upgraded Giant Metrewave Radio Telescope to measure scintillation arc properties in six bright canonical pulsars with simultaneous dual frequency coverage. These observations at frequencies from 300 to 750 MHz allowed for detailed analysis of arc evolution across frequency and epoch. We perform more robust determinations of arc curvature and scattering delay frequency-dependence than allowed by single-frequency-band-per-epoch measurements, which we find to agree with theory and previous literature. We report the discovery of a strong correlation between arc asymmetry and arc curvature, potentially indicating a link between scattering screen distance and refraction strength or the effect of asymmetric distribution of scattering material on a scattering screen. The inclusion of a 155 minute observation allowed us to resolve the scale of scintillation variations on short timescales, which we find to be directly tied to the amount of ISM sampled over the observation. Some of our pulsars showed either consistent or emerging asymmetries in arc curvature, indicating instances of refraction across their lines of sight. The presence of significant features in various pulsars, such as multiple scintillation arcs in PSR J1136+1551 and flat arclets in PSR J1509+5531, that have been found in previous works, were also sufficiently detected. Possible evidence for a timescale over which a given scattering screen dominates signal propagation was found by tracking visible scintillation arcs in each epoch in PSR J1136+1551. The interesting pulsar science accomplished with this upgraded telescope shows strong promise for important future work in pulsar astronomy.

Despite the success of galaxy-scale strong gravitational lens studies with Hubble-quality imaging, the number of well-studied strong lenses remains small. As a result, robust comparisons of the lens models to theoretical predictions are difficult. This motivates our application of automated Bayesian lens modeling methods to observations from public data releases of overlapping large ground-based imaging and spectroscopic surveys: Kilo-Degree Survey (KiDS) and Galaxy and Mass Assembly (GAMA), respectively. We use the open-source lens modeling software PyAutoLens to perform our analysis. We demonstrate the feasibility of strong lens modeling with large-survey data at lower resolution as a complementary avenue to studies that utilize more time-consuming and expensive observations of individual lenses at higher resolution. We discuss advantages and challenges, with special consideration given to determining background source redshifts from single-aperture spectra and to disentangling foreground lens and background source light. High uncertainties in the best-fit parameters for the models due to the limits of optical resolution in ground-based observatories and the small sample size can be improved with future study. We give broadly applicable recommendations for future efforts, and with proper application this approach could yield measurements in the quantities needed for robust statistical inference.

The new observational data concerning distribution, excitation, and kinematics of the ionized gas in the giant early-type disk galaxy NGC 2655 obtained at the 6m telescope of the Special Astrophysical Observatory (SAO RAS) and at the 2.5m telescope of the Caucasian Mountain Observatory of the Sternberg Astronomical Institute (CMO SAI MSU) are presented in this work. The joint analysis of these and earlier spectral observations has allowed us to make a conclusion about multiple nature of the gas in NGC 2655. Together with a proper large gaseous disk experiencing regular circular rotation in the equatorial plane of the stellar potential of the galaxy for billions years, we observe also remnants of a merged small satellite having striked the central part of NGC 2655 almost vertically for some 10 million years ago.

Hubble Space Telescope's (HST) Space Telescope Imaging Spectrograph (STIS) targeted 556 stars in a long-running program called Next Generation Spectral Library (NGSL) via proposals GO9088, GO9786, GO10222, and GO13776. Exposures through three low resolution gratings provide wavelength coverage from 0.2 $< \lambda <$ 1 $\mu$m at $\lambda/\Delta\lambda\sim$ 1000, providing unique coverage in the ultraviolet (UV). The UV grating (G230LB) scatters red light and this results in unwanted flux that becomes especially troubling for cool stars. We applied scattered light corrections based on \cite{2022stis.rept....5W} and flux corrections arising from pointing errors relative to the center of the 0\farcs2 slit. We present 514 fully reduced spectra, fluxed, dereddened, and cross-correlated to zero velocity. Because of the broad spectral range, we can simultaneously study H$\alpha$ and Mg II $\lambda$2800, indicators of chromospheric activity. Their behaviours are decoupled. Besides three cool dwarfs and one giant with mild flares in H$\alpha$, only Be stars show strong H$\alpha$ emission. Mg2800 emission, however, strongly anti-correlates with temperature such that warm stars show absorption and stars cooler than $5000 \: \! \rm{K}$ universally show chromospheric emission regardless of dwarf/giant status or metallicity. Transformed to Mg2800 flux emerging from the stellar surface, we find a correlation with temperature with approximately symmetric astrophysical scatter, in contrast to other workers who find a basal level with asymmetric scatter to strong values. Unsurprisingly, we confirm that Mg2800 activity is variable.

The bright A-type metallic-line star o Peg was reported in the early 1990s to have a surface magnetic field of ~2kG by analyzing the widths and strengths of spectral lines. In respect that those old studies were of rather empirical or approximate nature and the quality of observational data was not sufficient, this problem has been newly reinvestigated based on physically more rigorous simulations of line flux profiles, along with the observed equivalent widths (W) and full-widths at half-maximum (h) of 198 Fe I and 182 Fe II lines measured from the high-quality spectra. Given the Fe abundance derived from the conventional analysis, theoretical W and h values calculated for various sets of parameters were compared with the observed ones, which lead to the following conclusion regarding <H> (mean field strength). (1) An analysis of W yielded <H>~1-1.5kG from Fe II lines with the microturbulence of vt~1.5km/s. (2) A comparison of h resulted in <H>~1.5-2kG as well as the projected rotational velocity of vsini~5km/s. (3) Accordingly, the existence of mean magnetic field on the order of <H>~1-2kG in o Peg was confirmed, which is almost consistent with the consequence of the previous work.

Context. High latitude molecular clouds (hereafter HLMCs) permit the study of interstellar gas dynamics and astrochemistry with good accuracy due to their proximity, generally clear lines of sight, and lack of internal star-forming activity which can heavily modify the physical context. MBM 40, one of the nearest HLMCs, has been extensively studied, making it a superb target to infer and study the dust-to-gas mixing ratio (DGMR). Aims. The mixing of dust and gas in the interstellar medium remains a fundamental issue to keep track of astrochemistry evolution and molecular abundances. Accounting for both molecular and atomic gas is difficult because $H_2$ is not directly observable and HI spectra always show different dynamical profiles blended together which are not directly correlated with the cloud. We used two independent strategies to infer the molecular and atomic gas column densities and compute the dust-to-gas mixing ratio. Methods. We combined $HI$ 21 cm and $^{12}CO$ line observations with the IRAS 100 $\mu$m image to infer the dust-to-gas mixing ratio within the cloud. The cloud 21 cm profile was extracted using a hybrid Gaussian decomposition where $^{12}CO$ was used to deduce the total molecular hydrogen column density. Infrared images were used to calculate the dust emission. Results. The dust-to-gas mixing ratio is nearly uniform within the cloud as outlined by the hairpin structure. The total hydrogen column density and 100 $\mu$m emissivity are linearly correlated over a range in $N(H_{tot})$ of one order of magnitude.

We presented the multiwavelength analysis of a heavily obscured active galactic nucleus (AGN) in NGC 449. We first constructed a broadband X-ray spectrum using the latest NuSTAR and XMM-Newton data. Its column density ($\simeq 10^{24} \rm{cm}^{-2}$) and photon index ($\Gamma\simeq 2.4$) were reliably obtained by analyzing the broadband X-ray spectrum. However, the scattering fraction and the intrinsic X-ray luminosity could not be well constrained. Combined with the information obtained from the mid-infrared (mid-IR) spectrum and spectral energy distribution (SED) fitting, we derived its intrinsic X-ray luminosity ($\simeq 8.54\times 10^{42} \ \rm{erg\ s}^{-1}$) and scattering fraction ($f_{\rm{scat}}\simeq 0.26\%$). In addition, we also derived the following results: (1). The mass accretion rate of central AGN is about $2.54 \times 10^{-2} \rm{M}_\odot\ \rm{yr}^{-1}$, and the Eddington ratio is $8.39\times 10^{-2}$; (2). The torus of this AGN has a high gas-to-dust ratio ($N_{\rm H}/A_{\rm V}=8.40\times 10^{22}\ \rm{cm}^{-2}\ \rm{mag}^{-1}$); (3). The host galaxy and the central AGN are both in the early stage of co-evolution.

Millisecond magnetars produced in the center of dying massive stars are one prominent model to power gamma-ray bursts (GRBs). Their detailed nature, however, remains unsolved. To explore the effects of the initial mass, rotation, mass loss, and metallicity of the progenitor stars of $10-30~M_\odot$ on the formation and properties of the protomagnetars, we evolve over 150 single star models from the pre-main-sequence to core collapse by using the stellar evolution code MESA. We find that all of the fast-rotating stars become Wolf-Rayet stars. The final stellar, helium and carbon-oxygen core masses roughly increase with increasing initial mass, and decrease moderately with increasing initial rotation rate. We illustrate the effects of these intrinsic signatures on the hydrogen and helium envelopes and the metallicity. We then discuss the progenitors of the different types of supernovae. Furthermore, we find that the compactness parameter remains a nonmonotonic function of the initial mass and initial velocity when the effects of different metallicity and wind mass loss are considered. More importantly, we present the estimated period, magnetic field strength, and masses of protomagnetars in all cases. The typical rotational energy of these millisecond magnetars is sufficient to power long-duration GRBs.

Primordial non-Gaussianity has set strong constraints on models of the early universe. Studies have shown that Loop Quantum Cosmology (LQC), which is an attempt to extend inflationary scenario to planck scales, leads to a strongly scale dependent and oscillatory non-Gaussianity. In particular, the non-Gaussianity function $f_{_{\rm NL}}(k_1,\, k_2,\, k_3)$ generated in LQC, though similar to that generated during slow roll inflation at small scales, is highly scale dependent and oscillatory at large wavelengths. In this work, we investigate the imprints of such a primordial bispectrum in the bispectrum of Cosmic Microwave Background (CMB). Inspired by earlier works, we propose an analytical template for the primordial bispectrum in LQC and compute the corresponding reduced bispectra of temperature and electric polarisation and their three-point cross-correlations. We show that CMB bispectra generated in LQC is consistent with the observations from Planck. We conclude with a discussion of our results and its implications to LQC.

In this work, we introduce a quantitative methodology to define what is the main trunk and what are the significant branches of a minimum spanning tree (MST). We apply it to the pulsar tree, i.e. the MST of the pulsar population constructed upon a Euclidean distance over the pulsar's intrinsic properties. Our method makes use of the betweenness centrality estimator, as well as of non-parametric tests to establish the distinct character of the defined branches. Armed with these concepts, we study how the pulsar population has evolved throughout history, and analyze how to judge whether a new class of pulsars appears in new data, future surveys, or new incarnations of pulsar catalogs.

Multiple populations (MPs), characterized by variations in light elemental abundances, have been found in stellar clusters in the Milky Way, Magellanic Clouds, as well as several other dwarf galaxies. Based on a large amount of observations, mass has been suggested to be a key parameter affecting the presence and appearance of MPs in stellar clusters. To further investigate the existence of MPs in low-mass clusters and explore the mass threshold for MP formation, we carried out a project studying the stellar population composition in several low-mass Galactic globular clusters. Here we present our study on the cluster Eridanus. With blue-UV low-resolution spectra obtained with the OSIRIS/Multi-object spectrograph on the Gran Telescopio Canarias, we computed the spectral indices of CH and CN for the sample giant stars, and derived their carbon and nitrogen abundances using model spectra. A significant dispersion in the initial surface abundance of nitrogen was found in the sample, indicating the existence of MPs in Eridanus. Inspecting the age-initial mass distribution of in-situ clusters with MPs, we find a slight trend that initial mass increases with increasing age, and the lowest initial mass of log Minitial ~4.98 and 5.26 are found at the young and old end, respectively, which might provide a rough reference for the mass threshold for clusters to form MPs. However, more observations of clusters with low initial masses are still necessary before any firm conclusion can be drawn.

We present a study of the distribution of galaxies along the radius of 157 groups and clusters of galaxies (200~km~s$^{-1}$ < $\sigma$ < 1100~km~s$^{-1}$) of the local Universe (0.01 < $z$ < 0.1). We introduced a new boundary of galaxy systems and identified it with the splashback radius $R_{sp}$. We also identified the central region of galaxy systems with a radius of $R_c$. These radii are defined by the observed integrated distribution of the total number of galaxies depending on the squared distance from the center of the groups/clusters coinciding, as a rule, with the brightest galaxy. We show that the radius $R_{sp}$ is proportional to the $R_{200c}$ (radius of the virialized region of a galaxy cluster) and to the radius of the central region $R_c$ with a slope close to 1. Among the obtained dependences of the radii on X-ray luminosity, the $\log R_{sp}$ - $\log L_X$ relation has the lowest scatter. We measured $<R_{sp}>$ = $1.67\pm0.05$~Mpc for the total sample, $<R_{sp}>$ = $1.14\pm0.14$~Mpc for galaxy groups with $\sigma \leq$ 400~km~s$^{-1}$, $<R_{sp}>$ = $2.00\pm0.20$~Mpc for galaxy clusters with $\sigma$ > 400~km~s$^{-1}$. We found the average ratio of the radii $R_{sp}/R_{200c} = 1.40\pm0.02$ or $R_{sp}/R_{200m} = 0.88\pm0.02$.}

We consider a conventional $\alpha-\Omega$-dynamo model with meridional circulation that exhibits typical features of the solar dynamo, including a Hale cycle period of around 20 years and a reasonable shape of the butterfly diagram. With regard to recent ideas of a tidal synchronization of the solar cycle, we complement this model by an additional time-periodic $\alpha$-term that is localized in the tachocline region. It is shown that amplitudes of some dm/s are sufficient for this $\alpha$-term to become capable of entraining the underlying dynamo. We argue that such amplitudes of $\alpha$ may indeed be realistic, since velocities in the range of m/s are reachable, e.g., for tidally excited magneto-Rossby waves.

Interplanetary shocks are large-scale heliospheric structures often caused by eruptive phenomena at the Sun, and represent one of the main sources of energetic particles. Several interplanetary shock crossings by spacecraft at $1$ AU have revealed enhanced energetic-ion fluxes that extend far upstream of the shock. Surprisingly, in some shock events, ion fluxes with energies between $100$ keV and about $2$ MeV acquire similar values (which we refer to as ``overlapped'' fluxes), corresponding to flat energy spectra in that range. In contrast, closer to the shock, the fluxes are observed to depend on energy. In this work, we analyze three interplanetary shock-related energetic particle events observed by the Advanced Composition Explorer spacecraft where flat ion energy spectra were observed upstream of the shock. We interpret these observations via a velocity filter mechanism for particles in a given energy range. This reveals that low energy particles tend to be confined to the shock front and cannot easily propagate upstream, while high energy particles can. The velocity filter mechanism has been corroborated from observations of particle flux anisotropy by the Solid-State Telescope of Wind/3DP.

This document summarises the UK astronomy community's science and technology priorities for funding and investments in the coming decades, following a series of national community consultations by the Astronomy Advisory Panel of the Science and Technology Facilities Council (STFC). The facility remit of STFC is ground-based so the infrastructure recommendations are necessarily also ground-based, but the report also recognises the importance of STFC-funded technology development for, and science exploitation of, the ESA science program including but not limited to X-ray, gamma-ray and multimessenger astronomy.

We report on multi-wavelength observations of the tidal disruption event (TDE) candidate eRASSt J074426.3+291606 (J0744), located in the nucleus of a previously quiescent galaxy at $z=0.0396$. J0744 was first detected as a new, ultra-soft X-ray source (photon index $\sim 4$) during the second SRG/eROSITA All-Sky Survey (eRASS2), where it had brightened in the 0.3--2~keV band by a factor of more than $\sim$160 relative to an archival 3$\sigma$ upper limit inferred from a serendipitous Chandra pointing in 2011. The transient was also independently found in the optical by the Zwicky Transient Factory (ZTF), with the eRASS2 detection occurring only $\sim$20 days after the peak optical brightness, suggesting that the accretion disc formed promptly in this TDE. Continued X-ray monitoring over the following $\sim$400 days by eROSITA, NICER XTI and Swift XRT showed a net decline by a factor of $\sim$100, albeit with large amplitude X-ray variability where the system fades, and then rebrightens, in the 0.3--2~keV band by a factor $\sim$50 during an 80 day period. Contemporaneous Swift UVOT observations during this extreme X-ray variability reveal a relatively smooth decline, which persists over $\sim$400 days post-optical peak. The peak observed optical luminosity (absolute $g$-band magnitude $\sim -16.8$ mag) from this transient makes J0744 the faintest optically-detected TDE observed to date. However, contrasting the known set of `faint and fast' TDEs, the optical emission from J0744 decays slowly (exponential decay timescale $\sim$120~days), making J0744 the first member of a potential new class of `faint and slow' TDEs.

The ROSAT-selected tidal disruption event (TDE) candidate RX J133157.6-324319.7 (J1331), was detected in 1993 as a bright (0.2-2 keV flux of $(1.0 \pm 0.1) \times 10^{-12}$ erg s$^{-1}$ cm$^{-2}$), ultra-soft ($kT=0.11 \pm 0.03$ keV) X-ray flare from a quiescent galaxy ($z=0.05189$). During its fifth All-Sky survey (eRASS5) in 2022, SRG/eROSITA detected the repeated flaring of J1331, where it had rebrightened to an observed 0.2-2 keV flux of $(6.0 \pm 0.7) \times 10^{-13}$ erg s$^{-1}$ cm$^{-2}$, with spectral properties ($kT=0.115 \pm 0.007$ keV) consistent with the ROSAT-observed flare $\sim$30 years earlier. In this work, we report on X-ray, UV, optical, and radio observations of this system. During a pointed XMM observation $\sim$17 days after the eRASS5 detection, J1331 was not detected in the 0.2-2 keV band, constraining the 0.2-2 keV flux to have decayed by a factor of $\gtrsim$40 over this period. Given the extremely low probability ($\sim5\times 10^{-6}$) of observing two independent full TDEs from the same galaxy over a 30 year period, we consider the variability seen in J1331 to be likely caused by two partial TDEs involving a star on an elliptical orbit around a black hole. J1331-like flares show faster rise and decay timescales ($\mathcal{O}(\mathrm{days})$) compared to standard TDE candidates, with neglible ongoing accretion at late times post-disruption between outbursts.

The nucleosynthetic isotope dichotomy between carbonaceous and non-carbonaceous meteorites has been interpreted as evidence for spatial separation and coexistence of two distinct planet-forming reservoirs for several million years in the solar protoplanetary disk. Rapid formation of Jupiter's core within one million years after CAIs has been suggested as a potential mechanism for spatial and temporal separation. In this scenario, Jupiter's core would open a gap in the disk and trap inwards-drifting dust grains in the pressure bump at the outer edge of the gap, separating the inner and outer disk materials from each other. We performed simulations of dust particles in a protoplanetary disk with a gap opened by an early formed Jupiter core, including dust growth and fragmentation as well as dust transport using the dust evolution software DustPy. Our numerical experiments indicate that particles trapped in the outer edge of the gap rapidly fragment and are transported through the gap, contaminating the inner disk with outer disk materials on a timescale that is inconsistent with the meteoritic record. This suggests that other processes must have initiated or at least contributed to the isotopic separation between the inner and outer Solar System.

Neutron stars that show X-ray and $\gamma$-ray pulsed emission must, somewhere in the magnetosphere, generate electron-positron pairs. Such pairs are also required for radio emission, but then why do a number of these sources appear radio quiet? Here, we carried out a deep radio search towards four such neutron stars that are isolated X-ray/$\gamma$-ray pulsars but for which no radio pulsations have been detected yet. These sources are 1RXS J141256.0+792204 (Calvera), PSR J1958+2846, PSR J1932+1916 and SGR J1907+0919. Searching at lower radio frequencies, where the radio beam is thought to be wider, increases the chances of detecting these sources, compared to the earlier higher-frequency searches. We thus carried a search for periodic and single-pulse radio emission with the LOFAR radio telescope at 150 MHz. We used the known periods, and searched a wide range of dispersion measures, as the distances are not well constrained. We did not detect pulsed emission from any of the four sources. However, we put very constraining upper limits on the radio flux density at 150 MHz, of $\lesssim$1.4 mJy.

We present a smoothed density-corrected $V_{\rm max}$ technique for building a random catalog for property-dependent galaxy clustering estimation. This approach is essentially based on the density-corrected $V_{\rm max}$ method of Cole(2011), with three improvements to the original method. To validate the improved method, we generate two sets of flux-limited samples from two independent mock catalogs with different $k+e$ corrections. By comparing the two-point correlation functions, our results demonstrate that the random catalog created by the smoothed density-corrected $V_{\rm max}$ approach provides a more accurate and precise measurement for both sets of mock samples than the commonly used $V_{\rm max}$ method and redshift shuffled method. For flux-limited samples and color-dependent subsamples, the accuracy of the projected correlation function is well constrained within $1\%$ on the scale $0.07 h^{-1}\rm Mpc$ to $30 h^{-1}\rm Mpc$. The accuracy of the redshift-space correlation function is less than $2\%$ as well. Currently, it is the only approach that holds promise for achieving the high-accuracy goal of clustering measures for next-generation surveys.

We apply a population-orbit superposition method to 16 galaxies in the Fornax cluster observed with MUSE/VLT in the context of the Fornax3D project. By fitting the luminosity distribution, stellar kinematics, and age and metallicity maps simultaneously, we obtained the internal stellar orbit distribution, as well as the age and metallicity distribution of stars on different orbits for each galaxy. Based on the model, we decompose each galaxy into a dynamically cold disk (orbital circularity $\lambda_z\ge0.8$) and a dynamically hot non-disk component (orbital circularity $\lambda_z<0.8$), and obtain the surface-brightness, age, and metallicity radial profiles of each component. The galaxy infall time into the cluster is strongly correlated with galaxy cold-disk age with older cold disks in ancient infallers. We quantify the infall time $t_{\rm infall}$ of each galaxy with its cold-disk age using a correlation calibrated with TNG50 cosmological simulations. For galaxies in the Fornax cluster, we found that the luminosity fraction of cold disk in galaxies with $t_{\rm infall}>8$ Gyr are a factor of $\sim 4$ lower than in more recent infallers while controlling for total stellar mass. Nine of the 16 galaxies have spatially extended cold disks, and most of them show positive or zero age gradients; stars in the inner disk are $\sim 2-5$ Gyr younger than that in the outer disk, in contrast to the expectation of inside-out growth. Our results indicate that the assembly of cold disks in galaxies is strongly affected by their infall into clusters, by either removal of gas in outer regions or even tidally stripping or heating part of the pre-existing disks. Star formation in outer disks can stop quickly after the galaxy falls into the cluster, while star formation in the inner disks can last for a few Gyrs more, building the positive age gradient measured in cold disks.

Two major challenges of contemporary cosmology are the Hubble tension and the cosmic dipole tension. At the crossroad of these, we investigate the impact of peculiar velocities on estimations of the Hubble constant from time-delay cosmography. We quantify the bias on the inference of the Hubble constant due to peculiar velocities of the lens, the source and of the observer. The former two, which may cancel from one system to another, affect the determination of the angular diameter distances in the time-delay formula, and reconstructed quantities like the angle to the source, via a lens model. On the other hand, the peculiar velocity of the observer, which is a debated quantity in the context of the cosmic dipole tension, systematically affects observed angles through aberration, redshifts, angular diameter distance and reconstructed quantities. We compute in detail the effect of these peculiar velocities on the inference of the Hubble constant to linear order in the peculiar velocities for the seven lenses of the H0LiCOW/TDCOSMO collaboration. The bias generated by the observer's peculiar velocity alone can reach $1.15\%$ for the lenses which are well aligned with it. This results in a $0.25 \%$ bias for the seven combined lenses. Assuming a typical peculiar velocity of $300$ km s$^{-1}$ for the lens and the source galaxies, these add an additional random uncertainty, which can reach $1\%$ for an individual lens, but reduces to $0.24\%$ for the full TDCOSMO sample. The picture may change if peculiar velocities turn out to be larger than expected. Any time-delay cosmography program which aims for percent precision on the Hubble constant may need to take this source of systematic bias into account. This is especially so for future ground-based surveys which cover a fraction of the celestial sphere that is well aligned with the observer's peculiar velocity.

One of the sources of solar energetic particle (SEP) events is shocks that are driven by fast coronal mass ejections (CMEs). They can accelerate SEPs up to relativistic energies and are attributed to the largest SEP events. New studies suggest that CME-driven shocks can potentially accelerate electrons to MeV energies in the vicinity of the Sun. We focus on relativistic electrons associated with strong IP shocks between 2007 and 2019 to determine whether the shocks can keep accelerating such electrons up to 1 AU distance. We have analyzed High Energy Telescope (HET) observations aboard the STEREO spacecraft of potential electron energetic storm particle (ESP) events, characterized by intensity time series that peak at the time of, or close to, the associated CME-driven shock crossing. We present a new filtering method to assess the statistical significance of particle intensity increases and apply it to MeV electron observations in the vicinity of interplanetary shocks. We identified 27 candidate events by visual inspection from a STEREO in-situ shock list. Our method identified nine clear cases, where a significant increase of MeV electrons was found in association with a shock. Typically, the highest statistical significance was observed in the highest HET energy channel of electrons. All nine cases were associated with shocks driven by interplanetary CMEs that showed large transit speeds, in excess of 900 km/s. In several cases multiple shocks were observed within one day of the shock related to the electron increase. Although electron ESP events at MeV energies are found to be rare at 1 AU our filtering method is not designed to identify a potential interplanetary shock contribution from distances closer to the Sun. Future observations taken during closer approaches to the Sun will likely provide clarity on interplanetary shock acceleration of electrons.

Synthetic observations produced from radiative magnetohydronamic simulations have predicted that higher polarization fractions in the quiet solar photosphere would be revealed by increasing the total integration time of observations at GREGOR resolutions. We present recently acquired disk centre observations of the Fe I $15648.5$ $\mathrm{\AA}$ line obtained with the GREGOR telescope equipped with the GRIS-IFU during excellent seeing conditions, showing exceptionally high polarization fractions. Our observation reveal an internetwork region with a majority ($>60\%$) of magnetised pixels displaying a clear transverse component of the magnetic field. This result is in stark contrast to previous disk-centre GRIS-IFU observations in this spectral line, which had predominantly vertical magnetic fields in the deep photosphere. At the same time, the median magnetic field strength is weaker than previous GRIS-IFU observations, indicating that the larger fraction of polarization signals cannot be explained by a more active target. We use the Stokes Inversion based on Response functions (SIR) code to analyse the data, performing over $45$ million inversions, and interrogate the impact of two conflicting approaches to the treatment of noise on the retrieval of the magnetic inclination and azimuth. We present several case studies of the zoo of magnetic features present in these data, including small-scale magnetic loops that seem to be embedded in a sea of magnetism, and serpentine fields, focusing on regions where full-vector spectropolarimetry has been achieved. We also present a new open-source Python 3 analysis tool, SIR Explorer (SIRE), that we use to examine the dynamics of these small-scale magnetic features.

Context. The Extreme Ultraviolet Imager (EUI), onboard Solar Orbiter consists of three telescopes: the two High Resolution Imagers in EUV (HRIEUV) and in Lyman-{\alpha} (HRILya), and the Full Sun Imager (FSI). Solar Orbiter/EUI started its Nominal Mission Phase on 2021 November 27. Aims. EUI images from the largest scales in the extended corona off limb, down to the smallest features at the base of the corona and chromosphere. EUI is therefore a key instrument for the connection science that is at the heart of the Solar Orbiter mission science goals. Methods. The highest resolution on the Sun is achieved when Solar Orbiter passes through the perihelion part of its orbit. On 2022 March 26, Solar Orbiter reached for the first time a distance to the Sun close to 0.3 au. No other coronal EUV imager has been this close to the Sun. Results. We review the EUI data sets obtained during the period 2022 March-April, when Solar Orbiter quickly moved from alignment with the Earth (2022 March 6), to perihelion (2022 March 26), to quadrature with the Earth (2022 March 29). We highlight the first observational results in these unique data sets and we report on the in-flight instrument performance. Conclusions. EUI has obtained the highest resolution images ever of the solar corona in the quiet Sun and polar coronal holes. Several active regions were imaged at unprecedented cadences and sequence durations. We identify in this paper a broad range of features that require deeper studies. Both FSI and HRIEUV operate at design specifications but HRILya suffered from performance issues near perihelion. We conclude emphasising the EUI open data policy and encouraging further detailed analysis of the events highlighted in this paper.

The $VRC$ light curves were regularly measured for the eclipsing binary NSVS 7453183 as a part of our long-term observational project for studying of low-mass eclipsing binaries with a short orbital period and surface activity. The TESS light curve solution in Phoebe results to the detached configuration, where the temperature of primary component was adopted to $T_1$ = 4300 K according to the SED approximation. It gives us $T_2 =$ 4080 $\pm$ 100 K for the secondary component. The spectral type of the primary component was estimated to be K6 and the photometric mass ratio was derived $q = 0.86$. We confirm presence of the third body in this system, a stellar companion with a minimal mass 0.33 M$_{\rm Sun}$ orbiting the eclipsing pair with a short period about 425 days, and propose the next, fourth body with a longer orbiting period of about 12 years, probably a brown dwarf with the minimal mass of 50 M$_{\rm Jup}$. The hierarchical structure ((1+1)+1)+1 of this quadruple system is assumed. Characteristics and temporal variations of the dark region on the surface of the primary component were estimated. The average migration speed of about 10 deg/month was found during years 2020-2022.

The radio polarization properties of the pulsar population are only superficially captured by the conventional picture of pulsar radio emission. We study the broadband polarization of 271 young radio pulsars, focusing particularly on circular polarization, using high quality observations made with the Ultra-Wideband Low receiver on Murriyang, the Parkes radio telescope. We seek to encapsulate polarization behaviour on a population scale by defining broad categories for frequency- and phase-dependent polarization evolution, studying the co-occurrences of these categorizations and comparing them with average polarization measurements and spin-down energy ($\dot{E}$). This work shows that deviations of the linear polarization position angle (PA) from the rotating vector model (RVM) are linked to the presence of circular polarization features and to frequency evolution of the polarization. Polarization fraction, circular polarization contribution and profile complexity all evolve with $\dot{E}$ across the population, with the profiles of high-$\dot{E}$ pulsars being simple and highly linearly polarized. The relationship between polarization fraction and circular contribution is also seen to evolve such that highly polarized profiles show less variation in circular contribution with frequency than less strongly polarized profiles. This evolution is seen both across the population and across frequency for individual sources. Understanding pulsar radio polarization requires detailed study of individual sources and collective understanding of population-level trends. For the former, we provide visualizations of their phase- and frequency-resolved polarization parameters. For the latter, we have highlighted the importance of including the impact of circular polarization and of $\dot{E}$.

The major satellites of Jupiter and Saturn are believed to have formed in circumplanetary discs, which orbit forming giant protoplanets. Gas and dust in CPDs have different distributions and affect each other by drag, which varies with grain size. Yet simulations of multiple dust grain sizes with separate dynamics have not been done before. We seek to assess how much dust of each grain size there is in circumplanetary discs. We run multifluid 3D hydrodynamical simulations including gas and four discrete grain sizes of dust from 1$\mu$m to 1mm, representing a continuous distribution. We consider a 1 $M_\mathrm{Jup}$ protoplanet embedded in a protoplanetary disc around a 1 $M_{\odot}$ star. Our results show a truncated MRN distribution at smaller grain sizes, which starts to tail off by $a=100\mu$m and is near zero at 1mm. Large dust grains, which hold most of the dust mass, have very inefficient accretion to the CPD, due to dust filtration. Therefore CPDs' dust masses must be small, with mass ratio ~ a few $\times 10^{-6}$ to the protoplanet. These masses and the corresponding millimetre opacities are in line with CPD fluxes observed to date.

Euclid will survey most of the accessible extragalactic sky with imaging and slitless spectroscopy observations, creating a unique spectroscopic catalog of galaxies with H$\alpha$ line in emission that will map the Universe from $z=0.9$ to $1.8$. With low expected statistical errors, the error budget will likely be dominated by systematic errors related to uncertainties in the data and modelling. I will discuss the strategy that has been proposed to mitigate the expected systematic effects and propagate the uncertainty of mitigation to cosmological parameter errobars.

Questions regarding energy dissipation in astrophysical jets are open to date, despite of numerous attempts to limit the diversity of models. Some of the most popular models assume that energy is transferred to particles via internal shocks, which develop as a consequence of non-uniform velocity of the jet matter. In this context, we study the structure and energy deposition of colliding plasma shells, focusing our attention on the case of initially inhomogeneous shells. This leads to formation of distorted (corrugated) shock fronts -- a setup that has recently been shown to revive particle acceleration in relativistic magnetized perpendicular shocks. Our studies show that the radiative power of the far downstream of non-relativistic magnetized perpendicular shocks is moderately enhanced with respect to the flat shock cases. Based on the decay rate of downstream magnetic field, we make predictions for multiwavelength polarization properties.

The next-generation Event Horizon Telescope (ngEHT) will provide us with the best opportunity to investigate supermassive black holes (SMBHs) at the highest possible resolution and sensitivity. With respect to the existing Event Horizon Telescope (EHT) array, the ngEHT will provide increased sensitivity and uv-coverage with the addition of new stations, wider frequency coverage (from 86 GHz to 345 GHz and higher), finer resolution (<15 micro-arcseconds), and better monitoring capabilities. This will offer a unique opportunity to deeply investigate the physics around SMBHs, such as the disk-jet connection, the mechanisms responsible for high-energy photon and neutrino events, the role of magnetic fields in shaping relativistic jets, as well as the nature of binary SMBH systems. In this white paper we describe some ngEHT science cases in the context of multi-wavelength studies and synergies.

We present measurements of the CO luminosity functions (LFs) and the evolution of the cosmic molecular gas density out to z~6 based on an 8.5 arcmin^2 spectral scan survey at 3mm of the iconic Hubble Deep Field North (HDF-N) observed with the NOrthern Extended Millimeter Array (NOEMA). We use matched filtering to search for line emission from galaxies and determine their redshift probability distributions exploiting the extensive multi-wavelength data for the HDF-N. We identify the 7 highest-fidelity sources as CO emitters at 1<z<6, including the well-known submillimeter galaxy HDF850.1 at z=5.18. Four high-fidelity 3mm continuum sources are all found to be radio galaxies at z<=1, plus HDF850.1. We constrain the CO LFs in the HDF-N out to z~6, including a first measurement of the CO(5-4) LF at <z>=5.0. The relatively large area and depth of the NOEMA HDF-N survey extends the existing luminosity functions at 1<z<4 above the knee, yielding a somewhat lower density by 0.15-0.4 dex at the overlap region for the CO(2-1) and CO(3-2) transitions, attributed to cosmic variance. We perform a joint analysis of the CO LFs in the HDF-N and Hubble Ultra Deep Field from ASPECS, finding that they can be well described by a single Schechter function. The evolution of the cosmic molecular gas density from a joint analysis is in good agreement with earlier determinations. This implies that the impact of cosmic field-to-field variance on the measurements is consistent with previous estimates, adding to the challenges for simulations that model galaxies from first principles.

Earth and other rocky objects in the inner Solar System are depleted in carbon compared to objects in the outer Solar System, the Sun, or the ISM. It is believed that this is a result of the selective removal of refractory carbon from primordial circumstellar material. In this work, we study the irreversible release of carbon into the gaseous environment via photolysis and pyrolysis of refractory carbonaceous material during the disk phase of the early Solar System. We analytically solve the one-dimensional advection equation and derive an explicit expression that describes the depletion of carbonaceous material in solids under the influence of radial and vertical transport. We find both depletion mechanisms individually fail to reproduce Solar System abundances under typical conditions. While radial transport only marginally restricts photodecomposition, it is the inefficient vertical transport that limits carbon depletion under these conditions. We show explicitly that an increase in the vertical mixing efficiency, and/or an increase in the directly irradiated disk volume, favors carbon depletion. Thermal decomposition requires a hot inner disk (> 500 K) beyond 3 AU to deplete the formation region of Earth and chondrites. We find FU Ori-type outbursts to produce these conditions such that moderately refractory compounds are depleted. However, such outbursts likely do not deplete the most refractory carbonaceous compounds beyond the innermost disk region. Hence, the refractory carbon abundance at 1 AU typically does not reach terrestrial levels. Nevertheless, under specific conditions, we find photolysis and pyrolysis combined to reproduce Solar System abundances.

COSINE-200 is the next phase experiment of the ongoing COSINE-100 that aims to unambiguously verify the annual modulation signals observed by the DAMA experiment and to reach the world competitive sensitivity on the low-mass dark matter search. To achieve the physics goal of the COSINE-200, the successful production of the low-background NaI(Tl) detectors is crucial and it must begin from the mass production of the ultra-low background NaI powder. A clean facility for mass-producing the pure-NaI powder has been constructed at the Center for Underground Physics (CUP) in Korea. Two years of operation determined efficient parameters of the mass purification and provided a total of 480 kg of the ultra-pure NaI powder in hand. The potassium concentration in the produced powders varied from 5.4 to 11 ppb, and the maximum production capacity of 35 kg per two weeks was achieved. Here, we report our operational practice with the mass purification of the NaI powder, which includes raw powder purification, recycling of the mother solution, and recovery of NaI from the residual melt that remained after crystal growth.

In the past 15 years, the number of known planets outside of our solar system has grown from about 200 to more than 5000. During that time, we have conducted one of the longest crowdfunding campaigns in history, a nonprofit adopt a star program that supports astronomy research. The program includes the targets of NASA space telescopes that are searching for planets around other stars, and it uses the proceeds to help determine the properties of those stars and their planetary systems. I summarize how this innovative program has evolved over the years and engaged the public worldwide to support an international team of astronomers.

We discuss the damping of inflationary gravitational waves (GW) that re-enter the horizon before or during an epoch, where the energy budget of the universe is dominated by an unstable right handed neutrino (RHN), whose out of equilibrium decay releases entropy. Starting from the minimal Standard Model extension, motivated by the observed neutrino mass scale, with nothing more than 3 RHN for the Seesaw mechanism, we discuss the conditions for high scale leptogenesis assuming a thermal initial population of RHN. We further address the associated production of potentially light non-thermal dark matter and a potential component of dark radiation from the same RHN decay. One of our main findings is that the frequency, above which the damping of the tensor modes is potentially observable, is completely determined by successful leptogenesis and a Davidson-Ibarra type bound to be at around $0.1\;\text{Hz}$. To quantify the detection prospects of this GW background for various proposed interferometers such as AEDGE, BBO, DECIGO, Einstein Telescope or LISA we compute the signal-to-noise ratio (SNR). This allows us to investigate the viable parameter space of our model, spanned by the mass of the decaying RHN $M_1 \gtrsim 2.4\times 10^8\;\text{GeV} \cdot \sqrt{2\times 10^{-7}\;\text{eV}/\tilde{m}_1}$ (for leptogenesis) and the effective neutrino mass parameterizing its decay width $\tilde{m}_1< 2.9\times 10^{-7}\;\text{eV}$ (for RHN matter domination). Thus gravitational wave astronomy is a novel way to probe both the Seesaw and the leptogenesis scale, which are completely inaccessible to laboratory experiments in high scale scenarios.