## Direct Method to Compute Doppler Beaming Factors in Binary Stars

Chuanjie Zheng, 1,2 Yang Huang, 2,1 Jifeng Liu, 1,2,3,4 Youjun Lu, 1,2 Henggeng Han, 1 Yuan Tan, 1,2 and Timothy C. Beers 5

1 National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, People's Republic of China; 2 School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China; 3 New Cornerstone Science Laboratory, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, People's Republic of China;

- 4

Institute for Frontiers in Astronomy and Astrophysics, Beijing Normal University, Beijing, 102206, People's Republic of China 5 Department of Physics and Astronomy and JINA Center for the Evolution of the Elements (JINA-CEE), University of Notre Dame, Notre Dame, IN 46556, USA

## ABSTRACT

The Doppler beaming effect, induced by the reflex motion of stars, introduces flux modulations and serves as an efficient method to photometrically determine mass functions for a large number of close binary systems, particularly those involving compact objects. In order to convert observed beamingflux variations into a radial-velocity curve, precise determination of the beaming factor is essential. Previously, this factor was calculated as a constant, assuming a power-law profile for stellar spectra. In this study, we present a novel approach to directly compute this factor. Our new method not only simplifies the computation, especially for blue bands and cool stars, but also enables us to evaluate whether the relationship between beaming flux and radial velocity can be accurately described as linear. We develop a python code and compute a comprehensive beaming-factor table for commonly used filter systems covering main-sequence, subgiant, and giant stars, as well as hot subdwarf and white dwarf stars. Both the code and our table are archived and publicly available at Zenodo.

Keywords: Binary stars (154), Radial velocity (1332), Computational methods (1965)

## 1. INTRODUCTION

Doppler beaming refers to the apparent brightness variations caused by reflex motions of stars in binary systems. In cases where the velocities of targets are far from the speed of light, this effect can be approximately described by a linear relation between the target's radial movement and the observed flux (Loeb &amp; Gaudi 2003; Zucker et al. 2007; Claret et al. 2020a). Thus, this effect can be utilized to derive radial-velocity curves for binary systems photometrically, which are otherwise traditionally measured through expensive and time-consuming spectroscopic observations. However, the modulations in light curves are too small to be detected by groundbased telescopes. Even considering close binary system with radial-velocity semi-amplitude K from a few tens to hundreds of km s -1 , the photometric variations are only on the order of K/c ∼ 10 -4 -10 -3 , where c is the velocity of light.

Corresponding author: Yang Huang

Thanks to great advances in space-based missions, including CoRoT (Baglin et al. 2009; Auvergne et al. 2009), Kepler (Borucki et al. 2010), and TESS (Ricker et al. 2015), continuous light curves with ultra-high precisions of 50-100 parts per million (ppm) are now available for a large sample of field stars. Doppler beaming effects have been clearly detected for binary systems from these precise light curves (Mazeh &amp; Faigler 2010; Bloemen et al. 2011a; Faigler &amp; Mazeh 2011). The largest published sample, 70 beaming binaries, has been described by Tal-Or et al. (2015), who applied a dedicated search algorithm for the light curves obtained by CoRoT. The sample size of such binary systems will increase steadily through the continuously released data from the ongoing TESS survey. Moreover, we can expect to discover distant beaming binaries even at the Galactic center using data from the upcoming China Space Station Telescope (CSST) bulge survey.

It is non-trivial to extract physical parameters of binary systems from the beaming effects by precise light curves alone. The Doppler-beaming factor, characterizing the flux variation caused by the radial motion of star in a binary, is also required to be precisely determined. In the field of high-energy astrophysics, this factor is usually expressed as δ 3 -α , by assuming the emission spectrum to following a power-law profile with an index of α (i.e., F ∝ ν α ), where δ is the Doppler factor 1 (Lind &amp; Blandford 1985; Dermer 1995). On top of that, a constant monochromatic Doppler beaming factor of 3 -α , characterising the linear relation between observed flux and radial velocity, is derived by Loeb &amp; Gaudi (2003) for binary or planet plus host-star systems, given an orbital speed much smaller than c in the Doppler factor. For a given filter, the photon-weighted bandpass-integrated beaming factor can be derived by convolving the filter's transmission curve with the synthetic spectrum of different types of star (e.g., Bloemen et al. 2011a). In this way, a systematic grid of beaming factors, for dozens of commonly used filters (such as in the SDSS, Kepler , TESS, and Gaia photometric systems), were calculated for white dwarfs, main-sequence, and giant stars by Claret et al. (2020a).

