Papers for Wednesday, Jun 30 2021

The LDN 1172/1174 cloud complex in the Cepheus Flare region presents a hub-filament structure with the reflection nebula, NGC 7023, illuminated by a Herbig Be star, HD 200775, which consists of the hub with a $\sim$5 pc long narrow filament attached to it. Formation of a sparse cluster of low- and intermediate-mass stars is presently taking place in the hub. The aim of this work is to map the magnetic field geometry of LDN 1172/1174 to understand the role played by the field lines in the formation of the molecular cloud. We made R-band polarization measurements of 249 stars projected on the entire LDN 1172/1174 cloud complex to map the geometry of the magnetic field of this region. The magnetic field geometry constructed from our R-band polarization measurements is found to be parallel to the elongated structure inferred from the column density distribution of the cloud produced using the Herschel images. Our R-band polarization measurements are found to be in good agreement with those obtained from Planck. There is evidence of a possible distortion of the magnetic fields toward the northwestern part of the cloud by HD 200775. The magnetic field strength is estimated as $\sim$30 $\mu$G. The estimated star formation rate (SFR)/mass of 2.0$\pm$1.3 \%Myr$^{-1}$ and 0.4$\pm$0.3 \%Myr$^{-1}$ for LDN 1172/1174 and the neighboring cloud complex, LDN 1147/1158, respectively, are found to be consistent with the mean SFR/mass found for the clouds with magnetic field orientations parallel and perpendicular to their elongated structures, respectively. These results support earlier findings that the clouds with magnetic field lines parallel to their long axes seem to have higher SFRs compared to those with the magnetic field orientation perpendicular to the cloud elongation.

We use hydrodynamical simulations of two Milky Way-mass galaxies to demonstrate the impact of cosmic-ray pressure on the kinematics of cool and warm circumgalactic gas. Consistent with previous studies, we find that cosmic-ray pressure can dominate over thermal pressure in the inner 50 kpc of the circumgalactic medium (CGM), creating an overall cooler CGM than that of similar galaxy simulations run without cosmic rays. We generate synthetic sightlines of the simulated galaxies' CGM and use Voigt profile fitting methods to extract ion column densities, Doppler-b parameters, and velocity centroids of individual absorbers. We directly compare these synthetic spectral line fits with HST/COS CGM absorption-line data analyses, which tend to show that metallic species with a wide range of ionization potential energies are often kinematically aligned. Compared to the Milky-Way simulation run without cosmic rays, the presence of cosmic-ray pressure in the inner CGM creates narrower OVI absorption features and broader SiIII absorption features, a quality which is more consistent with observational data. Additionally, because the cool gas is buoyant due to nonthermal cosmic-ray pressure support, the velocity centroids of both cool and warm gas tend to align in the simulated Milky Way with feedback from cosmic rays. Our study demonstrates that detailed, direct comparisons between simulations and observations, focused on gas kinematics, have the potential to reveal the dominant physical mechanisms that shape the CGM.

The vector Apodizing Phase Plate (vAPP) is a class of pupil plane coronagraph that enables high-contrast imaging by modifying the Point Spread Function (PSF) to create a dark hole of deep flux suppression adjacent to the PSF core. Here, we recover the known brown dwarf HR 2562 B using a vAPP coronagraph, in conjunction with the Magellan Adaptive Optics (MagAO) system, at a signal-to-noise of S/N = 3.04 in the lesser studied L-band regime. The data contained a mix of field and pupil-stabilised observations, hence we explored three different processing techniques to extract the companion, including Flipped Differential Imaging (FDI), a newly devised Principal Component Analysis (PCA)-based method for vAPP data. Despite the partial field-stabilisation, the companion is recovered sufficiently to measure a 3.94 $\mu$m narrow-band contrast of (3.05$\pm$1.00) $\times$ 10$^{-4}$ ($\Delta$m$_{3.94 {\mu}m}$ = 8.79$\pm$0.36 mag). Combined with archival GPI and SPHERE observations, our atmospheric modelling indicates a spectral type at the L/T transition with mass M = 29$\pm$15 M$_{\text{Jup}}$, consistent with literature results. However, effective temperature and surface gravity vary significantly depending on the wavebands considered (1200$\leq$T$_{\text{eff}}$(K)$\leq$1700 and 4.0$\leq$log(g)(dex)$\leq$5.0), reflecting the challenges of modelling objects at the L/T transition. Observations between 2.4-3.2 $\mu$m will be more effective in distinguishing cooler brown dwarfs due to the onset of absorption bands in this region. We explain that instrumental scattered light and wind-driven halo can be detrimental to FDI+PCA and thus must be sufficiently mitigated to use this processing technique. We thus demonstrate the potential of vAPP coronagraphs in the characterisation of high-contrast substellar companions, even in sub-optimal conditions, and provide new, complementary photometry of HR 2562 B.

We derive a new model for neutrino-plasma interactions in an expanding universe that incorporates the collective effects of the neutrinos on the plasma constituents. We start from the kinetic description of a multi-species plasma in the flat Friedmann-Robertson-Walker metric, where the particles are coupled to neutrinos through the charged- and neutral-current forms of the weak interaction. We then derive the fluid equations and specialize our model to (a) the lepton epoch, where we consider a pair electron-positron plasma interacting with electron (anti-)neutrinos, and (b) after the electron-positron annihilation, where we model an electron-proton plasma and take the limit of slow ions and inertia-less electrons to obtain a set of neutrino-electron magnetohydrodynamics (NEMHD) equations. In both models, the dynamics of the plasma is affected by the neutrino motion through a ponderomotive force and, as a result, new terms appear in the induction equation that can act as a source for magnetic field generation in the early universe. A brief discussion on the possible applications of our model is proposed.

We take advantage of the publicly available LEGA-C spectroscopic survey to measure the stellar population properties of 140 individual massive and passive galaxies at $z\sim0.7$. We develop and publicly release PyLick, a flexible python code to measure UV to near-IR spectral indices. With PyLick we study the H/K ratio as a new diagnostic based on the pseudo-Lick CaII H and K indices, and find that a cut in ${\rm H/K}<1.1$ can be used jointly with other criteria to select (or verify the purity of) samples of passive galaxies. By combining photometric and spectroscopic criteria, we select a reliable sample of passively evolving galaxies. We constrain single-burst stellar ages, metallicities $\mathrm{[Z/H]}$, and $\mathrm{[\alpha/Fe]}$ with an optimized set of Lick indices, exploring in detail the robustness of our measurement against different combinations. Even without imposing cosmological priors, the derived ages follow a clear trend compatible with the expected cosmological aging of the Universe. We observe no significant redshift evolution for the metal abundance with respect to the values derived at $z=0$, with median $\mathrm{[Z/H]}=0.08\pm0.18$ and $\mathrm{[\alpha/Fe]}=0.13\pm0.11$. Finally, we find positive correlations between $\log\mathrm{age}$, $\mathrm{[Z/H]}$, $\mathrm{[\alpha/Fe]}$ and the stellar velocity dispersion, with slopes of ($0.48\pm0.14$), ($0.26\pm0.17$), and ($0.23\pm0.11$), respectively; the small scatter of $<0.2$ dex points to rather homogeneous and short star formation histories. Overall, these results confirm and extend low-redshift findings of a mass-downsizing evolution. This work further strengthens the possibility of selecting pure samples of passive galaxies to be exploited reliably as cosmic chronometers to place independent cosmological constraints.

We study scalar-tensor-tensor and tensor-scalar-scalar three point cross correlations generated by the dynamics of a transiently rolling spectator axion-$\rm U(1)$ gauge field model during inflation. In this framework, tensor and scalar fluctuations are sourced by gauge fields at the non-linear level due to gravitational interactions, providing a chiral background of gravitational waves while keeping the level of scalar fluctuations at the observationally viable levels at CMB scales. We show that the gravitational couplings between the observable sector and gauge fields can also mediate strong correlations between scalar and tensor fluctuations, generating an amplitude for the mixed type three-point functions that is parametrically larger compared to the single-field realizations of inflation. As the amplification of the gauge field sources are localized around the time of horizon exit, the resulting mixed bispectra are peaked close to the equilateral configurations. The shape dependence, along with the scale dependence and the parity violating nature of the mixed bispectra can serve as a distinguishing feature of the underlying axion-gauge field dynamics and suggest a careful investigation of their signatures on the CMB observables including cross correlations between temperature T and E,B polarization modes.

Improving our knowledge of the global properties of the Milky Way (MW) is critical for connecting the detailed measurements only possible within our own galaxy to our understanding of the broader galaxy population. To do this we utilise Gaussian Process Regression (GPR), a method for making predictions from sparsely-sampled, multi-dimensional datasets that can capture both local and large-scale trends. We train GPR models to map from galaxy properties that are well-measured for both the MW and external galaxies to broadband fluxes in a wide variety of photometric bands. The galaxy properties we use to predict photometric characteristics are stellar mass, apparent axis ratio, star formation rate, bulge-to-total ratio, disk scale length, and bar vote fraction. We use these models to estimate the global UV (GALEX $FUV/NUV$), optical (SDSS $ugriz$) and IR (2MASS $JHKs$ and WISE $W1/W2/W3/W4$) photometric properties for the MW as they would be measured from outside, resulting in a full UV-to-IR spectral energy distribution (SED). We show for the first time that the MW must lie in the star-forming regime in standard UV and IR diagnostic diagrams, in contrast to its position in the green valley in optical colour-mass diagrams. This is characteristic of the population of red spiral galaxies, suggesting that the MW may be a member of that class. Although each GPR model only predicts one band at a time, we find that the resulting MW UV-IR SED is in good agreement with SEDs of local spirals with characteristics broadly similar to the MW, suggesting that we can combine the individual band fluxes with confidence. Our MW UV-IR SED will serve as a valuable tool for reconstructing the Milky Way's star formation history via the same tools used for external galaxies, allowing for comparisons of results from $\textit{in situ}$ measurements to those from the methods used for extra-galactic objects.

We study cool neutral gas traced by NaD absorption in 140 local ($\rm z<0.1)$ early-type ``red geyser'' galaxies. These galaxies show unique signatures in spatially-resolved strong-line emission maps that have been interpreted as large-scale active galactic nuclei driven ionized winds. To investigate the possible fuel source for these winds, we examine the abundance and kinematics of cool gas ($\rm T \sim 100-1000 K$) inferred from Na I D absorption in red geysers and matched control samples drawn from SDSS-IV MaNGA. We find that red geysers host greater amounts of NaD-associated material. Substantial cool gas components are detected in more than $\rm 50 \%$ of red geysers (compared to 25\% of the control sample) going up to 78$\%$ for radio-detected red geysers. Our key result is that cool gas in red geysers is predominantly infalling. Among our 30 radio-detected red geysers, 86$\%$ show receding NaD absorption velocities (with respect to the systemic velocity) between $\rm 40 - 50~km~s^{-1}$. We verify this result by stacking NaD profiles across each sample which confirms the presence of infalling NaD velocities within red geysers ( $\sim\rm 40~km~s^{-1}$) with no velocity offsets detected in the control samples. Interpreting our observations as signatures of inflowing cool neutral clouds, we derive an approximate mass inflow rate of $\rm \dot{M}_{in} \sim 0.1 M_{\odot} yr^{-1}$, similar to that expected from minor merging and internal recycling. Some red geysers show much higher rates ($\rm \dot{M}_{in} \sim 5 M_{\odot} yr^{-1}$) that may indicate an ongoing accretion event.

