Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Astrometry”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Keck/NIRC2 L’-band Imaging of Jovian-mass Accreting Protoplanets around PDS 70

We present L’-band imaging of the PDS 70 planetary system with Keck/NIRC2 using the new infrared pyramid wave front sensor. We detected both PDS 70 b and c in our images, as well as the front rim of the circumstellar disk. After subtracting off a model of the disk, we measured the astrometry and photometry of both planets. Placing priors based on the dynamics of the system, we estimated PDS 70 b to have a semimajor axis of 20{sub −4}{sup +3} au and PDS 70 c to have a semimajor axis of 34{sub −6}{sup +12} au (95% credible interval). We fit the spectral energy distribution (SED) of both planets. For PDS 70 b, we were able to place better constraints on the red half of its SED than previous studies and inferred the radius of the photosphere to be 2–3 R {sub Jup}. The SED of PDS 70 c is less well constrained, with a range of total luminosities spanning an order of magnitude. With our inferred radii and luminosities, we used evolutionary models of accreting protoplanets to derive a mass of PDS 70 b between 2 and 4 M {sub Jup} and a mean mass accretion rate between 3 × 10{sup −7} and 8 × 10{sup −7} M {sub Jup}/yr. For PDS 70 c, we computed a mass between 1 and 3 M {sub Jup} and mean mass accretion rate between 1 × 10{sup −7} and 5 × 10{sup −7} M {sub Jup}/yr. The mass accretion rates imply dust accretion timescales short enough to hide strong molecular absorption features in both planets’ SEDs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Nearby Quintuple Systems κ Tucanae and ξ Scorpii

Architecture and parameters of two wide, nearby hierarchical systems containing five solar-type stars each, κ Tuc and ξ Sco, are studied. Using Gaia astrometry and photometry, masses are determined from visual orbits and isochrones, and effective temperatures from spectra or colors. Both systems are ∼2 Gyr old. Their spatial motion corresponds to a young disk but does not match any known kinematic group. Internal proper motions relative to the center of mass and radial velocities show that wide ∼8 kau outer pairs are bound. No correlation between orbit orientations in the inner subsystems is observed. All masses except one are confined to the narrow range from 0.8 to 1.5 solar. Strongly correlated masses and wide orbits can be explained if those systems formed by fragmentation in relative isolation and their components accreted gas from a common source, as expected in a hierarchical collapse. Young moving groups could be formed in similar environments, and many of them contain high-order hierarchies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Survey for New Stars and Brown Dwarfs in the Ophiuchus Star-forming Complex

We have performed a survey for new members of the Ophiuchus cloud complex using high-precision astrometry from the second data release of Gaia, proper motions measured with multi-epoch images from the Spitzer Space Telescope, and color–magnitude diagrams constructed with photometry from various sources. Through spectroscopy of candidates selected with those data, we have identified 155 new young stars. Based on available measurements of kinematics, we classify 102, 47, and 6 of those stars as members of Ophiuchus, Upper Sco, and other populations in Sco–Cen, respectively. We have also assessed the membership of all other stars in the vicinity of Ophiuchus that have spectroscopic evidence of youth from previous studies, arriving at a catalog of 373 adopted members of the cloud complex. For those adopted members, we have compiled mid-infrared photometry from Spitzer and the Wide-field Infrared Survey Explorer and have used mid-infrared colors to identify and classify circumstellar disks. We find that 210 of the members show evidence of disks, including 48 disks that are in advanced stages of evolution. Finally, we have estimated the relative median ages of the populations near the Ophiuchus clouds and the surrounding Upper Sco association using absolute K-band magnitudes (M {sub K}) based on Gaia parallaxes. If we adopt an age 10 Myr for Upper Sco, then the relative values of M {sub K} imply median ages of ∼2 Myr for L1689 and embedded stars in L1688, 3–4 Myr for low-extinction stars near L1688, and ∼6 Myr for the group containing ρ Oph.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Infrared Observations of 2I/Borisov near Perihelion

