Papers for Thursday, Apr 07 2022

We explore the impact of diffusive cosmic rays (CRs) on the evolution of the interstellar medium (ISM) under varying assumptions of supernova explosion environment. In practice, we systematically vary the relative fractions of supernovae (SN) occurring in star-forming high-density gas and those occurring in random locations decoupled from star-forming gas to account for SN from run-away stars or explosions in regions that have been cleared by prior SN, stellar winds, or radiation. We explore various mixed models by adjusting these fractions relative to each other. We find that in the simple system of a periodic stratified gas layer the ISM structure will evolve to one of two solutions: a "peak driving" state where warm gas is volume filling or a "thermal runaway" state where hot gas is volume filling. CR pressure and transport are important factors that strongly influence the solution state the ISM reaches and have the ability to flip the ISM between solutions. Observable signatures such as gamma ray emission and HI gas are explored. We find that gamma ray luminosity from pion decay is largely consistent with observations for a range of model parameters. The thickness of the HI gas layer may be too compact, however, this may be due to a large cold neutral fraction of midplane gas. The volume fraction of hot gas evolves to stable states in both solutions, but neither settles to a Milky Way-like configuration, suggesting that additional physics which is omitted here (e.g. a cosmological circum-galactic medium, radiation transport, or spectrally resolved and spatially varying CR transport) may be required.

An extensive catalog of spatially-resolved galaxy rotation curves and multi-band optical light profiles for 1752 observed spiral galaxies is assembled to explore the drivers of diversity in galaxy structural parameters, rotation curve shapes, and stellar mass profiles. Similar data were extracted from the NIHAO galaxy simulations to identify any differences between observations and simulations. Several parameters, including the inner slope "S" of a rotation curve (RC), were tested for diversity. Two distinct populations are found in observed and simulated galaxies; (i) blue, low mass spirals with stellar mass M* < 10^9.3 Msol and roughly constant "S", and (ii) redder, more massive and more diverse spirals with rapidly increasing "S". In all cases, the value of "S" seems equally contributed by the baryonic and non-baryonic (dark) matter. Diversity is shown to increase mildly with mass. Numerical simulations reproduce well most baryon-dominated galaxy parameter distributions, such as the inner stellar mass profile slope and baryonic scaling relations, but they struggle to match the full diversity of observed galaxy rotation curves (through "S") and most dark-matter-dominated parameters. To reproduce observations, the error broadening of the simulation's intrinsic spread of RC metrics would have to be tripled. The differences in various projections of observed and simulated scaling relations may reflect limitations of current sub-grid physics models to fully capture the complex nature of galaxies. For instance, AGNs are shown to have a significant effect on the shapes of simulated RCs. The inclusion of AGN feedback brings simulated and observed inner RC shapes into closer agreement.

We study the relationship between the morphology and star formation history (SFH) of 361 quiescent galaxies (QGs) at redshift $\langle z_{obs}\rangle\approx 2$, with stellar mass $\log M_*\ge10.3$, selected with the UVJ technique. Taking advantage of panchromatic photometry covering the rest-frame UV-to-NIR spectral range ($\approx40$ bands), we reconstruct the non-parametric SFH of the galaxies with the fully Bayesian SED fitting code Prospector. We find that the half-light radius $R_e$, observed at $z_{obs}$, depends on the formation redshift of the galaxies, $z_{form}$, and that this relationship depends on stellar mass. At $\log M_*<11$, the relationship is consistent with $R_e\propto(1+z_{form})^{-1}$, in line with the expectation that the galaxies' central density depends on the cosmic density at the time of their formation, i.e. the "progenitor effect". At $\log M_*>11$, the relationship between $R_e$ and $z_{form}$ flattens, suggesting that mergers become increasingly important for the size growth of more massive galaxies after they quenched. We also find that the relationship between $z_{form}$ and galaxy compactness similarly depends on stellar mass. While no clear trend is observed for QGs with $\log M_*>11$, lower-mass QGs that formed earlier, i.e. with larger $z_{form}$, have larger central stellar mass surface densities, both within the $R_e$ ($\Sigma_e$) and central 1 kpc ($\Sigma_{1kpc}$), and also larger $M_{1kpc}/M_*$, the fractional mass within the central 1 kpc. These trends between $z_{form}$ and compactness, however, essentially disappear, if the progenitor effect is removed by normalizing the stellar density with the cosmic density at $z_{form}$. Our findings highlight the importance of reconstructing the SFH of galaxies before attempting to infer their intrinsic structural evolution.

Outflows play a central role in galaxy evolution shaping the properties of galaxies. Understanding outflows and their effects in low luminosity AGNs, such as LINERs, is essential (e.g. they are a numerous AGN population in the local Universe). We obtained VLT/MUSE and GTC/MEGARA optical IFS-data for NGC1052, the prototypical LINER. The stars are distributed in a dynamically hot disc, with a centrally peaked velocity dispersion map and large observed velocity amplitudes. The ionised gas, probed by the primary component is detected up to $\sim$30arcsec ($\sim$3.3 kpc) mostly in the polar direction with blue and red velocities ($\mid$V$\mid$$<$250 km/s). The velocity dispersion map shows a notable enhancement ($\sigma$$>$90 km/s) crossing the galaxy along the major axis of rotation in the central 10arcsec. The secondary component has a bipolar morphology, velocity dispersion larger than 150 km/s and velocities up to 660 km/s. A third component is detected but not spatially resolved. The maps of the NaD absorption indicate optically thick neutral gas with a velocity field consistent with a slow rotating disc ($\Delta$V = 77$\pm$12 km/s) but the velocity dispersion map is off-centred without any counterpart in the flux map. We found evidence of an ionised gas outflow with mass of 1.6$\pm$0.6 $\times$ 10$^{5}$ Msun, and mass rate of 0.4$\pm$0.2 Msun/yr. The outflow is propagating in a cocoon of gas with enhanced turbulence and might be triggering the onset of kpc-scale buoyant bubbles (polar emission). Taking into account the energy and kinetic power of the outflow (1.3$\pm$0.9 $\times$ 10$^{53}$ erg and 8.8$\pm$3.5 $\times$ 10$^{40}$ erg/s, respectively) as well as its alignment with both the jet and the cocoon, and that the gas is collisionally ionised, we consider that the outflow is jet-powered, although some contribution from the AGN is possible.

