1. INTRODUCTION
Apsidal motion in detached eclipsing binaries with eccentric orbits is one of the methods used to examine the effects of the tidal distortion on the dynamical evolution of the binary and to test stellar evolution models. Orbital elements and absolute parameters of the components of such systems are required to estimate the internal structure constants (ISCs) of the components which are tests of stellar evolutionary models. Based on this motivation, we have examined the V990 Her binary system, which has not been studied in detail before.
V990Her (HD169888, GSC0158101786)(α
2000 =18h25m56s.71, δ
2000 = 21◦36'21''
.34) is an early type (A0, Cannon
& Pickering, 1993) eclipsing binary system, which was classified as
an Algol type by Kazarovets et al. (1999).
Later the system was listed in the eccentric eclipsing binary systems catalogue
published by Otero et al. (2006) reporting
the light elements as T
0(HJD) = 2448048.755, P = 8
d
.19329. McDonald et al.
(2012) reported the effective temperature T
eff = 7514 K and luminosity
In the following section (§2), we present information about the photometric and spectroscopic observations and reduction procedures of the data used in our analysis of the V990 Her system. The analysis of the data is given in § 3. In § 4 and § 5 we present the results and the discussion.
2. OBSERVATIONS AND DATA REDUCTION
2.1. Photometric Data
In this study, both new observational and survey data (available to the public) were used for the photometric analysis.
Time series photometric observations of V990 Her were made between 2018-2019
using the 0.6 m automated telescope (T60) of the TÜBİTAK National Observatory
(TUG), which is equipped with a 2048×2048 pixels FLI ProLine 3041-UV CCD camera
with a pixel size of 15×15µm2. The Akdeniz
University 0.25 m telescope (AUT25) which is equipped with a 2184×1472 pixels
QSI632ws CCD camera with a pixel size of 6.8
Additional photometric data of the system have been collected from four sources, the Hipparcos and Tycho Catalogues (Esa 1997), the All Sky Automated Survey (ASAS-3) (Pojmanski & Maciejewski 2004) Database, the Multi-site All-Sky CAmeRA (MASCARA) (Talens et al. 2017) Survey and the QuickLook Pipeline (QLP) light curves (LCs) observed by the Transiting Exoplanet Survey Satellite (TESS) and produced by the MIT QuickLook Pipeline (QLP) (Huang et al. 2020). To obtain the times of minima to be used in the O − C analysis in § 3.1 we used the MASCARA LCs published by Burggraaff et al. (2018) and the TESS QLP LCs (Sector 26) which are available at the Mikulski Archive for Space Telescopes (MAST) as a High Level Science Product. The photometric data taken from Hipparcos photometry and ASAS-3 were not sufficient to specify the minimum depth of the system. Therefore, we prefer to use only our T60 and AUT25 photometric data for the LCs analysis (see § 3.3 and § 3.6).
2.2. Spectroscopic Data
Spectroscopic data were collected by using the eShel spectrograph, which has a resolving power of 12000 covering the 4045-8100 Å wavelength range in 27 orders that is connected to the 60 cm telescope (UBT60) of Akdeniz University with a 20-m fiber. More information on the UBT60 telescope system is given by Bakıs, et al. (2020). We collected 29 spectra of V990 Her on different observing nights during June-September 2017 and August-October 2018 seasons (see Table 1). In Table 1, the columns show the order of the observing night, mid-observing time of the spectra in heliocentric Julian date, total exposure time, S/N ratio value around 5500 Å and the orbital phase calculated with the ephemeris of the primary minimum given in equation (1).
No | Date | Julian Date | Exp. Time |
|
Orbital Phase |
---|---|---|---|---|---|
(YYYYMMDD) | (HJD) | (s) | |||
1 | 20170516 | 2457890.3826 | 3000 | 40 | 0.178 |
2 | 20170606 | 2457911.4226 | 3000 | 40 | 0.746 |
3 | 20170611 | 2457916.4894 | 3000 | 110 | 0.365 |
4 | 20170612 | 2457917.5082 | 3000 | 110 | 0.489 |
5 | 20170622 | 2457927.4824 | 3000 | 40 | 0.706 |
6 | 20170702 | 2457937.4004 | 3000 | 110 | 0.917 |
7 | 20170703 | 2457938.4406 | 3000 | 50 | 0.044 |
8 | 20170704 | 2457939.5348 | 3000 | 100 | 0.177 |
9 | 20170709 | 2457944.4475 | 3000 | 100 | 0.777 |
10 | 20170711 | 2457946.4449 | 3000 | 135 | 0.021 |
11 | 20170719 | 2457954.4968 | 3000 | 100 | 0.003 |
12 | 20170720 | 2457955.3306 | 3000 | 110 | 0.105 |
13 | 20170723 | 2457958.4583 | 4200 | 140 | 0.487 |
14 | 20170802 | 2457968.3904 | 3000 | 140 | 0.699 |
15 | 20170918 | 2458015.3474 | 6600 | 50 | 0.430 |
16 | 20170920 | 2458017.3004 | 6000 | 50 | 0.669 |
17 | 20170928 | 2458025.2668 | 7200 | 40 | 0.641 |
18 | 20171011 | 2458038.2554 | 6000 | 50 | 0.226 |
19 | 20180828 | 2458359.3891 | 3000 | 110 | 0.421 |
20 | 20180829 | 2458360.3098 | 3000 | 110 | 0.533 |
21 | 20180913 | 2458375.3204 | 3000 | 70 | 0.365 |
22 | 20180917 | 2458379.3137 | 3000 | 80 | 0.853 |
23 | 20180918 | 2458380.3664 | 3000 | 110 | 0.981 |
24 | 20180919 | 2458381.3540 | 3000 | 100 | 0.102 |
25 | 20180922 | 2458384.2969 | 3000 | 110 | 0.461 |
26 | 20181002 | 2458394.3723 | 3000 | 40 | 0.691 |
27 | 20181003 | 2458395.3006 | 3000 | 40 | 0.804 |
28 | 20181005 | 2458397.2865 | 3000 | 140 | 0.046 |
29 | 20181007 | 2458399.2858 | 3000 | 125 | 0.290 |
*S/N ratio refers to the wavelength region around 5500 Å.
IRAF was used for the reduction of the spectroscopic data and for the measurement of the components’ radial velocities (RVs). In each observing night, bias, dark frames and flat spectra from a tungsten lamp were obtained for image reduction. Order aperture extractions and elimination of scattered light through the orders were made by using the IRAF échelle package. The wavelength calibration was performed with thorium-argon lamp spectra which were taken before and after each object spectrum.
3. ANALYSIS
3.1. Improved Ephemerides and Apsidal Motion
Before attempting the photometric and spectroscopic analysis, the available ephemerides of the system given by Otero et al. (2006) were improved. Some photometric minimum data of V990 Her were published by Kreiner (2004), which were used, together with the newly obtained minima times, for the O−C analysis. We measured two secondary and one primary times of minima from the MASCARA LCs and two primary and one secondary times of minima from the TESS QLP LCs. The measured original times of minima in units of BJD from TESS LCs were converted to HJD unit by using online applets4 Eastman et al. 2010). Additionally, we obtained three primary and one secondary times of minima during our photometric observations. The Kwee & van Woerden (1956) method was used to determine the new times of minima and their corresponding errors for our observations. The obtained times of minima are listed in Table 2 in line with the literature; their uncertainty, epoch and O−C values computed with the linear light elements in equation (1), type of minima (“pri” for primary and “sec” for secondary minimum), observational method (“ccd” for charge-coupled device observations) are also listed; the references are given in the last column.
JD (Hel.) 2400000+ | Error |
|
|
Type | Method | Reference |
---|---|---|---|---|---|---|
48048.7550 | 0.0 | 0.0030 | pri | ccd | Otero et al. (2006) | |
48052.401 | 0.5000 | -0.4477 | sec | ccd | Otero et al. (2006) | |
53177.7065 | 0.0016 | 626.0 | -0.0601 | pri | ccd | Zakrzewski,B. (SuperWASP) unpublished |
53746.7457 | 0.0024 | 695.5 | -0.4477 | sec | ccd | Zakrzewski,B. (ASAS-3) unpublished |
53751.2885 | 0.0027 | 696.0 | -0.0100 | pri | ccd | Zakrzewski,B. (ASAS-3)unpublished |
57192.4960 | 0.001 | 1116.0 | 0.0056 | pri | ccd | This study (MASCARA) |
57212.5300 | 0.002 | 1118.5 | -0.4437 | sec | ccd | This study (MASCARA) |
57253.4970 | 0.001 | 1123.5 | -0.4433 | sec | ccd | This study (MASCARA) |
58249.4300 | 0.010 | 1245.0 | 0.0021 | pri | ccd | This study (AUT25) |
58290.4000 | 0.01 | 1250.0 | 0.0055 | pri | ccd | This study (AUT25) |
58634.5000 | 0.01 | 1292.0 | -0.0137 | pri | ccd | This study (AUT25) |
58646.3500 | 0.01 | 1293.5 | -0.4537 | sec | ccd | This study (AUT25) |
59011.4045 | 0.0005 | 1338.0 | -0.0016 | pri | ccd | This study (TESS) |
59019.5977 | 0.0006 | 1339.0 | -0.0017 | pri | ccd | This study (TESS) |
59023.2540 | 0.001 | 1339.5 | -0.4421 | sec | ccd | This study (TESS) |
* Not used in the apsidal motion analysis.
When we performed a linear O − C fit to all the primary and secondary minimum times given in Table 2 by using the unweighted least squares method, we found that the O−C differences between the observed and the calculated times of minima regarding to the primary minimum behave in opposite phases compared to those of the secondary minima. Then we applied the unweighted least squares fits to the primary and secondary eclipsing times separately and obtained the new linear ephemerides as follows (numbers in parentheses indicate errors in the last digit, and E is the number of cycles):
This result is a clear sign that the system shows apsidal motion. Therefore, we decided to perform an apsidal motion analysis using the times of minima (except for 2453177.7065 due to its large deviation from the expected location in the O − C diagram) given in Table 2. To model the observed O − Cs we used the equation given by Giménez & Bastero (1995) (equation 15 in their paper). We applied a differential correction method to represent the corrections to the input values of the parameters, which are selected as adjustable, and used the unweighted least squares method to obtain them. Unfortunately, when we selected all five parameters as adjustable, we were unable to get a convergent solution. The main reason for this may be that the range covered by the times of minima available is short compared to the apsidal period of the system.
In order to reduce the number of the adjustable parameters and thus obtain an acceptable solution, we decided to determine the orbital eccentricity from the shape of the system’s LC. In a LC of an eccentric eclipsing binary system, ecosw and esinw, the combinations of orbital eccentricity (e) with the longitude (w) of the periastron, are related with the amount of displacement of the secondary minimum with respect to the primary minimum and with the durations of the minima, respectively.
For relatively small eccentricities, these relationships are given as ecosw = πδ/(1 + csc2 i) and esinw = (D 2−D 1)/(D 2+D 1), where δ is the amount of displacement of the secondary minimum with respect to the primary minimum, i is the orbital inclination and D 1 and D 2 are the duration of the primary and secondary minimum, respectively. As a result, from the LCs (see Figure 7) of the system we obtained the values ecosw = −0.0867 and esinw = −0.0899 which yield the orbital eccentricity as e = 0.125.
We managed to obtain a solution by keeping this value of eccentricity constant. The results are given in Table 3 associated with their internal errors. In Table 3, the parameters T 0 , P, e, w 0 and w˙ denote the initial epoch, orbital period, eccentricity, longitude of the periastron point at the initial epoch, and apsidal motion rate, respectively. The observed O − C variations (symbols) of V990 Her are shown in Figure 1 concerning the theoretical curves (solid lines) computed with the apsidal motion elements given in Table 3.
Parameter | Value | Error |
---|---|---|
|
2448048.513 | 0.004 |
|
8.193315 | 0.000003 |
|
0.125 | (fixed) |
|
226.2 | 0.8 |
|
0.00055 | 0.00011 |
According to the results given in Table 3, V990 Her has an apsidal motion with a period of U = 14683 ± 2936 years. Due to the relatively long apsidal period compared with the time span of about 30 years which is covered by the observations, the apsidal motion parameters have larger errors. Therefore, the results given in Table 3 should be considered as a preliminary ones. Thus, for V990 Her more observations, which will be performed in the future, are essential.
3.2. Spectroscopic Orbital Parameters
MgII, SiII, FeI, FeII and HI lines can be clearly seen in the spectrum of V990 Her. RVs were obtained from the orders containing the relatively sharp lines of MgII (4481Å), FeI (4957Å), FeI (5616Å) and HI (6563Å). The studied regions are shown in Figure 2. From top to bottom along the y-axis in Figure 2, the normalized spectra are plotted in accordance with the decreasing HJD times given in Table 1.
Gaussian profiles were applied to the central parts of the lines to measure the RVs and then a barrycentric correction was applied to correct the RVs. Measured RVs were listed in Appendices A, B, C and D. The least squares method has been applied to fit the RV equation to the measured RVs. We used the orbital period derived in equation (1) and fixed it during the iterations. The measured RVs and fitted curves (solid red lines) are shown in Figure 3 for each selected line separately. The spectroscopic orbital parameters, such as the longitude of the periastron, ω, orbital eccentricity, e, mass ratio, q = K 1 /K 2 and velocity semi-amplitude of the components, K 1,2, were obtained from different lines and their arithmetic mean values are given in Table 4.
3.3. Photometric Elements
The previously mentioned three photometric data sets obtained in V RI bands were used for the photometric analysis of V990 Her to determine the photometric elements of the system and the light contributions of the components at different phases, which were used to model the component spectra in § 3.4 and § 3.5. The LC solutions were performed with the program phoebe 1.0 (Prša & Zwitter 2005), which is an extension of the wd program (Wilson & Devinney 1971).
The improved light elements given in equation (1) were used in the LCs analysis. The effective temperature of V990 Her was given as T eff = 7514 K by McDonald et al. (2012). However, we obtained a more reliable temperature for the primary component as T eff1 = 8000 K directly from the spectrum of the system (see § 3.5). The temperature for the primary component T eff1, the improved period given in § 3.1, the eccentricity e, and the mass ratio q given in Table 4 were fixed during the LC solution. Convergence was allowed for the orbital inclination i, the longitude of the periastron point ω, the dimensionless surface potentials of the components Ω1 and Ω2, the effective temperature of the secondary component T eff,2 , and the relative monochromatic luminosity of the primary component l 1 /l tot. The bolometric albedos A 1,2 and gravitational darkening coefficients g 1,2 were adopted as 1.0, as suggested for early-type stars. The bolometric limbdarkening coefficients x bol1,2 were computed automatically by phoebe 1.0 with the internal tables of the programme which were generated according to the models given by van Hamme (1993). The results of the solution are listed in Table 5.
Element | Primary | System | Secondary |
---|---|---|---|
|
8.193314 (fixed) | ||
|
|
||
|
0.88 (fixed) | ||
|
0.09(fixed) | ||
|
|
||
|
|
||
|
8000 (fixed) |
|
|
|
|
|
|
|
0.198 | 0.199 | |
|
0.549 | 0.547 | |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.777 |
3.4. Disentangling the Components Spectrum
One of the advantages of having numerous spectra taken at different phases is that they allow to obtain the spectrum of the components separately. Thus, we used the korel program provided by Hadrava (1995) which needs approximate orbital elements and the light contributions of the components at the phases in which the individual spectrum were taken. The orbital parameters and light contributions of the components obtained in § 3.2 and § 3.3 were used as initial values for the korel solution. The orbital period P given in equation (1) was used and fixed during the iterations. The korel solution gives the RVs (see Table 6) of the components and the orbital parameters (see Table 7) of the binary system. The korel solution agrees very well with those obtained from individual lines (see Table 4). The disentangled spectra in the H α region for both components can be seen in Figure 4.
No | HJD | Phase | Primary |
|
Secondary |
- |
---|---|---|---|---|---|---|
(+2400000) | km s |
km s |
km s |
km s |
||
1 | 57890.3826 | 0.178 | -113.42 | -4.06 | 51.26 | -5.99 |
2 | 57911.4226 | 0.746 | 32.72 | -5.56 | -119.83 | -9.04 |
3 | 57916.4894 | 0.365 | -84.21 | -4.91 | 14.56 | -8.48 |
4 | 57917.5082 | 0.489 | - | - | - | - |
5 | 57927.4824 | 0.706 | 38.85 | -0.81 | -112.49 | -0.13 |
6 | 57937.4004 | 0.917 | -2.22 | 1.07 | -65.03 | -1.55 |
7 | 57938.4406 | 0.044 | -62.92 | -3.32 | 1.15 | 0.55 |
8 | 57939.5348 | 0.177 | -111.31 | -2.14 | 51.65 | -5.38 |
9 | 57944.4475 | 0.777 | 38.67 | 3.79 | -107.78 | -0.85 |
10 | 57946.4449 | 0.021 | - | - | - | - |
11 | 57954.4968 | 0.003 | - | - | - | - |
12 | 57955.3306 | 0.105 | -85.53 | 1.04 | 32.20 | 0.88 |
13 | 57958.4583 | 0.487 | - | - | - | - |
14 | 57968.3904 | 0.699 | 43.63 | 4.10 | -107.86 | 4.35 |
15 | 58015.3474 | 0.430 | - | - | - | - |
16 | 58017.3004 | 0.669 | 63.50 | 25.82 | -83.33 | 26.77 |
17 | 58025.2668 | 0.641 | 57.11 | 22.97 | -103.97 | 2.10 |
18 | 58038.2554 | 0.226 | -77.16 | 37.56 | 86.45 | 23.10 |
19 | 58359.3891 | 0.421 | - | - | - | - |
20 | 58360.3098 | 0.533 | 2.17 | -1.44 | -71.89 | -0.57 |
21 | 58375.3204 | 0.365 | - | - | - | - |
22 | 58379.3137 | 0.853 | 20.40 | 2.04 | -89.44 | -1.31 |
23 | 58380.3664 | 0.981 | - | - | - | - |
24 | 58381.3540 | 0.102 | -83.37 | 1.92 | 32.56 | 2.71 |
25 | 58384.2969 | 0.461 | - | - | - | - |
26 | 58394.3723 | 0.691 | 39.51 | 0.27 | -110.91 | 0.97 |
27 | 58395.3006 | 0.804 | 29.90 | -0.34 | -100.61 | 1.02 |
28 | 58397.2865 | 0.046 | -58.81 | 2.11 | 8.22 | 6.12 |
29 | 58399.2858 | 0.290 | -113.89 | -6.95 | 45.34 | -9.15 |
Element | Value |
---|---|
|
- |
|
- |
|
- |
|
- |
|
- |
|
- |
3.5. Atmosphere Modelling
We have a spectrum of the system in the middle of the secondary minimum of the LC
where
We know that the components are not early type since there are no He I lines in the composite or disentangled spectra, which means that their temperature is no higher than ≈ 10000 K, and therefore, the local thermodynamic equilibrium assumption could be used to construct the atmosphere models. atlas9 atmosphere models computed by Castelli & Kurucz (2004) and synthe program of Kurucz (1993) were used for the generation of the atmospheric models and the synthetic spectra for the components. In the atlas9 code, Z=0.03 and ξ = 2 km s−1 were assumed for the metallicity and micro turbulence velocity, respectively. This method has been applied to the wavelength regions where the dominant lines reside. In Figure 5, the observational data and the synthetic spectra are given for four different wavelength regions as examples.
The atmospheric model parameters of the primary component, such as the temperature, T eff1=8000 ± 200 K, the surface gravity, logg 1 = 4.0 ± 0.1 and the rotational velocity v rot1 =28 ± 5 km s−1 were determined.
Moreover, we have a spectrum at orbital phase φ = 0.0. Using a similar approach, the atmospheric model parameters of the secondary component were estimated as T eff2=7500±200 K, log g 2=4.0±0.1 and v rot2 =25±5 km s−1. At three different orbital phases (φ = 0.699, 0.533, 0.290), the synthetic spectra of the components were combined using their light contributions and their specific RVs. The observed spectra taken at the indicated maximum phases and the combined theoretical spectra are shown together in Figure 6.
3.6. Simultaneous Radial Velocity and Light Curves Solution
The RVs of the components obtained from the korel solution and the multi-colour LCs of the system were analyzed simultaneously with phoebe 1.0. The initial orbital parameters were used from the averaged values given in Tables 4 and 7. We used the improved light elements given in equation (1). The orbital period P and the effective temperature of the primary component T eff1 were fixed as in § 3.3, the other orbital and photometric parameters were all adjusted in the analyses.
The results of the simultaneous solution are listed in Table 8. The observed LCs in VRI bands and the radial velocity curves (RVCs) of the components are shown along with the theoretical models in Figure 7. The Roche geometry of V990 Her at phases φ=0.0, 0.25 and 0.45 can be seen in Figure 8.
Element | Primary | System | Secondary |
---|---|---|---|
|
8.193314 (fixed) | ||
|
|
||
|
|
||
|
|
||
|
|
||
|
|
||
|
|
||
|
|
||
|
8000 (fixed) |
|
|
|
|
|
|
|
0.198 | 0.199 | |
|
0.549 | 0.547 | |
|
|
0.425 | |
|
|
0.436 | |
|
|
0.447 | |
|
|
|
|
|
0.800 |
4. RESULTS
4.1. Physical Properties and Evolutionary Status
Using the parameters given in Tables 7 and
8, the astrophysical parameters of the components and the distance of V990 Her
were derived and are given in Table 9.
The masses of the components (2.01
Element | Unit | Primary | Secondary |
---|---|---|---|
|
|
2.01 |
1.83 |
|
|
2.22 |
2.12 |
|
K | 8000 |
7570 |
|
|
1.25 |
1.12 |
|
mag | 1.6 |
2.0 |
log |
cgs | 4.06 |
4.04 |
|
km s |
28 |
25 |
|
km s |
13.6 |
13.1 |
|
pc | 212 |
Bolometric corrections (BC) for the primary and secondary components were calculated as 0.027 and 0.034 mag using T eff-BC tables of Flower (1996). The bolometric magnitudes M bol given in Table 9 and BCs were used to derive the absolute magnitudes of the components as M V1 = 1 m .60 ± 0.11 and M V2 = 1 m .93 ± 0.12, respectively. By using the V J = V T −0.09×(B T −V T ) relation and Tycho magnitudes (B T = 7 m .950, V T = 7 m .710) given by Esa (1997), we derived the total visual magnitude of the system as m V = 7 m .688. To obtain the E(B − V ) color excess and the distance of the system we used the total u 0 g 0 r 0 i 0 z 0 magnitudes of the system given by Pickles & Depagne (2010) in the following transformation relationships of Smith et al. (2002):
Using the light contribution of the components given in Table8 together with the
4.2. Apsidal Motion and Internal Structure Constant
To test the observed preliminary results concerning the apsidal motion of V990 Her given in § 3.1, we can use the results of previous section. As is known, the total apsidal motion in a binary system occurs mainly for two contributions: (1) a classical part which is due to tidal interactions of the oblate components, and (2) a relativistic contribution which is a direct result of the General Relativity Theory.
The relativistic contribution
which was first derived by Levi-Civita (1937). Here, m 1,2 are the masses of the components in units of solar mass, e is the orbital eccentricity and P is the orbital period in units of days. By using the related parameters given in Tables 3 and 9 we obtained w˙ rel = 0◦ .00033 cycle−1 and so w˙ cl = 0◦ .00022 cycle−1.
On the other hand, the relation for the classical part of the apsidal motion rate was given by Kopal (1978) as follows
where k 2i (i=1, 2) are internal structure constants (ISCs) of the components, while c 2i are as follows
Here
Using the apsidal and absolute parameters given in Tables 1 and4 we
obtained the observed average ISC for the system as
To obtain the theoretical ISCs of the components we used the stellar models given
by Claret & Gimenez (1992), which
assume X=0.649, Z=0.03 and take mass loss and
core overshooting into account. By the interpolation we obtained log
L=1.30 and 1.13
5. DISCUSSION
The absolute parameters of the eccentric binary V990 Her has been obtained with a
good precision;
Different E(B −V ) color excess values of 0 m .166 (Jayasinghe et al. 2018), 0 m .11 (Gontcharov & Mosenkov 2017) and 0 m .066 (Kervella et al. 2019) are also available in the literature. According to these values, we estimated the distance of the system as 172 ± 17, 186 ± 18 and 198 ± 19 pc, respectively. Our photometric distance (d = 212 ± 21 pc), estimated in § 4.1, is closer to the one given by Gaia than the distances determined using the above color excesses.
The apsidal motion in V990 Her puts it into a special place and allows us to estimate the structure of the component stars. However, the time span of the available data we have is a very small portion (≈ 0.2%) of the apsidal cycle. This situation creates problems in the apsidal motion analysis of the system and makes it difficult to determine the apsidal motion rate and orbital eccentricity together. For now, based on our findings we can say that the orbital eccentricity of the system has a value between 0.09 and 0.13. More accurate results can be obtained with more observations to be made in the future. Lastly, the apsidal motion of V990 Her has a period of about 14700 years and is predominantly (60%) due to relativity.