Papers for Monday, Jan 18 2021

We aim to demonstrate the scientific potential of the Gaia Early Data Release 3 (EDR3) for the study of the Milky Way structure and evolution. We used astrometric positions, proper motions, parallaxes, and photometry from EDR3 to select different populations and components and to calculate the distances and velocities in the direction of the anticentre. We explore the disturbances of the current disc, the spatial and kinematical distributions of early accreted versus in-situ stars, the structures in the outer parts of the disc, and the orbits of open clusters Berkeley 29 and Saurer 1. We find that: i) the dynamics of the Galactic disc are very complex with vertical asymmetries, and new correlations, including a bimodality with disc stars with large angular momentum moving vertically upwards from below the plane, and disc stars with slightly lower angular momentum moving preferentially downwards; ii) we resolve the kinematic substructure (diagonal ridges) in the outer parts of the disc for the first time; iii) the red sequence that has been associated with the proto-Galactic disc that was present at the time of the merger with Gaia-Enceladus-Sausage is currently radially concentrated up to around 14 kpc, while the blue sequence that has been associated with debris of the satellite extends beyond that; iv) there are density structures in the outer disc, both above and below the plane, most probably related to Monoceros, the Anticentre Stream, and TriAnd, for which the Gaia data allow an exhaustive selection of candidate member stars and dynamical study; and v) the open clusters Berkeley~29 and Saurer~1, despite being located at large distances from the Galactic centre, are on nearly circular disc-like orbits. We demonstrate how, once again, the Gaia are crucial for our understanding of the different pieces of our Galaxy and their connection to its global structure and history.

In the context of the INFN project PANDORA_Gr3 (Plasma for Astrophysics, Nuclear Decays Observation and Radiation for Archaeometry) and of multi-messenger astronomy, we propose a feasibility study for in-laboratory plasma's opacity investigation, in an environment resembling thermodynamic conditions typical of the ejecta of compact binary mergers containing at least a neutron star. We aim to advance knowledge on the physics of kilonovae, the electromagnetic transients following a merger, which are relevant for the study of the origin of heavy nuclei in the Universe produced via r-process nucleosynthesis. In this paper, we present preliminary results of numerical simulations for some physics cases considered in the light of a possible experimental setup for future in-laboratory opacity spectroscopic measurements.

We present the serendipitous detection of the two main OH maser lines at 1667 and 1665 MHz associated with IRAS 10597+5926 at z = 0.19612 in the untargeted Apertif Wide-area Extragalactic Survey (AWES), and the subsequent measurement of the OH 1612 MHz satellite line in the same source. With a total OH luminosity of log(L/L_Sun) = 3.90 +/- 0.03, IRAS 10597+5926 is the fourth brightest OH megamaser (OHM) known. We measure a lower limit for the 1667/1612 ratio of R_1612 > 45.9 which is the highest limiting ratio measured for the 1612 MHz OH satellite line to date. OH satellite line measurements provide a potentially valuable constraint by which to compare detailed models of OH maser pumping mechanisms. Optical imaging shows the galaxy is likely a late-stage merger. Based on published infrared and far ultraviolet fluxes, we find that the galaxy is an ultra luminous infrared galaxy (ULIRG) with log(L_TIR/L_Sun) = 12.24, undergoing a star burst with an estimated star formation rate of 179 +/- 40 M_Sun/yr. These host galaxy properties are consistent with the physical conditions responsible for very bright OHM emission. Finally, we provide an update on the predicted number of OH masers that may be found in AWES, and estimate the total number of OH masers that will be detected in each of the individual main and satellite OH 18 cm lines.

We study how half-mass radii, central mass densities ($\Sigma_1$), and color gradients change as galaxies evolve. We classify $\sim7,000$ galaxies into sixteen groups with similar spectral shapes; each group represents a different evolutionary stage. We find that different galaxy types populate different regions of both size-mass and $\Sigma_1$-mass space. The nine star-forming groups lie along the integrated star-forming $\Sigma_1$-mass relation. However, these star-forming groups form steep parallel relations in the size-mass plane, with slopes similar to the quiescent size-mass relation. These steep slopes can be explained as a transformation of the star-forming $\Sigma_1$-mass relation and its scatter. We identify three types of transitional galaxies. Green valley and post-starburst galaxies are similarly compact at $z>1.5$; however, their distinct color gradients indicate that the two populations represent different pathways to quenching. Post-starburst galaxies have flat color gradients and compact structures, consistent with a fast quenching pathway which requires structural change and operates primarily at high redshift. Green valley galaxies have negative color gradients, and are both larger and more numerous towards lower redshift. These galaxies are consistent with slow quenching without significant structural change. We find that dusty star-forming galaxies at $z\gtrsim2$ are very compact, and may represent the "burst" before post-starburst galaxies; at $z\lesssim2$, dusty star-forming galaxies are extended and have shallow color gradients consistent with slow quenching. Our results suggest that star-forming galaxies grow gradually up the \sig-mass relation until (a) they naturally reach the high $\Sigma_1$ values required for quiescence, or (b) a compaction-type event rapidly increases their $\Sigma_1$.

