Papers for Friday, Feb 04 2022

The diffuse supernova neutrino background (DSNB) is a powerful future tool to constrain core-collapse explosion mechanisms without observation of a nearby event, and the corresponding signal has been calculated for a variety of collapse models. For Supernova (SN) 1987A, a peculiar double neutrino burst was detected, but models for the double collapse have never been studied in the DSNB context. Here, we fill this gap and compare the DSNB signal expected in the Standard Collapse (SC) and the Double Collapse (DC) models in various future detectors, including Hyper-Kamiokande, JUNO, DUNE and the Large Baksan Neutrino Telescope (LBNT). We calculate the spectra of diffuse neutrinos and antineutrinos in the DC model and determine the rate of registered events as a function of energy of the detected particle, taking into account detector parameters. For each detector, we estimate the corresponding uncertainties and the background and compare the signals expected for the SC and DC models. We conclude that the combination of DUNE and LBNT data will have the highest sensitivity to discriminate between the SC and DC models.

We investigate Dark Energy by associating it with vacuum energy or Cosmological constant ${\Lambda}$ which is taken to be dynamic in nature. Our approach is phenomenological and falls within the domain of variable-$\Lambda$ Cosmology. However, motivated by quantum theory of metastable vacuum decay, we proposed a new phenomenological decay law of $\Lambda$(t) where $\Lambda$(t) is a superposition of constant and variable components viz. $\Lambda$(t) = $\Lambda_{C}$ + $\Lambda_{v}$ which is indicated by the word $"$hybrid dynamic$"$ in the title. By taking a simplified two-fluid scenario with the Universe consisting of Dark Energy and another major component, we found the solutions for three particular phenomenological expressions and made a parametrization of the model in terms of dilution parameter (the dilution parameter has been defined in the text as the exponent of scale factor in the expression of density of the other major component, representing the dilution of the component with the expansion of Universe in the presence of dynamic Dark Energy). For pressureless Dust and dynamic Dark Energy Universe, we found the present-day matter density ($\Omega_{m0}$) and dilution parameter (u) to be $\Omega_{m0}$ = 0.29 $\pm$ 0.03, u = 2.90 $\pm$ 0.54 at 1 $\sigma$ by analysing 580 supernova from Union 2.1 catalogue. The physical features of the model in regard to scale factor evolution, deceleration parameter, cosmic age has also been studied and parallels have been drawn with $\Lambda$CDM model. The status of Cosmological problems in the model has also been checked which showed that the model solves the Cosmological Constant Problem but the Coincidence problem still exists in the model.

Recent work has shown that UV-luminous reionization-era galaxies often exhibit strong Lyman-alpha emission despite being situated at redshifts where the IGM is thought to be substantially neutral. It has been argued that this enhanced Ly$\alpha$ transmission reflects the presence of massive galaxies in overdense regions which power large ionized bubbles. An alternative explanation is that massive galaxies shift more of their Ly$\alpha$ profile to large velocities (relative to the systemic redshift) where the IGM damping wing absorption is reduced. Such a mass-dependent trend is seen at lower redshifts, but whether one exists at $z\sim7$ remains unclear owing to the small number of existing systemic redshift measurements in the reionization era. This is now changing with the emergence of [CII]-based redshifts from ALMA. Here we report MMT/Binospec Ly$\alpha$ spectroscopy of eight UV-bright ($\mathrm{M_{UV}}^{}\sim-22$) galaxies at $z\simeq7$ selected from the ALMA REBELS survey. We detect Ly$\alpha$ in 4 of 8 galaxies and use the [CII] systemic redshifts to investigate the Ly$\alpha$ velocity profiles. The Ly$\alpha$ lines are significantly redshifted from systemic (average velocity offset=223 km/s) and broad (FWHM$\approx$300$-$650 km/s), with two sources showing emission extending to $\approx$750 km/s. We find that the broadest Ly$\alpha$ profiles are associated with the largest [CII] line widths, suggesting a potential link between the Ly$\alpha$ FWHM and the dynamical mass. Since Ly$\alpha$ photons at high velocities transmit efficiently through the $z=7$ IGM, our data suggest that velocity profiles play a significant role in boosting the Ly$\alpha$ visibility of the most UV-luminous reionization-era galaxies.

We study the stellar mass-to-halo mass relation at $z=0$ in 30 Milky Way-like systems down to the ultra-faint ($M_* < 10^5 M_\odot$) regime using the semi-analytic model A-SLOTH. A new model allows us to follow star formation and the stochastic stellar feedback from individually sampled Pop II stars. Our fiducial model produces consistent results with the stellar mass-to-halo mass relation derived from abundance matching and the observed cumulative stellar mass function above the observational completeness. We find a plateau in the stellar mass-to-halo mass relation in the ultra-faint regime. The stellar mass of this plateau tells us how many stars formed before supernovae occur and regulate further star formation, which is determined by the Pop~II star formation efficiency. We also find that the number of luminous satellites increases rapidly as $M_*$ decreases until $M_* \approx 10^4 M_\odot$. Finally, we find that the relative streaming velocity between baryons and dark matter at high redshift is important in determining the number of ultra-faint dwarf galaxies at $z=0$. The new model in A-SLOTH provides a framework to study the stellar properties and the formation history of metal-poor stars in Milky Way and its satellites.

Any viable cosmological framework has to match the observed proportion of early- and late-type galaxies. In this contribution, we focus on the distribution of galaxy morphological types in the standard model of cosmology (Lambda cold dark matter, $\Lambda$CDM). Using the latest state-of-the-art cosmological $\Lambda$CDM simulations known as Illustris, IllustrisTNG, and EAGLE, we calculate the intrinsic and sky-projected aspect ratio distribution of the stars in subhalos with stellar mass $M_* > 10^{10}\,M_\odot$ at redshift $z=0$. There is a significant deficit of intrinsically thin disk galaxies, which however comprise most of the locally observed galaxy population. Consequently, the sky-projected aspect ratio distribution produced by these $\Lambda$CDM simulations disagrees with the Galaxy And Mass Assembly (GAMA) survey and Sloan Digital Sky Survey at $\geq 12.52\sigma$ (TNG50-1) and $\geq 14.82\sigma$ (EAGLE50) confidence. The deficit of intrinsically thin galaxies could be due to a much less hierarchical merger-driven build-up of observed galaxies than is given by the $\Lambda$CDM framework. It might also arise from the implemented sub-grid models, or from the limited resolution of the above-mentioned hydrodynamical simulations. We estimate that an $8^5$ times better mass resolution realization than TNG50-1 would reduce the tension with GAMA to the $5.58\sigma$ level. Finally, we show that galaxies with fewer major mergers have a somewhat thinner aspect ratio distribution. Given also the high expected frequency of minor mergers in $\Lambda$CDM, the problem may be due to minor mergers. In this case, the angular momentum problem could be alleviated in Milgromian dynamics (MOND) because of a reduced merger frequency arising from the absence of dynamical friction between extended dark matter halos.

We produce several public void catalogs using a volume-limited subsample of the Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Using new implementations of three different void-finding algorithms, VoidFinder and two ZOBOV-based algorithms (VIDE and REVOLVER), we identify 1159, 534, and 518 cosmic voids with radii >10 Mpc/h, respectively, out to a redshift of z = 0.114. We compute effective radii and centers for all voids and find none with an effective radius >54 Mpc/h. The median void effective radius is 15-17 Mpc/h for all three algorithms. We extract and discuss several properties of the void populations, including radial density profiles, the volume fraction of the catalog contained within voids, and the fraction of galaxies contained within voids. Using 64 mock galaxy catalogs created from the Horizon Run 4 N-body simulation, we compare simulated and observed void properties and find good agreement between the SDSS~DR7 and mock catalog results.

ALMA observations have revealed the presence of dust in the first generations of galaxies in the Universe. However, the dust temperature $T_d$ remains mostly unconstrained due to the few available FIR continuum data at redshift $z>5$. This introduces large uncertainties in several properties of high-$z$ galaxies, namely their dust masses, infrared luminosities, and obscured fraction of star formation. Using a new method based on simultaneous [CII] 158$\mu$m line and underlying dust continuum measurements, we derive $T_ d$ in the continuum and [CII] detected $z\approx 7$ galaxies in the ALMA Large Project REBELS sample. We find $39\ \mathrm{K} < T_d < 58\ \mathrm{K}$, and dust masses in the narrow range $M_d = (0.9-3.6)\times 10^7 M_{\odot}$. These results allow us to extend for the first time the reported $T_d(z)$ relation into the Epoch of Reionization. We produce a new physical model that explains the increasing $T_ d(z)$ trend with the decrease of gas depletion time, $t_{dep}=M_g/\mathrm{SFR}$, induced by the higher cosmological accretion rate at early times; this hypothesis yields $T_d \propto (1+z)^{0.4}$. The model also explains the observed $T_d$ scatter at a fixed redshift. We find that dust is warmer in obscured sources, as a larger obscuration results in more efficient dust heating. For UV-transparent (obscured) galaxies, $T_d$ only depends on the gas column density (metallicity), $T_d \propto N_H^{1/6}$ ($T_d \propto Z^{-1/6}$). REBELS galaxies are on average relatively transparent, with effective gas column densities around $N_H \simeq (0.03-1)\times 10^{21} \mathrm{cm}^{-2}$. We predict that other high-$z$ galaxies (e.g. MACS0416-Y1, A2744-YD4), with estimated $T_d \gg 60$ K, are significantly obscured, low-metallicity systems. In fact $T_d$ is higher in metal-poor systems due to their smaller dust content, which for fixed $L_{ IR}$ results in warmer temperatures.

