Papers for Wednesday, Oct 06 2021

We determine radial velocities and mass flow rates in a sample of 54 local spiral galaxies by modelling high-resolution and high-sensitivity data of the atomic hydrogen emission line. We found that, although radial inflow motions seem to be slightly preferred over outflow motions, their magnitude is generally small. Most galaxies show radial flows of only a few km/s throughout their HI disks, either inwards or outwards, without any clear increase in magnitude in the outermost regions, as we would expect for continuous radial accretion. Gas mass flow rates for most galaxies are less than 1 M$_\odot$/yr. Over the entire sample, we estimated an average inflow rate of 0.3 M$_\odot$/yr outside the optical disk and of 0.1 M$_\odot$/yr in the outskirts of the HI disks. These inflow rates are about 5-10 times smaller than the average star formation rate of 1.4 M$_\odot$/yr. Our study suggests that there is no clear evidence for systematic radial accretion inflows that alone could feed and sustain the star formation process in the inner regions of local spiral galaxies at its current rate.

The polarization of the 21-cm radiation from the epoch of reionization arises from Thomson scattering of 21-cm photons from free electrons and provides information that complements that from the intensity fluctuation. Previous work showed that a direct detection of this signal will be difficult, and hinted that the signal might be enhanced via correlation with other tracers. Here, we discuss the cross-correlation between the cosmic microwave background (CMB) polarization and the 21-cm polarization. We treat reionization using an analytical model with parameters calibrated by semi-numerical simulations. We then derive the cross-correlation angular power spectrum using the total-angular-momentum formalism. We also provide a noise analysis to test against two closely related, but subtly different, null hypotheses. First, we assume no reionization as a null hypothesis, and determine how well this null hypothesis could be ruled out by an observed 21cm-CMB polarization correlation. Second, we determine how well the null hypothesis of no 21-cm polarization can be ruled out by seeking the cross-correlation, assuming reionization is established from the CMB. We find that the first question could be answered by a synergy of ambitious next-generation 21-cm and CMB missions, whereas the second question will still remain out of reach.

Astrometry -- the precise measurement of positions and motions of celestial objects -- has emerged as a promising avenue for characterizing the dark matter population in our Galaxy. By leveraging recent advances in simulation-based inference and neural network architectures, we introduce a novel method to search for global dark matter-induced gravitational lensing signatures in astrometric datasets. Our method based on neural likelihood-ratio estimation shows significantly enhanced sensitivity to a cold dark matter population and more favorable scaling with measurement noise compared to existing approaches based on two-point correlation statistics, establishing machine learning as a powerful tool for characterizing dark matter using astrometric data.

The Wide-Field Infrared Transient Explorer (WINTER) is a new 1 $\text{deg}^2$ seeing-limited time-domain survey instrument designed for dedicated near-infrared follow-up of kilonovae from binary neutron star (BNS) and neutron star-black hole mergers. WINTER will observe in the near-infrared Y, J, and short-H bands (0.9-1.7 microns, to $\text{J}_{AB}=21$ magnitudes) on a dedicated 1-meter telescope at Palomar Observatory. To date, most prompt kilonova follow-up has been in optical wavelengths; however, near-infrared emission fades slower and depends less on geometry and viewing angle than optical emission. We present an end-to-end simulation of a follow-up campaign during the fourth observing run (O4) of the LIGO, Virgo, and KAGRA interferometers, including simulating 625 BNS mergers, their detection in gravitational waves, low-latency and full parameter estimation skymaps, and a suite of kilonova lightcurves from two different model grids. We predict up to five new kilonovae independently discovered by WINTER during O4, given a realistic BNS merger rate. Using a larger grid of kilonova parameters, we find that kilonova emission is $\approx$2 times longer lived and red kilonovae are detected $\approx$1.5 times further in the infrared than the optical. For 90% localization areas smaller than 150 (450) $\rm{deg}^{2}$, WINTER will be sensitive to more than 10% of the kilonova model grid out to 350 (200) Mpc. We develop a generalized toolkit to create an optimal BNS follow-up strategy with any electromagnetic telescope and present WINTER's observing strategy with this framework. This toolkit, all simulated gravitational-wave events, and skymaps are made available for use by the community.

A transformative science case for the proposed next-generation Very Large Array (ngVLA) is resolving the surfaces of nearby stars, both spatially and temporally, enabled by the combination of milliarcsecond-scale resolution and unprecedented sensitivity to thermal radio emission. In a previous study, we demonstrated the feasibility of stellar imaging with simulated observations of nearby stars, using both traditional CLEAN techniques and newly developed regularized maximum likelihood (RML) imaging methods for image reconstruction. In this memo, we present a continued study of stellar imaging with the ngVLA, evaluating the imaging capability of the Revision D (henceforth Rev D) Main Array configuration compared to the previous Revision C (henceforth Rev C) configuration. We find that the Rev D configuration, with more uniform coverage and better circular symmetry, improves the synthesized beam, resulting in better CLEAN reconstructions of simulated images of evolved stars with complex morphology, especially with robust weighting. However, the highly non-Gaussian nature of the synthesized beam still persists with both robust and natural weightings in the Rev D configuration and continues to limit the image fidelity of image reconstructions with non-uniform weighting. The RML methods show stable performance that is resilient to different array configurations with image quality comparable to or better than CLEAN methods in the presented simulation, consistent with our previous work. Our simulation results suggest that the Rev D configuration will provide a better deconvolution beam compared with the Rev C configuration, which would enhance the imaging capability for non-uniform weighting, and they continue to demonstrate that RML methods are an attractive choice, even for the improved array configuration.

We use Hubble space telescope data of 45 nearby star-forming galaxies to investigate properties of Lyman-alpha (Ly$\alpha$) halos, Ly$\alpha$ morphology, and the star-forming characteristics of galaxies. We study how the morphology of Ly$\alpha$ emission is related to other Ly$\alpha$ observables. Furthermore, we study the interdependencies of Ly$\alpha$ morphological quantities. We studied the spatial extent of Ly$\alpha$ using surface brightness profiles in i) using circular apertures and ii) within faint Ly$\alpha$ isophotes. We measured the average intensity and the size of the regions with a high star formation rate density. The morphology of the galaxies was quantified by computing centroid position, axis ratio, and position angle in the Ly$\alpha$, ultraviolet continuum, and I band maps. We found that galaxies with more extended star-forming regions possess larger Ly$\alpha$ halos. Furthermore, galaxies with more elongated Ly$\alpha$ morphology are also more extended in Ly$\alpha$. Our data suggest that Ly$\alpha$ bright galaxies appear rounder in their Ly$\alpha$ morphology, and there is less of a contribution from their Ly$\alpha$ halo to their overall luminosity. We compared our results with studies at high redshift and found that whilst the Ly$\alpha$ extent in the inner regions of the galaxies in our sample are similar to the high$-z$ Ly$\alpha$ emitters (LAEs), Ly$\alpha$ halos are more extended in high$-z$ LAEs. Our analysis suggests that the Ly$\alpha$ morphology affects the measurement of other observable quantities concerning Ly$\alpha$ emission, and some of the conclusions drawn from high redshift LAEs might be biased towards galaxies with specific Ly$\alpha$ shapes. In particular, faint Ly$\alpha$ emitters have larger Ly$\alpha$ scale lengths and halo fractions. This implies that faint Ly$\alpha$ emitters are harder to detect at high redshift than previously believed.

