Radial Velocities Of Stars Research Articles

Measuring precise radial velocities on individual spectral lines

Context. Radial velocities (RVs) of stars contain both the Doppler reflex motion of potential planetary companions and the drowning and sometimes imitating effect of stellar activity. To separate the two, previous efforts have sought proxies that only trace the activity signals, yet the sub-meter-per-second floor required for the detection of Earth-like planets remains difficult to break. Aims. In this work, we analyze a sample of 12 G- to early M-type stars in order to investigate the feasibility of detecting a differential effect of stellar activity with photospheric depth, as traced by the spectral line-forming temperature, for observations with different sampling and noise levels. Methods. We computed the average line formation temperature for each point in the observed wavelength grids using the spectral synthesis code PySME. The final line selection was curated to exclude blended and poorly synthesized lines. We thereafter computed the convective blueshift (CB) of the line cores of our master spectra (composed of the stacked individual spectra of each star). Finally, we extracted RV time series for certain intervals of formation temperature using a template-matching approach. Results. We find the CB to follow a linear relation with the formation temperature of the line cores, and the CB slope to be steeper with increasing effective temperature. For the RV time series derived for different intervals of formation temperature, we find the RVs of line parts formed at higher temperatures, close to the spectral continuum, to be generally correlated with the S index, and the RVs of line parts formed at cooler temperatures, close to the spectral line cores, to be generally anti-correlated, especially for stars with low noise levels and significant variations over their magnetic cycles. Conclusions. RVs of line parts formed in the coolest 25% of the line-forming temperature range appear to be a strong tracer of stellar activity over the magnetic cycle for several stars. By detrending the total RV time series with a multi-linear combination of residuals of RVs measured at different temperature ranges and the S index, the RV scatter can be decreased to a greater extent than with the S index alone.

Determination of the Parameters of a Nonlinear Kinematic Galactic Rotation Model Based on the Proper Motions and Radial Velocities of Stars from the Gaia DR3 Catalogue

Determination of the Parameters of a Nonlinear Kinematic Galactic Rotation Model Based on the Proper Motions and Radial Velocities of Stars from the Gaia DR3 Catalogue

GaiaData Release 3

Context.The secondGaiadata release, DR2, contained radial velocities of stars with effective temperatures up toTeff = 6900 K. The third data release,GaiaDR3, extends this up toTeff = 14 500 K.Aims.We derive the radial velocities for hot stars (i.e., in theTeff = 6900 − 14 500 K range) from data obtained with the Radial Velocity Spectrometer (RVS) on boardGaia.Methods.The radial velocities were determined by the standard technique of measuring the Doppler shift of a template spectrum that was compared to the observed spectrum. The RVS wavelength range is very limited. The proximity to and systematic blueward offset of the calcium infrared triplet to the hydrogen Paschen lines in hot stars can result in a systematic offset in radial velocity. For the hot stars, we developed a specific code to improve the selection of the template spectrum, thereby avoiding this systematic offset.Results.With the improved code, and with the correction we propose to the DR3 archive radial velocities, we obtain values that agree with reference values to within 3 km s−1(in median). Because of the required S/N for applying the improved code, the hot star radial velocities in DR3 are mostly limited to stars with a magnitude in the RVS wavelength band ≤12 mag.

Kinematics of stellar substructures in the small magellanic cloud

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

The Variability and Radial Velocity of Planetary Nebula Central Stars

The extremely accurate estimates of stellar variability and radial velocity in the Gaia Data Release 3 (Gaia DR3) have enabled us to examine the close binarity and radial velocity (RV) of central stars (CSs) of planetary nebulae (PNe). This study is twofold: (1) searching for new close binary CS candidates to better understand how binarity affects the formation and evolution of PNe; and (2) extending the sample size of known RVs of PNe in order to understand their kinematics and the dynamics of the Milky Way. As a target sample, we used all true, possible, and likely PNe available in the literature. Then, we looked for their matched Gaia DR3 sources that provide measurements of variability and RV. As a result, we detected the first large collection of trustworthy photometric variability of 26 symbiotic stars and 82 CSs. In this CS group, there are 24 sources already classified as true close binary CSs in the literature. Hence, we discovered 58 new close binary CS candidates. This close binary (CB) sample represents more than half of what is currently available in the literature. In addition, we identified the radial velocities for 51 PNe. To our knowledge, 24 of these were measured for the first time. The RV measurements predicted by Gaia, based on the Doppler shift of the CS absorption lines, and those derived from nebular emission lines, show satisfactory agreement except for a few extremely high-velocity PNe.

