KIC 7385478: An eclipsing binary with a {\gamma} Doradus component

TL;DR

This study analyzes the eclipsing binary KIC 7385478, revealing a primary ba0Doradus pulsator with ba0Doradus pulsations and secondary star activity, combining spectroscopic and Kepler photometric data.

Contribution

It provides the first detailed spectroscopic and photometric analysis of KIC 7385478, identifying ba0Doradus pulsations and characterizing the system's stellar parameters.

Findings

Primary component is a ba0Doradus pulsator.

Detected ba0Doradus pulsations with a period of ba00.5 days.

Secondary star activity influences residual light variations.

Abstract

We present spectroscopic and photometric analysis of the eclipsing binary KIC\,7385478. We find that the system is formed by F1V + K4III-IV components. Combining results from analysis of spectroscopic data and $Kepler$ photometry, we calculate masses and radii of the primary and the secondary components as M$_{1}$ = 1.71 $\pm$ 0.08 \Msun, M$_{2}$ = 0.37 $\pm$ 0.04 \Msun~ and R$_{1}$ = 1.59 $\pm$ 0.03\Rsun, R$_{2}$ = 1.90 $\pm$ 0.03\Rsun, respectively. Position of the primary component in HR diagram is in the region of $\gamma$ Doradus type pulsators and residuals from light curve modeling exhibit additional light variation with a dominant period of $\sim$0.5 day. These are clear evidences of the $\gamma$ Doradus type pulsations on the primary component. We also observe occasional increase in amplitude of the residuals, where the orbital period becomes the most dominant period. These may be attributed to the cool star activity originating from the secondary component.

Tables

Table 1. Table 1: Log of spectroscopic observations together with measured radial velocities and their corresponding standard errors ( σ 𝜎 \sigma ) in km s -1 .

HJD	Orbital	Exposure	Primary		Secondary
(24 00000+)	Phase	time (s)	Vr	$\sigma$	Vr	$\sigma$
56845.4821	0.2399	3600	-54.5	5.2	159.0	19.9
56845.5248	0.2657	3600	-46.8	7.0	158.2	23.9
56846.4677	0.8352	3600	10.2	6.3	-177.1	22.5
56887.4623	0.5984	3200	8.1	6.8	-128.1	38.6
57591.5291	0.8948	3600	11.0	10.2	-136.2	41.6
57600.5477	0.3425	3600	-46.8	10.9	132.6	32.4
57601.3798	0.8452	3600	24.3	13.2	-152.0	38.1
57617.2875	0.4543	3600	-27.3	4.7	—	—
57672.3540	0.7174	3600	24.4	8.4	-191.5	25.3

Table 2. Table 2: Spectroscopic orbital elements of KIC 7385478. M 1 M{{}_{1}} and M 2 M{{}_{2}} denote the masses of the primary and secondary component, respectively, while M 𝑀 M shows the total mass of the system.

Parameter	Value
$P_{\mathrm{orb}}$ (days)	1.655473 (fixed)
$T_0$ (HJD2454+)	954.534784 (fixed)
$\gamma$ (km s-1)	-16.2$\pm$0.8
$K_1$ (km s-1)	38.3$\pm$2.7
$K_2$ (km s-1)	178.8$\pm$2.7
$e$	0 (fixed)
$a\sin i$ (R⊙)	7.10$\pm$0.12
$M{\sin }^{3}i$ (M⊙)	1.754$\pm$0.061
Mass ratio ($q=M{}_2{}^{}/M{}_1$)	0.21$\pm$0.02
fit rms (km s-1)	6.5

Table 3. Table 3: Light curve modeling results of KIC 7385478. ⟨ r 1 ⟩ delimited-⟨⟩ subscript 𝑟 1 \langle r_{1}\rangle and ⟨ r 2 ⟩ delimited-⟨⟩ subscript 𝑟 2 \langle r_{2}\rangle denote mean fractional radii of the primary and the secondary components, respectively. Spot parameters are given in the order of co-latitude ( θ 𝜃 \theta ), longitude ( φ 𝜑 \varphi ), radius ( r s ​ p ​ o ​ t subscript 𝑟 𝑠 𝑝 𝑜 𝑡 r_{spot} ) and temperature factor ( T ​ F 𝑇 𝐹 TF ). Internal errors of the adjusted parameters are given in parentheses for the last digits. Asterix symbols in the table denote fixed value for the corresponding parameters.

Parameter	Value
$q$	0.21*
$T_1\left(K\right)$	7000*
$g_1$, $g_2$	0.32*, 0.270(5)
$A_1$, $A_2$	0.5*, 0.699(3)
$F_1$ = $F_2$	1.0*
phase shift	0.00160(2)
$i{(}^{\circ })$	70.966(8)
$T_2\left(K\right)$	4293(150)
${\Omega }_1$	4.9582(43)
${\Omega }_2$	2.2574*
$L_1$/($L_1$+$L_2)$	0.8590(6)
$x{}_1{}^{}bol,x{}_2{}^{}bol$	0.471*, 0.531*
$x{}_1{}^{},x{}_2$	0.462*, 0.344(8)
$⟨r_1⟩,⟨r_2⟩$	0.2114(2), 0.2527*
$\theta {(}^{\circ })$	98(6)
$\phi {(}^{\circ })$	93(1)
$r_{spot}{(}^{\circ })$	13(1)
$TF$	1.022(2)
Model rms	1.1 $\times$ 10-4

Table 4. Table 4: Absolute physical properties of KIC 7385478. Error of each parameter is given in paranthesis for the last digits.

Parameter	Primary	Secondary
Spectral Type	F1V	K4 III-IV
Mass (M⊙)	1.71(8)	0.37(4)
Radius (R⊙)	1.59(3)	1.90(3)
Log $L/L_{\odot }$	0.737(52)	0.043(82)
log $g$ (cgs)	4.269(7)	3.444(31)
$M_{bol}$ (mag)	2.91(13)	4.64(21)