However, all of these calculations are based on the assumption that stellar spectra can be well-represented by a power-law profile. In reality, stellar spectra often deviate significantly from this profile. For instance, the blue range of optical spectra is heavily absorbed due to the blanketing effects of metallic lines. Additionally, the entire optical spectra of late-type stars are dominated by molecular bands. These factors may lead to deviations from a constant beaming factor, or at least add complexities to its calculation. To overcome these limitations, we compute the beaming flux at various radial velocities, instead of calculating a constant factor based on the linear relation derived by assuming a power-law stellar spectrum. Using this novel method, we can first determine whether a constant beaming factor adequately describes the relation between beaming flux and radial velocity. If it does indeed exhibit linearity, as assumed in previous studies, this factor can be calculated more straightforwardly using this new approach.

This paper is organised as follows. Section 2 describes our method in detail. Section 3 presents our main results and compares them to previous studies. Finally, a summary is given in Section 4.

1 δ is defined as √ 1 -( v/c ) 2 1+ v r /c .

## 2. THE DOPPLER BEAMING FACTOR

According to relativistic theory, the ratio between the flux F ν ( ) and the cube of the frequency ν 3 should remain constant across different reference frames (see Equation 4.111 in Rybicki &amp; Lightman 1986).

By adopting a relativistic Doppler shift ν = Dν 0 , we obtain:

$$F ( \nu ) _ { \text{obs} } = F ( \nu / D ) _ { 0 } D ^ { 3 } \, \text{ \quad \ \ } ( 1 )$$

Here, F ν ( ) obs refers to the observed flux and the subscript 0 denotes the corresponding quantities in the absence of the relative motion. D is given by:

$$D = \frac { 1 } { \gamma ( 1 + \frac { v _ { c } } { c } ) },$$

and:

$$\gamma = \frac { 1 } { \sqrt { 1 - \frac { v ^ { 2 } } { c ^ { 2 } } } } \,.$$

Here, v r refers to the stellar radial velocity and v refers to stellar space velocity. The monochromatic beaming factor B , as a function of v r , can be derived at a given frequency ν i :

$$B ( v _ { r } ) = \frac { \partial F ( \nu _ { i } ) _ { \text{obs} } } { \partial \frac { v _ { r } } { c } }. \text{ \quad \quad } ( 4 )$$

## 2.1. The Linear Approximation: A Constant Beaming Factor

To search for a binary or a planet plus host-star systems using beaming effects (Loeb &amp; Gaudi 2003), a further refinement of Equation 1 is required, based on two assumptions: 1) the low-speed approximation, where the stellar space velocities v are far from the speed of light; and 2) the stellar spectra are represented by power-law profiles: F ∝ ν α . Under these two assumptions, D is then equal to 1 -v r c ; thus, Equation 1 can be re-written as:

$$F ( \nu ) _ { o b s } = F ( \nu ) _ { 0 } \left ( 1 - ( 3 - \alpha ) \frac { v _ { r } } { c } \right ). \quad \ \ ( 5 )$$

Here, the beaming flux exhibits a linear dependency on the radial velocity v r ; thus, the beaming factor becomes a constant value of ( α -3), as defined in Equation 4. We note that this factor differs from that of Loeb &amp; Gaudi (2003) by a minus sign, simply because they define redshift in the opposite direction from the usual convention.

Naturally, under the power-law approximation, Equation 5 can also be written in wavelength space as:

$$F ( \lambda ) _ { \text{obs} } = F ( \lambda ) _ { 0 } \left ( 1 - B ( \lambda ) \frac { v _ { r } } { c } \right ), \quad \ \ ( 6 )$$

and:

$$B ( \lambda ) = 5 + \frac { d \ln F ( \lambda ) _ { 0 } } { d \ln \lambda }. \quad \quad ( 7 )$$

Figure 1. Synthetic B -band beaming flux at different radial velocities for two binary systems with circular orbits ( e = 0) and semi-amplitude radial velocity K = 300kms -1 . In the left panel, the primary is an extremely metal-poor Solar-type star with T eff = 5600K, log g = 4.0, [Fe/H] = -3.0, and E B ( -V ) = 0 0. . In the right panel, the primary is a metal-rich cool giant star with T eff = 2600K, log g = 2.0, [Fe/H] = +0 5, and moderate reddening, . E B ( -V ) = 0 3. . The red solid lines in both upper panels represent the best linear fit to the synthetic-beaming flux at the B -band, as a function of radial velocity variation of the binary system. The bottom panels show the relative-fitting residuals, defined as the ratio between the fit difference O -C and the maximum flux variation ∆ F (see text for definition) of the linear fits. The red, green, and blue dotted lines in the lower panels mark ratios of ± 1%, ± 2% and ± 5%, respectively. The corresponding deviation values of the D 1 , D 2 , and D 5 indexes (see text for definition) are labelled at the bottom of each upper panel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01033v2-figure-c6864785-8d65-44ca-8735-18e945c722a4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01033v2-figure-f69169b1-838a-4001-9efa-84f4cbd0d425.jpg)

