Search

We present a novel method for automatically detecting and characterising semi-resolved star clusters: clusters where the observational point-spread function (PSF) is smaller than the cluster's radius, but larger than the separations between individual stars. We apply our method to a 1.77 deg$^2$ field located in the Large Magellanic Cloud (LMC) using the VISTA survey of the Magellanic Clouds (VMC), which surveyed the LMC in the $YJK_\text{s}$ bands. Our approach first models the position-dependent PSF to detect and remove point sources from deep $K_\text{s}$ images; this leaves behind extended objects such as star clusters and background galaxies. We then analyse the isophotes of these extended objects to characterise their properties, perform integrated photometry, and finally remove any spurious objects this procedure identifies. We demonstrate our approach in practice on a deep VMC $K_\text{s}$ tile that contains the most active star-forming regions in the LMC: 30 Doradus, N158, N159, and N160. We select this tile because it is the most challenging for automated techniques due both to crowding and nebular emission. We detect 682 candidate star clusters, with an estimated contamination rate of 13% from background galaxies and chance blends of physically unrelated stars. We compare our candidates to publicly available James Webb Space Telescope data and find that at least 80% of our detections appear to be star clusters., Comment: 26 pages, 17 figures, submitted to MNRAS

Analysing stellar parameters and abundances from nearly one million Gaia DR3 Radial Velocity Spectrometer (RVS) spectra poses challenges due to the limited spectral coverage (restricted to the infrared Ca II triplet) and variable signal-to-noise ratios of the data. To address this, we use The Cannon, a data-driven method, to transfer stellar parameters and abundances from the GALAH Data Release 4 (DR4; R ~ 28,000) catalogue to the lower resolution Gaia DR3 RVS spectra (R ~ 11,500). Our model, trained on 14,484 common targets, predicts parameters such as Teff, log g, and [Fe/H], along with several other elements across approximately 800,000 Gaia RVS spectra. We utilise stars from open and globular clusters present in the Gaia RVS catalogue to validate our predicted mean [Fe/H] with high precision (~0.02-0.10 dex). Additionally, we recover the bimodal distribution of [Ti/Fe] versus [Fe/H], reflecting the high and low alpha-components of Milky Way disk stars, demonstrating The Cannon's capability for accurate stellar abundance determination from medium-resolution Gaia RVS spectra. The methodologies and resultant catalogue presented in this work highlight the remarkable potential of the RVS dataset, which by the end of the Gaia mission will comprise spectra of over 200 million stars., Comment: Submitted to MNRAS, 16 pages, 15 figures

During the first few centuries CE, the centre of the known world gradually shifted from Alexandria to Constantinople. Combined with a societal shift from pagan beliefs to Christian doctrines, Antiquity gave way to the Byzantine era. While Western Europe entered an extended period of intellectual decline, Constantinople developed into a rich cultural crossroads between East and West. Yet, Byzantine scholarship in astronomy and geography continued to rely heavily on their ancient Greek heritage, and particularly on Ptolemy's Geography. Unfortunately, Ptolemy's choices for his geographic coordinate system resulted in inherent and significant distortions of and inaccuracies in maps centred on the Byzantine Empire. This comprehensive review of Byzantine geographic achievements -- supported by a review of astronomical developments pertaining to position determination on Earth -- aims to demonstrate why and how, when Constantinople fell to the Turks in 1453 and the Ottoman Empire commenced, Byzantine astronomers had become the central axis in an extensive network of Christians, Muslims and Jews. Their influence remained significant well into the Ottoman era, particularly in the context of geographical applications., Comment: Accepted for publication in De Medio Aevo, special issue: Costache D. (ed.), Scientific Interests and Technological Innovation in Byzantium: Interdisciplinary Perspectives; 20 pages, incl. 9 figures

Stellar mergers and accretion events have been crucial in shaping the evolution of the Milky Way (MW). These events have been dynamically identified and chemically characterised using red giants and main-sequence stars. RR Lyrae (RRL) variables can play a crucial role in tracing the early formation of the MW since they are ubiquitous, old (t$\ge$10 Gyr) low-mass stars and accurate distance indicators. We exploited Data Release 3 of the GALAH survey to identify 78 field RRLs suitable for chemical analysis. Using synthetic spectra calculations, we determined atmospheric parameters and abundances of Fe, Mg, Ca, Y, and Ba. Most of our stars exhibit halo-like chemical compositions, with an iron peak around [Fe/H]$\approx -$1.40, and enhanced Ca and Mg content. Notably, we discovered a metal-rich tail, with [Fe/H] values ranging from $-$1 to approximately solar metallicity. This sub-group includes almost ~1/4 of the sample, it is characterised by thin disc kinematics and displays sub-solar $\alpha$-element abundances, marginally consistent with the majority of the MW stars. Surprisingly, they differ distinctly from typical MW disc stars in terms of the s-process elements Y and Ba. We took advantage of similar data available in the literature and built a total sample of 535 field RRLs for which we estimated kinematical and dynamical properties. We found that metal-rich RRLs (1/3 of the sample) likely represent an old component of the MW thin disc. We also detected RRLs with retrograde orbits and provided preliminary associations with the Gaia-Sausage-Enceladus, Helmi, Sequoia, Sagittarius, and Thamnos stellar streams., Comment: Accepted for publication in MNRAS. 29 pages, 20 figures

Quasi-stellar objects (QSOs) are a basis for an absolute reference system for astrometric studies. There is a need for creating such system behind nearby galaxies, to facilitate the measuring of the proper motions of these galaxies. However, the foreground contamination from the galaxies themselves is a problem for the QSO identification. We search for new QSOs behind both Magellanic Clouds, the Magellanic Bridge, and the Magellanic Stream. We identify QSO candidates with a combination of near-infrared colors and variability criteria from the public ESO Visual and Infrared Survey Telescope for Astronomy (VISTA) Magellanic Clouds (VMC) survey. We confirm their nature from broad emission lines with low-resolution optical spectroscopy. We confirmed the QSO nature of 136 objects. They are distributed as follows: 12 behind the LMC, 37 behind the SMC, 63 behind the Bridge, and 24 behind the Stream. The QSOs span a redshift range from z~0.1 to z~2.9. A comparison of our quasar selection with the Quaia quasar catalog, based on Gaia low-resolution spectra, yields a selection and confirmation success rate of 6-19%, depending on whether the quality of the photometry, the magnitude ranges and the colors are considered. Our candidate list is rather incomplete, but the objects in it are likely to be confirmed as quasars with ~90% probability. Finally, we report a list of 3609 objects across the entire VMC survey that match our color and variability selection criteria; only 1249 of them have Gaia counterparts. Our combined infrared color and variability criteria for QSO selection prove to be efficient - ~90% of the observed candidates are bona fide QSOs and allow to generate a list of new high-probability quasar candidates., Comment: 20 pages, 9 figures, 5 tables

