Asteroid mass estimation with the robust adaptive Metropolis algorithm

Context. The bulk density of an asteroid informs us about its interior structure and composition. To constrain the bulk density, one needs an estimated mass of the asteroid. The mass is estimated by analyzing an asteroid's gravitational interaction with another object, such as another asteroid during a close encounter. An estimate for the mass has typically been obtained with linearized least-squares methods, despite the fact that this family of methods is not able to properly describe non-Gaussian parameter distributions. In addition, the uncertainties reported for asteroid masses in the literature are sometimes inconsistent with each other and are suspected to be unrealistically low.Aims. We aim to present a Markov-chain Monte Carlo (MCMC) algorithm for the asteroid mass estimation problem based on asteroid-asteroid close encounters. We verify that our algorithm works correctly by applying it to synthetic data sets. We use astrometry available through the Minor Planet Center to estimate masses for a select few example cases and compare our results with results reported in the literature.Methods. Our mass-estimation method is based on the robust adaptive Metropolis algorithm that has been implemented into the OpenOrb asteroid orbit computation software. Our method has the built-in capability to analyze multiple perturbing asteroids and test asteroids simultaneously.Results. We find that our mass estimates for the synthetic data sets are fully consistent with the ground truth. The nominal masses for real example cases typically agree with the literature but tend to have greater uncertainties than what is reported in recent literature. Possible reasons for this include different astrometric data sets and weights, different test asteroids, different force models or different algorithms. For (16) Psyche, the target of NASA's Psyche mission, our maximum likelihood mass is approximately 55% of what is reported in the literature. Such a low mass would imply that the bulk density is significantly lower than previously expected and thus disagrees with the theory of (16) Psyche being the metallic core of a protoplanet. We do, however, note that masses reported in recent literature remain within our 3-sigma limits.Results. The new MCMC mass-estimation algorithm performs as expected, but a rigorous comparison with results from a least-squares algorithm with the exact same data set remains to be done. The matters of uncertainties in comparison with other algorithms and correlations of observations also warrant further investigation.Peer reviewe

Masses, bulk densities, and macroporosities of asteroids (15) Eunomia, (29) Amphitrite, (52) Europa, and (445) Edna based on Gaia astrometry

Context. Gaia Data Release 2 (DR2) includes milliarcsecond-accuracy astrometry for 14 099 asteroids. One of the main expected scientific applications of these data is asteroid mass estimation via the modeling of perturbations during asteroid-asteroid encounters. Aims. We explore the practical impact of the Gaia astrometry of asteroids for the purpose of asteroid mass and orbit estimation by estimating the masses individually for four large asteroids. We use various combinations of Gaia astrometry and/or Earth-based astrometry so as to determine the impact of Gaia on the estimated masses. By utilizing published information about estimated volumes and meteorite analogs, we also derive estimates for bulk densities and macroporosities. Methods. We apply a Markov chain Monte Carlo (MCMC) algorithm for asteroid mass and orbit estimation by modeling asteroidasteroid close encounters to four separate large asteroids in an attempt to estimate their masses based on multiple simultaneously studied close encounters with multiple test asteroids. In order to validate our algorithm and data treatment, we apply the MCMC algorithm to pure orbit determination for the main-belt asteroid (367) Amicitia and compare the residuals to previously published ones. In addition, we attempt to estimate a mass for (445) Edna with Gaia astrometry alone based on its close encounter with (1764) Cogshall. Results. In the case of the orbit of (367) Amicitia, we find a solution that improves on the previously published solution. The study of (445) Edna reveals that mass estimation with DR2 astrometry alone is unfeasible and that it must be combined with astrometry from other sources to achieve meaningful results. We find that a combination of DR2 and Earth-based astrometry results in dramatically reduced uncertainties and, by extension, significantly improved results in comparison to those computed based on less accurate Earth-based astrometry alone. Conclusions. Our mass estimation algorithm works well with a combination of Gaia DR2 and Earth-based astrometry and provides very impressive results with significantly reduced uncertainties. We note that the DR2 has a caveat in that many asteroids suitable for mass-estimation purposes are not included in the relatively small sample. This limits the number of asteroids to which mass estimation can be applied. However, this issue will largely be corrected with the forthcoming third Gaia data release, which is expected to allow for a wave of numerous accurate mass estimates for a wide range of asteroids.Peer reviewe

