arxiv: v1 [astro-ph.sr] 17 Jun 2016

Size: px

Start display at page:

Download "arxiv: v1 [astro-ph.sr] 17 Jun 2016"

Kelly Black
5 years ago
Views:

1 Research in Astron. Astrophys Vol. X No. XX, Research in Astronomy and Astrophysics arxiv: v1 [astro-ph.sr] 17 Jun 2016 Physical and Geometrical Parameters of CVBS XI: COU1511 (HIP12552) Mashhoor A. Al-Wardat 1, M. H. El-Mahameed 2, Nihad A. Yusuf 3,4, Awni M. Khasawneh 5, Suhail G. Masda 6 1 Physics Department, Al al-bayt University, PO Box: , Mafraq, Jordan; mwardat@aabu.edu.jo 2 Department of Physics, Al-Hussein Bin Talal University, Maan, Jordan 3 Department of Physics, Yarmouk University, Irbid, 21163, Jordan 4 Princess Sumaya University for Technology, Amman, Jordan 5 Royal Jordanian Geographic Center, Amman, Jordan 6 Astrophysikalisches Institut und Universitäts-Sternwarte Jena, FSU Jena, Jena, Germany Received 2016 May 4; accepted 2016 June 8 Abstract Model atmospheres of the close visual binary star COU1511 (HIP12552) are constructed using grids of Kuruz s blanketed models to build the individual synthetic SEDs for both components. These synthetic SED s are combined together for the entire system and compared with the observational one following Al-Wardat s complex method for analyzing close visual binary stars. The entire observational spectral energy distribution (SED) of the system is used as a reference for comparison between synthetic SED and the observed one. The parameters of both components are derived as:teff a = 6180±50 K, Teff b = 5865±70 K, logg a = 4.35±0.12, logg b = 4.45±0.14,R a = 1.262±0.08R, R b = 1.006±0.07R, L a = 2.09±0.10L, L b = 1.08±0.12L, with spectral types F8 and G1 for both components (a,b) respectively, and age of 3.0 ± 0.9 Gy. A modified orbit of the system is built and the masses of the two components are calculated as M a = 1.17±0.11M,M b = 1.06±0.10M. Key words: binaries; visual stars; fundamental parameters stars; COU1511 (Hip12552) 1 INTRODUCTION Recent surveys of the sky showed that more than 50% of the galactic stellar systems are binaries, which raises their importance in understanding the formation and evolution of the galaxy. This role in determining precisely different stellar parameters gives the study of binary stars a special importance. The case is a bit complicated in the case of close visual binary stars (CVBS), which are not resolved as binaries by inspecting limited images, but can be resolved in space based observations, or by using modern techniques of ground-based observations, like speckle interferometry and adaptive optics. In addition to that, Docobo et al. (2001) pointed out that the study of orbital motion of visual and interferometric pairs remains an important astronomical discipline. The visual binaries are the main key source of information about stellar masses and distances, and they define practically our understanding of stellar physical properties especially for the lower part of the main sequence stars. For now, hundreds of CVBS with periods on the order of 10 years or less are routinely observed by different groups of high resolution techniques around the world. This has helped in determining the