As the first interstellar comet, 2I/Borisov provides a unique opportunity to study the surface composition of a comet from another stellar system, particularly whether it has water ice. In order to investigate the nature of 2I/Borisov, we conducted infrared observations close to perihelion. The water ice, if present, is expected to be revealed by absorption features at 1.5 and 2 micron. We therefore used FLAMINGOS-2 mounted on the Gemini south telescope, to carry out deep imaging on 2019 November 30 UT and spectroscopy on 2019 December 7 UT. At first glance, our imaging did not reveal an apparent coma or a cometary tail. This is due to the bright sky background and our short exposure times. Nevertheless we were able to put an upper limit of the size of the nucleus, as well as provide high-precision astrometry that can be used to investigate nongravitational acceleration in the future. Our infrared spectra showed a negative slope, contrary to the results by Yang et al. It is not unheard of for a comet to show a negative slope, and to progressively exhibit a spectrum with decreasing slope. Possible causes of the decreasing slope are an increase in water ice and/or decrease in dust size. Given the fact that our observations were carried out close to perihelion, it is likely that both factors contribute to the decreasing, negative slope of the infrared spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

The Most Metal-poor Stars in the Inner Bulge

The bulge is the oldest component of the Milky Way. Since numerous simulations of Milky Way formation have predicted that the oldest stars at a given metallicity are found on tightly bound orbits, the Galaxy’s oldest stars are likely metal-poor stars in the inner bulge with small apocenters (i.e., R {sub apo} ≲ 4 kpc). In the past, stars with these properties have been impossible to find due to extreme reddening and extinction along the line of sight to the inner bulge. We have used the mid-infrared metal-poor star selection of Schlaufman & Casey (2014) on Spitzer/Galactic Legacy Infrared Mid-Plane Survey Extraordinaire data to overcome these problems and target candidate inner bulge metal-poor giants for moderate-resolution spectroscopy with Anglo-Australian Telescope/AAOmega. We used those data to select three confirmed metal-poor giants ([Fe/H] = −3.15, −2.56, −2.03) for follow-up high-resolution Magellan/Magellan Inamori Kyocera Echelle spectroscopy. A comprehensive orbit analysis using Gaia DR2 astrometry and our measured radial velocities confirms that these stars are tightly bound inner bulge stars. We determine the elemental abundances of each star and find high titanium and iron-peak abundances relative to iron in our most metal-poor star. We propose that the distinct abundance signature we detect is a product of nucleosynthesis in the Chandrasekhar-mass thermonuclear supernova of a CO white dwarf accreting from a helium star with a delay time of about 10 Myr. Even though chemical evolution is expected to occur quickly in the bulge, the intense star formation in the core of the nascent Milky Way was apparently able to produce at least one Chandrasekhar-mass thermonuclear supernova progenitor before chemical evolution advanced beyond [Fe/H]∼−3.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

MIRC-X: A Highly Sensitive Six-telescope Interferometric Imager at the CHARA Array

Michigan InfraRed Combiner-eXeter (MIRC-X) is a new highly sensitive six-telescope interferometric imager installed at the CHARA Array that provides an angular resolution equivalent of up to a 330 m diameter baseline telescope in J- and H-band wavelengths ((λ/(2B))∼0.6 mas). We upgraded the original Michigan InfraRed Combiner (MIRC) instrument to improve sensitivity and wavelength coverage in two phases. First, a revolutionary sub-electron noise and fast-frame-rate C-RED ONE camera based on an SAPHIRA detector was installed. Second, a new-generation beam combiner was designed and commissioned to (i) maximize sensitivity, (ii) extend the wavelength coverage to J band, and (iii) enable polarization observations. A low-latency and fast-frame-rate control software enables high-efficiency observations and fringe tracking for the forthcoming instruments of the CHARA Array. Since mid-2017, MIRC-X has been offered to the community and has demonstrated best-case H-band sensitivity down to 8.2 correlated magnitude. MIRC-X uses single-mode fibers to coherently combine the light from six telescopes simultaneously with an image-plane combination scheme and delivers a visibility precision better than 1%, and closure phase precision better than 1°. MIRC-X aims at (i) imaging protoplanetary disks, (ii) detecting exoplanets with precise astrometry, and (iii) imaging stellar surfaces and starspots at an unprecedented angular resolution in the near-infrared. In this paper, we present the instrument design, installation, operation, and on-sky results, and demonstrate the imaging capability of MIRC-X on the binary system ι Peg. The purpose of this paper is to provide a solid reference for studies based on MIRC-X data and to inspire future instruments in optical interferometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