Halo models and halo occupation distributions (HODs) are important tools to model the galaxy and matter distribution. We present and assess a new method for constraining the parameters of HODs using the gravitational lensing shear around galaxy pairs, galaxy-galaxy-galaxy-lensing (G3L). In contrast to galaxy-galaxy-lensing, G3L is sensitive to correlations between the per-halo numbers of galaxies from different populations. We use G3L to probe these correlations and test the default hypothesis that they are negligible. We derive a halo model for G3L and validate it with realistic mock data from the Millennium Simulation and a semi-analytic galaxy model. Then, we analyse public data from the Kilo-Degree Survey (KiDS), the VISTA Infrared Kilo-Degree Galaxy Survey (VIKING) and data from the Galaxy And Mass Assembly Survey (GAMA) to infer the HODs of galaxies at $z<0.5$ in five different stellar mass bins between $10^{8.5}h^{-2} M_\odot$ and $10^{11.5}h^{-2} M_\odot$ and two colours (red and blue), as well as correlations between satellite numbers. The analysis recovers the true HODs in the simulated data within the $68\%$ credibility range. The inferred HODs vary significantly with colour and stellar mass. There is also strong evidence ($>3\sigma$) for correlations, increasing with halo mass, between the numbers of red and blue satellites and galaxies with stellar masses below $10^{10} \Msun. Possible causes of these correlations are the selection of similar galaxies in different samples, the survey flux limit, or physical mechanisms like a fixed ratio between the satellite numbers of distinct populations. The decorrelation for halos with smaller masses is probably an effect of shot noise by low-occupancy halos. The inferred HODs can be used to complement galaxy-galaxy-lensing or galaxy clustering HOD studies or as input to cosmological analyses and improved mock galaxy catalogues.

We present FORS2@VLT follow-up photometry of YMCA-1, a recently discovered stellar system located 13\degr~from the Large Magellanic Cloud (LMC) center. The deep color-magnitude diagram (CMD) reveals a well-defined main sequence (MS) and a handful of stars in the post-MS evolutionary phases. We analyse the YMCA-1 CMD by means of the automated isochrone matching package {\tt ASteCA} and model its radial density profile with a Plummer function. We find that YMCA-1 is an old ($11.7^{+1.7}_{-1.3}$~Gyr), metal-intermediate ([Fe/H] $\simeq -1.12^{+0.21}_{-0.13}$~dex), compact (r$_{\rm h} = 3.5 \pm 0.5$ pc), low-mass (M $= 10^{2.45 \pm 0.02} M_{\odot}$) and low-luminosity (M$_V = -0.47 \pm 0.57$~mag) stellar system. The estimated distance modulus ($\mu_0 = 18.72^{+0.15}_{-0.17}$~mag), corresponding to about 55~kpc, suggests that YMCA-1 is associated to the LMC, but we cannot discard the scenario in which it is a Milky Way satellite. The structural parameters of YMCA-1 are remarkably different compared with those of the 15 known old LMC globular clusters. In particular, it resides in a transition region of the M$_V$-r$_h$ plane, in between the ultra-faint dwarf galaxies and the classical old clusters, and close to SMASH-1, another faint stellar system recently discovered in the LMC surroundings.

Analysis of interstellar absorption lines observed in high-resolution {\em HST} spectra of nearby stars provides temperatures, turbulent velocities, and kinetic properties of warm interstellar clouds. Previous studies identified 15 warm partially ionized clouds within about 10~pc of the Sun and measured their mean thermal and kinematic properties. A new analysis of 100 interstellar velocity components reveals a wide range of temperatures and turbulent velocities within the Local Interstellar Cloud (LIC) and other nearby clouds. These variations appear to be random with Gaussian distributions. We find no trends of these properties with stellar distance, angle from the Galactic Center, the main source of EUV radiation (the star $\epsilon$~CMa), the center of the LIC, or the direction of inflowing interstellar matter into the heliosphere. The spatial scale for temperature variations in the LIC is likely smaller than 5,100~au, a distance that the Sun will traverse in 1,000 years. Essentially all velocity components align with known warm clouds. We find that within 4~pc of the Sun, space is completely filled with partially ionized clouds, but at larger distances space is only partially filled with partially ionized clouds, indicating that fully ionized inter-cloud gas fills the voids. We find that the neutral hydrogen number density in the LIC and likely other warm clouds in the CLIC is about 0.10~cm$^{-3}$ rather than the 0.20~cm$^{-3}$ density that may be representative of only the immediate environment of the LIC. The 3,000--12,000~K temperature range for the gas is consistent with the predictions of theoretical models of the WNM and WIM, but the high degree of inhomogeneity within clouds argues against simple theoretical models. Finally, we find evidence for a shock in the sight line to the star AD~Leo.

The primary difficulty in understanding the sources and processes that powered cosmic reionization is that it is not possible to directly probe the ionizing Lyman Continuum (LyC) radiation at that epoch as those photons have been absorbed by the intervening neutral hydrogen in the IGM on their way to us. It is therefore imperative to build a model to accurately predict LyC emission using other properties of galaxies in the reionization era. In recent years, studies have shown that the LyC emission from galaxies may be correlated to their Lya emission. Here, we study this correlation by analyzing thousands of galaxies at high-z in the SPHINX cosmological simulation. We post-process these galaxies with the Lya radiative transfer code RASCAS and analyze the Lya - LyC connection. We find that the Lya and LyC luminosities are strongly correlated with each other, although with dispersion. There is a positive correlation between Lya and LyC escape fractions in the brightest Lya emitters (>$10^{41}$ erg/s), similar to the recent observational studies. However, when we also include fainter Lya emitters (LAEs), the correlation disappears, which suggests that the observed relationship may be driven by selection effects. We also find that bright LAEs are dominant contributors to reionization ($> 10^{40}$ erg/s galaxies contribute $> 90\%$ of LyC emission). Finally, we build predictive models using multivariate linear regression where we use the physical and the Lya properties of simulated galaxies to predict their intrinsic and escaping LyC luminosities with a high degree of accuracy. We find that the most important galaxy properties to predict the escaping LyC luminosity of a galaxy are its escaping Lya luminosity, gas mass, gas metallicity, and SFR. These models can be very useful to predict LyC emissions from galaxies and can help us identify the sources of reionization.

