Search

Context. The Villafranca project is combining Gaia data with ground-based surveys to analyze Galactic stellar groups (clusters, associations, or parts thereof) with OB stars. Aims. We want to analyze the poorly studied cluster Stock 18 within the Villafranca project, as it is a very young stellar cluster with a symmetrical and compact H II region around it, Sh 2-170, so it is likely to provide insights into the structure and dynamics of such objects at an early stage of their evolution. Methods. We used Gaia astrometry, photometry, spectrophotometry, and variability data as well as ground-based spectroscopy and imaging to determine the characteristics of Stock 18. We used these data to analyze its core, massive star population, extinction, distance, membership, internal dynamics, density profile, IMF, stellar variability, and Galactic location. Results. Stock 18 is a very young (∼1.0 Ma) cluster located at a distance of 2.91 ± 0.10 kpc and is dominated by the GLS 13 370 system, whose primary (Aa) is an O9 V star. We propose that Stock 18 was in a very compact state (∼0.1 pc) about 1.0 Ma ago and that most massive stars were ejected at that time without significantly affecting the less massive stars as a result of multi-body dynamical interactions. Different age estimates also point toward an age close to 1.0 Ma, indicating that the dynamical interactions took place very shortly after massive star formation. Well-defined expanding stellar clusters have been observed before, but none are as young as this one. If we include all of the stars, the initial mass function is top heavy, but if we discard the ejected ones, it becomes nearly canonical. Therefore, this is another example (in addition to the previous one we found – the Bermuda cluster) of (a) a very young cluster with an already evolved present day mass function (b) that has significantly contributed to the future population of free-floating compact objects. If confirmed in more clusters, the number of such compact objects may be higher in the Milky Way than previously thought. Stock 18 has a variable extinction with an average value of R5495 higher than the canonical one of 3.1. We have discovered a new visual component (Ab) in the GLS 13 370 system. The cluster is above our Galactic mid-plane, likely as a result of the Galactic warp, and it has a distinct motion with respect to its surrounding old population, which is possibly an influence of the Perseus spiral arm. [ABSTRACT FROM AUTHOR]

We present a hierarchical Bayesian pipeline, BP3M, that measures positions, parallaxes, and proper motions (PMs) for cross-matched sources between Hubble Space Telescope (HST) images and Gaia—even for sparse fields (N * < 10 per image)—expanding from the recent GaiaHub tool. This technique uses Gaia-measured astrometry as priors to predict the locations of sources in HST images, and is therefore able to put the HST images onto a global reference frame without the use of background galaxies/QSOs. Testing our publicly available code in the Fornax and Draco dwarf spheroidal galaxies, we measure PMs that are a median of 8–13 times more precise than Gaia DR3 alone for 20.5 < G < 21 mag. We are able to explore the effect of observation strategies on BP3M astrometry using synthetic data, finding an optimal strategy to improve parallax and position precision at no cost to the PM uncertainty. Using 1619 HST images in the sparse COSMOS field (median nine Gaia sources per HST image), we measure BP3M PMs for 2640 unique sources in the 16 < G < 21.5 mag range, 25% of which have no Gaia PMs; the median BP3M PM uncertainty for 20.25 < G < 20.75 mag sources is 0.44 mas yr−1 compared to 1.03 mas yr−1 from Gaia, while the median BP3M PM uncertainty for sources without Gaia-measured PMs (20.75 < G < 21.5 mag) is 1.16 mas yr−1. The statistics that underpin the BP3M pipeline are a generalized way of combining position measurements from different images, epochs, and telescopes, which allows information to be shared between surveys and archives to achieve higher astrometric precision than that from each catalog alone. [ABSTRACT FROM AUTHOR]

We analyze four epochs of Hubble Space Telescope imaging over 18 yr for the Draco dwarf spheroidal galaxy. We measure precise proper motions for hundreds of stars and combine these with existing line-of-sight (LOS) velocities. This provides the first radially resolved 3D velocity dispersion profiles for any dwarf galaxy. These constrain the intrinsic velocity anisotropy and resolve the mass–anisotropy degeneracy. We solve the Jeans equations in oblate axisymmetric geometry to infer the mass profile. We find the velocity dispersion to be radially anisotropic along the symmetry axis and tangentially anisotropic in the equatorial plane, with a globally averaged value β B ¯ = − 0.20 − 0.53 + 0.28 , (where 1 – β B ≡ 〈 v tan 2 〉 / 〈 v rad 2 〉 in 3D). The logarithmic dark matter (DM) density slope over the observed radial range, Γdark, is − 0.83 − 0.37 + 0.32 , consistent with the inner cusp predicted in ΛCDM cosmology. As expected given Draco's low mass and ancient star formation history, it does not appear to have been dissolved by baryonic processes. We rule out cores larger than 487, 717, and 942 pc at 1 σ, 2 σ, and 3 σ confidence, respectively, thus imposing important constraints on the self-interacting DM cross section. Spherical models yield biased estimates for both the velocity anisotropy and the inferred slope. The circular velocity at our outermost data point (900 pc) is 24.19 − 2.97 + 6.31 km s − 1 . We infer a dynamical distance of 75.37 − 4.00 + 4.73 kpc and show that Draco has a modest LOS rotation, with v / σ = 0.22 ± 0.09 . Our results provide a new stringent test of the so-called "cusp–core" problem that can be readily extended to other dwarfs. [ABSTRACT FROM AUTHOR]

Aims. We use near-infrared, ground-based data from the VISTA Variables in the Via Lactea (VVV) survey to indirectly extend the astrometry provided by the Gaia catalog to objects in heavily extinct regions toward the Galactic bulge and plane that are beyond Gaia's reach. Methods. We made use of state-of-the-art techniques developed for high-precision astrometry and photometry with the Hubble Space Telescope to process the VVV data. We employed empirical, spatially variable, effective point spread functions and local transformations to mitigate the effects of systematic errors, like residual geometric distortion and image motion, and to improve measurements in crowded fields and for faint stars. We also anchored our astrometry to the absolute reference frame of Gaia Data Release 3. Results. We measure between 20 and 60 times more sources than Gaia in the region surrounding the Galactic center, obtaining a single-exposure precision of about 12 mas and a proper-motion precision of better than 1 mas yr−1 for bright, unsaturated sources. Our astrometry provides an extension of Gaia into the Galactic center. We publicly release the astro-photometric catalogs of the two VVV fields considered in this work, which contain a total of ~3.5 million sources. Our catalogs cover ~3 sq. deg, about 0.5% of the entire VVV survey area. [ABSTRACT FROM AUTHOR]

Core-collapse supernovae (CCSNe) are considered the primary magnetar formation channel, with 15 magnetars associated with supernova remnants (SNRs). A large fraction of these should occur in massive stellar binaries that are disrupted by the explosion, meaning that |$\sim 45~{{\ \rm per\ cent}}$| of magnetars should be nearby high-velocity stars. Here, we conduct a multiwavelength search for unbound stars, magnetar binaries, and SNR shells using public optical (uvgrizy bands), infrared (J, H, K , and Ks bands), and radio (888 MHz, 1.4 GHz, and 3 GHz) catalogues. We use Monte Carlo analyses of candidates to estimate the probability of association with a given magnetar based on their proximity, distance, proper motion, and magnitude. In addition to recovering a proposed magnetar binary, a proposed unbound binary, and 13 of 15 magnetar SNRs, we identify two new candidate unbound systems: an OB star from the Gaia catalogue we associate with SGR J1822.3−1606, and an X-ray pulsar we associate with 3XMM J185246.6 + 003317. Using a Markov Chain Monte Carlo simulation that assumes all magnetars descend from CCSNe, we constrain the fraction of magnetars with unbound companions to |$5\lesssim f_u \lesssim 24~{{\ \rm per\ cent}}$|⁠ , which disagrees with neutron star population synthesis results. Alternate formation channels are unlikely to wholly account for the lack of unbound binaries as this would require |$31\lesssim f_{nc} \lesssim 66~{{\ \rm per\ cent}}$| of magnetars to descend from such channels. Our results support a high fraction (⁠|$48\lesssim f_m \lesssim 86~{{\ \rm per\ cent}}$|⁠) of pre-CCSN mergers, which can amplify fossil magnetic fields to preferentially form magnetars. [ABSTRACT FROM AUTHOR]