The Millimeter-wave Intensity Mapping Experiment (mmIME) recently reported a detection of excess spatial fluctuations at a wavelength of 3 mm, which can be attributed to unresolved emission of several CO rotational transitions between $z\sim1-5$. We study the implications of this data for the high-redshift interstellar medium using a suite of state-of-the-art semianalytic simulations which have successfully reproduced many other sub-millimeter line observations across the relevant redshift range. We find that the semianalytic predictions are mildly in tension with the mmIME result, with a predicted CO power $\sim3.5\sigma$ below what was observed. We explore some simple modifications to the models which could resolve this tension. Increasing the molecular gas abundance at the relevant redshifts to $\sim10^8\ M_\odot\ \rm{Mpc}^{-3}$, a value well above that obtained from directly imaged sources, would resolve the discrepancy, as would assuming a CO-$H_2$ conversion factor $\alpha_{\rm{CO}}$ of $\sim1.5\ M_{\odot}$ K$^{-1}$ $(\rm{km}/\rm{s})^{-1}$ pc$^{2}$, a value somewhat lower than is commonly assumed. We go on to demonstrate that these conclusions are quite sensitive to the detailed assumptions of our simulations, highlighting the need for more careful modeling efforts as more intensity mapping data become available.

We present results from a comprehensive study of ultrafast outflows (UFOs) detected in a sample of fourteen quasars, twelve of which are gravitationally lensed, in a redshift range of 1.41-3.91, near the peak of the AGN and star formation activity. New XMM-Newton observations are presented for six of them which were selected to be lensed and contain a narrow absorption line (NAL) in their UV spectra. Another lensed quasar was added to the sample, albeit already studied, because it was not searched for UFOs. The remaining seven quasars of our sample are known to contain UFOs. The main goals of our study are to infer the outflow properties of high-z quasars, constrain their outflow induced feedback, study the relationship between the outflow properties and the properties of the ionizing source, and compare these results to those of nearby AGN. Our study adds six new detections ( > 99% confidence) of UFOs at z > 1.4, almost doubling the current number of cases. Based on our survey of six quasars selected to contain a NAL and observed with XMM-Newton, the coexistence of intrinsic UV NALs and UFOs is found to be significant in > 83% of these quasars suggesting a link between multiphase AGN feedback properties of the meso- and micro-scale. The kinematic luminosities of the UFOs of our high-z sample are large compared to their bolometric luminosities (median of L_K/L_Bol ~ 50%). This suggests they provide efficient feedback to influence the evolution of their host galaxies and that magnetic driving may be a significant contributor to their acceleration.

We report on the results of an XMM-Newton observation of the Supergiant Fast X-ray Transient (SFXT) IGR J08408-4503 performed in June 2020. The source is composed by a compact object (likely a neutron star) orbiting around an O8.5Ib-II(f)p star, LM Vel. The X-ray light curve shows a very low level of emission, punctuated by a single, faint flare. Analysis of spectra measured during the flare and during quiescence is performed. The quiescent state shows a continuum spectrum well deconvolved to three spectral models: two components are from a collisionally-ionized plasma (with temperatures kT1=0.24 keV and kT2=0.76 keV), together with a power law model (photon index of 2.55), dominating above 2 keV. The X-ray flux emitted at this lowest level is 3.2$\times10^{-13}$ erg/cm2/s (0.5-10 keV, corrected for the interstellar absorption), implying an X-ray luminosity of 1.85$\times10^{32}$ erg/s (at 2.2 kpc). The two temperature collisionally-ionized plasma is intrinsic to the stellar wind of the donor star, while the power law can be interpreted as emission due to residual, low level accretion onto the compact object. The X-ray luminosity contributed by the power law component only, in the lowest state, is (4.8$\pm{1.4})\times10^{31}$ erg/s, the lowest quiescent luminosity detected from the compact object in an SFXT. Thanks to this very faint X-ray state caught by XMM-Newton, X-ray emission from the wind of the donor star LM Vel could be well-established and studied in detail for the first time, as well as a very low level of accretion onto the compact object. The residual accretion rate onto the compact object in IGR J08408-4503 can be interpreted as the Bohm diffusion of (possibly magnetized) plasma entering the neutron star magnetosphere at low Bondi capture rates from the supergiant donor wind at the quasi-spherical radiation-driven settling accretion stage.

A rotation curve inequality that holds for spherically symmetric mass distributions is derived, and tested against the SPARC galaxy rotation curves dataset. We identify several Galaxies, eg NGC7793 and UGC05253, which are candidates for hosting non-spherical dark matter structures that could be detected by more precise measurements.

Gamma-ray bursts (GRBs) are ultra-relativistic collimated outflows, which emit synchrotron radiation throughout the entire electromagnetic spectrum when they interact with their environment. This afterglow emission enables us to probe the dynamics of relativistic blast waves, the microphysics of shock acceleration, and environments of GRBs. We perform Bayesian inference on a sample of GRB afterglow data sets consisting of 22 long GRBs and 4 short GRBs, using the afterglow model "scalefit", which is based on 2D relativistic hydrodynamic simulations. We make use of Gaussian processes to account for systematic deviations in the data sets, which allows us to obtain robust estimates for the model parameters. We present the inferred parameters for the sample of GRBs, and make comparisons between short GRBs and long GRBs in constant-density and stellar-wind-like environments. We find that in almost all respects such as energy and opening angle, short and long GRBs are statistically the same. Short GRBs however have a markedly lower prompt $\gamma$-ray emission efficiency than long GRBs. We also find that for long GRBs in ISM-like ambient media there is a significant anti-correlation between the fraction of thermal energy in the magnetic fields, $\epsilon_B$, and the beaming corrected kinetic energy. Furthermore, we find no evidence that the mass-loss rates of the progenitor stars are lower than those of typical Wolf-Rayet stars.

Rejecting cosmic rays (CRs) is essential for scientific interpretation of CCD-captured data, but detecting CRs in single-exposure images has remained challenging. Conventional CR-detection algorithms require tuning multiple parameters experimentally making it hard to automate across different instruments or observation requests. Recent work using deep learning to train CR-detection models has demonstrated promising results. However, instrument-specific models suffer from performance loss on images from ground-based facilities not included in the training data. In this work, we present Cosmic-CoNN, a deep-learning framework designed to produce generic CR-detection models. We build a large, diverse ground-based CR dataset leveraging thousands of images from the Las Cumbres Observatory global telescope network to produce a generic CR-detection model which achieves a 99.91% true-positive detection rate and maintains over 96.40% true-positive rates on unseen data from Gemini GMOS-N/S, with a false-positive rate of 0.01%. Apart from the open-source framework and dataset, we also build a suite of tools including console commands, a web-based application, and Python APIs to make automatic, robust CR detection widely accessible by the community of astronomers.

The scaling relation between stellar mass ($M_{*}$) and physical effective radius ($r_{e}$) has been well-studied using wide spectroscopic surveys. However, these surveys suffer from severe surface brightness incompleteness in the dwarf galaxy regime, where the relation is poorly constrained. In this study, I use a Bayesian empirical model to constrain the power-law exponent $\beta$ of the $M_{*}$-$r_{e}$ relation for late-type dwarfs ($10^{7}$$\leq$$M_{*}$/$M_{\odot}$$\leq$$10^{9}$) using a sample of 188 isolated low surface brightness (LSB) galaxies, accounting for observational incompleteness. Surprisingly, the best-fitting model ($\beta$=0.40$\pm$0.07) indicates that the relation is significantly steeper than would be expected from extrapolating canonical models into the dwarf galaxy regime. Nevertheless, the best fitting $M_{*}$-$r_{e}$ relation closely follows the distribution of known dwarf galaxies. These results indicate that extrapolated canonical models over-predict the number of large dwarf (i.e. LSB) galaxies, including ultra-diffuse galaxies (UDGs), explaining why they are over-produced by some semi-analytic models. The best-fitting model also constrains the power-law exponent of the physical size distribution of UDGs to $n\mathrm{[dex^{-1}]}\propto$$~r_{e}^{3.54\pm0.33}$, consistent to within 1$\sigma$ of the corresponding value in cluster environments and with the theoretical scenario in which UDGs occupy the high-spin tail of the normal dwarf galaxy population.

In this work we characterise the properties of the object SDSS J020536.84-081424.7, an extended nebular region with projected extension of $14 \times 14$ kpc$^{2}$ in the line of sight of the ETG Mrk 1172, using unprecedented spectroscopic data from MUSE. We perform a spatially resolved stellar population synthesis and estimate the stellar mass for both Mrk 1172 ($1 \times 10^{11} M_{\odot}$) and our object of study ($3 \times 10^{9} M_{\odot}$). While the stellar content of Mrk 1172 is dominated by an old ($\sim 10$ Gyr) stellar population, the extended nebular emission has its light dominated by young to intermediate age populations (from $\sim 100$ Myr to $\sim 1$ Gyr) and presents strong emission lines such as: H${\beta}$, [O III] ${\lambda}{\lambda}$4959,5007, H${\alpha}$, [N II] ${\lambda}{\lambda}$6549,6585 and [S II] ${\lambda}{\lambda}$6717,6732. Using these emission lines we find that it is metal-poor (with $Z \sim$ 1/3 $Z_{\odot}$, comparable to the LMC) and is actively forming stars ($0.70$ M$_{\odot}$ yr$^{-1}$), especially in a few bright clumpy knots that are readily visible in H${\alpha}$. The object has an ionised gas mass $\geq 3.8 \times 10^{5}$ M$_{\odot}$. Moreover, the motion of the gas is well described by a gas in circular orbit in the plane of a disk and is being affected by interaction with Mrk 1172. We conclude that SDSS J020536.84-081424.7 is most likely a dwarf irregular galaxy (dIGal).

The distribution of hot interstellar medium in early-type galaxies bears the imprint of the various astrophysical processes it underwent during its evolution. The X-ray observations of these galaxies have identified various structural features related to AGN and stellar feedback and environmental effects such as merging and sloshing. In our XMM-Newton Galaxy Atlas (NGA) project, we analyze archival observations of 38 ETGs, utilizing the high sensitivity and large field of view of XMM-Newton to construct spatially resolved 2D spectral maps of the hot gas halos. To illustrate our NGA data products in conjunction with the Chandra Galaxy Atlas (Kim et al. 2019), we describe two distinct galaxies - NGC 4636 and NGC 1550, in detail. We discuss their evolutionary history with a particular focus on the asymmetric distribution of metal-enriched, low-entropy gas caused by sloshing and AGN- driven uplift. We will release the NGA data products to a dedicated website, which users can download to perform further analyses.

