Phase-Resolved Cyclotron Spectroscopy of Polars

Ryan Campbell

In this thesis we use phase-resolved cyclotron spectroscopy to study polars. Polars are a subset of cataclysmic variables where the primary WD is highly magnetic. In this case, the accretion flow is constrained along the magnetic field lines and eventually deposited on the WD, where the accreting material interacts with the atmosphere, forming a standing hydrodynamic shock at a location termed the accretion region, and emitting cyclotron radiation. Due to its nature, cyclotron radiation from polars falls at either optical or NIR wavelengths. While a substantial amount of optical cyclotron spectra have been published on polars, the NIR remains relatively unstudied. In this thesis, we present NIR spectroscopy for fifteen polars. Additionally, while a single cyclotron spectrum is needed to constrain cyclotron spectra, phase-resolved spectroscopy allows for a more in-depth analysis of the shock structure and the geometry of the accretion region. Of the fifteen polars observed, eight yielded spectra of adequate quality to be modeled in this manner: EF Eri, EQ Cet, AN UMa, VV Pup, AM Her, ST LMi, MR Ser, and MQ Dra.

Initially, we used the industry standard "Constant Lambda (CL)" code to model each object. The code is fast, but produces only globally averaged values of the salient shock parameters: B - the magnetic field strength, kT - the plasma temperature, logLambda - the "size parameter" of the accretion column, and theta- the viewing angle between the observer and the magnetic field. For each object we present CL models for our NIR phase-resolved cyclotron spectra. Subsequently, we use a more advanced "Structured-Shock" code built Fischer & Beuermann (2001) ("F&B") to remodel three objects: EQ Cet, MQ Dra, and EF Eri. The F&B code allows for input of more physical parameters and crucially does ray tracing through a simulated one-dimensional accretion column. To determine the outgoing spectrum, temperature and velocity profiles are needed to reconstruct the characteristics of the plasma at each location. However, the program to compute these profiles has been lost. Thus, we reconstructed these profiles. Finally, we compare the results of the CL and F&B codes to determine when the extra complexity and significantly longer computational times of F&B modeling are necessary to understand these systems.