Context. The nuclear stellar disk has been the most prolific star-forming region in the Milky Way over the past ∼30 million years. Notably, the cumulative mass of the three clusters currently found in the nuclear stellar disk, the Quintuplet, the Arches, and the Nuclear clusters, amounts to just 10% of the total anticipated mass of young stars that formed in this period. This discrepancy, known as the missing cluster problem, is attributed to factors such as high stellar density and tidal forces. Traces of dissolving clusters may exist as comoving groups of stars, providing insights into the star formation history of the region. Recently, a new cluster candidate associated with an HII region was reported through the analysis of kinematic data Aims. Our aim is to determine whether the young and massive stellar objects in the region share proper motion, positions in the plane of the sky, and line-of-sight distances. We use reddening as a proxy for the distances. Methods. We reduced and analyzed integral field spectroscopy data from the KMOS instrument at the ESO VLT to locate possible massive young stellar objects in the field. Then, we identified young massive stars with astrophotometric data from the two different catalogs to analyze their extinction and kinematics. Results. We present a group of young stellar objects that share velocities, are close together in the plane of the sky, and are located at a similar depth in the nuclear stellar disk. Conclusions. The results presented here offer valuable insights into the missing clusters problem. They indicate that not all young massive stars in the Galactic center form in isolation; some of them seem to be the remnants of dissolved clusters or stellar associations. [ABSTRACT FROM AUTHOR]

Context. Hypervelocity stars (HVSs) unbound to the Galaxy can be formed with extreme stellar interactions, for example close encounters with supermassive black holes or in massive star clusters, supernova explosions in binary systems, or the stripping of dwarf galaxies. Observational evidence comes from measurements of radial velocities (RVs) of objects crossing the outer Galactic halo and of tangential velocities based on high proper motions (HPMs) and distances of relatively nearby stars. Aims. I searched for new nearby HVS candidates and reviewed known objects using their Gaia astrometric measurements. Methods. Candidates were selected with significant Gaia parallaxes of >0.1mas, proper motions of >20 mas yr−1, and computed galactocentric tangential velocities vtan_g>500 km s−1. The DR2 and DR3 samples of several thousand HVS candidates were studied with respect to their proper motions, sky distribution, number of observations, location in crowded fields, colour-magnitude diagrams, selection effects with magnitude, and RVs in DR3. The 72 most extreme (vtan_g>700 km s−1) and nearest (within 4kpc) DR3 HVS candidates were investigated with respect to detected close neighbours, flags, and astrometric quality parameters of objects of similar magnitudes in DR3. The quality checks involved HPM objects in a global comparison and all objects in the vicinity of each target. Results. Spurious HPMs in the Galactic centre region led to false HVS interpretations in Gaia DR2 and are still present in DR3, although to a lesser extent. Otherwise there is good agreement between the HPMs of HVS candidates in DR2 and DR3. However, HVS candidates selected from DR2 tend to have larger parallaxes, and hence lower tangential velocities in DR3. Most DR3 RVs are much lower than the tangential velocities, indicating that the DR3 HVS candidates are still affected by underestimated parallaxes. None of the 72 extreme nearby DR3 HVS candidates, including three D6 stars, passed all the quality checks. Their tangential velocities may turn out to be lower, but at least some of them still appear unbound to the Galaxy. [ABSTRACT FROM AUTHOR]

This work presents the results of a kinematic analysis of the Galaxy that uses a new model as applied to the newest available Gaia data. We carry out the Taylor decomposition of the velocity field up to second order for 18 million high luminosity stars (i.e. OBAF-type stars, giants, and subgiants) from the Gaia DR3 data. We determine the components of mean stellar velocities and their first and second partial derivatives (relative to cylindrical coordinates) for more than 28 thousand points in the plane of our Galaxy. We estimate Oort's constants A, B, C , and K and other kinematics parameters and map them as a function of Galactocentric coordinates. The values found confirm the results of our previous works and are in excellent agreement with those obtained by other authors in the solar neighbourhood. In addition, the introduction of second order partial derivatives of the stellar velocity field allows us to determine the values of the vertical gradient of the Galaxy azimuthal, radial, and vertical velocities. Also, we determine the mean of the Galaxy rotation curve for Galactocentric distances from 4 to 18 kpc by averaging Galactic azimuths in the range −30° < θ <  + 30° about the direction Galactic Centre – Sun – Galactic anticentre. Maps of the velocity components and of their partial derivatives with respect to coordinates within 10 kpc of the Sun reveal complex substructures, which provide clear evidence of non-axisymmetric features of the Galaxy. Finally, we show evidence of differences in the Northern and Southern hemispheres stellar velocity fields. [ABSTRACT FROM AUTHOR]

We have determined positions, proper motions, and parallaxes of 77 millisecond pulsars (MSPs) from ∼3 yr of MeerKAT radio telescope observations. Our timing and noise analyses enable us to measure 35 significant parallaxes (12 of them for the first time) and 69 significant proper motions. Eight pulsars near the ecliptic have an accurate proper motion in ecliptic longitude only. PSR J0955−6150 has a good upper limit on its very small proper motion (<0.4 mas yr−1). We used pulsars with accurate parallaxes to study the MSP velocities. This yields 39 MSP transverse velocities, and combined with MSPs in the literature (excluding those in Globular Clusters) we analyse 66 MSPs in total. We find that MSPs have, on average, much lower velocities than normal pulsars, with a mean transverse velocity of only 78(8) km s−1 (MSPs) compared with 246(21) km s−1 (normal pulsars). We found no statistical differences between the velocity distributions of isolated and binary MSPs. From Galactocentric cylindrical velocities of the MSPs, we derive 3D velocity dispersions of σρ, σϕ, σ z  = 63(11), 48(8), 19(3) km s−1. We measure a mean asymmetric drift with amplitude 38(11) km s−1, consistent with expectation for MSPs, given their velocity dispersions and ages. The MSP velocity distribution is consistent with binary evolution models that predict very few MSPs with velocities >300 km s−1 and a mild anticorrelation of transverse velocity with orbital period. [ABSTRACT FROM AUTHOR]

We search for high-velocity stars in the inner region of the Galactic bulge using a selected sample of red clump stars. Some of those stars might be considered hypervelocity stars (HVSs). Even though the HVSs ejection relies on an interaction with the supermassive black hole (SMBH) at the centre of the Galaxy, there are no confirmed detections of HVSs in the inner region of our Galaxy. With the detection of HVSs, ejection mechanism models can be constrained by exploring the stellar dynamics in the Galactic centre through a recent stellar interaction with the SMBH. Based on a previously developed methodology by our group, we searched with a sample of preliminary data from version 2 of the Vista Variables in the Via Lactea (VVV) Infrared Astrometric Catalogue (VIRAC2) and Gaia DR3 data, including accurate optical and near-infrared proper motions. This search resulted in a sample of 46 stars with transverse velocities larger than the local escape velocity within the Galactic bulge, of which four are prime candidate HVSs with high-proper motions consistent with being ejections from the Galactic centre. Adding to that, we studied a sample of reddened stars without a Gaia DR3 counterpart and found 481 stars with transverse velocities larger than the local escape velocity, from which 65 stars have proper motions pointing out of the Galactic centre and are candidate HVSs. In total, we found 69 candidate HVSs pointing away from the Galactic centre with transverse velocities larger than the local escape velocity. [ABSTRACT FROM AUTHOR]

Hydrogen-deficient stars include the cool R CrB variable (RCBs) and hydrogen-deficient carbon (HdCs) giants through extreme helium stars (EHes) to the very hot helium-rich subdwarfs (He-sdO and O(He) stars) and white dwarfs. With surfaces rich in helium, nitrogen, and carbon, their origins have been identified with the merger of two white dwarfs. Using Gaia to focus on the EHes, we aim to identify progenitor populations and test the evolution models. Gaia DR3 measurements and ground-based radial velocities have been used to compute Galactic orbits using galpy. Each orbit has been classified by population; EHe stars are found in all of the thin disc, thick disc, halo, and bulge, as are RCB, HdC, and He-sdO stars. Spectral energy distributions were constructed for all EHes, to provide angular diameters, and hence radii and luminosities. The EHes fall into two luminosity groups divided at |$L \approx 2500 \, {\rm L_{\odot }}$|⁠. This supports theory for the origin of EHes, and is the strongest confirmation so far in terms of luminosity. The lower luminosity EHes correspond well with the post-merger evolution of a double helium white dwarf binary. Likewise, the higher luminosity EHes match the post-merger evolution of a carbon/oxygen plus helium white dwarf binary. In terms of parent populations, current models predict that double white dwarf mergers should occur in all Galactic populations, but favour mergers arising from recent star formation (i.e. thin disc), whereas the statistics favour an older epoch (i.e. thick disc). [ABSTRACT FROM AUTHOR]

