MNRAS 000,1–15 (2020) Preprint 17 September 2021 Compiled using MNRAS LATEX style file v3.0 The intermediate polar cataclysmic variable GK Persei 120 years after the nova explosion: a first dynamical mass study A. Álvarez-Hernández1,2¢, M. A. P. Torres1,2, P. Rodríguez-Gil1,2, T. Shahbaz1,2, G. C. Anupama3, K. D. Gazeas4, M. Pavana3,5, A. Raj6, P. Hakala7, G. Stone8, S. Gomez9, P. G. Jonker10,11, J.-J. Ren12, G. Cannizzaro10,11, I. Pastor-Marazuela13,14, W. Goff15, J. M. Corral-Santana16, R. Sabo17 1Instituto de Astrofísica de Canarias, E-38205 La Laguna, Tenerife, Spain 2Departamento de Astrofísica, Universidad de La Laguna, E-38206 La Laguna, Tenerife, Spain 3Indian Institute of Astrophysics, 560 034 Bangalore, India 4Section of Astrophysics, Astronomy and Mechanics, Department of Physics, National and Kapodistrian University of Athens, GR-15784 Zografos, Athens, Greece 5Department of Physics, Pondicherry University, Puducherry 605014, India 6Dept. of Physics & Astrophysics, University Road, University Enclave, Delhi 110007, India 7Finnish Centre for Astronomy with ESO (FINCA), Quantum, University of Turku, FI-20014, Turku, Finland 8First Light Observatory Systems, 9 Wildflower Way, Santa Fe, NM 87506 USA 9Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA 10SRON, Netherlands Institute for Space Research, Sorbonnelaan 2, 3584 CA, Utrecht, The Netherlands 11Department of Astrophysics/ IMAPP, Radboud University, Heyendaalseweg 135,6525 AJ, Nijmegen, The Netherlands 12CAS Key Laboratory of Space Astronomy and Technology, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, People’s Republic of China 13Anton Pannekoek Institute, University of Amsterdam, Postbus 94249, 1090 GE Amsterdam, The Netherlands 14ASTRON, the Netherlands Institute for Radio Astronomy, Oude Hoogeveensedijk 4, 7991 PDDwingeloo, The Netherlands 15American Association of Variable Star Observers, 49 Bay State Rd. Cambridge, MA 02138 16European Southern Observatory, Alonso de Córdova 3107, Vitacura, Casilla 19001, Santiago de Chile, Chile 17American Association of Variable Star Observers, 1344 Post Dr. Bozeman, MT 59715 Accepted XXX. Received YYY; in original form ZZZ ABSTRACT We present a dynamical study of the intermediate polar and dwarf nova cataclysmic variable GK Per (Nova Persei 1901) based on a multi-site optical spectroscopy and '-band photometry campaign. The radial velocity curve of the evolved donor star has −1 a semi-amplitude 2 = 126.4 ± 0.9 km s and an orbital period % = 1.996872 ± 0.000009 d. We refine the projected rotational −1 velocity of the donor star to Erot sin 8 = 52 ± 2 km s which, together with 2, provides a donor star to white dwarf mass ratio @ = "2/"1 = 0.38 ± 0.03. We also determine the orbital inclination of the system by modelling the phase-folded ellipsoidal = ◦ ± ◦ = ¸0.16 = ¸0.07 light curve and obtain 8 67 5 . The resulting dynamical masses are "1 1.03−0.11 M and "2 0.39−0.06 M at 68 per cent confidence level. The white dwarf dynamical mass is compared with estimates obtained by modelling the decline light curve of the 1901 nova event and X-ray spectroscopy. The best matching mass estimates come from the nova light curve models and an X-ray data analysis that uses the ratio between the Alfvén radius in quiescence and during dwarf nova outburst. Key words: accretion, accretion discs – binaries: close – novae, cataclysmic variables – stars: individual: GK Per (Nova Persei arXiv:2107.06932v2 [astro-ph.SR] 16 Sep 2021 1901) 1 INTRODUCTION part of the accretion flow. Polars are the most extreme magnetic CVs: the strong WD magnetic field (퐵 & 107 G) prevents the formation Cataclysmic variables (CVs) are binary systems where a non- of an accretion disc and forces the transferred material to follow the degenerate star fills its Roche lobe and transfers matter towards an field lines onto one or both poles of the WD (Chanmugam & Wagner accreting white dwarf (WD; Kraft 1964, see Warner 1995 and refer- 1977; see Cropper 1990 for a review). In contrast, in intermediate ences therein). For a weakly magnetic WD, the mass from the donor polars (IPs) the magnetic field is only able to take control over the star is accreted onto the surface via an accretion disc. In magnetic transferred plasma in close proximity to the WD (Patterson 1994). In CVs, however, the magnetic field is strong enough to dominate at least these systems, the accretion disc is truncated at a certain radius from the WD and the disc accretion flow is funneled from there to its mag- netic poles along the field lines. A remarkable difference between ¢ E-mail: [email protected] © 2020 The Authors 2 A. Álvarez-Hernández et al. both types of magnetic CV is the degree of synchronization of the by Hachisu & Kato(2007) and Shara et al.(2018), respectively. ¸0.34 WD spin with the orbit: in IPs the spin period is usually significantly In addition, estimates of the WD mass ranging from 0 52−0.16 to shorter than the orbital period, while both periods are nearly equal 1.24 ± 0.10 M have been obtained through modelling of X-ray for most polars (see e.g. Norton et al. 2004). spectra (see Section4 for details). GK Per was discovered as a nova on 1901 February 22 by Scottish In this work, we present the first dynamical study of GK Per that amateur astronomer Thomas David Anderson (Williams 1901). It yields reliable masses for the WD and the donor star. The paper is peaked at a visual apparent magnitude of 0.2 mag. After years of structured as follows: in Section2 we describe the time-resolved opti- irregular fluctuations in brightness with amplitudes up to 1.5 mag cal photometry and spectroscopy observations and the data reduction. and several dozens of days duration, in 1948 it reached a quiescence From the analysis of the absorption lines of the donor star we ob- state (<+ ' 13 mag) and started to show 1 − 3-mag dwarf nova tain its radial velocity curve (Section 3.1), constrain its spectral type outbursts that typically last 50 d and recur about every three years (Section 3.2) and determine its rotational broadening (Section 3.3). (Hudec 1981; Bianchini et al. 1982; Sabbadin & Bianchini 1983; In Section 3.4 we present the modelling of the '-band light curve, Šimon 2002). which provides the orbital inclination for the first time. The stellar Crampton et al.(1986) reported a binary orbital period of nearly dynamical masses are then obtained using the measured quantities. 2 d and Watson et al.(1985) found a WD spin period of 351 s These are discussed in Section4, where we also compare the WD in the modulation of the hard X-ray emission, thus confirming the dynamical mass with the available estimates from X-ray spectral fit- IP nature of GK Per. The spin period also modulates the *-band ting and modelling of the nova light curve. Finally, we draw our flux (Patterson 1991) and the equivalent width of the emission lines conclusions in Section5. at optical wavelengths (Garlick et al. 1994; Reinsch 1994). 퐽-band circular polarimetry of GK Per in quiescence is consistent with a null detection (Stockman et al. 1992). However, the intensity of the 2 OBSERVATIONS AND DATA REDUCTION WD magnetic field is estimated at 퐵 ∼ 105 G from X-ray spectral modelling (Wada et al. 2018). The 1.9968-d orbital period of GK Per makes a given orbital phase The geometrical, kinematic and physical properties of the nova occur 4.6 min earlier every next orbital cycle. In addition, during shell in GK Per as well as its interaction with its surroundings a typical 10-h observing night only 20 per cent of the orbit can be have been studied in detail at different frequencies (Seaquist et al. covered. Thus, the fact that the orbital period is close to an integer 1989; Scott et al. 1994; Anupama & Kantharia 2005; Liimets et al. number of days precludes ground-based observers at a single location 2012, and references therein). In particular, far-infrared observations from achieving entire photometric coverage of the orbit in contem- showed that the nova shell is embedded in an ancient, possibly bipolar poraneous nights. This makes the light curve strongly susceptible to planetary nebula centred on the binary and extending ≈ 17 arcmin to aperiodic night-to-night accretion variability. To overcome this dif- the NW and SE (Dougherty et al. 1996). At the time of its discovery, ficulty and thus obtain the full ellipsoidal modulation produced by this nebula was interpreted as being the remnant of the binary com- the donor star, we performed multi-site photometry between 2017 mon envelope phase (Bode et al. 1987). However, ejecta presumably December and 2018 January. We also took multi-site spectroscopy from the WD progenitor star expelled during a second asymptotic during 2017–2019 to improve the light curve modelling and obtain giant branch phase (and thus a second common envelope event), a full dynamical determination of the system parameters. triggered by a period of high mass transfer rate from the donor star In this Section, we describe all the collected data sets. Tables1 −7 −1 19 −1 (¡ 3 × 10 M yr = 1.9 × 10 g s ), was proposed as a more and2 summarize the spectroscopic and photometric observations, likely origin for the nebula (Dougherty et al. 1996). respectively. Note that we have adopted orbital phase 0 as the moment Several spectroscopic classifications of the donor star in GK Per of inferior conjunction of the donor star. have been reported by different authors: K2 V-IVp (Kraft 1964; Gallagher & Oinas 1974), K0 III-IV (Crampton et al. 1986), K2–3 V (Reinsch 1994) and K1 IV (Morales-Rueda et al.