The Sun and sun-like stars commonly host the multi-million-Kelvin coronae and the 10,000-Kelvin chromospheres. These extremely hot gases generate X-ray and Extreme Ultraviolet emissions that may impact the erosion and chemistry of (exo)planetary atmospheres, influencing the climate and conditions of habitability. However, the mechanism of coronal and chromospheric heating is still poorly understood. While the magnetic field most probably plays a key role in driving and transporting energy from the stellar surface upwards, it is not clear if the atmospheric heating mechanisms of the Sun and active sun-like stars can be described in a unified manner. To this end, we report on a systematic survey of the responses of solar and stellar atmospheres to surface magnetic flux over a wide range of temperatures. By analyzing 10 years of multi-wavelength synoptic observations of the Sun, we reveal that the irradiance and magnetic flux show power-law relations with an exponent decreasing from above- to sub-unity as the temperature decreases from the corona to the chromosphere. Moreover, this trend indicating the efficiency of atmospheric heating can be extended to sun-like stars. We also discover that the power-law exponent has a solar cycle dependence, where it becomes smallest at activity maximum, probably due to the saturation of atmospheric heating. Our study provides observational evidence that the mechanism of atmospheric heating is universal among the Sun and sun-like stars, regardless of age or activity.

We present optical spectroscopy observations of 145 high-mass pre-main sequence candidates from the catalogue of Vioque et al. (2020). From these, we provide evidence for the Herbig nature of 128 sources. This increases the number of known objects of the class by $\sim50\%$. We determine the stellar parameters of these sources using the spectra and Gaia EDR3 data. The new sources are well distributed in mass and age, with 23 sources between $4$-$8$ M$_{\odot}$ and 32 sources above $8$ M$_{\odot}$. Accretion rates are inferred from H$\alpha$ and H$\beta$ luminosities for 104 of the new Herbigs. These accretion rates, combined with previous similar estimates, allow us to analyze the accretion properties of Herbig stars using the largest sample ever considered. We provide further support to the existence of a break in accretion properties at $\sim3$-$4$ M$_{\odot}$, which was already reported for the previously known Herbig stars. We re-estimate the potential break in accretion properties to be at $3.87^{+0.38}_{-0.96}$ M$_{\odot}$. As observed for the previously known Herbig stars, the sample of new Herbig stars independently suggests intense inner-disk photoevaporation for sources with masses above $\sim7$ M$_{\odot}$. These observations provide robust observational support to the accuracy of the Vioque et al. (2020) catalogue of Herbig candidates.

MS 0735.6+7421 ($z = 0.216$) is a massive cool core galaxy cluster hosting one of the most powerful active galactic nuclei (AGN) outbursts known. The radio jets of the AGN have carved out an unusually large pair of X-ray cavities, each reaching a diameter of $200$ kpc. This makes MS 0735.6+7421 a unique case to investigate active galactic nuclei feedback processes, as well as other cluster astrophysics at radio wavelengths. We present new low-radio-frequency observations of MS 0735.6+7421 taken with the Karl G. Jansky Very Large Array (VLA): 5 hours of P-band ($224-480$ MHz) and 5 hours of L-band ($1-2$ GHz) observations, both in C configuration. Our VLA P-band ($224-480$ MHz) observations reveal the presence of a new diffuse radio component reaching a scale of $\sim$ $900$ kpc in the direction of the jets and of $\sim$ $500$ kpc in the direction perpendicular to the jets. This component is centered on the cluster core and has a radio power scaled at $1.4$ GHz of $P_{1.4\text{ GHz}} = (4\pm2)\times 10^{24}$ WHz$^{-1}$. Its properties are consistent with those expected from a radio mini-halo as seen in other massive cool core clusters, although it may also be associated with radio plasma that has diffused out of the X-ray cavities. Observations at higher spatial resolution are needed to fully characterize the properties and nature of this component. We also suggest that if radio mini-halos originate from jetted activity, we may be witnessing the early stages of this process.

The XMM Cluster Survey (XCS) have developed a new Python module, X-ray: Generate and Analyse (XGA) to provide interactive and automated analyses of X-ray emitting sources observed by the XMM-Newton space telescope. XGA only requires that a set of cleaned, processed, event lists has been created, and (optionally) that a source detector has generated region lists for the observations. XGA is centered around the concept of making all available data easily accessible and analysable. The user provides information (e.g. RA, Dec, redshift) on the source they wish to investigate, and XGA locates relevant observations and generates required data products. This allows the user to quickly and easily complete common analyses using all relevant observations, thus being left free to focus on extracting the maximum scientific gain. XGA is centered around source and sample classes, which represent different types of X-ray emitting astrophysical objects and have different properties and methods relevant to that type of object. XGA also contains product classes, which provide interfaces to X-ray data products or information (such as radial profiles) that have been derived from them, with built-in methods for fitting, analysis, and visualisation. XGA can fit models to spectra (both global and annular) with XSPEC, measuring spectral properties such as temperature, photon index, and luminosity. XGA can also measure radial profiles of density and temperature for galaxy clusters, allowing the measurement of gas and total mass profiles of galaxy clusters. In the future, we will add support for X-ray telescopes other than XMM (e.g. Chandra, eROSITA), as well as the ability to perform multi-mission joint analyses. With the advent of new X-ray observatories such as eROSITA, XRISM, ATHENA, and Lynx, it is the perfect time for a new, open-source, software package that is open for anyone to use and scrutinise.

We report the discovery and characterization of TOI-1759~b, a temperate (400 K) sub-Neptune-sized exoplanet orbiting the M~dwarf TOI-1759 (TIC 408636441). TOI-1759 b was observed by TESS to transit on sectors 16, 17 and 24, with only one transit observed per sector, creating an ambiguity on the orbital period of the planet candidate. Ground-based photometric observations, combined with radial-velocity measurements obtained with the CARMENES spectrograph, confirm an actual period of $18.85019 \pm 0.00014$ d. A joint analysis of all available photometry and radial velocities reveal a radius of $3.17 \pm 0.10\,R_\oplus$ and a mass of $10.8 \pm 1.5\,M_\oplus$. Combining this with the stellar properties derived for TOI-1759 ($R_\star = 0.597 \pm 0.015\,R_\odot$; $M_\star = 0.606 \pm 0.020\,M_\odot$; $T_{\textrm{eff}} = 4065 \pm 51$ K), we compute a transmission spectroscopic metric (TSM) value of over 80 for the planet, making it a good target for transmission spectroscopy studies. TOI-1759 b is among the top five temperate, small exoplanets ($T_\textrm{eq} < 500$ K, $R_p < 4 \,R_\oplus$) with the highest TSM discovered to date. Two additional signals with periods of 80 d and $>$ 200 d seem to be present in our radial velocities. While our data suggest both could arise from stellar activity, the later signal's source and periodicity are hard to pinpoint given the $\sim 200$ d baseline of our radial-velocity campaign with CARMENES. Longer baseline radial-velocity campaigns should be performed in order to unveil the true nature of this long period signal.

Recent sub-arcsecond resolution surveys of the dust continuum emission from nearby protoplanetary disks showed a strong correlation between the sizes and luminosities of the disks. We aim to explain the origin of the (sub-)millimeter size-luminosity relation (SLR) between the $68\%$ effective radius ($r_{eff}$) of disks with their continuum luminosity ($L_{mm}$), with models of gas and dust evolution in a simple viscous accretion disk and radiative transfer calculations. We use a large grid of models ($10^{5}$ simulations) with and without planetary gaps, varying the initial conditions of the key parameters. We calculate the disk continuum emission and the effective radius for all models as a function of time. By selecting those simulations that continuously follow the SLR, we can derive constraints on the input parameters of the models. We confirm previous results that models of smooth disks in the radial drift regime are compatible with the observed SLR ($L_{mm}\propto r_{eff}^{2}$) but only smooth disks cannot be the reality. We show that the SLR is more widely populated if planets are present. However they tend to follow a different relation than smooth disks, potentially implying that a mixture of smooth and sub-structured disks are present in the observed sample. We derive a SLR ($L_{mm}\propto r_{eff}^{5/4}$) for disks with strong sub-structure. To be compatible with the SLR, models need to have an initially high disk mass ($\geq 2.5 \cdot 10^{-2}M_{\star}$) and low turbulence-parameter $\alpha$ values ($\leq 10^{-3}$). Furthermore, we find that the grain composition and porosity drastically affects the evolution of disks on the size-luminosity diagram where relatively compact grains that include amorphous carbon are favoured. Moreover, a uniformly optically thick disk with high albedo ($0.9$) that follows the SLR cannot be formed from an evolutionary procedure.