Globular clusters (GCs) have been proposed as promising sites for discovering intermediate-mass black holes (IMBHs), offering the possibility to gain crucial insights into the formation and evolution of these elusive objects. The Galactic GC 47 Tucanae (also known as NGC 104) has been suggested as a potential IMBH host, yet previous studies have yielded conflicting results. Therefore, we present a set of self-consistent dynamical models based on distribution functions (DFs) that depend on action integrals to assess the presence (or absence) of an IMBH in 47 Tucanae. Leveraging the state-of-the-art Multi Unit Spectroscopic Explorer (MUSE) and Hubble Space Telescope (HST) data, we analyzed the three-dimensional (3D) kinematics of the cluster's central regions, fitting individual star velocities down to the sub-arcsec scale (approximately 10−2 pc). According to our analysis, the inner kinematics of 47 Tucanae is incompatible with a central BH more massive than 578 M⊙ (at 3σ). This is the most stringent upper limit placed thus far on the mass of a putative IMBH in 47 Tucanae via a dynamical study. [ABSTRACT FROM AUTHOR]

We investigate the spatial structure and dynamical properties of the young (∼60 Myr), richly populated open cluster NGC 6649 based on reliable cluster members from Gaia -DR3 data. The density-based spatial clustering of applications with noise is used to estimate membership probabilities for 25045 sample stars in the field of the cluster. More than 1000 likely cluster members down to G ∼ 20 mag can be identified based on the membership probabilities. The cluster members of NGC 6649 extend up to a radial distance of |${\sim} 21\,\mathrm{ arcmin}$| (12 pc), which is much larger than previously thought. We estimate a core radius of |$R_{\mathrm{ c}}=2\overset{^{\prime }}{.}38\pm 0\overset{^{\prime }}{.}06$| (1.4 pc) and a tidal radius of |$R_{\mathrm{ t}}=35\overset{^{\prime }}{.}85\pm 6\overset{^{\prime }}{.}64$| (20.9 pc) for the cluster, indicating that the cluster has formed a dense core and an extended halo. We find that mass segregation has taken place within the cluster. We also detect two small and sparse clusters (C1 and C2) around NGC 6649. C1 has recently been reported, but C2 is detected for the first time. The Gaussian mixture model clustering method is used to identify 184 and 60 members for C1 and C2, respectively. We find that NGC 6649, C1, and C2 are located at similar distances. Moreover, we estimate a small spatial separation of ∼11.6 pc between C1 and NGC 6649, which is close to or slightly smaller than the detectable radius of NGC 6649. This indicates that NGC 6649 and C1 are physically connected, and there may exist strong tidal interactions and mass transfer between the two clusters. C2 shows a significant lack of faint members (G > 17 mag), probably due to star evaporation and tidal stripping. Our results show that NGC 6649 may be the dominant cluster in a triple-cluster system (NGC 6649+C1+C2). [ABSTRACT FROM AUTHOR]

Be stars make up almost 20 per cent of the B star population, and are rapidly rotating stars surrounded by a disc; however the origin of this rotation remains unclear. Mass transfer within close binaries provides the leading hypothesis, with previous detections of stripped companions to Be stars supporting this. Here, we exploit the exquisite astrometric precision of Gaia to carry out the largest to date comparative study into the binarity of matched samples of nearby B and Be stars from the Bright Star Catalogue. By utilizing new 'proper motion anomaly' values, derived from Gaia DR2 and DR3 astrometric data alongside previous values calculated using Hipparcos and Gaia data, and the Gaia -provided RUWE, we demonstrate that we can identify unresolved binaries down to separations of 0.02 arcsec. Using these measures, we find that the binary fractions of B and Be stars are similar between 0.04 and 10 arcsec, but the Be binary fraction is significantly lower than that of the B stars for separations below 0.04 arcsec. As the separation range of these 'missing' binaries is too large for mass transfer, and stripped companions are not retrieved by these measures, we suggest the companions migrate inwards via binary hardening within a triple system. This confirms statistically for the first time the hypothesis that binary interaction causes the Be phenomenon, with migration causing the dearth of Be binaries between 0.02 and 0.04 arcsec. Furthermore, we suggest that triplicity plays a vital role in this migration, and thus in the formation of Be stars as a whole. [ABSTRACT FROM AUTHOR]

We fit various colour–magnitude diagrams (CMDs) of the Galactic globular clusters NGC 6397 and NGC 6809 (M55) by isochrones from the Dartmouth Stellar Evolution Database (DSED) and Bag of Stellar Tracks and Isochrones (BaSTI) for α–enhanced [α/Fe] = +0.4. For the CMDs, we use data sets from Hubble Space Telescope, Gaia, Visible and Infrared Survey Telescope for Astronomy , and other sources utilizing 32 and 23 photometric filters for NGC 6397 and NGC 6809, respectively, from the ultraviolet to mid-infrared. We obtain the following characteristics for NGC 6397 and NGC 6809, respectively: metallicities [Fe/H] = −1.84 ± 0.02 ± 0.1 and −1.78 ± 0.02 ± 0.1 (statistic and systematic uncertainties); distances 2.45 ± 0.02 ± 0.06 and 5.24 ± 0.02 ± 0.18 kpc; ages 12.9 ± 0.1 ± 0.8 and 13.0 ± 0.1 ± 0.8 Gyr; reddenings E (B − V) = 0.178 ± 0.006 ± 0.01 and 0.118 ± 0.004 ± 0.01 mag; extinctions A V = 0.59 ± 0.01 ± 0.02 and 0.37 ± 0.01 ± 0.04 mag; and extinction-to-reddening ratio |$R_\mathrm{V}=3.32^{+0.32}_{-0.28}$| and |$3.16^{+0.66}_{-0.56}$|⁠. Our estimates agree with most estimates from the literature. BaSTI gives systematically higher [Fe/H] and lower reddenings than DSED. Despite nearly the same metallicity, age, and helium enrichment, these clusters show a considerable horizontal branch (HB) morphology difference, which must therefore be described by another parameter. This parameter must predominantly explain why the least massive HB stars (0.58–0.63 solar masses) are only found within NGC 6809. Probably they have been lost by the core-collapse cluster NGC 6397 during its dynamical evolution and mass segregation. In contrast, NGC 6809 has a very low central concentration and, hence, did not undergo this process. [ABSTRACT FROM AUTHOR]

To study the dust dynamics in the dust-trapping vortices in the protoplanetary disk around MWC 758, we analyzed the 1.3 mm continuum images of the MWC 758 disk obtained with the Atacama Large Millimeter/submillimeter Array in 2017 and 2021. We detect proper motions of 22 mas and 24 mas in the two dust clumps at radii of 0.″32 and 0.″54 in the disk on the plane of the sky, respectively. On the assumption that the dust clumps are located in the disk midplane, the velocities of the observed proper motions along the azimuthal direction of the inner and outer dust clumps are sub- and super-Keplerian, respectively, and both have angular velocities corresponding to the Keplerian angular velocity at a radius of 0.″46 ± 0.″04. This deviation from the Keplerian motion is not expected in the conventional theory of vortices formed by the Rossby wave instability. The observed non-Keplerian proper motions of the dust clumps are unlikely due to the disk warp and eccentricity, nor are they associated with any predicted planets. The two dust clumps are likely spatially coincident with the infrared spirals. In addition, we detect the changes in the intensity profiles of the dust clumps over the 4 yr span. Therefore, we suggest that the observed proper motions are possibly due to changes in the density distributions in the dust clumps caused by their interaction with the spirals in the disk.

