Search

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]

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.

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]

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]

Calvera (1RXS J141256.0+792204) is a pulsar of characteristic age 285 kyr at a high Galactic latitude of b = +37°, detected only in soft thermal X-rays. We measure a new and precise proper motion for Calvera using Chandra High Resolution Camera observations obtained 10 yr apart. We also derive a new phase-connected ephemeris using 6 yr of NICER data, including the astrometric position and proper motion as fixed parameters in the timing solution. Calvera is located near the center of a faint, circular radio ring that was recently discovered by LOFAR and confirmed as a supernova remnant (SNR) by the detection of γ -ray emission with Fermi Large Area Telescope. The proper motion of 78.5 ± 2.9 mas yr ^−1 at position angle 241.°3 ± 2.°2 (in Galactic coordinates) points away from the center of the ring, a result which differs markedly from a previous low-significance measurement, and greatly simplifies the interpretation of the SNR/pulsar association. It argues that the supernova indeed birthed Calvera

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.

We verify candidate hypervelocity red clump stars located in the Galactic bulge that were selected based on the VVV and the Gaia DR2 data by Luna et al. To do so, we analyze data from the OGLE-IV survey: difference images and astrometric time series. We have data for 30 stars out of 34 hypervelocity candidates. We confirmed the high proper motion of only one of these stars and find out that it is a nearby one, hence, not a hypervelocity star. To sum up, we do not confirm the candidate stars as hypervelocity ones. Hence, we disprove the production rate of hypervelocity red clump stars by the central supermassive black hole provided by Luna et al.

The Cosmological Advanced Survey Telescope for Optical and ultraviolet Research (CASTOR) is a proposed Canadian-led 1 m class space telescope that will carry out ultraviolet and blue optical wide-field imaging, spectroscopy, and photometry. CASTOR will provide an essential bridge in the post-Hubble era, preventing a protracted UV-optical gap in space astronomy and enabling an enormous range of discovery opportunities from the solar system to the nature of the cosmos, in conjunction with the other great wide-field observatories of the next decade (e.g., Euclid, Roman, Vera Rubin). FORECASTOR (Finding Optics Requirements and Exposure times for CASTOR) will supply a coordinated suite of mission-planning tools that will serve as the one-stop shop for proposal preparation, data reduction, and analysis for the CASTOR mission. We present the first of these tools: a pixel-based, user-friendly, extensible, multi-mission exposure time calculator built in Python, including a modern browser-based graphical user interface that updates in real time. We then provide several illustrative examples of FORECASTOR’s use that advance the design of planned legacy surveys for the CASTOR mission: a search for the most massive white dwarfs in the Magellanic Clouds, a study of the frequency of flaring activity in M stars and their distribution and impacts on habitability of exoplanets, mapping the proper motions of faint stars in the Milky Way, wide and deep galaxy surveys, and time-domain studies of active galactic nuclei.

The (sub)millimeter wavelengths (86–690 GHz) very long baseline interferometry will provide ∼5–40 μ as angular resolution, ∼10 mJy baseline sensitivity, and ∼1 μ as yr ^−1 proper-motion precision, which can directly detect supermassive black hole binary (SMBHB) systems by imaging the two visible sources and tracking their relative motions. Such a way exhibits an advantage compared to indirect detect methods of observing periodic signals in motion and light curves, which are difficult to confirm from competing models. Moreover, tracking relative motion at (sub)millimeter wavelengths is more reliable, as there is a negligible offset between the emission region and the black hole center. In this way, it is unnecessary to correct the black hole location by a prior of jet morphology as it would be required at longer wavelengths. We extend the formalism developed in D’Orazio & Loeb (2018) to link the observations with the orbital evolution of SMBHBs from the ≲10 kpc dynamical friction stages to the ≲0.01 pc gravitational radiation stages, and estimate the detectable numbers of SMBHBs. By assuming 5% of active galactic nuclei holding SMBHBs, we find that the number of detectable SMBHBs with redshift z ≤ 0.5 and mass M ≤ 10 ^11 M _⊙ is about 20. Such a detection relies heavily on proper-motion precision and sensitivity. Furthermore, we propose that the simultaneous multifrequency technique plays a key role in meeting the observational requirements.