The high-redshift intergalactic medium (IGM) and the primeval galaxy population are rapidly becoming the new frontier of extra-galactic astronomy. We investigate the IGM properties and their connection to galaxies at $z\geq5.5$ under different assumptions for the ionizing photon escape and the nature of dark matter, employing our novel THESAN radiation-hydrodynamical simulation suite, designed to provide a comprehensive picture of the emergence of galaxies in a full reionization context. Our simulations have realistic `late' reionization histories, match available constraints on global IGM properties and reproduce the recently-observed rapid evolution of the mean free path of ionizing photons. We additionally examine high-z Lyman-$\alpha$ transmission. The optical depth evolution is consistent with data, and its distribution suggests an even-later reionization than simulated, although with a strong sensitivity to the source model. We show that the effects of these two unknowns can be disentangled by characterising the spectral shape and separation of Lyman-$\alpha$ transmission regions, opening up the possibility to observationally constrain both. For the first time in simulations, THESAN reproduces the modulation of the Lyman-$\alpha$ flux as a function of galaxy distance, demonstrating the power of coupling a realistic galaxy formation model with proper radiation-hydrodynamics. We find this feature to be extremely sensitive on the timing of reionization, while being relatively insensitive to the source model. Overall, THESAN produces a realistic IGM and galaxy population, providing a robust framework for future analysis of the high-z Universe.

Diffuse radio emission at the centre of galaxy clusters has been observed both in merging clusters on scales of Mpc, called giant radio haloes, and in relaxed systems with a cool-core on smaller scales, named mini haloes. Giant radio haloes and mini haloes are thought to be distinct classes of sources. However, recent observations have revealed the presence of diffuse radio emission on Mpc scales in clusters that do not show strong dynamical activity. RX J1720.1+2638 is a cool-core cluster, presenting both a bright central mini halo and a fainter diffuse, steep-spectrum emission extending beyond the cluster core that resembles giant radio halo emission. In this paper, we present new observations performed with the LOFAR Low Band Antennas (LBA) at 54 MHz. These observations, combined with data at higher frequencies, allow us to constrain the spectral properties of the radio emission. The large-scale emission presents an ultra-steep spectrum with $\alpha_{54}^{144}\sim3.2$. The radio emission inside and outside the cluster core have strictly different properties, as there is a net change in spectral index and they follow different radio-X-ray surface brightness correlations. We argue that the large-scale diffuse emission is generated by particles re-acceleration after a minor merger. While for the central mini halo we suggest that it could be generated by secondary electrons and positrons from hadronic interactions of relativistic nuclei with the dense cool-core gas, as an alternative to re-acceleration models.

The photoevaporation of protoplanetary discs by nearby massive stars present in their birth cluster plays a vital role in their evolution. Previous modelling assumes that the disc behaves like a classical Keplerian accretion disc out to a radius where the photoevaporative outflow is launched. There is then an abrupt change in the angular velocity profile, and the outflow is modelled by forcing the fluid parcels to conserve their specific angular momenta. Instead, we model externally photoevaporating discs using the slim disc formalism. The slim disc approach self consistently includes the advection of radial and angular momentum as well as angular momentum redistribution by internal viscous torques. Our resulting models produce a smooth transition from a rotationally supported Keplerian disc to a photoevaporative driven outflow, where this transition typically occurs over ~4-5 scale heights. The penetration of UV photons predominately sets the radius of the transition and the viscosity's strength plays a minor role. By studying the entrainment of dust particles in the outflow, we find a rapid change in the dust size and surface density distribution in the transition region due to the steep gas density gradients present. This rapid change in the dust properties leaves a potentially observable signature in the continuum spectral index of the disc at mm wavelengths. Using the slim disc formalism in future evolutionary calculations will reveal how both the gas and dust evolves in their outer regions and the observable imprints of the external photoevaporation process.

Tidal Disruption Events (TDEs) occur when stars pass close to supermassive black holes, and have long been predicted to emit cosmic rays and neutrinos. Recently the TDE AT2109dsg was identified in spatial and temporal coincidence with high-energy neutrino IC191001A as part of the Zwicky Transient Facility (ZTF) neutrino follow-up program, providing the first direct observational evidence supporting these objects as multi-messenger sources. In these proceedings, I will place the recent results of our ZTF neutrino follow-up program into the broader context of developments in TDE and neutrino astronomy.

It is well established from cosmological simulations that dark matter haloes are not precisely self-similar and an additional parameter, beyond their concentration, is required to accurately describe their spherically-averaged mass density profiles. We present, for the first time, a model to consistently predict both halo concentration, $c$, and this additional `shape' parameter, $\alpha$, for a halo of given mass and redshift for a specified cosmology. Following recent studies, we recast the dependency on mass, redshift, and cosmology to a dependence on `peak height'. We show that, when adopting the standard definition of peak height, which employs the so-called spherical top hat (STH) window function, the concentration--peak height relation has a strong residual dependence on cosmology (i.e., it is not uniquely determined by peak height), whereas the $\alpha$--peak height relation is approximately universal when employing the STH window function. Given the freedom in the choice of window function, we explore a simple modification of the STH function, constraining its form so that it produces universal relations for concentration and $\alpha$ as a function of peak height using a large suite of cosmological simulations. It is found that universal relations for the two density profile parameters can indeed be derived and that these parameters are set by the linear power spectrum, $P(k)$, filtered on different scales. We show that the results of this work generalise to any (reasonable) combination of $P(k)$ and background expansion history, $H(z)$, resulting in accurate predictions of the density profiles of dark matter haloes for a wide range of cosmologies.

The high incidence rate of the O VI $\lambda\lambda$1032,1038 absorption around low-redshift, $\sim$$L^*$ star-forming galaxies has generated interest in studies of the circumgalactic medium. We use the high-resolution EAGLE cosmological simulation to analyze the circumgalactic O VI gas around $z\approx0.3$ star-forming galaxies. Motivated by the limitation that observations do not reveal where the gas lies along the line-of-sight, we compare the O VI measurements produced by gas within fixed distances around galaxies and by gas selected using line-of-sight velocity cuts commonly adopted by observers. We show that gas selected by a velocity cut of $\pm300$ km s$^{-1}$ or $\pm500$ km s$^{-1}$ produces a higher O VI column density, a flatter column density profile, and a higher covering fraction compared to gas within one, two, or three times the virial radius ($r_\mathrm{vir}$) of galaxies. The discrepancy increases with impact parameter and worsens for lower mass galaxies. For example, compared to the gas within 2$r_\mathrm{vir}$, identifying the gas using velocity cuts of 200-500 km s$^{-1}$ increases the O VI column density by 0.2 dex (0.1 dex) at 1$r_\mathrm{vir}$ to over 0.75 dex (0.7 dex) at $\approx2$$r_\mathrm{vir}$ for galaxies with stellar masses of $10^{9}$-$10^{9.5}$ $\rm M_\odot$ ($10^{10}$-$10^{10.5}$ $\rm M_\odot$). We furthermore estimate that excluding O VI outside $r_\mathrm{vir}$ decreases the circumgalactic oxygen mass measured by Tumlinson et al. (2011) by over 50%. Our results demonstrate that gas at large line-of-sight separations but selected by conventional velocity windows has significant effects on the O VI measurements and may not be observationally distinguishable from gas near the galaxies.

Gravitational wave detectors are starting to reveal the redshift evolution of the binary black hole (BBH) merger rate, $R_{\mathrm{BBH}}(z)$. We make predictions for $R_{\mathrm{BBH}}(z)$ as a function of black hole mass for systems originating from isolated binaries. To this end, we investigate correlations between the delay time and black hole mass by means of the suite of binary population synthesis simulations, COMPAS. We distinguish two channels: the common envelope (CE), and the stable Roche-lobe overflow (RLOF) channel, characterised by whether the system has experienced a common envelope or not. We find that the CE channel preferentially produces BHs with masses below about $30\rm{M}_{\odot}$ and short delay times ($t_{\rm delay} \lesssim 1$Gyr), while the stable RLOF channel primarily forms systems with BH masses above $30\rm{M}_{\odot}$ and long delay times ($t_{\rm delay} \gtrsim 1$Gyr). We provide a new fit for the metallicity specific star-formation rate density based on the Illustris TNG simulations, and use this to convert the delay time distributions into a prediction of $R_{\mathrm{BBH}}(z)$. This leads to a distinct redshift evolution of $R_{\mathrm{BBH}}(z)$ for high and low primary BH masses. We furthermore find that, at high redshift, $R_{\mathrm{BBH}}(z)$ is dominated by the CE channel, while at low redshift it contains a large contribution ($\sim 40\%$) from the stable RLOF channel. Our results predict that, for increasing redshifts, BBHs with component masses above $30\rm{M}_{\odot}$ will become increasingly scarce relative to less massive BBH systems. Evidence of this distinct evolution of $R_{\mathrm{BBH}}(z)$ for different BH masses can be tested with future detectors.

