We have implemented prescriptions for modelling pulsars in the rapid binary population synthesis code COMPAS. We perform a detailed analysis of the double neutron star (DNS) population, accounting for radio survey selection effects. The surface magnetic field decay timescale (${\sim}1000$\,Myr) and mass scale (${\sim}0.02$\,M$_\odot$) are the dominant uncertainties in our model. Mass accretion during common envelope evolution plays a non-trivial role in recycling pulsars. We find a best-fit model that is in broad agreement with the observed Galactic DNS population. Though the pulsar parameters (period and period derivative) are strongly biased by radio selection effects, the observed orbital parameters (orbital period and eccentricity) closely represent the intrinsic distributions. The number of radio observable DNSs in the Milky Way at present is about 2500 in our model, corresponding to approximately 10\% of the predicted total number of DNSs in the galaxy. Using our model calibrated to the Galactic DNS population, we make predictions for DNS mergers observed in gravitational waves. The DNS chirp mass distribution varies from 1.1M$_\odot$ to 2.1M$_\odot$ and the median is found to be 1.14\,M$_\mathrm{\odot}$. The expected effective spin $\chi_\mathrm{eff}$ for isolated DNSs is $\lesssim$0.03 from our model. We predict that 34\% of the current Galactic isolated DNSs will merge within a Hubble time, and have a median total mass of 2.7\,M$_\mathrm{\odot}$. Finally, we discuss implications for fast radio bursts and post-merger remnant gravitational-waves.

Context. Moderately metal-poor inner bulge globular clusters are relics of a generation of long-lived stars that formed in the early Galaxy. Terzan 9, projected at 4d 12 from the Galactic center, is among the most central globular clusters in the Milky Way, showing an orbit which remains confined to the inner 1 kpc. Aims. Our aim is the derivation of the cluster's metallicity, together with an accurate measurement of the mean radial velocity. In the literature, metallicities in the range between have been estimated for Terzan 9 based on color-magnitude diagrams and CaII triplet (CaT) lines. Aims. Our aim is the derivation of the cluster's metallicity, together with an accurate measurement of the mean radial velocity. In the literature, metallicities in the range between -2.0 and -1.0 have been estimated for Terzan 9 based on color-magnitude diagrams and CaII triplet (CaT) lines. Methods. Given its compactness, Terzan 9 was observed using the Multi Unit Spectroscopic Explorer (MUSE) at the Very Large Telescope. The extraction of spectra from several hundreds of individual stars allowed us to derive their radial velocities, metallicities, and [Mg/Fe]. The spectra obtained with MUSE were analysed through full spectrum fitting using the ETOILE code. Results. We obtained a mean metallicity of [Fe/H] -1.10 0.15, a heliocentric radial velocity of vhr = 58.1 1.1 km/s , and a magnesium-to-iron [Mg/Fe] = 0.27 0.03. The metallicity-derived character of Terzan 9 sets it among the family of the moderately metal-poor Blue Horizontal Branch clusters HP 1, NGC 6558, and NGC 6522.

One of the most exciting scientific challenges is to detect Earth-like planets in the habitable zones of other stars in the galaxy and search for evidence of life. The ability to observe and characterise dozens of potentially Earth-like planets now lies within the realm of possibility due to rapid advances in key space and imaging technologies. The associated challenge of directly imaging very faint planets in orbit around nearby very bright stars is now well understood, with the key instrumentation also being perfected and developed. Such advances will allow us to develop large transformative telescopes, covering a broad UV-optical-IR spectral range, which can carry out the detailed research programmes designed to answer the questions we wish to answer: Carry out high contrast imaging surveys of nearby stars to search for planets within their habitable zones. Characterise the planets detected to determine masses and radii from photometric measurements. Through spectroscopic studies of their atmospheres and surfaces, search for habitability indicators and for signs of an environment that has been modified by the presence of life. Active studies of potential missions have been underway for a number of years. The latest of these is the Large UV Optical IR space telescope (LUVOIR), one of four flagship mission studies commissioned by NASA in support of the 2020 US Decadal Survey. LUVOIR, if selected, will be of interest to a wide scientific community and will be the only telescope capable of searching for and characterizing a sufficient number of exoEarths to provide a meaningful answer to the question - Are we alone?. This paper is a submission to the European Space Agency Voyage 2050 call for white papers outlining the case for an ESA contribution to a Large UVOIR telescope.

The polarization state of starlight reflected by a planetary atmosphere uniquely reveals coverage, particle size, and composition of aerosols as well as changing cloud patterns. It is not possible to obtain a comparable level of detailed from flux-only observations. Furthermore, polarization observations can probe the atmosphere of planets independently of the orbital geometry (i.e., transiting and non-transiting planets). We show that a high-resolution spectropolarimeter with a broad wavelength coverage, particularly if attached to a large space telescope, would enable simultaneous study of the polarimetric exoplanet properties of the continuum and to look for and characterize the polarimetric signal due to scattering from single molecules.

