We present high time resolution optical photometric data of the polar VV Puppis obtained simultaneously in three filters (u', HeII $\lambda$4686, r') with the ULTRACAM camera mounted at the ESO-VLT telescope. An analysis of a long 50 ks XMM-Newton observation of the source, retrieved from the database, is also provided. Quasi-periodic oscillations (QPOs) are clearly detected in the optical during the source bright phase intervals when the accreting pole is visible, confirming the association of the QPOs with the basis of the accretion column. QPOs are detected in the three filters at a mean frequency of $\sim$ 0.7 Hz with a similar amplitude $\sim$ 1\%. Mean orbitally-averaged power spectra during the bright phase show a rather broad excess with a quality factor Q= $\nu$/$\Delta \nu$ = 5-7 but smaller data segments commonly show a much higher coherency with Q up to 30. The XMM (0.5--10 keV) observation provides the first accurate estimation of the hard X-ray component with a high kT $\sim$ 40 keV temperature and confirms the high EUV-soft/hard ratio in the range of (4--15) for VV Pup. The detailed X-ray orbital light curve displays a short $\Delta \phi \simeq 0.05$ ingress into self-eclipse of the active pole, indicative of a accretion shock height of $\sim$ 75 km. No significant X-ray QPOs are detected with an amplitude upper limit of $\sim$30\% in the range (0.1--5) Hz. Detailed hydrodynamical numerical simulations of the post-shock accretion region with parameters consistent with VV Pup demonstrate that the expected frequencies from radiative instability are identical for X-rays and optical regime at values $\nu$ $\sim$ (40--70) Hz, more than one order magnitude higher than observed. This confirms previous statements suggesting that present instability models are unable to explain the full QPO characteristics within the parameters commonly known for polars.

A new target design is presented to model high-energy radiative accretion shocks in polars. In this paper, we present the experimental results obtained on the GEKKO XII laser facility for the POLAR project. The experimental results are compared with 2D FCI2 simulations to characterize the dynamics and the structure of plasma flow before and after the collision. The good agreement between simulations and experimental data confirm the formation of a reverse shock where cooling losses start modifying the post-shock region. With the multi-material structure of the target, a hydrodynamic collimation is exhibited and a radiative structure coupled with the reverse shock is highlighted in both experimental data and simulations. The flexibility on the laser energy produced on GEKKO XII, allowed us to produce high-velocity flows and study new and interesting radiation hydrodynamic regimes between those obtained on the LULI2000 and Orion laser facilities.

Some polars show quasi-periodic oscillations (QPO) in their optical light curves which have been interpreted as the result of shock oscillations driven by the cooling instability. Although numerical simulations can recover this physics, they wrongly predict QPOs in the X-ray luminosity and have also failed to reproduce the observed frequencies, at least for the limited range of parameters explored so far. Given the uncertainties on the observed polar parameters, it is still unclear whether simulations can reproduce the observations. The aim of this work is to study QPOs covering all relevant polars showing QPOs. We perform numerical simulations including gravity, cyclotron and bremsstrahlung radiative losses, for a wide range of polar parameters, and compare our results with the astronomical data using synthetic X-ray and optical luminosities.We show that shock oscillations are the result of complex shock dynamics triggered by the interplay of two radiative instabilities. The secondary shock forms at the acoustic horizon in the post-shock region in agreement with our estimates from steady-state solutions. We also demonstrate that the secondary shock is essential to sustain the accretion shock oscillations at the average height predicted by our steady-state accretion model. Finally, in spite of the large explored parameter space, matching the observed QPO parameters requires a combination of parameters inconsistent with the observed ones. This difficulty highlights the limits of one-dimensional simulations, suggesting that multi-dimensional effects are needed to understand the non-linear dynamics of accretion columns in polars and the origins of QPOs.

Quasi-periodic oscillations (QPOs) are observed in the optical flux of some polars with typical periods of 1 to 3 s but none have been observed yet in X-rays where a significant part of the accreting energy is released. QPOs are expected and predicted from shock oscillations. Most of the polars have been observed by the XMM-Newton satellite. We made use of the homogeneous set of observations of the polars by XMM-Newton to search for the presence of QPOs in the (0.5-10 keV) energy range and to set significant upper limits for the brightest X-ray polars. We extracted high time-resolution X-ray light curves by taking advantage of the 0.07 sec resolution of the EPIC-PN camera. Among the 65 polars observed with XMM-Newton from 1998 to 2012, a sample of 24 sources was selected on the basis of their counting rate in the PN instrument to secure significant limits. We searched for QPOs using Fast Fourier Transform (FFT) methods and defined limits of detection using statistical tools. Among the sample surveyed, none shows QPOs at a significant level. Upper limits to the fractional flux in QPOs range from 7% to 71%. These negative results are compared to the detailed theoretical predictions of numerical simulations based on a 2D hydrodynamical code presented in Paper II. Cooling instabilities in the accretion column are expected to produce shock quasi-oscillations with a maximum amplitude reaching ~ 40% in the bremsstrahlung (0.5-10 keV) X-ray emission and ~ 20% in the optical cyclotron emission. The absence of X-ray QPOs imposes an upper limit of ~ (5-10) g.cm-2.s-1 on the specific accretion rate but this condition is found inconsistent with the value required to account for the amplitudes and frequencies of the observed optical QPOs. This contradiction outlines probable shortcomings with the shock instability model.

Magnetic cataclysmic variables are close binary systems containing a strongly magnetized white dwarf that accretes matter coming from an M-dwarf companion. High-energy radiation coming from those objects is emitted from the accretion column close to the white dwarf photosphere at the impact region. Its properties depend on the characteristics of the white dwarf and an accurate accretion column model allows the properties of the binary system to be inferred, such as the white dwarf mass, its magnetic field, and the accretion rate. We study the temporal and spectral behaviour of the accretion region and use the tools we developed to accurately connect the simulation results to the X-ray and optical astronomical observations. The radiation hydrodynamics code Hades was adapted to simulate this specific accretion phenomena. Classical approaches were used to model the radiative losses of the two main radiative processes: bremsstrahlung and cyclotron. The oscillation frequencies and amplitudes in the X-ray and optical domains are studied to compare those numerical results to observational ones. Different dimensional formulae were developed to complete the numerical evaluations. The complete characterization of the emitting region is described for the two main radiative regimes: when only the bremsstrahlung losses and when both cyclotron and bremsstrahlung losses are considered. The effect of the non-linear cooling in- stability regime on the accretion column behaviour is analysed. Variation in luminosity on short timescales (~ 1 s quasi-periodic oscillations) is an expected consequence of this specific dynamic. The importance of secondary shock instability on the quasi-periodic oscillation phenomenon is discussed. The stabilization effect of the cyclotron process is confirmed by our numerical simulations, as well as the power distribution in the various modes of oscillation.