We present a detection of 21-cm emission from large-scale structure (LSS) between redshift 0.78 and 1.43 made with the Canadian Hydrogen Intensity Mapping Experiment (CHIME). Radio observations acquired over 102 nights are used to construct maps which are foreground filtered and stacked on the angular and spectral locations of luminous red galaxies (LRG), emission line galaxies (ELG), and quasars (QSO) from the eBOSS clustering catalogs. We find decisive evidence for a detection when stacking on all three tracers of LSS, with the logarithm of the Bayes Factor equal to 18.9 (LRG), 10.8 (ELG), and 56.3 (QSO). An alternative frequentist interpretation, based on the likelihood-ratio test, yields a detection significance of $7.1\sigma$ (LRG), $5.7\sigma$ (ELG), and $11.1\sigma$ (QSO). These are the first 21-cm intensity mapping measurements made with an interferometer. We constrain the effective clustering amplitude of neutral hydrogen (HI), defined as $\mathcal{A}_{\rm HI}\equiv 10^{3}\,\Omega_\mathrm{HI}\left(b_\mathrm{HI}+\langle\,f\mu^{2}\rangle\right)$, where $\Omega_\mathrm{HI}$ is the cosmic abundance of HI, $b_\mathrm{HI}$ is the linear bias of HI, and $\langle\,f\mu^{2}\rangle=0.552$ encodes the effect of redshift-space distortions at linear order. We find $\mathcal{A}_\mathrm{HI}=1.51^{+3.60}_{-0.97}$ for LRGs $(z=0.84)$, $\mathcal{A}_\mathrm{HI}=6.76^{+9.04}_{-3.79}$ for ELGs $(z=0.96)$, and $\mathcal{A}_\mathrm{HI}=1.68^{+1.10}_{-0.67}$ for QSOs $(z=1.20)$, with constraints limited by modeling uncertainties at nonlinear scales. We are also sensitive to bias in the spectroscopic redshifts of each tracer, and find a non-zero bias $\Delta\,v= -66 \pm 20 \mathrm{km/s}$ for the QSOs. We split the QSO catalog into three redshift bins and have a decisive detection in each, with the upper bin at $z=1.30$ producing the highest redshift 21-cm intensity mapping measurement thus far.

Cosmic rays (CRs) are thought to play an important role in galaxy evolution. We study their effect when coupled to other important sources of feedback, namely supernovae and stellar radiation. Using the RAMSES-RT code, we perform the first radiation-magnetohydrodynamics simulations of isolated disc galaxies with and without CRs. We study galaxies embedded in dark matter haloes of $10^{10}$, $10^{11}$ and $10^{12}\, \rm M_{\odot}$ with a maximum resolution of $9 \,\rm pc$. We find that CRs reduce star formation rate in our two dwarf galaxies by a factor 2, with decreasing efficiency with increasing galaxy mass. They increase significantly the outflow mass loading factor in all our galaxies and make the outflows colder. We study the impact of the CR diffusion coefficient, exploring values from $\kappa = 10^{27}$ to $\rm 3\times 10^{29}\, cm^2\, s^{-1}$. With lower $\kappa$, CRs remain confined for longer on small scales and are consequently efficient in suppressing star formation, whereas a higher diffusion coefficient reduces the effect on star formation and increases the generation of cold outflows. Finally, we compare CR feedback to a calibrated 'strong' supernova feedback model known to sufficiently regulate star formation in high-redshift cosmological simulations. We find that CR feedback is not sufficiently strong to replace this strong supernova feedback. As they tend to smooth out the ISM and fill it with denser gas, CRs also lower the escape fraction of Lyman continuum photons from galaxies.

Efforts are underway to use high-precision timing of pulsars in order to detect low-frequency gravitational waves. A limit to this technique is the timing noise generated by dispersion in the plasma along the line of sight to the pulsar, including the solar wind. The effects due to the solar wind vary with time, influenced by the change in solar activity on different time scales, ranging up to $\sim 11$ years for a solar cycle. The solar wind contribution depends strongly on the angle between the pulsar line of sight and the solar disk, and is a dominant effect at small separations. Although solar wind models to mitigate these effects do exist, they do not account for all the effects of the solar wind and its temporal changes. Since low-frequency pulsar observations are most sensitive to these dispersive delays, they are most suited to test the efficacy of these models and identify alternative approaches. Here, we investigate the efficacy of some solar wind models commonly used in pulsar timing using long-term, high-cadence data on 6 pulsars taken with the Long Wavelength Array, and compare them with an operational solar wind model. Our results show that stationary models of the solar wind correction are insufficient to achieve the timing noise desired by pulsar timing experiments, and we need to use non-stationary models, which are informed by other solar wind observations, to obtain accurate timing residuals.

We report the detection and characterization of the transiting sub-Neptune TOI-1759 b, using photometric time-series from TESS and near infrared spectropolarimetric data from SPIRou on the CFHT. TOI-1759 b orbits a moderately active M0V star with an orbital period of $18.849975\pm0.000006$ d, and we measure a planetary radius and mass of $3.06\pm0.22$ R$_\oplus$ and $6.8\pm2.0$ M$_\oplus$. Radial velocities were extracted from the SPIRou spectra using both the CCF and the LBL methods, optimizing the velocity measurements in the near infrared domain. We analyzed the broadband SED of the star and the high-resolution SPIRou spectra to constrain the stellar parameters and thus improve the accuracy of the derived planet parameters. A LSD analysis of the SPIRou Stokes $V$ polarized spectra detects Zeeman signatures in TOI-1759. We model the rotational modulation of the magnetic stellar activity using a GP regression with a quasi-periodic covariance function, and find a rotation period of $35.65^{+0.17}_{-0.15}$ d. We reconstruct the large-scale surface magnetic field of the star using ZDI, which gives a predominantly poloidal field with a mean strength of $18\pm4$ G. Finally, we perform a joint Bayesian MCMC analysis of the TESS photometry and SPIRou RVs to optimally constrain the system parameters. At $0.1176\pm0.0013$ au from the star, the planet receives $6.4$ times the bolometric flux incident on Earth, and its equilibrium temperature is estimated at $433\pm14$ K. TOI-1759 b is a likely gas-dominated sub-Neptune with an expected high rate of photoevaporation. Therefore, it is an interesting target to search for neutral hydrogen escape, which may provide important constraints on the planetary formation mechanisms responsible for the observed sub-Neptune radius desert.

Pulsar radio emission may be generated in pair discharges which fill the pulsar magnetosphere with plasma as an accelerating electric field is screened by freshly created pairs. In this Letter we develop a simplified analytic theory for the screening of the electric field in these pair discharges and use it to estimate total radio luminosity and spectrum. The discharge has three stages. First, the electric field is screened for the first time and starts to oscillate. Next, a nonlinear phase occurs. In this phase, the amplitude of the electric field experiences strong damping because the field dramatically changes the momenta of newly created pairs. This strong damping ceases, and the system enters a final linear phase, when the electric field can no longer dramatically change pair momenta. Applied to pulsars, this theory may explain several aspects of radio emission, including the observed luminosity, $L_{\rm{rad}} \sim 10^{28} \rm{erg} \, \rm{s}^{-1}$, and the observed spectrum, $S_\omega \sim \omega^{-1.4 \pm 1.0} $.

Modern observational surveys allow us to probe the phase space distribution function (DF) of wide binaries in the Solar neighbourhood. This DF exhibits non-trivial features, in particular a superthermal distribution of eccentricities for semimajor axes $a\gtrsim 10^3$AU. To interpret such features we must first understand how the binary DF is affected by dynamical perturbations, which typically fall into two classes: (i) stochastic kicks from passing stars, molecular clouds, etc. and (ii) secular torques from the Galactic tide. Here we isolate effect (ii) and calculate the time-asymptotic, phase-mixed DF for an ensemble of wide binaries under quadrupole-order tides. For binaries wide enough that the phase-mixing assumption is valid, none of our results depend explicitly on semimajor axes, masses, etc. We show that unless the initial DF is both isotropic in binary orientation and thermal in eccentricity, then the final phase-mixed DF is always both anisotropic and non-thermal. However, the only way to produce a superthermal DF under phase mixing is for the initial DF to itself be superthermal.

A methodology to provide the polarization angle requirements for different sets of detectors, at a given frequency of a CMB polarization experiment, is presented. The uncertainties in the polarization angle of each detector set are related to a given bias on the tensor-to-scalar ratio $r$ parameter. The approach is grounded in using a linear combination of the detector sets to obtain the CMB polarization signal. In addition, assuming that the uncertainties on the polarization angle are in the small angle limit (lower than a few degrees), it is possible to derive analytic expressions to establish the requirements. The methodology also accounts for possible correlations among detectors, that may originate from the optics, wafers, etc. The approach is applied to the LiteBIRD space mission. We show that, for the most restrictive case (i.e., full correlation of the polarization angle systematics among detector sets), the requirements on the polarization angle uncertainties are of around 1 arcmin at the most sensitive frequency bands (i.e., $\approx 150$ GHz) and of few tens of arcmin at the lowest (i.e., $\approx 40$ GHz) and highest (i.e., $\approx 400$ GHz) observational bands. Conversely, for the least restrictive case (i.e., no correlation of the polarization angle systematics among detector sets), the requirements are $\approx 3$ times %smaller less restrictive than for the previous scenario. At the global and the telescope levels, polarization angle knowledge of a few arcmins is sufficient for correlated global systematic errors and can be relaxed by a factor of two for fully uncorrelated errors in detector polarization angle. The reported uncertainty levels are needed in order to have the bias on $r$ due to systematics below the limit established by the LiteBIRD collaboration.