NGC 6558 is a bulge globular cluster with a blue horizontal branch (BHB), combined with a metallicity of [Fe/H]~-1.0. It is similar to HP 1 and NGC 6522, which could be among the oldest objects in the Galaxy. Element abundances in these clusters could reveal the nature of the first supernovae. We aim to carry out detailed spectroscopic analysis for four red giants of NGC 6558, in order to derive the abundances of the light elements C, N, O, Na, Al, the alpha-elements Mg, Si, Ca, Ti, and the heavy elements Y, Ba, and Eu. High-resolution spectra of four stars with FLAMES-UVES@VLT UT2-Kueyen were analysed. Spectroscopic parameter-derivation was based on excitation and ionization equilibrium of FeI and FeII. This analysis results in a metallicity of [Fe/H] = -1.17+-0.10 for NGC 6558. We find the expected alpha-element enhancements in O and Mg with [O/Fe]=+0.40, [Mg/Fe]=+0.33, and low enhancements in Si and Ca. Ti has a moderate enhancement of [Ti/Fe]=+0.22. The r-element Eu appears very enhanced with a mean value of [Eu/Fe]=+0.63. Ba appears to have a solar abundance ratio relative to Fe. NGC 6558 shows an abundance pattern that could be typical of the oldest inner bulge globular clusters, together with the pattern in the similar clusters NGC 6522 and HP 1. They show low abundances of the odd-Z elements Na and Al, and of the explosive nucleosynthesis alpha-elements Si, Ca, and Ti. The hydrostatic burning alpha-elements O and Mg are normally enhanced as expected in old stars enriched with yields from core-collapse supernovae, and the iron-peak elements Mn, Cu, Zn show low abundances, which is expected for Mn and Cu, but not for Zn. Finally, the cluster trio NGC 6558, NGC 6522, and HP 1 show similar abundance patterns.

Context: Detecting moons around exoplanets is a major goal of current and future observatories. Moons are suspected to influence rocky exoplanet habitability, and gaseous exoplanets in stellar habitable zones could harbour abundant and diverse moons to target in the search for extraterrestrial habitats. Exomoons will contribute to exoplanetary signals but are virtually undetectable with current methods. Aims: We identify and analyse traces of exomoons in the temporal variation of total and polarised fluxes of starlight reflected by an Earth-like exoplanet and its spatially unresolved moon across all phase angles, with both orbits viewed in an edge-on geometry. Methods: We compute the total and linearly polarised fluxes, and the degree of linear polarization P of starlight that is reflected by the exoplanet with its moon along their orbits, accounting for the temporal variation of the visibility of the planetary and lunar disks, and including effects of mutual transits and mutual eclipses. Our computations pertain to a wavelength of 450 nm. Results: Total flux F shows regular dips due to planetary and lunar transits and eclipses. Polarization P shows regular peaks due to planetary transits and lunar eclipses, and P can increase and/or slightly decrease during lunar transits and planetary eclipses. Changes in F and P will depend on the radii of the planet and moon, on their reflective properties, and their orbits, and are about one magnitude smaller than the smooth background signals. The typical duration of a transit or an eclipse is a few hours. Conclusions: Traces of an exomoon due to planetary and lunar transits and eclipses show up in F and P of sunlight reflected by planet-moon systems and could be searched for in exoplanet flux and/or polarisation phase functions.

The circular polarization of light that planets reflect is often neglected because it is very small compared to the linear polarization. It could, however, provide information on a planet's atmosphere and surface, and on the presence of life, because homochiral molecules that are the building blocks of life on Earth are known to reflect circularly polarized light. We compute $P_c$, the degree of circular polarization, for light that is reflected by rocky (exo)planets with liquid water or sulfuric acid solution clouds, both spatially resolved across the planetary disk and, for planets with patchy clouds, integrated across the planetary disk, for various planetary phase angles $\alpha$. The optical thickness and vertical distribution of the atmospheric gas and clouds, the size parameter and refractive index of the cloud particles, and $\alpha$ all influence $P_c$. Spatially resolved, $P_c$ varies between $\pm 0.20\%$ (the sign indicates the polarization direction). Only for small gas optical thicknesses above the clouds do significant sign changes (related to cloud particle properties) across the planets' hemispheres occur. For patchy clouds, the disk--integrated $P_c$ is typically smaller than $\pm 0.025\%$, with maximums for $\alpha$ between $40^\circ$ and $70^\circ$, and $120^\circ$ to $140^\circ$. As expected, the disk--integrated $P_c$ is virtually zero at $\alpha=0^\circ$ and 180$^\circ$. The disk--integrated $P_c$ is also very small at $\alpha \approx 100^\circ$. Measuring circular polarization signals appears to be challenging with current technology. The small atmospheric circular polarization signal could, however, allow the detection of circular polarization due to homochiral molecules. Confirmation of the detectability of such signals requires better knowledge of the strength of circular polarization signals of biological sources.