Observations of the dynamics of solar coronal structures are necessary to investigate space weather phenomena and global heating of the corona. The profiles of high temperature lines emitted by the hot plasma are usually integrated by narrow band filters or recorded by classical spectroscopy. We present in this paper details of a new transportable instrument (under construction) for imaging spectroscopy: the Solar Line Emission Dopplerometer (SLED). It uses the Multi-channel Subtractive Double Pass (MSDP) technique, which combines the advantages of both filters and narrow slit spectrographs, i.e. high temporal, spatial and spectral resolutions. The SLED will measure at high cadence (1 Hz) the line-of-sight velocities (Doppler shifts) of hot coronal loops, in the forbidden lines of FeX 637. nm and FeXIV 530.3 nm. It will follow the dynamics of fast evolving events of solar activity such as flares or Coronal Mass Ejections (CMEs), and also study coronal heating by short period waves. Observations will be performed with the coronagraph at the Lomnicky Stit Observatory (LSO, in Slovakia) or during total eclipses. The SLED will also observe the dynamics of solar prominences in Halpha 656.3 nm or HeD3 587.6 nm lines when mounted on the Bialkow coronagraph (near Wroclaw, Poland). It is fully compatible with polarimetric measurements by various techniques.

We present a spectroscopic and photometric calibration to derive effective temperatures $T_{\mathrm{eff}}$ and metallicities [Fe/H] for M dwarfs, based on a Principal Component Analysis of 147 spectral indices measured off moderate resolution $R \sim 11\,000$), high S/N ($>100$) spectra in the $\lambda \lambda$ 8390-8834 region, plus the J$-$H color. Internal uncertainties, estimated by the residuals, are 81 K and 0.12 dex, respectively, for $T_{\mathrm{eff}}$ and [Fe/H], the calibrations being valid for 3050 K $< T_{\mathrm{eff}} <$ 4100 K and $-$0.45 $<$ [Fe/H] $<$ $+$0.50 dex. The PCA calibration is a competitive model-independent method to derive $T_{\mathrm{eff}}$ and [Fe/H] for large samples of M dwarfs, well suited to the available database of far-red spectra. The median uncertainties are 105 K and 0.23 dex for $T_{\mathrm{eff}}$ and [Fe/H], respectively, estimated by Monte Carlo simulations. We compare our values to other works based on photometric and spectroscopic techniques and find median differences 75 $\pm$ 273 K and 0.02 $\pm$ 0.31 dex for $T_{\mathrm{eff}}$ and [Fe/H], respectively, achieving good accuracy but relatively low precision. We find considerable disagreement in the literature between atmospheric parameters for stars in common. We use the new calibration to derive $T_{\mathrm{eff}}$ and [Fe/H] for 178 K7-M5 dwarfs, many previously unstudied. Our metallicity distribution function for nearby M dwarfs peaks at [Fe/H]$\sim$-0.10 dex, in good agreement with the RAVE distribution for GK dwarfs. We present radial velocities (internal precision 1.4 km/s) for 99 objects without previous measurements. The kinematics of the sample shows it to be fully dominated by thin/thick disk stars, excepting the well-known high-velocity Kapteyn's star.

Large-scale beamforming with radio interferometers has the potential to revolutionize the science done with pulsars and fast radio bursts by improving the survey efficiency for these sources. We describe a wide-field beamformer for the MeerKAT radio telescope and outline strategies to optimally design such surveys. A software implementation of these techniques, ${\rm M{\small OSAIC}}$ is introduced and its application in the MeerKAT telescope is presented. We show initial results using the beamformer by observing a globular cluster to track several pulsars simultaneously and demonstrate the source localization capability of this observation.

Methanol, one of the simplest complex organic molecules in the Interstellar Medium (ISM), has been shown to be present and extended in cold environments such as starless cores. We aim at studying methanol emission across several starless cores and investigate the physical conditions at which methanol starts to be efficiently formed, as well as how the physical structure of the cores and their surrounding environment affect its distribution. Methanol and C$^{18}$O emission lines at 3 mm have been observed with the IRAM 30m telescope within the large program "Gas phase Elemental abundances in Molecular CloudS" (GEMS) towards 66 positions across 12 starless cores in the Taurus Molecular Cloud. A non-LTE radiative transfer code was used to compute the column densities in all positions. We then used state-of-the-art chemical models to reproduce our observations. We have computed N(CH$_3$OH)/N(C$^{18}$O) column density ratios for all the observed offsets, and two different behaviours can be recognised: the cores where the ratio peaks at the dust peak, and the cores where the ratio peaks with a slight offset with respect to the dust peak ($\sim $10000 AU). We suggest that the cause of this behaviour is the irradiation on the cores due to protostars nearby which accelerate energetic particles along their outflows. The chemical models, which do not take into account irradiation variations, can reproduce fairly well the overall observed column density of methanol, but cannot reproduce the two different radial profiles observed. We confirm the substantial effect of the environment onto the distribution of methanol in starless cores. We suggest that the clumpy medium generated by protostellar outflows might cause a more efficient penetration of the interstellar radiation field in the molecular cloud and have an impact on the distribution of methanol in starless cores.

Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available.

The IceCube neutrino telescope discovered PeV-energy neutrinos originating beyond our Galaxy with an energy flux that is comparable to that of GeV-energy gamma rays and EeV-energy cosmic rays. These neutrinos provide the only unobstructed view of the cosmic accelerators that power the highest energy radiation reaching us from the universe. We will review the results from IceCube's first decade of operations, emphasizing the measurement of the diffuse multiflavored neutrino flux from the universe and the identification of the supermassive black hole TXS 0506+056 as a source of cosmic neutrinos and, therefore, cosmic rays. We will speculate on the lessons learned for multimessenger astronomy, among them that extragalactic neutrino sources may be a relatively small subset of the cosmic accelerators observed in high-energy gamma rays and that these may be gamma-ray-obscured at the times that they emit neutrinos.

We demonstrate that using up to seven stellar abundance ratios can place observational constraints on the star formation histories (SFHs) of Local Group dSphs, using the Sculptor dSph as a test case. We use a one-zone chemical evolution model to fit the overall abundance patterns of $\alpha$ elements (which probe the core-collapse supernovae that occur shortly after star formation), $s$-process elements (which probe AGB nucleosynthesis at intermediate delay times), and iron-peak elements (which probe delayed Type Ia supernovae). Our best-fit model indicates that Sculptor dSph has an ancient SFH, consistent with previous estimates from deep photometry. However, we derive a total star formation duration of $\sim0.9$ Gyr, which is shorter than photometrically-derived SFHs. We explore the effect of various model assumptions on our measurement and find that modifications to these assumptions still produce relatively short SFHs of duration $\lesssim1.4$ Gyr. Our model is also able to compare sets of predicted nucleosynthetic yields for supernovae and AGB stars, and can provide insight into the nucleosynthesis of individual elements in Sculptor dSph. We find that observed [Mn/Fe] and [Ni/Fe] trends are most consistent with sub-$M_{\mathrm{Ch}}$ Type Ia supernova models, and that a combination of "prompt" (delay times similar to core-collapse supernovae) and "delayed" (minimum delay times $\gtrsim50$ Myr) $r$-process events may be required to reproduce observed [Ba/Mg] and [Eu/Mg] trends.