Table 5. Table 5: Extracted frequencies in multi-frequency analysis. N, F, A, P and SNR means number, frequency (in c/d), amplitude (in mmag), phase and signal – to –noise ratio, respectively.

N	F	A	P	SNR
1	2.0252	0.470	0.2618	275.8
2	1.9428	0.315	0.6111	210.5
3	0.6035	0.183	0.7206	105.5
4	0.6023	0.169	0.5600	84.7
5	0.6563	0.143	0.5643	76.2
6	0.6055	0.127	0.1400	71.5
7	1.8052	0.124	0.0554	91.7
8	0.7138	0.115	0.7896	69.5
9	2.4208	0.102	0.4827	87.2
10	0.6048	0.090	0.3848	62.7
11	0.6062	0.079	0.5564	50.9
12	1.2087	0.058	0.1224	50.2
13	1.2119	0.064	0.5436	48.1
14	0.6008	0.065	0.7873	42.2
15	0.4780	0.061	0.4596	42.2
16	2.1178	0.060	0.4606	54.5
17	3.7537	0.058	0.3023	67.9
18	0.6110	0.062	0.5936	38.9
19	0.6002	0.054	0.2745	39.7
20	0.6084	0.045	0.1297	38.7
21	0.6028	0.050	0.7349	36.8
22	1.2041	0.051	0.3854	35.8
23	0.5983	0.034	0.3999	29.4
24	1.2051	0.036	0.5672	32.1
25	0.8446	0.036	0.8880	27.6
26	1.2110	0.036	0.6077	28.4
27	1.8109	0.033	0.6449	29.9
28	0.6041	0.035	0.0754	25.2
29	0.6075	0.034	0.5545	26.5
30	1.2099	0.032	0.5354	27.3
31	0.5956	0.030	0.6500	22.3
32	1.1569	0.026	0.9567	23.0
33	0.5993	0.030	0.3687	23.2
34	1.6469	0.024	0.7514	23.8
35	0.5918	0.025	0.1310	19.6
36	0.5970	0.021	0.9181	20.0
37	0.0008	0.024	0.8344	18.6
38	1.4211	0.022	0.8007	20.5
39	2.8693	0.022	0.3247	24.7
40	2.6815	0.022	0.0867	23.9
41	0.6091	0.026	0.5645	18.4
42	0.6101	0.022	0.8878	20.5
43	1.2078	0.020	0.8653	20.3
44	0.0129	0.021	0.2401	15.2
45	1.8117	0.020	0.2413	20.2
46	1.3322	0.019	0.7038	18.6
47	0.6197	0.021	0.8756	16.6
48	3.8374	0.020	0.9090	25.4
49	1.1790	0.018	0.4249	17.4
50	0.6662	0.017	0.4237	15.9
51	3.0771	0.019	0.1901	22.3
52	0.5902	0.016	0.0661	15.7
53	0.6144	0.017	0.9196	15.4
54	0.6162	0.016	0.5804	15.6
55	1.2030	0.018	0.0640	17.3
56	0.0441	0.018	0.9902	13.9
57	0.0103	0.025	0.7827	13.7
58	2.4687	0.018	0.5315	18.9
59	4.4414	0.017	0.7160	24.4
60	0.8171	0.018	0.0733	15.3
61	2.7390	0.018	0.4742	19.9
62	0.6216	0.018	0.4458	15.1
63	1.2849	0.017	0.1149	16.3
64	4.6194	0.016	0.5229	23.7
65	2.3371	0.016	0.3393	17.2
66	1.1995	0.018	0.3142	14.7
67	1.2005	0.018	0.6592	17.3
68	2.6293	0.014	0.3300	17.5
69	0.4964	0.013	0.1267	13.2
70	4.8330	0.013	0.6804	22.5
71	3.0200	0.014	0.3884	18.4
72	0.0309	0.014	0.8959	11.9
73	1.2133	0.014	0.6766	14.5
74	0.6301	0.011	0.8103	13.1
75	3.1345	0.014	0.0059	18.4

Table 6. Table 5: Continued.

N	F	A	P	SNR
226	0.0354	0.006	0.4586	7.6
227	0.0938	0.007	0.1194	7.1
228	1.1982	0.007	0.5649	8.6
229	0.9672	0.006	0.0567	8.3
230	2.8896	0.007	0.9090	10.8
231	1.2020	0.009	0.0116	9.4
232	3.1509	0.007	0.9775	11.2
233	0.5869	0.009	0.8703	8.1
234	1.2355	0.007	0.5330	8.6
235	3.0208	0.007	0.4495	11.1
236	2.7926	0.006	0.8003	10.5
237	0.3717	0.007	0.5017	7.2
238	1.2140	0.008	0.0111	8.7
239	0.0222	0.009	0.8167	6.8
240	0.5759	0.007	0.5043	7.7
241	0.5772	0.006	0.1883	8.2
242	0.0739	0.006	0.3209	6.7
243	0.0574	0.005	0.2890	7.2
244	0.0077	0.006	0.9595	7.0
245	4.4674	0.006	0.2380	12.8
246	0.9993	0.006	0.4145	8.1
247	0.0322	0.007	0.6979	6.5
248	1.1968	0.006	0.3491	8.2
249	0.6349	0.004	0.1773	7.5
250	0.0787	0.006	0.2070	6.8
251	4.6394	0.006	0.1738	12.8
252	2.0074	0.006	0.1696	9.1
253	0.6599	0.006	0.8010	7.3
254	1.2186	0.006	0.7455	8.1
255	2.4178	0.007	0.2341	10.0
256	0.1532	0.006	0.5155	6.4
257	3.6261	0.008	0.0635	11.5
258	1.8128	0.006	0.4256	9.2
259	0.0178	0.009	0.9468	7.1
260	0.1111	0.007	0.2124	6.6
261	0.0609	0.009	0.2394	6.6
262	0.1258	0.009	0.3570	6.6
263	1.2226	0.008	0.0595	8.1
264	0.6411	0.006	0.8622	7.2
265	0.6323	0.008	0.4660	7.4
266	3.0235	0.007	0.6443	10.0
267	4.8395	0.006	0.7207	12.3
268	3.0726	0.005	0.2607	9.9
269	0.0494	0.007	0.2219	6.7
270	0.0755	0.009	0.6666	6.8
271	0.5363	0.004	0.8838	6.9
272	1.0775	0.006	0.5053	7.6
273	4.2214	0.006	0.7295	11.4
274	1.9485	0.006	0.3600	8.5
275	2.0246	0.007	0.7458	8.6
276	0.5568	0.007	0.1948	6.7
277	0.5676	0.005	0.3004	6.8
278	0.4949	0.006	0.3021	6.7
279	4.8282	0.006	0.0358	11.8
280	4.8343	0.007	0.7484	12.1
281	0.1558	0.006	0.2316	6.2
282	0.1548	0.007	0.8164	7.1
283	0.1661	0.006	0.3987	6.2
284	0.5587	0.006	0.9196	7.1
285	0.0972	0.005	0.5152	6.0
286	0.0644	0.005	0.9109	6.2
287	1.9416	0.006	0.6954	8.5
288	1.1844	0.006	0.0778	7.4
289	1.2148	0.007	0.6590	7.7
290	2.4615	0.005	0.7649	8.7
291	1.1953	0.005	0.1796	7.5
292	0.1286	0.006	0.4343	6.2
293	2.4120	0.006	0.6780	8.6
294	2.4219	0.006	0.3595	9.1
295	0.6379	0.004	0.7334	6.8
296	4.2270	0.006	0.2032	10.9
297	0.1249	0.005	0.5003	6.2
298	4.8304	0.005	0.5935	11.5
299	3.6203	0.006	0.6814	10.0
300	1.3331	0.006	0.5969	7.6