In addition to being the most magnetic objects in the known Universe, magnetars are the only objects observed to generate fast-radio-burst-like emissions. The formation mechanism of magnetars is still highly debated and may potentially be probed with the magnetar velocity distribution. We carried out a 3 yr long astrometric campaign on Swift J1818.0−1607, the fastest-spinning magnetar, using the Very Long Baseline Array. After applying the phase-calibrating 1D interpolation strategy, we obtained a small proper motion of 8.5 mas yr ^−1 mag and a parallax of 0.12 ± 0.02 mas (uncertainties at 1 σ confidence throughout the Letter) for Swift J1818.0−1607. The latter is the second magnetar parallax and is among the smallest neutron star parallaxes ever determined. From the parallax, we derived the distance ${9.4}_{-1.6}^{+2.0}$ kpc, which locates Swift J1818.0−1607 at the far side of the Galactic central region. Combined with the distance, the small proper motion leads to a transverse peculiar velocity ${v}_{\perp }={48}_{-16}^{+50}$ km s ^−1 —a new lower limit to magnetar v _⊥ . Incorporating previous v _⊥ estimates of seven other magnetars, we acquired ${v}_{\perp }={149}_{-68}^{+132}$ km s ^−1 for the sample of astrometrically studied magnetars, corresponding to the three-dimensional space velocity $\sim {190}_{-87}^{+168}$ km s ^−1 , smaller than the average level of young pulsars. Additionally, we found that the magnetar velocity sample does not follow the unimodal young pulsar velocity distribution reported by Hobbs et al. at >2 σ confidence, while loosely agreeing with more recent bimodal young pulsar velocity distributions derived from relatively small samples of quality astrometric determinations.

We use the Subaru Hyper Suprime-Cam (HSC) data to study structural parameters and systemic proper motion of the Sextans dwarf spheroidal galaxy at the heliocentric distance of 86 kpc, which is one of the most important targets for studies of dark matter nature and galaxy formation physics. Thanks to the superb image quality and wide area coverage, the HSC data enable a secure selection of member star candidates based on the colour–magnitude cut, yielding about 10 000 member candidates at magnitudes down to i ∼ 24. We use a likelihood analysis of the two-dimensional distribution of stars to estimate the structural parameters of Sextans taking into account the contamination of foreground halo stars in the Milky Way, and find that the member star distribution is well fitted by an elliptical King profile with ellipticity ϵ ≃ 0.25 and the core and tidal radii of R c = (368.4 ± 8.5) pc and R t = (2.54 ± 0.046) kpc, respectively. Then using the two HSC data sets of 2.66 yr time baseline on average, we find the systemic proper motions of Sextans to be (μα, μδ) = (−0.448 ± 0.075, 0.058 ± 0.078) mas yr−1, which are consistent with some of the previous works using the Gaia data of relatively bright member stars in Sextans. Thus, our results give a demonstration that ground-based, large-aperture telescope data that cover a wide solid angle of the sky and have a long time baseline, such as the upcoming data from Legacy Survey of Space and Time (LSST), can be used to study systemic proper motions of dwarf galaxies. [ABSTRACT FROM AUTHOR]

We compile a catalogue of low-mass and high-mass X-ray binaries, some recently reported binaries that likely host a neutron star (NS) or a black hole (BH), and binary pulsars (a pulsar and a non-degenerated companion) that have measured systemic radial velocities (γ). Using Gaia and radio proper motions together with γ, we integrate their Galactic orbits and infer their post-supernova (post-SN) 3D peculiar velocities (⁠|$v_\mathrm{pec }^{z=0}$| at Galactic plane crossing); these velocities bear imprints of the natal kicks that compact objects received at birth. With the sample totalling 85 objects, we model the overall distribution of |$v_\mathrm{pec }^{z=0}$| and find a two-component Maxwellian distribution with a low- (⁠|$\sigma _v \approx 21\, \mathrm{km s^{-1}}$|⁠) and a high-velocity (⁠|$\sigma _v \approx 107\, \mathrm{km s^{-1}}$|⁠) component. A further comparison between distributions of binary subgroups suggests that binaries hosting high-mass donors/luminous companions mostly have |$v_\mathrm{pec }^{z=0}\lesssim 100\, \mathrm{km s^{-1}}$|⁠ , while binaries with low-mass companions exhibit a broader distribution that extends up to |$\sim 400\, \mathrm{km s^{-1}}$|⁠. We also find significant anticorrelations of |$v_\mathrm{pec }^{z=0}$| with binary total mass (M tot) and orbital period (P orb), at over 99 per cent confidence. Specifically, our fit suggests |$v_\mathrm{pec }^{z=0}\propto M_\mathrm{tot}^{-0.5}$| and |$v_\mathrm{pec }^{z=0}\propto P_\mathrm{orb}^{-0.2}$|⁠. Discussions are presented on possible interpretation of the correlations in the context of kinematics and possible biases. The sample should enable a range of follow-up studies on compact object binary kinematics and evolution. [ABSTRACT FROM AUTHOR]

To accurately link the radio and optical Celestial Reference Frames (CRFs) at optical bright end, i.e. with GaiaG -band magnitude |$\lesssim$| 13, increasing number and improving sky distribution of radio stars with accurate astrometric parameters from both Very Long Baseline Interferometry (VLBI) and Gaia measurements are mandatory. We selected two radio stars HD 199178 and AR Lacertae as the target for a pilot program for the frame link, using the Very Long Baseline Array at 15 GHz at six epochs spanning about 1 yr, to measure their astrometric parameters. The measured parallax of HD 199178 is 8.949 ± 0.059 mas and the proper motion is μ αcos δ = 26.393 ± 0.093 and μ δ = −0.950 ± 0.083 mas yr−1, while the parallax of AR Lac is 23.459 ± 0.094 mas and the proper motion is μ αcos δ = −51.906 ± 0.138 and μ δ = 46.732 ± 0.131 mas yr−1. Our VLBI measured astrometric parameters have accuracies about 4–5 times better than the corresponding historic VLBI measurements and comparable accuracies with those from Gaia , validating the feasibility of frame link using radio stars. With the updated astrometric parameters for these two stars, there is a ∼25 per cent reduction of the uncertainties on the Y -axis for both orientation and spin parameters. [ABSTRACT FROM AUTHOR]

This article describes the dynamical behavior of any linear control system on the group of proper motions SE(2). It characterizes the controllability property and the control sets of these systems. [ABSTRACT FROM AUTHOR]

Astrometric surveys can be used to constrain the stochastic gravitational wave background (SGWB) at very low frequencies. We use proper motion data provided by Gaia DR3 to fit a generic dipole+quadrupole field. We analyse several quasar-based data sets and discuss their purity and idoneity to set constraints on gravitational waves. For the cleanest data set, we derive an upper bound on the (frequency-integrated) energy density of the SGWB |$h_{70}^2\Omega _{\rm GW}\lesssim 0.087$| for 4.2 × 10−18 Hz ≲ f ≲ 1.1 × 10−8 Hz. We also reanalyse previous VLBI-based data to set the constraint |$h_{70}^2\Omega _{\rm GW}\lesssim 0.024$| for 5.8 × 10−18 Hz ≲ f ≲ 1.4 × 10−9 Hz under the same formalism, standing as the best astrometric constraint on GWs. Based on our results, we discuss the potential of future Gaia data releases to impose tighter constraints. [ABSTRACT FROM AUTHOR]

Observations and theories reveal that star-forming regions eject their members during the early evolution of the cluster. We present an initial search for stars ejected from the young cluster IC 348 based on an analysis of Gaia EDR3 astrometry of previously identified members. We identify six candidate walkaway stars with projected velocities of 5–8 km s−1. The candidate with properties most consistent with a walkaway star, V701 Per, was likely ejected from a small group to the west of IC 348. Three candidates are currently located or have proper motions that trace back to the center of IC 348, though their classification as walkaway stars may also be affected by large astrometric errors. Two other candidate walkaway stars were likely born in other subgroups around IC 348 within the Perseus star-forming region. Since our search criteria are restricted to previously identified members of IC 348, we did not identify any stars that were ejected at high velocities and traveled far from the cluster. Our study suggests that ejections have potentially occurred in both the central region of IC 348 and other subgroups in Perseus. Extending the search range will lead to more candidate ejections and constrain the dynamic properties and evolution of IC 348 by comparing with simulations. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the spatial structure and dynamical state of the poorly studied old (∼ 4 Gyr) open cluster Trumpler 19 based on reliable cluster members from G a i a -DR3. The DBSCAN clustering algorithm is used to estimate membership probabilities and select likely cluster members in a normalized 5D parametric space. We identify 859 likely cluster members down to G ∼ 20 mag in the field of the cluster. We estimate a cluster radius of R c l ∼ 18 arcmin (13.1 pc) based on the radial distribution of the cluster members. We find that the cluster is deficient in faint and low-mass stars (G > 18.5 mag) and strong mass segregation has taken place within the cluster. We estimate a core radius of R c = 3.0 ± 0.2 arcmin (2.2 ± 0.1 pc) and a tidal radius of R t = 30.7 ± 9.0 arcmin (22.3 ± 6.5 pc) for the cluster. We estimate a concentration parameter of log (R t / R c) ∼ 1.0 , indicating that Trumpler 19 has formed a clear core-halo structure due to dynamical evolution. We investigate the dynamical state of the cluster using 26 blue stragglers (BSs) identified in the cluster. We find that the BSs extend to a radius of ∼ 10.5 arcmin (3.5 R c ), and 16 of them (∼ 62 % ) are located in the core radius of the cluster. We also find that the BSs are significantly more concentrated than the red giant branch (RGB) and main-sequence (MS) stars, indicating that Trumpler 19 is a dynamically old cluster. [ABSTRACT FROM AUTHOR]