The Solar Neighborhood L: Spectroscopic Discovery of K Dwarfs Younger Than 1 Gyr and New Binaries within 30 pc

As part of a comprehensive effort to characterize the nearest stars, the CHIRON echelle spectrograph on the CTIO/SMARTS 1.5 m telescope is being used to acquire high-resolution (R = 80,000) spectra of K dwarfs within 50 pc. This paper provides spectral details about 35 K dwarfs from five benchmark sets with estimated ages spanning 20 Myr–5.7 Gyr. Four spectral age and activity indicators are tested, three of which aligned with the estimated ages of the benchmark groups—the Na i doublet (5889.95 and 5895.92 Å), the Hα line (6562.8 Å), and the Li i resonance line (6707.8 Å). The benchmark stars are then used to evaluate seven field K dwarfs exhibiting variable radial velocities for which initial CHIRON data did not show obvious companions. Two of these stars are estimated to be younger than 700 Myr, while one exhibits stellar activity unusual for older K-dwarf field stars and is possibly young. The four remaining stars turn out to be spectroscopic binaries, two of which are being reported here for the first time with orbital periods found using CHIRON data. Spectral analysis of the combined sample of 42 benchmark and variable radial velocity stars indicates temperatures ranging from 3900 to 5300 K and metallicities from −0.4 < [Fe/H] < +0.2. We also determine for main-sequence K dwarfs. Ultimately, this study will target several thousand of the nearest K dwarfs and provide results that will serve present and future studies of stellar astrophysics and exoplanet habitability.

APOGEE-2S Mg–Al anti-correlation of the metal-poor globular cluster NGC 2298

We present detailed elemental abundances and radial velocities of stars in the metal-poor globular cluster (GC) NGC 2298, based on near-infrared high-resolution (R ∼ 22 500) spectra of 12 members obtained during the second phase of the Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) at Las Campanas Observatory as part of the seventeenth Data Release (DR 17) of the Sloan Digital Sky Survey IV (SDSS-IV). We employed the Brussels Automatic Code for Characterizing High accuracy Spectra (BACCHUS) software to investigate abundances for a variety of species including α elements (Mg, Si, and Ca), the odd-Z element Al, and iron-peak elements (Fe and Ni) located in the innermost regions of NGC 2298. We find a mean and median metallicity [Fe/H] = −1.76 and −1.75, respectively, with a star-to-star spread of 0.14 dex, which is compatible with the internal measurement errors. Thus, we find no evidence for an intrinsic [Fe/H] abundance spread in NGC 2298. The typical α-element enrichment in NGC 2298 is overabundant relative to the Sun, and it follows the trend of other metal-poor GCs. We confirm the existence of an Al-enhanced population in this cluster, which is clearly anti-correlated with Mg, indicating the prevalence of the multiple-population phenomenon in NGC 2298.

Central kinematics of the Galactic globular cluster M80

ABSTRACT We use spectra observed with the integral-field spectrograph Multi Unit Spectroscopic Explorer (MUSE) to reveal the central kinematics of the Galactic globular cluster Messier 80 (M80, NGC 6093). Using observations obtained with the recently commissioned narrow-field mode of MUSE, we are able to analyse 932 stars in the central 7.5 arcsec by 7.5 arcsec of the cluster for which no useful spectra previously existed. Mean radial velocities of individual stars derived from the spectra are compared to predictions from axisymmetric Jeans models, resulting in radial profiles of the velocity dispersion, the rotation amplitude, and the mass-to-light ratio. The new data allow us to search for an intermediate-mass black hole (IMBH) in the centre of the cluster. Our Jeans model finds two similarly probable solutions around different dynamical cluster centres. The first solution has a centre close to the photometric estimates available in the literature and does not need an IMBH to fit the observed kinematics. The second solution contains a location of the cluster centre that is offset by about 2.4 arcsec from the first one and it needs an IMBH mass of $4600^{+1700}_{-1400}~\text{M}_\odot {}$. N-body models support the existence of an IMBH in this cluster with a mass of up to 6000 M⊙ in this cluster, although models without an IMBH provide a better fit to the observed surface brightness profile. They further indicate that the cluster has lost nearly all stellar-mass black holes. We further discuss the detection of two potential high-velocity stars with radial velocities of 80–90 $\text{km}\, \text{s}^{-1}$ relative to the cluster mean.

Unveiling the Architecture of a Pulsar–Binary Black Hole Triple System with Pulsar Arrival Time Analysis