Table 7. Table 5: Continued.

N	F	A	P	SNR
451	0.2924	0.005	0.9228	5.1
452	4.2316	0.005	0.9844	8.8
453	4.2304	0.004	0.4969	9.4
454	4.3576	0.004	0.8228	8.8
455	1.8041	0.004	0.9325	6.7
456	2.2187	0.004	0.9185	6.8
457	4.8310	0.004	0.1935	9.3
458	1.2346	0.004	0.5441	6.1
459	1.8602	0.004	0.0884	6.3
460	3.5830	0.004	0.8596	7.7
461	0.6746	0.004	0.0340	5.4
462	0.5685	0.005	0.8404	6.0
463	0.5419	0.004	0.7589	5.6
464	0.2362	0.004	0.7916	5.0
465	0.1931	0.003	0.8427	5.0
466	3.5256	0.004	0.5599	7.6
467	1.1938	0.004	0.1118	5.9
468	0.0705	0.007	0.8385	5.6
469	0.6673	0.006	0.7820	5.6
470	0.6685	0.005	0.2409	5.8
471	0.4626	0.004	0.8015	5.1
472	4.8361	0.004	0.9403	9.3
473	4.9617	0.004	0.5742	8.9
474	0.0945	0.006	0.3966	5.3
475	1.9475	0.004	0.4560	6.5
476	0.6667	0.006	0.1897	6.1
477	0.6985	0.004	0.1118	5.4
478	0.6803	0.004	0.6363	5.4
479	0.0303	0.006	0.8319	5.3
480	0.7038	0.004	0.1086	5.3
481	0.3299	0.004	0.7127	4.9
482	0.3033	0.004	0.6330	5.0
483	0.1611	0.004	0.0336	4.9
484	0.2035	0.004	0.1517	4.8
485	0.2694	0.004	0.9749	4.8
486	3.9009	0.004	0.1089	7.7
487	0.1649	0.004	0.5565	4.9
488	1.1698	0.003	0.7949	5.4
489	1.2301	0.004	0.3618	5.6
490	1.5136	0.003	0.8284	5.6
491	0.5192	0.004	0.6906	4.9
492	0.5221	0.003	0.0884	5.0
493	0.3680	0.003	0.2712	4.8
494	0.4213	0.003	0.8097	5.0
495	0.1854	0.004	0.4820	4.7
496	0.4602	0.003	0.5508	5.0
497	0.2642	0.004	0.2734	4.6
498	0.0652	0.008	0.0828	5.4
499	0.0732	0.004	0.9467	4.8
500	0.0555	0.007	0.3539	5.0
501	0.1031	0.004	0.4012	5.0
502	0.1074	0.005	0.6009	4.7
503	4.0991	0.003	0.6757	7.6
504	0.5535	0.003	0.4260	5.0
505	2.8701	0.003	0.3971	6.6
506	1.9436	0.004	0.2756	6.1
507	4.8374	0.004	0.7261	9.0
508	4.0504	0.003	0.7416	7.6
509	3.0118	0.003	0.3390	6.7
510	0.5441	0.003	0.0740	5.0
511	0.4287	0.004	0.6468	4.8
512	0.7063	0.003	0.3598	5.0
513	2.4231	0.004	0.6310	6.4
514	1.2313	0.004	0.4308	5.5
515	1.8756	0.003	0.3359	5.8
516	1.1806	0.004	0.1709	5.4
517	0.1280	0.005	0.5063	4.8
518	0.1297	0.005	0.3037	5.5
519	0.1899	0.004	0.9423	4.8
520	0.0457	0.005	0.1781	5.0
521	2.0281	0.003	0.1735	5.9
522	0.1244	0.005	0.9761	4.8
523	0.1475	0.004	0.1391	4.5
524	0.1595	0.004	0.5191	4.8
525	0.5283	0.004	0.5056	5.0

Table 8. Table 5: Continued.