2 2 Al-Wardat et al. orbital parameters and magnitude differences for some of these CVBS. However, this is not sufficient to determine the individual physical parameters of the components of the system. Al-Wardat s method for analyzing CVBS (Al-Wardat 2012) offers a complementary solution for this problem by implementing differential photometry, spectrophotometry, atmospheric modeling, and orbital solution in accurate determination of different physical and geometrical parameters of this category of stars. The method was successfully applied to several solar type and subgaint binary systems such as Hip70973, Hip72479 (Al-Wardat 2012), Hip689 (Al-Wardat et al. 2014), and Hip11253 (Al- Wardat & Widyan 2009). As a consequence of the previous work, this paper (the XI in its series) presents the analysis of the nearby solar type CVBS COU1511 (Hip12552), with a modification to its parallax. Table 1 shows the basic data of the system taken from SIMBAD and NACA/IPAC catalogues, and Table 2 shows data from Hipparcos and Tycho catalogues (ESA 1997), while Table 3 shows the magnitude difference of the system along with filters used to observe them. Table 1 Basic data of the system Hip12552 Reference α h 41 m 28 ṣ 88 SIMBAD δ Tyc HD Sp. Typ. G0 - E(B V) 0.076±0.002 NASA/IPAC A v 0 ṃ 24 NASA/IPAC Table 2 Data of Hipparcos and Tycho Catalogues Hip12552 Source of data V J(Hip) 8 ṃ 51 Hipparcos (B V) J(Hip) 0 ṃ 60± π Hip (mas) 9.69± B T 9 ṃ 24±0.016 Tycho V T 8 ṃ 59± π Tyc (mas) 14.7± π Hip (mas) 11.07±1.07 New Hipparcos (van Leeuwen 2007) Table 3 Magnitude difference between the components of the system Hip12552, along with the filters used to obtain the observations. m Filter (λ/ λ) References 0 ṃ 65± nm/ ṃ 75± nm/ ṃ 88± 562nm/ ṃ 86± 692nm/ Balega et al. (2007), 2 Docobo et al. (2006), 3 Horch et al. (2011).

3 Parameters of CVBS XI: COU1511 (HIP12552) 3 2 ATMOSPHERIC MODELING The observational SED of the system Hip12552 obtained by Al-Wardat (2002) is used as reference for the comparison with synthetic SED. Using m V = 8 ṃ 51 (see Table 2), m= 0 ṃ 76±0.03 (as the average of all m using the different filters for V-band only ( nm), see Table 3), and Hipparcos trigonometric parallax (π = ± 1.07 mas), the individual and absolute magnitudes of both components (a, b) of the system are calculated using the following relations: F a F b = m (1) M v = m v +5 5log(d) A v (2) to get m a v = 8ṃ 95±0.02, m b v = 9ṃ 71±0.05, and M a v = 4ṃ 07±0.18, M b v = 4ṃ 83±0.19, where the extinction valuea v was taken from Table 1. To calculate the preliminary input parameters used to build the atmospheric modelling, we use the bolometric magnitudes, the luminosities from Lang (1992), and Gray (2005) with the following relations: log(r/r ) = 0.5log(L/L ) 2log(T eff /T ), (3) logg = log(m/m ) 2log(R/R )+4.43, (4) to estimate the effective temperatures and gravity acceleration. These values for the effective temperature and gravity acceleration allow us to construct the model atmosphere for each component using grids of Kuruz s line blanketed models (ATLAS9) (Kurucz 1994). Here we usedt = 5777K andm bol = 4ṃ 75 in all calculations. The total energy flux from a binary star is due to the net luminosity of the components a, and b located at a distance d from the Earth. The total energy flux may be written as: Rearranging equ. 5 gives F λ d 2 = H a λ R2 a +Hb λ R2 b, (5) F λ = (R a /d) 2 (H a λ +H b λ (R b /R a ) 2 ), (6) where R a and R b are the radii of the primary and secondary components of the system in solar units, H a λ andhb λ are the fluxes at the surface of the star andf λ is the flux for the entire SED of the system. Many attempts were made to achieve the best fit (Fig. 1) between the observed flux of Al-Wardat (2002) and the total computed one using the iteration method of different sets of parameters. The best fit is found using the following set of parameters: T a eff = 6180±50K,Tb eff = 5865±70K, logg a = 4.35±0.12,logg b = 4.45±0.14, R a = 1.262±0.08R,R b = 1.006±0.07R Using equ. 3 the luminosities are calculated yielding the following values: L a = 2.09 ± 0.10L and L b = 1.08 ± 0.12L. Using tables of Gray (2005) or the Sp T eff. empirical relation from Lang (1992), the spectral types of the components (a, b) of the system are F8 and G2 respectively.