We present two new GW Vir-type pulsating white dwarf stars, TIC\,0403800675 (WD\,J115727.68-280349.64) and TIC\,1989122424 (WD J211738.38-552801.18) discovered in the Transiting Exoplanet Survey Satellite (TESS) photometric data. For both stars, the TESS light curves reveal the presence of oscillations with periods in a narrow range between 400 and 410\,s, which are associated with typical gravity ($g$)-modes. Follow-up ground-based spectroscopy shows that both stars have similar effective temperature ($T_\mathrm{eff} = 110,000 \pm 10,000$\,K) and surface gravity ($\log g = 7.5 \pm 0.5$), but different He/C composition (mass fractions): He\,=\,0.75 and C\,=\,0.25 for TIC\,0403800675, and He\,=\,0.50 and C\,=\,0.50 for TIC\,1989122424. By performing a fit to their spectral energy distributions, we found for both stars radii and luminosities of $R=0.019\pm0.002\,R_\odot$ and $\log(L/L_\odot)=1.68^{+0.15}_{-0.24}$, respectively. By employing evolutionary tracks of PG~1159 stars, we find the masses of both stars to be $0.56\pm0.18 M_{\odot}$ from the $\log g$-$T_\mathrm{eff}$ diagram and $0.60^{+0.11}_{-0.09} M_{\odot}$ from the Hertzsprung Russell diagram.

How molecular clouds fragment and create the dense structures which go on to form stars is an open question. We investigate the relative importance of different energy terms (kinetic, thermal, magnetic, and gravity - both self-gravity and tidal forces) for the formation and evolution of molecular clouds and their sub-structures based on the SILCC-Zoom simulations. These simulations follow the self-consistent formation of cold molecular clouds down to scales of 0.1 pc from the diffuse supernova-driven interstellar medium in a stratified galactic disc. We study the time evolution of seven molecular clouds (five with magnetic fields and two without) for 1.5-2 Myr. Using a dendrogram, we identify hierarchical 3D sub-structures inside the clouds with the aim to understand their dynamics and distinguish between the theories of gravo-turbulent fragmentation and global hierarchical collapse. The virial analysis shows that the dense gas is indeed dominated by the interplay of gravity and turbulence, while magnetic fields and thermal pressure are only important for fluffy, atomic structures. Over time, gravitationally bound sub-structures emerge from a marginally bound medium (viral ratio $1 \leq \alpha_{\rm vir}^{\rm vol} <2$) as a result of large-scale supernova-driven inflows rather than global collapse. A detailed tidal analysis shows that the tidal tensor is highly anisotropic. Yet the tidal forces are generally not strong enough to disrupt either large-scale or dense sub-structures but cause their deformation. By comparing tidal and crossing time scales, we find that tidal forces do not seem to be the main driver of turbulence within the molecular clouds.

In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exciting prospects for exploiting OHMs to probe the cosmic history of gas-rich mergers. Using observations for the Looking At the Distant Universe with the MeerKAT Array (LADUMA) deep HI survey, we report the first untargeted detection of an OHM at $z > 0.5$, LADUMA J033046.20$-$275518.1 (nicknamed "Nkalakatha"). The host system, WISEA J033046.26$-$275518.3, is an infrared-luminous radio galaxy whose optical redshift $z \approx 0.52$ confirms the MeerKAT emission line detection as OH at a redshift $z_{\rm OH} = 0.5225 \pm 0.0001$ rather than HI at lower redshift. The detected spectral line has 18.4$\sigma$ peak significance, a width of $459 \pm 59\,{\rm km\,s^{-1}}$, and an integrated luminosity of $(6.31 \pm 0.18\,{\rm [statistical]}\,\pm 0.31\,{\rm [systematic]}) \times 10^3\,L_\odot$, placing it among the most luminous OHMs known. The galaxy's far-infrared luminosity $L_{\rm FIR} = (1.576 \pm 0.013) \times 10^{12}\,L_\odot$ marks it as an ultra-luminous infrared galaxy; its ratio of OH and infrared luminosities is similar to those for lower-redshift OHMs. A comparison between optical and OH redshifts offers a slight indication of an OH outflow. This detection represents the first step towards a systematic exploitation of OHMs as a tracer of galaxy growth at high redshifts.