We predict the stellar mass -- halo mass (SMHM) relationship for dwarf galaxies, using simulated galaxies with peak halo masses of M$_{\rm peak} = 10^{11}$ M$_{\odot}$ down into the ultra-faint dwarf range to M$_{\rm peak} =$ 10$^7$ M$_{\odot}$. Our simulated dwarfs have stellar masses of M$_{\rm star} = $ 790 M$_{\odot}$ to $8.2 \times 10^8$ M$_{\odot}$, with corresponding $V$-band magnitudes from $-2$ to $-18.5$. For M$_{\rm peak} > 10^{10}$ M$_{\odot}$, the simulated SMHM relationship agrees with literature determinations, including exhibiting a small scatter of 0.3 dex. However, the scatter in the SMHM relation increases for lower-mass halos. We first present results for well-resolved halos that contain a simulated stellar population, but recognize that whether a halo hosts a galaxy is inherently mass resolution dependent. We thus adopt a probabilistic model to populate "dark" halos below our resolution limit to predict an "intrinsic" slope and scatter for the SMHM relation. We fit linearly growing log-normal scatter in stellar mass, which grows to more than 1 dex at M$_{\rm peak}$ $=$ 10$^8$ M$_{\odot}$. At the faintest end of the SMHM relation probed by our simulations, a galaxy cannot be assigned a unique halo mass based solely on its luminosity. Instead, we provide a formula to stochastically populate low-mass halos following our results. Finally, we show that our growing log-normal scatter steepens the faint-end slope of the predicted stellar mass function.

We present numerical simulations of misaligned discs around a spinning black hole covering a range of parameters. Previous simulations have shown that discs that are strongly warped by a forced precession -- in this case the Lense-Thirring effect from the spinning black hole -- can break apart into discrete discs or rings that can behave quasi-independently for short timescales. With the simulations we present here, we confirm that thin and highly inclined discs are more susceptible to disc tearing than thicker or low inclination discs, and we show that lower values of the disc viscosity parameter lead to instability at lower warp amplitudes. This is consistent with detailed stability analysis of the warped disc equations. We find that the growth rates of the instability seen in the numerical simulations are similar across a broad range of parameters, and are of the same order as the predicted growth rates. However, we did not find the expected trend of growth rates with viscosity parameter. This may indicate that the growth rates are affected by numerical resolution, or that the wavelength of the fastest growing mode is a function of local disc parameters. Finally, we also find that disc tearing can occur for discs with a viscosity parameter that is higher than predicted by a local stability analysis of the warped disc equations. In this case, the instability manifests differently producing large changes in the disc tilt locally in the disc, rather than the large changes in disc twist that typically occur in lower viscosity discs.