Observational Hertzsprung-Russell diagrams

Context. Gaia Data Release 2 provides high-precision astrometry and three-band photometry for about 1.3 billion sources over the full sky. The precision, accuracy, and homogeneity of both astrometry and photometry are unprecedented. Aims. We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. Methods. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. Results. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Conclusions. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.Peer reviewe

Gaia Early Data Release 3 : The celestial reference frame (Gaia-CRF3)

Publisher Copyright: ©Context. Gaia-CRF3 is the celestial reference frame for positions and proper motions in the third release of data from the Gaia mission, Gaia DR3 (and for the early third release, Gaia EDR3, which contains identical astrometric results). The reference frame is defined by the positions and proper motions at epoch 2016.0 for a specific set of extragalactic sources in the (E)DR3 catalogue. Aims. We describe the construction of Gaia-CRF3 and its properties in terms of the distributions in magnitude, colour, and astrometric quality. Methods. Compact extragalactic sources in Gaia DR3 were identified by positional cross-matching with 17 external catalogues of quasi-stellar objects (QSO) and active galactic nuclei (AGN), followed by astrometric filtering designed to remove stellar contaminants. Selecting a clean sample was favoured over including a higher number of extragalactic sources. For the final sample, the random and systematic errors in the proper motions are analysed, as well as the radio-optical offsets in position for sources in the third realisation of the International Celestial Reference Frame (ICRF3). Results. Gaia-CRF3 comprises about 1.6 million QSO-like sources, of which 1.2 million have five-parameter astrometric solutions in Gaia DR3 and 0.4 million have six-parameter solutions. The sources span the magnitude range G = 13-21 with a peak density at 20.6 mag, at which the typical positional uncertainty is about 1 mas. The proper motions show systematic errors on the level of 12 μas yr-1 on angular scales greater than 15 deg. For the 3142 optical counterparts of ICRF3 sources in the S/X frequency bands, the median offset from the radio positions is about 0.5 mas, but it exceeds 4 mas in either coordinate for 127 sources. We outline the future of Gaia-CRF in the next Gaia data releases. Appendices give further details on the external catalogues used, how to extract information about the Gaia-CRF3 sources, potential (Galactic) confusion sources, and the estimation of the spin and orientation of an astrometric solution.Peer reviewe

Gaia Data Release 2 Variable stars in the colour-absolute magnitude diagram

Context. The ESA Gaia mission provides a unique time-domain survey for more than 1.6 billion sources with G less than or similar to 21 mag. Aims. We showcase stellar variability in the Galactic colour-absolute magnitude diagram (CaMD). We focus on pulsating, eruptive, and cataclysmic variables, as well as on stars that exhibit variability that is due to rotation and eclipses. Methods. We describe the locations of variable star classes, variable object fractions, and typical variability amplitudes throughout the CaMD and show how variability-related changes in colour and brightness induce "motions". To do this, we use 22 months of calibrated photometric, spectro-photometric, and astrometric Gaia data of stars with a significant parallax. To ensure that a large variety of variable star classes populate the CaMD, we crossmatched Gaia sources with known variable stars. We also used the statistics and variability detection modules of the Gaia variability pipeline. Corrections for interstellar extinction are not implemented in this article. Results. Gaia enables the first investigation of Galactic variable star populations in the CaMD on a similar, if not larger, scale as was previously done in the Magellanic Clouds. Although the observed colours are not corrected for reddening, distinct regions are visible in which variable stars occur. We determine variable star fractions to within the current detection thresholds of Gaia. Finally, we report the most complete description of variability-induced motion within the CaMD to date. Conclusions. Gaia enables novel insights into variability phenomena for an unprecedented number of stars, which will benefit the understanding of stellar astrophysics. The CaMD of Galactic variable stars provides crucial information on physical origins of variability in a way that has previously only been accessible for Galactic star clusters or external galaxies. Future Gaia data releases will enable significant improvements over this preview by providing longer time series, more accurate astrometry, and additional data types (time series BP and RP spectra, RVS spectra, and radial velocities), all for much larger samples of stars.Peer reviewe