The presence of young stellar populations in the Large Magellanic Cloud cluster NGC 1783 has caught significant attention, with suggestions ranging from it being a genuine secondary stellar generation to a population of blue straggler stars or simply contamination from background stars. Thanks to multi-epoch observations with the Hubble Space Telescope, proper motions for stars within the field of NGC 1783 have been derived, thus allowing accurate cluster membership determination. Here, we report that the younger stars within NGC 1783 indeed belong to the cluster, and their spatial distribution is more extended compared to the bulk of the older stellar population, consistent with previous studies. Through N-body simulations, we demonstrate that the observed characteristics of the younger stars cannot be explained solely by blue straggler stars in the context of the isolated dynamical evolution of NGC 1783. Instead, accretion of the external, low-mass stellar system can better account for both the inverse spatial concentration and the radial velocity isotropy of the younger stars. We propose that NGC 1783 may have accreted external stars from low-mass stellar systems, resulting in a mixture of external younger stars and blue straggler stars from the older bulk population, thereby accounting for the characteristics of the younger sequence., Comment: 11 pages, 6 figures; accepted for publication in ApJ

The ability to measure precise and accurate stellar effective temperatures ($T_{\rm{eff}}$) and surface gravities ($\log(g)$) is essential in determining accurate and precise abundances of chemical elements in stars. Measuring $\log(g)$ from isochrones fitted to colour-magnitude diagrams of open clusters is significantly more accurate and precise compared to spectroscopic $\log(g)$. By determining the ranges of ages, metallicity, and extinction of isochrones that fit the colour-magnitude diagram, we constructed a joint probability distribution of $T_{\rm{eff}}$ and $\log(g)$. The joint photometric probability shows the complex correlations between $T_{\rm{eff}}$ and $\log(g)$, which depend on the evolutionary stage of the star. We show that by using this photometric prior while fitting spectra, we can acquire more precise spectroscopic stellar parameters and abundances of chemical elements. This reveals higher-order abundance trends in open clusters like traces of atomic diffusion. We used photometry and astrometry provided by the \textit{Gaia} DR3 catalogue, Padova isochrones, and Galactic Archaeology with HERMES (GALAH) DR4 spectra. We analysed the spectra of 1979 stars in nine open clusters, using MCMC to fit the spectroscopic abundances of 26 elements, $T_{\rm{eff}}$, $\log(g)$, $v_{\rm{mic}}$, and $v_{\rm{broad}}$. We found that using photometric priors improves the accuracy of abundances and $\log(g)$, which enables us to view higher-order trends of abundances caused by atomic diffusion in M67 and Ruprecht 147., Comment: 47 Pages. 14 figures in the main body. Accepted by MNRAS on the 8th of February 2024

Type II Cepheids (T2C) are less frequently used counterparts of classical Cepheids which provide the primary calibration of the distance ladder for measuring $H_0$ in the local Universe. In the era of the Hubble Tension, T2C variables with the RR Lyrae stars (RRL) and the tip of the red giant branch (TRGB) can potentially provide classical Cepheid independent calibration of the cosmic distance ladder. Our goal is to provide an absolute calibration of the Period-Luminosity, Period-Luminosity-Color and Period-Wesenheit relations(PL,PLC and PW) of T2Cs in the Large Magellanic Cloud (LMC). We exploited time-series photometry in the near-infrared (NIR) bands for a sample of more than 320 T2Cs in the Magellanic Clouds (MC). These observations were acquired during 2009-2018 in the context of the VMC ESO public survey (The VISTA near-infrared YJKs survey of the Magellanic System). The NIR photometry was supplemented with well-sampled optical light curves and accurate pulsation periods from the OGLE IV survey and the Gaia mission. We used the best-quality NIR light curves to generate custom templates for modelling sparsely sampled light curves in YJKs bands; in turn, we derived accurate and precise intensity-averaged mean magnitudes and pulsation amplitudes of 339 T2Cs in the MC. We used optical and NIR mean magnitudes to derive PL/PLC/PW relations for T2Cs in multiple bands, which were calibrated with the geometric distance to the LMC and with the Gaia parallaxes. We used our new empirical calibrations of PL/PW relations to obtain distances to 22 T2C-host Galactic globular clusters, which were found to be systematically smaller by 0.1 mag and 0.03-0.06 mag compared with the literature. A better agreement is found between our distances and those based on RRLs in globular clusters, providing strong support for using these population II stars with the TRGB for future distance scale studies.

The split main sequences found in the colour-magnitude diagrams of star clusters younger than ~600 Myr are suggested to be caused by the dichotomy of stellar rotation rates of upper main-sequence stars. Tidal interactions have been suggested as a possible explanation of the dichotomy of the stellar rotation rates. This hypothesis proposes that the slow rotation rates of stars along the split main sequences are caused by tidal interactions in binaries. To test this scenario, we measured the variations in the radial velocities of slowly rotating stars along the split main sequence of the young Galactic cluster NGC 2422 (~90 Myr) using spectra obtained at multiple epochs with the Canada-France-Hawai'i Telescope. Our results show that most slowly rotating stars are not radial-velocity variables. Using the theory of dynamical tides, we find that the binary separations necessary to fully or partially synchronise our spectroscopic targets, on time-scales shorter than the cluster age, predict much larger radial velocity variations across multiple-epoch observations, or a much larger radial velocity dispersion at a single epoch, than the observed values. This indicates that tidal interactions are not the dominant mechanism to form slowly rotating stars along the split main sequences. As the observations of the rotation velocity distribution among B- and A-type stars in binaries of larger separations hint at a much stronger effect of braking with age, we discuss the consequences of relaxing the constraints of the dynamical tides theory., Comment: 14 pages, 10 figures, 2 tables, accepted for publication in MNRAS