N	F	A	P	SNR
676	0.3005	0.003	0.5771	4.2
677	1.2216	0.003	0.0587	5.2
678	3.9971	0.003	0.4223	6.1
679	1.8259	0.003	0.3384	5.0
680	1.8076	0.003	0.8685	5.2
681	2.8270	0.003	0.5419	5.4
682	0.4640	0.002	0.4254	4.3
683	0.4567	0.003	0.5376	4.4
684	0.0664	0.004	0.5410	4.8
685	4.8297	0.003	0.6131	6.8
686	0.0549	0.004	0.5179	4.8
687	0.0425	0.004	0.6134	4.5
688	0.3584	0.003	0.0050	4.3
689	0.1907	0.003	0.9363	4.6
690	0.1162	0.003	0.8436	4.8
691	0.3911	0.003	0.3221	4.3
692	0.6542	0.003	0.7953	4.6
693	0.1395	0.003	0.6188	4.3
694	0.2577	0.003	0.0831	4.2
695	0.1088	0.003	0.7113	4.6
696	2.0526	0.003	0.3691	5.1
697	1.9362	0.003	0.8072	5.0
698	1.3772	0.003	0.8655	4.7
699	2.4126	0.003	0.1311	5.4
700	0.6570	0.003	0.2109	4.9
701	0.5151	0.003	0.4084	4.4
702	0.8439	0.003	0.7323	4.3
703	1.0644	0.003	0.8552	4.5
704	1.8472	0.003	0.5657	4.9
705	4.2238	0.003	0.9221	6.5
706	0.9428	0.003	0.3426	4.3
707	0.6134	0.004	0.3359	5.1
708	0.5909	0.003	0.1009	4.9
709	0.3504	0.003	0.7099	4.0
710	0.3447	0.003	0.9244	4.2
711	0.6714	0.004	0.3836	4.5
712	0.6435	0.003	0.1370	5.3
713	0.5451	0.003	0.1752	4.8
714	0.6725	0.003	0.7879	4.7
715	0.6979	0.003	0.8107	4.7
716	0.7168	0.003	0.7007	4.4
717	0.6623	0.003	0.1877	4.9
718	3.9300	0.003	0.8226	5.9
719	1.2402	0.003	0.0255	4.7
720	1.2586	0.003	0.6966	4.7
721	0.5260	0.003	0.7856	4.5
722	2.8341	0.002	0.3379	5.2
723	0.2627	0.003	0.1238	4.1
724	0.5121	0.003	0.0645	4.4
725	0.7406	0.003	0.0371	4.5
726	0.5621	0.003	0.6052	4.7
727	0.4930	0.003	0.8863	4.3
728	1.8162	0.003	0.0778	5.1
729	0.8539	0.002	0.1619	4.3
730	1.2447	0.003	0.8448	4.7
731	3.9228	0.002	0.3068	5.7
732	1.3589	0.002	0.2463	4.5
733	3.2856	0.002	0.2462	5.3
734	3.8896	0.002	0.5791	5.7
735	0.2824	0.002	0.4703	4.0

Full text

1 Introduction

High precision optical photometry from space telescopes is a milestone for modeling of eclipsing binary systems, where current models can not reproduce very small amplitude variations adequately. Even though, reasonable light curve modeling of these systems provides not only absolute physical properties, but also reveals additional small amplitude light variations beside the eclipse events and binarity effects (i.e. reflection, gravity darkening, etc…). This leads us to find accurate position of the components of these system on HR diagram and evaluate their variation nature more accurately.

There is a few indicators to reveal the internal layers of the stars. One of them is the stellar pulsation. Unfortunately the pulsations can not be observable for each star. On the other hand, the initial analysis and some studies in the literature, such as Uytterhoeven et al., (2011), indicated that KIC 7385478 is one of the candidates for the eclipsing binaries with pulsating component. Moreover, the pulsation behaviour seen in the stars is very important pattern to understand the stellar itself and its evolution. According to the observations lasting as long as several decades indicate that there are several type pulsating stars such as Cepheid, $\gamma$ Doradus, and $\delta$ Scuti type pulsating stars in the Instability Strip in the Hertzsprung-Russell diagram, especially on the main sequence. All these types are separated by their locations in the Instability Strip from each other. Analysing the pulsation frequencies, which is generally known as stellar seismology called asteroseismology, the physical processes behind both the pulsating nature and stellar interiors can be revealed. This is why the pulsating stars have an important role to understanding stellar evolution (Cunha et al., 2007; Aerts et al., 2010).

KIC 7385478 is an interesting eclipsing binary whose out-of-eclipse variation does not seem regular. The system was identified as a variable star in ASAS catalog (ASAS J195058+4259.8) for the first time (Pigulski et al. 2009) with a variation period of 1d.6551 and $V-I$ color of 0m.661. Variation amplitude in $V$ and $I$ filters are given as 0m.12 and 0m.19, respectively. Later, Slawson et al. (2011) compiled the $Kepler$ eclipsing binary catalog111http://keplerebs.villanova.edu/, including KIC 7385478. Extracted high precision light curves, ephemeris, period and effective temperature estimation for most of the system are provided in the catalog. The catalog provides effective the temperature of KIC 7385478 as 6477 $K$. Pinsonneault et al. (2012) revised the temperature to 6735 $K$. More recently, Armstrong et al. (2014) estimated the effective temperatures of the primary and the secondary components of KIC 7385478 as 6346 $K$ and 4719 $K$, respectively, which are based on spectral energy distribution fitting. The ephemeris and the period were revised by Borkovits et al. (2016), who provided eclipse timing analysis of the system. Their analysis revealed the orbit of a third body physically bound to the eclipsing binary and they estimated the minimum mass of the third body as 0.27M*⊙*. Beyond these studies, there is no comprehensive analysis of the system published so far.

In this study, we focus on spectroscopic and photometric modeling of the system. In the next section we give summary of photometric observations, and spectroscopic observations including reduction process. Section 3 comprises spectroscopic, orbital and light curve modeling together with calculated physical properties and evolutionary status of the system. In addition we further give analysis on residuals from light curve modeling in order to investigate out-of-eclipse light variations. In the final section, we summarize and discuss our findings.

2 Observations and data reductions

2.1 Kepler photometry