## 2.2. Applications to Widely Used Filter Systems

## 2.3. A New Method: Back to the Definition

To compute beaming factors for the most widely used photometric systems, Bloemen et al. (2011a) and Claret et al. (2020a) modified the monochromatic factor B λ ( ) to the photon-weighted bandpass-integrated beaming factor ¯ , B incorporating filter transmission curves T λ ( ) and synthetic spectra f ( λ ):

$$\bar { B } = \frac { \int _ { \lambda _ { 1 } } ^ { \lambda _ { 2 } } T ( \lambda ) \lambda B ( \lambda ) f ( \lambda ) d \lambda } { \int _ { \lambda _ { 1 } } ^ { \lambda _ { 2 } } T ( \lambda ) \lambda f ( \lambda ) d \lambda }, \quad \quad ( 8 ) \quad \begin{pmatrix} h \\ t \\ a \\ d \end{pmatrix}$$

where λ 1 and λ 2 denote the lower and upper wavelengths for a specified filter.

Following Equation 8, Claret et al. (2020a) presents the first publicly available beaming-factor table for hot white dwarf stars (DA, DB, and DBA), main-sequence stars, and giant stars in commonly adopted broad-band systems, including the SDSS, UBVRI , HiPERCAM, Kepler , TESS, and Gaia bands. The white dwarf atmospheric models from Claret et al. (2020b) and models from Claret et al. (2020a) for main-sequence and giant stars are adopted in their calculations. We note that their calculations didn't consider the reddening effect, which cannot be ignored for disk and bulge stars.

The constant beaming factor defined in previous studies is based on the power-law approximation, which actually deviates significantly from real stellar spectra, especially in the blue range of optical spectra. This leads us to consider whether or not a linear relationship between observed flux and radial velocity is suitable for all types of stars observed in various filter systems, which has never been tested. We thus propose a new method to calculate the beaming flux at various radial velocities and then compute beaming factors based on the direct definition given in Equation 4.

In comparison to the previous approach, the new technique has two advantages: 1) it facilitates the examination of whether a constant beaming factor adequately describes the relation between the observed flux and radial velocity; and 2) the calculation is straightforward and applicable to any kind of spectra; for example, those regions heavily affected by metallic absorption or even the entire optical bands covered by molecular bands for cool stars. In the following subsections, we provide a detailed introduction to the new method for computing the beaming factor, as well as a method to check its validity as a constant.

Figure 2. Predicted power-law profile inferred from the beaming factors in this work (red-solid line) and Claret et al. (2020a) (cyan-dashed line). We note that both profiles are shifted in arbitrary amounts to match the synthetic spectra from ATLAS9 with fixed log g = 5 0, Solar metallicity, and effective temperature labeled in the bottom-left corner. . The upper three panels show results for the U -band, and the bottom three panels are for the TESS band. Transmission curves for both bands are overplotted with gray-solid lines.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01033v2-figure-61deff73-d407-4696-a2c4-0db9d88bf2d2.jpg)

Table 1. Synthetic Spectral Libraries Adopted in this Work

| Library   | Reference                  | Resolving power   | T eff (K)         | log g        | Metallicity        |
|-----------|----------------------------|-------------------|-------------------|--------------|--------------------|
| PHOENIX a | Husser et al. (2013)       | 500,000           | [3000, 12,000]    | [0.0, 6.0]   | [ - 4 . 0, +1.0] e |
| ATLAS9 b  | Castelli & Kurucz (2003)   | 150-300           | [3000, 50,000]    | [0.0, 5.0]   | [ - 4 . 0, +0.5] e |
| TMAP c    | H¨gelmeyer u et al. (2007) | ≥ 10,000          | [50,000, 190,000] | [5.0, 9.0]   | [0.0, +1.0] f      |
| TLUSTY d  | Lanz & Hubeny (2003, 2007) | 1800 g            | [15,000, 55,000]  | [1.75, 4.75] | [0.0, 2.0] h       |

- a The average steps in T eff , log g , and metallicity are 100 K, 0.5 dex, and 0.5 dex, respectively.
- b The average steps in T eff , log g , and metallicity are 250 K, 0.5 dex, and 0.5 dex, respectively.
- c The average steps in T eff , log g , and metallicity are 10,000 K, 0.5 dex, and 0.1 dex, respectively.
- d The average steps in T eff and log g are 1000 K and 0.25 dex, respectively. The metallicity grid is represented by 0.0, 0.1, 0.2, 0.5, 1.0, and 2.0, respectively.