PyMieDAP (the Python Mie Doubling-Adding Programme) is a Python--based tool for computing the total, linearly, and circularly polarized fluxes of incident unpolarized sun- or starlight that is reflected by, respectively, Solar System planets or moons, or exoplanets at a range of wavelengths. The radiative transfer computations are based on an adding--doubling Fortran algorithm and fully include polarization for all orders of scattering. The model (exo)planets are described by a model atmosphere composed of a stack of homogeneous layers containing gas and/or aerosol and/or cloud particles bounded below by an isotropically, depolarizing surface (that is optionally black). The reflected light can be computed spatially--resolved and/or disk--integrated. Spatially--resolved signals are mostly representative for observations of Solar System planets (or moons), while disk--integrated signals are mostly representative for exoplanet observations. PyMieDAP is modular and flexible, and allows users to adapt and optimize the code according to their needs. PyMieDAP keeps options open for connections with external programs and for future additions and extensions. In this paper, we describe the radiative transfer algorithm that PyMieDAP is based on and the code's principal functionalities. And we provide benchmark results of PyMieDAP that can be used for testing its installation and for comparison with other codes. PyMieDAP is available online under the GNU GPL license at http://gitlab.com/loic.cg.rossi/pymiedap

We present orbit analysis for a sample of 8 inner bulge globular clusters, together with one reference halo object. We used proper motion values derived from long time base CCD data. Orbits are integrated in both an axisymmetric model and a model including the Galactic bar potential. The inclusion of the bar proved to be essential for the description of the dynamical behavior of the clusters. We use the Monte Carlo scheme to construct the initial conditions for each cluster, taking into account the uncertainties in the kinematical data and distances. The sample clusters show typically maximum height to the Galactic plane below 1.5 kpc, and develop rather eccentric orbits. Seven of the bulge sample clusters share the orbital properties of the bar/bulge, having perigalactic and apogalatic distances, and maximum vertical excursion from the Galactic plane inside the bar region. NGC 6540 instead shows a completely different orbital behaviour, having a dynamical signature of the thick-disc. Both prograde and prograde-retrograde orbits with respect to the direction of the Galactic rotation were revealed, which might characterize a chaotic behaviour.

Flight height is a fundamental parameter for correcting the gamma signal produced by terrestrial radionuclides measured during airborne surveys. The frontiers of radiometric measurements with UAV require light and accurate altimeters flying at some 10 m from the ground. We equipped an aircraft with seven altimetric sensors (three low-cost GNSS receivers, one inertial measurement unit, one radar altimeter and two barometers) and analyzed $\sim$ 3 h of data collected over the sea in the (35-2194) m altitude range. At low altitudes (H $<$ 70 m) radar and barometric altimeters provide the best performances, while GNSS data are used only for barometer calibration as they are affected by a large noise due to the multipath from the sea. The $\sim$ 1 m median standard deviation at 50 m altitude affects the estimation of the ground radioisotope abundances with an uncertainty less than 1.3%. The GNSS double-difference post-processing enhanced significantly the data quality for H $>$ 80 m in terms of both altitude median standard deviation and agreement between the reconstructed and measured GPS antennas distances. Flying at 100 m the estimated uncertainty on the ground total activity due to the uncertainty on the flight height is of the order of 2%.

Context. Clouds have already been detected in exoplanetary atmospheres. They play crucial roles in a planet's atmosphere and climate and can also create ambiguities in the determination of atmospheric parameters such as trace gas mixing ratios. Knowledge of cloud properties is required when assessing the habitability of a planet. Aims. We aim to show that various types of cloud cover such as polar cusps, subsolar clouds, and patchy clouds on Earth-like exoplanets can be distinguished from each other using the polarization and flux of light that is reflected by the planet. Methods. We have computed the flux and polarization of reflected starlight for different types of (liquid water) cloud covers on Earth-like model planets using the adding-doubling method, that fully includes multiple scattering and polarization. Variations in cloud-top altitudes and planet-wide cloud cover percentages were taken into account. Results. We find that the different types of cloud cover (polar cusps, subsolar clouds, and patchy clouds) can be distinguished from each other and that the percentage of cloud cover can be estimated within 10%. Conclusions. Using our proposed observational strategy, one should be able to determine basic orbital parameters of a planet such as orbital inclination and estimate cloud coverage with reduced ambiguities from the planet's polarization signals along its orbit.