OGLE-2018-BLG-1269Lb: A Jovian Planet with a Bright I = 16 Host

We report the discovery of a planet in the microlensing event OGLE-2018-BLG-1269 with a planet–host mass ratio q ∼ 6 × 10{sup −4}, i.e., 0.6 times smaller than the Jupiter/Sun mass ratio. Combined with the Gaia parallax and proper motion, a strong one-dimensional constraint on the microlens parallax vector allows us to significantly reduce the uncertainties of lens physical parameters. A Bayesian analysis that ignores any information about light from the host yields that the planet is a cold giant (M{sub 2}=0.69{sub −0.22}{sup +0.44} M{sub J}) orbiting a Sun-like star (M{sub 1}=1.13{sub −0.35}{sup +0.72} M{sub ⊙}) at a distance of D{sub L}=2.56{sub −0.62}{sup +0.92} kpc. The projected planet–host separation is a{sub ⊥}=4.61{sub −1.17}{sup +1.70} au. Using Gaia astrometry, we show that the blended light lies ≲12 mas from the host and therefore must be either the host star or a stellar companion to the host. An isochrone analysis favors the former possibility at >99.6%. The host is therefore a subgiant. For host metallicities in the range of 0.0⩽[Fe/H]⩽+0.3, the host and planet masses are then in the range of 1.16⩽M{sub 1}/M{sub ⊙}⩽1.38 and 0.74⩽M{sub 2}/M{sub J}⩽0.89, respectively. Low host metallicities are excluded. The brightness and proximity of the lens make the event a strong candidate for spectroscopic follow-up both to test the microlensing solution and to further characterize the system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Dynamical Mass Estimates of the β Pictoris Planetary System through Gaussian Process Stellar Activity Modeling

Nearly 15 yr of radial velocity (RV) monitoring and direct imaging enabled the detection of two giant planets orbiting the young, nearby star β Pictoris. The δ Scuti pulsations of the star, which overwhelm planetary signals, need to be carefully suppressed. In this work, we independently revisit the analysis of the RV data following a different approach than available in the literature to model the activity of the star. We show that a Gaussian process (GP) with a stochastically driven damped harmonic oscillator kernel can model the δ Scuti pulsations. It provides similar results to parametric models but with a simpler framework, using only three hyperparameters. It also enables us to model poorly sampled RV data that were excluded from previous analyses, hence extending the RV baseline by nearly five years. Altogether, the orbit and mass of both planets can be constrained from RV only, which was not possible with the parametric modeling. To characterize the system more accurately, we also perform a joint fit of all available relative astrometry and RV data. Our orbital solutions for β Pic b favor a low eccentricity of 0.029$_{−0.024}^{+0.061}$ and a relatively short period of 21.1$_{−0.8}^{+2.0}$ yr. The orbit of β Pic c is eccentric with 0.206$_{−0.063}^{+0.074}$ with a period of 3.36 ± 0.03 yr. We find model-independent masses of 11.7 ± 1.4 and 8.5 ± 0.5 M Jup for β Pic b and c, respectively, assuming coplanarity. The mass of β Pic b is consistent with the hottest start evolutionary models, at an age of 25 ± 3 Myr. A direct detection of β Pic c would provide a second calibration measurement in a coeval system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Two-epoch Orbit Estimation for Wide Binaries Resolved in Hipparcos and Gaia

The Hipparcos catalog and its Double and Multiple System Annex (DMSA) lists 4099 components with individual proper motions and coordinates on the epoch 1991.25. Many of these long-period binary stars are also present in Gaia Data Release 2 (DR2). Using the available relative positions and proper motions separated by 24.25 yr, the equations of relative orbital motion can be solved for the two epoch eccentric anomalies, orbital period, and eccentricity. This method employs elimination of the linear Thiele–Innes unknowns and nonlinear optimization of the remaining condition equations. The quality of these solutions is compromised by the insufficient condition and modest precision of the Hipparcos astrometric data, as revealed by Monte Carlo simulations with artificially perturbed data points. The presence of multiple systems and optical pairs can also perturb the results. Limited experiments with artificial data indicate that useful estimates can be obtained with a 25 yr epoch difference for wide binaries with orbital periods up to ~500 yr. The prospects of this method dramatically improve with the proposed next-generation space astrometry missions such as Gaia-NIR and Theia, especially when additional conditions are included from astrometric or spectroscopic measurements. An ancillary catalog of cross-identification and astrometric information for 1295 double-star pairs cross-matched in Gaia DR2 and Hipparcos is also published.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