In our recent catalogue of BY Draconis (BY Dra) variables based on Zwicky Transient Facility data, we found traces of a period gap in the period-colour diagram. We combined our BY Dra database with catalogues from the {\sl Kepler} and K2 surveys, revealing a prominent period gap. Here, we use this combined ZTF-{\sl Kepler}-K2 data set to investigate the origin of the period gap observed for BY Dra stars using chromospheric activity indices. We use low- and medium-resolution spectra from the LAMOST Data Release 7 to derive magnetic activity indices for the Ca {\sc ii} H and K and H$\alpha$ emission lines. We find a strong dependence of chromospheric activity on both stellar mass and rotation period. For partially convective K-M-type stars, the activity decreases steeply up to an age of $\sim$700-1000 Myr, subsequently evolving to the type of low-level saturation associated with spin-down stallation. In contrast, F-G-type stars with thinner convective envelopes exhibit constant activity with increasing age. We suspect that the observed steep decrease for partially convective stars is driven by core-envelope coupling. This mechanism reduces differential rotation at the core-envelope transition, hence leading to decreased magnetic activity. Moreover, we derive activity indices for previously known star clusters and find similar trends as regards their activity levels as a function of age. In particular, very low-level activity is observed around the location of the period gap. Therefore, we conclude that the period gap, defined by the non-detection of variable sources, is driven by a minimum in chromospheric activity.

We present a kinematic analysis of the Small Magellanic Cloud using 3700 spectra extracted from the European Southern Observatory archive. We used data from Gaia and near-infrared photometry to select stellar populations and discard Galactic foreground stars. The sample includes main-sequence, red giant branch and red clump stars, observed with the Fibre Large Array Multi Wavelength Spectrograph. The spectra have a resolving power lambda/Delta(lambda) from 6500 to 38000. We derive radial velocities by employing a full spectrum fitting method using a penalised pixel fitting routine. We obtain a mean radial velocity for the galaxy of 159+/-2 km/s, with a velocity dispersion of 33+/-2 km/s. Our velocities agree with literature estimates for similar (young or old) stellar populations. The radial velocity of stars in the Wing and bar-like structure differ as a consequence of the dynamical interaction with the Large Magellanic Cloud. The higher radial velocity of young main-sequence stars in the bar compared to that of supergiants can be attributed to star formation around 40 Myr ago from gas already influenced by tidal stripping. Similarly, young main-sequence stars in the northern part of the bar, resulting from a prominent episode 25 Myr ago, have a higher radial velocity than stars in the southern part. Radial velocity differences between the northern and southern bar over densities are also traced by giant stars. They are corroborated by studies of the cold gas and proper motion indicating stretching/tidal stripping of the galaxy., Comment: 18 pages, 13 figures, accepted for publication in MNRAS

We discuss David Bohm's dual contributions as a physicist and thinker. First, de Grijs introduces Bohm's universe, with an emphasis on the physical quest that led Bohm to the elaboration of an original cosmology at the nexus of science and philosophy. Next, Costache takes his cue from de Grijs' explorations by highlighting the affinity between Bohm's scientific cosmology and patristic ideas that are central to the Orthodox worldview. It is our hope that this approach will stir the interest of Bohm scholars in the Orthodox worldview and also lead Orthodox theologians to nurture an appreciation for Bohm's cosmology., Comment: 19 pages; accepted for publication in "Theology and Science"

As the European maritime powers expanded their reach beyond north Atlantic coastal waters to distant lands as far away as the East Indies, access to a practical means of maritime navigation in the southern hemisphere became imperative. The first few voyages undertaken by the Dutch East India Company and its predecessor explicitly aimed at compiling star charts and constellations that were only visible south of the Equator, as practical navigation aids. The oldest known star atlas of southern constellations was published in 1603 by Frederick de Houtman. Controversies have plagued de Houtman's astronomical credentials from their inception, however, with contemporaries variously attributing the early southern star charts to Pieter Dirkszoon Keyser, de Houtman, or even to their tutor Petrus Plancius. The balance of available evidence suggests that Keyser initially led the astronomical observing campaign, ably assisted by de Houtman. Upon Keyser's untimely death, de Houtman embraced a leading role in compiling astronomical observations for maritime navigation purposes, whereas Plancius most probably led the delineation of the 12 new southern constellations that soon became part and parcel of the nautical consciousness., Comment: To be published in Essays on Astronomical History and Heritage: A Tribute to Wayne Orchiston on his 80th Birthday, eds Steven Gullberg and Peter Robertson, Springer (2023)

RR Lyrae variables are excellent population II distance indicators thanks to their well-defined period-luminosity relations (PLRs) at infrared wavelengths. We present results of near-infrared (NIR) monitoring of Galactic globular clusters to empirically quantify the metallicity dependence of NIR PLRs for RR Lyrae variables. Our sample includes homogeneous, accurate, and precise photometric data for 964 RR Lyrae variables in 11 globular clusters covering a large metallicity range ($\Delta\textrm{[Fe/H]}\sim2$~dex). We derive $JHK_s$ band period-luminosity-metallicity (PLZ) and period-Wesenheit-metallicity (PWZ) relations anchored using 346 Milky Way field RR Lyrae stars with {\it Gaia} parallaxes, and simultaneously solved for independent distances to globular clusters. We find a significant metallicity dependence of $\sim0.2$~mag/dex in $JHK_s$ band PLZ and PWZ relations for RR Lyrae stars independent of the adopted metallicity scale. The metallicity coefficients and the zero-points of the empirical PLZ and PWZ relations are in excellent agreement with the predictions from the horizontal branch evolution and pulsation models. Furthermore, RR Lyrae based distances to our sample of globular clusters are also statistically consistent with other independent measurements in the literature. Our recommended empirical $JHK_s$ band PLZ relations are also provided for RR Lyrae based distance measurements., Comment: 10 pages, 4 figures, abridged abstract, accepted in the ApJ Letters