There is no ”$true$” photometric filters in photometer of the $Kepler$ space telescope, hence $Kepler$ photometry contains no color information. However, response function of the photometer covers a very broad wavelength range, which is between 4100 Å and 9100 Å and this allows to collect photons from a large part of the optical spectrum and increases the photometric precision, thus enables to detect light variations with a sub-milimag amplitude, such as planet transit and small amplitude oscillations. Photometric measurements are collected in short cadence and long cadence mode, which have typical exposure times of 58.89 seconds and 29.4 minutes, respectively. In this study, we use long cadence (29.4 min) data of KIC 7385478 available at $Kepler$ eclipsing binary catalog. We consider detrended and normalized fluxes (Prša et al. 2011; Slawson et al. 2011) in the catalogue. However, observations from quarter 12 and 13 are not included in the catalogue data file, hence we downloaded data files of these missing quarters from MAST222http://archive.stsci.edu/kepler/ data archive center, in fits format. In order to evaluate these data together with the catalogue data, we consider simple aperture photometry (SAP) measurements in the fits data files of quarter 12 and 13. Then we detrend and normalize SAP fluxes as described in Slawson et al. (2011). The final data set covers $\sim$4 years of time span with 65 722 data points in total. MAST archive reports 0.9% contamination level in the measurements, which indicates negligible contribution to the measured fluxes of KIC 7385478, if any.

2.2 Spectroscopy

We carried out spectroscopic observations of KIC 7385478 by 1.5 m Russian – Turkish telescope equipped with Turkish Faint Object Spectrograph Camera (TFOSC) at Tubitak National Observatory333http://www.tug.tubitak.gov.tr/rtt150_tfosc.php. Spectra were recorded on a back illuminated 2048 $\times$ 2048 pixels CCD camera with a pixel size of 15 $\times$ 15 $\mu m^{2}$. Observed spectra were obtained in échelle mode, which provides an effective resolution of R = $\lambda/\Delta\lambda$ = 2500 , which indicates $\Delta\lambda$ value of 2.6 Å around 6500 Å region. The spectra covers usable wavelength range between 3900 – 9100 Å  in 11 orders.

We obtained nine optical spectra of KIC 7385478 between 2014 and 2016 observing seasons. Exposure time of observations are 3200 s and 3600 s depending on atmospheric seeing conditions. Signal – to – noise ratio (SNR) of observed are between 80 and 145, which are estimated via photon statistic . We further obtained high SNR spectra of $\iota$ Psc (HD 222368, F7V, $v_{r}=5.656$ km s*-1*), and HD 184499 (G0V, $v_{r}=-166.1$ km s*-1*) to use as radial velocity template and spectroscopic comparison.

We used standard IRAF444The Image Reduction and Analysis Facility is hosted by the National Optical Astronomy Observatories in Tucson, Arizona at URL iraf.noao.edu. packages and tasks to reduce all observations. In each observing run, several bias and halogen lamp (flat field) frames were obtained as well as comparison lamp (Fe-Ar) images, just before or after the target star observation. In the beginning of the reduction process, master bias frame obtained from nightly taken 8-10 bias frames were subtracted from all object, Fe-Ar comparison lamp and halogen lamp frames. Then, normalized master flat field frame was produced from bias corrected halogen lamp frames and all target and Fe-Ar spectra were divided by the normalized flat field frame. Cosmic rays removal and scattered light corrections were applied to the bias and flat field corrected images. The spectra from reduced frames were extracted with IRAF task ”$apall$” under noao.imred.echelle package. Fe-Ar images were used for wavelength calibration and finally wavelength calibrated target star spectra were normalized to the unity by using cubic spline functions.

3 Analysis

3.1 Spectral type

Instead of the published effective temperature estimations in the literature, we use the advantage of having medium resolution spectra of the system to estimate the spectral type and effective temperature, which would be more accurate.

Our preliminary light curve analysis indicates that the secondary component makes $\sim$15% contribution to the total light of the system, which means the contribution of the secondary component in the spectra taken around secondary minimum would be negligible in our resolution. We have a good spectrum taken at orbital phase $\sim$0.45 (see Table 1) where the signal of the secondary component fairly diminishes in the resolution of our spectrum, so that we mainly observe the spectrum of the primary component. We adopt this spectrum as $reference~{}spectrum$ for the primary component and compare it with the observed spectra of $\iota$ Psc and HD 184499. The comparison shows that the primary seems hotter than the both standard stars, therefore, we switch to synthetic spectrum fitting method. In general, we use the latest version of python framework $iSpec$ (Blanco-Cuaresma et al. 2014), which provides easy use of various synthetic spectrum calculation codes. Among these codes, we adopt SPECTRUM555http://www.appstate.edu/$\sim$grayro/spectrum/spectrum.html code (Gray & Corbally 1994) and calculate synthetic spectra by using ATLAS-9 (Castelli & Kurucz 2004) model atmospheres in conjunction with line list from the third version of the Vienna atomic line database ($VALD3$) (Ryabchikova et al. 2015). Grids of model atmospheres are taken from the temperature range between 6000 K and 7500 K in steps of 250 K and gravity (log $g$) range between 4.5 and 3.5. In all calculations, we adopt 2 km s*-1* of microturbulence velocity and solar metallicity. Calculated spectra are convolved with a Gaussian line-spread function to match the resolution of the reference spectrum and this convolution is done either by auxiliary program $smooth$ provided by the SPECTRUM code or by built-in function of $iSpec$ code. However, we do not consider the rotational broadening and other broadening mechanisms, which could be ignored due to the relatively large instrumental broadening in our spectra. The models with the temperature 7000 K and log $g$ values 4.5 and 4.0 provide very close matches to the reference spectrum. However, calculated physical properties of the system show that the log $g$ value of the primary component is 4.27 (see Section 3.4). We also use synthetic spectrum fitting routine in $iSpec$, which mainly adopts method of minimizing $\chi^{2}$ value of the difference between the reference spectrum and synthetic spectra. For this purpose, we choose wavelength region of 4750 – 5700 Å and executed fitting routine by fixing log $g$ value and choosing different starting parameters for the temperature and metallicity. This method confirmed primary star temperature as 7000 K and solar metallicity assumption. This temperature corresponds to the $F1$ spectral type according to the calibration given by Gray (2005). Considering calculated error of the temperature from $\chi^{2}$ minimization method, temperature steps in ATLAS-9 grids and current resolution of the observed spectra, uncertainty of the temperature is estimated to be 200 K. We plot the observed reference spectrum and the best matched synthetic spectrum in Figure 1, for three different regions.