An Asymmetric Eclipse Seen toward the Pre-main-sequence Binary System V928 Tau

K2 observations of the weak-lined T Tauri binary V928 Tau A and B show the detection of a single, asymmetric eclipse, which may be due to a previously unknown substellar companion eclipsing one component of the binary with an orbital period >66 days. Over an interval of about 9 hr, one component of the binary dims by around 60%, returning to its normal brightness about 5 hr later. From modeling of the eclipse shape, we find evidence that the eclipsing companion may be surrounded by a disk or a vast ring system. The modeled disk has a radius of 0.9923 ± 0.0005 R∗, with an inclination of 56.°78 ± 0.°03, a tilt of 41.°22 ± 0.°05, an impact parameter of -0.2506 ± 0.0002 R∗, and an opacity of 1.00. The occulting disk must also move at a transverse velocity of 6.637 ± 0.002 R∗ day -1 , which, depending on whether it orbits V928 Tau A or B, corresponds to approximately 73.53 or 69.26 km s -1 . A search in ground-based archival data reveals additional dimming events, some of which suggest periodicity, but no unambiguous period associated with the eclipse observed by K2. We present a new epoch of astrometry that is used to further refine the orbit of the binary, presenting a new lower bound of 67 yr, and constraints on the possible orbital periods of the eclipsing companion. The binary is also separated by 18' (~2250 au) from the lower-mass CFHT-BD-Tau 7, which is likely associated with V928 Tau A and B. We also present new high-dispersion optical spectroscopy that we use to characterize the unresolved stellar binary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Statistics of Wide Pre-main-sequence Binaries in the Orion OB1 Association

Statistics of low-mass pre-main-sequence binaries in the Orion OB1 association with separations ranging from 0."6 to 20" (220 to 7400 au at 370 pc) are studied using images from the VISTA Orion mini survey and astrometry from Gaia. The input sample based on the CVSO catalog contains 1137 stars of K and M spectral types (masses between 0.3 and 0.9 M ⊙ ), 1021 of which are considered to be association members. There are 135 physical binary companions to these stars with mass ratios above ~0.13. The average companion fraction is 0.09 ± 0.01 over 1.2 decades in separation, slightly less than, but still consistent with, the field. We found a difference between the Ori OB1a and OB1b groups, the latter being richer in binaries by a factor of 1.6 ± 0.3. No overall dependence of the wide-binary frequency on the observed underlying stellar density is found, although in the Ori OB1a off-cloud population, these binaries seem to avoid dense clusters. The multiplicity rates in Ori OB1 and in sparse regions like Taurus differ significantly, hinting that binaries in the field may originate from a mixture of diverse populations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

ARMADA. I. Triple Companions Detected in B-type Binaries α Del and ν Gem

Ground-based optical long-baseline interferometry has the power to measure the orbits of close binary systems at ∼10 μas precision. This precision makes it possible to detect “wobbles” in the binary motion due to the gravitational pull from additional short-period companions. We started the ARrangement for Micro-Arcsecond Differential Astrometry (ARMADA) survey with the Michigan Infra-Red Combiner (MIRC)/MIRC-X instrument at the Center for High Angular Resoloution Astronomy (CHARA) array for the purpose of detecting giant planets and stellar companions orbiting individual stars in binary systems. We describe our observations for the survey, and introduce the wavelength calibration scheme that delivers precision at the tens of microarcseconds level for <02 binaries. We test our instrument performance on a known triple system, κ Peg, and show that our survey is delivering a factor of 10 better precision than previous similar surveys. We present astrometric detections of tertiary components to two B-type binaries: a 30 day companion to α Del, and a 50 day companion to ν Gem. We also collected radial velocity data for α Del with the Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We are able to measure the orbits and masses of all three components in these systems. We find that the previously published radial velocity orbit for the inner pair of ν Gem is not consistent with our visual orbit. The precision achieved for these orbits suggests that our ARMADA survey will be successful at discovering new compact triple systems to A/B-type binary systems, leading to better statistics of hierarchical system architectures and formation history.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Constraining the Nature of the PDS 70 Protoplanets with VLTI/GRAVITY