The interaction of Galactic-Centre (GC) super bubbles (GSB) with the gaseous disc and halo of the Milky Way is investigated using radio continuum, X-ray, HI and CO line surveys. The radio North Polar Spur (NPS) constitutes the brightest eastern ridge of GSB, brightening towards the galactic plane and reaching $ l = 22\deg, \ b = + 2\deg$ at the sharpest end, where it intersects the tangential direction of the 3-kpc expanding ring and crater. Examination of the spur ridges reveals that the entire GSB, including the NPS and its counter spurs, constitutes a GC-symmetrical $\Omega /$\rotatebox[origin=c]{180}{$\Omega$} shape. The thickness and gas density of the HI and CO discs are shown to increase sharply from the inside (lower longitude) to the outside of the 3-kpc crater. Formation of crater is explained by the sweeping of the upper layer of disc gas by the shock wave from the GC by the explosion $ \sim 10 $ My ago with the emitted energy of several $10 ^ {55} $ ergs. Based on the discussion, a unified view on the structure and formation mechanism of GSB is presented.

From Swift monitoring of a sample of active galactic nuclei (AGN) we found a transient X-ray obscuration event in Seyfert-1 galaxy NGC 3227, and thus triggered our joint XMM-Newton, NuSTAR, and Hubble Space Telescope (HST) observations to study this event. Here in the first paper of our series we present the broadband continuum modelling of the spectral energy distribution (SED) for NGC 3227, extending from near infrared (NIR) to hard X-rays. We use our new spectra taken with XMM-Newton, NuSTAR, and HST/COS in 2019, together with archival unobscured XMM-Newton, NuSTAR, and HST/STIS data, in order to disentangle various spectral components of NGC 3227 and recover the underlying continuum. We find the observed NIR-optical-UV continuum is explained well by an accretion disk blackbody component (Tmax = 10 eV), which is internally reddened by E(B-V) = 0.45 with a Small Magellanic Cloud (SMC) extinction law. We derive the inner radius (12 Rg) and the accretion rate (0.1 solar mass per year) of the disk by modelling the thermal disk emission. The internal reddening in NGC 3227 is most likely associated with outflows from the dusty AGN torus. In addition, an unreddened continuum component is also evident, which likely arises from scattered radiation, associated with the extended narrow-line region (NLR) of NGC 3227. The extreme ultraviolet (EUV) continuum, and the 'soft X-ray excess', can be explained with a 'warm Comptonisation' component. The hard X-rays are consistent with a power-law and a neutral reflection component. The intrinsic bolometric luminosity of the AGN in NGC 3227 is about 2.2e+43 erg/s in 2019, corresponding to 3% Eddington luminosity. Our continuum modelling of the new triggered data of NGC 3227 requires the presence of a new obscuring gas with column density NH = 5e+22 cm^-2, partially covering the X-ray source (Cf = 0.6).

We measure the small-scale clustering of the Data Release 16 extended Baryon Oscillation Spectroscopic Survey Luminous Red Galaxy sample, corrected for fibre-collisions using Pairwise Inverse Probability weights, which give unbiased clustering measurements on all scales. We fit to the monopole and quadrupole moments and to the projected correlation function over the separation range $7-60\,h^{-1}$Mpc with a model based on the Aemulus cosmological emulator to measure the growth rate of cosmic structure, parameterized by $f\sigma_8$. We obtain a measurement of $f\sigma_8(z=0.737)=0.408\pm0.038$, which is $1.4\sigma$ lower than the value expected from 2018 Planck data for a flat $\Lambda$CDM model, and is more consistent with recent weak-lensing measurements. The level of precision achieved is 1.7 times better than more standard measurements made using only the large-scale modes of the same sample. We also fit to the data using the full range of scales $0.1-60\,h^{-1}$Mpc modelled by the Aemulus cosmological emulator and find a $4.5\sigma$ tension in the amplitude of the halo velocity field with the Planck+$\Lambda$CDM model, driven by a mismatch on the non-linear scales. This may not be cosmological in origin, and could be due to a breakdown in the Halo Occupation Distribution model used in the emulator. Finally, we perform a robust analysis of possible sources of systematics, including the effects of redshift uncertainty and incompleteness due to target selection that were not included in previous analyses fitting to clustering measurements on small scales.

We present new 3 mm observations of the ionized gas toward the nuclear starburst in the nearby (D ~ 3.5 Mpc) galaxy NGC 253. With ALMA, we detect emission from the H40-alpha and He40-alpha lines in the central 200 pc of this galaxy on spatial scales of ~4 pc. The recombination line emission primarily originates from a population of approximately a dozen embedded super star clusters in the early stages of formation. We find that emission from these clusters is characterized by electron temperatures ranging from 7000-10000 K and measure an average singly-ionized helium abundance <Y+> = 0.25 +/- 0.06, both of which are consistent with values measured for HII regions in the center of the Milky Way. We also report the discovery of unusually broad-linewidth recombination line emission originating from seven of the embedded clusters. We suggest that these clusters contribute to the launching of the large-scale hot wind observed to emanate from the central starburst. Finally, we use the measured recombination line fluxes to improve the characterization of overall embedded cluster properties, including the distribution of cluster masses and the fractional contribution of the clustered star formation to the total starburst, which we estimate is at least 50%.

Two extremely red main-belt asteroids: 203 Pompeja and 269 Justita, were identified from combined visible and near-infrared spectroscopic observations collected at the IRTF and SAO observatories. These two asteroids have a redder spectral slope than any other D-type body, which are the reddest objects in the asteroid belt, and similar to RR and IR-class objects found in the outer Solar System among trans-Neptunian objects and Centaurs. Spectroscopic results suggest the presence of complex organic materials on the surface layer of these asteroids, implying that they could have formed in the vicinity of Neptune and been transplanted to the main belt region during a phase of planetary migration. 203 Pompeia is the only very red asteroid known so far among the ~250 bodies with diameter larger than 110 km (i.e. presumably structurally intact) found in the asteroid belt. These discoveries add another piece of evidence that the main asteroid belt hosts a population of bodies that were formed in the outskirt of the Solar System.

We present the first systematic search for GHz frequency radio emission from directly imaged exoplanets using Very Large Array (VLA) observations of sufficient angular resolution to separate the planets from their host stars. We obtained results for five systems and eight exoplanets located at $\lesssim 50$ pc, through new observations (Ross 458, GU Psc, and 51 Eri) and archival data (GJ 504 and HR 8799). We do not detect radio emission from any of the exoplanets, with $3\sigma$ luminosity upper limits of $(0.9-23)\times10^{21}$ erg s$^{-1}$. These limits are comparable to the level of radio emission detected in several ultracool dwarfs, including T dwarfs, whose masses are only a factor of two times higher than those of the directly-imaged exoplanets. Despite the lack of detections in this pilot study, we highlight the need for continued GHz frequency radio observations of nearby exoplanets at $\mu$Jy-level sensitivity.

IW And stars are a recently recognized subgroup of dwarf novae which are characterized by (often repetitive) slowly rising standstills terminated by brightening, but the exact mechanism for this variation is not yet identified. We have identified BO Cet, which had been considered as a novalike cataclysmic variable, as a new member of IW And stars based on the behavior in 2019-2020. In addition to this, the object showed dwarf nova-type outbursts in 2020-2021, and superhumps having a period 7.8% longer than the orbital one developed at least during one long outburst. This object has been confirmed as an SU UMa-type dwarf nova with an exceptionally long orbital period (0.1398 d). BO Cet is thus the first cataclysmic variable showing both SU UMa-type and IW And-type features. We obtained a mass ratio (q) of 0.31-0.34 from the superhumps in the growing phase (stage A superhumps). At this q, the radius of the 3:1 resonance, responsible for tidal instability and superhumps, and the tidal truncation radius are very similar. We interpret that in some occasions this object showed IW And-type variation when the disk size was not large enough, but that the radius of the 3:1 resonance could be reached as the result of thermal instability. We also discuss that there are SU UMa-type dwarf novae above q=0.30, which is above the previously considered limit (q~0.25) derived from numerical simulations and that this is possible since the radius of the 3:1 resonance is inside the tidal truncation radius. We constrained the mass of the white dwarf larger than 1.0Msol, which may be responsible for the IW And-type behavior and the observed strength of the He II emission. The exact reason, however, why this object is unique in that it shows both SU UMa-type and IW And-type features is still unsolved.

We revisit in this work a model for repeating Fast Radio Bursts based of the release of energy provoked by the magnetic field dynamics affecting a magnetar's crust. We address the basic needs of such a model by solving the propagation approximately, and quantify the energetics and the radiation by bunches of charges in the so-called {\it charge starved} region in the magnetosphere. The (almost) simultaneous emission of newly detected X-rays from SGR 1935+2154 is tentatively associated to a reconnection behind the propagation. The strength of $f$-mode gravitational radiation excited by the event is quantified, and more detailed studies of the non-linear (spiky) soliton solutions suggested.

We performed photometric and spectroscopic investigations of NSVS 5029961 for the first time. The new BV(RI)$_c$-band light curves were obtained with the 1.0-m telescope at Weihai Observatory of Shandong University. Applying the Wilson-Devinney program, we found that NSVS 5029961 is an A-subtype shallow contact binary with extremely low mass ratio (q = 0.1515, f = 19.1\%). Six spectra have been obtained by LAMOST, and many chromospheric activity emission line indicators were detected in the spectra, revealing that the target exhibits strong chromospheric activity. We calculated the absolute parameters with the photometric solutions and Gaia distance, and estimated the initial masses of the two components and the age of the binary. The evolutionary status was discussed by using the mass-radius and mass-luminosity diagrams. The result shows the primary component is a little evolved star and the secondary component has evolved away from the main sequence. The formation and evolution investigations of NSVS 5029661 indicate that it may have evolved from a detached binary with short period and low mass ratio by angular momentum loss via magnetic braking and case A mass transfer, and is in a stable contact stage at present.

This paper continues the authors' previous work and presents a coplanar averaged ellipsoid-ellipsoid model of synchronous binary asteroid system (BAS) plus thermal and tidal effects. Using this model, we analyze the breakup mechanism of the synchronous BAS. Different from the classical spin-orbit coupling model which neglects the rotational motion's influence on the orbital motion, our model considers simultaneously the orbital motion and the rotational motions. Our findings are following. (1) Stable region of the secondary's synchronous state is mainly up to the secondary's shape. The primary's shape has little influence on it. (2) The stable region shrinks continuously with the increasing value of the secondary's shape parameter $a_B/b_B$. Beyond the value of $a_B/b_B=\sqrt{2}$, the planar stable region for the secondary's synchronous rotation is small but not zero. (3) Considering the BYORP torque, our model shows agreement with the 1-degree of freedom adiabatic invariance theory in the outwards migration process, but an obvious difference in the inwards migration process. In particular, our studies show that the so-called 'long-term' stable equilibrium between the BYORP torque and the tidal torque is never a real equilibrium state, although the binary asteroid system can be captured in this state for quite a long time. (4) In case that the primary's angular velocity gradually reduces due to the YORP effect, the secondary's synchronous state may be broken when the primary's rotational motion crosses some major spin-orbit resonances.