Among Type Ia supernova remnants (SNRs), Tycho’s SNR has been considered as a typical object from the viewpoints of its spectroscopic, morphological, and environmental properties. A recent reanalysis of Chandra data showed that its forward shock is experiencing a substantial deceleration since around 2007, which suggests recent shock interactions with a dense medium as a consequence of a cavity-wall environment inside a molecular cloud. Such a nonuniform environment can be linked back to the nature and activities of its progenitor. In this study, we perform hydrodynamic simulations to characterize Tycho’s cavity-wall environment using the latest multiepoch proper motion measurements of the forward shock. A range of parameters for the environment is explored in the hydrodynamic models to fit with the observational data for each azimuthal region. Our results show that a wind-like cavity with ρ ( r ) ∝ r ^−2 reconciles with the latest data better than a uniform medium with a constant density. In addition, our best-fit model favors an anisotropic wind with an azimuthally varying wind parameter. The overall result indicates a mass-loss rate which is unusually high for the conventional single-degenerate explosion scenario.

We report the results of careful astrometric measurements of the cannonball pulsar J0002+6216 carried out over 3 yr using the High Sensitivity Array. We significantly refine the proper motion to μ = 35.3 ± 0.6 mas yr ^−1 and place new constraints on the distance, with the overall effect of lowering the velocity and increasing the inferred age to 47.60 ± 0.80 kyr. Although the pulsar is brought more in line with the standard natal kick distribution, this new velocity has implications for the morphology of the pulsar wind nebula that surrounds it, the density of the interstellar medium through which it travels, and the age of the supernova remnant (CTB 1) from which it originates.

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.

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.

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.

We develop and disseminate effective point-spread functions and geometric-distortion solutions for high-precision astrometry and photometry with the JWST NIRISS instrument. We correct field dependencies and detector effects, and assess the quality and the temporal stability of the calibrations. As a scientific application and validation, we study the proper motion (PM) kinematics of stars in the JWST calibration field near the Large Magellanic Cloud (LMC) center, comparing to a first-epoch Hubble Space Telescope (HST) archival catalog with a 16 yr baseline. For stars with G ∼ 20, the median PM uncertainty is ∼13 μ as yr ^−1 (3.1 km s ^−1 ), better than Gaia DR3 typically achieves for its very best-measured stars. We kinematically detect the known star cluster OGLE-CL LMC 407, measure its absolute PM for the first time, and show how this differs from other LMC populations. The inferred cluster dispersion sets an upper limit of 24 μ as yr ^−1 (5.6 km s ^−1 ) on systematic uncertainties. Red-giant-branch stars have a velocity dispersion of 33.8 ± 0.6 km s ^−1 , while younger blue populations have a narrower velocity distribution, but with a significant kinematical substructure. We discuss how this relates to the larger velocity dispersions inferred from Gaia DR3. These results establish JWST as capable of state-of-the-art astrometry, building on the extensive legacy of HST. This is the first paper in a series by our JWST Telescope Scientist Team, in which we will use Guaranteed Time Observations to study the PM kinematics of various stellar systems in the Local Group.

Our understanding of the kinematic properties of multiple stellar populations (mPOPs) in Galactic globular clusters (GCs) is still limited compared to what we know about their chemical and photometric characteristics. Such limitation arises from the lack of a comprehensive observational investigation of this topic. Here we present the first homogeneous kinematic analysis of mPOPs in 56 GCs based on high-precision proper motions computed with Hubble Space Telescope data. We focused on red-giant-branch stars, for which the mPOP tagging is clearer, and measured the velocity dispersion of stars belonging to first (1G) and second generations (2G). We find that 1G stars are generally kinematically isotropic even at the half-light radius, whereas 2G stars are isotropic at the center and become radially anisotropic before the half-light radius. The radial anisotropy is induced by a lower tangential velocity dispersion of 2G stars with respect to the 1G population, while the radial component of the motion is comparable. We also show possible evidence that the kinematic properties of mPOPs are affected by the Galactic tidal field, corroborating previous observational and theoretical results suggesting a relation between the strength of the external tidal field and some properties of mPOPs. Although limited to the GCs’ central regions, our analysis leads to new insights into the mPOP phenomenon, and provides the motivation for future observational studies of the internal kinematics of mPOPs.