We present K-band interferometric observations of the PDS 70 protoplanets along with their host star using VLTI/GRAVITY. We obtained K-band spectra and 100 μas precision astrometry of both PDS 70 b and c in two epochs, as well as spatially resolving the hot inner disk around the star. Rejecting unstable orbits, we found a nonzero eccentricity for PDS 70 b of 0.17 ± 0.06, a near-circular orbit for PDS 70 c, and an orbital configuration that is consistent with the planets migrating into a 2:1 mean motion resonance. Enforcing dynamical stability, we obtained a 95% upper limit on the mass of PDS 70 b of 10 M {sub Jup}, while the mass of PDS 70 c was unconstrained. The GRAVITY K-band spectra rules out pure blackbody models for the photospheres of both planets. Instead, the models with the most support from the data are planetary atmospheres that are dusty, but the nature of the dust is unclear. Any circumplanetary dust around these planets is not well constrained by the planets’ 1–5 μm spectral energy distributions (SEDs) and requires longer wavelength data to probe with SED analysis. However with VLTI/GRAVITY, we made the first observations of a circumplanetary environment with sub-astronomical-unit spatial resolution, placing an upper limit of 0.3 au on the size of a bright disk around PDS 70 b.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Astrometric Membership Tests for the Zinn–Newell–Gibson UV-bright Stars in Galactic Globular Clusters

In 1972, Zinn, Newell, & Gibson (ZNG) published a list of 156 candidate UV-bright stars they had found in 27 Galactic globular clusters (GCs), based on photographs in the U and V bands. UV-bright stars lie above the horizontal branch (HB) and blueward of the asymptotic giant branch (AGB) and red giant branch in the clusters’ color–magnitude diagrams. They are in rapid evolutionary phases—if they are members and not unrelated bright foreground stars. The ZNG list has inspired numerous follow-up studies, aimed at understanding late stages of stellar evolution. However, the ZNG candidates were presented only in finding charts, and celestial coordinates were not given. Using my own collection of CCD frames in u and V, I have identified all of the ZNG objects, and have assembled their coordinates, parallaxes, and proper motions from the recent Gaia Early Data Release 3 (EDR3). Based on the Gaia astrometry, I have determined which objects are probable cluster members (45% of the sample). For the members, using photometry from EDR3, I have assigned the stars to various evolutionary stages, including luminous post-AGB stars, and stars above the HB. I point out several ZNG stars of special interest that have still, to my knowledge, never been studied in detail. This study is an adjunct to a forthcoming survey of the Galactic GCs in the uBVI photometric system, designed for detection of low-gravity stars with large Balmer discontinuities.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Functional Principal Component Analysis of Radio–Optical Reference Frame Tie

The Gaia optical reference frame is intrinsically undefined with respect to global orientation and spin, so it needs to be anchored in the radio-based International Celestial Reference Frame (ICRF) to provide a referenced and quasi-inertial celestial coordinate system. The link between the two fundamental frames is realized through two samples of distant extragalactic sources, mostly active galactic nuclei and quasars, but only the smaller sample of radio-loud ICRF sources with optical counterparts is available to determine the mutual orientation. The robustness of this link can be mathematically formulated in the framework of functional principal component analysis using a set of vector spherical harmonics to represent the differences in celestial positions of the common objects. The weakest eigenvectors are computed, which describe the greatest deficiency of the link. The deficient or poorly determined terms are specific vector fields on the sphere that carry the largest errors of absolute astrometry using Gaia in reference to the ICRF. This analysis provides guidelines for the future development of the ICRF maximizing the accuracy of the link over the entire celestial sphere. A measure of the robustness of a least-squares solution, which can be applied to any linear model fitting problem, is introduced to help discriminate between reference frame tie models of different degrees.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

DELVE-ing into the Jet: A Thin Stellar Stream on a Retrograde Orbit at 30 kpc