e Here, metallicity is taken to be [Fe/H].

f Here, metallicity refers to the mass ratios between H and He.

- g The resolving power drops to 450 for metal-free cases.

h Here, metallicity refers to Z/Z ⊙ .

## 2.3.1. Synthetic Beaming Flux and Beaming Factor

First, using theoretical stellar spectra, the monochromatic beaming flux at a given radial velocity can be determined by Equation 1:

$$F ( \lambda ) _ { \text{beam} } = F ( D \lambda ) _ { 0 } D ^ { 5 }.$$

For filters with transmission curve T λ ( ) and targets with color excess E B ( -V ), the synthetic beaming flux

¯ F obs observed by energy-counting detectors can be computed with:

$$\bar { F } _ { o b s } = \frac { \int _ { \lambda _ { 1 } } ^ { \lambda _ { 2 } } F ( \lambda ) _ { \text{beam} } R ( \lambda ) T ( \lambda ) d \lambda } { \int _ { \lambda _ { 1 } } ^ { \lambda _ { 2 } } T ( \lambda ) d \lambda }, \quad \ ( 1 0 )$$

where R λ ( ) is the reddening term:

$$R ( \lambda ) = 1 0 ^ { - \frac { k _ { \lambda } E ( B - V ) } { 2. 5 } }, \text{ \quad \ \ } ( 1 1 )$$

where k λ represents the reddening coefficient predicted from the R V = A V E B ( -V ) = 3 1, the Fitzpatrick extinction . law (Fitzpatrick et al. 2019).

For a specific beaming binary with a circular orbit e = 0, a semi-amplitude K of a radial velocity curve and a reddening value of E B ( -V ), the beaming factor of a given filter with transmission curve T λ ( ) can be directly calculated by Equation 4:

$$B _ { X } \left ( K, E ( B - V ) \right ) = \frac { \partial \bar { F } _ { o b s } } { \partial \frac { v _ { r } } { c } } ( - K \leq v _ { r } - v _ { s y s } \leq K ), \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \text{--} \colon \dots \text{ } v \colon \dots \colon \dots \colon \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \colon \dots \end{array} \text{the}$$

where X indicates the name of a given filter and ¯ F obs can be computed from Equation 10.

̸

The method can also work for eccentric orbits ( e = 0), by revising the upper and lower boundaries of radial velocity using the binary orbital parameters (Paddock 1913). A dedicated Python code, BeamingFactor , has been developed to calculate the beaming factor based on Equation 12. By default, the code calculates ¯ F obs from -K to K in steps of 1% K , or 2 5 km s . -1 when K ≥ 250 kms -1 . The beaming factor is then determined from the slope of a linear fit between ¯ F obs and v r . The code accommodates both the well-known filter systems integrated within the code and any passbands provided by the user. For synthetic spectra, this code can directly read spectra files form several widely used libraries: PHOENIX (Husser et al. 2013), ATLAS9 (Castelli &amp; Kurucz 2003), and TAMP (H¨ ugelmeyer et al. 2007), generated with the codes PHOENIX (Hauschildt &amp; Baron 1999), ATLAS9 (Kurucz 1970, 2005), and TAMP (Werner &amp; Dreizler 1999; Werner et al. 2003; Rauch &amp; Deetjen 2003), respectively. Additionally, users can upload their own spectral libraries. The code is made publicly available at Zenodo (Zheng 2024), and a version under continuous development is available at https: //github.com/shift-method/beamingfactor.

## 2.3.2. Beyond the Linear Approximation

In contrast to Claret et al. (2020a), our procedure offers the capability to assess whether a constant beaming factor sufficiently captures the relationship between the observed flux and radial velocity. In other words, our method can work effectively even if the relationship is not linear, as assumed in previous studies. To accomplish this, we define the D 5 index, which counts the fraction of data points deviating from the linear fit by δK K ∆ , with F δK = 5% K . Here ∆ F is the maximum flux change within an orbit: ¯ F obs ( v r = v sys -K ) -¯ F obs ( v r = v sys + K ). Similarly, we can define D 1 and D 2 indexes for different types of studies. If the deviation from linearity is significant, let's say D 5 = 10% (users can define this fraction according to their own requirements), the constant beaming factor will fail to accurately convert the observed flux variation to radial-velocity variation.