3.2 Radial velocities and spectroscopic orbit

We calculate radial velocities of the system by cross-correlating each observed spectrum with a template spectrum using $fxcor$ task in IRAF (Tonry & Davis 1979). Here, we adopt $\iota$ Psc as the template since it provides the most similar spectrum to KIC 7385478 and obtained with the same instrumental setup. We consider absorption lines in échelle orders 3, 4, 5 and 6, which cover 5000 - 6800 Å, except strongly blended lines and broad lines, such as Hα and Na I D lines. In Figure 2, we show cross-correlation functions of two spectra obtained around orbital quadratures.

We give log of observations in Table 1, together with measured radial velocities and their standard errors. Orbital phases in the table are calculated via ephemeris and period given in Borkovits et al. (2016), which we also adopt for further analysis. Investigating Table  1, one can easily notice the large standard error values for the measured radial velocities of secondary component. This is primarily caused by the small contribution of the secondary component to the total spectrum of the system, which can be evaluated from the cross correlation peaks of components given in Figure 2. We note that 3600 s of exposure time corresponds to 0.025 orbital phase step and estimated radial velocity shift during this exposure time is about 8 km s*-1*at the orbital quadratures for the secondary component, and much smaller for the primary component. These shifts are negligible compared to the standard errors of the measured radial velocities, therefore we can safely ignore the velocity shift due to the long exposure time.

Preliminary light curve examination of KIC 7385478 indicates no evidence for eccentric orbit, therefore we calculate spectroscopic orbital elements of the system by assuming circular orbit. In addition, we fixed the orbital period and ephemeris to the values given by (Borkovits et al. 2016). We use a simple python script written by us, which applies differential corrections by least squares method to all radial velocities, as described in Aitken (1935). We tabulate calculated spectroscopic orbital elements in Table 2 and plot measured radial velocities together with theoretical spectroscopic orbit and residuals from solution in Figure 3. Residuals, which belongs to the primary component scatters over zero level, while the residuals of the secondary component occupy the sub-zero level and this may be interpreted as if the whole fit could be improved further. However, if one consider the total rms of the fit given in Table 2 (last row) as one $\sigma$, then the scatter of residuals is inside $\sim$1.5 $\sigma$ level and indicates that any further improvement to the orbital solution has no statistical significance, but only cause slight changes in stellar parameters, which would still stay inside corresponding standard error.

3.3 Light curve modeling

In case of KIC 7385478, light curve modeling would not be practical by considering 65 722 long cadence data points, therefore we first phase the whole data with respect to the orbital period and then calculate binned light curve with a phase step of 0.002 by using freely available fortran code $lcbin$666http://www.astro.keele.ac.uk/$\sim$jkt/codes.html$\#$lcbin written by John Southworth. We use the phase binned light curve for modeling the eclipsing binary. The first step of the light curve modeling is to find geometric and physical parameters of the system, based on phase binned data and the second step is to calculate theoretical light curve with the best geometric and physical parameters and subtract the theoretical model from the whole long cadence data in order to inspect possible out of eclipse variations. We plot the phase binned light curve in Figure 4, panel $a$.

By the quick inspection of the phase binned light curve, one may see that the system is on a circular orbit and there is a possibility of semi-detached configuration for the system. It may also be noticed that light level around 0.75 phase is slightly higher than the light level around 0.25 phase (see Figure 4, panel $b$), which means this effect is dominant through 4 years of $Kepler$ observations. Difference in the depth of the primary and secondary minima is another important feature of the light curve, which indicates large difference between the temperatures of the components. Shape of the light curve at out of eclipse phases indicates variations beside the geometric eclipse events, such as distortion in geometric shape of the components, reflection, spots, and even oscillations as we will focus on in the next section.

We use 2015 version of the Wilson-Devinney code (Wilson & Devinney 1971; Wilson & Van Hamme 2014) for light curve modeling. Thanks to our spectroscopic observations, we have already determined the two most critical parameters of the light curve modeling process, i.e. effective temperature of the primary component and the mass ratio of the system. We fixed these two parameters during the modeling. Since the effective temperature of the primary component is at critical location where the convective outer envelope is very thin or almost becomes radiative, we carried out solutions by setting gravity darkening ($g_{1}$) and albedo ($A_{1}$) values to 0.32 and 0.5, respectively (for convective envelopes) and setting both parameters to 1.0 (for radiative envelopes). In both cases, there is almost negligible difference between solutions with radiative envelope and convective envelope assumption, where convective envelope assumption leads to a slightly lower residuals, thus we continue with 0.32 and 0.5 values for $g_{1}$ and $A_{1}$, respectively. We adopt $g_{2}$ and $A_{2}$ values as 0.32 and 0.5 for the secondary component, which are typical for convective envelopes. Both light curve and spectroscopic orbit solution indicate circular orbit for the system, therefore we assume synchronous rotation for the components, which is proper for circular orbits, and fix rotation parameter ($F$) of each component to 1.0. Here, $F$ is defined as the ratio of the axial rotation rate to the orbital rate. Linear limb darkening coefficients ($x{{}_{1}},x{{}_{2}}$) of the components are adopted from van Hamme (1993).