Previous studies have shown significant differences in the enhancement of the star-formation rate (SFR) and the star-formation efficiency (SFE=SFR/M_mol) between spiral-spiral and spiral-elliptical mergers. In order to shed light on the physical mechanism of these differences, we present NOEMA observations of the molecular gas distribution and kinematics (linear resolutions of ~ 2kpc) in two representative close major-merger star-forming pairs: the spiral-elliptical pair Arp142 and the spiral-spiral pair Arp238. The CO in Arp142 is widely distributed over a highly distorted disk without any nuclear concentration, and an off-centric ring-like structure is discovered in channel maps. The SFE varies significantly within Arp142, with a starburst region (Region 1) near the eastern tip of the distorted disk showing an SFE ~0.3 dex above the mean of the control sample of isolated galaxies, and the SFE of the main disk (Region 4) 0.43 dex lower than the mean of the control sample. In contrast, the CO emission in Arp238 is detected only in two compact sources at the galactic centers. Compared to the control sample, Arp238-E shows an SFE enhancement of more than 1 dex whereas Arp238-W has an enhancement of ~0.7 dex. We suggest that the extended CO distribution and the large SFE variation in Arp142 are due to an expanding large-scale ring triggered by a recent high-speed head-on collision between the spiral galaxy and the elliptical galaxy, and the compact CO sources with high SFEs in Arp238 are associated with nuclear starbursts induced by gravitational tidal torques in a low-speed coplanar interaction.

Globular Clusters (GCs) are now well known to almost universally show multiple popu-lations (MPs). The HST UV Legacy Survey of a large number of Galactic GCs in UV filters optimized to explore MPs finds that a small fraction of GCs, termed Type II, also display more complex, anomalous behavior. Several well-studied Type II GCs show intrinsic Fe abundance variations, suggesting that the other, less well-studied, Type II GCs should also exhibit similar behavior. Our aim is to perform the first detailed metallicity analysis of NGC 1261, an intermediate mass Type II GC, in order to determine if this object shows an intrinsic Fe variation. We determined the Fe abundance in eight red giant members using Magellan-MIKE and UVES-FLAMES high-resolution, high S/N spectroscopy. The full range of [Fe/H] for the entire sample from the spectra is from -1.05 to -1.43 dexwith an observed spread sigma_obs=0.133 dex. Compared with the total internal error of Sigma_tot=0.06,this indicates a significant intrinsic metallicity spread of Sigma_int=0.119 dex. We found a very similar variation in [Fe/H] using an independent method to derive the atmospheric parameters based on near-IR photometry. More importantly, the mean metallicity of the five presumed normal metallicity stars is -1.37+/-0.02, while that of the three presumed anomalous/highmetallicity stars is -1.18+/-0.09. This difference is significant at the $\pm$2.4Sigma level. We find indications from existing data of other Type II GCs that several of them presumedto have real metallicity spreads may in fact posses none. The minimum mass required for a GC to acquire an intrinsic Fe spread appears to be $\pm$10^5 Msun. We find no strong correlation betwee nmass and metallicity variation for Type II GCs. The metallicity spread is also independent of the fraction of anomalous stars within the Type II GCs and of GC origin.

Hot massive stars present strong stellar winds that are driven by absorption, scattering and re\-emission of photons by the ions of the atmosphere (\textit{line-driven winds}). A better comprehension of this phenomenon, and a more accurate calculation of hydrodynamics and radiative acceleration is required to reduce the number of free parameters in spectral fitting, to determine accurate wind parameters such as mass-loss rates and velocity profiles. We use the non-LTE model-atmosphere code CMFGEN to numerically solve the radiative transfer equation in the stellar atmosphere and to calculate the radiative acceleration $g_\text{rad}(r)$. Under the assumption that the radiative acceleration depends only on the radial coordinate, we solve analytically the equation of motion by means of the Lambert $W$-function. An iterative procedure between the solution of the radiative transfer and the equation of motion is executed in order to obtain a final self-consistent velocity field that is no longer based on any $\beta$-law. We apply the Lambert-procedure to three O supergiant stars ($\zeta$-Puppis, HD~165763 and $\alpha$-Cam) and discuss the Lambert-solutions for the velocity profiles. It is found that, even without recalculation of the mass-loss rate, the Lambert-procedure allows the calculation of consistent velocity profiles that reduce the number of free parameters when a spectral fitting using CMFGEN is performed. Synthetic spectra calculated from our Lambert-solutions show significant differences compared to the initial $\beta$-law CMFGEN models. The results indicate the importance of consistent velocity profile calculation in the CMFGEN code and its usage in a fitting procedure and interpretation of observed spectra.

The majority of observed pixels on the Transiting Exoplanet Survey Satellite (TESS) are delivered in the form of full frame images (FFI). However, the FFIs contain systematic effects such as pointing jitter and scattered light from the Earth and Moon that must be removed before downstream analysis. We present unpopular, an open-source Python package to de-trend TESS FFI light curves based on the causal pixel model method. Under the assumption that shared flux variations across multiple distant pixels are likely to be systematics, unpopular removes these common (i.e., popular) trends by modeling the systematics in a given pixel's light curve as a linear combination of light curves from many other distant pixels. To prevent overfitting we employ ridge regression and a train-and-test framework where the data points being de-trended are separated from those used to obtain the model coefficients. We also allow for simultaneous fitting with a polynomial model to capture any long-term astrophysical trends. We validate our method by de-trending different sources (e.g., supernova, tidal disruption event, exoplanet-hosting star, fast rotating star) and comparing our light curves to those obtained by other pipelines when appropriate. We also show that unpopular is able to preserve sector-length astrophysical signals, allowing for the extraction of multi-sector light curves from the FFI data. The unpopular source code and tutorials are freely available online.

We present optical follow-up observations with the DDOTI telescope of gravitational-wave events detected during the Advanced LIGO and Advanced Virgo O3 observing run. DDOTI is capable of responding to an alert in a few minutes, has an instantaneous field of about 69 deg$^{2}$, and obtains $10\sigma$ upper limits of $w_{\rm lim}=18.5$ to 20.5 AB mag in 1000~s of exposure, depending on the conditions. We observed 54\% (26 out of 48) of the unretracted gravitational-wave alerts and did not find any electromagnetic counterparts. We compare our upper limits to various possible counterparts: the kilonova AT~2017gfo, models of radioactive- and magnetar-powered kilonovae, short gamma-ray burst afterglows, and AGN flares. Although the large positional uncertainties of GW sources do not allow us to place strong constraints during O3, DDOTI observations of well-localized GW events in O4 and beyond could meaningfully constrain models of compact binary mergers. We show that DDOTI is able to detect kilonovae similar to AT~2017gfo up to about 200~Mpc and magnetar-powered kilonovae up to 1~Gpc. We calculate that nearby ($\lesssim$200 Mpc) afterglows have a high chance ($\approx$70\%) to be detected by rapid ($\lesssim$3 hours) DDOTI observations if observed on-axis, whereas off-axis afterglows are unlikely to be seen. Finally, we suggest that long-term monitoring of massive BBH events with DDOTI could confirm or rule out late AGN flares associated with these events.

We perform a zero-$\beta$ magnetohydrodynamic simulation for the C7.7 class flare initiated at 01:18 UT on 2011 June 21 using the Message Passing Interface Adaptive Mesh Refinement Versatile Advection Code (MPI-AMRVAC). The initial condition for the simulation involves a flux rope which we realize through the regularized Biot-Savart laws, whose parameters are constrained by observations from the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO) and the Extreme Ultraviolet Imager (EUVI) on the twin Solar Terrestrial Relations Observatory (STEREO). This data-constrained initial state is then relaxed to a force-free state by the magneto-frictional module in MPI-AMRVAC. The further time-evolving simulation results reproduce the eruption characteristics obtained by SDO/AIA 94 A, 304 A, and STEREO/EUVI 304 A observations fairly well. The simulated flux rope possesses similar eruption direction, height range, and velocities to the observations. Especially, the two phases of slow evolution and fast eruption are reproduced by varying the density distribution in light of the filament material draining process. Our data-constrained simulations also show other advantages, such as a large field of view (about 0.76 solar radii). We study the twist of the magnetic flux rope and the decay index of the overlying field, and find that in this event, both the magnetic strapping force and the magnetic tension force are sufficiently weaker than the magnetic hoop force, thus allowing the successful eruption of the flux rope. We also find that the anomalous resistivity is necessary in keeping the correct morphology of the erupting flux rope.

Understanding the cluster population of clusters of galaxies is of the utmost importance for using cluster samples in both astrophysical and cosmological studies. We present an in-depth analysis of the X-ray morphological parameters of the galaxy clusters and groups detected in the eROSITA Final Equatorial-Depth Survey (eFEDS). We study the eROSITA X-ray imaging data for a sample of 325 clusters and groups that are significantly detected in the eFEDS field. We characterize their dynamical properties by measuring a number of dynamical estimators: concentration, central density, cuspiness, ellipticity, power-ratios, photon asymmetry, and Gini coefficient. The galaxy clusters and groups detected in eFEDS, covering a luminosity range of more than three orders of magnitude and large redshift range out to 1.2 provide an ideal sample for studying the redshift and luminosity evolution of the morphological parameters and characterization of the underlying dynamical state of the sample. Based on these measurements we construct a new dynamical indicator, relaxation score, for all the clusters in the sample. We find no evidence for bimodality in the distribution of morphological parameters of our clusters, rather we observe a smooth transition from the cool-core to non-cool-core and from relaxed to disturbed states. A significant evolution in redshift and luminosity is also observed in the morphological parameters examined in this study after carefully taking into account the selection effects. We determine that our eFEDS-selected cluster sample, differently than ROSAT-based cluster samples, is not biased toward relaxed clusters, but contains a similar fraction of disturbed as SZ surveys.

Using TESS we are doing a systematic study of outbursting AM~CVn systems to place some limits on the current outbursts models. We present the TESS light curve (LC) for 9 AM~CVns showing both superoutbursts (SO) and normal outbursts (NO). The continuous coverage of the outbursts with TESS allows us to place stringent limits on the duration and structures of the SO and the NO. We present evidence that in at least some of the systems enhanced mass transfer (EMT) has to be taken into account to explain the observed LC of the SO and rebrighthening phase after the SO. For others, the colour evolution from simultaneous observations in $g$ and $r$ with ZTF differs from previously reported color evolution of longer period AM~CVns where EMT is responsible for the SO. We also find that due to the lack of sufficiently high cadence coverage many of the duration might have been overestimated in previous ground-based surveys and report the SO duration for 6 AM~CVns. We also found that precursors are a common feature of SO in AM~CVns and are seen in the LC of 5 of the 6 reported SO. Finally with the 10-minute and 2-minute cadence LC from TESS also allowed us to find two new candidates orbital periods of AM~CVns, both of which are in reasonably good agreement with the predictions for their periods based on their past outburst histories.