Figure 1 presents two examples, both asssuming a binary system in a circular orbit with semi-amplitude K = 300kms -1 . In the left panel, the primary of the binary is a Solar-type extremely metal-poor star ( T eff = 5600K, log g = 4 0, . and [Fe/H] = -3 0) . with zero reddening, E B ( -V ) = 0. The beaming flux at the B -band exhibits a perfectly linear correlation with radial velocity, and the D 5 index is zero. However, in the right panel, clear deviations from a linear fit are evident between the beaming flux in the B -band and radial velocity for a metal-rich cool giant ( T eff = 2600K, log g = 2 0, and [Fe/H] = +0 5) with moderate redden-. . ing, E B ( -V ) = 0 3. . This deviation is a natural result of the strong absorption within the B band, which causes the spectra to deviate significantly from a powerlaw profile. The D 5 index is as high as 21% and D 2 reaches as high as 69%. The radial velocity-dependent beaming factors are required to convert the observed flux variations into a radial-velocity curve. This can be done by our code, if needed, while previous methods cannot achieve this.

As stated earlier, our method is more straightforward for computing the beaming factor than the former method based on the power-law assumption for stellar spectra. This is particularly true for the blue bands and cool stars, whose spectra are dominated by absorption features. Figure 2 shows the comparison between the real synthetic spectra and spectral indexes inferred from the beaming factors computed from our method and from Claret et al. (2020a), which assume a powerlaw profile for stellar spectra. As expected, the spectral slopes predicted by our beaming factor are in excellent agreement with the synthetic ones across all effectivetemperature ranges (i.e., T eff from 4000 K to 50,000 K) in both the blue and optical bands. On the contrary, the slopes given by the beaming factors from Claret et al. (2020a) exhibit significant deviations in the U -band spectra across all effective-temperature ranges and moderate deviations even for cooler temperatures (i.e., T eff = 4000K) in the optical band. These comparisons clearly demonstrate the superiority of our method for calculating the beaming factor.

## 3. BEAMING FACTOR TABLES

Similar to Claret et al. (2020a), a comprehensive table of beaming factors is calculated using the BeamingFactor code developed in this study. The factors are computed for binaries with circular orbits and a semi-amplitude K = 300 km s -1 , which are applicable to most known beaming binaries (Bloemen et al.

2011b; Tal-Or et al. 2015; Pelisoli et al. 2021; Schaffenroth et al. 2022, 2023). The table also includes D 1 , D 2 and D 5 to help evaluate linearity and precision. The results are partly summarized graphically in Figures 3 and 4. We note that the results also hold very well for binaries with smaller or larger values of K . However, if one is interested in very blue passbands (e.g., U ) or cool stars, it is better to calculate them individually using the BeamingFactor code with the exact values of the binary orbital parameters.

## 3.1. Parameter Coverage

The spectral libraries we adopted are PHOENIX (Husser et al. 2013), ATLAS9 (Castelli &amp; Kurucz 2003), TLUSTY OSTAR2002 and BSTAR2006 (hereafter TLUSTY ; Lanz &amp; Hubeny 2003, 2007), and TAMP (H¨ ugelmeyer et al. 2007). Table 1 summarizes the resolving power and atmospheric parameter coverages of the four libraries. Our computations incorporate all available spectra files, so the factor tables follow the same parameter grids as these synthetic spectra libraries. In short, both PHOENIX and ATLAS9 work very well for normal main-sequence, subgiant, and giant stars over wide metallicity ranges, while ATLAS9 provides spectra for hot stars with effective temperature as high as 50,000 K. TLUSTY is dedicated for OB stars with T eff ranging from 15,000 K to 55,000 K. For white dwarf stars and hot subdwarf stars, the TMAP libraries are recommended, as they computed an extended atmosphere grid of H+He composition models with different H/He abundance ratios. For more details about these four libraries, we refer the interested reader to Husser et al. (2013), Castelli &amp; Kurucz (2003), Lanz &amp; Hubeny (2003, 2007), and H¨gelmeyer et al. (2007). u

Different choices of reddening values 2 are also considered for calculating the beaming factors. In total, we calculated beaming factors for six filter systems, including SDSS ( ugriz ), Johnson-Cousins ( UBVRI ), Gaia ( G,G BP , G RP ), Kepler , TESS, and CSST ( NUV and ugrizy ).

## 3.2. Results and Comparison

The results for beaming factors are made publicly available at Zenodo (Zheng 2024). For the listed filter systems, the beaming factors vary between 1 and 15. Consequently, the Doppler-beaming signal is typically on the order of a few thousandths or less when K is around 100 km s -1 . This tiny signal demands highly precise photometry to detect, as it is much smaller compared to other effects such as ellipsoidal modulation and