Since the hypervelocity stars were discovered in the Milky Way, various mechanisms have been proposed to explain these runaway stars. Up to now, however, the dominant ejected mechanism of hypervelocity stars is still unclear. As the largest and closest face-on satellite galaxy of the Milky Way, the Large Magellanic Cloud (LMC) serves as a good target, allowing us to study this issue. Based on the high-precision astrometric parameters provided by Gaia Data Release 3, we researched the proper motions of 3119 massive O–B2-type stars in the LMC and identified 98 (∼3%) as having the ability to escape from the LMC, with more than 50% confidence. Furthermore, by investigating the characteristics of the identified massive stars and adopting a regression analysis, we find that the dynamic ejection scenario might be the dominant mechanism of the massive hypervelocity runaway stars in the LMC. [ABSTRACT FROM AUTHOR]

Mrk 1018 is a nearby changing-look active galactic nucleus (AGN) that has oscillated between spectral Type 1.9 and Type 1 over a period of 40 yr. Recently, a recoiling supermassive black hole (rSMBH) scenario has been proposed to explain the spectral and flux variability observed in this AGN. Detections of rSMBHs are important for understanding the processes by which SMBH binaries merge and how rSMBHs influence their galactic environment through feedback mechanisms. However, conclusive identification of any rSMBHs has remained elusive to date. In this paper, we present an analysis of 6.5 yr of multifrequency Very Long Baseline Array monitoring of Mrk 1018. We find that the radio emission is compact down to 2.4 pc, and it displays flux density and spectral variability over the length of our campaign, typical of a flat-spectrum radio core. We observe proper motion in RA of the radio core at −36.4 ± 8.6 μ as yr−1 (4.2 σ), or 0.10 c ± 0.02 c at the redshift of Mrk 1018. No significant proper motion is found in DEC (31.3 ± 25.1 μ as yr−1). We discuss possible physical mechanisms driving the proper motion, including an rSMBH. We conclude that the apparent velocity we observe of the VLBI radio core is too high to reconcile with theoretical predictions of rSMBH velocities and that the proper motion is most likely dominated by an unresolved, outflowing jet component. Future observations may yet reveal the true nature of Mrk 1018. However, our observations are not able to confirm it as a true rSMBH. [ABSTRACT FROM AUTHOR]

Recent studies of nearby globular clusters have discovered excess dark mass in their cores, apparently in an extended distribution, and simulations indicate that this mass is composed mostly of white dwarfs (respectively stellar-mass black holes) in clusters that are core collapsed (respectively with a flatter core). We perform mass-anisotropy modelling of the closest globular cluster, M4, with intermediate slope for the inner stellar density. We use proper motion data from Gaia Early Data Release 3 (EDR3) and from observations by the Hubble Space Telescope. We extract the mass profile employing Bayesian Jeans modelling, and check our fits with realistic mock data. Our analyses return isotropic motions in the cluster core and tangential motions (β ≈ −0.4 ± 0.1) in the outskirts. We also robustly measure a dark central mass of roughly |$800\pm 300 \, \rm M_\odot$| , but it is not possible to distinguish between a point-like source, such as an intermediate-mass black hole (IMBH), and a dark population of stellar remnants of extent |${\approx} 0.016\, {\rm pc} \simeq 3300\, {\rm au}$|. However, when removing a high-velocity star from the cluster centre, the same mass excess is found, but more extended (⁠|${\sim} 0.034\, {\rm pc} \approx 7000\, {\rm au}$|). We use Monte Carlo N -body models of M4 to interpret the second outcome, and find that our excess mass is not sufficiently extended to be confidently associated with a dark population of remnants. Finally, we discuss the feasibility of these two scenarios (i.e. IMBH versus remnants), and propose new observations that could help to better grasp the complex dynamics in M4's core. [ABSTRACT FROM AUTHOR]

The application of very long baseline interferometry (VLBI) to the search for extraterrestrial intelligence (SETI) has been limited to date, despite the technique offering many advantages over traditional single-dish SETI observations. In order to further develop interferometry for SETI, we used the European VLBI network (EVN) at 21 cm to observe potential secondary phase calibrators in the Kepler field. Unfortunately, no secondary calibrators were detected. However, a VLBA primary calibrator in the field, J1926+4441, offset only ∼1.88 arcmin from a nearby exoplanet Kepler-111 b, was correlated with high temporal |$\left(0.25 \, \rm {s}\right)$| and spectral |$\left(16384 \times 488\, \rm {Hz \ channels}\right)$| resolution. During the analysis of the high-resolution data, we identified a spectral feature that was present in both the auto and cross-correlation data with a central frequency of 1420.424 ± 0.0002 MHz and a width of 0.25 MHz. We demonstrate that the feature in the cross-correlations is an artefact in the data, associated with a significant increase in each telescope's noise figure due to the presence of H  i in the beam. This would typically go unnoticed in data correlated with standard spectral resolution. We flag (excluded from the subsequent analysis) these channels and phase rotate the data to the location of Kepler-111 b aided by the Gaia catalogue and search for signals with |$\rm {SNR}\gt 7$|. At the time of our observations, we detect no transmitters with an equivalent isotropically radiated power ≳4 × 1015 W. [ABSTRACT FROM AUTHOR]

The majority of massive stars are born in binaries, and most unbind upon the first supernova. With precise proper motion surveys such as Gaia , it is possible to trace back the motion of stars in the vicinity of young remnants to search for ejected companions. Establishing the fraction of remnants with an ejected companion, and the photometric and kinematic properties of these stars, offers unique insight into supernova progenitor systems. In this paper, we employ binary population synthesis to produce kinematic and photometric predictions for ejected secondary stars. We demonstrate that the unbound neutron star velocity distribution from supernovae in binaries closely traces the input kicks. Therefore, the observed distribution of neutron star velocities should be representative of their natal kicks. We evaluate the probability for any given filter, magnitude limit, minimum measurable proper motion (as a function of magnitude), temporal baseline, distance, and extinction that an unbound companion can be associated with a remnant. We compare our predictions with results from previous companion searches, and demonstrate that the current sample of stars ejected by the supernova of their companion can be increased by a factor of ∼5–10 with Gaia data release 3. Further progress in this area is achievable by leveraging the absolute astrometric precision of Gaia , and by obtaining multiple epochs of deep, high resolution near-infrared imaging with the Hubble Space Telescope ,  JWST and next-generation wide-field near-infrared observatories such as Euclid or the Nancy Grace Roman Space Telescope. [ABSTRACT FROM AUTHOR]

This paper is an extension of the study by Khamitov et al. (2022) with regard to the catalog and the astrophysical interpretation of the imitation of significant proper motions in active galactic nuclei (AGNs) and quasars based on data from the Gaia space observatory. We present a sample of SRG/eROSITA X-ray sources in the eastern Galactic hemisphere () having significant proper motions in the Gaia EDR3 measurements with the confirmed extragalactic nature of the objects. The catalog consists of 248 extragalactic sources with spectroscopically measured redshifts. The catalog includes all of the objects available in the SIMBAD database and coincident with the identified optical counterpart within 0.5 arcsec. Eighteen sources with spectroscopically measured redshifts from observations with the Russian–Turkish 1.5-m telescope RTT-150 (Khamitov et al. 2022) have been additionally included in the catalog. The sources in the catalog are AGNs of various types (Sy1, Sy2, LINER), quasars, radio galaxies, and star-forming galaxies. The imitation of significant proper motions can be explained (by the VIM effect previously known in astrometry) by the presence of transient events on the line of sight in the vicinity of AGNs and quasars (within the Gaia optical resolution element). Among such astrophysical events are supernova outbursts, tidal disruption events in binary AGNs, the variability of high-mass supergiants, the presence of OB associations against the background of AGNs with a variable brightness, etc. The model of outbursts with a fast rise–exponential decay profile allows the variable positional parameters of most sources observed in Gaia to be described. This approach can be used as an independent way of detecting transient events in the vicinity of AGNs (on scales of several hundred parsecs in the plane of the sky) based on data from the SRG/eROSITA catalogs of X-ray sources and the optical Gaia catalog. [ABSTRACT FROM AUTHOR]