Earth-like, potentially habitable exoplanets are prime targets in the search for extraterrestrial life. Information about their atmosphere and surface can be derived by analyzing light of the parent star reflected by the planet. We investigate the influence of the surface albedo $A_{\rm s}$, the optical thickness $b_{\rm cloud}$ and altitude of water clouds, and the mixing ratio $\eta$ of biosignature O$_2$ on the strength of the O$_2$ A-band (around 760 nm) in flux and polarization spectra of starlight reflected by Earth-like exoplanets. Our computations for horizontally homogeneous planets show that small mixing ratios ($\eta$ < 0.4) will yield moderately deep bands in flux and moderate to small band strengths in polarization, and that clouds will usually decrease the band depth in flux and the band strength in polarization. However, cloud influence will be strongly dependent on their properties such as optical thickness, top altitude, particle phase, coverage fraction, horizontal distribution. Depending on the surface albedo, and cloud properties, different O$_2$ mixing ratios $\eta$ can give similar absorption band depths in flux and band strengths in polarization, in particular if the clouds have moderate to high optical thicknesses. Measuring both the flux and the polarization is essential to reduce the degeneracies, although it will not solve them, in particular not for horizontally inhomogeneous planets. Observations at a wide range of phase angles and with a high temporal resolution could help to derive cloud properties and, once those are known, the mixing ratio of O$_2$ or any other absorbing gas.

The evolution of star clusters is largely affected by the tidal field generated by the host galaxy. It is thus in principle expected that under the assumption of an "universal" initial cluster mass function the properties of the evolved present-day mass function of star cluster systems should show a dependency on the properties of the galactic environment in which they evolve. To explore this expectation a sophisticated model of the tidal field is required in order to study the evolution of star cluster systems in realistic galaxies. Along these lines, in the present work we first describe a method developed for coupling $N$-body simulations of galaxies and star clusters. We then generate a database of galaxy models along the Hubble sequence and calibrate evolutionary equations to the results of direct $N$-body simulations of star clusters in order to predict the clusters' mass evolution as function of the galactic environment. We finally apply our methods to explore the properties of evolved "universal" initial cluster mass functions and any dependence on the host galaxy morphology and mass distribution. The preliminary results show that an initial power-law distribution of the masses "universally" evolves into a log-normal distribution, with the properties correlated with the stellar mass and stellar mass density density of the host galaxy.

We computed proper motions of a selected sample of globular clusters projected on the central bulge, employing CCD images gathered along the last 25 years at the ESO-NTT, ESO-Danish and HST telescopes. We presented a method to derive their proper motions, and a set of coordinate transformations to obtain 3D Galactic velocity vectors of the clusters. We analysed 10 globular clusters, namely Terzan 1, Terzan 2, Terzan 4, Terzan 9, NGC 6522, NGC 6558, NGC 6540, AL~3,ESO456--SC38 and Palomar 6. For comparison purposes we also studied the outer bulge cluster NGC 6652. We discuss the general properties of the proper-motion-cleaned Colour-Magnitude Diagrams, derived for the first time for most of them. A general conclusion is that the inner bulge globular clusters have clearly lower transverse motions (and spatial velocities) than halo clusters, and appear to be trapped in the bulge bar.

We propose an evolutionary model to describe the dynamical evolution of star cluster systems in tidal fields, in which we calibrated the parametric equations defining the model by running direct N-body simulations of star clusters with a wide range of initial masses and set of orbital parameters, living within the external tidal field generated by a disc-like galaxy. We derived a new method to solve numerically the evolutionary equations, allowing us to infer constraints on the mass of a star cluster from its age, present-day mass, orbital parameters and external gravitational potential. The result has been applied to the metal-rich subsample of Galactic globular clusters, being a good representation of a disc-bulge population. We reconstructed the initial mass function of these objects from the present-day mass function, finding that a lognormal distribution is well preserved during the evolution of the globular cluster system. The evolution of a power-law initial mass function has been evaluated, confirming that it transforms into a lognormal distribution of the cluster masses within an Hubble time. Our results are consistent with a formation scenario in which metal-rich Galactic globular clusters formed from giant molecular clouds in high-pressure regions during the early phases of the evolution of the Galactic disc and bulge.

The Sun Watcher with Active Pixels and Image Processing (SWAP) is an EUV solar telescope on board ESA's Project for Onboard Autonomy 2 (PROBA2) mission launched on 2 November 2009. SWAP has a spectral bandpass centered on 17.4 nm and provides images of the low solar corona over a 54x54 arcmin field-of-view with 3.2 arcsec pixels and an imaging cadence of about two minutes. SWAP is designed to monitor all space-weather-relevant events and features in the low solar corona. Given the limited resources of the PROBA2 microsatellite, the SWAP telescope is designed with various innovative technologies, including an off-axis optical design and a CMOS-APS detector. This article provides reference documentation for users of the SWAP image data.