Abstract A large number of binary black holes (BBHs) with longer orbital periods are supposed to exist as progenitors of BBH mergers recently discovered with gravitational wave (GW) detectors. In our previous papers, we proposed to search for such BBHs in triple systems through the radial velocity modulation of the tertiary orbiting star. If the tertiary is a pulsar, high-precision and high-cadence observations of its arrival time enable an unambiguous characterization of the pulsar–BBH triples located at several kiloparsecs, which are inaccessible with the radial velocity of stars. The present paper shows that such inner BBHs can be identified through the short-term Rømer delay modulation, on the order of 10 ms for our fiducial case, a triple consisting of a 20 M ⊙ BBH and 1.4 M ⊙ pulsar with P in = 10 days and P out = 100 days. If the relativistic time delays are measured as well, one can determine basically all the orbital parameters of the triple. For instance, this method is applicable to inner BBHs of down to ∼1 hr orbital periods if the orbital period of the tertiary pulsar is around several days. Inner BBHs with ≲1 hr orbital period emit the GW detectable by future space-based GW missions, including LISA, DECIGO, and BBO, and very short inner BBHs with subsecond orbital period can be even probed by the existing ground-based GW detectors. Therefore, our proposed methodology provides a complementary technique to search for inner BBHs in triples, if they exist at all, in the near future.

Internal kinematics of M10 and M71

ABSTRACT We present a study of the internal kinematics of two globular clusters, M10 (NGC 6254) and M71 (NGC 6838), using individual radial velocity (RV) measurements obtained from observations using the Hydra multiobject spectrograph on the WIYN 3.5 m telescope. We measured 120 RVs for stars in M10, of which 107 were determined to be cluster members. In M71, we measured 82 RVs and determined 78 of those measurements belonged to cluster members. Using the cluster members, we determine a mean RV of 75.9 ± 4.0 (s.d.) km s−1 and −22.9 ± 2.2 (s.d.) km s−1 for M10 and M71, respectively. We combined the Hydra RV measurements with literature samples and performed a line-of-sight rotational analysis on both clusters. Our analysis has not revealed a statistically significant rotation in either of these clusters with the exception of the inner region (10–117 arcsec) of M10 for which we find hints of a marginally significant rotation with amplitude Vrot = 1.14 ± 0.18 km s−1. For M10, we calculate a central velocity dispersion of σ0 = 5.44 ± 0.61 km s−1, which gives a ratio of the amplitude of rotation to the central velocity dispersion Vrot/σ0 = 0.21 ± 0.04. We also explored the rotation of the multiple stellar populations identified in M10 and M71 and found rotation (or lack thereof) in each population consistent with each other and the cluster global rotation signatures.

Investigation of the Kinematics of Stars from the Gaia Data Release 2 with Radial Velocities Catalogue Using Scalar and Vector Spherical Harmonics

This paper is a continuation of our work. The proper motions and radial velocities of stars from the Gaia DR2 with RV catalogue have been decomposed. We have confirmed that the Ogorodnikov-Milne model is consistent with the observational data and found the kinematic components that are not described by this model. The beyond-the-model harmonics have been partially identified with the nonlinear terms of the extended Oort model and the parameters of the second-order kinematic model.

Precise Radial Velocities of Cool Low-mass Stars with iSHELL

Abstract The coolest dwarf stars are intrinsically faint at visible wavelengths and exhibit rotationally modulated stellar activity from spots and plages. It is advantageous to observe these stars at near-infrared (NIR) wavelengths (1–2.5 μm) where they emit the bulk of their bolometric luminosity and are most quiescent. In this work, we describe our methodology and results in obtaining precise radial velocity (RV) measurements of low-mass stars using K-band spectra taken with the R ∼ 80,000 iSHELL spectrograph and the NASA Infrared Telescope Facility using a methane isotopologue gas cell in the calibration unit. Our novel analysis pipeline extracts RVs by minimizing the rms of the residuals between the observed spectrum and a forward model. The model accounts for the gas cell, tellurics, blaze function, multiple sources of quasi-sinusoidal fringing, and line spread function of the spectrograph. The stellar template is derived iteratively using the target observations themselves through averaging barycenter-shifted residuals. We have demonstrated 5 m s−1 precision over one-year timescales for the M4 dwarf Barnard’s Star and K dwarf 61 Cygni A, and 3 m s−1 over a month for the M2 dwarf GJ 15 A. This work demonstrates the potential for iSHELL to determine dynamical masses for candidate exoplanets discovered with the NASA Transiting Exoplanet Survey Satellite mission, and to search for exoplanets orbiting moderately active and/or young K &amp; M dwarfs.