Current observational facilities have yet to conclusively detect $10^3 - 10^4 M_{\odot}$ intermediate mass black holes (IMBHs) that fill in the evolutionary gap between early universe seed black holes and $z \sim 0$ supermassive black holes. Dwarf galaxies present an opportunity to reveal active IMBHs amidst persistent star formation. We introduce photoionization simulations tailored to address key physical uncertainties: coincident vs. non-coincident mixing of IMBH and starlight excitation, open vs. closed surrounding gas cloud geometries, and different AGN SED shapes. We examine possible AGN emission line diagnostics in the optical and mid-IR, and find that the diagnostics are often degenerate with respect to the investigated physical uncertainties. In spite of these setbacks, and in contrast to recent work, we are able to show that [O III]/H$\beta$ typically remains bright for dwarf AGN powered by IMBHs down to $10^3 M_{\odot}$. Dwarf AGN are predicted to have inconsistent star-forming and Seyfert/LINER classifications using the most common optical diagnostics. In the mid-IR, [O IV] 25.9$\mu$m and [Ar II] 6.98$\mu$m are less sensitive to physical uncertainties than are optical diagnostics. Based on these emission lines, we provide several mid-IR emission line diagnostic diagrams with demarcations for separating starbursts and AGN with varying levels of activity. The diagrams are valid over a wide range of ionization parameters and metallicities out to $z\sim0.1$, so will prove useful for future JWST observations of local dwarf AGN in the search for IMBHs. We make our photoionization simulation suite freely available.

Asteroid impacts with the Earth may have played an essential role in the emergence of life on Earth through their creation of favorable niches for life, changes to the atmosphere and delivery of water. Consequently, we suggest two potential requirements for life in an exoplanetary system: first, that the system has an asteroid belt, and second, that there is a mechanism to drive asteroids to impact the terrestrial habitable planet. Since in the solar system, the $\nu_6$ secular resonance has been shown to have been important in driving these impacts, we explore how the masses and locations of two giant planets determine the location and strength of this secular resonance. Examining observed exoplanetary systems with two giant planets, we find that a secular resonance within the asteroid belt region may not be uncommon. Hence the solar system is somewhat special, but the degree of fine-tuning that may be necessary for the emergence of life is not excessive. Finally, with $n$-body simulations, we show that when the two giant planets are close to the 2:1 mean motion resonance, the asteroid belt is unstable but this does not lead to increased asteroid delivery.

Water clouds are expected to form on Y dwarfs and giant planets with equilibrium temperatures near or below that of Earth, drastically altering their atmospheric compositions and their albedos and thermal emission spectra. Here we use the 1D Community Aerosol and Radiation Model for Atmospheres (CARMA) to investigate the microphysics of water clouds on cool substellar worlds to constrain their typical particle sizes and vertical extent, taking into consideration nucleation and condensation, which have not been considered in detail for water clouds in H/He atmospheres. We compute a small grid of Y dwarf and temperate giant exoplanet atmosphere models with water clouds forming through homogeneous nucleation and heterogeneous nucleation on cloud condensation nuclei composed of meteoritic dust, organic photochemical hazes, and upwelled potassium chloride cloud particles. We present comparisons with the Ackerman & Marley parameterization of cloud physics to extract the optimal sedimentation efficiency parameter (f$_{sed}$) using Virga. We find that no Virga model replicates the CARMA water clouds exactly and that a transition in f$_{sed}$ occurs from the base of the cloud to the cloud top. Furthermore, we generate simulated thermal emission and geometric albedo spectra and find large, wavelength-dependent differences between the CARMA and Virga models, with different gas absorption bands reacting differently to the different cloud distributions and particularly large differences in the M band. Therefore, constraining the vertically-dependent properties of water clouds will be essential to estimating the gas abundances in these atmospheres.

To quantitatively study the driving mechanisms of magnetospheric convection in the magnetotail lobes on a global scale, we utilize data from the ARTEMIS spacecraft in the deep tail and the Cluster spacecraft in the near tail. Previous work demonstrated that, in the lobes near the Moon, we can estimate the convection by utilizing ARTEMIS measurements of lunar ions velocity. In this paper, we analyze these datasets with machine learning models to determine what upstream factors drive the lobe convection in different magnetotail regions and thereby understand the mechanisms that control the dynamics of the tail lobes. Our results show that the correlations between the predicted and test convection velocities for the machine learning models (> 0.75) are much better than those of the multiple linear regression model (~ 0.23 - 0.43). The systematic analysis reveals that the IMF and magnetospheric activity play an important role in influencing plasma convection in the global magnetotail lobes.

Strongly lensed type Ia supernovae (SNe Ia) are expected to have some advantages in measuring time delays of multiple images, and so they have a great potential to be developed into a powerful late-universe cosmological probe. In this paper, we simulate a sample of lensed SNe Ia with time-delay measurements in the era of the Legacy Survey of Space and Time (LSST). Based on the distance sum rule, we use lensed SNe Ia to implement model-independent constraints on the Hubble constant $H_0$ and cosmic curvature parameter $\Omega_K$ in the late universe. We find that if 20 lensed SNe Ia could be observed, the constraint on $H_{0}$ is better than the measurement by the SH0ES collaboration. When the event number of lensed SNe Ia increases to 100, the constraint precision of $H_{0}$ is comparable with the result from \emph{Planck} 2018 data. Considering 200 lensed SNe Ia events as the optimistic estimation, we obtain $\Delta H_0=0.33$ $\rm km\ s^{-1}\ Mpc^{-1}$ and $\Delta\Omega_K=0.053$. In addition, we also simulate lensed quasars in different scenarios to make a comparison and we find that they are still a useful cosmological probe even though the constraint precision from them is much less than that obtained from lensed SNe Ia. In the era of LSST, the measurements of time delay from both lensed SNe Ia and lensed quasars are expected to yield the results of $\Delta H_0=0.26 ~\rm km\ s^{-1}\ Mpc^{-1}$ and $\Delta\Omega_K=0.044$.

To reveal natures of the dark matter (DM) particles, a gamma-ray signal produced in annihilation processes of DM into the standard model particles has been one of the major probes. The cross-correlation between highly DM dominated structures, such as local dwarf galaxies, and observed photons in the direction of the structures, has been explored and provided stringent constraints on the annihilation rate. In our previous work, we have shown that it is sufficient to know the distance distribution of the galaxy sample and individual distance measurement is not required for constraining the annihilation rate. In this work, we apply the method to low-surface brightness galaxies (LSBGs) with unknown individual redshifts provided from the Dark Energy Survey (DES) Year 3 data. With all DES-LSBGs of about 24,000 objects, we find that the upper limits of the cross section for bb channel with 95% C.L. is 3*10^-25 cm^3/s at DM mass of 100 GeV. To be more conservative, we remove about 7000 LSBGs within 1 degree from resolved gamma-ray point sources and then the constraint becomes 30 % weaker than the one with all samples in all DM mass ranges.

We report the results from the broad-band X-ray monitoring of the new Galactic black hole candidate MAXI J1803$-$298 with the MAXI/GSC and Swift/BAT during its outburst. After the discovery on 2021 May 1, the soft X-ray flux below 10 keV rapidly increased for $\sim 10$ days and then have been gradually decreasing over 5 months. At the brightest phase, the source exhibited the state transition from the low/hard state to the high/soft state via the intermediate state. The broad-band X-ray spectrum during the outburst was well described with a disk blackbody plus its thermal or non-thermal Comptonization. Before the transition the source spectrum was described by a thermal Comptonization component with a photon index of $\sim 1.7$ and an electron temperature of $\sim 30$ keV, whereas a strong disk blackbody component was observed after the transition. The spectral properties in these periods are consistent with the low/hard state and the high/soft state, respectively. A sudden flux drop with a few days duration, unassociated with a significant change in the hardness ratio, was found in the intermediate state. A possible cause of this variation is that the mass accretion rate rapidly increased at the disk transition, which induced a strong Compton-thick outflow and scattered out the X-ray flux. Assuming a non-spinning black hole, we estimated a black hole mass of MAXI J1803$-$298 as $5.8 \pm 0.4~(\cos i/\cos 70^\circ)^{-1/2} (D/8~\mathrm{kpc})~M_\odot$ (where $i$ and $D$ are the inclination angle and the distance) from the inner disk radius obtained in the high/soft state.

The observable characteristics and subsequent evolution of young stellar populations is dominated by their massive stars. As our understanding of those massive stars and the factors affecting their evolution improves, so our interpretation of distant, unresolved stellar systems can also advance. As observations increasingly probe the distant Universe, and the rare low metallicity starbursts nearby, so the opportunity arises for these two fields to complement one another, and lead to an improved conception of both stars and galaxies. Here we review the current state of the art in modelling of massive star dominated stellar populations, and discuss their applications and implications for interpreting the distant Universe. Our principle findings include: - Binary evolutionary pathways must be included to understand the stellar populations in early galaxies. - Observations constraining the extreme ultraviolet spectrum of early galaxies are showing that current models are incomplete. The best current guess is that some form of accretion onto compact remnants is required. - The evolution and fates of very massive stars, of the order of 100Msun and above, may be key to fully understand aspects of early galaxies.