Accretion discs around black holes power some of the most luminous objects in the Universe. Discs that are misaligned to the black hole spin can become warped over time by Lense-Thirring precession. Recent work has shown that strongly warped discs can become unstable, causing the disc to break into discrete rings producing a more dynamic and variable accretion flow. In a companion paper, we present numerical simulations of this instability and the resulting dynamics. In this paper, we discuss the implications of this dynamics for accreting black hole systems, with particular focus on the variability of Active Galactic Nuclei (AGN). We discuss the timescales on which variability might manifest, and the impact of the observer orientation with respect to the black hole spin axis. When the disc warp is unstable near the inner edge of the disc, we find quasi periodic behaviour of the inner disc which may explain the recent quasi periodic eruptions observed in, for example, the Seyfert 2 galaxy GSN 069 and in the galactic nucleus of RX J1301.9+2747. These eruptions are thought to be similar to the `heartbeat' modes observed in some X-ray binaries (e.g. GRS 1915+105 and IGR J17091-3624). When the instability manifests at larger radii in the disc, we find that the central accretion rate can vary on timescales that may be commensurate with, e.g., changing-look AGN. We therefore suggest that some of the variability properties of accreting black hole systems may be explained by the disc being significantly warped, leading to disc tearing.

Recent observations of the periodic Fast Radio Burst source 180916.J0158+65 (FRB 180916) find small linear polarization position angle swings during and between bursts, with a burst activity window that becomes both narrower and earlier at higher frequencies. Although the observed chromatic activity window disfavors models of periodicity in FRB 180916 driven by the occultation of a neutron star by the optically-thick wind from a stellar companion, the connection to theories where periodicity arises from the motion of a bursting magnetar remains unclear. In this paper, we show how altitude-dependent radio emission from magnetar curvature radiation, with bursts emitted from regions which are asymmetric with respect to the magnetic dipole axis, can lead to burst activity windows and polarization consistent with the recent observations. In particular, the fact that bursts arrive systematically earlier at higher frequencies disfavors theories where the FRB periodicity arises from forced precession of a magnetar by a companion or fallback disk, but is consistent with theories where periodicity originates from a slowly-rotating or freely-precessing magnetar. Several observational tests are proposed to verify/differentiate between the remaining theories, and pin-down which theory explains the periodicity in FRB 180916.

We present results of two-year photometric monitoring of Comet C/2011 J2 (LINEAR) that spans the time period from February 2013 through December 2014, before and after perihelion passage. The observations were conducted with broadband R filter. Analysis of observations of Comet C/2011 J2 (LINEAR) allows estimating the nucleus radius as function of geometric albedo and phase-darkening coefficient. Furthermore, our observations showed split of the comet. Fragment (B) additional to the primary C/2011 J2 nucleus was unambiguously detected; relative velocity is estimated to be near 0.17 ''/day. We derive the Afrho parameter and estimate the dust production rate in Comet C/2011 J2 (LINEAR) over the entire run of observations. We found a noticeable increase in Afrho parameter between September 18, 2014 and November 5, 2014, epoch when the nucleus presumably got split. On September 28, 2014, we observed a transit of the 16-magnitude star (USNO-A2 1275-18299027) near nucleus of the Comet C/2011 J2 (LINEAR). We retrieve the optical depth of the coma 0.034. The filling factor f that corresponds to such optical depth is in good quantitative agreement with the value that can be derived from the Afrho parameter under reasonable assumption on geometric albedo of cometary dust.

In this era of Gaia and ALMA, dynamical stellar mass measurements provide benchmarks that are independent of observations of stellar characteristics and their uncertainties. These benchmarks can then be used to validate and improve stellar evolutionary models, which can lead to both imprecise and inaccurate mass predictions for pre-main-sequence, low-mass stars. We present the dynamical stellar masses derived from disks around three M-stars (FP Tau, J0432+1827, and J1100-7619) using ALMA observations of $^{12}$CO (J=2--1) and $^{13}$CO (J=2--1) emission. These are the first dynamical stellar mass measurements for J0432+1827 and J1100-7619 and the most precise measurement for FP Tau. Fiducial stellar evolutionary model tracks, which do not include any treatment of magnetic activity, agree with the dynamical measurement of J0432+1827 but underpredict the mass by $\sim$60\% for FP Tau and $\sim$80\% for J1100-7619. Possible explanations for the underpredictions include inaccurate assumptions of stellar effective temperature, undetected binarity for J1100-7619, and that fiducial stellar evolutionary models are not complex enough to represent these stars. In the former case, the stellar effective temperatures would need to be increased by $\sim$40K to $\sim$340K to reconcile the fiducial model predictions with the dynamically-measured masses. In the latter case, we show that the dynamical masses can be reproduced using results from stellar evolutionary models with starspots, which incorporate fractional starspot coverage to represent the manifestation of magnetic activity. Folding in low-mass M-stars from the literature and assuming that the stellar effective temperatures are imprecise but accurate, we find tentative evidence of a relationship between fractional starspot coverage and observed effective temperature for these young, cool stars.

We report the discovery of two unsaturated organic species, trans-(E)-cyanovinylacetylene and vinylcyanoacetylene, using the second data release of the GOTHAM deep survey towards TMC-1 with the 100 m Green Bank Telescope. For both detections, we performed velocity stacking and matched filter analyses using Markov chain Monte Carlo simulations, and for trans-(E)-cyanovinylacetylene, three rotational lines were observed at low signal-to-noise (${\sim}$3$\sigma$). From this analysis, we derive column densities of $2\times10^{11}$ and $3\times10^{11}$ cm$^{-2}$ for vinylcyanoacetylene and trans-(E)-cyanovinylacetylene, respectively, and an upper limit of $<2\times10^{11}$ cm$^{-2}$ for trans-(Z)-cyanovinylacetylene. Comparisons with G3//B3LYP semi-empirical thermochemical calculations indicate abundances of the [H$_3$C$_5$N}] isomers are not consistent with their thermodynamic stability, and instead their abundances are mainly driven by dynamics. We provide discussion into how these species may be formed in TMC-1, with reference to related species like vinyl cyanide (CH$_2$=CHCN). As part of this discussion, we performed the same analysis for ethyl cyanide (CH$_3$CH$_2$CN), the hydrogenation product of CH$_2$=CHCN. This analysis provides evidence -- at 4.17$\sigma$ significance -- an upper limit to the column density of $<4\times10^{11}$ cm$^{-2}$; an order of magnitude lower than previous upper limits towards this source.

In this work, we propose a cosmological model-independent and non-local method to constrain the Hubble Constant $H_0$. Inspired by the quasi cosmological model-independent and $H_0$-free properties of the `shifted' Hubble diagram of HII galaxies (HIIGx) defined by Wei et al. (2016), we joint analyze it with the parametric type Ia supernova (SN Ia) Hubble diagram (e.g. the joint-lightcurves-analysis sample, JLA) and get a Bayesian Inference of Hubble constant, $H_0 = 71 \pm 20\ \mathrm{km\ s^{-1}\ Mpc^{-1}}$. Although with large uncertainty, we find that $H_0$ is only strongly degenerate with the B-band absolute magnitude ($M_B$) of SN Ia but almost independent on other nuisance parameters. Therefore the accuracy can be simultaneously improved by a tight constraint of $M_B$ through a cosmological and $H_0$ independent way. This method can be extended further to get more-literally non-local results of $H_0$ by using other Hubble diagrams at higher redshifts.