Sagittarius A* (Sgr A*) is a strong, compact radio source believed to be powered by a supermassive black hole at the galactic center. Extinction by dust and gas in the galactic plane prevents observing it optically, but its position and proper motion have previously been estimated using radio interferometry. We present new VLBI absolute astrometry measurements of its precise position and proper motion in the frame of the third realization of the International Celestial Reference Frame (ICRF3). The observations used were made at 52 epochs on the VLBA at K band (24 GHz) between 2006 June and 2022 August. We find the proper motion of Sgr A* to be −3.128 ± 0.042 mas yr ^−1 in R.A. and −5.584 ± 0.075 mas yr ^−1 in decl., or 6.400 ± 0.073 mas yr ^−1 at a position angle of 209.°26 ± 0.°51. We also find its J2000 ICRF3 coordinates at the 2015.0 proper motion epoch to be 17 ^h 45 ^m 40.ˢ034047 ± 0.ˢ000018, −29°00′28.″21601 ± 0.″00044. In galactic coordinates, Sgr A* shows proper motion of −6.396 ± 0.071 mas yr ^−1 in galactic longitude and −0.239 ± 0.045 mas yr ^−1 in galactic latitude, indicating solar motion of 248.0 ± 2.8 km s ^−1 in the galactic plane and 9.3 ± 1.9 km s ^−1 toward the north galactic pole.

This paper presents improved constraints on the low-mass stellar initial mass function (IMF) of the Boötes I (Boo I) ultrafaint dwarf galaxy, based on our analysis of recent deep imaging from the Hubble Space Telescope. The identification of candidate stellar members of Boo I in the photometric catalog produced from these data was achieved using a Bayesian approach, informed by complementary archival imaging data for the Hubble Ultra Deep Field. Additionally, the existence of earlier-epoch data for the fields in Boo I allowed us to derive proper motions for a subset of the sources and thus identify and remove likely Milky Way stars. We were also able to determine the absolute proper motion of Boo I, and our result is in agreement with, but completely independent of, the measurement(s) by Gaia. The best-fitting parameter values of three different forms of the low-mass IMF were then obtained through forward modeling of the color–magnitude data for likely Boo I member stars within an approximate Bayesian computation Markov Chain Monte Carlo algorithm. The best-fitting single power-law IMF slope is $\alpha =-{1.95}_{-0.28}^{+0.32}$ , while the best-fitting broken power-law slopes are ${\alpha }_{1}=-{1.67}_{-0.57}^{+0.48}$ and ${\alpha }_{2}=-{2.57}_{-1.04}^{+0.93}$ . The best-fitting lognormal characteristic mass and width parameters are ${M}_{c}={0.17}_{-0.11}^{+0.05}{{ \mathcal M }}_{\odot }$ and $\sigma ={0.49}_{-0.20}^{+0.13}$ . These broken power-law and lognormal IMF parameters for Boo I are consistent with published results for the stars within the Milky Way, and thus it is plausible that Boötes I and the Milky Way are populated by the same stellar IMF.