We analyzed 12-year Fermi Large Area Telescope $\gamma$-ray data in the inner Galaxy centered at (l=30$^{\deg}$, b=0$^{\deg}$.) and (l=330$^{\deg}$, b=0$^{\deg}$). We found significant hardening of the spectrum of the diffuse $\gamma$-ray emission in these regions as previously reported. We further deduced that the diffuse $\gamma$ rays can be divided into two components from the likelihood analysis. One component is associated with the total gas column density and reveals a soft spectrum, while the other is associated with the HII gas and presents a hard spectrum. Assuming the diffuse $\gamma$-ray emissions are mainly produced through the interaction between cosmic rays (CRs) and the ambient gas, these two components are produced by the CR populations with spectral indices of 2.8 ("soft") and 2.3 ("hard"), respectively. We argue that the hard CR population may come from the vicinity of the CR accelerators. The soft CR population has a similar spectral shape and density as measured in the solar neighborhood, which implies a uniform CR "sea" with a similar density and spectral shape in the Galaxy.

We report the detection of flaring events in NGC 4395 ULX1, a nearby ultraluminous X-ray source(ULX), for the first time, using recent XMM-NEWTON observations. The flaring episodes are spectrally harder than the steady emission periods, resulting in higher fractional variability in the high energy regime. A thin Keplerian and a slim accretion disk provide the best-fit continuum for XMM-NEWTON spectra. In all observations, the presence of a broad Gaussian emission feature around ~0.9 keV suggests a strong wind outflow in this ULX. The flaring spectra correspond to higher slim disk temperatures due to higher mass accretion rate under an advection-dominated accretion scenario. The luminosity-temperature profile in different flux states is consistent with the theoretical predictions for a slim accretion disk in the case of super-Eddington accretion onto a stellar-mass compact object. The unperturbed absorption column and wind outflow line emission during flaring events suggest that the site of origin of these flares is mainly the advection flow in the inner region of the accretion disk.

Global models of planet formation tend to begin with an initial set of planetary embryos for the sake simplicity. While this approach gives valuable insights on the evolution of the initial embryos, the initial distribution itself is a bold assumption. Limiting oneself to an initial distribution may neglect essential physics that precedes, or follows said initial distribution. We wish to investigate the effect of dynamic planetary embryo formation on the formation of planetary systems. The presented framework begins with an initial disk of gas, dust and pebbles. The disk evolution, the formation of planetesimals and the formation of planetary embryos is modeled consistently. Embryos then grow by pebble, planetesimal and eventually gas accretion. Planet disk interactions and N-body dynamics with other simultaneously growing embryos is included in the framework. We show that the formation of planets can occur in multiple consecutive phases. Earlier generations grow massive by pebble accretion but are subject to fast type I migration and thus accretion to the star. The later generations of embryos that form grow to much smaller masses by planetesimal accretion, as the amount of pebbles in the disk has vanished. The formation history of planetary systems may be far more complex than an initial distribution of embryos could reflect. The dynamic formation of planetary embryos needs to be considered in global models of planet formation to allow for a complete picture of the systems evolution.

The shape of the ambient circumnuclear medium (ACM) density profile can probe the history of accretion onto the central supermassive black hole in galaxies and the circumnuclear environment. However, due to the limitation of the instrument resolution, the density profiles of the ACM for most of galaxies remain largely unknown. In this work, we propose a novel method to measure the ACM density profile of active galactic nucleus (AGN) by the equilibrium between the radiation pressure on the warm absorbers (WAs, a type of AGN outflows) and the drag pressure from the ACM. We study the correlation between the outflow velocity and ionization parameter of WAs in each of the five Seyfert 1 galaxies (NGC 3227, NGC 3783, NGC 4051, NGC 4593, and NGC 5548), inferring that the density profile of the ACM is between n\propto r^-1.7 and n \propto r^-2.15 (n is number density and r is distance) from 0.01 pc to pc scales in these five AGNs. Our results indicate that the ACM density profile in Seyfert 1 galaxies is steeper than the prediction by the spherically symmetric Bondi accretion model and the simulated results of the hot accretion flow, but more in line with the prediction by the standard thin disk model.

Black hole low mass X-ray binaries (BH LMXBs) can launch powerful outflows in the form of discrete ejecta. Observing the entire trajectory of these ejecta allows us to model their motion with great accuracy, and this is essential for measuring their physical properties. In particular, observing the final deceleration phase, often poorly sampled, is fundamental to obtain a reliable estimate of the jet's energy. During its 2019/2020 outburst, the BH LMXB MAXI J1348$-$630 launched a single-sided radio-emitting jet that was detected at large scales after a strong deceleration due to the interaction with the interstellar medium (ISM). We successfully modelled the jet motion with a dynamical external shock model, which allowed us to constrain the jet initial Lorentz factor $\Gamma_0 = 1.85^{+0.15}_{-0.12}$, inclination angle $\theta = 29.3_{-3.2}^{+2.7}$ deg and ejection date $t_{\rm ej} = 21.5_{-3.0}^{+1.8}$ (MJD $-$ $58500$). Under simple assumptions on the jet opening angle and on the external ISM density, we find that the jet has a large initial kinetic energy $E_0 = 4.6^{+20.0}_{-3.4} \times 10^{46}$ erg, far greater than what commonly measured for LMXBs from the jet's synchrotron emission. This implies that discrete ejecta radiate away only a small fraction of their total energy, which is instead transferred to the environment. The jet power estimate is larger than the simultaneous available accretion power, and we present several options to mitigate this discrepancy. We infer that MAXI J1348$-$630 is likely embedded in an ISM cavity with internal density $n = 0.0010^{+0.0005}_{-0.0003}$ cm$^{-3}$ and radius $R_{\rm c} = 0.61^{+0.11}_{-0.09}$ pc, which could have been produced by the system's previous activity, as proposed for other BH LMXBs.

We performed NH$_3\ (J,K)=(1,1),(2,2),$ and $(3,3)$ mapping observations toward the Galactic massive star-forming region Sh 2-255 and Sh 2-257 using the Nobeyama 45-m telescope as a part of the KAGONMA (KAgoshima Galactic Object survey with the Nobeyama 45-metre telescope by Mapping in Ammonia lines) project. NH$_3$ (1,1) has an intensity peak at the cluster S255 N, is distributed over 3 pc $\times$ 2 pc and is located between two HII regions. The kinetic temperature derived from the NH$_3 (2,2)/(1,1)$ ratio was $\sim 35$ K near the massive cluster S255 IR. These clusters also show emission with a large line width of $\sim$ 3-4 km s$^{-1}$. Based on the reported data we suggest that NH$_3$ gas in these regions is affected by stellar feedback from embedded YSO clusters in S255 IR and S255 N. We also detected NH$_3$ (1,1) emission in a region west of the main gas clump at the location of a concentration of Class II YSOs adjacent to the HII regions Sh 2-254. The presence of Class II YSOs implies $\sim$ 2 Myr of star formation, younger than Sh 2-254 ($\sim 5$ Myr), thus we suggest that star formation in the western region could be influenced by the older HII region Sh 2-254.

Here we present an improved algorithm to model the serpentinization process in planetesimals in the early Solar system. Although it is hypothesized that serpentinization-like reactions played an important role in the thermal evolution of planetesimals, few and restricted models are available in this topic. These processes may be important as the materials involved were abundant in these objects. Our model is based on the model by (Gobi & Kereszturi 2017), and contains improvements in the consideration of heat capacities and lithospheric pressure, and in the calculation of the amount of interfacial water. Comparison of our results with previous calculations shows that there are significant differences in the e.g. the serpentinization time -- the time necessary to consume most of the reactants at specific initial conditions -- or the amount of heat produced by this process. In a simple application we show that in icy bodies, under some realistic conditions, below the melting point of water ice, serpentinization reaction using interfacial water may be able to proceed and eventually push the local temperature above the melting point to start a 'runaway' serpentinization. According to our calculations in objects with radii R $\gtrsim$ 200 km serpentinization might have quickly reformed nearly the whole interior of these bodies in the early Solar system.

In order to extract the precise physical information encoded in the gravitational and electromagnetic signals from powerful neutron-star merger events, we need to include as much of the relevant physics as possible in our numerical simulations. This presents a severe challenge, given that many of the involved parameters are poorly constrained. In this paper we focus on the role of nuclear reactions. Combining a theoretical discussion with an analysis connecting to state-of-the-art simulations, we outline multiple arguments that lead to a reactive system being described in terms of a bulk viscosity. The results demonstrate that in order to properly account for nuclear reactions, future simulations must be able to handle different regimes where rather different assumptions/approximations are appropriate. We also touch upon the link to models based on the large-eddy-strategy required to capture turbulence.

We analyse Gaia EDR3 and re-calibrated HST proper motion data from the core-collapsed and non core-collapsed globular clusters NGC 6397 and NGC 3201, respectively, with the Bayesian mass-orbit modelling code MAMPOSSt-PM. We use Bayesian evidence and realistic mock data sets constructed with AGAMA to select between different mass models. In both clusters, the velocities are consistent with isotropy within the extent of our data. We robustly detect a dark central mass (DCM) of roughly 1000 solar masses in both clusters. Our MAMPOSSt-PM fits strongly prefer an extended DCM in NGC 6397, while only presenting a mild preference for it in NGC 3201, with respective sizes of a roughly one and a few per cent of the cluster effective radius. We explore the astrophysics behind our results with the CMC Monte Carlo N-body code, whose snapshots best matching the phase space observations lead to similar values for the mass and size of the DCM. The internal kinematics are thus consistent with a population of hundreds of massive white dwarfs in NGC 6397, and roughly 100 segregated stellar-mass black holes in NGC 3201, as previously found with CMC. Such analyses confirm the accuracy of both mass-orbit modelling and Monte Carlo N-body techniques, which together provide more robust predictions on the DCM of globular clusters (core-collapsed or not). This opens possibilities to understand a vast range of interesting astrophysical phenomena in clusters, such as fast radio bursts, compact object mergers, and gravitational waves.