An Adaptive secondary mirror (ASM) allows for the integration of adaptive optics (AO) into the telescope itself. Adaptive secondary mirrors, based on hybrid variable reluctance (HVR) actuator technology, developed by TNO, provide a promising path to telescope-integrated AO. HVR actuators have the advantage of allowing mirrors that are stiffer, more power efficient, and potentially less complex than similar, voice-coil based ASM's. We are exploring the application of this technology via a laboratory testbed that will validate the technical approach. In parallel, we are developing conceptual designs for ASMs at several telescopes including the Automated Planet Finder Telescope (APF) and for Keck Observatory. An ASM for APF has the potential to double the light through the slit for radial velocity measurements, and dramatically improved the image stability. An ASM for WMKO enables ground layer AO correction and lower background infrared AO observations, and provides for more flexible deployment of instruments via the ability to adjust the location of the Cassegrain focus.

Advancements in high-efficiency hybrid variable reluctance (HVR) actuators are an enabling technology for building the next generation of large-format deformable mirrors, including adaptive secondary mirrors. The Netherlands Organization for Applied Scientific Research (TNO) has developed a new style of hybrid variable reluctance actuator that requires approximately seventy-five times less power to operate as compared to the traditional style of voice-coil actuators. We present the initial performance results from laboratory testing of TNO's latest 19-actuator prototype deformable mirror, FLASH. We report the actuator cross-coupling, linearity, hysteresis, natural shape flattening, and drift as measured with a Zygo interferometer and a set of four capacitive sensors. We also present results of the dynamic performance of the FLASH on sub-millisecond timescales to estimate the limits of this technology for high-contrast imaging adaptive optics. We confirm that this technology has strong potential for use in on-sky adaptive secondary mirrors without the need for active cooling.

In light of the most recent cosmological observations, we provide new updated constraints on the slow-roll inflation in different extended scenarios beyond the $\Lambda\rm{CDM}$ cosmological model. Along with the usual six parameters, we simultaneously vary different combinations of additional parameters, including the running of the scalar spectral index $\alpha_s$, its running of running $\beta_s$, the tensor amplitude $r$ and the spatial curvature $\Omega_k$. From the Planck 2018 data, we find no evidence for a scalar running or a running of running, while analyzing the Atacama Cosmology Telescope data combined with WMAP 9-years observations data we find a preference for non-zero $\alpha_s$ and $\beta_s$ at the level of 2.9$\sigma$ and 2.7$\sigma$, respectively. Anyway, this preference is reduced when the tensor amplitude can vary in the model or $\beta_s$ is fixed to zero. The upper bound on $r$ is only slightly affected by the additional parameters while the differences in the datasets can remarkably change the compatibility among the different inflationary models, sometimes leading to discordant conclusions.