We measure the proper motion of the pulsar PSR J1745-2900 relative to the Galactic center massive black hole, Sgr A*, using the Very Long Baseline Array (VLBA). The pulsar has a transverse velocity of 236 ± 11 km s-1 at position angle 22 ± 2 deg east of north at a projected separation of 0.097 pc from Sgr A*. Given the unknown radial velocity, this transverse velocity measurement does not conclusively prove that the pulsar is bound to Sgr A*; however, the probability of chance alignment is very small. We do show that the velocity and position are consistent with a bound orbit originating in the clockwise disk of massive stars orbiting Sgr A* and a natal velocity kick of 500 km s-1. An origin among the isotropic stellar cluster is possible but less probable. If the pulsar remains radio-bright, multiyear astrometry of PSR J1745-2900 can detect its acceleration and determine the full three-dimensional orbit. We also demonstrate that PSR J1745-2900 exhibits the same angular broadening as Sgr A* over a wavelength range of 3.6 cm to 0.7 cm, further confirming that the two sources share the same interstellar scattering properties. Finally, we place the first limits on the presence of a wavelength-dependent shift in the position of Sgr A*, i.e., the core shift, one of the expected properties of optically thick jet emission. Our results for PSR J1745-2900 support the hypothesis that Galactic center pulsars will originate from the stellar disk and deepen the mystery regarding the small number of detected Galactic center pulsars. [ABSTRACT FROM AUTHOR]

We present the BANYAN All-Sky Survey (BASS) catalog, consisting of 228 new late-type (M4-L6) candidate members of nearby young moving groups (YMGs) with an expected false-positive rate of ~13%. This sample includes 79 new candidate young brown dwarfs and 22 planetary-mass objects. These candidates were identified through the first systematic all-sky survey for late-type low-mass stars and brown dwarfs in YMGs. We cross-matched the Two Micron All Sky Survey and AllWISE catalogs outside of the galactic plane to build a sample of 98,970 potential ≥M5 dwarfs in the solar neighborhood and calculated their proper motions with typical precisions of 5-15 mas yr-1. We selected highly probable candidate members of several YMGs from this sample using the Bayesian Analysis for Nearby Young AssociatioNs II tool (BANYAN II). We used the most probable statistical distances inferred from BANYAN II to estimate the spectral type and mass of these candidate YMG members. We used this unique sample to show tentative signs of mass segregation in the AB Doradus moving group and the Tucana-Horologium and Columba associations. The BASS sample has already been successful in identifying several new young brown dwarfs in earlier publications, and will be of great interest in studying the initial mass function of YMGs and for the search of exoplanets by direct imaging; the input sample of potential close-by ≥M5 dwarfs will be useful to study the kinematics of low-mass stars and brown dwarfs and search for new proper motion pairs. [ABSTRACT FROM AUTHOR]

We identify proper motion objects in the Hubble Ultra Deep Field (UDF) using the optical data from the original UDF program in 2004 and the near-infrared data from the 128 orbit UDF 2012 campaign. There are 12 sources brighter than I = 27 mag that display >3σ significant proper motions. We do not find any proper motion objects fainter than this magnitude limit. Combining optical and near-infrared photometry, we model the spectral energy distribution of each point-source using stellar templates and state-of-the-art white dwarf models. For I ⩽ 27 mag, we identify 23 stars with K0-M6 spectral types and two faint blue objects that are clearly old, thick disk white dwarfs. We measure a thick disk white dwarf space density of 0.1-1.7 × 10-3 pc-3 from these two objects. There are no halo white dwarfs in the UDF down to I = 27 mag. Combining the Hubble Deep Field North, South, and the UDF data, we do not see any evidence for dark matter in the form of faint halo white dwarfs, and the observed population of white dwarfs can be explained with the standard Galactic models. [ABSTRACT FROM AUTHOR]

Using multi-epoch observations of the Stripe 82 region from the Sloan Digital Sky Survey (SDSS), we measure precise statistical proper motions of the stars in the Sagittarius (Sgr) stellar stream. The multi-band photometry and SDSS radial velocities allow us to efficiently select Sgr members and thus enhance the proper-motion precision to ∼0.1mas yr-1.Wemeasure separately the proper motion of a photometrically selected sample of themain-sequence turn-off stars, as well as spectroscopically selected Sgr giants. The data allow us to determine the proper motion separately for the two Sgr streams in the south found in Koposov et al. Together with the precise velocities from SDSS, our proper motions provide exquisite constraints of the three-dimensional motions of the stars in the Sgr streams [ABSTRACT FROM AUTHOR]