Context: Properly modelling scattering by interstellar dust grains requires a good characterisation of the scattering phase function. The Henyey-Greenstein phase function has become the standard for describing anisotropic scattering by dust grains, but it is a poor representation of the real scattering phase function outside the optical range. Aims: We investigate alternatives for the Henyey-Greenstein phase function that would allow the scattering properties of dust grains to be described. Our goal is to find a balance between realism and complexity: the scattering phase function should be flexible enough to provide an accurate fit to the scattering properties of dust grains over a wide wavelength range, and it should be simple enough to be easy to handle, especially in the context of radiative transfer calculations. Methods: We fit various analytical phase functions to the scattering phase function corresponding to the BARE-GR-S model, one of the most popular and commonly adopted models for interstellar dust. We weigh the accuracy of the fit against the number of free parameters in the analytical phase functions. Results: We confirm that the Henyey-Greenstein phase functions poorly describe scattering by dust grains, particularly at ultraviolet (UV) wavelengths, with relative differences of up to 50%. The Draine phase function alleviates this problem at near-infrared (NIR) wavelengths, but not in the UV. The two-term Reynolds-McCormick phase function, recently advocated in the context of light scattering in nanoscale materials and aquatic media, describes the BARE-GR-S data very well, but its five free parameters are degenerate. We propose a simpler phase function, the two-term ultraspherical-2 (TTU2) phase function, that also provides an excellent fit to the BARE-GR-S phase function over the entire UV-NIR wavelength range. (Abridged)

High quality low-frequency radio surveys have the promise of advancing our understanding of many important topics in astrophysics, including the life cycle of active galactic nuclei (AGN), particle acceleration processes in jets, the history of star formation, and exoplanet magnetospheres. Currently leading low-frequency surveys reach an angular resolution of a few arcseconds. However, this resolution is not yet sufficient to study the more compact and distant sources in detail. Sub-arcsecond resolution is therefore the next milestone in advancing these fields. The biggest challenge at low radio frequencies is the ionosphere. If not adequately corrected for, ionospheric seeing blurs the images to arcsecond or even arcminute scales. Additionally, the required image size to map the degree-scale field of view of low-frequency radio telescopes at this resolution is far greater than what typical soft- and hardware is currently capable of handling. Here we present for the first time (to the best of our knowledge) widefield sub-arcsecond imaging at low radio frequencies. We derive ionospheric corrections in a few dozen individual directions and apply those during imaging efficiently using a recently developed imaging algorithm (arXiv:1407.1943, arXiv:1909.07226). We demonstrate our method by applying it to an eight hour observation of the International LOw Frequency ARray (LOFAR) Telescope (ILT) (arXiv:1305.3550). Doing so we have made a sensitive $7.4\ \mathrm{deg}^2$ $144\ \mathrm{MHz}$ map at a resolution of $0.3''$ reaching $25\ \mu\mathrm{Jy\ beam}^{-1}$ near the phase centre. The estimated $250,000$ core hours used to produce this image, fit comfortably in the budget of available computing facilities. This result will enable future mapping of the entire northern low-frequency sky at sub-arcsecond resolution.

Turbulence is expected to generate characteristic velocity and density fluctuations in the interstellar medium (ISM). Using HI survey data, we distinguish these contributions as the basis of comparisons with theoretical and magnetohydrodynamical (MHD) simulations. We used HI4PI multiphase observations and Gaussian components representing three HI phases, the cold, warm, and unstable lukewarm medium (CNM, WNM and LNM, respectively) to deduce characteristic fluctuations in turbulent density and velocity fields. We applied the velocity decomposition algorithm (VDA) for separating such fluctuations in the position-position-velocity (PPV) space. The VDA extends the velocity channel analysis (VCA) by Lazarian and Pogosyan and predicts that turbulent velocity and density fields are statistically uncorrelated. Applying the VDA decomposition to the observational data (here, HI4PI) yields - contrary to our expectations - a significant correlation between velocity and density fields. All HI phases contribute to this correlation. The correlations in the velocity wings suffer from attenuation caused by uncorrelated noise. Both VCA and VDA predict that dispersions from fluctuations in narrow velocity channels scale in proportion to the average intensity. Observed brightness temperature fluctuations follow, however, a square-root scaling - as expected for a sum of normally distributed random sources. Fluctuations in HI channel maps at high spatial frequencies are observed to be dominated by density structures, which is in opposition to the results from MHD simulations, where small-scale structures are predicted to have been generated by velocity caustics. The VCA predictions and VDA results on MHD simulations are not compatible with the HI observations. The observed turbulent velocity and density fields in the ISM are not statistically uncorrelated, but they do reveal significant interrelations.

The dynamical evolution of the Prometheus and Pandora pair of satellites is chaotic, with a short 3.3 years Lyapunov time. It is known that the anti-alignment of the apses line of Prometheus and Pandora, which occurs every 6.2 years, is a critical configuration that amplifies their chaotic dynamical evolution. However, the mutual interaction between Prometheus and Pandora is not enough to explain the longitudinal lags observed by the Hubble Space Telescope. The main goal of the current work is to identify the main contributors to the chaotic dynamical evolution of the Prometheus-Pandora pair beyond themselves. Therefore, in this work, we first explore the sensibility of this dynamical system to understand it numerically and then build numerical experiments to reach our goals. We identified that almost all major satellites of the Saturn system play a significant role in the evolution of Prometheus' and Pandora's orbits.

We present a spectroscopic analysis of Type I superluminous supernova (SLSN-I), SN 2018bsz. While it closely resembles SLSNe-I, the multi-component H$\alpha$ line appearing at $\sim30$ d post-maximum is the most atypical. The H$\alpha$ is characterised by two emission components, one at $+3000$ km/s and a second at $-7500$ km/s, with a third, near-zero velocity component appearing after a delay. The blue and central components can be described by Gaussian profiles of intermediate width, but the red component is significantly broader and Lorentzian. The blue component evolves towards lower velocity before fading at $100$ d post-peak, concurrently with a light curve break. Multi-component profiles are observed in other hydrogen lines including Pa$\beta$, and in lines of Ca II and He I. Spectropolarimetry obtained before (10.2 d) and after (38.4 d) the appearance of the H lines show a large shift on the Stokes $Q$ -- $U$ plane consistent with SN 2018bsz undergoing radical changes in its geometry. Assuming the SN is almost unpolarised at 10.2 d, the continuum polarisation at 38.4 d reaches $P \sim1.8\%$ implying a highly asymmetric configuration. We propose that the observed evolution of SN 2018bsz can be explained by highly aspherical CSM. After the SN explosion, the CSM is quickly overtaken by the ejecta, but as the photosphere starts to recede, the different CSM regions re-emerge producing the peculiar line profiles. Based on the first appearance of H$\alpha$, we can constrain the distance of the CSM to be less than $430$ AU, or even lower ($<87$ AU) if the pre-peak plateau is related to an eruption that created the CSM. The presence of CSM has been inferred for other SLSNe-I. However, it is not clear whether the rare properties of SN 2018bsz can be generalised for SLSNe-I or whether they are the result of an uncommon evolutionary path, possibly involving a binary companion.

Observations of the SNR Cassiopeia A (Cas A) show asymmetries in the reverse shock that cannot be explained by models describing a remnant expanding through a spherically symmetric wind of the progenitor star. We investigate whether a past interaction of Cas A with a massive asymmetric shell of the circumstellar medium can account for the observed asymmetries. We performed 3D MHD simulations that describe the remnant evolution from the SN to its interaction with a circumstellar shell. The initial conditions are provided by a 3D neutrino-driven SN model whose morphology resembles Cas A. We explored the parameter space of the shell, searching for a set of parameters able to produce reverse shock asymmetries at the age of 350 years analogous to those observed in Cas A. The interaction of the remnant with the shell can match the observed reverse shock asymmetries if the shell was asymmetric with the densest portion in the nearside to the northwest (NW). According to our models, the shell was thin with radius 1.5 pc. The reverse shock shows the following asymmetries at the age of Cas A: i) it moves inward in the observer frame in the NW region, while it moves outward in other regions; ii) the geometric center of the reverse shock is offset to the NW by 0.1 pc from the geometric center of the forward shock; iii) the reverse shock in the NW region has enhanced nonthermal emission because, there, the ejecta enter the reverse shock with a higher relative velocity (between 4000 and 7000 km/s) than in other regions (below 2000 km/s). Our findings suggest the interaction of Cas A with an asymmetric circumstellar shell between 180 and 240 years after the SN event. We suggest that the shell was, most likely, the result of a massive eruption from the progenitor star that occurred about 10^5 years prior to core-collapse. We estimate a total mass of the shell of approximately 2.6 Msun.