The ultraviolet (UV) to sub-millimetre (submm) spectral energy distribution of galaxies can be roughly divided into two sections: the stellar emission (attenuated by dust) at UV to near-infrared wavelengths and dust emission at longer wavelengths. In Dobbels et al. (2020), we show that these two sections are strongly related, and we can predict the global dust properties from the integrated UV to mid-infrared emission with the help of machine learning techniques. We investigate if these machine learning techniques can also be extended to resolved scales. Our aim is to predict resolved maps of the specific dust luminosity, specific dust mass, and dust temperature starting from a set of surface brightness images from UV to mid-infrared wavelengths. We used a selection of nearby galaxies retrieved from the DustPedia sample, in addition to M31 and M33. These were convolved and resampled to a range of pixel sizes, ranging from 150 pc to 3 kpc. We trained a random forest model which considers each pixel individually. We find that the predictions work well on resolved scales, with the dust mass and temperature having a similar root mean square error as on global scales (0.32 dex and 3.15 K on 18" scales respectively), and the dust luminosity being noticeably better (0.11 dex). We find no significant dependence on the pixel scale. Predictions on individual galaxies can be biased, and we find that about two-thirds of the scatter can be attributed to scatter between galaxies (rather than within galaxies). A machine learning approach can be used to create dust maps, with its resolution being only limited to the input bands, thus achieving a higher resolution than Herschel. These dust maps can be used to improve global estimates of dust properties, they can lead to a better estimate of dust attenuation, and they can be used as a constraint on cosmological simulations that trace dust.

Dust and gas in protoplanetary disks dissipate as central stars evolve. In order to estimate the dust dissipation timescales in the protoplanetary disks, we stacked the WISE 12, 22, and the AKARI 90 $\mu$m survey images of known T Tauri stars and derived the average fluxes, well below the survey flux limit in the 90 $\mu$m band. We classified 4,783 T Tauri stars into three age groups, which are young ($<$2 Myr), mid-age (2-6 Myr), and old ($>$6 Myr) groups, and stacked the WISE 12 and 22 and the AKARI 90 $\mu$m images in each group. The photometry of the stacked image shows the flux decay timescales of 1.4$\pm$0.2, 1.38$\pm$0.05, and 1.4$^{+0.6}_{-0.5}$ Myr in the 12, 22, and 90 $\mu$m bands, respectively. In optically thin disks with one-solar luminosity central stars, the 12 and 22 $\mu$m fluxes are attributed to the emission from the intermediate ($\sim$ 1 au) region and the 90 $\mu$m flux corresponds to that from the outer ($\sim$ 10 au) region in the disk. We hence conclude the dust dissipation timescale is $\tau_{\rm med,dust}\sim$1.4 Myr in the intermediate disks and is $\tau_{\rm outer, dust}=$1.4$^{+0.6}_{-0.5}$ Myr in the outer disks. The dust-dissipation time difference between the outer and intermediate disks is $\Delta\tau_{\rm dust} = \tau_{\rm outer,dust} -\tau_{\rm med,dust}=0.0^{+0.6}_{-0.5}$ Myr, indicating that the dust in the intermediate and outer disks dissipates on almost the same timescale.

We present the results from an analysis of deep Herschel far-infrared observations of the edge-on disk galaxy NGC 3079. The PSF-cleaned PACS images at 100 and 160 um display a 25 kpc x 25 kpc X-shape structure centered on the nucleus that is similar in extent and orientation to that seen in Halpha, X-rays, and the far-ultraviolet. One of the dusty filaments making up this structure is detected in the SPIRE 250 um map out to ~25 kpc from the nucleus. The match between the far-infrared filaments and those detected at other wavelengths suggests that the dusty material has been lifted out of the disk by the same large-scale galactic wind that has produced the other structures in this object. A closer look at the central 10 kpc x 10 kpc region provides additional support for this scenario. The dust temperatures traced by the 100-to-160 um flux ratios in this region are enhanced within a biconical region centered on the active galactic nucleus, aligned along the minor axis of the galaxy, and coincident with the well-known double-lobed cm-wave radio structure and Halpha-X-ray nuclear superbubbles. PACS imaging spectroscopy of the inner 6-kpc region reveals broad [C II] 158 um emission line profiles and OH 79 um absorption features along the minor axis of the galaxy with widths well in excess of those expected from beam smearing of the disk rotational motion. This provides compelling evidence that the cool material traced by the [C II] and OH features directly interacts with the nuclear ionized and relativistic outflows traced by the Halpha, X-ray, and radio emission.

Amines, in particular primary amines (R-NH$_2$) are closely related to the primordial synthesis of amino acids since they share the same structural backbone. However, only limited number of amines has been identified in the ISM which prevents us from studying their chemistry as well as their relation to pre-biotic species that could lead to the emergence of life. In this letter, we report the first interstellar detection of vinylamine (C$_2$H$_3$NH$_2$) and tentative detection of ethylamine (C$_2$H$_5$NH$_2$) towards the Galactic Centre cloud G+0.693-0.027. The derived abundance with respect to H$_2$ is (3.3$\pm$0.4)$\times$10$^{-10}$ and (1.9$\pm$0.6)$\times$10$^{-10}$, respectively. The inferred abundance ratios of C$_2$H$_3$NH$_2$ and C$_2$H$_5$NH$_2$ with respect to methylamine (CH$_3$NH$_2$) are $\sim$0.02 and $\sim$0.008 respectively. The derived abundance of C$_2$H$_3$NH$_2$, C$_2$H$_5$NH$_2$ and several other NH$_2$-bearing species are compared to those obtained towards high-mass and low-mass star-forming regions. Based on recent chemical and laboratory studies, possible chemical routes for the interstellar synthesis of C$_2$H$_3$NH$_2$ and C$_2$H$_5$NH$_2$ are discussed.

We present Quokka, a new subcycling-in-time, block-structured adaptive mesh refinement (AMR) radiation hydrodynamics code optimised for graphics processing units (GPUs). Quokka solves the equations of hydrodynamics with the piecewise parabolic method (PPM) in a method-of-lines formulation, and handles radiative transfer via the variable Eddington tensor (VET) radiation moment equations with a local closure. In order to maximise GPU performance, we combine explicit-in-time evolution of the radiation moment equations with the reduced speed-of-light approximation. We show results for a wide range of test problems for hydrodynamics, radiation, and coupled radiation hydrodynamics. On uniform grids in 3D, we achieve a peak of 93 million hydrodynamic updates per second per GPU, and 22 million radiation hydrodynamic updates per second per GPU. For radiation hydrodynamics problems on uniform grids in 3D, our code also scales from 4 GPUs to 256 GPUs with an efficiency of 80 percent. The code is publicly released under an open-source license on GitHub.

We present orbits for 24 binaries in the field of open cluster NGC 2516 (~150 Myr) and 13 binaries in the field of open cluster NGC 2422 (~130 Myr) using results from a multi-year radial velocity survey of the cluster cores. Six of these systems are double-lined spectroscopic binaries (SB2s). We fit these RV variable systems with orvara, a MCMC-based fitting program that models Keplerian orbits. We use precise stellar parallaxes and proper motions from Gaia EDR3 to determine cluster membership. We impose a barycentric radial velocity prior on all cluster members; this significantly improves our orbital constraints. Two of our systems have periods between 5 and 15 days, the critical window in which tides efficiently damp orbital eccentricity. These binaries should be included in future analyses of circularization across similarly-aged clusters. We also find a relatively flat distribution of binary mass ratios, consistent with previous work. With the inclusion of TESS lightcurves for all available targets, we identity target 378-036252 as a new eclipsing binary. We also identify a field star whose secondary has a mass in the brown dwarf range, as well as two cluster members whose RVs suggest the presence of an additional companion. Our orbital fits will help constrain the binary fraction and binary properties across stellar age and across stellar environment.

Interplanetary Coronal Mass Ejection (ICME) shocks are known to accelerate particles and contribute significantly to Solar Energetic Particle (SEP) events. We have performed Magnetohydrodynamic-Particle in Cell (MHD-PIC) simulations of ICME shocks to understand the acceleration mechanism. These shocks vary in Alfv\'enic Mach numbers as well as in magnetic field orientations (parallel \& quasi-perpendicular). We find that Diffusive Shock Acceleration (DSA) plays a significant role in accelerating particles in a parallel ICME shock. In contrast, Shock Drift Acceleration (SDA) plays a pivotal role in a quasi-perpendicular shock. High-Mach shocks are seen to accelerate particles more efficiently. Our simulations suggest that background turbulence and local particle velocity distribution around the shock can indirectly hint at the acceleration mechanism. Our results also point towards a few possible \textit{in situ} observations that could validate our understanding of the topic.

Context: In this paper we present a validation scheme to investigate the quality of coronal magnetic field models, which is based upon comparisons with observational data from multiple sources. Aims: Many of these coronal models may use a range of initial parameters that produce a large number of physically reasonable field configurations. However, that does not mean that these results are reliable and comply with the observations. With an appropriate validation scheme the quality of a coronal model can be assessed. Methods: The validation scheme is developed on the example of the EUropean Heliospheric FORecasting Information Asset (EUHFORIA) coronal model. For observational comparison we use EUV and white-light data to detect coronal features on the surface (open magnetic field areas) and off-limb (streamer and loop) structures from multiple perspectives (Earth view and the Solar Terrestrial Relations Observatory - STEREO). The validation scheme can be applied to any coronal model that produces magnetic field line topology. Results: We show its applicability by using that validation scheme on a large set of model configurations, which can be efficiently reduced to an ideal set of parameters that matches best with observational data. Conclusions: We conclude that by using a combined empirical visual classification with a mathematical scheme of topology metrics a very efficient and rather objective quality assessment for coronal models can be performed.

We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from the 1.3 mm continuum visibility profiles and the radial temperature profiles are estimated from H2CO rotation temperature maps. Molecular column densities are derived with the line fitting tool XCLASS. The physical and chemical properties are combined by applying the physical-chemical model MUSCLE in order to constrain the chemical timescales of a few line-rich cores. The morphology and spatial correlations of the molecular emission are analyzed using the HOG method. The mid-infrared data show that both regions contain a cluster of young stellar objects. Bipolar molecular outflows are observed in the CO 2-1 transition toward the strong mm cores indicating protostellar activity. We find strong molecular emission of SO, SiO, H2CO, and CH3OH in locations which are not associated with the mm cores. These shocked knots can be either associated with the bipolar outflows or, in the case of ISOSS J23053+5953, with a colliding flow that creates a large shocked region between the mm cores. The mean chemical timescale of the cores is lower (20 000 yr) compared to that of the sources of the more evolved CORE sample (60 000 yr). With the HOG method, we find that the spatial emission of species tracing the extended emission and of shock-tracing molecules are well correlated within transitions of these groups.

We present 97 new high-quality strong lensing candidates found in the final $\sim 350\,\rm deg^2$, that completed the full $\sim 1350\,\rm deg^2$ area of the Kilo-Degree Survey (KiDS). Together with our previous findings, the final list of high-quality candidates from KiDS sums up to 268 systems. The new sample is assembled using a new Convolutional Neural Network (CNN) classifier applied to $r$-band (best seeing) and $g,~r,~i$ color-composited images separately. This optimizes the complementarity of the morphology and color information on the identification of strong lensing candidates. We apply the new classifiers to a sample of luminous red galaxies (LRGs) and a sample of bright galaxies (BGs) and select candidates that received a high probability to be a lens from the CNN ($P_{\rm CNN}$). In particular, setting $P_{\rm CNN}>0.8$ for the LRGs, the $1$-band CNN predicts 1213 candidates, while the $3$-band classifier yields 1299 candidates, with only $\sim$30\% overlap. For the BGs, in order to minimize the false positives, we adopt a more conservative threshold, $P_{\rm CNN} >0.9$, for both CNN classifiers. This results in 3740 newly selected objects. The candidates from the two samples are visually inspected by 7 co-authors to finally select 97 "high-quality" lens candidates which received mean scores larger than 6 (on a scale from 0 to 10). We finally discuss the effect of the seeing on the accuracy of CNN classification and possible avenues to increase the efficiency of multi-band classifiers, in preparation of next-generation surveys from ground and space.

[abridged] Measurements of relative isotope abundances can provide unique insights into the formation and evolution histories of celestial bodies. The five stable isotopes of titanium are used to study the early history of the solar system and constrain Galactic chemical models. The minor isotopes of titanium are relatively abundant compared to those of other elements, making them more accessible for challenging observations. We assessed the feasibility of performing titanium isotope measurements in exoplanet atmospheres, and in particular, whether processing techniques used for high-resolution spectroscopy affect the derived isotope ratios. We used an archival high-dispersion CARMENES spectrum of the M-dwarf GJ 1002 as a proxy for an exoplanet observed at very high signal-to-noise. Spectral retrievals using petitRADTRANS models were performed on both narrow (7045-7090 {\AA}) and wide (7045-7500 {\AA}) wavelength regions, resulting in isotope ratios and uncertainties. These retrievals were repeated on the spectrum with its continuum removed to mimic typical high-dispersion exoplanet observations. The relative abundances of all minor Ti isotopes are found to be slightly enhanced compared to terrestrial values. Loss of continuum information from broadband spectral filtering has little effect on the isotope ratios. The CARMENES spectrum was subsequently degraded by adding varying levels of Gaussian noise to estimate the signal-to-noise requirements for future exoplanet atmospheric observations. For the wide wavelength range, a spectrum with signal-to-noise of 5 is required to determine the isotope ratios with relative errors $\lesssim$10%. Super Jupiters at large angular separations from their host star are the most accessible exoplanets, requiring about an hour of observing time on 8-meter-class telescopes, and less than a minute of observing time with the future Extremely Large Telescope.

A wide range of mechanisms have been put forward to explain the quenching of star formation in galaxies with cosmic time, however, the true balance of responsible mechanisms remains unknown. The identification and study of galaxies that have shut down their star formation on different timescales might elucidate which mechanisms dominate at different epochs and masses. Here we study the population of rapidly quenched galaxies (RQGs) in the SIMBA cosmological hydrodynamic simulation at $0.5<z<2$, comparing directly to observational post-starburst galaxies in the UKIDSS Ultra Deep Survey via their colour distributions and mass functions. We find that the fraction of quiescent galaxies that are rapidly quenched in SIMBA is 59% (or 48% in terms of stellar mass), which is higher than observed. A similar "downsizing" of RQGs is observed in both SIMBA and the UDS, with RQGs at higher redshift having a higher average mass. However, SIMBA produces too many RQGs at $1<z_q<1.5$ and too few low mass RQGs at $0.5<z_q<1$. The precise colour distribution of SIMBA galaxies compared to the observations also indicates various inconsistencies in star formation and chemical enrichment histories, including an absence of short, intense starbursts. Our results will help inform the next generation of galaxy evolution models, particularly with respect to the quenching mechanisms employed.

We are compiling a new list of gamma-ray jetted active galactic nuclei (AGN), starting from the fourth catalog of point sources of the Fermi Large Area Telescope (LAT). Our aim is to prepare a list of jetted AGN with known redshifts and classifications to be used to calibrate jet power. We searched in the available literature for all the published optical spectra and multiwavelength studies useful to characterize the sources. We found new, missed, or even forgotten information leading to a substantial change in the redshift values and classification of many sources. We present here the preliminary results of this analysis and some statistics based on the gamma-ray sources with right ascension within the interval $0^{\rm h}-12^{\rm h}$ (J2000). Although flat-spectrum radio quasars and BL Lac objects are still the dominant populations, there is a significant increase in the number of other objects, such as misaligned AGN, narrow-line Seyfert 1 galaxies, and Seyfert galaxies. We also introduced two new classes of objects: changing-look AGN and ambiguous sources. About one third of the sources remain unclassified.

Multiband photometric transit observations or low-resolution spectroscopy (spectro-photometry) are normally used to retrieve the broadband transmission spectra of transiting exoplanets in order to assess the chemical composition of their atmospheres. In this paper, we present an alternative approach for recovering the broadband transmission spectra using chromatic Doppler tomography based on physical modeling through the SOAP tool: CHOCOLATE (CHrOmatiC line prOfiLe tomogrAphy TEchnique). To validate the method and examine its performance, we use observational data recently obtained with the ESPRESSO instrument to retrieve the transmission spectra of the archetypal hot Jupiter HD 209458b. Our findings indicate that the recovered transmission spectrum is in good agreement with the results presented in previous studies, which used different methodologies to extract the spectrum, achieving similar precision. We explored several atmospheric models and inferred from spectral retrieval that a model containing H2O and NH3 is the preferred scenario. The CHOCOLATE methodology is particularly interesting for future studies of exoplanets around young and active stars or moderate to fast rotating stars, considering SOAP's ability to model stellar active regions and the fact that the rotational broadening of spectral lines favors its application. Furthermore, CHOCOLATE will allow the broad transmission spectrum of a planet to be retrieved using high signal-to-noise ratio, high-resolution spectroscopy with the next generation of Extremely Large Telescopes (ELTs), where low-resolution spectroscopy will not always be accessible.

The eastern region of the Small Magellanic Cloud (SMC) is found to have a foreground stellar sub-structure, which is identified as a distance bimodality (12 kpc apart) in the previous studies using Red Clump (RC) stars. Interestingly, studies of Red giant branch (RGB) stars in the eastern SMC indicate a bimodal radial velocity (RV) distribution. In this study, we investigate the connection between these two bimodal distributions to better understand the nature and origin of the foreground stellar sub-structure in the eastern SMC. We use the Gaia EDR3 astrometric data and archival RV data of RGB stars for this study. We found a bimodal RV distribution of RGB stars (separated by 35 - 45 km/s) in the eastern and south-western (SW) outer regions. The observed proper motion values of the lower and higher RV RGB components in the eastern regions are similar to those of the foreground and main-body RC stars respectively. This suggests that the two RGB populations in the eastern region are separated by a similar distance as those of the RC stars, and the RGB stars in the lower RV component are part of the foreground sub-structure. Based on the differences in the distance and RV of the two components, we estimated an approximate time of formation of this sub-structure as 307+/-65 Myr ago. This is comparable with the values predicted by simulations for the recent epoch of tidal interaction between the Magellanic Clouds. Comparison of the observed properties of RGB stars, in the outer SW region, with N-body simulations shows that the higher RV component in the SW region is at a farther distance than the main body, indicating the presence of a stellar Counter-Bridge in the SW region of the SMC.

Aims. We study the X-ray emission of the galactic supernova remnant (SNR) G18.95-1.1 with the eROSITA telescope on board the Spectrum Rentgen Gamma (SRG) orbital observatory. In addition to the pulsar wind nebula that was previously identified and examined by ASCA and Chandra, we study the X-ray spectra of the bright SNR ridge, which is resolved into a few bright clumps. Methods. The wide field of view and the large collecting area in the 0.2-2.3 keV energy range of SRG/eROSITA allowed us to perform spatially resolved spectroscopy of G18.95-1.1. Results. The X-ray ridge of G18.95-1.1 is asymmetric, indicating either supernova ejecta asymmetry or their interaction with a cloud. The X-ray dim northern regions outside the pulsar wind nebula can be described by a thin thermal plasma emission with a temperature ~0.3 keV and a solar composition. The X-ray spectra of a few bright clumps located along the southern ridge may be satisfactorily approximated by a single thermal component of the Si-rich ejecta at the collisional ionization equilibrium with a temperature of about 0.3 keV. The bright ridge can be alternatively fit with a single component that is not dominated by equilibrium ejecta with T ~ 0.6 keV. The high ratio of the derived Si/O abundances indicates that the ejecta originated in deep layers of the progenitor star. The plasma composition of a southern Si-rich clump and the bright ridge are similar to what was earlier found in the Vela shrapnel A and G.

We present the results of photometric reverberation mapping observations on the changing look active galactic nucleus Mrk 590 at z = 0.026. The observations were carried out from July to December, 2018 using broad band B, R and narrow band H{\alpha} and S II filters. B-band traces the continuum emission from the accretion disk, R-band encompasses both the continuum emission from the accretion disk and the redshifted H{\alpha} line from the broad line region (BLR), S II band contains the redshifted H{\alpha} emission and the H{\alpha} band traces the continuum emission underneath the S II band. All the light curves showed strong variation with a fractional root-mean-square variation of 0.132 (+/-) 0.001 in B-band and 0.321 (+/-) 0.001 in H{\alpha} line. From cross-correlation function analysis, we obtained a delayed response of H{\alpha} line emission to the optical B-band continuum emission of 21.44(+1.49/-2.11) days in the rest-frame of the source, corresponding to a linear size of the BLR of 0.018 pc. This is consistent with previous estimates using H{\beta}. By combining the BLR size with the H{\alpha} line full width at half maximum of 6478 (+/-) 240 km/s measured from a single-epoch spectrum obtained with the Subaru telescope, we derived a black hole mass of 1.96 (+0.15 / -0.21) X 10^8 Msun.

In large scale maps of the Galactic disc, the Carina Arm stands out as a clear spiral feature, hosting prominent star clusters and associations rich in massive stars. We study the proper motions of 4199 O and early B most likely in the far Carina Arm, at distances mainly in excess of 4 kpc from the Sun, within the sky region, $282^{\circ} < \ell < 294^{\circ}$ and $-3^{\circ} < b < +1^{\circ}$ (Galactic coordinates). The sample is constructed by extending an existing blue-selected catalogue, and cross-matching with Gaia EDR3 astrometry. The observed pattern of proper motions is modulated into a saw-tooth pattern, with full amplitude approaching 1 mas yr$^{-1}$, recurring roughly every 2--3 degrees of longitude (200--300 pc at the median OB-star distance of 5.8 kpc). Kinematic perturbation of underlying circular rotation is most likely present. The data also reveal a moving group containing $>50$ OB stars at $\ell \sim 286^{\circ}$, $b \sim -1^{\circ}.4$ behind the main run of the far arm. An analysis of relative proper motions is performed that yields an incidence of runaway O stars of at least 10\% (potentially $>20$\% when full space motions become available). To map where runaways have run away from, we set up simulations for the region that assume linear trajectories and test for trajectory impact parameter in order to identify likely ejection hot spots. We find the method currently gives good results for times of flight of up to $\sim$4 Myr. It shows convincingly that only NGC 3603 and Westerlund 2 have ejected OB stars in significant numbers. Indeed, both clusters have experienced intense spells of ejection between 0.6--0.9 and 0.5--0.8 Myr ago, respectively.

Knowledge of X-ray shock and radio relic connection in merging galaxy clusters has been greatly extended in terms of both observation and theory over the last decade. ZwCl 2341+0000 is a double-relic merging galaxy cluster; previous studies have shown that half of the S relic is associated with an X-ray surface brightness discontinuity, while the other half not. The discontinuity was believed to be a shock front. Therefore, it is a mysterious case of an only partial shock-relic connection. By using the 206.5 ks deep Chandra observations, we aim to investigate the nature of the S discontinuity. Meanwhile, we aim to explore new morphological and thermodynamical features. In addition, we utilize the GMRT and JVLA images to compute radio spectral index (SI) maps. In the deep observations, the previously reported S surface brightness discontinuity is better described as a sharp change in slope or as a kink, which is likely contributed by the disrupted core of the S subcluster. The radio SI maps show spectral flattening at the SE edge of the S relic, suggesting that the location of the shock front is 640 kpc away from the kink. We update the radio shock Mach number to be $2.2\pm0.1$ and $2.4\pm0.4$ for the S and N radio relics, respectively, based on the injection SI. We also put a 3 sigma lower limit on the X-ray Mach number of the S shock to be >1.6. Meanwhile, the deep observations reveal that the N subcluster is in a perfect cone shape, with a ~400 kpc linear cold front on each side. This type of conic subcluster has been predicted by simulations but is observed here for the first time. It represents a transition stage between a blunt-body cold front and a slingshot cold front. Strikingly, we found a 400 kpc long gas trail attached to the apex of the cone, which could be due to the gas stripping. In addition, an over-pressured hot region is found in the SW flank of the cluster.

We investigate the structural properties, distribution and abundance of LCDM dark matter subhaloes using the Phi-4096 and Uchuu suite of N-body cosmological simulations. Thanks to the combination of their large volume, high mass resolution and superb statistics, we are able to quantify -- for the first time consistently over more than seven decades in ratio of subhalo-to-host-halo mass -- dependencies of subhalo properties with mass, maximum circular velocity, Vmax, host halo mass and distance to host halo centre. We also dissect the evolution of these dependencies over cosmic time. We provide accurate fits for the subhalo mass and velocity functions, both exhibiting decreasing power-law slopes in the expected range of values and with no significant dependence on redshift. We also find subhalo abundance to depend weakly on host halo mass. We explore the distribution of subhaloes within their hosts and its evolution over cosmic time for subhaloes located as deep as ~0.1 per cent of the host virial radius. Subhalo structural properties are codified via a concentration parameter, cV, that does not depend on any specific, pre-defined density profile and relies only on Vmax. We derive the cV-Vmax relation in the range 7-1500 km/s and find an important dependence on distance of the subhalo to the host halo centre, as already described in Molin\'e et al. (2017). Interestingly, we also find subhaloes of the same mass to be significantly more concentrated into more massive hosts. Finally, we investigate the redshift evolution of cV, and provide accurate fits that take into account all mentioned dependencies. Our results offer an unprecedented detailed characterization of the subhalo population, consistent over a wide range of subhalo and host halo masses, as well as cosmic times. Our work enables precision work in any future research involving dark matter halo substructure.

Context: Open clusters (OCs) are widely used as test particles to investigate a variety of astrophysical phenomena, from stellar evolution to Galactic evolution. Gaia and the complementary massive spectroscopic surveys are providing an unprecedented wealth of information about these systems. Aims: The Open Cluster Chemical Abundances from Spanish Observatories (OCCASO) survey aims to complement all this work by determining OCs accurate radial velocities and chemical abundances from high-resolution, R$\geq$60\,000, spectra. Methods: Radial velocities have been obtained by cross-correlating the observed spectra with a library of synthetic spectra which covers from early M to A spectral types. Results: We provide radial velocities for 336 stars, including several Gaia Benchmark Stars and objects belonging to 51 open clusters. The internal uncertainties of the derived radial velocities go from 10 m/s to 21 m/s as a function of the instrumental configuration used. The derived radial velocities, together with the Gaia proper motions, have been used to investigate the cluster membership of the observed stars. After this careful membership analysis, we obtain average velocities for 47 open clusters. To our knowledge, this is the first radial velocity determination for 5 of these clusters. Finally, the radial velocities, proper motions, distances and ages have been used to investigate the kinematics of the observed clusters and in the integration of their orbits.

Major galaxy mergers can trigger nuclear activities and are responsible for high-luminosity quasi-stellar objects /active galactic nuclei (QSOs/AGNs). In certain circumstances, such mergers may cause dual active galactic nuclei (dAGN) phenomenon. This study investigates dAGN triggering and evolution of massive black holes (MBHs) during the merging processes using hydrodynamic code GADGET-2 to simulate several gas-rich major mergers at redshift $z=2$ and $3$, respectively. Results reveal that gas-rich major mergers can trigger significant nuclear activities after the second and third pericentric passages and the formation of dAGN with significant time duration ($\sim 10 - 390$\,Myr). During the merging processes, galactic bulge evolves with time because of the rapid star formation in each (or both) galactic centers and initial mixing of stars in galactic disks due to violent relaxation. MBHs grow substantially due to accretion and finally merge into a bigger black hole. The growth of galactic bulges and corresponding increases of its velocity dispersions predate the growth of MBHs in the dAGN stages. The MBHs in these stages deviate below the relation between MBH mass and bulge mass (or velocity dispersion), and they revert to the relation after the final mergers due to the significant accretion that occurs mostly at a separation less than a few kpc. Then, the two MBHs merge with each other.

The Radio Arc is a system of organized non-thermal filaments (NTFs) located within the Galactic Center (GC) region of the Milky Way. Recent observations of the Radio Arc NTFs revealed a magnetic field which alternates between being parallel and rotated with respect to the orientation of the filaments. This pattern is in stark contrast to the predominantly parallel magnetic field orientations observed in other GC NTFs. To help elucidate the origin of this pattern, we analyze spectro-polarimetric data of the Radio Arc NTFs using an Australian Telescope Compact Array data set covering the continuous frequency range from $\sim$4 to 11 GHz at a spectral resolution of 2 MHz. We fit depolarization models to the spectral polarization data to characterize Faraday effects along the line-of-sight. We assess whether structures local to the Radio Arc NTFs may contribute to the unusual magnetic field orientation. External Faraday effects are identified as the most likely origin of the rotation observed for the Radio Arc NTFs; however, internal Faraday effects are also found to be likely in regions of parallel magnetic field. The increased likelihood of internal Faraday effects in parallel magnetic field regions may be attributed to the effects of structures local to the GC. One such structure could be the Radio Shell local to the Radio Arc NTFs. Future studies are needed to determine whether this alternating magnetic field pattern is present in other multi-stranded NTFs, or is a unique property resulting from the complex interstellar region local to the Radio Arc NTFs.

An overwhelming majority of galactic spiral arms trail with respect to the rotation of the galaxy, though a small sample of leading spiral arms has been observed. The formation of these leading spirals is not well understood. Here we show, using collisionless $N$-body simulations, that a barred disc galaxy in a retrograde dark matter halo can produce long-lived ($\sim3$ Gyr) leading spiral arms. Due to the strong resonant coupling of the disc to the halo, the bar slows rapidly and spiral perturbations are forced ahead of the bar. We predict that such a system, if observed, will also host a dark matter wake oriented perpendicular to the stellar bar. More generally, we propose that any mechanism that rapidly decelerates the stellar bar will allow leading spiral arms to flourish.

Super-soft X-ray sources were established as a heterogeneous class of objects from observations of the Large Magellanic Cloud (LMC). We have searched for new sources of this class in the X-ray images obtained from the XMM-Newton survey of the LMC and additional archival observations. We first selected candidates by visual inspection of the image, and screened out artefacts which can mimic super-soft X-ray sources as well as bright foreground stars which create optical loading on the CCD image. We finally obtained 4 new super-soft X-ray sources for which we performed detailed X-ray timing and spectral analysis and searched for possible optical counterparts to identify their nature. XMMU J050452.0-683909 is identified as the central star of the planetary nebula SMP LMC21 in the LMC. We suggest XMMU J051854.8-695601 and XMMU J050815.1-691832 as new soft intermediate polars from the nature of their X-ray spectrum. Their estimated absorption-corrected luminosities and the blackbody radii indicate that they are located in our Galaxy, rather than the LMC. We discovered coherent pulsations of 497 s from XMMU J044626.6-692011 which indicates a magnetic cataclysmic variable nature of the source. The location of XMMU J044626.6-692011 in the LMC or our Galaxy is less clear. It could either be a white dwarf in the LMC with nuclear burning on its surface near the Eddington limit, or another soft intermediate polar in our Galaxy. The discovery of new super-soft X-ray sources makes a significant contribution to the known population in our own Galaxy. An observed higher density of sources in the direction of the Magellanic Clouds can likely be explained by the relatively low Galactic column density in their direction as well as a large number of existing observations sensitive at low X-ray energies.

One of the grand challenges in heliophysics is the characterisation of coronal mass ejection (CME) magnetic structure and evolution from eruption at the Sun through heliospheric propagation. At present, the main difficulties are related to the lack of direct measurements of the coronal magnetic fields and the lack of 3D in-situ measurements of the CME body in interplanetary space. Nevertheless, the evolution of a CME magnetic structure can be followed using a combination of multi-point remote-sensing observations and multi-spacecraft in-situ measurements as well as modelling. Accordingly, we present in this work the analysis of two CMEs that erupted from the Sun on 28 April 2012. We follow their eruption and early evolution using remote-sensing data, finding indications of CME--CME interaction, and then analyse their interplanetary counterpart(s) using in-situ measurements at Venus, Earth, and Saturn. We observe a seemingly single flux rope at all locations, but find possible signatures of interaction at Earth, where high-cadence plasma data are available. Reconstructions of the in-situ flux ropes provide almost identical results at Venus and Earth but show greater discrepancies at Saturn, suggesting that the CME was highly distorted and/or that further interaction with nearby solar wind structures took place before 10 AU. This work highlights the difficulties in connecting structures from the Sun to the outer heliosphere and demonstrates the importance of multi-spacecraft studies to achieve a deeper understanding of the magnetic configuration of CMEs.

We present a cosmological solution to the electroweak hierarchy problem. After discussing general features of cosmological approaches to naturalness, we extend the Standard Model with two light scalars very weakly coupled to the Higgs and present the mechanism, which we recently introduced in a companion paper to explain jointly the electroweak hierarchy and the strong-CP problem. In this work we show that this solution can be decoupled from the strong-CP problem and discuss its possible implementations and phenomenology. The mechanism works with any standard inflationary sector, it does not require weak-scale inflation or a large number of e-folds, and does not introduce ambiguities related to eternal inflation. The cutoff of the theory can be as large as the Planck scale, both for the Cosmological Constant and for the Higgs sector. Reproducing the observed dark matter relic density fixes the couplings of the two new scalars to the Standard Model, offering a target to future axion or fifth force searches. Depending on the specific interaction of the scalars with the Standard Model, the mechanism either yields rich phenomenology at colliders or provides a novel joint solution to the strong-CP problem. We highlight what predictions are common to most realizations of cosmological selection of the weak scale and will allow to test this general framework in the near future.

This paper tries to obtain a simple picture of several aspects of the mode structure in relativistic non-equilibrium thermodynamics. Its pedagogical focus is on the relation between long-wavelength perturbation modes of the causal M\"{u}ller-Israel-Stewart (MIS) theory and those of the traditional Eckart theory. Principally, this issue was clarified in a series of papers by Hiscock and Lindblom (see [8-10]). Here, I put together some essential features which do not require the entire formalism of the complete theory.

Rapid sky localization of gravitational wave sources is crucial to enable prompt electromagnetic follow-ups. In this article, we present a novel semi-analytical approach for sky localization of gravitational waves from compact binary coalescences. We use the Bayesian framework with an analytical approximation to the prior distributions for a given astrophysical model. We derive a semi-analytical solution to the posterior distribution of source directions. This method only requires one-fold numerical integral that marginalizes over the merger time, compared to the five-fold numerical integration otherwise needed in the Bayesian localization method. The performance of the method is demonstrated using a set of binary neutron stars (BNS) injections on Gaussian noise using LIGO-Virgo's design and O2 sensitivity. We find the median of 90% confidence area in O2 sensitivity to be $\mathcal{O}(10^2) ~\mathrm{deg}^2$, comparable to that of the existing LIGO-Virgo online localization method Bayestar and parameter estimation toolkit LALInference. In the end, we apply this method to localize the BNS event GW170817 and find the 50% (90%) confidence region of 11 $\mathrm{deg}^2$ (50 $\mathrm{deg}^2$). The detected optical counterpart of GW170817 resides within our 50% confidence area.

While biologists have not yet reached a consensus on the definition of life, homochirality - the specific molecular handedness of biomolecules - is a phenomenon only produced by life. The unraveling of its origin requires interdisciplinary research, by exploring fundamental physics, chemistry, astrophysics and biology. Here, we consider the origin of biological homochirality in the context of astrophysics and particle physics. The weak force, one of the fundamental forces operating in nature, is parity-violating. On Earth, at ground level, most of our cosmic radiation dose comes from polarized muons formed in a decay involving the weak force. We discuss how the magnetic polarization is transmitted in cosmic showers within several different environments which are prime targets in the search for the origin of life. We consider how this polarization could have induced a biological preference for one type of chirality over the other, and discuss the implications for the search of life in other worlds.

We discuss stiffening of dense matter in two color QCD (QC$_2$D) where hadrons are mesons and diquark baryons. We study two models which describe a transition of matter from the Bose-Einstein-Condensation regime at low density to the Bardeen-Cooper-Schrieffer regime at high density. The first model is based on coherent states of diquarks, and the second is the Nambu-Jona-Lasinio model with diquark pairing terms. We particularly focus on how quark states are occupied as baryon density increases. We find that, due to the occupied quark levels, the ideal gas picture of diquarks breaks down at density significantly less than the density where baryon cores overlap. The saturation of quark states at low momenta stiffens equations of state. We also study the effects of interactions which depend on the quark occupation probability. We argue that equations of state become very stiff when the bulk part of the quark Fermi sea has the effective repulsion but the Fermi surface enjoys the attractive correlations. This disparity for different momentum domains is possible due to the strong channel dependence in gluon exchanges with momentum transfer of $0.2-1$ GeV. These concepts can be transferred from QC$_2$D to QCD in any numbers of colors.