We start analysis with detached configuration, therefore we consider inclination of the orbit ($i$), temperature of the secondary component ($T_{2}$), dimensionless omega potentials of the primary and the secondary component ($\Omega_{1}$, $\Omega_{2}$), luminosity of the primary component ($L_{1}$) as adjustable parameters. In addition, we observe a general phase shift of 0.002 in the phase binned light curve, which possibly arises from a shift in the adopted ephemeris due to the third body (Borkovits et al. 2016), hence we also adjust phase shift during analysis. In a few iterations, we observe that the $\Omega_{2}$ value jiggles around the inner critical potential value and secondary component has entirely filled its Roche-lobe in the corresponding Roche geometry. Then we switch to the semi-detached configuration and fixed $\Omega_{2}$. In this case, we consider distorted shape of the secondary component and adopt $g_{2}$, $A_{2}$ and $x{{}_{2}}$ as adjustable parameters during iterations. After a few iterations, we achieve statistically the best parameter set, however, we still observe that the residuals from the solution exhibit additional wave-like variation through an orbital cycle.

During the iterations, we noticed that the contribution of the secondary component to the total light is not more than 15%, therefore we expect negligible contribution from the secondary component to the wave-like variation in residuals. The most possible explanation could be that the Roche-lobe filled secondary star transfers its own mass to the primary component through the inner Lagrange point, L1, of the system. The transferred mass possibly hits directly to the photosphere of the primary component without forming an accretion disk, thus forms a local region warmer than the surrounding photosphere. Therefore, we consider this hot region as a bright spot on the primary component and add a single bright spot into our light curve model. At first step, we adopt spot longitude and radius as adjustable parameters together with eclipsing binary parameters and fixed the spot co-latitude and temperature factor to 45*∘* and 1.03, respectively. After these parameters are adjusted, we fix the spot longitude and radius, and adjust spot co-latitude and temperature factor. Until we achieve the best solution, we adjust different number of spot parameters simultaneously with a different combinations. When we reach the best solution, we adopt all parameters (spot and eclipsing binary) adjustable and run a single iteration in order to obtain statistical uncertainties. We tabulate our results in Table 3. Note that we do not give the internal error of the $T_{2}$ since it is unrealistically small ($\sim$ 1 K), therefore we adopt the uncertainty of $T_{1}$ estimated in Section 3.1 as the uncertainty of $T_{2}$. In Figure 4 we over plot the best fit models of spotless and spotted solutions (panel $a$ and $b$) and show residuals from both solutions in panel $c$ and $d$ in the same figure.

3.4 Physical properties and evolutionary status

Combining results from spectroscopic orbital solution and light curve modeling, we calculate absolute physical parameters of the system listed in Table 4. Inspecting the parameters, we immediately see that the secondary component has a very large radius compared to its mass which causes lower gravity compared to a typical main sequence star. A typical main sequence star has a radius of $\sim$0.41 R*⊙*(Gray 2005). Adopting temperature calibration of Gray (2005), we estimate the spectral type of the secondary as K4 III-IV. All these indicate that the lower mass secondary component has already evolved off the main sequence, which seems contradictory to the stellar evolution theory. However, we know that the secondary component has already filled its Roche-lobe and there must be mass transfer from secondary to the primary via inner Lagrange point, L1. The situation could be explained as the secondary component was actually the more massive component in the system and as it evolved off the main sequence, it filled its Roche-lobe and started a mass transfer to the less massive component (in our case, the primary component). We observe the effect of the mass transfer as hot spot on the primary component. Continuous mass transfer in time finally reversed the mass ratio of the system, therefore the final configuration became like a more massive main sequence star and a less massive sub-giant star. This is typical scenario adopted for Algols, which are well known to have reverted their mass ratio because of Roche lobe overflow.

In Figure 5, we show the position of the primary component of KIC 7385478 on Hertzsprung-Russel diagram, together with some eclipsing binaries with known $\gamma$ Doradus type pulsating components, i.e. VZ CVn (Ibanoǧlu et al. 2007), CoRot 102918586 (Maceroni et al. 2013), KIC 11285625 (Debosscher et al. 2013), KIC 9851142 (Çakırlı 2015), KIC 9851944 (Guo et al. 2016), V2653 Oph (Çakırlı & Ibanoglu 2016), CoRot 100866999 (Chapellier & Mathias 2013). Our target is located in the middle of the $\gamma$ Doradus instability strip given by Warner et al. (2003), which takes us to the possibility of $\gamma$ Doradus type pulsation on the primary component.

3.5 The out-of-eclipse variations

In order to inspect the out-of-eclipse variations, we first construct a theoretical light curve by using the parameters in Table 3, then subtract it from whole long cadence data and obtain residuals. In this process, we first divide the whole long cadence data into subsets where each subset covers a single orbital cycle. Then, we run differential correction program of the Wilson-Devinney code by only adjusting ephemeris reference time of the related subset and keeping all remaining parameters fixed. This process eliminates any shift in the ephemeris reference time due to the light time travel effect caused by the third body and gives correct residuals.

In Figure 6, we plot the whole residuals (panel $a$), and a sample of residuals covering a time span of a month (panel $b$ and $c$). One may easily notice that removing eclipsing binary model from the data unravels a clear variation with a dominant period of $\sim$ 0.5 day (panel $b$) with variable amplitude, which suggests the possibility of $\gamma$ Doradus type pulsation (Kaye et al. 1999) on the primary component, with a beat period of about 12 days. In addition, we observe occasional increase in amplitude of residuals (panel $c$), where the most dominant period becomes the orbital period, however by keeping 0.5 day variation as small humps and pits through an orbital cycle.

We apply multi-frequency analysis to the residual data using $pysca$ software package (Herzberg & Glogowski 2014) to investigate these variations. In case of continuous and long time series photometry, $pysca$ is very practical to automatically extract significant frequencies above a defined SNR limit. We start with Fourier analysis of the data for the frequency between 0 and 24.498 cycle/day (c/d), where the 24.498 denotes the nyquist frequency. In Figure 7, we show amplitude spectrum of the residuals (panel $a$). We observe that the dominant frequencies are located in lower frequency region ($f$ $<$ 5 c/d) and consider this region for frequency extraction process. In this process, the most dominant frequency in the amplitude spectrum is determined. Then its amplitude and phase are calculated via $A_{i}\Sigma sin[2\pi(f_{i}t+\phi_{i})]$, where t is the time of the corresponding measurement, while $A_{i}$, $f_{i}$ and $\phi_{i}$ show the amplitude, frequency and phase of the $i$th frequency, respectively. Next, this frequency is removed from the data and the same process is repeated for the remaining ”prewhitened” residuals. We adopt criteria of Breger et al. (1993), which puts a lower SNR limit of 4 for a frequency to be accepted as significant. Uncertainties of the extracted frequencies are estimated as $\sim$7$\times 10^{-3}$ d*-1*, which is determined by using the Rayleigh criterion.