In this work we report the discovery and analysis of six new compact triply eclipsing triple star systems found with the TESS mission: TICs 37743815, 42565581, 54060695, 178010808, 242132789, and 456194776. All of these exhibit distinct third body eclipses where the inner eclipsing binary (EB) occults the third (`tertiary') star, or vice versa. We utilized the TESS photometry, archival photometric data, and available archival spectral energy distribution curves (SED) to solve for the properties of all three stars, as well as many of the orbital elements. We describe in detail our SED fits, search of the archival data for the outer orbital period, and the final global photodynamical analyses. From these analyses we find that all six systems are coplanar to within $0^\circ$ - $5^\circ$, and are viewed nearly edge on (i.e., within a couple of degrees). The outer orbital periods and eccentricities of the six systems are {$P_{\rm out}$ (days), $e$}: {68.7, 0.36}, {123, 0.16}, {60.7, 0.01}, {69.0, 0.29}, {41.5, 0.01}, {93.9, 0.29}, respectively, in the order the sources are listed above. The masses of all 12 EB stars were in the range of 0.7-1.8 M$_\odot$ and were situated near the main sequence. By contrast, the masses and radii of the tertiary stars ranged from 1.5-2.3 M$_\odot$ and 2.9-12 R$_\odot$, respectively. We use this information to estimate the occurrence rate of compact flat triple systems.

Context. The close approach of the near-Earth asteroid (99942) Apophis to Earth in 2029 will provide a unique opportunity to examine how the physical properties of the asteroid could be changed due to the Earth's gravitational perturbation. As a result, the Republic of Korea is planning a rendezvous mission to Apophis. Aims. Our aim was to use photometric data from the apparitions in 2020-2021 to refine the shape model and spin state of Apophis. Methods. Using thirty-six 1 to 2-m class ground-based telescopes and the Transiting Exoplanet Survey Satellite, we performed a photometric observation campaign throughout the 2020-2021 apparition. The convex shape model and spin state were refined using the light-curve inversion method. Results. According to our best-fit model, Apophis is rotating in a short axis mode with rotation and precession periods of 264.178 hours and 27.38547 hours, respectively. The angular momentum vector orientation of Apophis was found as (275$^\circ$, -85$^\circ$) in the ecliptic coordinate system. The ratio of the dynamic moments of inertia of this asteroid was fitted to $I_a:I_b:I_c=0.64:0.97:1$, which corresponds to an elongated prolate ellipsoid. These findings regarding the spin state and shape model could be used to not only design the space mission scenario but also investigate the impact of the Earth's tidal force during close encounters.

Detecting the cosmic 21-cm signal during the Epoch of Reionisation and Cosmic Dawn will reveal insights into the properties of the first galaxies and advance cosmological parameter estimation. Until recently, the primary focus for astrophysical parameter inference from the 21-cm signal centred on the power spectrum (PS). However, the cosmic 21-cm signal is highly non-Gaussian rendering the PS sub-optimal for characterising the cosmic signal. In this work, we introduce a new technique to analyse the non-Gaussian information in images of the 21-cm signal called the Wavelet Scattering Transform (WST). This approach closely mirrors that of convolutional neural networks with the added advantage of not requiring tuning or training of a neural network. Instead, it compresses the 2D spatial information into a set of coefficients making it easier to interpret while also providing a robust statistical description of the non-Gaussian information contained in the cosmic 21-cm signal. First, we explore the application of the WST to mock 21-cm images to gain valuable physical insights by comparing to the known behaviour from the 21-cm PS. Then we quantitatively explore the WST applied to the 21-cm signal by extracting astrophysical parameter constraints using Fisher Matrices from a realistic 1000 hr mock observation with the Square Kilometre Array. We find that: (i) the WST applied only to 2D images can outperform the 3D spherically averaged 21-cm PS, (ii) the excision of foreground contaminated modes can degrade the constraining power by a factor of ~1.5-2 with the WST and (iii) higher cadences between the 21-cm images can further improve the constraining power.

We studied Romanov V48, which had been considered as a possible polar below the period minimum (orbital period of 0.0420991 d). We analyzed the publicly available Zwicky Transient Facility data and found that this object is an intermediate polar with an orbital period of 0.102312 d (in the period gap) and a spin period of 0.0420991 d. The amplitude of the spin variation was very large (0.6 mag) and the profile has been confirmed to vary with the beat period of 0.071534 d. The ratio between the spin and orbital periods was large and the system resembles the intermediate polar DW Cnc below the period gap. Infrared emission from Romanov V48, however, could not be explained by radiation from the secondary. The infrared excess and the shape of the spectral energy distribution resembled those of polars, suggesting that the emission mechanism in Romanov V48 is similar to those of polars. Romanov V48 may be an intermediate object between intermediate polars and polars.

Understanding the origin of merging binary black holes is currently one of the most pressing quests in astrophysics. We show that if isolated binary evolution dominates the formation mechanism of merging binary black holes, one should expect a correlation between the effective spin parameter, $\chi_\mathrm{eff}$, and the redshift of the merger, $z$, of binary black holes. This correlation comes from tidal spin-up systems preferentially forming and merging at higher redshifts due to the combination of weaker orbital expansion from low metallicity stars given their reduced wind mass loss rate, delayed expansion and have smaller maximal radii during the supergiant phase compared to stars at higher metallicity. As a result, these tightly bound systems merge with short inspiral times. Given our fiducial model of isolated binary evolution, we show that the origin of a $\chi_\mathrm{eff}-z$ correlation in the detectable LIGO--Virgo binary black hole population is different from the intrinsic population, which will become accessible only in the future by third-generation gravitational-wave detectors such as Einstein Telescope and Cosmic Explorer. Finally, we compare our model predictions with population predictions based on the current catalog of binary black hole mergers and find that current data favor a positive correlation of $\chi_\mathrm{eff}-z$ as predicted by our model of isolated binary evolution.

In the context of the IRAM 30m Large Program GEMS, we present a study of thioformaldehyde in several starless cores located in star-forming filaments of Taurus, Perseus, and Orion. We investigate the influence of the environmental conditions on the abundances of these molecules in the cores, and the effect of time evolution. We have modelled the observed lines of H2CS, HDCS, and D2CS using the radiative transfer code RADEX. We have also used the chemical code Nautilus to model the evolution of these species depending on the characteristics of the starless cores. We derive column densities and abundances for all the cores. We also derive deuterium fractionation ratios, Dfrac, to determine and compare the evolutionary stage between different parts of each star-forming region. Our results indicate that the north region of the B213 filament in Taurus is more evolved than the south, while the north-eastern part of Perseus presents an earlier evolutionary stage than the south-western zone. Model results also show that Dfrac decreases with the cosmic-ray ionisation rate, while it increases with density and with the degree of sulphur depletion. In particular, we only reproduce the observations when the initial sulphur abundance in the starless cores is at least one order of magnitude lower than the solar elemental sulphur abundance. The progressive increase in HDCS/H2CS and D2CS/H2CS with time makes these ratios powerful tools for deriving the chemical evolutionary stage of starless cores. However, they cannot be used to derive the temperature of these regions, since both ratios present a similar evolution at two different temperature ranges (7-11 K and 15-19 K). Regarding chemistry, (deuterated) thioformaldehyde is mainly formed through gas-phase reactions (double-replacement and neutral-neutral displacement reactions), while surface chemistry plays an important role as a destruction mechanism.

Cassiopeia A (Cas A) supernova remnant shows strong radiation from radio to gamma-ray bands. The mechanism of gamma-ray radiation in Cas A and its possible contribution to PeV cosmic rays are still under debate. The X-ray imaging reveals an asymmetric profile of Cas A, suggesting the existence of a jet-like structure. In this work, we propose an asymmetrical model for Cas A, consisting of a fast moving jet-like structure and a slowly expanding isotropic shell. This model can account for the multi-wavelength spectra of Cas A, especially for the power-law hard X-ray spectrum from $\sim$ 60 to 220 keV. The GeV to TeV emission from Cas A should be contributed by both hadronic and leptonic processes. Moreover, the jet-like structure may produce a gamma-ray flux of $\sim 10^{-13}\rm erg\ cm^{-2}\ s^{-1}$ at $\sim 100$ TeV, to be examined by LHAASO and CTA.

The specific angular momenta ($j\equiv J/M$) of stars ($j_{\star}$), gas ($j_{\mathrm{gas}}$), baryons as a whole ($j_{\mathrm{b}}$) and dark matter haloes ($j_{\mathrm{h}}$) contain clues of vital importance about how galaxies form and evolve. Using one of the largest samples of disc galaxies (S0-BCD) with high-quality rotation curves and near-infrared surface photometry, we perform a detailed comparative analysis of $j$ that stretches across a variety of galaxy properties. Our analysis imposes tight constraints on the "retained" fractions of specific angular momentum ($j_{\star}/j_{\mathrm{h}}$, $j_{\mathrm{HI}}/j_{\mathrm{h}}$ and $j_{\mathrm{b}}/j_{\mathrm{h}}$), as well as on their systematic trends with mass fraction and galaxy morphology, thus on how well specific angular momentum is conserved in the process of disc galaxy formation and evolution. Besides, our analysis demonstrates how challenging it is to characterize barred galaxies from a gravitational instability point of view. This is true not only for the popular Efstathiou, Lake & Negroponte (1982) bar instability criterion, which fails to separate barred from non-barred galaxies in about 55% of the cases, but also for the mass-weighted Toomre (1964) parameter of atomic gas, $\langle Q_{\mathrm{HI}}\rangle$, which succeeds in separating barred from non-barred galaxies, but only in a statistical sense.

X-ray images of galaxy clusters often show disturbed structures that are indication of cluster mergers. As a complementary to our previous work on dynamical state of 964 clusters observed by the Chandra, we process the X-ray images for 1308 clusters from XMM-Newton archival data, together with the images of 22 clusters newly released by the Chandra, and evaluate their dynamical state from these X-ray images. The concentration index $c$, the centroid shift $\omega$, the power ratio $P_3/P_0$ are calculated in circular regions with a certain radius of 500 kpc, and the morphology index $\delta$ is estimated within elliptical regions adaptive to the cluster size and shape. In addition, the dynamical parameters for 42 clusters previously estimated from Chandra images are upgraded based on the newly available redshifts. Good consistence is found between dynamical parameters derived from XMM-Newton and Chandra images for the overlapped sample of clusters in the two datasets. The dependence of mass scaling relations on dynamical state is shown by using data of 388 clusters.

LAGO, the Latin American Giant Observatory, is an extended cosmic ray observatory, consisting of a wide network of water Cherenkov detectors located in 10 countries. With different altitudes and geomagnetic rigidity cutoffs, their geographic distribution, combined with the new electronics for control, atmospheric sensing and data acquisition, allows the realisation of diverse astrophysics studies at a regional scale. It is an observatory designed, built and operated by the LAGO Collaboration, a non-centralised alliance of 30 institutions from 11 countries. While LAGO has access to different computational frameworks, it lacks standardised computational mechanisms to fully grasp its cooperative approach. The European Commission is fostering initiatives aligned to LAGO objectives, especially to enable Open Science and its long-term sustainability. This work introduces the adaptation of LAGO to this paradigm within the EOSC-Synergy project, focusing on the simulations of the expected astrophysical signatures at detectors deployed at the LAGO sites around the World.

We report high-precision observations of the linear polarization of the F1$\,$III star $\theta$ Scorpii. The polarization has a wavelength dependence of the form expected for a rapid rotator, but with an amplitude several times larger than seen in otherwise similar main-sequence stars. This confirms the expectation that lower-gravity stars should have stronger rotational-polarization signatures as a consequence of the density dependence of the ratio of scattering to absorption opacities. By modelling the polarization, together with additional observational constraints (incorporating a revised analysis of Hipparcos astrometry, which clarifies the system's binary status), we determine a set of precise stellar parameters, including a rotation rate $\omega\, (= \Omega/\Omega_{\rm c})\ge 0.94$, polar gravity $\log{g_p} = 2.091 ^{+0.042}_{-0.039}$ (dex cgs), mass $3.10 ^{+0.37}_{-0.32}$ solar masses, and luminosity $\log(L/Lsun) =3.149^{+0.041}_{-0.028}$. These values are incompatible with evolutionary models of single rotating stars, with the star rotating too rapidly for its evolutionary stage, and being undermassive for its luminosity. We conclude that $\theta$ Sco A is most probably the product of a binary merger.

Faraday tomography of radio polarimetric data below 200 MHz from LOFAR are providing us with a new perspective on the diffuse and magnetized interstellar medium (ISM). Of particular interest is the discovery of Faraday-rotated synchrotron polarization associated with neutral gas, as traced by atomic hydrogen (HI) and dust. Here we present the first in-depth numerical study of these LOFAR results. We produce and analyze comprehensive synthetic observations of low-frequency synchrotron polarization from MHD simulations of colliding super shells in the multiphase ISM. We define five distinct gas phases over more than four orders of magnitude in gas temperature and density, ranging from hot, and warm, fully ionized gas to cold neutral medium. We focus on the contribution of each gas phase to synthetic observations of both rotation measure and synchrotron polarized intensity below 200 MHz. We also investigate the link between the latter and synthetic observations of optically thin HI gas. We find that, not only the fully ionized gas but also the warm partially ionized and neutral phases strongly contribute to the total rotation measure and polarized intensity. However, the contribution of each phase to the observables strongly depends on the choice of integration axis and the orientation of the mean magnetic field with respect to the shell-collision axis. Strong correlation between HI synthetic data and synchrotron polarized intensity, reminiscent of LOFAR results, is obtained with lines of sight perpendicular to the mean magnetic field direction. Our study suggests that multiphase modelling of MHD processes is needed in order to interpret observations of the radio sky at low frequency. This work is a first step toward understanding the complexity of low-frequency synchrotron emission that will be soon revolutionized by large-scale surveys with LOFAR and the SKA.

In the neutral hydrogen (HI) intensity mapping (IM) survey, the foreground contamination on the cosmological signals is extremely severe, and the systematic effects caused by radio telescopes themselves further aggravate the difficulties in subtracting foreground. In this work, we investigate whether the deep learning method, concretely the 3D U-Net algorithm here, can play a crucial role in foreground subtraction when considering the systematic effect caused by the telescope primary beam. We consider two beam models, i.e., the Gaussian beam model as a simple case and the Cosine beam model as a sophisticated case. To make a comparison, and also a combination, with the U-Net method, we also employ the traditional principal component analysis (PCA) method in the foreground subtraction. We find that, in the case of the Gaussian beam, both the PCA and U-Net methods can effectively clean the foreground, but in the case of the Cosine beam, U-Net performs much better than PCA in cleaning the foreground. In order to show how well the PCA and U-Net methods can recover the HI signals, we also derive the angular power spectra, as well as the 2D power spectrum of HI after performing the foreground subtractions. It is found that, in the case of Gaussian beam, the concordance with the original HI map using U-Net is better than that using PCA by $27.4\%$, and in the case of Cosine beam, the concordance using U-Net is better than that using PCA by $144.7\%$. Therefore, the U-Net based foreground subtraction can efficiently eliminate the telescope primary beam effect and shed new light on recovering the HI power spectrum for future HI IM experiments.

Relativistic models of light propagation adopted for high-precision astrometry are based on the parametrised post-Newtonian formalism, which provides a framework for examining the effects of a hypothetical violation of general relativity on astrometric data. Astrometric observations are strongly affected by the post-Newtonian parameter $\gamma$ describing the strength of gravitational light deflection. We study both analytically and numerically how a deviation in the PPN parameter $\gamma$ from unity, which is the value predicted by general relativity, affects the parallax estimations in Gaia-like astrometry. Changes in the observable quantities produced by a small variation in PPN $\gamma$ were calculated analytically. We then considered how such variations of the observables are reflected in the parallax estimations, and we performed numerical simulations to check the theoretical predictions. A variation in the PPN $\gamma$ results in a global shift of parallaxes and we present a formula describing the parallax bias in terms of the satellite barycentric distance, the angle between the spin axis and the direction to the Sun, and the PPN $\gamma$ uncertainty. Numerical simulations of the astrometric solutions confirm the theoretical result. The up-to-date estimation of PPN $\gamma$ suggests that a corresponding contribution to the Gaia parallax zero point unlikely exceeds 0.2 $\mu$as. The numerical simulations indicate that the parallax shift is strongly dependent on ecliptic latitude. It is argued that this effect is due to an asymmetry in the Gaia scanning law and this conclusion is fully validated by additional simulations with a reversed direction of the precession of the spin axis around the direction to the Sun.

We present observations of the carbon-rich protoplanetary nebula (PPN) CRL 2688 made with the Institut de Radioastronomie Millimetrique (IRAM) 30 m telescope in the 3mm and 2mm bands. In total, 196 transition lines belonging to 38 molecular species and isotopologues are detected, among which, to our best knowledge, 153 transition lines and 13 species are the first report for this object. Additionally, in order to contribute to future research, we have collected observational data on the molecular lines of CRL 2688 from the literature and compiled them into a single unified catalog. We find that the molecular abundance of CRL 2688 cannot be explained by the standard model of a circumstellar envelope. The implications of metal-bearing molecules on circumstellar chemistry are discussed.

Long secondary periods (LSPs), observed in a third of pulsating red giant and supergiant stars, are the only unexplained type of large-amplitude stellar variability known at this time. Numerous authors have explored various scenarios for the origin of LSPs, but were unable to give a final solution to this problem. We present known properties of LSP variables and show new results proving that the physical mechanism responsible for LSPs is binarity. Namely, the LSP light changes are due to the presence of a dusty cloud orbiting the red giant together with a brown-dwarf companion and obscuring the star once per orbit. In this scenario, the low-mass companion is a former planet that accreted a significant amount of mass from the envelope of its host star and grew into a brown dwarf.

Stochastic gravitational waves can be induced from the primordial curvature perturbations generated during inflation, through scalar-tensor mode coupling at the second order of cosmological perturbation theory. Here we discuss a model of warm inflation in which large curvature perturbations are generated at the small scales because of inflaton dissipation. These overdense perturbations then collapse at later epoch to form primordial black holes, as was studied in our earlier work (Ref. \cite{Arya:2019wck}), and therefore may also act as a source to the second-order tensor perturbations. In this study, we calculate the spectrum of these secondary gravitational waves from our warm inflationary model. We find that our model leads to a generation of scalar induced gravitational waves (SIGW) over a frequency range ($1-10^6$) Hz. Further, we discuss the detection possibilities of these SIGWs, taking in account the sensitivity of different ongoing and future gravitational wave experiments.

We simulate the collapse of a turbulent gas cloud with the particle-based code Gadget2. We choose two sub-clouds, formed by those particles located around the centers $\vec{r}_L$ and $\vec{r}_R$ and within the radii $r_L$ and $r_R$, respectively. A translational velocity $\vec{v}_L$ or $\vec{v}_R$ is added, so that the sub-clouds move towards each other to collide. The radius and pre-collision velocity of the sub-clouds are chosen to be un-equal, and both head-on and oblique collisions are considered. The simulations are all calibrated to have the same total mass and the initial energy ratio $\alpha=0.16$, which is defined as the ratio of thermal energy to gravitational energy. We compare low-$\beta$ models to a high-$\beta$ models, where $\beta$ is defined as the ratio of kinetic energy to gravitational energy. Finally, we consider the turbulent cloud to be under the gravitational influence of an object located far enough, in order to approximate the tidal effects by means of an azimuthal velocity $V_{\rm cir}$ added to the cloud particles in addition to the translational and turbulent velocities mentioned above. We compare a low-$V_{\rm cir}$ model with a high-$V_{\rm cir}$ one.

We investigate the kinematics of the molecular gas in a sample of seven edge-on (i>60 degrees) galaxies identified as hosting large-scale outflows of ionised gas, using ALMA CO(1-0) observations at ~1 kpc resolution. We build on Hogarth et al. (2021), where we find that molecular gas is more centrally concentrated in galaxies which host winds than in control objects. We perform full 3-dimensional kinematic modelling with multiple combinations of kinematic components, allowing us to infer whether these objects share any similarities in their molecular gas structure. We use modelling to pinpoint the kinematic centre of each galaxy, in order to interpret their minor- and major-axis position velocity diagrams (PVDs). From the PVDs, we find that the bulk of the molecular gas in our galaxies is dynamically cold, tracing the rotation of the stellar disc, but with minor flux asymmetries. Most notably, we find evidence of radial gas motion in some of our minor-axis PVDs. In our highest S/N object, we include bi-symmetric radial flow in our kinematic model, and find (via the Bayesian Information Criterion) that the presence of radial gas motion is strongly favoured. This may provide a mechanism by which molecular gas and star formation are centrally concentrated, enabling the launch of massive ionised gas winds.

The spectrograph ESPRESSO obtained updated limits on the variation of the fine-structure constant, $\alpha$, through measurements along the line of sight of a bright quasar with unprecedented precision and accuracy. These impose new constraints on cosmological models with a varying $\alpha$. We assume such a model where the electromagnetic sector is coupled to a scalar field dark energy responsible for the current acceleration of the Universe. We parametrise the variation of $\alpha$ with two extra parameters, one defining the cosmological evolution of the quintessence component and the other fixing the coupling with the electromagnetic field. The objective of this work is to constrain these parameters with both astrophysical and local probes. We also make a comparative analysis on how each data probe may constrain our parametrisation. We perform a Bayesian analysis by comparing the predictions of the model with observations. The astrophysical datasets are composed of quasars spectra measurements, including the latest ESPRESSO data point, as well as Planck observations of the cosmic microwave background. We combine them with local results from atomic clocks and the MICROSCOPE experiment. The constraints placed on the quintessence parameter are consistent with a null variation of the field, i.e. compatible with a $\Lambda$CDM cosmology. The constraints on the coupling to the electromagnetic sector are dominated by the E\"otv\"os parameter local bound. The ESPRESSO data point is extremely important to study the cosmological evolution of $\alpha$ as it probes an interval of redshift not accessible to other observations. However, for this particular model, current available data favours a null variation of $\alpha$ resulting mostly from the strong MICROSCOPE limits.

The isomerism of molecules in the interstellar medium and the mechanisms behind it are essential questions in the chemistry of organic molecules in space. In particular, for the simple formic and thioformic acids, the low temperatures found in molecular clouds indicate that cis-trans isomerization in the gas-phase must be impeded. Reactions happening on top of interstellar dust grains may explain the isomer interconversion at low temperatures. We studied the isomerization processes of formic and thioformic acid susceptible to happen on the surface of interstellar dust grains and initiated by H abstraction reactions. Similarly, deuterium enrichment of the acids can occur by the same mechanism. Our objective is to shed light on both topics to increase our understanding of key precursors of organic molecules in space.

Luminosity, which is the total amount of radiant energy emitted by an object, is one of the most critical quantities in astrophysics for characterizing stars. Equally important is the temporal evolution of a star's luminosity because of its intimate connection with the stellar energy budget, large-scale convective motion, and heat storage in the stellar interior. Here, we model the solar luminosity by extending a semi-empirical total solar irradiance (TSI) model that uses solar-surface magnetism to reconstruct solar irradiance over the entire 4{\pi} solid angle around the Sun. This model was constrained by comparing its output to the irradiance in the Earth's direction with the measured TSI. Comparing the solar luminosity to the TSI on timescales from days to for cycles 23 and 24, we find poor agreement on short timescales (< solar rotation). On longer timescales, however, we find good agreement between the luminosity model and the TSI, which suggests that the extrapolation of luminosities to multi-cycle timescales based on TSI reconstructions may be possible. We show that the solar luminosity is not constant but varies in phase with the solar cycle. This variation has an amplitude of 0.14% from minimum to maximum for solar cycle 23. Considering the energetics in the solar convection zone, it is therefore obvious that a steady-state input from the radiative zone at the solar minimum level would lead to a gradual reduction in the energy content in the convection zone over multi-century timescales. We show that the luminosity at the base of the convection zone should be approximately 0.032% higher than that at the solar surface during solar minimum to maintain net energy equilibrium through the solar cycle. These different energy-input scenarios place constraints on the long-term evolution of the total solar irradiance and its impact on the solar forcing of climate variability.

Massive neutrinos suppress the growth of structure under their free-streaming scales. The effect is most prominent on small scales where the widely-used two-point statistics can no longer capture the full information. In this work, we study the signatures massive neutrinos leave on large-scale structure (LSS) as revealed by its morphological properties, which are fully described by $4$ Minkowski functionals (MFs), and quantify the constraints on the summed neutrino mass $M_{\nu}$ from the MFs, by using publicly available N-body simulations. We find the MFs provide important complementary information, and give tighter constraints on $M_{\nu}$ than the power spectrum. Specifically, depending on whether massive neutrinos are included in the density field (the `m' field) or not (the `cb' field), we find the constraint on $M_{\nu}$ from the MFs with a smoothing scale of $R_G=5 h^{-1}$Mpc is $48$ or $4$ times better than that from the power spectrum. When the MFs are combined with the power spectrum, they can improve the constraint on $M_{\nu}$ from the latter by a factor of 63 for the `m' field and 5 for the `cb' field. Notably, when the `m' field is used, the constraint on $M_{\nu}$ from the MFs can reach $0.0177$eV with a volume of $1(h^{-1}\rm Gpc)^3$, while the combination of the MFs and power spectrum can tighten this constraint to be $0.0133$eV, a $4.5\sigma$ significance on detecting the minimum sum of the neutrino masses. For the `m' field, we also find the $\sigma_8$ and $M_{\nu}$ degeneracy is broken with the MFs, leading to stronger constraints on all 6 cosmological parameters considered in this work than the power spectrum.

The Large High Altitude Air Shower Observatory~(LHAASO) recently reported the detection of gamma-ray emissions with energies up to $1.1~\textrm{PeV}$ from the Crab Nebula. Using the absence of vacuum Cherenkov effect by inverse-Compton electrons, we improve previous bounds to linear-order Lorentz invariance violation (LV) in the dispersion relations of electrons by $10^{4}$ times. We show that the LV effect on electrons is severely constrained, compatible with certain type of LV as expected by some models of quantum gravity~(QG), such as the string/D-brane inspired space-time foam. We argue that such models are supported by the Crab Nebula constraints from the LHAASO observations, as well as various LV phenomenologies for photons to date.

We construct $\alpha$-attractor versions of hybrid inflation models. In these models, the potential of the inflaton field $\varphi$ is uplifted by the potential of the second field $\chi$. This uplifting ends due to a tachyonic instability with respect to the field $\chi$, which appears when $\varphi$ becomes smaller than some critical value $\varphi_{c}$. In the large $N$ limit, these models have the standard universal $\alpha$-attractor predictions. In particular, $n_{s }= 1- {2 \over N}$ for the exponential attractors. However, in some special cases the large $N$ limit is reached only beyond the horizon, for $N \gtrsim 60$. This may change predictions for the cosmological observations. For any fixed $N$, in the limit of large uplift $V_{\rm up}$, or in the limit of large $\varphi_{c}$, we find another attractor prediction, $ n_s = 1$. By changing the parameters $V_{\rm up}$ and $\varphi_{c}$ one can continuously interpolate between the two attractor predictions $n_{s }= 1- {2 \over N}$ and $n_{s} = 1$.

Within quantum-gravity approaches and beyond, different mechanisms for singularity resolution in black holes exist. Under a set of assumptions that we spell out in detail, these mechanisms leave their imprint in shadow images of spherically symmetric black holes. We find that even current EHT accuracy is sufficient to place nontrivial constraints on the scale of new physics within one modified spacetime, if the EHT measurement of M87* is combined with an independent measurement of the black-hole mass. In other spacetimes, increased accuracy is required that the next-generation EHT may deliver. We show how the combination of $n=1$ and $n=2$ photon rings is a powerful probe of the spacetime geometry of regular black holes, even when considering astrophysical uncertainties in accretion disks. Further, we generate images containing a localized emission region, inspired by the idea of hotspots in accretion flows. Finally, we investigate the photon-ring structure of a horizonless object, which is characterized by either two or no photon spheres. We show how photon rings annihilate each other, when there is no photon sphere in the spacetime.

String theory enjoys an elevated role among quantum gravity theories, since it seems to be the most consistent UV completion of general relativity and the Standard Model. However, it is hard to verify the existence of this underlying theory on terrestrial accelerators. One way to probe string theory is to study its imprints on the low-energy effective inflationary Lagrangian, which are quantified in terms of high energy correction terms. It is highly likely, thus, to find higher order curvature terms combined with string moduli, that is scalar fields, since both these types of interactions and matter fields appear in string theory. In this work we aim to stress the probability that the inflationary dynamics are controlled by the synergy of scalar fields and higher order curvature terms. Specifically, we shall consider a well motivated quantum corrected canonical scalar field theory, with the quantum corrections being of the $\mathcal{R}^2$ type. The reason for choosing minimally coupled scalar theory is basically because if scalar fields are evaluated in their vacuum configuration, they will either be minimally coupled or conformally coupled. Here we choose the former case, and the whole study shall be performed in the string frame (Jordan frame), in contrast to similar studies in the literature where the Einstein frame two scalar theory is considered. We derive the field equations of the quantum-corrected theory at leading order and we present the form the slow-roll indices obtain for the quantum corrected theory. We exemplify our theoretical framework by using the quadratic inflation model, and as we show, the $\mathcal{R}^2$ quantum corrected quadratic inflation model produces a viable inflationary phenomenology, in contrast with the simple quadratic inflation model.

The Neoproterozoic Earth experienced at least two global-scale glaciations termed Snowball Earth events. 'Cap carbonates' were widely deposited after the events, but controversy surrounds their origin. Here, we apply the novel $\delta^{44/40}$Ca-$\delta^{88/86}$Sr multi-proxy to two Marinoan (ca. 635 Ma) cap carbonate sequences from Namibia and show that the rocks archive primary environmental signals deriving from a combination of seawater-glacial meltwater mixing and kinetic isotope effects. In an outer platform section, dolostone $\delta^{44/40}$Ca and $\delta^{88/86}$Sr values define a line predicted for kinetic mass-dependent isotope fractionation. This dolostone mostly precipitated from meltwater. Moreover, stratigraphically higher samples exhibiting the fastest precipitation rates correlate with elevated 87Sr/86Sr ratios, consistent with long-held expectations that a rapid deglacial weathering pulse forced cap carbonate formation. An inner-platform dolostone shows greater effects from water-mass mixing but still reveals that precipitation rates increased up-section. Overlying limestones show the greatest Ca and Sr contributions from seawater. Amplification of local coastal processes during global ice sheet collapse offers a simple but sufficient proposition to explain the Ca isotope heterogeneity of cap carbonates. Detection of kinetic isotope effects in the rock record provides a basis for developing the $\delta^{44/40}$Ca-$\delta^{88/86}$Sr multi-proxy as an indicator of saturation state and $p$CO$_2$.

The empirical law of Gladstone-Dale is insufficient for high-precision studies using the refractivity of a gas: this is not exactly proportional to its density, and the gas may not be properly described as perfect. An optical Mariotte temperature allows making a comparative analysis of the results given by various authors. The effect of hygrometry on the refractivity at visible wavelengths is historically traced and its small effect on the astronomical refraction angle numerically shown. Finally at infrared and radio wavelengths, the effects of the humidity in the lower atmosphere can be strong; as for the ionosphere, its curvature plays an essential role for the astronomical refraction angle unlike in the visible.

The fundamental quadrature governing light rays in a spherically symmetrical medium is first recalled. A rigorous discussion of some qualitative properties of its solutions follows, using the Young-Kattawar diagram which leads to a geometric formulation of the ray curvature. The case of an optical duct is deepened, analyzing transfer curves for different positions of the observer with respect to the duct. New analytical expressions for their wavelength dependence are derived, and their numerical consequences are coherent with computer simulations.

Neutron star mergers are very violent events involving extreme physical processes: dynamic, strong-field gravity, large magnetic field, very hot, dense matter, and the copious production of neutrinos. Accurate modeling of such a system and its associated multi-messenger signals, such as gravitational waves, short gamma ray burst, and kilonova, requires the inclusion of all these processes, and is increasingly important in light of advancements in multi-messenger astronomy generally, and in gravitational wave astronomy in particular (such as the development of third-generation detectors). Several general relativistic codes have been incorporating some of these elements with different levels of realism. Here, we extend our code MHDuet, which can perform large eddy simulations of magnetohydrodynamics to help capture the magnetic field amplification during the merger, and to allow for realistic equations of state and neutrino cooling via a leakage scheme. We perform several tests involving isolated and binary neutron stars demonstrating the accuracy of the code.

We investigate the long-term relaxation of one-dimensional (${1D}$) self-gravitating systems, using both kinetic theory and $N$-body simulations. We consider thermal and Plummer equilibria, with and without collective effects. All combinations are found to be in clear agreement with respect to the Balescu-Lenard and Landau predictions for the diffusion coefficients. Interestingly, collective effects reduce the diffusion by a factor ${\sim 10}$. The predicted flux for Plummer equilibrium matches the measured one, which is a remarkable validation of kinetic theory. We also report on a situation of quasi kinetic blocking for the same equilibrium.

The stochastic spectral expansion method offers a simple framework for calculations in de Sitter spacetimes. We show how to extend its reach to metastable vacuum states, both in the case when the potential is bounded from below, and when it is unbounded from below and therefore no stable vacuum state exists. In both cases, the decay rate of the metastable vacuum is given by the lowest non-zero eigenvalue associated to the Fokker-Planck equation. We show how the corresponding eigenfunction determines the field probability distribution which can be used to compute correlation functions and other observables in the metastable vacuum state.