Detection of a gamma-ray flare from the high-redshift blazar DA 193

High-redshift ($z>2$) blazars are the most powerful members of the blazar family. Yet, only a handful of them have both X-ray and $\gamma$-ray detection, thereby making it difficult to characterize the energetics of the most luminous jets. Here, we report, for the first time, the Fermi-Large Area Telescope detection of the significant $\gamma$-ray emission from the high-redshift blazar DA 193 ($z=2.363$). Its time-averaged $\gamma$-ray spectrum is soft ($\gamma$-ray photon index = $2.9\pm0.1$) and together with a relatively flat hard X-ray spectrum (14$-$195 keV photon index = $1.5\pm0.4$), DA 193 presents a case to study a typical high-redshift blazar with inverse Compton peak being located at MeV energies. An intense GeV flare was observed from this object in the first week of 2018 January, a phenomenon rarely observed from high-redshift sources. What makes this event a rare one is the observation of an extremely hard $\gamma$-ray spectrum (photon index = $1.7\pm0.2$), which is somewhat unexpected since high-redshift blazars typically exhibit a steep falling spectrum at GeV energies. The results of our multi-frequency campaign, including both space- (Fermi, NuSTAR, and Swift) and ground-based (Steward and Nordic Optical Telescope) observatories, are presented and this peculiar $\gamma$-ray flare is studied within the framework of a single-zone leptonic emission scenario.Comment: 14 pages, 9 figures, 2 tables, to appear in the Astrophysical Journa

Time-series analysis of long-term photometry of BM Canum Venaticorum

Long-term photometry is commonly used to monitor chromospheric activity of late-type stars. We study standard Johnson differential V photometry of the RS CVn binary BM Canum Venaticorum (BM CVn) spanning over a quarter of a century. Our main aims are to determine the activity cycles, the rate of surface differential rotation, and the rotation period of the active longitudes of BM CVn. The continuous period search (CPS) algorithm is applied to the photometry. The changes of the mean and amplitude of the light curves are used to search for activity cycles. The rotation period changes give an estimate of the rate of surface differential rotation. The Kuiper method is applied to the epochs of the primary and secondary minima to search for active longitudes. The photometry reveals the presence of a stable mean light curve (MLC) connected to the orbital period P-orb=20.(d)6252 of this binary. We remove this MLC from the original V magnitudes, which gives us the corrected V magnitudes. These two samples of Vand Vdata are analyzed separately with CPS. The fraction of unreliable CPS models decreases when the MLC is removed. The same significant activity cycle of approximately 12.5 years is detected in both V and V samples. The estimate for the surface differential rotation coefficient, k >= 0.10, is the same for both samples, but the number of unrealistic period estimates decreases after removing the MLC. The same active longitude period of P-al=20.(d)511 +/- 0.(d)005 is detected in the V and V magnitudes. This long-term regularity in the epochs of primary and secondary minima of the light curves is not caused by the MLC. On the contrary, the MLC hampers the detection of active longitudes.Peer reviewe