A relativistic jet has been produced in the single well-localised binary neutron star (BNS) merger detected to date in gravitational waves (GWs), and the local rates of BNS mergers and short gamma-ray bursts are of the same order of magnitude. This suggests that jet formation is not a rare outcome for BNS mergers, and we show that this intuition can be turned into a quantitative constraint: at least about $1/3$ of GW-detected BNS mergers, and at least about $1/5$ of all BNS mergers, should produce a successful jet (90\% credible level). Whether a jet is launched depends on the properties of the merger remnant and of the surrounding accretion disc, which in turn are a function of the progenitor binary masses and equation of state (EoS). The fraction of jets in the population therefore carries information about the binary component mass distribution and EoS. Under the assumption that a jet can only be produced by a black hole remnant surrounded by a non-negligible accretion disc, we show how the jet fraction can be used to place a joint constraint on the space of BNS component mass distributions and EoS. The result points to a broad mass distribution, with particularly strong support for masses in the $1.3-1.6\,\mathrm{M_\odot}$ range. The constraints on the EoS are shallow, but we show how they will tighten as the knowledge on the jet fraction and mass distribution improve. We also discuss how to extend the method to include future BNS events, with possibly uncertain jet associations.

Some of the earliest work on neutrino astronomy was accomplished by a class of underground detectors primarily designed for particle physics goals . These detectors used inexpensive water to obtain the large masses needed to observe the very low interaction rates expected from neutrinos. They exploited the relatively large light attenuation length and the index of refraction of the water to get a very inexpensive cost per thousand tons of detector. The results obtained from these pioneering neutrino detectors have included real time observation of solar neutrinos, supernova neutrinos, and atmospheric neutrinos. Searches for neutrino point sources, dark matter and primordial magnetic monopoles were also made using them.

Asymmetric debris discs have been found around stars other than the Sun; asymmetries are sometimes attributed to perturbations induced by unseen planets. The presence or absence of asymmetries in our own trans-Neptunian belt remains controversial. The study of sensitive tracers in a sample of objects relatively free from the perturbations exerted by the four known giant planets and most stellar flybys may put an end to this debate. The analysis of the distribution of the mutual nodal distances of the known extreme trans-Neptunian objects (ETNOs) that measure how close two orbits may get to each other could be such a game changer. Here, we use a sample of 51 ETNOs together with random shufflings of this sample and two unbiased scattered-disc orbital models to confirm a statistically significant (62 sigma) asymmetry between the shortest mutual ascending and descending nodal distances as well as the existence of multiple highly improbably (p < 0.0002) correlated pairs of orbits with mutual nodal distances as low as 0.2 au at 152 au from the Solar system's barycentre or 1.3 au at 339 au. We conclude that these findings fit best with the notion that trans-Plutonian planets exist.

We present a method to derive the logarithmic extinction coefficient in optical wavelengths using the emission lines of HeI. Using this procedure we can avoid selection biases when studying regions with different surface brightness and we can obtain better measurements of temperature lines, for example [SIII]6312.

We review the present methodology for the use of Sulfur as global metallicity tracer in galaxies, which allows performing a complete abundance analysis using mainly the red to near infrared spectral region, and extending the range of directly derived abundances up to 5 times the S solar photospheric value. The empirical calibration of Sulfur via the S23 parameter is also reviewed.

We present the wide field slitless spectroscopy mode of the NIRISS instrument on the James Webb Space Telescope. This mode employs two orthogonal low-resolution (resolving power $\approx 150$) grisms in combination with a set of six blocking filters in the wavelength range 0.8 to $2.3\,\mu$m to provide a spectrum of almost every source across the field-of-view. When combined with the low background, high sensitivity and high spatial resolution afforded by the telescope, this mode will enable unprecedented studies of the structure and evolution of distant galaxies. We describe the performance of the as-built hardware relevant to this mode and expected imaging and spectroscopic sensitivity. We discuss operational and calibration procedures to obtain the highest quality data. As examples of the observing mode usage, we present details of two planned Guaranteed Time Observations programs: The Canadian NIRISS Unbiased Cluster Survey (CANUCS) and The NIRISS Survey for Young Brown Dwarfs and Rogue Planets.

Despite helium abundance (AHe = nH/nHe) is ~ 8 % at the solar photospheric/chromospheric heights, AHe can be found to exceed 8% in interplanetary coronal mass ejections (ICMEs) on many occasions. Although various factors like interplanetary shocks, chromospheric evaporation and "sludge removal" have been separately invoked in the past to address the AHe enhancements in ICMEs, none of these processes could explain the variability of AHe in ICMEs comprehensively. Based on extensive analysis of 275 ICME events, we show that there is a solar activity variation of ICME averaged AHe values. The investigation also reveals that the first ionization potential effect as well as coronal temperature are not the major contributing factors for AHe enhancements in ICMEs. Investigation on concurrent solar flares and ICME events for 63 cases reveals that chromospheric evaporation in tandem with gravitational settling determine the AHe enhancements and variabilities beyond 8% in ICMEs. While chromospheric evaporation releases the helium from chromosphere into the corona, the gravitationally settled helium is thrown out during the ICMEs.We show that the intensity and timing of the preceding flares from the same active region from where the CME erupts are important factors to understand the AHe enhancements in ICMEs.