Young open clusters are ideal laboratories to understand star formation process. We present deep UBV I and Halpha photometry for the young open cluster IC 1590 in the center of the H II region NGC 281. Early-type members are selected from UBV photometric diagrams, and low-mass pre-main sequence (PMS) members are identified by using Halpha photometry. In addition, the published X-ray source list and Gaia astrometric data are also used to isolate probable members. A total of 408 stars are selected as members. The mean reddening obtained from early-type members is <E(B-V) = 0.40 +/- 0.06 (s.d.). We confirm the abnormal extinction law for the intracluster medium. The distance modulus to the cluster determined from the zero-age main-sequence fitting method is 12.3 +/- 0.2 mag (d = 2.88 +/- 0.28 kpc), which is consistent with the distance d = 2.70 ^+0.24 _-0.20 kpc from the recent Gaia parallaxes. We also estimate the ages and masses of individual members by means of stellar evolutionary models. The mode of the age of PMS stars is about 0.8 Myr. The initial mass function of IC 1590 is derived. It appears a steeper shape (Gamma = -1.49 +/- 0.14) than that of the Salpeter/Kroupa initial mass function for the high mass regime (m > 1 M_sun). The signature of mass segregation is detected from the difference in the slopes of the initial mass functions for the inner (r < 2.'5) and outer region of this cluster. We finally discuss the star formation history in NGC 281.

We present a fully nonlinear and relativistically covariant expression for the observed galaxy density contrast. Building on a null tetrad tailored to the cosmological observer's past light cone, we find a decomposition of the nonlinear galaxy over-density into manifestly gauge-invariant quantities, each of which has a clear physical interpretation as a cosmological observable. This ensures that the monopole of the galaxy over-density field is properly accounted for. We anticipate that this decomposition will be useful for future work on nonlinearities in galaxy number counts, for example, deriving the relativistic expression for the galaxy bispectrum. We then specialise our results to conformal Newtonian gauge, with a Hubble parameter either defined globally or measured locally, illustrating the significance of the different contributions to the observed monopole of the galaxy density.

For the first time, simultaneous broadband spectral and timing study of the atoll source 4U 1705-44 was performed using AstroSat Soft X-ray Telescope (SXT) and Large Area X-ray Proportional Counter (LAXPC) data (0.8-70 keV). Based on the HID, the source was in the soft banana state during these observations. Spectral modeling was performed using the full reflection framework and an inner disk radii of 14 Rg was obtained. A hard powerlaw tail was noticed in the soft state and hot component fluxes and varying powerlaw indices point towards a varying corona/sub-keplerian flow. Based on the spectral fits the boundary layer radius and magnetospheric radius were constrained to be $\sim$ 14-18 km and $\sim$ 9-19 km respectively. Cross- Correlation Function studies were performed between the 0.8-3 keV soft SXT lightcurve and 10-20 keV hard LAXPC lightcurve and correlated and anticorrelated lags were found, which was used to constrain the coronal height to 0.6-20 km (\b{eta}=0.1). Since the inner disk radius is not varying during the observations, we conclude that the detected lags are possibly caused by a varying structure of corona/boundary layer in the inner region of the accretion disk. Based on the observations, a geometrical model is proposed for explaining the detected lags in the atoll source 4U 1705-44.

The invariance of the speed of light in the distant universe has profound significance for fundamental physics. In this paper, we propose a new model-independent method to test the invariance of the speed of light $c$ at different redshifts by combining the strong gravitational lensing (SGL) systems and the observations of type-Ia supernovae (SNe Ia). All the quantities used to test the deviation of $c$ come from the direct observations, and the absolute magnitudes of SNe Ia need not to be calibrated. Our results show that the speed of light in the distant universe is no obvious deviation from the constant value $c_0$ within the uncertainty based on current observations. Moreover, we conclude that the currently compiled SGL and SNe Ia Pantheon samples may achieve much higher precision $\Delta c/c\sim10^{-2}$ for the deviation of $c$ than all previously considered approaches. The forthcoming data from the Legacy Survey of Space and Time and Wide-Field InfraRed Space Telescope will achieve more stringent testing for deviation of the SOL (at the level of $\Delta c/c \sim10^{-3}$) by using our model-independent method. Finally, we discuss the potential ways in which our technique might be improved, focusing on the treatment of possible sources of systematic uncertainties.

Various solutions of the kinetic equation for the equilibrium of a gravitating sphere of uniform density with a quadratic gravitational potential and a linear dependence of gravitational force on radius are examined. New analytic solutions are obtained for a uniform sphere with a hollow spherical volume and central mass inside the sphere. Solutions are also obtained for an arbitrary number of spherical layers with the same density, but with different equilibrium distribution functions.

The search for extraterrestrial intelligence (SETI) has been conducted for nearly 60 years. A Dyson Sphere, a spherical structure that surrounds a star and transports its radiative energy outward as an energy source for an advanced civilisation, is one of the main targets of SETI. In this study, we discuss whether building a Dyson Sphere around a black hole is effective. We consider six energy sources: (i) the cosmic microwave background, (ii) the Hawking radiation, (iii) an accretion disk, (iv) Bondi accretion, (v) a corona, and (vi) relativistic jets. To develop future civilisations (for example, a Type II civilisation), $4\times10^{26}\,{\rm W}$($1\,{\rm L_{\odot}}$) is expected to be needed. Among (iii) to (vi), the largest luminosity can be collected from an accretion disk, reaching $10^{5}\,{\rm L_{\odot}}$, enough to maintain a Type II civilisation. Moreover, if a Dyson Sphere collects not only the electromagnetic radiation but also other types of energy (e.g., kinetic energy) from the jets, the total collected energy would be approximately 5 times larger. Considering the emission from a Dyson Sphere, our results show that the Dyson Sphere around a stellar-mass black hole in the Milky Way ($10\,\rm kpc$ away from us) is detectable in the ultraviolet$(\rm 10-400\,{\rm nm)}$, optical$(\rm 400-760\,{\rm nm)}$, near-infrared($\rm 760\,{\rm nm}-5\,{\rm \mu m}$), and mid-infrared($\rm 5-40\,{\rm \mu m}$) wavelengths via the waste heat radiation using current telescopes such as Galaxy Evolution Explorer Ultraviolet Sky Surveys. Performing model fitting to observed spectral energy distributions and measuring the variability of radial velocity may help us to identify these possible artificial structures.

In the effort to define the main features a successful solution to the Hubble tension should have, growing evidence has emerged pointing to the need for models able to modify the expansion (and possibly thermal) history of the universe both prior and after recombination, and that could thereby restore the overall concordance between early- and late-time observations without introducing nor worsening other cosmological tensions in the process. In precisely this spirit, here we consider a multi-interacting dark energy model with two complementary interaction channels: one with the dark matter and one with the photons. The former most significantly affects the matter dominated epoch as well as the very late universe, and has been shown to be able to significantly resolve the $S_8$ tension. The latter has been introduced in this work to allow for a larger and natural time dependence of the evolution of the universe since it mostly impacts the radiation dominated epoch as well as the temperature scaling of the photons, as extensively explained in the text. As a result, considering data from Planck+BAO+Pantheon+KV450+DES (which can be combined since the $S_8$ value is lowered even neglecting the weak lensing data) we find that the significance of the $H_0$ tension only reduces to about 3.5$\sigma$, while that of the $S_8$ tension falls below the 2$\sigma$ level. Also, the statistical analysis we perform strongly favors $\Lambda$CDM in all considered cases. Overall, we conclude that the specific multi-interacting dark energy model considered here, despite its broad generality and very rich cosmological phenomenology, cannot successfully restore the aforementioned overall concordance between early- and late-time observations.

We present archival Giant Metrewave Radio Telescope (GMRT) observations of two exoplanetary systems, $\tau$ Bo\"otis, and 55 Cancri, at 610 MHz and 150 MHz, respectively. Theoretical models predict these systems to have some of the highest expected flux densities at radio wavelengths. Both $\tau$ Bo\"otis and 55 Cancri have been previously observed at low frequency ($\sim$ 30 MHz) with Low-Frequency Array (LOFAR) (Turner et al. 2020). $\tau$ Bo\"otis shows tentative signatures of circularly polarized emission at 30 MHz, while no emission was detected from 55 Cancri. We do not detect radio emission from both the systems, but the GMRT observations set $3\sigma$ upper limits of 0.6 mJy at 610 MHz for $\tau$ Bo\"otis and 4.6 mJy at 150 MHz for 55 Cancri. The sensitivity achieved at 610 MHz in these observations is comparable to some of the deepest images of an exoplanet field.

We present the prospects for the pre-merger detection and localization of binary neutron star mergers with third generation gravitational-wave observatories. We consider a wide variety of gravitational-wave networks which may be operating in the 2030's and beyond; these networks include up to two Cosmic Explorer sites, the Einstein Telescope, and continued observation with the existing second generation ground-based detectors. For a fiducial merger rate of 300 Gpc$^{-3}$yr$^{-1}$, we find that the Einstein Telescope on its own is able to detect 6 (2) sources per year 5 (30) minutes before merger and provide a localization of $<10~\textrm{deg}^2$. A single Cosmic Explorer would detect but be unable to localize sources on its own. A two-detector Cosmic Explorer network, however, would detect 22 (0.4) mergers per year using the same criteria. A full three-detector network with the operation of dual Cosmic Explorers and the Einstein Telescope would allow for $<1~\textrm{deg}^2$ source localization at 5 minutes before merger for $\sim7$ sources per year. Given the dramatic increase in localization and detection capabilities, third generation observatories will enable the regular observation of the prompt emission of mergers by a broad array of observatories including gamma-ray, x-ray, and optical telescopes. Moreover, sub-degree localizations minutes before merger, combined with narrow-field-of-view high-energy telescopes, could strongly constrain the high-energy pre-merger emission models proposed in the last decade.

We present here up-to-date neutrino mass limits exploiting the most recent cosmological data sets. By making use of the Cosmic Microwave Background temperature fluctuation and polarization measurements, Supernovae Ia luminosity distances, Baryon Acoustic Oscillation observations and determinations of the growth rate parameter, we are able to set the most constraining bound to date, $\sum m_\nu<0.09$~eV at $95\%$~CL. This very tight limit is obtained without the assumption of any prior on the value of the Hubble constant and highly compromises the viability of the inverted mass ordering as the underlying neutrino mass pattern in nature. The results obtained here further strengthen the case for very large multitracer spectroscopic surveys as unique laboratories for cosmological relics, such as neutrinos: that would be the case of the Dark Energy Spectroscopic Instrument (DESI) survey and of the Euclid mission.

We use Gemini Near-Infrared Integral Field Spectrograph (NIFS) observations of the inner 3.5$\times$3.5 kpc$^2$ of the radio galaxy Cygnus A to map the gas excitation and kinematics at a spatial resolution of 200 pc. The emission of the ionised gas shows a biconical morphology, with half-opening angle of 45$^\circ$ and oriented along the position angle of the radio jet. Coronal line emission is seen within the cone, up to 1.75 kpc from the nucleus, with higher ionisation gas observed in the easterly side. The H$_2$ and [Fe II] emission lines are consistent with excitation by the central AGN, with some contribution of shocks to the southwest of the nucleus. The gas visual extinction and electron density are larger than those from optical-based measurements, consistent with the fact that near-IR observations penetrate deeply into the gas emission structure, probing denser and more obscured regions. The gas kinematics shows two components: (i) a rotating disc with kinematic position angle of $\Psi_0=21^\circ\pm2^\circ$, seen both in ionised and molecular gas, and (ii) outflows with velocities of up to 600 kms$^{-1}$ observed within the ionisation cone in ionised gas and restricted to inner 0.5 arcsec in molecular gas. The mass outflow rate in ionised gas is in the range $\sim100-280 {\rm M_\odot yr^{-1}}$ and the kinetic power of the outflow corresponds to 0.3-3.3 per cent of the AGN bolometric luminosity, indicating that the outflows in Cygnus A may be effective in suppressing star formation.

The streaming instability provides an efficient way of overcoming the growth barriers in the initial stages of the planet formation process. Considering the realistic case of a particle size distribution, the dynamics of the system is altered compared to the outcome of single size models. In order to understand the outcome of the multi-species streaming instability in detail, we perform a large parameter study in terms of particle number, particle size distribution, particle size range, initial metallicity and initial particle scale height. We study vertically stratified systems and determine the metallicity threshold for filament formation. We compare these with a system where the initial particle distribution is unstratified and find that its evolution follows that of its stratified counterpart. We find that change in particle number does not result in significant variation in the efficiency and timing of filament formation. We also see that there is no clear trend for how varying the size distribution in combination with particle size range changes the outcome of the multi-species streaming instability. Finally, we find that an initial metallicity of $Z_{\rm{init}}=0.005$ and $Z_{\rm{init}}=0.01$ both result in similar critical metallicity values for the start of filament formation. Our results show that the inclusion of a particle size distribution into streaming instability simulations, while changing the dynamics as compared to mono-disperse systems, does not result in overall unfavorable conditions for solid growth. We attribute the sub-dominant role of multiple species to the high-density conditions in the midplane, conditions under which also linear stability analysis predict little difference between single and multiple species.

There is a wide consensus that type Ia supernovae (SN Ia) originate from the thermonuclear explosion of CO white dwarfs (WD), with the lack of hydrogen in the observed spectra as a distinctive feature. Here, we present SN 2016jae, which was classified as a Type Ia SN from a spectrum obtained soon after the discovery. The SN reached a B-band peak of -17.93 +- 0.34 mag, followed by a fast luminosity decline with sBV 0.56 +- 0.06 and inferred Dm15(B) of 1.88 +- 0.10 mag. Overall, the SN appears as a "transitional" event between "normal" SNe Ia and very dim SNe Ia such as 91bg-like SNe. Its peculiarity is that two late-time spectra taken at +84 and +142 days after the peak show a narrow line of Halpha (with full width at half-maximum of ~650 and 1000 kms-1, respectively). This is the third low-luminosity and fast-declining Type Ia SN after SN 2018cqj/ATLAS18qtd and SN 2018fhw/ASASSN-18tb, found in the 100IAS survey that shows resolved narrow Halpha line in emission in their nebular-phase spectra. We argue that the nebular Halpha emission originates in an expanding hydrogen-rich shell (with velocity < 1000 kms-1). The hydrogen shell velocity is too high to be produced during a common envelope phase, while it may be consistent with some material stripped from an H-rich companion star in a single-degenerate progenitor system. However, the derived mass of this stripped hydrogen is ~0.002-0.003 Msun, which is much less than that expected (>0.1 Msun) for standard models for these scenarios. Another plausible sequence of events is a weak SN ejecta interaction with a H-shell ejected by optically thick winds or a nova-like eruption on the C/O WD progenitor some years before the supernova explosion.

The mode switching in Cepheids is studied using the methods of the nonlinear theory of stellar pulsation, depending on the main sequence mass $M_0$ and the abundance of elements heavier than helium $Z$. The grid of evolutionary and hydrodynamic models of core--helium burning Cepheids is represented by 30 evolutionary sequences of stars with initial masses $5.7M_\odot\le M_0\le 7.2M_\odot$ and $Z=0.014$, 0.018, 0.022. For considered values of $Z$ the periods of the fundamental mode and the first overtone at the oscillation mode switching are shown to depend on the mean density of the stellar matter. The upper limit of the period of the first overtone decreases witn increasing $Z$ from $\approx 6.9$ day for $Z=0.014$ to $\approx 4.1$ day for $Z=0.022$. The theoretical period--radius relation is independent of $Z$ and agrees well (within 2.5\%) with recent measurements of Cepheid radii based on the Baade--Wesselink method. The fundamental parameters of the short--period Cepheid CG Cas were derived with application of observational estimates of the period and the rate of period change. This star is shown to be the first--overtone pulsator.

Ultra-hot Jupiters (UHJs) present excellent targets for atmospheric characterisation. Their hot dayside temperatures (T $\gtrsim$ 2200 K) strongly suppress the formation of condensates, leading to clear and highly-inflated atmospheres extremely conducive to transmission spectroscopy. Recent studies using optical high-resolution spectra have discovered a plethora of neutral and ionised atomic species in UHJs, placing constraints on their atmospheric structure and composition. Our recent work has presented a search for molecular features and detection of Fe I in the UHJ WASP-121b using VLT/UVES transmission spectroscopy. Here, we present a systematic search for atomic species in its atmosphere using cross-correlation methods. In a single transit, we uncover potential signals of 17 atomic species which we investigate further, categorising 5 as strong detections, 3 as tentative detections, and 9 as weak signals worthy of further exploration. We confirm previous detections of Cr I, V I, Ca I, K I and exospheric H I and Ca II made with HARPS and ESPRESSO, and independently re-recover our previous detection of Fe I at 8.8 $\sigma$ using both the blue and red arms of the UVES data. We also add a novel detection of Sc II at 4.2 $\sigma$. Our results further demonstrate the richness of UHJs for optical high-resolution spectroscopy.

The Gaia space astrometry mission is measuring accurate distances and space motions of more than two billion stars throughout our Galaxy and beyond. This is a first look at how Gaia is contributing to fundamental physics, and in particular to our understanding of dark matter, for which a few examples are given from the current literature. One of our goals is to illustrate how deep, and often surprising, insight into very diverse areas of fundamental physics can be extracted from this new and enormous high-accuracy stellar data set. In this spirit, we finish by suggesting a search for a connection between stellar activity, dark matter streams, and planetary configuration in nearby exoplanetary systems, as has been tentatively proposed in the case of the solar system.

We present a detailed analysis of an eruptive event that occurred on early 2019 March 8 in active region AR 12734, to which we refer as the International Women's day event. The event under study is intriguing in several aspects: 1) low-coronal eruptive signatures come in ''pairs'' (a double-peak flare, two coronal dimmings, and two EUV waves); 2) although the event is characterized by a complete chain of eruptive signatures, the corresponding coronagraphic signatures are weak; 3) although the source region of the eruption is located close to the center of the solar disc and the eruption is thus presumably Earth-directed, heliospheric signatures are very weak with little Earth-impact. We analyze a number of multi-spacecraft and multi-instrument (both remote-sensing and in situ) observations, including Soft X-ray, (extreme-) ultraviolet (E)UV), radio and white-light emission, as well as plasma, magnetic field and particle measurements. We employ 3D NLFF modeling to investigate the coronal magnetic field configuration in and around the active region, the GCS model to make a 3D reconstruction of the CME geometry and the 3D MHD numerical model EUHFORIA to model the background state of the heliosphere. Our results indicate two subsequent eruptions of two systems of sheared and twisted magnetic fields, which merge already in the upper corona and start to evolve further out as a single entity. The large-scale magnetic field significantly influences both, the early and the interplanetary evolution of the structure. During the first eruption the stability of the overlying field was disrupted which enabled the second eruption. We find that during the propagation in the interplanetary space the large-scale magnetic field, i.e. , the location of heliospheric current sheet between the AR and the Earth likely influences propagation and the evolution of the erupted structure(s).

Our aim was to determine the initial Li content of two clusters of similar metallicity but very different ages, the old open cluster NGC 2243 and the metal-rich globular cluster NGC 104. We compared the lithium abundances derived for a large sample of stars (from the turn-off to the red giant branch) in each cluster. For NGC 2243 the Li abundances are from the catalogues released by the Gaia-ESO Public Spectroscopic Survey, while for NGC 104 we measured the Li abundance using FLAMES/GIRAFFE spectra, which include archival data and new observations. We took the initial Li of NGC 2243 to be the lithium measured in stars on the hot side of the Li dip. We used the difference between the initial abundances and the post first dredge-up Li values of NGC 2243, and by adding this amount to the post first dredge-up stars of NGC~104 we were able to infer the initial Li of this cluster. Moreover, we compared our observational results to the predictions of theoretical stellar models for the difference between the initial Li abundance and that after the first dredge-up. The initial lithium content of NGC 2243 was found to be A(Li)_i = 2.85dex by taking the average Li abundance measured from the five hottest stars with the highest lithium abundance. This value is 1.69 dex higher than the lithium abundance derived in post first dredge-up stars. By adding this number to the lithium abundance derived in the post first dredge-up stars in NGC~104, we infer a lower limit of its initial lithium content of A(Li)_i= 2.30dex. Stellar models predict similar values. Therefore, our result offers important insights for further theoretical developments.

In recent years, thanks to new facilities such as LOFAR capable of sensitive observations, much work has been done on the detection of stellar radio emission at low frequencies. Such emission has commonly been shown to be coherent emission, generally attributed to electron-cyclotron maser emission, and has usually been detected from main-sequence M dwarfs. Here we report the first detection of coherent emission at low frequencies from T Tauri stars, which are known to be associated with high levels of stellar activity. Using LOFAR, we have detected several bright radio bursts at 150 MHz from two weak-line T Tauri stars: KPNO-Tau 14 and LkCa 4. All of the bursts have high brightness temperatures ($10^{13} - 10^{14}\ \mathrm{K}$) and high circular polarization fractions (60 - 90 \%), indicating that they must be due to a coherent emission mechanism. This could be either plasma emission or electron-cyclotron maser (ECM) emission. Due to the exceptionally high brightness temperatures seen in at least one of the bursts ($\geq 10^{14}\ \mathrm{K}$), and the high circular polarization levels, it seems unlikely that plasma emission could be the source and so ECM is favoured as the most likely emission mechanism. Assuming this is the case, the required magnetic field in the emission regions would be 40 - 70 G. We determine that the most likely method of generating ECM emission is plasma co-rotation breakdown in the stellar magnetosphere. There remains the possibility, however, it could be due to an interaction with an orbiting exoplanet.

The tidal evolution of interacting binaries when the orbital period is short compared to the primary star's convective time scale is a problem of long-standing. Terquem (2021) has argued that, when this temporal ordering scheme is obeyed, the rate of energy transfer from tides to convection (denoted $D_R$) is given by the product of the averaged Reynolds stress associated with the tidal velocity and the mean shear associated with the convective flow. In a recent response, Barker and Astoul (2021, hereafter BA21) claim to show that $D_R$ (in this form) cannot contribute to tidal dissipation. Their analysis is based on a study of Boussinesq and anelastic models. Here, we demonstrate that BA21 misidentify the correct term responsible for energy transfer between tides and convection. As a consequence, their anelastic calculations do not prove that the $D_R$ formulation is invalidated as an energy-loss coupling between tides and convection. BA21 also carry out a calculation in the Boussinesq approximation. Here, their claim that $D_R$ once again does not contribute is based on boundary conditions that do not apply to any star or planet that radiates energy from its surface, which is a key dissipational process in the problem we consider.

We propose two model independent methods to obtain constraints on the transition and equivalence redshifts $z_{tr}$, $z_{eq}$. In particular, we consider $z_{tr}$ as the onset of cosmic acceleration, whereas $z_{eq}$ the redshift at which the densities of dark energy and pressureless matter are equated. With this prescription, we expand the Hubble and deceleration parameters up to two hierarchical orders and show a linear correlation between transition and equivalence, from which we propose exclusion plots where $z_{eq}$ is not allowed to span. To this end, we discuss how to build up cosmographic expansions in terms of $z_{tr}$ and compute the corresponding observable quantities directly fitting the luminosity and angular distances and the Hubble rate with cosmic data. We make our computations through Monte Carlo fits involving type Ia supernova, baryonic acoustic oscillation and Hubble most recent data catalogs. We show at $1\sigma$ confidence level the $\Lambda$CDM predictions on $z_{tr}$ and $z_{eq}$ are slightly confirmed, although at $2\sigma$ confidence level dark energy expectations cannot be excluded. Finally, we theoretically interpret our outcomes and discuss possible limitations of our overall approach.

Galactic cosmic-rays (GCRs) are thought to be accelerated in strong shocks induced by massive star winds and supernova explosions sweeping across the interstellar medium. But the phase of the interstellar medium from which the CRs are extracted has remained elusive up to now. Here, we study in detail the GCR source composition deduced from recent measurements by the AMS-02, Voyager 1 and SuperTIGER experiments to obtain information on the composition, ionisation state and dust content of the GCR source reservoirs. We show that the volatile elements of the CR material are mainly accelerated from a plasma of temperature above about 2 MK, which is typical of the hot medium found in galactic superbubbles energised by the activity of massive star winds and supernova explosions. Another GCR component, which is responsible for the overabundance of $^{22}$Ne, most likely arises from acceleration of massive star winds in their termination shocks. From the CR-related gamma-ray luminosity of the Milky Way, we estimate that the ion acceleration efficiency in both supernova shocks and wind termination shocks is of the order of $10^{-5}$. The GCR source composition also shows evidence for a preferential acceleration of refractory elements contained in interstellar dust. We suggest that the GCR refractories are also produced in superbubbles, from shock acceleration and subsequent sputtering of dust grains continuously incorporated into the hot plasma through thermal evaporation of embedded molecular clouds. Our model explains well the measured abundances of all primary and mostly primary CRs from H to Zr, including the overabundance of $^{22}$Ne.

We present an optical study of the strong lensing galaxy cluster MS 0440.5$+$0204 at $z=0.19593$, based on CFHT/MegaCam g', r'-photometry and GMOS/Gemini and CFHT/MOS/SIS spectroscopy in a broader area compared to previous works. We have determined new spectroscopic redshifts for the most prominent gravitational arcs surrounding the central galaxy in the cluster. The new redshifts and the information provided by the photometric catalog yield us to perform a detailed weak and strong lensing mass reconstruction of the cluster. The large number of member galaxies and the area covered by our observations allow to estimate more accurately the velocity dispersion and mass of cluster and examine in detail the nature of the cluster and surroundings structures. The dynamical mass is in good agreement with the mass inferred from the lensing analysis and X-ray estimates. About $\sim$68\% of the galaxies are located in the inner $\lesssim$0.86 h$^{-1}_{70}$ Mpc region of the cluster. The galaxy redshift distribution in the inner region of the cluster shows a complex structure with at least three sub-structures along the line-of-sight. Other sub-structures are also identified in the galaxy density map and in the weak lensing mass map. The member galaxies in the North-East overdensity are distributed in a filament between MS 0440.5$+$0204 and ZwCL 0441.1$+$0211 clusters, suggesting that these two structures might be connected. MS 0440$+$0204 appears to be dynamically active, with a cluster core that is likely experiencing a merging process and with other nearby groups at projected distances of $\lesssim$1 h$^{-1}_{70}$ Mpc that could be being accreted by the cluster.

The Multi-slit Solar Explorer (MUSE) is a proposed NASA MIDEX mission, currently in Phase A, composed of a multi-slit EUV spectrograph (in three narrow spectral bands centered around 171A, 284A, and 108A) and an EUV context imager (in two narrow passbands around 195A and 304A). MUSE will provide unprecedented spectral and imaging diagnostics of the solar corona at high spatial (<0.5 arcsec), and temporal resolution (down to ~0.5s) thanks to its innovative multi-slit design. By obtaining spectra in 4 bright EUV lines (Fe IX 171A , Fe XV 284A, Fe XIX-Fe XXI 108A) covering a wide range of transition region and coronal temperatures along 37 slits simultaneously, MUSE will for the first time be able to "freeze" (at a cadence as short as 10 seconds) with a spectroscopic raster the evolution of the dynamic coronal plasma over a wide range of scales: from the spatial scales on which energy is released (~0.5 arcsec) to the large-scale often active-region size (170 arcsec x 170 arcsec) atmospheric response. We use advanced numerical modeling to showcase how MUSE will constrain the properties of the solar atmosphere on the spatio-temporal scales (~0.5 arcsec, ~20 seconds) and large field-of-view on which various state-of-the-art models of the physical processes that drive coronal heating, solar flares and coronal mass ejections (CMEs) make distinguishing and testable predictions. We describe how the synergy between MUSE, the single-slit, high-resolution Solar-C EUVST spectrograph, and ground-based observatories (DKIST and others) can address how the solar atmosphere is energized, and the critical role MUSE plays because of the multi-scale nature of the physical processes involved. In this first paper, we focus on how comparisons between MUSE observations and theoretical models will significantly further our understanding of coronal heating mechanisms.

Current state-of-the-art spectrographs cannot resolve the fundamental spatial (sub-arcseconds) and temporal scales (less than a few tens of seconds) of the coronal dynamics of solar flares and eruptive phenomena. The highest resolution coronal data to date are based on imaging, which is blind to many of the processes that drive coronal energetics and dynamics. As shown by IRIS for the low solar atmosphere, we need high-resolution spectroscopic measurements with simultaneous imaging to understand the dominant processes. In this paper: (1) we introduce the Multi-slit Solar Explorer (MUSE), a spaceborne observatory to fill this observational gap by providing high-cadence (<20 s), sub-arcsecond resolution spectroscopic rasters over an active region size of the solar transition region and corona; (2) using advanced numerical models, we demonstrate the unique diagnostic capabilities of MUSE for exploring solar coronal dynamics, and for constraining and discriminating models of solar flares and eruptions; (3) we discuss the key contributions MUSE would make in addressing the science objectives of the Next Generation Solar Physics Mission (NGSPM), and how MUSE, the high-throughput EUV Solar Telescope (EUVST) and the Daniel K Inouye Solar Telescope (and other ground-based observatories) can operate as a distributed implementation of the NGSPM. This is a companion paper to De Pontieu et al. (2021a), which focuses on investigating coronal heating with MUSE.

Post-starburst galaxies are crucial to disentangling the effect of star formation and quenching on galaxy demographics. They comprise, however, a heterogeneous population of objects, described in numerous ways. To obtain a well-defined and uncontaminated sample, we take advantage of spatially resolved spectroscopy to construct an unambiguous sample of E + A galaxies - post-starburst systems with no observed ongoing star formation. Using data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) Survey, in the fourth generation of the Sloan Digital Sky Survey (SDSS-IV), we have identified 30 E + A galaxies that lie within the green valley of color-stellar mass space. We first identified E + A candidates by their central, single-fiber spectra and (u-r) color from SDSS DR15, and then further required each galaxy to exhibit E + A properties throughout the entirety of the system to three effective radii. We describe our selection criteria in detail, note common pitfalls in E + A identification, and introduce the basic characteristics of the sample. We will use this E + A sample, which has been assembled with stringent criteria and thus re-establishes a well-defined subpopulation within the broader category of post-starburst galaxies, to study the evolution of galaxies and their stellar populations in the time just after star formation within them is fully quenched.