4 4 Al-Wardat et al Log F (erg/cm 2.s.Å) Component A Component B Entire Flux (Synthetic) Entire Flux (Observed) Wavelength (Å) Fig. 1 Best fit between the entire observed spectrum (dotted line) (Al-Wardat 2002) and the synthetic entire SED (solid line) for the system Hip12552 using the following parameters Teff A = 6180±50 K, logg A = 4.35±0.10,R A = 1.262±0.09R,Teff B = 5865±70 K, log g B = 4.45±0.15,R B = 1.006±0.10R with d = 84.53±0.009 pc (π = mas). 3 ORBITAL ANALYSIS The orbit of the system is built using the positional measurements listed in Table 4, following Tokovinin s method (Tokovinin 1992). The modified orbital elements of the system along with those taken from the sixth interferometric catalogue are listed in Table 5. The table shows a good agreement between our estimated orbital period, P; eccentricity, e; semi-major axis, a; inclination, i; argument of periastron, ω; position angle of nodes, Ω; and time of primary minimum, T 0 with those previously reported results. 4 MASSES Using the estimated orbital element, the masses of the system and the corresponding errors are calculated using the following relations: M A +M B = ( a3 π 3 P 2) M (7) σ M M = (3 σ π π )2 +(3 σ a a )2 +(2 σ p p )2 (8)

5 Parameters of CVBS XI: COU1511 (HIP12552) 5 Table 4 Positional measurements of the system from the Fourth Interferometric Catalogue. Epoch θ,deg ρ,arcsec References * * * * * * * These points were modified by180 to achieve consistency with nearby points. 1 McAlister & Hendry (1982), 2 McAlister et al. (1987), 3 Hartkopf et al. (1994), 4 Balega et al. (1994), 5 ten Brummelaar et al. (2000), 6 Balega et al. (1999), 7 Hartkopf et al. (1997), 8 Balega et al. (2006), 9 Balega et al. (2013), 10 Hartkopf et al. (2008), 11 Balega et al. (2007), 12 Docobo et al. (2006), 13 Mason et al. (2011), 14 Gili & Prieur (2012), 15 Horch et al. (2011). The preliminary result obtained using the new Hipparcos trigonometric parallax (π = ± 1.07 mas) (van Leeuwen 2007) is M a + M b = 2.72 ± 0.75M, while it is 3.26 ± 1.18M when using Hipparcos trigonometric parallax (π = ± 0.2 mas, See Table 5). But depending on our analysis (Sec. 2), we achieved the best fit between the synthetic and observational entire SED using π = mas, this new parallax value gives a mass sum of M a + M b = 2.23 ± 0.57M, which fits better the positions of the two components on the evolutionary tracks as shown in Fig. 3. Table 5 Orbital elements of the system. Parameters This work Hartkopf & Mason (2001) Couteau (1996) P, yr ± T 0, yr ± e ± a, arcsec ± i, deg ± ω, deg ± Ω, deg ± M a +M b,m 2.72± ± ±1.18 Based on new Hipparcos trigonometric parallax (π = 11.07±1.07 mas) Based on the parallax estimated in this work (π = mas) Based on Hipparcos trigonometric parallax (π = 10.42±0.2 mas)

6 6 Al-Wardat et al. Fig. 2 Relative visual orbit of the system HIP12552 showing the epoch of the positional measurements. 5 SYNTHETIC PHOTOMETRY As a double-check for the best fit and to present a new synthetic photometrical data of the unseen individual components of the system, we apply the following relation (Maíz Apellániz 2006, 2007): Pp (λ)f λ,s (λ)λdλ m p [F λ,s (λ)] = 2.5log Pp (λ)f λ,r (λ)λdλ +ZP p, (9) to calculate total and individual synthetic magnitudes of the systems, where m p is the synthetic magnitude of the passbandp,p p (λ) is the dimensionless sensitivity function of the passbandp,f λ,s (λ) is the synthetic SED of the object and F λ,r (λ) is the SED of the reference star (Vega). Here the zero points (ZP p ) of Maíz Apellániz (2007) are adopted. Calculated synthetic magnitudes and color indices of the entire system and individual components of different photometrical systems are shown in Table 6. 6 RESULTS AND DISCUSSION The synthetic SEDs of the individual component and the system HIP12552 are built using atmospheric modeling and the visual magnitude difference between the two components along with the total observed SED. Least square fitting with weights inversely proportional to the squares of the positional measurement errors is used to modify the orbit of the system. So, the physical and geometrical parameters of HIP12552 are estimated. Fig. 1 shows the best fit of the total synthetic SED to the observed one.