Between the mid-eighteenth and mid-nineteenth centuries, long-haul oceanic voyages of exploration and discovery routinely carried astronomical tent observatories to support land-based longitude determinations using heavy and cumbersome astronomical regulators and transit telescopes. Following James Cook's deployment of a pilot tent observatory on his first voyage to the Pacific in 1768-1771, the tent design was altered by William Bayly for more convenient use on Cook's second and third voyages to the Pacific. Bayly's design became the standard structure of tent observatories assigned to shipboard astronomers during the Age of Sail. By the middle of the nineteenth century, a subtle shift in focus had occurred, with tent observatories now being deployed to observe specific celestial events (such as the 1882 Venus transit or a variety of eclipses), while longitude determinations increasingly relied on the novel, compact and improved box chronometers of the day. A further shift in the application of tent observatories occurred towards the end of the nineteenth century, when astronomical applications largely gave way to a renewed focus on meteorological measurements., Comment: 21 pages; Journal of Astronomical History and Heritage, in press (March 2023 issue)

By the twelfth century, northern European scholars gradually embraced Arabic innovations in science and technology. England naturally developed into a significant centre of the new learning in western Europe. Hereford, and specifically its cathedral school, played a particularly important role in the transition of English scholarship to the new learning. Hereford cathedral developed into a focal point for high-level scholarship, attracting numerous scholars from across the continent. Roger of Hereford stands out among his peers as an enlightened scholar who made more practical use than most of the full astronomical and astrological knowledge base available in England at the time. A significant body of recent scholarship focuses on twelfth-century ecclesiastical developments, including those relating to Roger of Hereford's Computus. However, much less scholarly emphasis is placed on Roger's astronomical calculations, particularly those which allowed him to establish an important reference meridian at Hereford. Those aspects are the focus of this paper., Comment: 17 pages, 5 figures; The Local Historian, in press (April 2023 issue; note that the version of record will have B/W figures)

Voyages of exploration often included astronomers among their crew to aid with maritime navigation. William Dawes, a British Marine who had been trained in practical astronomy, was assigned to the "First Fleet", a convoy of eleven ships that left England in May 1787 bound for Botany Bay (Sydney, Australia). Dawes was also expected to take measurements of the local gravitational acceleration, $g$, at any port of call by measuring the daily rate by which his Shelton pendulum clock differed from that at Greenwich, its calibration location. Although Dawes and Nevil Maskelyne, Britain's fifth Astronomer Royal, had planned to obtain clock-rate measurements in the Canary Islands, San Sebastian (Rio de Janeiro) and Table Bay, Captain Arthur Phillip, Commander of the First Fleet, only allowed Dawes to disembark the clock in Rio de Janeiro. Therefore, we have just one set of clock-rate measurements from the voyage, in addition to land-based measurements obtained in New South Wales. If gravity was the dominant factor affecting the clock's changing rate, Dawes' measurement of $-48.067$ sec per (sidereal) day obtained in Rio de Janeiro implies a local gravitational acceleration, $g = 9.7946$ m sec$^{-2}$. On the other hand, if we adopt the modern value, $g = 9.7878$ m s$^{-2}$, the implied daily decay rate is almost exactly 30 sec greater than Dawes' clock-rate determination, a difference that is well in excess of the prevailing uncertainties. This suggests that the pendulum's regulator nut may have been offset by a full turn, thus implying that our assumptions regarding the pendulum length may have to be revisited., Comment: 9 pages, 3 figures; Journal of Astronomical History and Heritage, in press (December 2022 issue)

In addition to the extended main-sequence turnoffs widely found in young and intermediate-age (~ 600 Myr-2 Gyr-old) star clusters, some younger clusters even exhibit split main sequences (MSs). Different stellar rotation rates are proposed to account for the bifurcated MS pattern, with red and blue MSs (rMS and bMS) populated by fast and slowly rotating stars, respectively. Using photometry from Gaia Early Data Release 3, we report a Galactic open cluster with a bifurcated MS, NGC 2422 ( ~ 90 Myr). We exclude the possibilities that the bifurcated MS pattern is caused by photometric noise or differential reddening. We aim to examine if stellar rotation can account for the split MSs. We use spectra observed with the Canada-France-Hawaii Telescope and the Southern African Large Telescope, and directly measured v sin i, the projected rotational velocities, for stars populating the bMS and rMS. We find that their v sin i values are weakly correlated with their loci in the color-magnitude diagram because of contamination caused by a large fraction of rMS stars with low projected rotational velocities. Based on the spectral energy distribution fitting method, we suggest that these slowly rotating stars at the rMS may hide a binary companion, which breaks the expected v sin i-color correlation. Future time-domain studies focusing on whether these slowly rotating stars are radial velocity variables are crucial to test the roles of stellar rotation and binarity in generating the split MSs., Comment: 20 pages, 12 figures, 1 table; Accepted for publication in The Astrophysical Journal (ApJ)

Multi-mode Cepheids pulsate simultaneously in more than one mode of oscillation. They provide an independent means to test stellar models and pulsation theories. They can also be used to derive metallicities. In recent years, the number of known multi-mode Cepheids has increased dramatically with the discovery of a large number of Galactic double-mode Cepheids. To date, 209 double-mode Cepheids have been detected in the Galactic bulge and disk, mostly based on the Optical Gravitational Lensing Experiment's (OGLE) catalog. In this paper, we conduct a comprehensive search for double-mode Cepheids in the northern sky based on Zwicky Transient Facility Data Release 5. We found 72 such objects in the Milky Way. The periods of the 30 sample objects already included in the OGLE catalog show excellent agreement with the OGLE periods. The period ratios of our new Cepheids are consistent with those of known double-mode Cepheids, as evidenced by their loci in the so-called `Petersen diagram'. Compared with OGLE, the completeness of our double-mode Cepheid sample is around 71\%. The much improved temporal sampling of the Zwicky Transient Facility offers significant scope to find more double-mode Cepheids, especially at the distribution's short-period end., Comment: 11 pages, 3 figures, 3 tables; accepted for publication in ApJS