Tracing the transition between the diffuse atomic interstellar medium (ISM) and cold, dense gas is crucial for deciphering the star formation cycle in galaxies. Here we present MACH, a new survey of cold neutral hydrogen (HI) absorption at $21\rm\,cm$ by the Karl G. Jansky Very Large Array. We target 42 bright background sources with $60<l<110^{\circ}$, $30<b<62^{\circ}$, significantly expanding the sample of publicly-available, sensitive $21\rm\,cm$ absorption outside the Galactic plane. With matching $21\rm\,cm$ emission data from the EBHIS survey, we measure the total column density and cold HI fraction, and quantify the properties of individual HI structures along each sightline via autonomous Gaussian decomposition. Combining the MACH sample with results from recent HI absorption surveys, we produce a robust characterization of the cool atomic medium at high and intermediate Galactic latitudes. We find that MACH HI has significantly smaller column density relative to samples at similar latitudes, and the detected cold HI structures have smaller line widths, temperatures and turbulent Mach numbers, suggesting that MACH probes a particularly quiescent region. Using all available observations, we compute the cumulative covering fraction ($c$) of cold HI at local velocities outside the disk: structures with $\tau>0.001$ are ubiquitous ($c\sim100\%$), whereas high optical depths ($\tau>1$) are extremely rare ($c\sim0\%$).