7 Parameters of CVBS XI: COU1511 (HIP12552) 7 Table 6 Synthetic magnitudes and color indices of the system. Sys. Filter Entire Comp. Comp. σ = ±0.03 a b Joh- U Cou. B V R U B B V V R Ström. u v b y u v v b b y Tycho B T V T B T V T Table 7 Comparison between the observational and synthetic magnitudes, colors and magnitude differences of the system Hip Observed a Synthetic b (This work) V J 8 ṃ 51 8 ṃ 51±0.03 B T 9 ṃ 24± ṃ 25±0.03 V T 8 ṃ 59± ṃ 58±0.03 (B V) J 0 ṃ 60± ṃ 60±0.03 m 0 ṃ 76 c ± ṃ 76±0.04 a See Table 2 b See Table 6 c As the average of all m using the differnt filters under V-band (see Table 3). Table 8 Parameters of the components of the system HIP Parameters Comp. a Comp. b T eff (K) 6180± ±70 Radius (R ) 1.262± ±0.07 logg 4.35± ±0.14 L(L ) 2.09± ±0.12 M bol 3 ṃ 95± ṃ 67±0.19 M V 4 ṃ 07± ṃ 83±0.19 Mass (M ) 1.17± ±0.10 Sp. Type F8 G2 Parallax (mas) ± 1.07 ( Ma+M b M ) 2.23±0.57 depending on the evolutionary tracks (Fig. 3), depending on the tables of Lang (1992), depending on the orbital solution. Table 7 shows a comparison between the observational and synthetic magnitudes, colors and magnitudes differences for the system HIP This gives a good indication of the reliability of the estimated parameters of the individual components of the system which are listed in Table 8.

8 8 Al-Wardat et al. 0.8 Log (L/L ) M 1.3M 1.2M Main Sequence HIP Component A Component B 1.1M 1M Log T eff Fig. 3 The system components on the evolutionary tracks of Girardi et al. (2000b). The positions of the system s components on the evolutionary tracks of Girardi et al. (2000a) (Fig. 3) show that both components with massesm A = 1.20 andm B = 1.09M belong to the main-sequence stars. And their positions on Girardi et al. (2000a) isochrones for low- and intermediate-mass stars of different metallicities and that of the solar composition [Z = 0.019, Y = 0.273] are shown in Figs 4 & 5, which give an age of the system around3.0±0.9 Gy. 7 CONCLUSION The CVBS COU1511 (HIP12552) is analyzed using Al-Wardat s complex method for analyzing close visual binary stars, which is based on combining magnitude difference measurements from speckle interferometry, entire spectral energy distribution (SED) from spectrophotometry, atmospheres modeling and orbital analysis to estimate the individual physical and geometrical parameters of the system. The entire and individual Johnson-Cousin UBVR, Strömgren uvby, and Tycho BV synthetic magnitudes and color indices of the system are calculated. A modified orbit and geometrical elements of the system are introduced and compared with earlier results. The positions of the two components on the evolutionary tracks and isochrones are shown, their spectral types are estimated as F8 and G1 respectively with the age of3.0±0.9 Gy. Acknowledgements This work made use of SAO/ NASA, SIMBAD database, Fourth Catalog of Interferometric Measurements of Binary Stars, IPAC data systems and CHORIZOS code of photometric and spectrophotometric data analysis. The authors thank Mrs. Kawther Al-Waqfi for her help in some orbital calculations.