2 The grid contains values of E B ( -V ): 0.0, 0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.40, 0.60, 0.80, 1.00, 1.25, and 1.5.

reflection, which can reach several tens of percent. Although its amplitude is small, the Doppler-beaming signal is distinguishable from these other effects, enabling clear separation. To facilitate this, a dedicated program named BEER was developed by Faigler &amp; Mazeh (2011) for the photometric detection and analysis of these signals. However, we note that another hard-to-distinguish factor is the presence of varying stellar spots, which are commonly found in M dwarf stars with strong magnetic activity.

In Figures 3 and 4, we present a comparison between the beaming factors computed in this study and those calculated by Claret et al. (2020a). Overall, for factors computed with ATLAS9 , as also adopted by Claret et al. (2020a), the results exhibit consistency, with one notable exception: for cool stars with T eff ≤ 7000 K in the U -band, the values reported by Claret et al. (2020a) are 20% to 40% lower than those obtained in our study. Such differences between factors will result in different values of the photometric semi-amplitude radial velocity, K . Even for hot stars with T eff ≥ 20 000 K, the devi-, ations remain at 10% in the U -band. This discrepancy is easily understood, since their method has problems capturing the real slope of the spectra in the blue region of stars, especially for those that are cooler than 8000 K (see Figure 2). Moreover, we note that reddening can cause variations of a few to ten percent in the derived beaming factors, especially for blue bands. This effect has now been properly accounted for in our tables, with different choices of E B ( -V ) ranging from 0.0 to 1.5 provided.

In contrast to the substantial discrepancy between the ATLAS9 results from this work and those of Claret et al. (2020a), the beaming factors computed with different models exhibit only minor differences, typically within 1%. However, significant discrepancies (as large as 20%) are detected for T eff higher than 30,000 K in V and TESS filters (see right panels of Figure 4). Such deviations may arise from two factors. First, the models use different assumptions. For example, ATLAS9 generates synthetic spectra assuming local thermodynamic equilibrium (LTE), whereas TLUSTY incorporates nonLTE (NLTE) effects. We thus recommend the beaming factors calculated by TLUSTY , given that NLTE effects are important when modeling the atmosphere of hotter stars. Second, modeling synthetic spectra of hotter stars may encounter numerical convergence issues, leading to unexpected artificial features in the resulting spectra. For example, the beaming factors for the TESS filter display a 'U' shape as a function of T eff when using ATLAS9 spectra, whereas the results from TLUSTY spectra show a decreasing trend, as expected. This unexpected 'U'

Figure 3. Beaming factors for the U , B V , , Kepler , Gaia G , and TESS filters, as a function of effective temperature ranging from 3500 K to 10,000 K. All results are computed with ATLAS9 (Castelli &amp; Kurucz 2003) and PHOENIX (Husser et al. 2013) synthetic spectra with Solar metallicity, log g = 4.5. The results from Claret et al. (2020a) are also included for comparison. Readers who wish to use beaming factors within the gray areas should carefully examine the corresponding spectra to ensure no artificial features are present. The decreasing trend observed in all filters except for U reflects the evolution of the spectral slope with T eff . The 'U' shape in the U filter indicates an additional influence from the Balmer jump within the U filter, which otherwise would increase the slope.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01033v2-figure-7094bde1-a662-48d3-9ecc-8feac3fcc192.jpg)

feature is actually caused by artifacts around 5881 ˚ A in the ATLAS9 spectra for T eff above 30,000 K (see Figure A1 in the appendix). We therefore recommend that readers visually inspect the spectra when using beaming factors for hotter stars with T eff ≥ 30 000 K. ,

## 4. SUMMARY

In this study, we present a novel approach for computing the beaming factor for binary systems based on its direct definition. This method is superior to previous approaches that relied on the assumption of a power-law profile for stellar spectra. Our method offers several advantages, particularly in two aspects: 1) calculating the factors straightforwardly, especially for blue bands and cool stars whose spectra are dominated by absorption features; and 2) evaluating whether a constant factor can accurately describe the relationship between beaming flux and radial velocity, which previous methods assumed.

puted for widely used photometric filter systems, including SDSS ( ugriz ), Johnson-Cousins ( UBVRI ), Gaia ( G,G BP , G RP ), Kepler , TESS, and CSST ( NUV and ugrizy ), covering main-sequence, subgiant, and giant stars, as well as hot subdwarf and white dwarf stars. Both the code and table are archived and public available at Zenodo (Zheng 2024). For the code under continuous development, please refer to https://github.com/ shift-method/beamingfactor.