We present a large VLT/MUSE mosaic (3.8 x 3.8 kpc) of the nearby spiral galaxy M83, with a spatial resolution ~20 pc. We obtained the kinematics of the stars and ionised gas, and compared them with molecular gas kinematics from ALMA CO(2-1). We separated the ionised gas into HII regions and diffuse ionised gas (DIG) and determined the fraction of Ha luminosity originating from the DIG (f_DIG). We observe that both stars and gas trace the galactic disk rotation, as well as a fast-rotating nuclear component, likely connected to secular processes driven by the galactic bar. In the gas kinematics, we observe a stream east of the nucleus, redshifted with respect to the disk. The stream is surrounded by an extended ionised gas region with enhanced velocity dispersion and a high ionisation state, which is largely consistent with being ionised by slow shocks. We interpret this feature as either the superposition of the disk and an extraplanar layer of DIG, or as a bar-driven inflow of shocked gas. A double Gaussian component fit to the Ha line also reveals the presence of a nuclear biconic structure whose axis of symmetry is perpendicular to the bar. The two cones appear blue- and redshifted along the line of sight and stand out for having an Ha emission separated by up to 200 km s-1 from that of the disk, and a high velocity dispersion ~80-200 km s-1. At the far end of the cones, we observe that the gas is consistent with being ionised by shocks. These features had never been observed before in M83; we postulate that they are tracing a starburst-driven outflow shocking into the surrounding ISM. Finally, we obtain f_DIG ~ 13% in our field of view. We inspect the emission of the HII regions and DIG in `BPT' diagrams, finding that in HII regions photoionisation accounts for 99.8% of the Ha flux, whereas the DIG has a mixed contribution from photoionisation (94.9%) and shocks (5.1%). [abridged]

We study the consequences of an enlarged 4 parameter dynamical dark energy (4pDE) equation of state using the latest Planck, BAO, and Pantheon supernovae data. This parameterization of the dark energy equation of state incorporates a generic non-linear monotonic evolution of the dark energy equation of state, where the four parameters are the early and the present value of the equation of state, the transition scale factor, and the sharpness of the transition. In this study, we use SH0ES $M_B$ prior and the KIDS/Viking $S_8$ prior while keeping the neutrino mass $\Sigma m_\nu$ as a free parameter. We show that in this case, the dynamical dark energy 4pDE model can bring down the Hubble tension to $\sim 2.5 \sigma$ level and the $S_8$ tension to $\sim 1.5 \sigma$ level when tested against Planck, BAO, and Pantheon supernovae data together. We also compare our results with the well-explored CPL model. We find that the present data can not constrain all the four dark energy equations of state parameters ensuring the fact that the present observations do not demand a complex non-linear multi-parameter evolution of the time-dependent DE equation of state. We also report that with SH0ES $M_B$ and KIDS/Viking $S_8$ prior 4pDE and CPL model favours a non-zero value for the neutrino mass parameter at the most at $\sim 1 \sigma$ level ($\Sigma m_\nu \sim 0.2 \pm 0.1$ eV).

We have analyzed the circumnuclear ring of the spiral galaxy NGC7742 in order to understand its formation and evolution. We have obtained gaseous abundances, characterized the interstellar medium of the clusters and studied the properties of the ionizing clusters. We have also implemented a new methodology using the red wavelength range of optical spectra, with the purpose of understanding how star formation evolves in high metallicity environments.

We present a new study of the Algol-type eclipsing binary system AS Eri based on the combination of the MOST and TESS light curves and a collection of very precise radial velocities obtained with the spectrographs HERMES operating at the Mercator telescope, La Palma, and TCES operating at the Alfred Jensch telescope, Tautenburg. The primary component is an A3 V-type pulsating, mass-accreting star. We fitted the light and velocity data with the package PHOEBE, and determined the best-fitting model adopting the configuration of a semi-detached system. The orbital period has been improved using a recent (O-C) analysis and the phase shift detected between both light curves to the value 2.6641496 $\pm$ 0.0000001 days. The absence of any cyclic variation in the (O-C) residuals confirms the long-term stability of the orbital period. Furthermore, we show that the models derived for each light curve separately entail small differences, e.g. in the temperature parameter T$_{\rm eff,2}$. The high quality of the new solutions is illustrated by the residuals. We obtained the following absolute component parameters: L$_1$ = 14.125~L$_{\odot}$, M$_1$ = 2.014~M$_{\odot}$, R$_1$ = 1.733~R$_{\odot}$, log g$_1$ = 4.264, L$_2$ = 4.345~L$_{\odot}$, M$_2$ = 0.211~M$_{\odot}$, R$_2$ = 2.19~R$_{\odot}$, log g$_2$ = 3.078~ with T$_{\rm eff,2}$/T$_{\rm eff,1}$ = 0.662 $\pm$ 0.002. Although the orbital period appears to be stable on the long term, we show that the light-curve shape is affected by a years-long modulation which is most probably due to the magnetic activity of the cool companion.

It has become clear that early career astrophysics researchers (doctoral researchers, post-docs, etc) have a very diverse appreciation of their career, with some declaring it the best job that you can have and others suffering from overwork, harrassment and stress from the precarity of their job, and associated difficulties. In order to establish how astrophysics researchers, primarily in France, experience their career, we sent out a survey to understand the impact that their job has on their well-being. 276 people responded to the survey. Whilst around half of the respondents expressed pleasure derived from their career, it is clear that many (early career) researchers are suffering due to overwork, with more than a quarter saying that they work in excess of 50 hours per week and 2\% in excess of 90 h per week. Almost 30\% professed to having suffered harrassment or discrimination in the course of their work. Further, whilst only 20\% had suffered mental health issues before starting their career in astrophysics, $\sim$45\% said that they suffered with mental health problems since starting in astrophysics. Here we provide results from the survey as well as possible avenues to explore and a list of recommendations to improve (early) careers in astrophysics.

Gravitational waves (GWs) generate oscillating electromagnetic effects in the vicinity of external electric and magnetic fields. We discuss this phenomenon with a particular focus on reinterpreting the results of axion haloscopes based on lumped-element detectors, which probe GWs in the 100 kHz-100 MHz range. Measurements from ABRACADABRA and SHAFT already place bounds on GWs, although the present strain sensitivity is weak. However, we demonstrate that the sensitivity scaling with the volume of such instruments is significant - faster than for axions - and so rapid progress will be made in the future. With no modifications, DMRadio-m$^3$ will have a GW strain sensitivity of $h \sim 10^{-20}$ at 200 MHz. A simple modification of the pickup loop used to readout the induced magnetic flux can parametrically enhance the GW sensitivity, particularly at lower frequencies.

Stars that pass close to the supermassive black holes located in the center of galaxies can be violently disrupted by tidal forces, leading to flares that are observed as bright transient events in sky surveys. The rate for these events to occur depends on the black hole spins, which in turn can be affected by ultra-light bosons due to superradiance. We perform a detailed analysis of these effects and show that searches for stellar tidal disruptions have a significant potential to uncover the existence of ultra-light bosons. In particular, we find that upcoming stellar tidal disruption rate measurements by the Vera Rubin Observatory's Legacy Survey of Space and Time can be used to either discover or rule out bosons with masses ranging from $10^{-20}$ to $10^{-18}$ eV. Our analysis also indicates that these measurements may be used to constrain a variety of supermassive black hole spin distributions and determine if close-to maximal spins are preferred.

The universal relaxation bound suggests that the relaxation times of perturbed thermodynamical systems is bounded from below by the simple time-times-temperature (TTT) quantum relation $\tau\times T\geq {{\hbar}\over{\pi}}$. It is known that some perturbation modes of near-extremal Kerr black holes in the regime $MT_{\text{BH}}/\hbar\ll m^{-2}$ are characterized by normalized relaxation times $\pi\tau\times T_{\text{BH}}/\hbar$ which, in the approach to the limit $MT_{\text{BH}}/\hbar\to0$, make infinitely many oscillations with a tiny constant amplitude around $1$ and therefore cannot be used directly to verify the validity of the TTT bound in the entire parameter space of the black-hole spacetime (Here $\{T_{\text{BH}},M\}$ are respectively the Bekenstein-Hawking temperature and the mass of the black hole, and $m$ is the azimuthal harmonic index of the linearized perturbation mode). In the present compact paper we explicitly prove that all rapidly-spinning Kerr black holes respect the TTT relaxation bound. In particular, using analytical techniques, it is proved that all black-hole perturbation modes in the complementary regime $m^{-1}\ll MT_{\text{BH}}/\hbar\ll1$ are characterized by relaxation times with the simple dimensionless property $\pi\tau\times T_{\text{BH}}/\hbar\geq1$.

A minimal model based on the Covariant Physical Couplings (CPC) framework for gravity is proposed. The CPC framework is based on the assumptions of a metric-compatible four-dimensional Riemannian manifold where a covariantly conserved stress-energy tensor acts as source of the field equations which are formally the same as Einstein field equations, but where the couplings $\{ G, c,\Lambda \}$ are allowed to vary simultaneously. The minimal CPC model takes $\Lambda$ as a genuine constant while $c$ and $G$ vary in an entangled way that is consistent with Bianchi identity and the aforementioned assumptions. The model is constrained using the most recent galaxy cluster gas mass fraction observational data. Our result indicates that the functions $c(z)$ and $G\left(z\right)=G_{0}\left(c/c_{0}\right)^{4}$ are compatible with constant couplings for the two different parameterizations of $c=c(z)$ adopted here.

Quasi-normal oscillation modes of neutron stars provide a means to probe their interior composition using gravitational wave astronomy. We compute the frequencies and damping times of composition-dependent core g-modes of neutron stars containing quark matter employing linearized perturbative equations of general relativity. We find that ignoring background metric perturbations due to the oscillating fluid, as in the Cowling approximation, underestimates the g-mode frequency by up to 10% for higher mass stars, depending on the parameters of the nuclear equation of state and how the mixed phase is constructed. The g-mode frequencies are well-described by a linear scaling with the central lepton (or combined lepton and quark) fraction for nucleonic (hybrid) stars. Our findings suggest that neutron stars with and without quarks are manifestly different with regards to their quasi-normal g-mode spectrum, and may thus be distinguished from one another in future observations of gravitational waves from merging neutron stars.

Measurements of cosmogenic particles at various locations and altitudes are becoming increasingly important in view of worldwide interests in rare signals for search of new physics. In this work, we report measurement of muon zenith angle distributions and integrated flux using a portable setup of four one-meter long liquid scintillator bars. Each scintillator bar is read out from both sides via photomultiplier tubes followed by an 8-channel Digitizer. We exploit energy deposition and excellent timing of scintillators to construct two dimensional tracks and hence angles of charged particles. We use liquid scintillators since they have an added advantage of pulse shape discrimination (PSD) which can be used for detecting muon induced particles. The energy deposition, time window of event and PSD cuts are used to reduce the random as well as correlated backgrounds. In addition, we propose three track quality parameters which are applied to obtain a clean muon spectrum. The zenith angle measurement is performed upto $60^\circ$. With our improved analysis, we demonstrate that a setup of 3 bars can also be used for quicker and precise measurements. The vertical muon flux measured is $66.70 \pm 0.36 \pm 1.50$ with $n=2.10 \pm 0.05 \pm 0.25 $ in $\cos^n \theta$ at the location of Mumbai, India ($19^{\circ}$N, $72.9^{\circ}$E) at Sea level with a muon momentum above $255$ GeV/$c$. The muon flux has dependence on various factors, the most prominents are latitudes, altitudes and momentum cut of muon so that portable setup like this can be a boon for such measurements at various locations.

The objective of the proposed MAQRO mission is to harness space for achieving long free-fall times, extreme vacuum, nano-gravity, and cryogenic temperatures to test the foundations of physics in macroscopic quantum experiments. This will result in the development of novel quantum sensors and a means to probe the foundations of quantum physics at the interface with gravity. Earlier studies showed that the proposal is feasible but that several critical challenges remain, and key technologies need to be developed. These new technologies will open up the potential for achieving additional science objectives. The proposed research campaign aims to advance the state of the art and to perform the first macroscopic quantum experiments in space. Experiments on the ground, in micro-gravity, and in space will drive the proposed research campaign during the current decade to enable the implementation of MAQRO within the subsequent decade.

Antimatter particles such as positrons and antiprotons abound in the cosmos. Much less common are light antinuclei, composed of antiprotons and antineutrons, which can be produced in our galaxy via high-energy cosmic-ray collisions with the interstellar medium or could also originate from the annihilation of the still undiscovered dark-matter particles. On Earth, the only way to produce and study antinuclei with high precision is to create them at high-energy particle accelerators like the Large Hadron Collider (LHC). Though the properties of elementary antiparticles have been studied in detail, knowledge of the interaction of light antinuclei with matter is rather limited. This work focuses on the determination of the disappearance probability of \ahe\ when it encounters matter particles and annihilates or disintegrates. The material of the ALICE detector at the LHC serves as a target to extract the inelastic cross section for \ahe\ in the momentum range of $1.17 \leq p < 10$ GeV/$c$. This inelastic cross section is measured for the first time and is used as an essential input to calculations of the transparency of our galaxy to the propagation of $^{3}\overline{\rm He}$ stemming from dark-matter decays and cosmic-ray interactions within the interstellar medium. A transparency of about 50% is estimated using the GALPROP program for a specific dark-matter profile and a standard set of propagation parameters. For cosmic-ray sources, the obtained transparency with the same propagation scheme varies with increasing $^{3}\overline{\rm He}$ momentum from 25% to 90%. The absolute uncertainties associated to the $^{3}\overline{\rm He}$ inelastic cross section measurements are of the order of 10%$-$15%. The reported results indicate that $^{3}\overline{\rm He}$ nuclei can travel long distances in the galaxy, and can be used to study cosmic-ray interactions and dark-matter decays.