We investigate the properties of the double sequences of the Milky Way discs visible in the [$\alpha$/Fe] vs [Fe/H] diagram. In the framework of Galactic formation and evolution, we discuss the complex relationships between age, metallicity, [$\alpha$/Fe], and the velocity components. We study stars with measured chemical, seismic and astrometric properties from the APOGEE survey, the Kepler and Gaia satellites, respectively. We separate the [$\alpha$/Fe]-[Fe/H] diagram into 3 stellar populations: the thin disc, high-$\alpha$ metal-poor thick disc and high-$\alpha$ metal-rich thick disc and characterise each of these in the age-chemo-kinematics parameter space. We compare results obtained from different APOGEE data releases and using two recent age determinations. We use the Besan\c{c}on Galaxy model (BGM) to highlight selection biases and mechanisms not included in the model. The thin disc exhibits a flat age-metallicity relation while [$\alpha$/Fe] increases with stellar age. We confirm no correlation between radial and vertical velocities with [Fe/H], [$\alpha$/Fe] and age for each stellar population. Considering both samples, V$_\varphi$ decreases with age for the thin disc, while it is in the opposite direction for the h$\alpha$mp thick disc. Although the age distribution of the h$\alpha$mr thick disc is very close to that of the h$\alpha$mp thick disc between 7 and 14 Gyr, its kinematics seems to follow the thin disc kinematics. This feature, not predicted by the hypotheses included in the BGM, suggests a different origin and history of this population. Finally, we show a maximum of the dispersion of the vertical velocity, $\sigma_Z$, with age for the h$\alpha$mp thick disc around 8 Gyr. The comparisons with the BGM simulations suggest a more complex chemo-dynamical scheme to explain this feature, most likely including mergers and radial migration effects