The distribution of young stars into OB associations has long been in need of updating. High-precision {\sl Gaia} early Data Release 3 astrometry, coupled with modern machine-learning methods, allows this to be done. We have compiled a well-defined sample which includes OB stars and young open clusters, in total comprising about 47,700 objects. To break the sample down into groupings resembling associations, we applied the HDBSCAN$^{*}$ clustering algorithm. We used a Monte Carlo method to estimate the kinematic ages of the resulting clusters and the Student's $t$-test to assess the significance of the linear correlations between proper motions and coordinates, indicating the presence of possible cluster expansion signatures. The ages of the majority of clusters demonstrating a general expansion at a 1$\sigma$ confidence level are several tens of Myr, which is in agreement with the expected ages of OB associations. We found 32 open clusters which turned out to be members of the resulting groupings; their ages are consistent with one another within the uncertainties. Comparison of the clusters thus obtained with the historical composition of OB associations in the literature shows a correspondence between their positions in the Galaxy but an apparent absence of good one-to-one stellar matches. Therefore, we suggest that the historical composition of OB associations needs to be revised., Comment: 12 pages, 14 figures

We present an extensive catalogue of BY Draconis (BY Dra)-type variables and their stellar parameters. BY Dra are main-sequence FGKM-type stars. They exhibit inhomogeneous starspots and bright faculae in their photospheres. These features are caused by stellar magnetic fields, which are carried along with the stellar disc through rotation and which produce gradual modulations in their light curves (LCs). Our main objective is to characterise the properties of BY Dra variables over a wide range of stellar masses, temperatures and rotation periods. A recent study categorised 84,697 BY Dra variables from Data Release 2 of the Zwicky Transient Facility based on their LCs. We have collected additional photometric data from multiple surveys and performed broad-band spectral energy distribution fits to estimate stellar parameters. We found that more than half of our sample objects are of K spectral type, covering an extensive range of stellar parameters in the low-mass regime (0.1-1.3 M$_{\odot}$ ). Compared with previous studies, most of the sources in our catalogue are rapid rotators, and so most of them must be young stars for which a spin-down has not yet occurred. We subdivided our catalogue based on convection zone depth and found that the photospheric activity index, $S_{\rm ph}$, is lower for higher effective temperatures, i.e., for thinner convective envelopes. We observe a broad range of photospheric magnetic activity for different spectral classes owing to the presence of stellar populations of different ages. We found a higher magnetically active fraction for K- than M-type stars., Comment: 12 pages, 13 figures

We present a map of the total intrinsic reddening across ~90 deg$^{2}$ of the Large Magellanic Cloud (LMC) derived using optical (ugriz) and near-infrared (IR; YJKs) spectral energy distributions (SEDs) of background galaxies. The reddening map is created from a sample of 222,752 early-type galaxies based on the LEPHARE $\chi^{2}$ minimisation SED-fitting routine. We find excellent agreement between the regions of enhanced intrinsic reddening across the central (4x4 deg$^2$) region of the LMC and the morphology of the low-level pervasive dust emission as traced by far-IR emission. In addition, we are able to distinguish smaller, isolated enhancements that are coincident with known star-forming regions and the clustering of young stars observed in morphology maps. The level of reddening associated with the molecular ridge south of 30 Doradus is, however, smaller than in the literature reddening maps. The reduced number of galaxies detected in this region, due to high extinction and crowding, may bias our results towards lower reddening values. Our map is consistent with maps derived from red clump stars and from the analysis of the star formation history across the LMC. This study represents one of the first large-scale categorisations of extragalactic sources behind the LMC and as such we provide the LEPHARE outputs for our full sample of ~2.5 million sources., Comment: 18 pages, 16 figures, accepted for publication in MNRAS

We perform a statistical clustering analysis of upper main-sequence stars in the Large Magellanic Cloud (LMC) using data from the Visible and Infrared Survey Telescope for Astronomy survey of the Magellanic Clouds. We map over 2500 young stellar structures at 15 significance levels across ~120 square degrees centred on the LMC. The structures have sizes ranging from a few parsecs to over 1 kpc. We find that the young structures follow power-law size and mass distributions. From the perimeter-area relation, we derive a perimeter-area dimension of 1.44+-0.20. From the mass-size relation and the size distribution, we derive two-dimensional fractal dimensions of 1.50+-0.10 and 1.61+-0.20, respectively. We find that the surface density distribution is well-represented by a lognormal distribution. We apply the Larson relation to estimate the velocity dispersions and crossing times of these structures. Our results indicate that the fractal nature of the young stellar structures has been inherited from the gas clouds from which they form and that this architecture is generated by supersonic turbulence. Our results also suggest that star formation in the LMC is scale-free from 10 pc to 700 pc., Comment: 20 pages, 17 figures, accepted for publication in MNRAS

Multiple stellar populations (MPs) representing star-to-star light-element abundance variations are common in nearly all ancient Galactic globular clusters. Here we provide the strongest evidence yet that the populous, ~ 1.7 Gyr-old Large Magellanic Cloud cluster NGC 2173 also exhibits light-element abundance variations. Thus, our results suggest that NGC 2173 is the youngest cluster for which MPs have been confirmed to date. Our conclusion is based on the distinct bifurcation at the tip of its red-giant branch in high-quality color--magnitude diagrams generated from Hubble Space Telescope imaging observations. Our results are further supported by a detailed analysis of 'pseudo-$UBI$' maps, which reveal clear evidence of a bimodality in the cluster's red-giant-branch color distribution. Young clusters in the Magellanic Clouds can provide critical insights into galaxy evolution histories. Our discovery of MPs in NGC 2173 suggests that ancient Galactic globular clusters and young massive clusters might share a common formation process., Comment: 8 pages, 5 figures