We have developed the Python code BeamingFactor to derive beaming factors for binary systems with known orbital parameters with any photometric bands. Using this code, we provide a table of beaming factors for users to quickly reference. This table is com-

Figure 4. Similar to Figure 3, but for effective temperature ranging from 20,000 K to 50,000 K and with synthetic spectra libraries TLUSTY (Lanz &amp; Hubeny 2003, 2007) and ATLAS9 (Castelli &amp; Kurucz 2003). The factors computed with ATLAS9 present a 'U' shape in the TESS filter, which originates from artificial spectral features near 5881 ˚ A (see Figure A1 in appendix for details). We thus recommend to visually check the spectra when computing beaming factors for hotter stars with T eff ≥ 30,000 K (see gray areas).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01033v2-figure-b9c4595e-bdf8-4941-89d2-201e340c6d7c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01033v2-figure-a3c351f3-ded2-41a5-9d59-65d8c8d44219.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01033v2-figure-3269a472-43ee-45eb-93f9-81e45c3e1922.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01033v2-figure-1f42efd8-fc27-484f-8de6-e6fe214d4be9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01033v2-figure-5b12d10c-dcbe-40c1-95f0-09a4b231d407.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01033v2-figure-db1e9e87-108e-455c-b1f2-dd0bf840e99c.jpg)

## ACKNOWLEDGEMENTS

We would like to thank the referee for helpful comments. This work acknowledges support from the National Natural Science Foundation of China (NSFC) for grants No. 11933004 and No. 11988101. YH acknowledges support from the National Key Basic R&amp;D Program of China via 2023YFA1608303 and the China Manned Space Project with No. CMS-CSST-2021-A08. JFL acknowledges support from the New Cornerstone Science Foundation through the New Corner- stone Investigator Program and the XPLORER PRIZE. T.C.B. acknowledges partial support for this work from grant PHY 14-30152; Physics Frontier Center/JINA Center for the Evolution of the Elements (JINA-CEE), and OISE-1927130: The International Research Network for Nuclear Astrophysics (IReNA), awarded by the US National Science Foundation.

## REFERENCES

Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&amp;A,

506, 411, doi: 10.1051/0004-6361/200810860

Baglin, A., Auvergne, M., Barge, P., et al. 2009, in Transiting Planets, ed. F. Pont, D. Sasselov, &amp; M. J. Holman, Vol. 253, 71-81, doi: 10.1017/S1743921308026252

| Bloemen, S., Marsh, T. R., Østensen, R. H., et al. 2011a, MNRAS, 410, 1787, doi: 10.1111/j.1365-2966.2010.17559.x -. 2011b, MNRAS, 410, 1787, doi: 10.1111/j.1365-2966.2010.17559.x                |
|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977, doi: 10.1126/science.1185402                                                                                                  |
| Castelli, F., & Kurucz, R. L. 2003, in Modelling of Stellar Atmospheres, ed. N. Piskunov, W. W. Weiss, & D. F. Gray, Vol. 210, A20,                                                                |
| Claret, A., Cukanovaite, E., Burdge, K., et al. 2020a, A&A, 641, A157, doi: 10.1051/0004-6361/202038436 -. 2020b, A&A, 634, A93,                                                                   |
| doi: 10.1051/0004-6361/201937326                                                                                                                                                                   |
| Dermer, C. D. 1995, ApJL, 446, L63, doi: 10.1086/187931 Faigler, S., & Mazeh, T. 2011, MNRAS, 415, 3921, doi: 10.1111/j.1365-2966.2011.19011.x                                                     |
| Fitzpatrick, E. L., Massa, D., Gordon, K. D., Bohlin, R., & Clayton, G. C. 2019, ApJ, 886, 108, doi: 10.3847/1538-4357/ab4c3a                                                                      |
| Hauschildt, P. H., & Baron, E. 1999, Journal of Computational and Applied Mathematics, 109, 41, doi: 10.48550/arXiv.astro-ph/9808182                                                               |
| Dwarfs, ed. R. Napiwotzki & M. R. Burleigh, 187, doi: 10.48550/arXiv.astro-ph/0610744 Husser, T. O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6, doi: 10.1051/0004-6361/201219058 |
| Kurucz, R. L. 1970, SAO Special Report, 309 -. 2005, Memorie della Societa Astronomica Italiana Supplementi, 8, 14                                                                                 |
| Lanz, T., & Hubeny, I. 2003, ApJS, 146, 417, doi: 10.1086/374373 -. 2007, ApJS, 169, 83, doi: 10.1086/511270                                                                                       |