The VISTA Variables in the Via Lactea (VVV) survey and its extension, have been monitoring about 560 square degrees of sky centred on the Galactic bulge and inner disc for nearly a decade. The photometric catalogue contains of order 10$^9$ sources monitored in the K$_s$ band down to 18 mag over hundreds of epochs from 2010-2019. Using these data we develop a decision tree classifier to identify microlensing events. As inputs to the tree, we extract a few physically motivated features as well as simple statistics ensuring a good fit to a microlensing model both on and off the event amplification. This produces a fast and efficient classifier trained on a set of simulated microlensing events and catacylsmic variables, together with flat baseline light curves randomly chosen from the VVV data. The classifier achieves 97 per cent accuracy in identifying simulated microlensing events in a validation set. We run the classifier over the VVV data set and then visually inspect the results, which produces a catalogue of 1,959 microlensing events. For these events, we provide the Einstein radius crossing time via a Bayesian analysis. The spatial dependence on recovery efficiency of our classifier is well characterised, and this allows us to compute spatially resolved completeness maps as a function of Einstein crossing time over the VVV footprint. We compare our approach to previous microlensing searches of the VVV. We highlight the importance of Bayesian fitting to determine the microlensing parameters for events with surveys like VVV with sparse data.

Significant galaxy mergers throughout cosmic time play a fundamental role in theories of galaxy evolution. The widespread usage of human classifiers to visually assess whether galaxies are in merging systems remains a fundamental component of many morphology studies. Studies that employ human classifiers usually construct a control sample, and rely on the assumption that the bias introduced by using humans will be evenly applied to all samples. In this work, we test this assumption and develop methods to correct for it. Using the standard binomial statistical methods employed in many morphology studies, we find that the merger fraction, error, and the significance of the difference between two samples are dependent on the intrinsic merger fraction of any given sample. We propose a method of quantifying merger biases of individual human classifiers and incorporate these biases into a full probabilistic model to determine the merger fraction and the probability of an individual galaxy being in a merger. Using 14 simulated human responses and accuracies, we are able to correctly label a galaxy as ''merger'' or ''isolated'' to within 1\% of the truth. Using 14 real human responses on a set of realistic mock galaxy simulation snapshots our model is able to recover the pre-coalesced merger fraction to within 10\%. Our method can not only increase the accuracy of studies probing the merger state of galaxies at cosmic noon, but also can be used to construct more accurate training sets in machine learning studies that use human classified data-sets.

Axions constituting dark matter (DM) are often considered to form a non-relativistic oscillating field. We explore bursts of relativistic axions from transient astrophysical sources, such as axion star explosions, where the sources are initially non-relativistic. For the QCD axion, bursts from collapsing axion stars lead to potentially detectable signals over a wide range of axion masses $10^{-15} \, \textrm{eV} \lesssim m \lesssim 10^{-7} \, \textrm{eV}$ in future experiments, such as ABRACADABRA, DMRadio and SHAFT. Unlike conventional cold axion DM searches, the sensitivity to axion bursts is not necessarily suppressed as $1/f$ for large decay constants $f$. The detection of axion bursts could provide new insights into the fundamental axion potential, which is challenging to probe otherwise. An ensemble of bursts in the distant past, in direct analogy with neutrinos, would give rise to a diffuse axion background distinct from the usual cold axion DM. Coincidence with other signatures, such as electromagnetic and gravitational-wave emission, would provide a new beyond-the-standard-model window into multi-messenger astronomy.

We present a novel mechanism which leads to the baryon asymmetry generation during the strong first order phase transition. If the bubble wall propagates with ultra-relativistic velocities, it has been shown that it can produce states much heavier than the scale of the transition and that those states are then out-of-equilibrium. In this paper, we show that the production mechanism can also induce CP-violation at one-loop level. We calculate those CP violating effects during the heavy particle production and show, that combined with baryon number violating interactions, those can lead to successful baryogenesis. Two models based on this mechanism are constructed and their phenomenology is discussed. Stochastic gravitational wave signals turn out to be generic signatures of this type of models.

It is shown that the slopes of the superhorizon hypermagnetic spectra produced by the variation of the gauge couplings are practically unaffected by the relative strength of the parity-breaking terms. A new method is proposed for the estimate of the gauge power spectra in the presence of pseudoscalar interactions during inflation. To corroborate the general results, various concrete examples are explicitly analyzed. Since the large-scale gauge spectra also determine the late-time magnetic fields it turns out that the pseudoscalar contributions have little impact on the magnetogenesis requirement. Conversely the parity-breaking terms crucially affect the gyrotropic spectra that may seed, in certain models, the baryon asymmetry of the universe. From the viewpoint of the effective field theory description of magnetogenesis scenarios these considerations hold generically for the whole class of inflationary models where the inflaton is not constrained by any underlying symmetry.

Quintessential inflation utilises a single scalar field to account for the observations of both cosmic inflation and dark energy. The requirements for modelling quintessential inflation are described and two explicit successful models are presented in the context of $\alpha$-attractors and Palatini modified gravity.

A dark matter search project needs and extremely low background radiation detector since the expected event rate of dark matter is less than a few events in one year in one tonne of the detector mass. The authors developed a highly radiopure NaI(Tl) crystal to search for dark matter. The best combination of the purification methods was developed, resulting $^{\mathrm{nat}}$K and $^{210}$Pb were less than 20 ppb and 5.7 $\mu$Bq/kg, respectively. The authors will construct a large volume detector system with high-purity NaI(Tl) crystals. The design and the performance of the prototype detector module will be reported in this article.

Theories of gravity that obey the Weak Equivalence Principle have the same Parametrised Post-Newtonian parameter $\gamma$ for all particles at all energies. The large Shapiro time delays of extragalactic sources allow us to put tight constraints on differences in $\gamma$ between photons of different frequencies from spectral lag data, since a non-zero $\Delta \gamma$ would result in a frequency-dependent arrival time. The majority of previous constraints have assumed that the Shapiro time delay is dominated by a few local massive objects, although this is a poor approximation for distant sources. In this work we consider the cosmological context of these sources by developing a source-by-source, Monte Carlo-based forward model for the Shapiro time delays by combining constrained realisations of the local density field using the BORG algorithm with unconstrained large-scale modes. Propagating uncertainties in the density field reconstruction and marginalising over an empirical model describing other contributions to the time delay, we use spectral lag data of Gamma Ray Bursts from the BATSE satellite to constrain $\Delta \gamma < 3.4 \times 10^{-15}$ at $1 \sigma$ confidence between photon energies of $25 {\rm \, keV}$ and $325 {\rm \, keV}$.

We present new equilibrium solutions of stationary models of magnetized thick disks (or tori) around Kerr black holes with synchronised scalar hair. The models reported here largely extend our previous results based on constant radial distributions of the specific angular momentum along the equatorial plane. We introduce a new way to prescribe the distribution of the disk's angular momentum based on a combination of two previous proposals and compute the angular momentum distribution outside the equatorial plane by resorting to the construction of von Zeipel cylinders. We find that the effect of the scalar hair on the black hole spacetime can yield significant differences in the disk morphology and properties compared to what is found if the spacetime is purely Kerr. Some of the tori built within the most extreme, background hairy black hole spacetime of our sample exhibit the appearance of two maxima in the gravitational energy density which impacts the radial profile distributions of the disk's thermodynamical quantities. The models reported in this paper can be used as initial data for numerical evolutions with GRMHD codes to study their stability properties. Moreover, they can be employed as illuminating sources to build shadows of Kerr black holes with scalar hair which might help further constrain the no-hair hypothesis as new observational data is collected.

During its two-year mission at comet 67P, Rosetta nearly continuously monitored the inner coma plasma environment for gas production rates varying over three orders of magnitude, at distances to the nucleus from a few to a few hundred km. To achieve the best possible measurements, cross-calibration of the plasma instruments is needed. We construct with two different physical models to cross-calibrate the electron density as measured by the Mutual Impedance Probe (MIP) to the ion current and spacecraft potential as measured by the Rosetta Langmuir Probe (LAP), the latter validated with the Ion Composition Analyser (ICA). We retrieve a continuous plasma density dataset for the entire cometary mission with a much improved dynamical range compared to any plasma instrument alone and, at times, improve the temporal resolution from 0.24-0.74~Hz to 57.8~Hz. The new density dataset is consistent with the existing MIP density dataset and covers long time periods where densities were too low to be measured by MIP. The physical model also yields, at 3~hour time resolution, ion flow speeds as well as a proxy for the solar EUV flux from the photoemission from the Langmuir Probes. We report on two independent mission-wide estimates of the ion flow speed which are consistent with the bulk H$_2$O$^+$ ion velocities as measured by ICA. We find the ion flow to consistently be much faster than the neutral gas over the entire mission, lending further evidence that the ions are collisionally decoupled from the neutrals in the coma. RPC measurements of ion speeds are therefore not consistent with the assumptions made in previously published plasma density models of the comet ionosphere at the start and end of the mission. Also, the measured EUV flux is perfectly consistent with independently derived values previously published from Johansson et al. (2017) and lends support for the conclusions drawn therein.

A large number of particle detectors employ liquid argon as their target material owing to its high scintillation yield and its ability to drift ionization charge for large distances. Scintillation light from argon is peaked at 128 nm and a wavelength shifter is required for its efficient detection. In this work we directly compare the light yield achieved in two identical liquid argon chambers, one of which is equipped with PolyEthylene Naphthalate (PEN) and the other with TetraPhenyl Butadiene (TPB) wavelength shifter. Both chambers are lined with enhanced specular reflectors and instrumented with SiPMs with a coverage fraction of approximately 1%, which represents a geometry comparable to the future large scale detectors. We measured the light yield of the PEN chamber to be 39.4$\pm$0.4(stat)$\pm$1.9(syst)% of the yield of the TPB chamber. Using a Monte Carlo simulation this result is used to extract the wavelength shifting efficiency of PEN relative to TPB equal to 47.2$\pm$5.7%. This result paves the way for the use of easily available PEN foils as a wavelength shifter, which can substantially simplify the construction of the future liquid argon detectors.

The thermodynamics-gravity conjecture reveals that one can derive the gravitational field equations by using the first law of thermodynamics and vice versa. Considering the entropy associated with the horizon in the form of non-extensive Tsallis entropy, $S\sim A^{\beta}$ here we first derive the corresponding gravitational field equations by applying the Clausius relation $\delta Q=T \delta S$ to the horizon. We then construct the Friedmann equations of Friedmann-Lema\^itre-Robertson-Walker (FLRW) universe based on Tsallis modified gravity (TMG). Moreover, in order to constrain the cosmological parameters of TMG model, we use observational data, including Planck cosmic microwave background (CMB), weak lensing, supernovae, baryon acoustic oscillations (BAO), and redshift-space distortions (RSD) data. Numerical results indicate that TMG model is more compatible with the low redshift observational data by predicting a lower value for the structure growth parameter $\sigma_8$ with respect to $\Lambda$CDM model. This implies that TMG model would slightly alleviate the $\sigma_8$ tension.