This process leads to 735 frequencies above SNR limit of 4. The most dominant two peaks are located at 2.0252 c/d and 1.9427 c/d, corresponding 0.4938 day 0.5147 day and we define these frequencies as $\gamma$ Doradus type pulsation frequencies. Third and fourth peaks are 0.6034 c/d and 0.6023 c/d, corresponding 1.6573 day and 1.6603 day, and these frequencies indicate the orbital frequency. Beyond the first four frequencies, we do not find any independent frequency but low and high order combination of pulsation and orbital frequencies. We list the extracted frequencies in Table 5.

In Figure 7, we also show close view around orbital frequency (panel $b$) and pulsation frequencies (panel $c$ and $d$). One may easily notice that peaks at pulsation frequencies are single and sharp peaks, while the peak at the orbital frequency is more shallow in amplitude and broader. This picture leads us to calculate amplitude spectrum of long cadence data for each $Kepler$ data quarter separately and investigate behavior of the dominant peaks in 4 years time range. We plot amplitude spectrum of each quarter separately in Figure 8. Amplitudes of the pulsation frequencies are almost stable in all quarters, while the amplitude of the orbital frequency clearly varies in time, which causes shallow and broader peak(s) at the orbital frequency in the amplitude spectrum of the 4 year long cadence data.

4 Summary and discussion

Photometric and spectroscopic analysis of KIC 7385478 shows that the system is a low mass ratio ($q$ = 0.21) eclipsing binary formed by an F1V primary, and a K4III-IV secondary components which entirely fills its Roche lobe. This means the system can be classified as semi-detached. Physical properties of the components suggest that the low mass secondary component has already evolved off the main sequence before the dwarf primary component, which is not expected in the scope of the stellar evolution theory. Continuous mass transfer between the components could explain this situation, which could change the roles (and masses) of the components in long time scales. If the mass transfer does not lead to a disk, a hot spot forms on the photosphere of the mass gaining component, where the streaming matter directly hits to its photosphere. Almost continuous $Kepler$ long cadence photometry through four years indicates a brightness level difference between orbital quadrature phases (i.e. 0.25 and 0.75 phases), which could be explained by this kind of a hot spot on the primary (mass gaining) component. At that point, the assumption of an impact system (no disc) can be justified by the mass ratio (0.22) and fractional radius of the primary (0.21) according to Lubow-Shu criterium (Lubow & Shu 1975) for disc formation. Furthermore, we do not observe any emission feature in Hα, indicates the absence of the disk and this is consistent with impact system assumption.

The position of the primary component on HR diagram corresponds to the middle of the region where $\gamma$ Doradus variables are located and suggests intrinsic light variations related to pulsation. When we subtract the eclipsing binary model from the long cadence data, we see clear signal which has a dominant period of $\sim$0.5 day indicating $\gamma$ Doradus type pulsations and confirms the intrinsic variation suggested by the position on HR diagram. However, we also observe occasional increase in the amplitude of the residuals where the dominant period becomes the orbital period, including $\sim$0.5 day variation as a smaller amplitude variation. Multiple frequency analysis of the residuals results in 735 frequencies, where the first two frequencies indicate $\gamma$ Doradus type pulsations, while the third and the fourth frequencies correspond to the orbital frequency. Most of the remaining 731 frequencies correspond to either the orbital frequency or its harmonics. This indicates an additional light variation with a period almost identical to the orbital one and this variation is usually suppressed by the light variation originating from pulsations.

Figure 7 and Figure 8 clearly show that the amplitude of the orbital frequency varies in time and this causes broad and shallow peak structure at the orbital frequency in amplitude spectrum of 4 year data, which is another evidence for the additional light variation mentioned above. Possible star spot activity originating from the cool secondary component could easily cause this kind variation, thus leads to many low and high order combinations of the pulsation and orbital frequencies in the amplitude spectrum, especially in period analysis of continuous long term data, which covers $\sim$4 years in case of KIC 7385478. It is also known that stellar magnetic activity may cause orbital period modulations via mechanism proposed by Applegate (1992). Another possibility for orbital period modulations is variable mass transfer rate in the donor. Since the typical time scales for Applegate mechanism and variable mass transfer rate are decades or longer, these are not comparable to the modulations observed in residual data of KIC 7385478, therefore are not likely.

Estimated spectral type and luminosity class of the secondary component provide support for the spot activity possibility. According to the eclipsing binary model, 14% contribution of the secondary component to the total light is expected at very broad band $Kepler$ filter, which means we may observe small amplitude light variation due to the possible spot activity of the secondary star. However, when we check the observed spectra of the system, we do not observe any emission feature in Ca II H & K lines, which are very sensitive to the chromospheric activity in cool stars. Considering the temperatures and radii of the components, one may easily conclude that the contribution of the secondary component to the total light around 3950 Å  is almost completely negligible, hence, considering only the spectroscopic indicators, we may not arrive at the conclusion on the existence of star spot activity in case of KIC 7385478. Therefore, we can still speculate that the secondary component might have star spot activity, which seems the most possible cause of the occasional amplitude increase in the light residuals.

Acknowledgments

We thank to TUBITAK for a partial support in using RTT150 (Russian-Turkish 1.5-m telescope in Antalya) with project number 14BRTT150-667. This paper includes data collected by the Kepler mission. Funding for the Kepler mission is provided by the NASA Science Mission Directorate. Some of the data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST). STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. Support for MAST for non-HST data is provided by the NASA Office of Space Science via grant NNX13AC07G and by other grants and contracts.