We have solved the Ogorodnikov–Milne stellar-kinematics equations in the Galactic rectangular coordinate system based on the total velocities for a special sample of stars with radial velocities from the final Gaia Data Release 3 catalogue. We have found the region of applicability of the linear model and the regions that it describes poorly. We have constructed a second-order model that takes into account the peculiarities of stellar kinematics more accurately and showed its applicability for stars at distances up to 5 kpc. [ABSTRACT FROM AUTHOR]

Several observational studies have shown that many Galactic globular clusters (GCs) are characterized by internal rotation. Theoretical studies of the dynamical evolution of rotating clusters have predicted that, during their long-term evolution, these stellar systems should develop a dependence of the rotational velocity around the cluster's centre on the mass of stars, with the internal rotation increasing for more massive stars. In this paper, we present the first observational evidence of the predicted rotation–mass trend. In our investigation, we exploited the Gaia Data Release 3 catalogue of three GCs: NGC 104 (47 Tuc), NGC 5139 (ω Cen), and NGC 5904 (M 5). We found clear evidence of a cluster rotation–mass relation in 47 Tuc and M 5, while in ω Cen, the dynamically youngest system among the three clusters studied here, no such trend was detected. [ABSTRACT FROM AUTHOR]

Abstract We present the first proper-motion (PM) measurements for the dwarf elliptical galaxies NGC 147 and NGC 185, two satellite galaxies of M31, using multiepoch Hubble Space Telescope (HST) imaging data with time baselines of ∼8 yr. For each galaxy, we take an error-weighted average of measurements from HST Advanced Camera for Surveys/WFC and Wide Field Camera 3/UVIS to determine the PMs. Our final results for the PMs are (μ W, μ N)N147 = (−0.0232, 0.0378) ± (0.0143, 0.0146) mas yr−1 for NGC 147 and (μ W, μ N)N185 = (−0.0242, 0.0058) ± (0.0141, 0.0147) mas yr−1 for NGC 185. The 2D direction of motion for NGC 147 about M31 is found to be aligned with its tidal tails. The 3D positions and velocities of both galaxies are transformed into a common M31-centric coordinate system to study the detailed orbital histories of the combined M31+NGC 147+NGC 185 system via numerical orbit integration. We find that NGC 147 (NGC 185) had its closest passage to M31 0.3–0.5 Gyr (≳1.6 Gyr) within the past 6 Gyr at distances of ∼70 kpc (70–260 kpc). The pericentric times of NGC 147/NGC 185 correlate qualitatively well with the presence/absence of tidal tails seen around the galaxies. Our PMs show that the orbital poles of NGC 147, as well as NGC 185, albeit to a lesser degree, agree within the uncertainties with the normal of the Great Plane of Andromeda (GPoA). These are the first measurements of the 3D angular momentum vector of any satellite identified as original GPoA members. Our results strengthen the hypothesis that the GPoA may be a dynamically coherent entity. We revisit previous claims that NGC 147 and NGC 185 are binary galaxies and conclude that it is very unlikely that the two galaxies were ever gravitationally bound to each other.

CatWISE is a program to catalog sources selected from combined WISE and NEOWISE all-sky survey data at 3.4 and 4.6 μm (W1 and W2). The CatWISE Preliminary Catalog consists of 900,849,014 sources measured in data collected from 2010 to 2016. This data set represents four times as many exposures and spans over 10 times as large a time baseline as that used for the AllWISE Catalog. CatWISE adapts AllWISE software to measure the sources in coadded images created from six-month subsets of these data, each representing one coverage of the inertial sky, or epoch. The catalog includes the measured motion of sources in eight epochs over the 6.5 yr span of the data. From comparison to Spitzer, signal-to-noise ratio = 5 limits in magnitudes in the Vega system are W1 = 17.67 and W2 = 16.47, compared to W1 = 16.96 and W2 = 16.02 for AllWISE. From comparison to Gaia, CatWISE positions have typical accuracies of 50 mas for stars at W1 = 10 mag and 275 mas for stars at W1 = 15.5 mag. Proper motions have typical accuracies of 10 mas yr-1 and 30 mas yr-1 for stars with these brightnesses, an order of magnitude better than from AllWISE. The catalog is available in the WISE/NEOWISE Enhanced and Contributed Products area of the NASA/IPAC Infrared Science Archive.

We present the results of determining the coordinates of the vertices of various stellar systems, the centroids of which are located in the Galactic plane. To do this, the positions, parallaxes, proper motions, and radial velocities of red giants and subgiants contained in the Gaia  DR3 catalogue have been used. When determining the components of the deformation velocity tensors in local coordinate systems, we found the coordinates of the vertices of the stellar systems under study. It turned out that there is a complex dependence of vertex deviations lxy in Galactocentric cylindrical (R, θ) and Galactic rectangular (X, Y) coordinates. Based on the approach proposed in this paper, heliocentric distances to vertices have been determined for the first time. The results obtained show that in addition to the fact that the angular coordinates of the Galactic Centre and the vertices of stellar systems do not coincide, their heliocentric distances do not coincide as well. This presumably indicates that there are structures in the Galaxy that noticeably affect its axisymmetry. [ABSTRACT FROM AUTHOR]

Understanding the natal kicks received by neutron stars (NSs) during formation is a critical component of modelling the evolution of massive binaries. Natal kicks are an integral input parameter for population synthesis codes, and have implications for the formation of double NS systems and their subsequent merger rates. However, many of the standard observational kick distributions that are used are obtained from samples created only from isolated NSs. Kick distributions derived in this way overestimate the intrinsic NS kick distribution. For NSs in binaries, we can only directly estimate the effect of the natal kick on the binary system, instead of the natal kick received by the NS itself. Here, for the first time, we present a binary kick distribution for NSs with low-mass companions. We compile a catalogue of 145 NSs in low-mass binaries with the best available constraints on proper motion, distance, and systemic radial velocity. For each binary, we use a three-dimensional approach to estimate its binary kick. We discuss the implications of these kicks on system formation, and provide a parametric model for the overall binary kick distribution, for use in future theoretical modelling work. We compare our results with other work on isolated NSs and NSs in binaries, finding that the NS kick distributions fit using only isolated pulsars underestimate the fraction of NSs that receive low kicks. We discuss the implications of our results on modelling double NS systems, and provide suggestions on how to use our results in future theoretical works. [ABSTRACT FROM AUTHOR]

Context. The nuclear stellar disc (NSD) is a flat dense stellar structure at the heart of the Milky Way. Recent work has shown that analogous structures are common in the nuclei of external spiral galaxies, where there is evidence of an age gradient that indicates that they form inside-out. However, the characterisation of the age of the NSD stellar population along the line of sight is still missing due to its extreme source crowding and the high interstellar extinction towards the Galactic centre. Aims. We aim to characterise the age of the stellar population at different average Galactocentric NSD radii to investigate for the first time the presence of an age gradient along the line of sight. Methods. We selected two groups of stars at different NSD radii via their different extinction and proper motion distribution. We analysed their stellar population by fitting their de-reddened Ks luminosity functions with a linear combination of theoretical models. Results. We find significant differences in the stellar population at different NSD radii, indicating the presence of an age gradient along the line of sight. Our sample from the closest edge of the NSD contains a significant fraction (∼40% of its total stellar mass) of intermediate-age stars (2–7 Gyr) that is not present in the sample from stars deeper inside the NSD, in which ∼90% of the stellar mass is older than 7 Gyr. Our results suggest that the NSD age distribution is similar to the one found in external galaxies and they imply that bar-driven processes observed in external galaxies are similarly at play in the Milky Way. [ABSTRACT FROM AUTHOR]