Gaia Focused Product Release: A catalogue of sources around quasars to search for strongly lensed quasars

Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of $0.18$" in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those expected for most lenses. Aims. We present the Data Processing and Analysis Consortium GravLens pipeline, which was built to analyse all Gaia detections around quasars and to cluster them into sources, thus producing a catalogue of secondary sources around each quasar. We analysed the resulting catalogue to produce scores that indicate source configurations that are compatible with strongly lensed quasars. Methods. GravLens uses the DBSCAN unsupervised clustering algorithm to detect sources around quasars. The resulting catalogue of multiplets is then analysed with several methods to identify potential gravitational lenses. We developed and applied an outlier scoring method, a comparison between the average BP and RP spectra of the components, and we also used an extremely randomised tree algorithm. These methods produce scores to identify the most probable configurations and to establish a list of lens candidates. Results. We analysed the environment of 3 760 032 quasars. A total of 4 760 920 sources, including the quasars, were found within 6" of the quasar positions. This list is given in the Gaia archive. In 87\% of cases, the quasar remains a single source, and in 501 385 cases neighbouring sources were detected. We propose a list of 381 lensed candidates, of which we identified 49 as the most promising. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the unique Gaia data for strong lensing studies further.Comment: 35 pages, 60 figures, accepted for publication by Astronomy and Astrophysic

Gaia Focused Product Release: Radial velocity time series of long-period variables

The third Gaia Data Release (DR3) provided photometric time series of more than 2 million long-period variable (LPV) candidates. Anticipating the publication of full radial-velocity (RV) in DR4, this Focused Product Release (FPR) provides RV time series for a selection of LPVs with high-quality observations. We describe the production and content of the Gaia catalog of LPV RV time series, and the methods used to compute variability parameters published in the Gaia FPR. Starting from the DR3 LPVs catalog, we applied filters to construct a sample of sources with high-quality RV measurements. We modeled their RV and photometric time series to derive their periods and amplitudes, and further refined the sample by requiring compatibility between the RV period and at least one of the $G$, $G_{\rm BP}$, or $G_{\rm RP}$ photometric periods. The catalog includes RV time series and variability parameters for 9\,614 sources in the magnitude range $6\lesssim G/{\rm mag}\lesssim 14$, including a flagged top-quality subsample of 6\,093 stars whose RV periods are fully compatible with the values derived from the $G$, $G_{\rm BP}$, and $G_{\rm RP}$ photometric time series. The RV time series contain a mean of 24 measurements per source taken unevenly over a duration of about three years. We identify the great most sources (88%) as genuine LPVs, with about half of them showing a pulsation period and the other half displaying a long secondary period. The remaining 12% consists of candidate ellipsoidal binaries. Quality checks against RVs available in the literature show excellent agreement. We provide illustrative examples and cautionary remarks. The publication of RV time series for almost 10\,000 LPVs constitutes, by far, the largest such database available to date in the literature. The availability of simultaneous photometric measurements gives a unique added value to the Gaia catalog (abridged)Comment: 36 pages, 38 figure

Gaia Early Data Release 3 Acceleration of the Solar System from Gaia astrometry

Context. Gaia Early Data Release 3 (Gaia EDR3) provides accurate astrometry for about 1.6 million compact (QSO-like) extragalactic sources, 1.2 million of which have the best-quality five-parameter astrometric solutions. Aims. The proper motions of QSO-like sources are used to reveal a systematic pattern due to the acceleration of the solar systembarycentre with respect to the rest frame of the Universe. Apart from being an important scientific result by itself, the acceleration measured in this way is a good quality indicator of the Gaia astrometric solution. Methods. Theeffect of the acceleration was obtained as a part of the general expansion of the vector field of proper motions in vector spherical harmonics (VSH). Various versions of the VSH fit and various subsets of the sources were tried and compared to get the most consistent result and a realistic estimate of its uncertainty. Additional tests with the Gaia astrometric solution were used to get a better idea of the possible systematic errors in the estimate. Results. Our best estimate of the acceleration based on Gaia EDR3 is (2.32 +/- 0.16) x 10(-10) m s(-2) (or 7.33 +/- 0.51 km s(-1) Myr-1) towards alpha = 269.1 degrees +/- 5.4 degrees, delta = -31.6 degrees +/- 4.1 degrees, corresponding to a proper motion amplitude of 5.05 +/- 0.35 mu as yr(-1). This is in good agreement with the acceleration expected from current models of the Galactic gravitational potential. We expect that future Gaia data releases will provide estimates of the acceleration with uncertainties substantially below 0.1 mu as yr(-1).Peer reviewe