Based on long baseline (5-7 years) multi-epoch HST/ACS photometry, used previously to measure the proper motion of M31, we present the proper motions (PMs) of 13 main-sequence Milky Way halo stars. The sample lies at an average distance of r ≃ 24 kpc from the Galactic center, with a root-mean-square spread of 6 kpc. At this distance, the median PM accuracy is 5 km s-1. We devise a maximum likelihood routine to determine the tangential velocity ellipsoid of the stellar halo. The velocity second moments in the directions of the Galactic (l, b) system are (v2 1 )1/2 = 123+29 -23 km s-1, and (v2 b )1/2 = 83+24 -16 km s-1. We combine these results with the known line-of-sight second moment, (v2 los )1/2 = 105 ± 5 km s-1, at this (r) to study the velocity anisotropy of the halo. We find approximate isotropy between the radial and tangential velocity distributions, with anisotropy parameter β = 0.0+0.2 -0.4. Our results suggest that the stellar halo velocity anisotropy out to r ∼ 30 kpc is less radially biased than solar neighborhood measurements. This is opposite to what is expected from violent relaxation, and may indicate the presence of a shell-type structure at r ∼ 24 kpc. With additional multi-epoch HST data, the method presented here has the ability to measure the transverse kinematics of the halo for more stars, and to larger distances. This can yield new improved constraints on the stellar halo formation mechanism, and the mass of the Milky Way. [ABSTRACT FROM AUTHOR]