The Australian Square Kilometre Array Pathfinder (ASKAP) collects images of the sky at radio wavelengths with an unprecedented field of view, combined with a high angular resolution and sub-millijansky sensitivities. The large quantity of data produced is used by the ASKAP Variables and Slow Transients (VAST) survey science project to study the dynamic radio sky. Efficient pipelines are vital in such research, where searches often form a `needle in a haystack' type of problem to solve. However, the existing pipelines developed among the radio-transient community are not suitable for the scale of ASKAP datasets. In this paper we provide a technical overview of the new "VAST Pipeline": a modern and scalable Python-based data pipeline for transient searches, using up-to-date dependencies and methods. The pipeline allows source association to be performed at scale using the Pandas DataFrame interface and the well-known Astropy crossmatch functions. The Dask Python framework is used to parallelise operations as well as scale them both vertically and horizontally, by means of a cluster of workers. A modern web interface for data exploration and querying has also been developed using the latest Django web framework combined with Bootstrap.

A compact green astro-comb with 43-GHz repetition rate is developed based on a Ti:Sapphire optical frequency comb (OFC) and a mode-selecting cavity. The OFC's large repetition rate of 1.6 GHz eases the requirements for the mode-selecting cavity. Unnecessary frequency-modes of the OFC are suppressed down to $5 \times 10^{-4}$ at 535 nm - 550 nm using a single mode-selecting cavity with 70-MHz linewidth. The radial velocity precision $\sigma \sim 1.4$ m/s is achieved at the High Dispersion Echelle Spectrosraph for the Okayama 188-cm telescope of the National Astronomical Observatory of Japan using our astro-comb. With further improvements of the mode-selecting cavity and removal of fiber modal noises, our system will provide a simple, compact, and precise astro-comb setup in visible wavelength region.

We present the "StoP" photometer-polarimeter (Stokes-polarimeter) used for observations with the 1-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences since the beginning of 2020. We describe the instrument and its parameters in observations performed in the photometric and polarimetric modes. We demonstrate the capabilities of the instrument through the polarimetry of the blazar S5 0716+714 and compare the results with those earlier obtained with the 6-m telescope.

Some observations such as those presented in Walker et al. show that the observed entropy profiles of the intra-cluster medium (ICM) deviate from the power-law prediction of adiabatic simulations. This implies that non-gravitational processes, which are absent in the simulations, may be important in the evolution of the ICM, and by quantifying the deviation, we may be able to estimate the feedback energy in the ICM and use it as a probe of the non-gravitational processes. To address this issue we calculate the ICM entropy profiles in a sample of 47 galaxy clusters and groups, which have been observed out to at least $\sim r_{500}$ with Chandra, XMM-Newton and/or Suzaku, by constructing a physical model to incorporate the effects of both gravity and non-gravitational processes to fit the observed gas temperature and surface brightness profiles via Bayesian statistics. After carefully evaluating the effects of systematic errors, we find that the gas entropy profiles derived with best-fit results of our model are consistent with the simulation-predicted power-law profile near the virial radius, while the flattened profiles reported previously can be explained by introducing the gas clumping effect, the existence of which is confirmed in 19 luminous targets in our sample. We calculate the total feedback energy per particle and find that it decreases from $\sim 10$ keV at the center to about zero at $\sim 0.35$$r_{200}$ and is consistent with zero outside $\sim 0.35$$r_{200}$, implying the upper limit of the feedback efficiency $\sim 0.02$ for the super-massive black holes hosted in the brightest cluster galaxies.

We have examined 13 pointed observations of the TeV emitting high synchrotron peak blazar PKS 2155-304, taken by the Suzaku satellite throughout its operational period. We found that the blazar showed large-amplitude intraday variabilities in the soft (0.8 - 1.5 keV) and the hard (1.5 - 8.0 keV) bands in the light curves. Spectral variability on intraday timescales is estimated using the hardness ratio. The blazar usually becomes harder when brighter and vice versa, following the typical behavior of high synchrotron peak blazars. The power spectral density (PSD) analyses of 11 out of 13 light curves in the total energy (0.8 - 8.0 keV) are found to be red-noise dominated, with power-law spectral indices that span a large range, from -2.81 to -0.88. Discrete correlation function analyses of all the 13 light curves between the soft and the hard bands show that they are well correlated and peak at, or very close to, zero lag. This indicates that the emissions in soft and hard bands are probably cospatial and emitted from the same population of leptons. Considering fluxes versus variability timescales, we found no correlation on intraday timescales, implying that X-ray emission from PKS 2155-304 is not dominated by simple changes in the Doppler factor. We briefly discuss the most likely emission mechanisms responsible for the observed flux and spectral variabilities and place constraints on magnetic field strength and Lorentz factors of the electrons emitting the X-rays in the most likely scenario.

The timing-based localization, which utilize the triangulation principle with the different arrival time of gamma-ray photons, with a fleet of Cubesats is a unique and powerful solution for the future all-sky gamma-ray observation, which is a key for identification of the electromagnetic counterpart of the gravitational wave sources. The Cubesats Applied for MEasuring and Localising Transients (CAMELOT) mission is now being promoted by the Hungarian and Japanese collaboration with a basic concept of the nine Cubesats constellations in low earth orbit. The simulation framework for estimation of the localization capability has been developed including orbital parameters, an algorithm to estimate the expected observed profile of gamma-ray photons, finding the peak of the cross-correlation function, and a statistical method to find a best-fit position and its uncertainty. It is revealed that a degree-scale localization uncertainty can be achieved by the CAMELOT mission concept for bright short gamma-ray bursts, which could be covered by future large field of view ground-based telescopes. The new approach utilizing machine-learning approach is also investigated to make the procedure automated for the future large scale constellations. The trained neural network with 10$^6$ simulated light curves generated by the artificial short burst templates successfully predicts the time-delay of the real light curve and achieves a comparable performance to the cross-correlation algorithm with full automated procedures.

We present a bright galaxy sample with accurate and precise photometric redshifts (photo-zs), selected using $ugriZYJHK_\mathrm{s}$ photometry from the Kilo-Degree Survey (KiDS) Data Release 4 (DR4). The highly pure and complete dataset is flux-limited at $r<20$ mag, covers $\sim1000$ deg$^2$, and contains about 1 million galaxies after artifact masking. We exploit the overlap with Galaxy And Mass Assembly (GAMA) spectroscopy as calibration to determine photo-zs with the supervised machine learning neural network algorithm implemented in the ANNz2 software. The photo-zs have mean error of $|\langle \delta z \rangle| \sim 5 \times 10^{-4}$ and low scatter (scaled mean absolute deviation of $\sim 0.018(1+z)$), both practically independent of the $r$-band magnitude and photo-z at $0.05 < z_\mathrm{phot} < 0.5$. Combined with the 9-band photometry, these allow us to estimate robust absolute magnitudes and stellar masses for the full sample. As a demonstration of the usefulness of these data we split the dataset into red and blue galaxies, use them as lenses and measure the weak gravitational lensing signal around them for five stellar mass bins. We fit a halo model to these high-precision measurements to constrain the stellar-mass--halo-mass relations for blue and red galaxies. We find that for high stellar mass ($M_\star>5\times 10^{11} M_\odot$), the red galaxies occupy dark matter halos that are much more massive than those occupied by blue galaxies with the same stellar mass. The data presented here will be publicly released via the KiDS webpage.

We explore constraints imposed by shear-driven instabilities on the acceleration of energetic particles in relativistic shearing flows. We show that shearing layers in large-scale AGN jets are likely to encompass a sizeable fraction ($\geq 0.1$) of the jet radius, requiring seed injection of GeV electrons for efficient acceleration. While the diffusion process may depend on pre-developed turbulence if injection occurs at higher energies, electron acceleration to PeV and proton acceleration to EeV energies appears possible within the constraints imposed by jet stability.

The operation of the future Cherenkov Telescope Array (CTA), the next generation of Imaging Atmospheric Cherenkov Telescopes (IACTs), as an Open Observatory requires data products and analysis tools to be accessible and usable by the scientific community. This requirement has impelled gamma-ray astronomers to develop open-source science tools and standardised data formats. The objective of this presentation is to provide a perspective on the current effort to define a prototypical format for high-level gamma-ray astronomical data. The initiative, started by the IACT community in 2016, has gained in the years 2018 and 2019 full recognition thanks to the first public data releases in such format by current-generation IACTs, followed by a series of papers employing it either to test the capabilities of open-source gamma-ray analysis tools or to showcase examples of reproducible, multi-instrument analysis.

We present the results from an Australian Square Kilometre Array Pathfinder search for radio variables on timescales of hours. We conducted an untargeted search over a 30 deg$^2$ field, with multiple 10-hour observations separated by days to months, at a central frequency of 945 MHz. We discovered six rapid scintillators from 15-minute model-subtracted images with sensitivity of $\sim 200\,\mu$Jy/beam; two of them are extreme intra-hour variables with modulation indices up to $\sim 40\%$ and timescales as short as tens of minutes. Five of the variables are in a linear arrangement on the sky with angular width $\sim 1$ arcmin and length $\sim 2$ degrees, revealing the existence of a huge plasma filament in front of them. We derived kinematic models of this plasma from the annual modulation of the scintillation rate of our sources, and we estimated its likely physical properties: a distance of $\sim 4$ pc and length of $\sim 0.1$ pc. The characteristics we observe for the scattering screen are incompatible with published suggestions for the origin of intra-hour variability leading us to propose a new picture in which the underlying phenomenon is a cold tidal stream. This is the first time that multiple scintillators have been detected behind the same plasma screen, giving direct insight into the geometry of the scattering medium responsible for enhanced scintillation.

Strong accretion shocks are expected to illuminate the warm-hot inter-galactic medium encompassed by the filaments of the cosmic web, through synchrotron radio emission. Given their high sensitivity, low-frequency large radio facilities may already be able to detect signatures of this extended radio emission from the region in between two close and massive galaxy clusters. In this work we exploit the non-detection of such diffuse emission by deep observations of two pairs of relatively close ($\simeq 10$ Mpc) and massive ($M_{500}\geq 10^{14}M_\odot$) galaxy clusters using the LOw-Frequency ARray (LOFAR). By combining the results from the two putative inter-cluster filaments, we derive new independent constraints on the median strength of inter-galactic magnetic fields: $B_{\rm Mpc}< 2.5\times 10^2\,\rm nG\,(95\%\, \rm CL)$. Based on cosmological simulations and assuming a primordial origin of the B-fields, these estimates can be used to limit the amplitude of primordial seed magnetic fields: $B_0\leq10\,\rm nG$. We advise the observation of similar cluster pairs as a powerful tool to set tight constraints on the amplitude of extragalactic magnetic fields.

The use of the full potential of stellar seismology is made difficult by the improper modeling of the upper-most layers of solar-like stars and their influence on the modeled frequencies. Our knowledge on these \emph{surface effects} has improved thanks to the use of 3D hydrodynamical simulations but the calculation of eigenfrequencies relies on empirical models for the description of the Lagrangian perturbation of turbulent pressure: the reduced-$\Gamma_1$ model (RGM) and the gas-$\Gamma_1$ model (GGM). Starting from the fully compressible turbulence equations, we derive both the GGM and RGM models using a closure to model the flux of turbulent kinetic energy. It is found that both models originate from two terms: the source of turbulent pressure due to compression produced by the oscillations and the divergence of the flux of turbulent pressure. It is also demonstrated that they are both compatible with the adiabatic approximation but also imply a number of questionable assumptions mainly regarding mode physics. Among others hypothesis, one has to neglect the Lagrangian perturbation of the dissipation of turbulent kinetic energy into heat and the Lagrangian perturbation of buoyancy work.

We explore the possibility of inferring the properties of the Galactic neutron star population through machine learning. In particular, in this paper we focus on their dynamical characteristics and show that an artificial neural network is able to estimate with high accuracy the parameters which control the current positions of a mock population of pulsars. For this purpose, we implement a simplified population-synthesis framework (where selection biases are neglected at this stage) and concentrate on the natal kick-velocity distribution and the distribution of birth distances from the Galactic plane. By varying these and evolving the pulsar trajectories in time, we generate a series of simulations that are used to train and validate a suitably structured convolutional neural network. We demonstrate that our network is able to recover the parameters governing the kick-velocity and Galactic height distribution with a mean relative error of about $10^{-2}$. We discuss the limitations of our idealized approach and study a toy problem to introduce selection effects in a phenomenological way by incorporating the observed proper motions of 216 isolated pulsars. Our analysis highlights that increasing the sample of pulsars with accurate proper motion measurements by a factor of $\sim$10, one of the future breakthroughs of the Square Kilometer Array, we might succeed in constraining the birth spatial and kick-velocity distribution of the neutron stars in the Milky Way with high precision through machine learning.

The Auger Surface Detector consists of a large array of water Cherenkov detector tanks each with a volume of 12,000 liters, for the detection of high energy cosmic rays. The accuracy in the measurement of the integrated signal amplitude of the detector unit has been studied using experimental air shower data. It can be described as a Poisson-like term with a normalization constant that depends on the zenith angle of the primary cosmic ray. This dependence reflects the increasing contribution to the signal of the muonic component of the shower, both due to the increasing muon/electromagnetic (e+- and gamma) ratio and muon track length with zenith angle.

2-aminooxazole (2AO), a N-heterocyclic molecule, has been proposed as an intermediate in prebiotic syntheses. It has been demonstrated that it can be synthesized from small molecules such as cyanamide and glycoaldehyde, which are present in interstellar space. The aim of this work is to provide infrared spectra, in the solid phase for conditions typical of astrophysical environments and to estimate its stability toward UV photons and cosmic rays. Infrared (4000-600 cm$^{-1}$) absorption spectra at 20 K, 180 K, and 300 K, IR band strengths, and room temperature UV (120-250 nm) absorption spectra are given for the first time for this species. Destruction cross-sections of 9.5 10$^{-18}$ cm$^2$ and 2 10$^{-16}$ cm$^2$ were found in the irradiation at 20 K of pure 2AO and 2AO:H$_2$O ices with UV (6.3-10.9 eV) photons or 5 keV electrons, respectively. These data were used for the estimate of half-life times for the molecule in different environments. It is estimated that 2AO could survive UV radiation and cosmic rays in the ice mantles of dense clouds beyond cloud collapse. In contrast, it would be very unstable at the surface of cold Solar System bodies like Kuiper belt objects, but the molecule could still survive within dust grain agglomerates or cometesimals.

The presence of A/F-type {\it Kepler} hybrid stars extending across the entire $\delta$ Sct-$\gamma$ Dor instability strips and beyond remains largely unexplained. In order to better understand these particular stars, we performed a multi-epoch spectroscopic study of 49 candidate A/F-type hybrid stars and one cool(er) hybrid object detected by the {\it Kepler} mission. We determined a lower limit of 27 % for the multiplicity fraction. For six spectroscopic systems, we also reported long-term variations of the time delays. For four systems, the time delay variations are fully coherent with those of the radial velocities and can be attributed to orbital motion. We aim to improve the orbital solutions for those systems with long orbital periods (order of 4-6 years) among the {\it Kepler} hybrid stars. The orbits are computed based on a simultaneous modelling of the RVs obtained with high-resolution spectrographs and the photometric time delays derived from time-dependent frequency analyses of the {\it Kepler} light curves. We refined the orbital solutions of four spectroscopic systems with A/F-type {\it Kepler} hybrid component stars: KIC 4480321, 5219533, 8975515 and KIC 9775454. Simultaneous modelling of both data types analysed together enabled us to improve the orbital solutions, obtain more robust and accurate information on the mass ratio, and identify the component with the short-period $\delta$ Sct-type pulsations. In several cases, we were also able to derive new constraints for the minimum component masses. From a search for regular frequency patterns in the high-frequency regime of the Fourier transforms of each system, we found no evidence of tidal splitting among the triple systems with close (inner) companions. However, some systems exhibit frequency spacings which can be explained by the mechanism of rotational splitting.

Astrophysical observations show complex organic molecules (COMs) in the gas phase of protoplanetary disks. X-rays emitted from the central young stellar object (YSO) that irradiate interstellar ices in the disk, followed by the ejection of molecules in the gas phase, are a possible route to explain the abundances observed in the cold regions. This process, known as X-ray photodesorption, needs to be quantified for methanol-containing ices. This paper I focuses on the case of X-ray photodesorption from pure methanol ices. We aim at experimentally measuring X-ray photodesorption yields of methanol and its photo-products from pure CH$_3$OH ices, and to shed light on the mechanisms responsible for the desorption process. We irradiated methanol ices at 15 K with X-rays in the 525 - 570 eV range. The release of species in the gas phase was monitored by quadrupole mass spectrometry, and photodesorption yields were derived. Under our experimental conditions, the CH$_3$OH X-ray photodesorption yield from pure methanol ice is 10$^{-2}$ molecule/photon at 564 eV. Photo-products such as CH$_4$, H$_2$CO, H$_2$O, CO$_2$ , and CO also desorb at increasing efficiency. X-ray photodesorption of larger COMs, which can be attributed to either ethanol, dimethyl ether, and/or formic acid, is also detected. The physical mechanisms at play are discussed and must likely involve the thermalization of Auger electrons in the ice, thus indicating that its composition plays an important role. Finally, we provide desorption yields applicable to protoplanetary disk environments for astrochemical models. The X-rays are shown to be a potential candidate to explain gas-phase abundances of methanol in disks. However, more relevant desorption yields derived from experiments on mixed ices are mandatory to properly support the role played by X-rays in nonthermal desorption of methanol (see paper II).

We present the first analysis of results from the SuperWASP Variable Stars Zooniverse project, which is aiming to classify 1.6 million phase-folded light curves of candidate stellar variables observed by the SuperWASP all sky survey with periods detected in the SuperWASP periodicity catalogue. The resultant data set currently contains $>$1 million classifications corresponding to $>$500,000 object-period combinations, provided by citizen scientist volunteers. Volunteer-classified light curves have $\sim$89 per cent accuracy for detached and semi-detached eclipsing binaries, but only $\sim$9 per cent accuracy for rotationally modulated variables, based on known objects. We demonstrate that this Zooniverse project will be valuable for both population studies of individual variable types and the identification of stellar variables for follow up. We present preliminary findings on various unique and extreme variables in this analysis, including long period contact binaries and binaries near the short-period cutoff, and we identify 301 previously unknown binaries and pulsators. We are now in the process of developing a web portal to enable other researchers to access the outputs of the SuperWASP Variable Stars project.

In this paper, we employ Machine Learning algorithms on multi-mission observations for the classification of accretion states of outbursting Black hole X-ray binaries for the first time. Archival data from RXTE, Swift, MAXI and AstroSat observatories are used to generate the hardness intensity diagrams (HIDs) for outbursts of the sources XTE J1859+226 (1999 outburst), GX 339-4 (2002, 2004, 2007 and 2010 outbursts), IGR J17091-3624 (2016 outburst), and MAXI J1535-571 (2017 outburst). Based on variation of X-ray flux, hardness ratios, presence of various types of Quasi-periodic Oscillations (QPOs), photon indices and disk temperature, we apply clustering algorithms like K-Means clustering and Hierarchical clustering to classify the accretion states (clusters) of each outburst. As multiple parameters are involved in the classification process, we show that clustering algorithms club together the observations of similar characteristics more efficiently than the `standard' method of classification. We also infer that K-Means clustering provides more reliable results than Hierarchical clustering. We demonstrate the importance of the classification based on machine learning by comparing it with results from `standard' classification.