A small star hosts a big planet

Exoplanets M dwarfs, the most common type of star, are low-mass objects that emit most of their faint light in the near-infrared, making it difficult to detect any orbiting exoplanets. Morales et al. have observed the nearby M dwarf GJ 3512 in the optical and near-infrared (see the Perspective by Laughlin). Periodic variations in the star's radial velocity show that it hosts a gas giant exoplanet on an eccentric orbit. The authors use simulations to show that such a large exoplanet around such a small star has implications for models of planet formation. Science , this issue p. [1441][1]; see also p. [1382][2] [1]: /lookup/doi/10.1126/science.aax3198 [2]: /lookup/doi/10.1126/science.aay7775

The 1D Relativistic Doppler Formula Is an Incorrect Approximation in Precise Radial Velocity Work

Stellar Doppler velocimetry determines a star's radial velocity $v_r$ via measurement of a redshift, $z$. At precisions below 10 m s$^{-1}$ conversion between the two quantities is complex, and care must be taken to properly account for the effects of relativity. One particular aspect of the problem that bears repeating is that the one-dimensional version of the relativistic Doppler formula, which does not distinguish between the motion of the source and the observer, is incorrect in this context, and indeed does not even provide the correct coefficient for variations in the second-order terms involving $\beta$. Nonetheless, it is often useful to report a redshift in the more familiar units of velocity without a rigorous calculation, and much code already exists that does this. In these cases it is important to clearly document which formula is being used, and I recommend simply (and explicitly) using the approximation $v_r \approx cz$. This choice is trivially inverted, does not misrepresent the degree of relativistic rigor that has been applied in translating between redshift and radial velocity, and is, I believe, the most commonly followed convention in astronomy and cosmology. I also briefly discuss the differences between the Wright & Eastman barycentric correction procedure and the Lindegren & Dravins barycentric radial velocity measure.

Physical properties of the fullerene C60-containing planetary nebula SaSt2-3★

We perform a detailed analysis of the fullerene C60-containing planetary nebula (PN) SaSt2-3 to investigate the physical properties of the central star (B0-1II) and nebula based on our own Subaru/HDS spectra and multiwavelength archival data. By assessing the stellar absorption, we derive the effective temperature, surface gravity, and photospheric abundances. For the first time, we report time variability of the central star's radial velocity, strongly indicating a binary central star. Comparison between the derived elemental abundances and those predicted values by asymptotic giant branch (AGB) star nucleosynthesis models indicates that the progenitor is a star with initial mass of ~1.25 Msun and metallicity Z = 0.001/alpha-element/Cl-rich ([alpha,Cl/Fe] ~ +0.3-0.4). We determine the distance (11.33 kpc) to be consistent with the post-AGB evolution of 1.25 Msun initial mass stars with Z = 0.001. Using the photoionisation model, we fully reproduce the derived quantities by adopting a cylindrically shaped nebula. We derive the mass fraction of the C-atoms present in atomic gas, graphite grain, and C60. The highest mass fraction of C60 (~0.19%) indicates that SaSt2-3 is the C60-richest PN amongst Galactic PNe. From comparison of stellar/nebular properties with other C60 PNe, we conclude that the C60 formation depends on the central star's properties and its surrounding environment (e.g., binary disc), rather than the amount of C-atoms produced during the AGB phase.

Radial Velocity Comparison of Gaia DR2 and RAVE DR5 Survey: A Systematic Offset in Radial Velocities among a Group of Highly Accurate Radial Velocity Stars

Abstract Here, we report a comparative study of radial velocity (RV) data of two major surveys: Gaia Data Release 2 and RAVE Data Release 5. We restricted the sample to stars with relatively accurate RVs ( km s−1 or ≤2%, and km s−1 or ≤2%). The difference between RV Gaia and RVRAVE for a majority of the sample follows a normal distribution with mean = 0.28 km s−1 and σ = 1.49 km s−1. However, we found a very small group of stars (≈0.08% of the total) for which the difference in RVs between the two surveys is significantly larger with an offset of −104.50 km s−1 with σ = 4.92 km s−1. Kinematics based on RV Gaia suggest that most of the group members belong to the Galactic thin disk, which agrees with the group’s metallicity range of −1.2 to +0.5 dex suggesting the offset in RV is probably due to RAVE velocity data for this particular group.

Bowen emission from Aquila X-1: evidence for multiple components and constraint on the accretion disc vertical structure