9 Parameters of CVBS XI: COU1511 (HIP12552) Gy Gy Log L/L HIP Component A Component B Log T eff Fig. 4 The systems components on the isochrones of low- and intermediate-mass, solar composition [Z=0.019,Y=0.273] stars of Girardi et al. (2000a). References Al-Wardat, M. A. 2002, Bull. Special Astrophys. Obs., 53, 58 3, 4 Al-Wardat, M. A. 2012, PASA, 29, Al-Wardat, M. A., Balega, Y. Y., Leushin, V. V., et al. 2014, Astrophysical Bulletin, 69, 58 2 Al-Wardat, M. A., & Widyan, H. 2009, Astrophysical Bulletin, 64, Balega, I. I., Balega, A. F., Maksimov, E. V., et al. 2006, Bull. Special Astrophys. Obs., 59, 20 5 Balega, I. I., Balega, Y. Y., Belkin, I. N., et al. 1994, A&AS, 105, Balega, I. I., Balega, Y. Y., Gasanova, L. T., et al. 2013, Astrophysical Bulletin, 68, 53 5 Balega, I. I., Balega, Y. Y., Maksimov, A. F., et al. 2007, Astrophysical Bulletin, 62, 339 2, 5 Balega, I. I., Balega, Y. Y., Maksimov, A. F., et al. 1999, A&AS, 140, Couteau, P. 1996, IAU Commission on Double Stars, 128, 1 5 Docobo, J. A., Tamazian, V. S., Balega, Y. Y., et al. 2001, A&A, 366, Docobo, J. A., Tamazian, V. S., Balega, Y. Y., & Melikian, N. D. 2006, AJ, 132, 994 2, 5 ESA. 1997, The Hipparcos and Tycho Catalogues (ESA) 2 Gili, R., & Prieur, J.-L. 2012, Astronomische Nachrichten, 333, Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000a, A&AS, 141, 371 8, 9, 10 Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000b, VizieR Online Data Catalog, 414, Gray, D. F. 2005, The Observation and Analysis of Stellar Photospheres, 505 3

10 10 Al-Wardat et al Log L/L HIP Component A Component B Log T eff Fig. 5 The systems components on the isochrones for low- and intermediate-mass stars of different metallicities of Girardi et al. (2000a). Hartkopf, W. I., & Mason, B. D. 2001, IAU Commission on Double Stars, 145, 1 5 Hartkopf, W. I., Mason, B. D., & Rafferty, T. J. 2008, AJ, 135, Hartkopf, W. I., McAlister, H. A., Mason, B. D., et al. 1994, AJ, 108, Hartkopf, W. I., McAlister, H. A., Mason, B. D., et al. 1997, AJ, 114, Horch, E. P., Gomez, S. C., Sherry, W. H., et al. 2011, AJ, 141, 45 2, 5 Kurucz, R. 1994, Solar abundance model atmospheres for 0,1,2,4,8 km/s. Kurucz CD-ROM No. 19. Cambridge, Mass.: Smithsonian Astrophysical Observatory, 1994., 19 3 Lang, K. R. 1992, Astrophysical Data I. Planets and Stars., 133 3, 7 Maíz Apellániz, J. 2006, AJ, 131, Maíz Apellániz, J. 2007, in Astronomical Society of the Pacific Conference Series, Vol. 364, The Future of Photometric, Spectrophotometric and Polarimetric Standardization, ed. C. Sterken (San Francisco: Astronomical Society of the Pacific), Mason, B. D., Hartkopf, W. I., Raghavan, D., et al. 2011, AJ, 142, McAlister, H. A., Hartkopf, W. I., Hutter, D. J., & Franz, O. G. 1987, AJ, 93, McAlister, H. A., & Hendry, E. M. 1982, ApJS, 49, ten Brummelaar, T., Mason, B. D., McAlister, H. A., et al. 2000, AJ, 119, Tokovinin, A. 1992, in Astronomical Society of the Pacific Conference Series, Vol. 32, IAU Colloq. 135: Complementary Approaches to Double and Multiple Star Research, ed. H. A. McAlister & W. I. Hartkopf, 573 4

11 van Leeuwen, F. 2007, A&A, 474, 653 2, 5 Parameters of CVBS XI: COU1511 (HIP12552) 11

Physical and geometrical parameters of CVBS XI: Cou 1511 (HIP 12552)

Research in Astronomy and Astrophysics PAPER Physical and geometrical parameters of CVBS XI: Cou 1511 (HIP 12552) Recent citations - Stellar parameters of the two binary systems: HIP 14075 and HIP 14230