We investigate strongly gravitationally lensed type II supernovae (LSNe II) for time-delay cosmography incorporating microlensing effects, which expands on previous microlensing studies of type Ia supernovae (SNe Ia). We use the radiative-transfer code ${\rm \small TARDIS}$ to recreate five spectra of the prototypical SN 1999em at different times within the plateau phase of the light curve. The microlensing-induced deformations of the spectra and light curves are calculated by placing the SN into magnification maps generated with the code ${\rm \small GERLUMPH}$. We study the impact of microlensing on the color curves and find that there is no strong influence on them during the investigated time interval of the plateau phase. The color curves are only weakly affected by microlensing due to the almost achromatic behavior of the intensity profiles. However, the lack of non-linear structure in the color curves makes time-delay measurements difficult given the possible presence of differential dust extinction. Therefore, we further investigate SN phase inference through spectral absorption lines under the influence of microlensing and Gaussian noise. As the spectral features shift to longer wavelengths with progressing time after explosion, the measured wavelength of a specific absorption line provides information on the epoch of the SN. The comparison between retrieved epochs of two observed lensing images then gives the time delay of the images. We find that the phase retrieval method using spectral features yields accurate delays with uncertainties $\small {\lesssim}$2 days, making it a promising approach.

Publications in astrophysics are nowadays mainly published and read in digitized formats. Astrophysical publications in both research and in popular outreach often use colorful representations of stars to indicate various stellar types, that is, different spectral types or effective temperatures. Computer generated and computer displayed imagery has become an integral part of stellar astrophysics communication. There is, however, no astrophysically motivated standard color palette for illustrative representations of stars and some stars are actually represented in misleading colors. We use pre-computed PHOENIX and TLUSTY stellar model spectra and convolve them with the three standard color matching functions for human color perception between 360$\,$nm and 830$\,$nm. The color matching functions represent the three sets of receptors in the eye that respond to red, green, and blue light. For a grid of main sequence stars with effective temperatures between 2300$\,$K and 55,000$\,$K of different metallicities we present the red-blue-green and hexadecimal color codes that can be used for digitized color representations of stars as if seen from space. We find significant deviations between the color codes of stars computed from stellar spectra and from a black body radiator of the same effective temperature. We illustrate the main sequence in the color wheel and demonstrate that there are no yellow, green, cyan, or purple stars. Red dwarf stars (spectral types M0V - M9V) actually look orange to the human eye. Old white dwarfs such as WD$\,$1856$+$534, host to a newly discovered transiting giant planet candidate, occur pale orange to the human eye, not white. Our freely available software can be used to generate color codes for any input spectrum such as those from planets, galaxies, quasars etc.