We present a detailed spectroscopic study of the optical counterpart of the neutron star X-ray transient Aquila X-1 during its 2011, 2013 and 2016 outbursts. We use 65 intermediate resolution GTC-10.4m spectra with the aim of detecting irradiation-induced Bowen blend emission from the donor star. While Gaussian fitting does not yield conclusive results, our full phase coverage allows us to exploit Doppler mapping techniques to independently constrain the donor star radial velocity. By using the component N III 4640.64/4641.84 A we measure Kem = 102 +-6 km/s. This highly significant detection ( >13 sigma) is fully compatible with the true companion star radial velocity obtained from near-infrared spectroscopy during quiescence. Combining these two velocities we determine, for the first time, the accretion disc opening angle and its associated error from direct spectroscopic measurements and detailed modelling, obtaining alpha = 15.5 +2.5 -5 deg. This value is consistent with theoretical work if significant X-ray irradiation is taken into account and is important in the light of recent observations of GX339-4, where discrepant results were obtained between the donor's intrinsic radial velocity and the Bowen-inferred value. We also discuss the limitations of the Bowen technique when complete phase coverage is not available.

Evidence for at least three planet candidates orbiting HD 20794

We explore the feasibility of detecting Earth analogs around Sun-like stars using the radial velocity method by investigating one of the largest radial velocities datasets for the one of the most stable radial-velocity stars HD20794. We proceed by disentangling the Keplerian signals from correlated noise and activity-induced variability. We diagnose the noise using the differences between radial velocities measured at different wavelength ranges, so-called "differential radial velocities". We apply this method to the radial velocities measured by HARPS, and identify four signals at 18, 89, 147 and 330 d. The two signals at periods of 18 and 89 d are previously reported and are better quantified in this work. The signal at a period of about 147 d is reported for the first time, and corresponds to a super-Earth with a minimum mass of 4.59 Earth mass located 0.51 AU from HD20794. We also find a significant signal at a period of about 330 d corresponding to a super-Earth or Neptune in the habitable zone. Since this signal is close to the annual sampling period and significant periodogram power in some noise proxies are found close to this signal, further observations and analyses are required to confirm it. The analyses of the eccentricity and consistency of signals provide weak evidence for the existence of the previously reported 43 d signal and a new signal at a period of about 11.9 d with a semi amplitude of 0.4 m/s. We find that the detection of a number of signals with radial velocity variations around 0.5\,m/s likely caused by low mass planet candidates demonstrates the important role of noise modeling in searching for Earth analogs.

On the mass of the Galactic star cluster NGC 4337

Only a small number of galactic open clusters survives for longer than few hundred million years. Longer lifetimes are routinely explained in term of larger initial masses, particularly quiet orbits, and off-plane birth-places. We derive in this work the actual mass of NGC 4337, one of the few open clusters in the Milky Way inner disk that managed to survive for about 1.5 Gyr. We derive its mass in two different ways. First, we exploit an unpublished photometric data set in the UBVI passbands to estimate -using star counts- the cluster luminosity profile, and luminosity and mass function, and hence its actual mass both from the luminosity profile and from the mass function.This data-set is also used to infer crucial cluster parameters, as the cluster half-mass radius and distance. Second, we make use of a large survey of cluster star radial velocities to derive dynamical estimates for the cluster mass. Under the assumption of virial equilibrium and neglecting the external gravitational field leads to values for the mass significantly larger than those obtained by mean of observed density distribution or with the mass function but still marginally compatible with the inferred values of the invisible mass in form of both low mass stars or remnants of high mass stars in the cluster. Finally, we derive the cluster initial mass by computing the mass loss experienced by the cluster during its lifetime, and adopting the various estimates of the actual mass.

PHOTOMETRIC PROPERTIES OF SELECTED ALGOL-TYPE BINARIES. IX. V548 CYGNI

ABSTRACT We present a new UV light curve of the Algol eclipsing binary V548 Cyg obtained with the Lunar Ultraviolet Telescope. We model the UV light curve together with two previously published (B and V) light curves, primary star radial velocities, and eclipse timings in a unified multi-data-type solution and determine orbital parameters and absolute dimensions. Timing residuals hint at the presence of a third star in the system. This star is possibly the source of the third light that is needed to obtain a good fit to each of the light curves simultaneously. The light–time oscillation in the timing residuals has a period of either or years. The third body orbit inclination would have to be low ( or , respectively) for the third star to have a mass of , which would be expected for a main-sequence star of color , as determined from the light curve solution. In an H–R diagram, the mass-gaining, primary component of V548 Cyg is located between the zero-age and terminal-age main sequence for solar composition stars, and close to the 0.4 Gyr isochrone.