Context. The Carina Nebula is one of the major massive star-forming regions in the Galaxy. Its relatively nearby distance (2.35 kpc) makes it an ideal laboratory for the study of massive star formation, structure, and evolution, both for individual stars and stellar systems. Thanks to the high-quality spectra provided by the Gaia-ESO survey and the LiLiMaRlin library, as well as Gaia EDR3 astrometry, a detailed and homogeneous spectroscopic characterization of its massive stellar content can be carried out. Aims. Our main objective is to spectroscopically characterize all massive members of the Carina Nebula in the Gaia-ESO survey footprint to provide an updated census of massive stars in the region and an updated estimate of the binary fraction of O stars. Methods. We performed accurate spectral classification using an interactive code that compares spectra with spectral libraries of OB standard stars, as well as line-based classic methods. We calculated membership using our own algorithm based on Gaia EDR3 astrometry. To check the correlation between the spectroscopic n-qualifier and the rotational velocity, we used a semi-automated tool for the line-broadening characterization of OB stars based on a combined Fourier transform and goodness-of-fit methodology. Results. The Gaia-ESO survey sample of massive OB stars in the Carina Nebula consists of 234 stars. The addition of brighter sources from the Galactic O-Star Spectroscopic Survey and additional sources from the literature allows us to create the most complete census of massive OB stars so far in the region. It contains a total of 316 stars, with 18 of them in the background and 4 in the foreground. Of the 294 stellar systems in Car OB1, 74 are of O type, 214 are of nonsupergiant B type, and 6 are of WR or nonO supergiant (II to Ia) spectral class. We identify 20 spectroscopic binary systems with an O-star primary, of which 6 are reported for the first time, and another 18 with a B-star primary, of which 13 are new detections. The average observed double-lined binary fraction of O-type stars in the surveyed region is 0.35, which represents a lower limit. We find a good correlation between the spectroscopic n-qualifier and the projected rotational velocity of the stars. The fraction of candidate runaways among the stars with and without the n-qualifier is 4.4% and 2.4%, respectively, although nonresolved double-lined binaries could be contaminating the sample of fast rotators. [ABSTRACT FROM AUTHOR]

Context. Star-forming regions are excellent benchmarks for testing and validating theories of star formation and stellar evolution. The Perseus star-forming region, being one of the youngest (< 10 Myr), closest (280−320 pc), and most studied in the literature, is a fundamental benchmark. Aims. We aim to study the membership, phase-space structure, mass, and energy (kinetic plus potential) distribution of the Perseus star-forming region using public catalogues (Gaia, APOGEE, 2MASS, and Pan-STARRS). Methods. We used Bayesian methodologies that account for extinction to identify the Perseus physical groups in the phase-space, retrieve their candidate members, derive their properties (age, mass, 3D positions, 3D velocities, and energy), and attempt to reconstruct their origin. Results. We identify 1052 candidate members in seven physical groups (one of them new) with ages between 3 and 10 Myr, dynamical super-virial states, and large fractions of energetically unbounded stars. Their mass distributions are broadly compatible with that of Chabrier for masses ≳0.1 M⊙ and do not show hints of over-abundance of low-mass stars in NGC 1333 with respect to IC 348. These groups' ages, spatial structure, and kinematics are compatible with at least three generations of stars. Future work is still needed to clarify if the formation of the youngest was triggered by the oldest. Conclusions. The exquisite Gaia data complemented with public archives and mined with comprehensive Bayesian methodologies allow us to identify 31% more members than previous studies, discover a new physical group (Gorgophone: 7 Myr, 191 members, and 145 M⊙), and confirm that the spatial, kinematic, and energy distributions of these groups support the hierarchical star formation scenario. [ABSTRACT FROM AUTHOR]

Context. The Carte du Ciel archive at the University of Helsinki enables us to see the sky as it was about 120 yr ago. The archive consists of single-exposure and triple-exposure plates between epochs 1896.8–1925.8. Aims. Our main aim is to find binary and multiple stars by combining Carte du Ciel and Gaia data. Methods. The plates were digitised with a commercial digital camera. We used Gaia data to calculate predicted coordinates of stars at the epoch of each plate. These stars were used as reference stars to fit astrometry for each plate, giving fitted coordinates for stars on the Carte du Ciel plates. If the predicted and fitted coordinates differed at a significant level, we classified the star as a non-single star, for which the proper motion values given in the Gaia catalogue can be unreliable. Results. We find that several astrometric quality indicators of Gaia indicate that the uncertainties of Gaia's single-star model fit are, in general, larger for our non-single-star candidates. The percentage of our non-single-star candidates, which are in the catalogues of known binary stars, is relatively low, ~10% at maximum. Conclusions. The combination of the Carte du Ciel and Gaia data can be used to identify candidates of non-single stars. We propose that the sources with a significant difference between the predicted and fitted coordinates are long-period binaries, although astrophysical and/or instrumental effects as origin for the coordinate difference cannot be excluded for individual cases. [ABSTRACT FROM AUTHOR]

The flattened, possibly co-rotating plane of satellite galaxies around Centaurus A, if more than a fortuitous alignment, adds to the pre-existing tension between the well-studied Milky Way and M31 planes and the Lambda cold dark matter model of structure formation. It was recently reported that the Centaurus A satellite plane (CASP) may be rotationally supported, but a further understanding of the system's kinematics is elusive in the absence of full 3D velocities. We constrain the transverse velocities of 27 satellites that would rotationally stabilize the Centaurus A plane, and classify the satellites by whether their possible orbits are consistent with the CASP. Five satellites are identified to be unlikely to participate in the plane, two of which are clearly non-members. Despite their previously reported line-of-sight velocity trend suggestive of a common co-rotating motion, 17 out of 22 potential CASP members are consistent with either orbital direction within both the full range of possible kinematics as well as when limiting orbits to those within the plane. On the other hand, disregarding the 5 off-plane satellites found to be inconsistent with CASP membership enhances the significance of the CASP's line-of-sight velocity trend fivefold. Our results are robust with different mass estimates of the Centaurus A halo, and the adoption of either spherical or triaxial NFW potentials. [ABSTRACT FROM AUTHOR]

Decadal changes in a nearby supernova remnant (SNR) were analyzed using a multiepoch maximum likelihood estimation (MLE) approach. To achieve greater accuracy in capturing the dynamics of SNRs, kinematic features and point-spread function effects were integrated into the MLE framework. Using Cassiopeia A as a representative example, data obtained by the Chandra X-ray Observatory in 2000, 2009, and 2019 were utilized. The proposed multiepoch MLE was qualitatively and quantitatively demonstrated to provide accurate estimates of various motions, including shock waves and faint features, across all regions. To investigate asymmetric structures, such as singular components that deviate from the direction of expansion, the MLE method was extended to combine multiple computational domains and classify kinematic properties using the k -means algorithm. This approach allowed for the mapping of different physical states onto the image, and one classified component was suggested to interact with circumstellar material by comparison with infrared observations from the James Webb Space Telescope. Thus, this technique will help quantify the dynamics of SNRs and discover their unique evolution.

If a planetary nebula (PN) is shown to be a member of a star cluster, we obtain important new constraints on the mass and chemical composition of the PN’s progenitor star, which cannot be determined for PNe in the field. Cluster membership can be tested by requiring the projected separation between the PN and cluster to be within the tidal radius of the cluster, and the objects to have nearly identical radial velocities (RVs) and interstellar extinctions, and nearly identical proper motions (PMs). In an earlier study, we used PMs to confirm that three PNe, which had already passed the other tests, are highly likely to be members of Galactic globular clusters (GCs). For a fourth object, the PN JaFu 1, which lies in the Galactic bulge near the GC Palomar 6 on the sky and has a similar RV, the available PM measurement gave equivocal results. We have now obtained new high-resolution images of the central star of JaFu 1 with the Hubble Space Telescope (HST), which, combined with archival HST frames taken 14 and 16 yr earlier, provide a high-precision PM. Unfortunately, we find that the PM of the central star differs from that of the cluster with high statistical significance. Nevertheless, JaFu 1 is of astrophysical interest because its nucleus appears to be a member of the rare class of “EGB 6-type” central stars, which are associated with compact emission-line knots.