We consider a short rollercoaster cosmology based on two stages of monodromy inflation separated by a stage of matter domination, generated after the early inflaton falls out of slow roll. If the first stage is controlled by a flat potential, $V \sim \phi^p$ with $p < 1$ and lasts ${\cal N} \sim 30 - 40$ efolds, the scalar and tensor perturbations at the largest scales will fit the CMB perfectly, and produce relic gravity waves with $0.02 \lesssim r \lesssim 0.06$, which can be tested by LiteBIRD and CMB-S4 experiments. If in addition the first inflaton is strongly coupled to a hidden sector $U(1)$, there will be an enhanced production of vector fluctuations near the end of the first stage of inflation. These modes convert rapidly to tensors during the short epoch of matter domination, and then get pushed to superhorizon scales by the second stage of inflation, lasting another $20-30$ efolds. This band of gravity waves is chiral, arrives today with wavelengths in the range of $10^8$ km, and with amplitudes greatly enhanced compared to the long wavelength CMB modes by vector sources. It is therefore accessible to LISA. Thus our model presents a rare early universe theory predicting several simultaneous signals testable by a broad range of gravity wave searches in the very near future.

In this paper, we study the dynamics of k-essence in loop quantum cosmology (LQC). The study indicates that the loop quantum gravity (LQG) effect plays a key role only in the early epoch of the universe and is diluted at the later stage. The fixed points in LQC are basically consistent with that in standard Friedmann-Robertson-Walker (FRW) cosmology. For most of the attractor solutions, the stability conditions in LQC are in agreement with that for the standard FRW universe. But for some special fixed point, more tighter constraints are imposed thanks to the LQG effect.

We investigate a relativistic adaptation of the Lattice Boltzmann Method that reproduces the equations of motion for a turbulent, two-dimensional, massless hydrodynamic system. The classical Lattice Boltzmann Method and its extension to relativistic fluid dynamics is described. The numeric formulation is evaluated using a zero-averaged stirring force introduced into the numerics to induce turbulence, and the flow characteristics produced are compared to properties of a classical turbulent hydrodynamic flow. The model can reasonably be expected to offer quantitative simulations of electron fluid flows in graphene or Kagome lattices.

Cosmological models of the early or late universe exhibit (quasi) de Sitter space-times with different stability properties. Considering models derived from string theory, the swampland program does not provide for now a definite characterisation of this stability. In this work we focus on de Sitter solutions of 10d type II supergravities, candidates for classical de Sitter string backgrounds: surprisingly, all known examples are unstable with $\eta_V < -1$. We aim at proving the existence of such a systematic tachyon, and getting formally a bound on the value of $\eta_V$. To that end, we develop three methods, giving us various sufficient conditions for having a tachyon upon assumptions, in analogy with de Sitter no-go theorems. We use this knowledge to find 10 new de Sitter solutions of type IIB supergravity, that have different tachyons, higher $\eta_V$ values and new 6d geometries.