The Large Magellanic Cloud (LMC) is the most luminous satellite galaxy of the Milky Way and owing to its companion, the Small Magellanic Cloud (SMC), represents an excellent laboratory to study the interaction of dwarf galaxies. The aim of this study is to investigate the kinematics of the outer regions of the LMC by using stellar proper motions to understand the impact of interactions, e.g. with the SMC about 250 Myr ago. {We calculate proper motions using multi-epoch $K_\mathrm{s}$-band images from the VISTA survey of the Magellanic Clouds system (VMC). Observations span a time baseline of 2$-$5 yr. We combine the VMC data with data from the Gaia early Data Release 3 and introduce a new method to distinguish between Magellanic and Milky Way stars based on a machine learning algorithm. This new technique enables a larger and cleaner sample selection of fainter sources as it reaches below the red clump of the LMC. We investigate the impact of the SMC on the rotational field of the LMC and find hints of stripped SMC debris. The south east region of the LMC shows a slow rotational speed compared to the overall rotation. $N$-body simulations suggest that this could be caused by a fraction of stripped SMC stars, located in that particular region, that move opposite to the expected rotation., Comment: 21 pages, 21 figures, accepted for publication in A&A

High-level Chinese cartographic developments predate European innovations by several centuries. Whereas European cartographic progress -- and in particular the search for a practical solution to the perennial "longitude problem" at sea -- was driven by persistent economic motivations, Chinese mapmaking efforts responded predominantly to administrative, cadastral and topographic needs. Nevertheless, contemporary Chinese scholars and navigators, to some extent aided by experienced Arab navigators and astronomers, developed independent means of longitude determination both on land and at sea, using a combination of astronomical observations and timekeeping devices that continued to operate adequately on pitching and rolling ships. Despite confusing and speculative accounts in the current literature and sometimes overt nationalistic rhetoric, Chinese technical capabilities applied to longitude determination at sea, while different in design from European advances owing to cultural and societal circumstances, were at least on a par with those of their European counterparts., Comment: 23 pages incl. 9 figures; Journal of Astronomical History and Heritage, in press (March 2022 issue)

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., Comment: Accepted for publication in MNRAS, 11 pages, 6 figures

Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing the galaxy's large-scale structure. We have developed a technique for producing high-resolution submillimeter images of galaxies of large angular size by using the telescope on the ground to determine the small-scale structure (the large Fourier components) and a space telescope (Herschel or Planck) to determine the large-scale structure (the small Fourier components). Using this technique, we are carrying out the HARP and SCUBA-2 High Resolution Terahertz Andromeda Galaxy Survey (HASHTAG), an international Large Program on the James Clerk Maxwell Telescope, with one aim being to produce the first high-fidelity high-resolution submillimeter images of Andromeda. In this paper, we describe the survey, the method we have developed for combining the space-based and ground-based data, and present the first HASHTAG images of Andromeda at 450 and 850um. We also have created a method to predict the CO(J=3-2) line flux across M31, which contaminates the 850um band. We find that while normally the contamination is below our sensitivity limit, the contamination can be significant (up to 28%) in a few of the brightest regions of the 10 kpc ring. We therefore also provide images with the predicted line emission removed., Comment: 26 pages, 19 figures. Submitted to ApJS June 2021, Accepted September 2021

We have derived high-spatial-resolution metallicity maps covering $\sim$105~deg$^2$ across the Large Magellanic Cloud (LMC) using near-infrared passbands from the VISTA Survey of the Magellanic Clouds. We attempt to understand the metallicity distribution and gradients of the LMC up to a radius of $\sim$ 6~kpc. We identify red giant branch (RGB) stars in spatially distinct $Y, (Y-K_{\rm s})$ colour-magnitude diagrams. In any of our selected subregions, the RGB slope is used as an indicator of the average metallicity, based on calibration to metallicity using spectroscopic data. The mean LMC metallicity is [Fe/H] = $-$0.42~dex ($\sigma$[Fe/H] = 0.04~dex). We find the bar to be mildly metal-rich compared with the outer disc, showing evidence of a shallow gradient in metallicity ($-0.008 \pm 0.001$ dex kpc$^{-1}$) from the galaxy's centre to a radius of 6~kpc. Our results suggest that the LMC's stellar bar is chemically similar to the bars found in large spiral galaxies. The LMC's radial metallicity gradient is asymmetric. It is metal-poor and flatter towards the southwest, in the direction of the Bridge. This hints at mixing and/or distortion of the spatial metallicity distribution, presumably caused by tidal interactions between the Magellanic Clouds., Comment: 13 pages, 15 figures, accepted for publication in the MNRAS

We derive the spatially-resolved star formation history (SFH) for a $96$ deg$^2$ area across the main body of the Large Magellanic Cloud (LMC), using the near-infrared photometry from the VISTA survey of the Magellanic Clouds (VMC). The data and analyses are characterised by a great degree of homogeneity and a low sensitivity to the interstellar extinction. 756 subregions of size $0.125$ deg$^2$ -- corresponding to projected sizes of about $296\times322\,\mathrm{pc}^{2}$ in the LMC -- are analysed. The resulting SFH maps, with typical resolution of $0.2$--$0.3$ dex in logarithm of age, reveal main features in the LMC disc at different ages: the patchy star formation at recent ages, the concentration of star formation on three spiral arms and on the Bar up to ages of $\sim\!1.6$ Gyr, and the wider and smoother distribution of older populations. The period of most intense star formation occurred roughly between 4 and 0.5 Gyr ago, at rates of $\sim\!0.3\,\mathrm{M}_{\odot}\mathrm{yr}^{-1}$. We compare young and old star formation rates with the observed numbers of RR Lyrae and Cepheids. We also derive a mean extinction and mean distance for every subregion, and the plane that best describes the spatial distribution of the mean distances. Our results cover an area about 50 per cent larger than the classical SFH maps derived from optical data by Harris & Zaritsky (2009). Main differences with respect to those maps are lower star formation rates at young ages, and a main peak of star formation being identified at ages slightly younger than $1$ Gyr., Comment: Accepted for publication in MNRAS

Angular momentum is a key property regulating star formation and evolution. However, the physics driving the distribution of the stellar rotation rates of early-type main-sequence stars is as yet poorly understood. Using our catalog of 40,034 early-type stars with homogeneous $v\sin i$ parameters, we review the statistical properties of their stellar rotation rates. We discuss the importance of possible contaminants, including binaries and chemically peculiar stars. Upon correction for projection effects and rectification of the error distribution, we derive the distributions of our sample's equatorial rotation velocities, which show a clear dependence on stellar mass. Stars with masses less than $2.5\ {M_\odot}$ exhibit a unimodal distribution, with the peak velocity ratio increasing as stellar mass increases. A bimodal rotation distribution, composed of two branches of slowly and rapidly rotating stars, emerges for more massive stars ($M>2.5\ {M_\odot}$). For stars more massive than $3.0\ {M_\odot}$, the gap between the bifurcated branches becomes prominent. For the first time, we find that metal-poor ([M/H] $< -0.2$ dex) stars only exhibit a single branch of slow rotators, while metal-rich ([M/H] $> 0.2$ dex) stars clearly show two branches. The difference could be attributed to unexpectedly high spin-down rates and/or in part strong magnetic fields in the metal-poor subsample., Comment: 18 pages, 9 figures. Accepted for publication in ApJ; revised references for section 3.2

We derive stellar parameters and abundances (`stellar labels') of 40,034 late-B and A-type main-sequence stars extracted from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope Medium Resolution Survey (LAMOST--MRS). The primary selection of our early-type sample was obtained from LAMOST Data Release 7 based on spectral line indices. We employed the Stellar LAbel Machine (SLAM) to derive their spectroscopic stellar parameters, drawing on Kurucz spectral synthesis models with 6000 K $< T_\mathrm{eff} <$ 15,000 K and $-1$ dex $< \mathrm{[M/H]} <$ 1 dex. For a signal-to-noise ratio of $\sim 60$, the cross-validated scatter is $\sim$75 K, 0.06 dex, 0.05 dex, and $\sim 3.5\,\mathrm{km\,s^{-1}}$ for $T_\mathrm{eff}$, $\log g$, [M/H], and $v\sin i$, respectively. A comparison with objects with prior, known stellar labels shows great consistency for all stellar parameters, except for $\log g$. Although this is an intrinsic caveat that comes from the MRS's narrow wavelength coverage, it only has a minor effect on estimates of the stellar rotation rates because of the decent spectral resolution and the profile-fitting method employed. The masses and ages of our early-type sample stars were inferred from non-rotating stellar evolution models. This paves the way for reviewing the properties of stellar rotation distributions as a function of stellar mass and age., Comment: 18 pages, 10 figures, 1 table. Accepted for publication in ApJS

Despite frequent references in modern reviews to a seventeenth-century Venetian longitude prize, only a single, circumstantial reference to the alleged prize is known from contemporary sources. Edward Harrison's scathing assessment of the conditions governing the award of an alleged Venetian longitude prize simultaneously disparages the rewards offered by the Dutch States General. However, the latter had long run its course by 1696, the year of the citation, thus rendering Harrison's reference unreliable. Whereas other longitude awards offered by the leading European maritime nations attracted applicants from far and wide, often accompanied by extensive, self-published pamphlets, the alleged Venetian prize does not seem to have been subject to similar hype. The alleged existence of seventeenth-century Venetian award is particularly curious, because the city's fortune was clearly in decline, and longitude determination on the open seas does not appear to have been a priority; the city's mariners already had access to excellent "portolan" charts. It is therefore recommended that authors refrain from referring to a potentially phantom Venetian longitude prize in the same context as the major sixteenth- to eighteenth-century European awards offered by the dominant sea-faring nations., Comment: 13 pages; Journal of Astronomical History and Heritage, in press (September 2021 issue)

Although governments across Europe had realised the need to incentivise the development of practically viable longitude solutions as early as the late-sixteenth century, the English government was late to the party. An sense of urgency among the scientific community and maritime navigators led to the establishment of a number of longitude awards by private donors. The first private British award was bequeathed in 1691 by Thomas Axe, parish clerk of Ottery St. Mary (Devon). Despite the absence of an expenses component and the onerous and costly nature of its terms and conditions, the Axe prize attracted a number of optimistic claimants. Although the award was never disbursed, it may have contributed to the instigation of the government-supported monetary reward associated with the British Longitude Act of 1714. It is likely that the conditions governing the British "Longitude Prize", specifically the required accuracy and the need for sea trials and of disclosure of a successful method's theoretical principles, can be traced back at least in part to the Axe Prize requirements., Comment: 11 pages; Journal of Astronomical History and Heritage, in press (September 2021 issue)

We study the morphology of the stellar periphery of the Magellanic Clouds in search of substructure using near-infrared imaging data from the VISTA Hemisphere Survey (VHS). Based on the selection of different stellar populations using the ($J-K_\mathrm{s}$, $K_\mathrm{s}$) colour-magnitude diagram, we confirm the presence of substructures related to the interaction history of the Clouds and find new substructures on the easter side of the LMC disc which may be owing to the influence of the Milky Way, and on the northern side of the SMC, which is probably associated to the ellipsoidal structure of the galaxy. We also study the luminosity function of red clump stars in the SMC and confirm the presence of a bi-modal distance distribution, in the form of a foreground population. We find that this bi-modality is still detectable in the eastern regions of the galaxy out to a 10 deg distance from its centre. Additionally, a background structure is detected in the North between 7 and 10 deg from the centre which might belong to the Counter Bridge, and a foreground structure is detected in the South between 6 and 8 deg from the centre which might be linked to the Old Bridge., Comment: 20 pages, 19 figures, accepted for publication in MNRAS

We independently determine the zero-point offset of the Gaia early Data Release-3 (EDR3) parallaxes based on $\sim 110,000$ W Ursae Majoris (EW)-type eclipsing binary systems. EWs cover almost the entire sky and are characterized by a relatively complete coverage in magnitude and color. They are an excellent proxy for Galactic main-sequence stars. We derive a $W1$-band Period-Luminosity relation with a distance accuracy of $7.4\%$, which we use to anchor the Gaia parallax zero-point. The final, global parallax offsets are $-28.6\pm0.6$ $\mu$as and $-25.4\pm4.0$ $\mu$as (before correction) and $4.2\pm0.5$ $\mu$as and $4.6\pm3.7$ $\mu$as (after correction) for the five- and six-parameter solutions, respectively. The total systematic uncertainty is $1.8$ $\mu$as. The spatial distribution of the parallax offsets shows that the bias in the corrected Gaia EDR3 parallaxes is less than 10 $\mu$as across $40\%$ of the sky. Only $15\%$ of the sky is characterized by a parallax offset greater than 30 $\mu$as. Thus, we have provided independent evidence that the parallax zero-point correction provided by the Gaia team significantly reduces the prevailing bias. Combined with literature data, we find that the overall Gaia EDR3 parallax offsets for Galactic stars are $[-20, -30]$ $\mu$as and 4-10 $\mu$as, respectively, before and after correction. For specific regions, an additional deviation of about 10 $\mu$as is found., Comment: 7 pages, 4 figures. Accepted for publication in The Astrophysical Journal Letters

The conventional picture of coeval, chemically homogeneous, populous star clusters -- known as `simple stellar populations' (SSPs) -- is a view of the past. Photometric and spectroscopic studies reveal that almost all ancient globular clusters in the Milky Way and our neighbouring galaxies exhibit star-to-star light-element abundance variations, typically known as 'multiple populations' (MPs). Here, we analyse photometric $\it Hubble$ $\it Space$ $\it Telescope$ observations of three young ($<$2 Gyr-old) Large and Small Magellanic Cloud clusters, NGC 411, NGC 1718 and NGC 2213. We measure the widths of their red-giant branches (RGBs). For NGC 411, we also use a pseudo-colour--magnitude diagram (pseudo-CMD) to assess its RGB for evidence of MPs. We compare the morphologies of the clusters' RGBs with artificially generated SSPs. We conclude that their RGBs do not show evidence of significant broadening beyond intrinsic photometric scatter, suggesting an absence of significant chemical abundance variations in our sample clusters. Specifically, for NGC 411, NGC 1718 and NGC 2213 we derive maximum helium-abundance variations of delta_Y=0.003$\pm$0.001 Y=0.300), 0.002$\pm$0.001 (Y=0.350) and 0.004$\pm$0.002 (Y=0.300), respectively. We determined an upper limit to the NGC 411 nitrogen-abundance variation of $\Delta$[N/Fe] = 0.3 dex; the available data for our other clusters do not allow us to determine useful upper limits. It thus appears that the transition from SSPs to MPs occurs at an age of ~2 Gyr, implying that age might play an important role in this transition. This raises the question as to whether this is indeed a fundamental minimum-age limit for the formation of MPs., Comment: 11 pages, 5 Figures

The color-magnitude diagrams (CMDs) of young star clusters show that, particularly at ultraviolet wavelengths, their upper main sequences (MSs) bifurcate into a sequence comprising the bulk population and a blue periphery. The spatial distribution of stars is crucial to understand the reasons for these distinct stellar populations. This study uses high-resolution photometric data obtained with the Hubble Space Telescope to study the spatial distributions of the stellar populations in seven Magellanic Cloud star clusters. The cumulative radial number fractions of blue stars within four clusters are strongly anti-correlated with those of the high-mass-ratio binaries in the bifurcated region, with negative Pearson coefficients < -0.7. Those clusters generally are young or in an early dynamical evolutionary stage. In addition, a supporting N-body simulation suggests the increasing percentage of blue-MS stars from the cluster centers to their outskirts may be associated with the dissolution of soft binaries. This study provides a different perspective to explore the MS bimodalities in young clusters and adds extra puzzles. A more comprehensive study combined with detailed simulations is needed in the future., Comment: 19 pages, 9 figures, 3 tables. Accepted for publication in ApJ

The impact of stellar rotation on the morphology of star cluster colour-magnitude diagrams is widely acknowledged. However, the physics driving the distribution of the equatorial rotation velocities of main-sequence turn-off (MSTO) stars is as yet poorly understood. Using Gaia Data Release 2 photometry and new Southern African Large Telescope medium-resolution spectroscopy, we analyse the intermediate-age ($\sim1\,$Gyr-old) Galactic open clusters NGC 3960, NGC 6134 and IC 4756 and develop a novel method to derive their stellar rotation distributions based on SYCLIST stellar rotation models. Combined with literature data for the open clusters NGC 5822 and NGC 2818, we find a tight correlation between the number ratio of slow rotators and the clusters' binary fractions. The blue-main-sequence stars in at least two of our clusters are more centrally concentrated than their red-main-sequence counterparts. The origin of the equatorial stellar rotation distribution and its evolution remains as yet unidentified. However, the observed correlation in our open cluster sample suggests a binary-driven formation mechanism., Comment: 9 pages, 7 figures, 1 table. Accepted for publication in MNRAS

We present a catalog of eclipsing binaries in the northern Galactic Plane from the Kiso Wide-Field Camera Intensive Survey of the Galactic Plane (KISOGP). We visually identified 7055 eclipsing binaries spread across $\sim$330 square degrees, including 4197 W Ursa Majoris/EW-, 1458 $\beta$ Lyrae/EB-, and 1400 Algol/EA-type eclipsing binaries. For all systems, $I$-band light curves were used to obtain accurate system parameters. We derived the distances and extinction values for the EW-type objects from their period--luminosity relation. We also obtained the structure of the thin disk from the distribution of our sample of eclipsing binary systems, combined with those of high-mass star-forming regions and Cepheid tracers. We found that the thin disk is inhomogeneous in number density as a function of Galactic longitude. Using this new set of distance tracers, we constrain the detailed structure of the thin disk. Finally, we report a global parallax zero-point offset of $ \Delta \pi=-42.1\pm1.9\mbox{(stat.)}\pm12.9\mbox{(syst.)}$ $\mu$as between our carefully calibrated EW-type eclipsing binary positions and those provided by Gaia Early Data Release 3. Implementation of the officially recommended parallax zero-point correction results in a significantly reduced offset. Additionally, we provide a photometric characterization of our EW-type eclipsing binaries that can be applied to further analyses., Comment: 17 pages, 9 figures, 2 tables. Accepted for publication in AJ

On 13 May 1787, a convict fleet of 11 ships left Portsmouth, England, on a 24,000 km, 8-month-long voyage to New South Wales. The voyage would take the "First Fleet" under Captain Arthur Phillip via Tenerife (Canary Islands), the port of Rio de Janeiro (Brazil), Table Bay at the southern extremity of the African continent and the southernmost cape of present-day Tasmania to their destination of Botany Bay. Given the navigation tools available at the time and the small size of the convoy's ships, their safe arrival within a few days of each other was a phenomenal achievement. This was particularly so, because they had not lost a single ship and only a relatively small number of crew and convicts. Phillip and his crew had only been able to ensure their success because of the presence of crew members who were highly proficient in practical astronomy, most notably Lieutenant William Dawes. We explore in detail his educational background and the events leading up to Dawes' appointment by the Board of Longitude as the convoy's dedicated astronomer-cum-Marine. In addition to Dawes, John Hunter, second captain of the convoy's flagship H.M.S. Sirius, Lieutenant William Bradley and Lieutenant Philip Gidley King were also experts in navigation and longitude determination, using both chronometers and "lunar distance" measurements. The historical record of the First Fleet's voyage is remarkably accurate, even by today's standards., Comment: 36 pages, including 19 figures; to appear in the Journal of Astronomical History and Heritage (March 2021 issue)