We present a hierarchical Bayesian pipeline, BP3M , that measures positions, parallaxes, and proper motions (PMs) for cross-matched sources between Hubble Space Telescope (HST) images and Gaia—even for sparse fields ( N _* < 10 per image)—expanding from the recent GaiaHub tool. This technique uses Gaia-measured astrometry as priors to predict the locations of sources in HST images, and is therefore able to put the HST images onto a global reference frame without the use of background galaxies/QSOs. Testing our publicly available code in the Fornax and Draco dwarf spheroidal galaxies, we measure PMs that are a median of 8–13 times more precise than Gaia DR3 alone for 20.5 < G < 21 mag. We are able to explore the effect of observation strategies on BP3M astrometry using synthetic data, finding an optimal strategy to improve parallax and position precision at no cost to the PM uncertainty. Using 1619 HST images in the sparse COSMOS field (median nine Gaia sources per HST image), we measure BP3M PMs for 2640 unique sources in the 16 < G < 21.5 mag range, 25% of which have no Gaia PMs; the median BP3M PM uncertainty for 20.25 < G < 20.75 mag sources is 0.44 mas yr ^−1 compared to 1.03 mas yr ^−1 from Gaia, while the median BP3M PM uncertainty for sources without Gaia-measured PMs (20.75 < G < 21.5 mag) is 1.16 mas yr ^−1 . The statistics that underpin the BP3M pipeline are a generalized way of combining position measurements from different images, epochs, and telescopes, which allows information to be shared between surveys and archives to achieve higher astrometric precision than that from each catalog alone.

We have determined the proper motions (PMs) of 12 dwarf galaxies in the Local Group (LG), ranging from the outer Milky Way (MW) halo to the edge of the LG. We used the Hubble Space Telescope (HST) as the first and Gaia as the second epoch using the GaiaHub software. For Leo A and Sag DIG, we also used multi-epoch HST measurements relative to background galaxies. Orbital histories derived using these PMs show that two-thirds of the galaxies in our sample are on first infall with >90% certainty. The observed star formation histories of these first-infall dwarfs are generally consistent with infalling dwarfs in simulations. The remaining four galaxies have crossed the virial radius of either the MW or M31. When we compare their star formation (SF) and orbital histories we find tentative agreement between the inferred pattern of SF with the timing of dynamical events in the orbital histories. For Leo I, SF activity rises as the dwarf crosses the MW’s virial radius, culminating in a burst of SF shortly before pericenter (≈1.7 Gyr ago). The SF then declines after pericenter, but with some smaller bursts before its recent quenching (≈0.3 Gyr ago). This shows that even small dwarfs like Leo I can hold onto gas reservoirs and avoid quenching for several gigayears after falling into their host, which is longer than generally found in simulations. Leo II, NGC 6822, and IC 10 are also qualitatively consistent with this SF pattern in relation to their orbit, but more tentatively due to larger uncertainties.

Omega Centauri ( ω Cen) is the most massive globular cluster of the Milky Way. It is thought to be the nucleus of an accreted dwarf galaxy because of its high mass and its complex stellar populations. To decipher its formation history and study its dynamics, we created the most comprehensive kinematic catalog for its inner region, by analyzing both archival and new Hubble Space Telescope (HST) data. Our catalog contains 1,395,781 proper-motion measurements out to the half-light radius of the cluster ( $\sim 5.0^{\prime} $ ) and down to m _F625W ≈ 25 mag. The typical baseline for our proper-motion measurements is 20 yr, leading to a median 1D proper motion precision of ∼11 μ as yr ^−1 for stars with m _F625W ≈ 18 mag, with even better precision (∼6.6 μ as yr ^−1 ) achieved in the extensively observed centermost ( $r\lt 1.5^{\prime} $ ) region. In addition to our astrometric measurements, we also obtained precise HST photometry in seven filters spanning from the ultraviolet to the near-infrared. This allows detailed color–magnitude diagram studies and separation of the multiple stellar populations of the cluster. In this work, we describe the data reduction used to obtain both the photometric and the proper-motion measurements. We also illustrate the creation and the content of our catalog, which is made publicly available. Finally, we present measurements of the plane-of-sky rotation of ω Cen in the previously unprobed inner few arcminutes and a precise measurement of the inclination, i = 43.°9 ± 1.°3.

The mergers of double neutron star (DNS) systems are believed to drive the majority of short γ -ray bursts (SGRBs), while also serving as production sites of heavy r -process elements. Despite being key to (i) confirming the nature of the extragalactic SGRBs, (ii) addressing the poorly understood r -process enrichment in the ultrafaint dwarf galaxies (UFDGs), and (iii) probing the formation process of DNS systems, the space velocity distribution of DNSs is still poorly constrained, due to the small number of DNSs with well-determined astrometry. In this work, we determine new proper motions and parallaxes of two Galactic DNSs, PSR J0509+3801 and PSR J1930−1852, using the Very Long Baseline Array, and we estimate the transverse velocities v _⊥ of all 11 isolated Galactic DNSs having proper-motion measurements in a consistent manner. Our correlation analysis reveals that the DNS v _⊥ is tentatively correlated with three parameters: spin period, orbital eccentricity, and companion mass. With the preliminary v _⊥ distribution, we obtain the following findings. First, the refined v _⊥ distribution is confirmed to agree with the observed displacements of the localized SGRBs from their host galaxy birth sites. Second, we estimate that around 11% and 25% of DNSs remain gravitationally bound to UFDGs with escape velocities of 15 and 25 km s ^−1 , respectively. Hence, the retained DNSs might indeed be responsible for the r -process enrichment confirmed so far in a few UFDGs. Finally, we discuss how a future ensemble of astrometrically determined DNSs may probe the multimodality of the v _⊥ distribution.

We analyze four epochs of Hubble Space Telescope imaging over 18 yr for the Draco dwarf spheroidal galaxy. We measure precise proper motions for hundreds of stars and combine these with existing line-of-sight (LOS) velocities. This provides the first radially resolved 3D velocity dispersion profiles for any dwarf galaxy. These constrain the intrinsic velocity anisotropy and resolve the mass–anisotropy degeneracy. We solve the Jeans equations in oblate axisymmetric geometry to infer the mass profile. We find the velocity dispersion to be radially anisotropic along the symmetry axis and tangentially anisotropic in the equatorial plane, with a globally averaged value $\bar{{\beta }_{{\rm{B}}}}=-{0.20}_{-0.53}^{+0.28}$ , (where 1 – ${\beta }_{{\rm{B}}}\equiv \langle {v}_{\tan }^{2}\rangle /\langle {v}_{\mathrm{rad}}^{2}\rangle $ in 3D). The logarithmic dark matter (DM) density slope over the observed radial range, Γ _dark , is $-{0.83}_{-0.37}^{+0.32}$ , consistent with the inner cusp predicted in ΛCDM cosmology. As expected given Draco’s low mass and ancient star formation history, it does not appear to have been dissolved by baryonic processes. We rule out cores larger than 487, 717, and 942 pc at 1 σ , 2 σ , and 3 σ confidence, respectively, thus imposing important constraints on the self-interacting DM cross section. Spherical models yield biased estimates for both the velocity anisotropy and the inferred slope. The circular velocity at our outermost data point (900 pc) is ${24.19}_{-2.97}^{+6.31}\,\mathrm{km}\,{{\rm{s}}}^{-1}$ . We infer a dynamical distance of ${75.37}_{-4.00}^{+4.73}$ kpc and show that Draco has a modest LOS rotation, with $\left\langle v/\sigma \right\rangle =0.22\pm 0.09$ . Our results provide a new stringent test of the so-called “cusp–core” problem that can be readily extended to other dwarfs.

A large number of candidate binary stars with apparent acceleration on the sky has emerged from analysis of astrometric data collected by the Hipparcos, Tycho-2, and Gaia space missions. Although the apparent acceleration can serve as a relatively reliable indicator of binarity, it provides scarce information about the orbital and physical parameters of the components. With an emphasis on the search for stellar-mass black holes and neutron stars hidden in binary systems, we start a broader effort to characterize the most promising candidates using follow-up ground-based observations. Accurate quantification of orbital and physical parameters of systems with dim or invisible companions requires combination of Hipparcos, Gaia, and precision spectroscopic measurements. In this paper, we review the necessary steps in this implementation and describe the improved Hipparcos–Gaia sample of long-term astrometric accelerations, which includes correction of sky-correlated systematic errors using the vector spherical decomposition method. As an example, we study one Hipparcos star with a large acceleration, HIP 44842, where the companion is revealed to be a normal main-sequence star.