| Lind, K. R., & Blandford, R. D. 1985, ApJ, 295, 358, doi: 10.1086/163380                                                                                                       |
|--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Loeb, A., & Gaudi, B. S. 2003, ApJL, 588, L117, doi: 10.1086/375551                                                                                                            |
| Mazeh, T., & Faigler, S. 2010, A&A, 521, L59, doi: 10.1051/0004-6361/201015550                                                                                                 |
| Paddock, G. F. 1913, PASP, 25, 208, doi: 10.1086/122237 Pelisoli, I., Neunteufel, P., Geier, S., et al. 2021, Nature Astronomy, 5, 1052, doi: 10.1038/s41550-021-01413-0       |
| Rauch, T., & Deetjen, J. L. 2003, in Astronomical Society of the Pacific Conference Series, Vol. 288, Stellar Atmosphere Modeling, ed. I. Hubeny, D. Mihalas, &                |
| Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015,                                                                                                                       |
| Journal of Astronomical Telescopes, Instruments, and Systems, 1, 014003, doi: 10.1117/1.JATIS.1.1.014003                                                                       |
| Rybicki, G. B., & Lightman, A. P. 1986, Radiative Processes in Astrophysics                                                                                                    |
| Schaffenroth, V., Barlow, B. N., Pelisoli, I., Geier, S., & Kupfer, T. 2023, A&A, 673, A90, doi: 10.1051/0004-6361/202244697                                                   |
| Schaffenroth, V., Pelisoli, I., Barlow, B. N., Geier, S., & Kupfer, T. 2022, A&A, 666, A182, doi: 10.1051/0004-6361/202244214                                                  |
| Tal-Or, L., Faigler, S., & Mazeh, T. 2015, A&A, 580, A21, doi: 10.1051/0004-6361/201526425                                                                                     |
| Werner, K., Deetjen, J. L., Dreizler, S., et al. 2003, in                                                                                                                      |
| Astronomical Society of the Pacific Conference Series, Vol. 288, Stellar Atmosphere Modeling, ed. I. Hubeny, D. Mihalas, & K. Werner, 31, doi: 10.48550/arXiv.astro-ph/0209535 |
| Werner, K., & Dreizler, S. 1999, Journal of Computational and Applied Mathematics, 109, 65, doi: 10.48550/arXiv.astro-ph/9906130                                               |
| Zheng, C. 2024, beamingfactor, 2.0.0, Zenodo,                                                                                                                                  |
| doi: 10.5281/zenodo.13049419                                                                                                                                                   |
| Zucker, S., Mazeh, T., & Alexander, T. 2007, ApJ, 670, 1326, doi: 10.1086/521389                                                                                               |

## APPENDIX

## A. SYNTHETIC SPECTRA COMPARISONS BETWEEN ATLAS9 AND TLUSTY

We present a comparison of spectra from ATLAS9 and TLUSTY in the wavelength range of 5600 to 6200 ˚ . As illustrated A in Figure A1, the absorption features near 5881 ˚ A in the ATLAS9 spectra exhibit abnormal profiles and show significant deviation from those in the TLUSTY spectra when T eff exceeds 30,000 K.

Figure A1. Normalized spectra from ATLAS9 (left panel) and TLUSTY (right panel) are shown for T eff of 23,000, 29,000, 35,000, and 40,000 K. All spectra have Solar abundances with a surface gravity log g of 4.5 dex. The resolving power of the TLUSTY spectra has been degraded to match that of ATLAS9 , which is approximately 300.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01033v2-figure-24f85c0d-b755-4907-93a8-af9a6bcfe7fd.jpg)
Chuanjie Zheng Yang Huang Jifeng Liu Youjun Lu Henggeng Han Yuan Tan Timothy C. Beers
Direct Method to Compute Doppler Beaming Factors in Binary Stars
The Doppler beaming effect, induced by the reflex motion of stars, introduces
flux modulations and serves as an efficient method to photometrically determine
mass functions for a large number of close binary systems, particularly those
involving compact objects. In order to convert observed beaming-flux variations
into a radial-velocity curve, precise determination of the beaming factor is
essential. Previously, this factor was calculated as a constant, assuming a
power-law profile for stellar spectra. In this study, we present a novel
approach to directly compute this factor. Our new method not only simplifies
the computation, especially for blue bands and cool stars, but also enables us
to evaluate whether the relationship between beaming flux and radial velocity
can be accurately described as linear. We develop a python code and compute a
comprehensive beaming-factor table for commonly used filter systems covering
main-sequence, subgiant, and giant stars, as well as hot subdwarf and white
dwarf stars. Both the code and our table are archived and publicly available at
http://doi.org/10.5281/zenodo.13049419.