We discover a late-type L dwarf, WISEP J060738.65+242953.4 (W0607+2429), by comparing the Wide-field Infrared Survey Explorer (WISE) preliminary data release to the Two Micron All Sky Survey (2MASS) in search of high proper motion objects (≳0".3 yr-1). W0607+2429 was found to have a proper motion of 0.57 ± 0".02 yr-1. Based on colors and color-color diagrams using 2MASS and Sloan Digital Sky Survey photometry, we estimate the spectral type (optical) to be L8 within a spectral sub-type. Based on the spectral type estimated we find W0607+2429 to have a distance of Due to image rights restrictions, multiple line equation(s) cannot be graphically displayed pc, making it one of only four very late L dwarfs within 10 pc, and the third closest L dwarf overall. This close L/T transition dwarf will play a pivotal role in resolving outstanding issues of condensate clouds of low-temperature atmospheres. [ABSTRACT FROM AUTHOR]

The Pan-STARRS1 survey has obtained multi-epoch imaging in five bands (Pan-STARRSI gPI, rPI, iPI, zPI, and yPI) on 12 "Medium-Deep fields," each of which spans a 3°73 circle. For the period between 2009 April and 2011 April these fields were observed 50-200 times. Using a reduced proper motion diagram, we have extracted a list of 47 white dwarf (WD) candidates whose Pan-STARRS 1 astrometry indicates a non-zero proper motion at the 6σ level, with a typical 1σ proper motion uncertainty of 10 mas yr-1. We also used astrometry from the Sloan Digital Sky Survey (when available) and USNO-B to assess our proper motion fits. None of the WD candidates exhibits evidence of statistically significant parallaxes, with a typical 1σ uncertainty of 8 mas. Twelve of these candidates are known WDs, including the high proper motion (1".7 yr-1) WD LHS 291. We confirm seven more objects as WDs through optical spectroscopy. Based on the Pan-STARRS1 colors, ten of the stars are likely to be cool WDs with 4170 K < Teff < 5000 K and cooling ages ≈ 8 Gyr. Based on their kinematics, these objects are likely thick disk WDs, but a subset may belong to the halo. Our current sample represents only a small fraction of the Pan-STARRS1 data. With continued coverage from the Medium-Deep Field Survey and the 3π survey, Pan-STARRS1 should find many more high proper motion WDs that are part of the old thick disk and halo. [ABSTRACT FROM AUTHOR]

We report astrometric results of phase-referencing very long baseline interferometry observations of 43 GHz SiO maser emission toward the red hypergiant VY Canis Majoris (VY CMa) using the Very Long Baseline Array (VLBA). We measured a trigonometric parallax of 0.83 ± 0.08 mas, corresponding to a distance of Due to image rights restrictions, multiple line equation(s) cannot be graphically displayed kpc. Compared to previous studies, the spatial distribution of SiO masers has changed dramatically, while its total extent remains similar. The internal motions of the maser spots are up to 1.4 mas yr-1, corresponding to 8 km s-1, and show a tendency for expansion. After modeling the expansion of maser spots, we derived an absolute proper motion for the central star of µx, = -2.8±0.2 and µy = 2.6 ± 0.2 mas yr-1 eastward and northward, respectively. Based on the maser distribution from the VLBA observations, and the relative position between the radio photosphere and the SiO maser emission at 43 GHz from the complementary Very Large Array observations, we estimate the absolute position of VY CMa at mean epoch 2006.53 to be &alphaJ2000 = 07h22m58s.3259 ± 0s.0007, δJ2000 = -25°46'03".063 ± 0".010. The position and proper motion of VY CMa from the VLBA observations differ significantly with values measured by the Hipparcos satellite. These discrepancies are most likely associated with inhomogeneities and dust scattering the optical light in the circumstellar envelope. The absolute proper motion measured with VLBA suggests that VY CMa may be drifting out of the giant molecular cloud to the east of it. [ABSTRACT FROM AUTHOR]

We present high-precision timing of five millisecond pulsars (MSPs) carried out for more than seven years; four pulsars are in binary systems and one is isolated. We are able to measure the pulsars' proper motions and derive an estimate for their space velocities. The measured two-dimensional velocities are in the range 70-210 km s-1, consistent with those measured for other MSPs. We also use all the available proper motion information for isolated and binary MSPs to update the known velocity distribution for these populations. As found by earlier works, we find that the velocity distribution of binary and isolated MSPs are indistinguishable with the current data. Four of the pulsars in our observing program are highly recycled with low-mass white dwarf companions and we are able to derive accurate binary parameters for these systems. For three of these binary systems, we are able to place initial constraints on the pulsar masses with best-fit values in the range 1.0-1.6 M⊙. The implications of the results presented here to our understanding of binary pulsar evolution are discussed. The updated parameters for the binary systems studied here, together with recently discovered similar systems, allowed us to update previous limits on the violation of the strong equivalence principle through the parameter |Δ| to 4.6 x 10-3 (95% confidence) and the violation of Lorentz invariance/momentum conservation through the parameter |α̂ to 5.5 x 10-20 (95% confidence). [ABSTRACT FROM AUTHOR]

We use 14 year baseline images obtained with the Wide Field and Planetary Camera 2 on board the Hubble Space Telescope (HST) to derive a proper motion for one of the Milky Way's most distant dwarf spheroidal companions, Leo II, relative to an extragalactic background reference frame. Astrometric measurements are performed in the effective point-spread function formalism using our own developed code. An astrometric reference grid is defined using 3224 stars that are members of Leo II and brighter than a magnitude of 25 in the F814W band. We identify 17 compact extragalactic sources, for which we measure a systemic proper motion relative to this stellar reference grid. We derive a proper motion [μα,μδ] = [+104±113,-33 ± 151] μas yr-1 for Leo II in the heliocentric reference frame. Though marginally detected, the proper motion yields constraints on the orbit of Leo II. Given a distance of d ≃ 230 kpc and a heliocentric radial velocity vr = +79 km s-1, and after subtraction of the solar motion, our measurement indicates a total orbital motion vG = 266.1 ± 128.7 km s-1 in the Galactocentric reference frame, with a radial component vrG = 21.5 ± 4.3 km s-1 and tangential component vtG = 265.2 ± 129.4 km s-1. The small radial component indicates that Leo II either has a low-eccentricity orbit or is currently close to perigalacticon or apogalacticon distance. We see evidence for systematic errors in the astrometry of the extragalactic sources which, while close to being point sources, are slightly resolved in the HST images. We argue that more extensive observations at later epochs will be necessary to better constrain the proper motion of Leo II. We provide a detailed catalog of the stellar and extragalactic sources identified in the HST data which should provide a solid early-epoch reference for future astrometric measurements. [ABSTRACT FROM AUTHOR]

We report the discovery of a substellar companion to the nearby solar-type star HD 46588 (F7V, 17.9 pc, τ ~ 3 Gyr). HD 46588 B was found through a survey for common proper motion companions to nearby stars using data from the Wide-field Infrared Survey Explorer (WISE) and the Two Micron All Sky Survey (2MASS). It has an angular separation of 79".2 from its primary, which corresponds to a projected physical separation of 1420 AU. We have measured a spectral type of L9 for this object based on near-infrared spectroscopy performed with TripleSpec at Palomar Observatory. We estimate a mass of Due to image rights restrictions, multiple line equation(s) cannot be graphically displayed. Mʘ from a comparison of its luminosity to the values predicted by theoretical evolutionary models for the age of the primary. Because of its companionship to a well-studied star, HD 46588 B is one of the few known brown dwarfs at the L/T transition for which both age and distance estimates are available. Thus, it offers new constraints on the properties of brown dwarfs during this brief evolutionary phase. The discovery of HD 46588 B also illustrates the value of the WISE for identifying brown dwarfs in the solar neighborhood via their proper motions. [ABSTRACT FROM AUTHOR]

With the advent of the Event Horizon Telescope (EHT), a millimeter/submillimeter very long baseline interferometer (VLBI), it has become possible to image a handful of black holes with sub-horizon resolutions. However, these images do not translate into microarcsecond absolute positions due to the lack of absolute phase information when an external phase reference is not used. Due to the short atmospheric coherence time at these wavelengths, nodding between the source and phase reference is impractical. However, here we suggest an alternative scheme which makes use of the fact that many of the VLBI stations within the EHT are arrays in their own right.With this we show that it should be possible to absolutely position the supermassive black holes at the centers of the Milky Way (Sgr A*) and M87 relative to nearby objects with precisions of roughly 1μas. This is sufficient to detect the perturbations to Sgr A*'s position resulting from interactions with the stars and stellar-mass black holes in the Galactic cusp on year timescales, and severely constrain the astrophysically relevant parameter space for an orbiting intermediate-mass black hole, implicated in some mechanisms for producing the young massive stars in the Galactic center. For M87, it allows the registering of millimeter images, in which the black hole may be identified by its silhouette against nearby emission, and existing larger-scale radio images, eliminating present ambiguities in the nature of the radio core and inclination, opening angle, and source of the radio jet. [ABSTRACT FROM AUTHOR]

The Pan-STARRS1 survey has obtained multi-epoch imaging in five bands (Pan-STARRS1 gP1, rP1, iP1, zP1, and yP1) on twelve “Medium Deep Fields”, each of which spans a 3.3 degree circle. For the period between Apr 2009 and Apr 2011 these fields were observed 50–200 times. Using a reduced proper motion diagram, we have extracted a list of 47 white dwarf (WD) candidates whose Pan-STARRS1 astrometry indicates a non-zero proper motion at the 6σ level, with a typical 1σ proper motion uncertainty of 10 mas/yr. We also used astrometry from SDSS (when available) and USNO-B to assess our proper motion fits. None of the WD candidates exhibits evidence of statistically significant parallaxes, with a typical 1σ uncertainty of 8 mas. Twelve of these candidates are known WDs, including the high proper motion (1.7′′ yr−1) WD LHS 291. We confirm three more objects as WDs through optical spectroscopy. Based on the Pan-STARRS1 colors, ten of the stars are likely to be cool WDs with 4170 K < Teff < 5000 K and cooling ages <9 Gyr. We classify these objects as likely thick disk WDs based on their kinematics. Our current sample represents only a small fraction of the Pan-STARRS1 data. With continued coverage from the Medium Deep Field Survey and the 3π survey, Pan-STARRS1 should find many more high proper motion WDs that are part of the old thick disk and halo., Astronomy