We perform a detailed photometric and astrometric analysis of stars in the Jet stream using data from the first data release of the DECam Local Volume Exploration Survey DR1 and Gaia EDR3. We discover that the stream extends over ∼ 29° on the sky (increasing the known length by 18°), which is comparable to the kinematically cold Phoenix, ATLAS, and GD-1 streams. Using blue horizontal branch stars, we resolve a distance gradient along the Jet stream of 0.2 kpc deg$^{−1}$, with distances ranging from D $_{⊙}$ ∼ 27–34 kpc. We use natural splines to simultaneously fit the stream track, width, and intensity to quantitatively characterize density variations in the Jet stream, including a large gap, and identify substructure off the main track of the stream. Furthermore, we report the first measurement of the proper motion of the Jet stream and find that it is well aligned with the stream track, suggesting the stream has likely not been significantly perturbed perpendicular to the line of sight. Finally, we fit the stream with a dynamical model and find that it is on a retrograde orbit, and is well fit by a gravitational potential including the Milky Way and Large Magellanic Cloud. These results indicate the Jet stream is an excellent candidate for future studies with deeper photometry, astrometry, and spectroscopy to study the potential of the Milky Way and probe perturbations from baryonic and dark matter substructure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Chemical Composition of Extreme-velocity Stars* †

Abstract Little is known about the origin of the fastest stars in the Galaxy. Our understanding of the chemical evolution history of the Milky Way and surrounding dwarf galaxies allows us to use the chemical composition of a star to investigate its origin and to say whether it was formed in situ or was accreted. However, the fastest stars, the hypervelocity stars, are young and massive and their chemical composition has not yet been analyzed. Though it is difficult to analyze the chemical composition of a massive young star, we are well versed in the analysis of late-type stars. We have used high-resolution ARCES/3.5 m Apache Point Observatory, MIKE/Magellan spectra to study the chemical details of 15 late-type hypervelocity star candidates. With Gaia EDR3 astrometry and spectroscopically determined radial velocities we found total velocities with a range of 274–520 km s −1 and mean value of 381 km s −1 . Therefore, our sample stars are not fast enough to be classified as hypervelocity stars, and are what is known as extreme-velocity stars. Our sample has a wide iron abundance range of −2.5 ≤ [Fe/H] ≤ −0.9. Their chemistry indicates that at least 50% of them are accreted extragalactic stars, with iron-peak elements consistent with prior enrichment by sub-Chandrasekhar mass Type Ia supernovae. Without indication of binary companions, their chemical abundances and orbital parameters indicate that they are the accelerated tidal debris of disrupted dwarf galaxies.

Reggiani, Henrique (ORCID:0000000165336179)↗
DOE OSTI2022

The Galactic Distribution of Phosphorus: A Survey of 163 Disk and Halo Stars*

Phosphorus (P) is a critical element for life on Earth, yet the cosmic production sites of P are relatively uncertain. To understand how P has evolved in the solar neighborhood, we measured abundances for 163 FGK stars over a range of –1.09 < [Fe/H] < 0.47 using observations from the Habitable-zone Planet Finder instrument on the Hobby–Eberly Telescope. Atmospheric parameters were calculated by fitting a combination of astrometry, photometry, and Fe I line equivalent widths. Phosphorus abundances were measured by matching synthetic spectra to a P I feature at 10529.52 Å. Our [P/Fe] ratios show that chemical evolution models generally underpredict P over the observed metallicity range. Additionally, we find that the [P/Fe] differs by ~0.1 dex between thin disk and thick disk stars that were identified with kinematics. The P abundances were compared with α-elements, iron-peak, odd-Z, and s-process elements, and we found that the evolution of P in the disk most strongly resembles that of the α-elements. We also find that molar P/C and N/C ratios for our sample match the scatter seen from other abundance studies. Finally, we measure a [P/Fe] = 0.09 ± 0.1 ratio in one low-α halo star and probable Gaia–Sausage–Enceladus member, an abundance ratio ~0.3–0.5 dex lower than the other Milky Way disk and halo stars at similar metallicities. Overall, we find that P is likely most significantly produced by massive stars in core-collapse supernovae, based on the largest P abundance survey to date.

79 ASTRONOMY AND ASTROPHYSICS↗