XMM-Newton observation of MV Lyr and the sandwiched model confirmation
A. Dobrotka J. -U. Ness S. Mineshige A. A. Nucita




firstpage–lastpage 2017


Maximal violation of n-locality inequalities in a star-shaped quantum network

    Fabio Sciarrino
    Accepted ???. Received ???; in original form December 30, 2023
===============================================================================






We present spectral and timing analyses of simultaneous X-ray and UV observations of the VY Scl system MV Lyr taken by , containing the longest continuous X-ray+UV light curve and highest signal-to-noise X-ray (EPIC) spectrum to date. The RGS spectrum displays emission lines plus continuum, confirming model approaches to be based on thermal plasma models. We test the sandwiched model based on fast variability that predicts a geometrically thick corona that surrounds an inner geometrically thin disc. The EPIC spectra are consistent with either a cooling flow model or a 2-T collisional plasma plus Fe emission lines in which the hotter component may be partially absorbed which would then originate in a central corona or a partially obscured boundary layer, respectively. The cooling flow model yields a lower mass accretion rate than expected during the bright state, suggesting an evaporated plasma with a low density, thus consistent with a corona. Timing analysis confirms the presence of a dominant break frequency around log(f/Hz) = -3 in the X-ray Power Density Spectrum (PDS) as in the optical PDS. The complex soft/hard X-ray light curve behaviour is consistent with a region close to the white dwarf where the hot component is generated. The soft component can be connected to an extended region. We find another break frequency around log(f/Hz) = -3.4 that is also detected by . We compared flares at different wavelengths and found that the peaks are simultaneous but the rise to maximum is delayed in X-rays with respect to UV.



accretion, accretion discs - turbulence - stars: individual: MV Lyr - novae, cataclysmic variables




§ INTRODUCTION



A great variety of objects such as cataclysmic variables (CVs), symbiotic systems, X-ray binaries or active galactic nuclei are powered by a common physical process: accretion. In binaries the process is based on mass loss from a companion star. The transported gas falls towards the central compact object, and in the absence of a strong magnetic field, an accretion disc forms. The central accretor can either be a white dwarf in the case of CVs and symbiotic systems or a neutron star or a black hole in the case of X-ray binaries (see e.g.  or  for a review).

The family of CVs is divided into several subclasses based on characteristic variability patterns. The most common are the dwarf novae showing quasiregular outbursts with durations of several days and appearing on a time scale of 10 - 100 days (see  for review). VY Scl systems spend most of their life time in a high state while they sporadically transition into a transient low state. This alternating behaviour can be explained by corresponding changes in accretion rate. The high state in VY Scl systems is stable for a relatively long time, while in a dwarf nova it is only temporary (outbursts). The shorter durations of high states in dwarf novae can be explained by the mass accretion rate being unstable in the framework of the disc instability model (see  for review). Meanwhile, in VY Scl systems, the mass accretion rate remains above the critical limit required for stability, explaining the longer duration of high states. The low state events are generated by a sudden drop of mass transfer from the secondary or even a total stop of mass transfer (, ).

Observations of accreting systems of various kinds suggest that the basic characteristic of accretion is fast stochastic variability (a.k.a. flickering), see, e.g. <cit.>, <cit.>, <cit.>. The remarkable observational similarities suggest that the same physical mechanism is responsible for flickering in all accretion systems (e.g. , , , ). Flickering has three basic observational characteristics; 1) linear correlation between variability amplitude and log-normally distributed flux (so called rms-flux relation) observed in a variety of accreting systems such as X-ray binaries or active galactic nuclei (), CVs (, ) and symbiotic systems (), 2) time lags where flares reach their maxima slightly earlier in the blue than in the red () and 3) red noise in power density spectra (PDS). The shape of such PDS can be a simple power law (see e.g. , ), a broken power law (two power law components) with a single break frequency in between (see e.g. , ) or a multicomponent PDS with several characteristic break frequencies (see e.g. , , , ).

For sufficiently detailed studies of the PDS, a high cadence and long, continuous light curves are needed. An ideal opportunity is offered by the  mission which allowed us to discover a multicomponent PDS in the optical waveband of the two CVs V1504 Cyg () and MV Lyr (). The former case shows two or three characteristic break frequencies, while for the latter, four components were found. It is believed that the accretion process from the companion star toward the white dwarf surface is structured and every single structure has its own flow characteristics behaving differently. Therefore, every single characteristic break frequency in a PDS can be a footprint of a separate accretion structure.

Different simulation techniques/models intended for identifying the sources of flickering (related to accretion structure) have been developed so far. A cellular-automaton model[A cellular automaton model consists of a regular grid of cells. For each cell, a set of adjacent cells is defined relative to the specified cell. An initial configuration is defined by assigning a state for each cell. Step by step, a new configuration is calculated, according to some rule that determines the new state of each cell in terms of the current state of the cell and the states of adjacent cells.] was proposed by <cit.>, where light fluctuations are produced by occasional flare-like events and a subsequent avalanche flow in the accretion disc atmosphere (see  for the original idea). Another cellular-automaton model was developed by <cit.>. It is based on changing the collection of magnetic flux tubes anchored in the disc, transporting angular momentum and driving accretion inhomogeneities. <cit.> developed a statistical model to simulate flickering based on the simple idea of angular momentum transport between two adjacent concentric rings in the accretion disc via discrete turbulent bodies. The method is based on a geometrically thin disc with a ratio of H/r < 0.01 (H is the disc scale height and r is the distance from the centre). Another way of the PDS reproduction was proposed by <cit.>. They derived an analytical expression for the fluctuating accretion rate in the disc, where fluctuations are generated at each radius on a local viscous time scale, and overall variability of the accretion rate in the innermost disc region is the product of all the fluctuations produced at all radii. Practically all mentioned models approve the basic idea of variations in the accretion rate that are produced at different disc radii (, , ), i.e. the mass accretion variability at the outer radii is propagating inside and is influencing the variability characteristics of the inner regions. Such process explains the mentioned linear rms-flux relation and log-normal flux distribution.

<cit.> applied the <cit.> method to fit the highest PDS break frequency log(f) = -3.01 ± 0.06 Hz detected in the optical PDS of the nova-like system MV Lyr (). The discs in CVs are believed to be geometrically thin and optically thick. During quiescence the mass accretion rate is low, resulting in a truncated disc (see  for review), i.e. a hole around the white dwarf is formed because of inefficient cooling[Cooling in an optically thin hot plasma in this case is based mainly on free-free transitions, therefore dependent on the square of the particle number density. The plasma energy is transported toward the disc by electron conduction, and a low concentration of electrons reduces cooling. The underlying disc matter thus heats and evaporates, which increases the local particle concentration until an equilibrium is reached.] generating evaporation of the matter (). Instead of the geometrically thin disc, a central geometrically thick (H/r no longer small compared to 1) optically thin corona forms due to evaporation. <cit.> determine a ratio H/r > 0.1 with a radius of 0.81^+0.20_-0.14× 10^10 cm and an α parameter of 0.6 - 1.0 (). The author interpreted such a disc as an expanded optically thin hot corona surrounding a geometrically thin standard disc, i.e. the so-called sandwiched model.

<cit.> applied a modelling method () to the study of MV Lyr. They find the geometrically thin disc to be responsible for the second lowest break frequency, while the lowest one can be generated by somehow enhanced activity of the outer disc rim. From their modelling, they also find a solution where part of the central disc can be responsible for the highest break frequency fitted by <cit.>. Although the derived disc radius of 10^10 cm and α = 0.8 - 1.0 are in agreement with <cit.>, the ratio H/r is completely different, i.e. H/r < 0.1 in <cit.> and H/r > 0.1 in <cit.>.

The explanation of the optical multicomponent PDS of the nova-like MV Lyr promises the most complex image of the accretion process in a CV so far. However, the character of the highest break frequency is still not perfectly clear because of two contradictory geometry conditions from numerical modelling, but the corona interpretation is very plausible. There is only one way to resolve this puzzle, i.e. a direct X-ray observation, because the hot optically thin corona is radiating in X-rays. If the corona interpretation is correct, the highest optical break frequency log(f) = -3.01 ± 0.06 Hz must be detected in the X-ray PDS as well.

While MV Lyr has been observed in X-rays before, none of them were sufficient for the required timing analysis. During the  all-sky survey, MV Lyr was observed three times for 722 s (in Oct. 1990), 20250 s (in Nov. 1992), and 2218 s (in May 1996) with the position-sensitive proportional counter (PSPC) being the first two observations coincident with high optical state periods and the last occurring during a low state. The PSPC 0.1-2.4 keV count rates were 0.079 ± 0.011, 0.069 ± 0.002 and < 0.0008 () corresponding (after re-scaling the assumed source distance from 320 pc to 500 pc) to luminosities of 4.87 × 10^31 erg s^-1, 3.87 × 10^31 erg s^-1, and < 1.20 × 10^30 erg s^-1, respectively. In contrast to our  observation, the PSPC bandpass does not cover the hard range above 2.4keV that contains a lot of critical information.

MV Lyr was also observed by  during both high (Obs ID 91443 lasting for 14240 s) and low state (Obs ID 32042, 3282 s). <cit.> concentrated on the former and found that the spectrum can be well fitted by a multi-temperature plasma emission model (CEVMKL in xspec). In order to test the hypothesis of the presence of any scattering effect of X-rays from a wind or an extended component, these authors realized that adding a power-law component increased the quality of the fit. In particular, they derived only a lower limit to the plasma temperature (kT > 21,keV), Γ≃ 0.82, α≃ 1.6 and hydrogen column density of 0.13 × 10^22 cm^-2. The 0.2-10 keV flux of 5.4 × 10^-12 erg s^-1 cm^-2 corresponds to a luminosity of 1.7 × 10^32 erg s^-1.

Using  data, <cit.> derived a 2 σ upper limit to the soft black-body component (if any) associated to the boundary layer finding kT < 6.6 eV. Taking the accretion rates derived by optical and UV bands, the standard disc models predict an optically thick boundary layer with temperature (13 - 33 eV) larger than that derived from the  data and luminosity (> 10^34 erg s^-1) much higher than that inferred by using the  data. Moreover, the X-ray luminosity in the 0.1 - 50 keV band of the source is ≃ 3.2 × 10^32 erg s^-1, so that its ratio to the disc luminosity (calculated from the UV and optical observations) is in the range 0.01 - 0.001 thus suggesting that the MV Lyr system in high state is characterized by an optically thin boundary layer associated with a low efficiency accretion model (an ADAF-like flow) and/or an X-ray corona on the inner disc and close to the WD.

<cit.> studied all  XRT spectra, thus the observation taken during low state and the same data studied by <cit.>. They found that a good fit is obtained (see their Table 5) with a two-component thermal plasma model but also with a thermal plasma to which a power-law is added. In the high state, MV Lyr is characterized by a 0.3 - 10 keV luminosity of 1.7 × 10^32 erg s^-1 (consistently with ) which reduces to 3.6 × 10^31 erg s^-1 in the low state.

In this paper we present timing and spectral analyses of our recent  observation of MV Lyr. The data are described in Section <ref>, the EPIC spectrum is analysed in Section <ref>, and timing analysis in Section <ref>. We discuss and summarize our results in Sections <ref> and <ref>, respectively.



§ OBSERVATIONS



Since the existing  and  X-ray observations either only cover a limited energy range or are too short for the sensitivity of timing analyses needed to put the  data into context, we have requested a 50-ks  observation. On 2015 September 6, this observation was realised during the more common high state. A 64.6-ks observation was taken under ObsID 0761040101 yielding 61.6, 63.3 and 63.5 ks for PN, MOS and RGS instruments, respectively. The Optical Monitor (OM) took 13 exposures of 2883s duration each in Image+Fast mode and the UVW1 filter inserted. The data were downloaded from the XMM-Newton Science Archive (XSA), and we obtained science products with the Science Analysis Software (SAS), version 14.0. We used the tool xmmextractor to re-generate calibrated events files from which in turn light curves were extracted for all instruments. Optimised extraction regions for the EPIC detectors were calculated with the tool eregionanalyse while for the Reflection Grating Spectrometer (RGS) and the optical monitor (OM), standard extraction regions were used by xmmextractor. We use the MOS light curves as a comparison while for the timing analysis we use only PN and OM light curves. The PN and OM light curves are shown in Fig. <ref> and the Julian date of the  observation in the long-term context of the AAVSO light curve in Fig. <ref>. For the spectral analysis we use both MOS and PN spectra while RGS spectra are only used for consistency checks.









§ SPECTRAL ANALYSIS



We have extracted X-ray spectra from the RGS and the EPIC cameras using the Science Analysis System (SAS) version 14.0. An RGS spectrum in flux units was obtained using rgsproc resulting in a merged spectrum from RGS1 and RGS2, making best use of all redundancies, especially filling chip gaps and bad pixels. The combined RGS spectrum is shown in Fig. <ref>. Prominent line transitions are marked with vertical lines and labels of the respective ions on top at the respective wavelengths which can be used as a check to see whether any of them are present. Clearly, the RGS spectrum is not sufficiently well exposed for a quantitative analysis, but it allows the conclusion that the spectrum consists of weak continuum below ∼ 17 Å (above 0.7 keV) with superimposed emission lines from oxygen, iron, and neon. Other lines, e.g., N vi may be present but can not be detected above the noise. The lines are too weak to derive any meaningful line ratios, and, e.g., density diagnostics with He-like triplets can not be applied. But we use the RGS as a consistency check as the high spectral resolution can potentially rule out spectral models, e.g., if they do not reproduce the emission lines that only the RGS can resolve.




For quantitative analysis, we therefore need to rely on low-resolution spectra from the three EPIC detectors PN and MOS1/2 over the energy range 0.1-10 keV. We merged MOS1 and MOS2 spectra with the SAS tool epicspeccombine while we use the PN spectrum separately in a simultaneous fit.






We use xspec to test spectral models against the combined MOS1+2 and PN spectra via simultaneous χ^2 minimization fitting. We take photoelectric absorption within the neutral interstellar medium plus potentially circumstellar material into account by the tbabs model (Tübingen absorption ). The only parameter of tbabs is the neutral hydrogen column density N_ H. As an independent estimate for the amount of interstellar absorption in the direction of MV Lyr we use the HEASARC NH tool[http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl]. Two values are computed from Galactic H i measurements in the Leiden/Argentine/Bonn (LAB) Survey and by <cit.> resulting in values of 5.35× 10^20 cm^-2 and 6.17× 10^20 cm^-2, respectively. During the fit to the X-ray spectra, we keep the parameter N_ H variable, starting with the interstellar value of 6 × 10^20 cm^-2, in the middle between the above values. If there are any measurable additional sources of circumstellar absorption, a higher value will result, while a lower value may be an indicator for using the wrong model.

In pursuit of finding evidence for either coronal or boundary layer emission, we test three spectral models. All three models are different realizations of the APEC model (), which basically assumes a collisional plasma in equilibrium, thus collisional ionizations and excitations being balanced by radiative recombination and de-excitations. This way a spectrum consisting of bremsstrahlung continuum (from free electrons) plus emission lines (from bound-bound transitions) is produced. The key parameter is the electron temperature T which is converted to a Maxwellian velocity distribution of electrons and ions. A low-density plasma is assumed with an arbitrary density of log(n_e)=1 (in cgs units), and that no photons are re-absorbed (thus optically thin). The intensity depends on the differential emission measure, basically the product of volume and density, V× n_e. For details we refer to <cit.>. An isothermal plasma is unlikely in nature, more likely is a broad distribution of temperatures. This can be approximated by a multi-temperature model or by integration over a continuous temperature distribution. A comparison of these two approaches is presented by <cit.>.
The APEC model assumes solar abundances by <cit.> that can be scaled by a factor to determine an overall metallicity. Meanwhile the VAPEC model allows modification of individual abundances. We used the VAPEC model iterating the Fe abundance if needed.


We first tested a single temperature model (1-T VAPEC) but were not able to reproduce the observation, thus the plasma is measurably not isothermal. We then defined a second temperature component (2-T VAPEC) with variable Fe abundance for each component and show the result in the top panel of Fig. <ref>, together with the combined MOS1+2 and PN spectra. The model from each temperature component is shown with dotted lines (for clarity only for MOS1+2). This model already reproduces the data fairly well with a value of χ^2_ red= 1.2 (at 3076 degrees of freedom). The Fe xxv lines are 6.7 keV can only be reproduced if the Fe abundance of the hot component is increased. Meanwhile, the abundance of the cooler component is solar, although with large errors. This means that the hot and cool emission could come from different regions with different metallicity. We added a Gaussian component at 6.4 keV with zero line width, thus only broadened by the instrument to determine whether there is any significant flux produced by Fe i. While some flux results, there is not really a line-like features. The inset of Fig. <ref> shows the PN spectrum around the Fe lines with the 2-T VAPEC model (blue line) and the additional emission in the Gaussian (light blue-shadings). The vertical dotted lines mark the energies of the four transitions listed in the left. Two nearby Fe xxv lines  appear to be resolved with a red-shifted and a blue-shifted component in the PN spectrum, but the model demonstrates that the energy resolution is not sufficient for such a conclusion.

The model parameters are given in the top part of Table <ref>. The values of emission measure, log(VEM) indicate that the hotter component dominates with a factor 10 over the cooler component. Both temperatures are typical of coronal plasma, and the 2-T VAPEC model can thus be considered representative of a coronal plasma with a large range of temperatures.

Closer inspection of the high-energy tail of the 2-T VAPEC model in Fig.  <ref> (and the inset) suggests that the model systematically underpredicts the continuum in the energy range ∼ 3-7.3 keV. This may be an indicator for the high-temperature component to be highly absorbed. This could occur if the high temperature component originates from the boundary layer, thus deeply embedded at the bottom of the accretion disc. To test this scenario, we modify the 2-T VAPEC model by introducing an additional absorber that acts only on the hot component while the main absorber acts on the entire emission. Since absorption within the accretion disc may not be uniform, we chose a partial absorber. The best-fit parameters are listed in the middle part of Table <ref> while the spectral model is compared with the data in the bottom part of Fig. <ref>. The lower temperature component is identical in the two models while the high-temperature component yields a lower value in the partial absorber model. The Fe abundance does not need to be modified from solar for this model, indicating that the abundance effect in the 2-T VAPEC model may not be real.
The overall value of N_ H is closer to the interstellar value in the partial absorber model than in the 2-T model where it is slightly higher. If ignoring the uncertainties in the interstellar value of N_ H, one might argue that the partial absorber model is more realistic, but the difference is marginal.

Finally, we test the cooling flow model that <cit.> have applied to CV spectra after noticing that their grating spectra seem consistent with a multitemperature thermal plasma with a relatively flat emission measure distribution. The flatness of this distribution can indicate an isobaric cooling flow, which assumes that the gas releases all of its energy in the form of optically thin radiation as it cools in a steady-state flow. The optically thin radiation can be modelled with the APEC model (see above), and the  xspec  model mkcflow assumes the flow as an interpolation between a minimum and a maximum APEC temperature. The mkcflow also has a flavour in which abundances can be modified individually but only for the full flow. The normalization parameter directly gives the total mass flow rate.

The mkcflow model has originally been developed for Galaxy Clusters, for which a parameter is needed to red-shift the model. In addition to red-shifting, the parameter z is used to get the flux from the model via distance/z. Owing to this dual purpose of z, the computation of the mass flow rate from the normalization fails for z=0, and for our Galactic source without the need for a red shift, we fix z at a small but non-zero value, z=5× 10^-8. This is equivalent to a Doppler shift of 0.015 km s^-1, which is negligible.

The best-fit mkcflow model is illustrated in Fig. <ref>, and the parameters are listed in the bottom part of Table <ref>. The inset of Fig. <ref> shows that the Fe lines are underestimated, even when leaving the Fe abundance free to vary. For the 2-T VAPEC model, a better reproduction of the Fe lines could be achieved because we could modify the abundances of the two components separately. The Fe abundance of the low-temperature component will be driven by the Fe L-shell lines around 1 keV while the Fe abundance of the high-temperature component will be driven by the K-shell lines at 6.7 keV. Since the Fe L-shell lines are much stronger, we believe the Fe abundance in the mkclow model is driven by the L-shell lines, leading to an overall lower Fe abundance and thus underprediction of the K-shell lines.
The mkcflow model also is a poor representation of the energy range 0.7-1 keV, however, the overall value of χ^2 is sufficiently small so this model is not necessarily unacceptable compared to the other two models.




Although the RGS spectra are not well-enough exposed for quantitative analyses, we can compare the results from the EPIC cameras with the RGS spectrum, and in Fig. <ref> we show the RGS1 and RGS2 count spectra in comparison with the 2-T VAPEC and mkcflow spectral models. We have not re-fitted nor re-normalized the models, thus only converted the models to equivalent count spectra by folding them through the spectral responses of the RGS. The main plot shows good agreement of the continuum (a result of the quality of the cross calibration) while the inset shows how much the models agree with selected emission lines. The strongest line, O viii (18.97 Å, bottom right inset) is well reproduced by both models, slightly better with the mkcflow model. The O vii triplet (21.6/21.8/22.1 Å, top right inset) is not at all reproduced, suggesting that both models may be missing some cooler plasma. The O vii forbidden line (22.1 Å) seems not detected while the intercombination line (21.8 Å) might be present. If this is real, the cooler plasma would have a density above ∼ 10^10 cm^-3, higher than typical coronal plasma and could thus come from the boundary layer.

In conclusion, the best fit to the spectra has been obtained with a cool optically thin plasma model plus a partially absorbed hotter component. This model can be interpreted as a coronal component with a temperature ∼ 1 keV (1.2× 10^7 K) and a hotter emission component of ∼ 6 keV (7× 10^7 K) from the boundary layer. In this picture the boundary layer is hidden behind semi-transparent (∼ 30% transparency) plasma with a hydrogen column density of N_ H∼ 4× 10^22 cm^-2 while absorption of the coronal plasma is consistent with the interstellar medium. In terms of emission measure, the boundary layer contributes 90%. However, the other models also lead to acceptable fits, also allowing the possibility that we are either seeing a multi-temperature coronal plasma or a cooling accretion flow, although the cooling flow model requires somewhat higher N_ H than interstellar.





§ TIMING ANALYSIS





 §.§ Power density spectra


The goal of the timing analysis is to search for characteristic break frequencies L_ i that have been detected in the optical PDS determined the from long-term light curve observed with  ().




For the timing analysis and PDS calculation we applied the Lomb-Scargle algorithm (). This method is particularly adequate for the OM data because it can handle gaps in the light curve. The high and low frequency limits within which the PDS can be studied are determined by the light curve characteristics. The high frequency end is usually limited by the white noise or PDS power rising to the Nyquist frequency. The lower frequency limit is usually determined by the duration of the light curve. The periodograms derived from the individual instruments are depicted in log-log scale in Fig. <ref>. The white noise starts to be noticeable around a frequency of log(f) = -2.4 Hz for the X-ray data, but for higher frequencies in the case of UV data. Therefore, we adopted the frequency of log(f) = -2.4 Hz as an upper limit for subsequent analysis.




Fig. <ref> reveals an apparent break frequency around L_1 in all instruments. Furthermore, the presence of the orbital frequency can clearly be seen in all PDSs except in the UV. While there is no obvious indication for other L_ i signals in X-rays, the OM data show some increased power around L_3 and L_4.

For a more detailed study we calculated the PDSs as we have done for the dwarf nova RU Peg (). The light curve is divided into several subsamples, and the corresponding periodograms in log-log scale of every subsample are averaged in order to get the PDS (). The main motivation of this method to average out random features in the individual periodograms, while only real intrinsic PDS features remain. Important is the basic rule, the more subsamples are used, the lower is the scatter in the PDS, while on the other hand, the studied frequency interval becomes narrower with shorter subsamples as the low frequency limit then increases. We subsequently binned the PDS into equally spaced bins in order to reduce the scatter even more and the mean values with the errors were fitted with a model. As bin error estimate we used the standard deviation because it describes the intrinsic PDS scatter. All PDS points within the bin were used for standard deviation calculation. The selection of binning frequency interval and number of divided light curve subsamples we based on χ^2_ red and on visual inspection of the binned PDS details. It was clear from the beginning that the binned PDS has a multicomponent shape. Therefore, for fitting we used a multicomponent broken power law consisting of 4 linear functions with three break frequencies. Finally, we divided the OM data into 5 and PN (as the most relevant for timing analysis) into 6 light curve subsamples, with higher binning resolution in the OM data (because of considerably larger count rate), in order to get a solution with χ^2_ red closer to 1. For fitting, we used the Gnuplot[http://www.gnuplot.info/] software, yielding best-fit values of break frequencies with standard errors calculated from the variance-covariance matrix.

It is worth noting, that we are comparing quantitatively very different data, i.e.  PDSs averaged from 5 one-day light curve segments, vs. 5/6 segments of  data with much shorter duration. Therefore, our investigation of the  observation is not adequate for detailed studies of low frequencies and thus we concentrate mainly on the break frequencies L_1 and L_2.

Resulting best cases with PDS, with binned data and fitted broken power laws are depicted in Fig. <ref>. The resulting fitted PDS parameters with standard errors are summarized in Table <ref>.






The break frequency L_1 can clearly be seen in both PN and OM PDSs. All fitted values of the searched break frequency (Table <ref>) agree well with the observed L_1 frequency (Table <ref>) within the errors. Furthermore, a plateau is seen between L_1 and L_2 in both PDSs. For OM data the power decrease toward lower frequencies starts at a break consistent with the  L_2 value. For the PN PDS this break is slightly higher, but all values (PN, OM,  L_2) agree within the errors, albeit PN and L_2 agree only marginally (the error extrema are equal). An additional power increase and possible break close to  L_3 value is seen in the OM data, although still speculative, but with no indication in PN PDS.










To investigate the low frequency part of the PDS, another PDS calculation with a smaller number of light curve subsamples is needed. However, this yields a much more scattered PDS with larger error bars. A lower frequency bin would reduce the errorbars, but also increases the scatter in the PDS, which is counterproductive. In general, this case appeared to not be adequate for a quantitative study and fitting.

In the last solution we concentrate on the high frequency part by increasing the number of light curve subsamples. This excludes low frequencies from the PDS, but details in the remaining PDS should be clearer with lower scatter, allowing shorter/finer frequency bins also for the PN PDS. Fig. <ref> shows a 10 subsample case where the presence of a break around L_2 is clearer in the PN PDS. The values of χ^2_ red are higher (2.0 in the OM PDS) or much higher (10.4 in the PN PDS) than 1 because of smaller errorbars as a natural result of more subsamples used for PDS calculation reducing the scatter.




Based on the results from the spectral analysis (Fig. <ref>), we created two light curves for a hard and a soft energy band in order to investigate whether the different radiation sources may show different variability patterns. The energy interval for the soft light curve extraction was chosen to roughly agree with the lower temperature plasma component (150 - 1500 eV) while the hard light curve was extracted from 1500 - 10000 eV. Non-averaged and non-binned PDSs are depicted in Fig. <ref>. Interestingly, the orbital modulation can only be seen in the soft band while it is absent in the hard X-rays. Averaged and binned PDSs with fits are depicted in Fig. <ref>. For direct comparison we keep the same binning resolution in both cases despite the fact that only one fit yields acceptable χ^2_ red. Both bands show a clear break frequency coincident with the  L_2 value, while the break is more pronounced in the hard band. The different PDS behaviour of both bands between L_1 and L_2 is probably the reason for the slightly higher value obtained from the full PN light curve, i.e. the integrated PDS is slightly deformed by the asymmetry in both bands, while the break agrees well with L_2. The fitted PDS parameter values are given in Table <ref>.







Finally we can conclude that except for the dominant break around the frequency L_1, the OM PDS shows an obvious L_2 and a possible L_3 break, also detected in  data, while the PN PDS only shows a break consistent with  L_2 value.



 §.§ PDS simulations


While the presence of the L_1 break frequency in  data is unambiguous, the second break at L_2 requires some more justification. For this purpose we simulated artificial light curves following the method by <cit.>. The principle is to use two gaussian-distributed random numbers following the input PDS, and use them as the real and imaginary part of the Fourier coefficient. This is done for every frequency and the required time series (synthetic light curve) is obtained by inverse Fourier transform. As input PDS we used PDS parameters from the fits to the real  data and the artificial light curves had the same duration and sampling as the original  data. We applied the same periodogram calculation, light curve subsample division and binning procedure as for the observations.

In Fig. <ref> we show various simulations compared to the PN binned PDS from the top panel of Fig. <ref>. We used three PDS models, i.e. a two-component model using two red noises before and after the break frequency at log(f/Hz) = -2.96, a three-component model, the same as the two component case but the lowest red noise is added, and finally a four-component model as shown by the fit in the top panel of Fig. <ref>. Therefore, only the latter model comprises the break frequency at log(f/Hz) = -3.28 (L_2 equivalent). The left column of Fig. <ref> shows the mean PDS calculated from 10000 simulations with standard deviations as assumed error. Apparently the four-component model is the best, although not perfect, i.e. the low-frequency part and the power depression below the L_2 break is slightly higher than the original fit used as input for the simulations. This suggests that the random process is important and not every simulated PDS from a light curve divided into 6 subsamples yields the L_2 break frequency. In the left column of Fig. <ref> we display examples of simulated PDSs close to the original fit. Apparently, the power depression below the L_2 break frequency can be the result of a random process without the L_2 break frequency inherently rooted in the PDS. However, these selected cases in the two and three-component models emerge with lower probability than in the four-component case. Therefore, we can conclude that it is not certain that the break frequency L_2 in the observed PN PDS is a real feature, but it is certainly more probable than that it is just a result of a random process.




The presented PDS simulations are ideal for PDS parameters uncertainty derivation directly from the randomness of the red noise. However, the lower break frequency is not present in all simulated PDSs and the fitting procedure did not converge properly in the majority of cases and direct access to the initial parameter estimate was required. Therefore, repeating the simulation process 10000 times, the uncertainty derivation can only be realised for the highest break frequency. We performed simple broken power law fits to the binned PDSs for frequencies higher than log(f/Hz) = -3.2. A Gaussian fit to the resulting histogram yields a mean value of -3.00 and 1-σ parameter of 0.12 as break frequency error. This uncertainty value is very close to the standard error of 0.14 from Table <ref>. This suggests that the standard errors derived by the Gnuplot software are good error estimates.



 §.§ Time delays


The UV and X-ray light curves shown in Fig. <ref> seem to be correlated and in this section we search for time lags between UV and X-rays using a cross-correlation function (CCF) as described in Section 2 of <cit.>. Because of the non-continuous nature of the OM light curve, we calculated a separate CCF for each of the 13 continuous OM light curves (Fig. <ref> in left panel). Subsequently we calculated a mean CCF which is depicted as a shaded area (the error of the mean is used as an estimate of the error in the CCF) in the right panel of Fig. <ref>.

The mean CCF has a maximum around zero time shift but is clearly asymmetric indicating that not all variations are simultaneous or constantly shifted[A constant time lag would result in a symmetric shape but with a shift of the peak.]. The same indication of the positive time lags (asymmetry or shift of the peak) is seen in almost all individual CCFs except 3 cases, where the indicated time lag is reverse (marked as thick lines in left panel of Fig. <ref>). We fitted the mean CCF with a single Lorentzian yielding a poor fit (Fig. <ref>, right panel, thin line), while a sum of two Lorentzians yielded an acceptable solution (Fig. <ref>, right panel, thick line). The exact values (time lags) of Lorentzian centers depends on the time lag interval over which we performed the fitting. Some time lags are summarized in Table <ref>[For time lags larger than 400 s the CCF is loosing the Lorentzian shape and the wings become constant.].







Clearly it is not straight forward to derive the time lag from the CCF, but it appears that part of the signal is (almost) simultaneous and part of the X-ray lags behind the UV. Therefore, we performed a more detailed frequency dependent time lag and coherence analysis following the method described by <cit.>. The two compared light curves must have the same sampling and number of data points. Therefore, we resampled the OM light curves by linear interpolation in order to synchronise them with the PN data[The difference is just a few seconds, which is (almost) negligible when comparing both the original and the resampled light curves.] and we filled the gaps with a linear function between two extreme data points of two consecutive OM observations. The coherence and time lags with 3-σ uncertainties are shown in Fig. <ref>. Clearly, there is a coherence peak around the L_1 break frequency, where the positive time lag, in spite of large error bars, confirms the conclusion from the CCF analysis. Also the time delay values are in good agreement with the ones derived from two-Lorentzian CCF fitting (t_ 2 value in Table. <ref>).




In Fig. <ref>, we visualize the behaviour of the PN lag. We selected some examples of UV rising before PN. In order to get comparable light curves, we rescale the selected light curve segments[Except the flare around 58 ks (Fig. <ref> and bottom right panel of Fig. <ref>) because there is a very strong peak in PN. Therefore, we choose the second highest peak.], transforming the count rates to get minima at 0 and maxima at 1. The results suggest that not the whole flare is offset toward earlier time, just the increasing branch shows this characteristic, while the decreasing parts of OM and PN flares are (almost) simultaneous. Therefore, the time delay derived from Figs. <ref> and <ref> is not the true time delay between the peaks, but it represents a time delay of the rising part only.




In order to study the asymmetry in more detail we analysed the time delay between the soft and hard X-ray bands. All corresponding frequency-dependent time lags are visualized in Fig. <ref>. It is clear that the soft X-ray emission lags behind the hard emission, while the soft emission lags behind the UV emission for frequencies between L_1 and L_2, with increasing time delay towards L_2. The delay between hard X-rays and UV is not as clear because the data are noisier and many values intersect the zero time lag within the error bars. It is possible that those data are (almost) simultaneous, but a PN lagging behind the OM is still possible mainly around the L_1 break frequency.



The situation is depicted in Fig. <ref> with two representative examples from Fig. <ref>. It is clear how the hard and soft PN count rates rise almost simultaneously, while the hard band is declining earlier. Therefore the mentioned asymmetry/lag between OM and PN is caused mainly by the soft X-ray band. Hard X-rays and UV can thus be more synchronized as suggestive from Fig. <ref>.




For a global study of the mentioned profile features, we analysed the flares in the same way as <cit.>, i.e. we averaged several flares in order to get a mean profile. Because of low photon statistics, only the 50 s binning was adequate for this analysis[Lower binning does not show smooth flares for the X-ray light curve and a central very thin spike is formed by the averaging instead of a flare.]. The peaks of the flares were defined as the local maxima with 10 data points before and 10 after the maximum in PN light curve, yielding 55 flares. The results are depicted in Fig. <ref> where the flares were divided by the underlying integrated area in order to make them comparable[Mainly because of OM flares which have much larger amplitude and count rate than the X-ray]. It is clear from the plots, that the OM is rising before the X-rays, while soft X-rays lag behind the hard band. The simultaneous rise of hard and soft X-rays as suggestive from Fig. <ref> is therefore not confirmed. The opposite seems to be present, but all lags are in agreement with the frequency-dependent analysis in Figs. <ref> and <ref>.




Finally, it seems that all peaks are simultaneous. Despite the fact, that lower binning is not adequate for this analysis because of low number of counts, we checked for this simultaneity and it appeared at 30 s binning. Lower binnings yield a small deviation, but not reliable because of X-ray average flare showing just a very narrow central spike[It is the result of a flare profile in the original light curve, where the flares are just an accumulation of individual photon spikes.].



§ RMS-FLUX RELATION


After the PDS study, we tested whether the  data satisfies another typical feature of the flickering in accreting systems, i.e. the rms-flux relation. For this purpose we used the OM light curve samples into 10 s bins, but for the PN case we needed larger bins of 50 s because of considerably lower count rate. The resulting rms-flux relation is shown in Fig. <ref>. The typical linear trend is shown by a linear fit to the data.






§ DISCUSSION



We present an analysis of  observations of the nova-like system MV Lyr in order to test the conclusion of <cit.> that the highest break frequency L_1 detected in optical  data () is generated by the central hot optically thin and geometrically thick corona. Such structure is radiating in X-rays, therefore the presence of the break frequency log(L_1) = -3.01 ± 0.06 Hz is expected in X-rays.



 §.§ Timing analysis


The presence of the searched break frequency around the  L_1 value can clearly be recognized in all forms of PDS we derived from the X-ray and UV data. The derived value is log(f) = -2.96 ± 0.14 Hz. Furthermore, we found significant indications of the presence of the second optical break frequency L_2 found in  data. This break frequency with a value of log(f) = -3.41 ± 0.10 Hz is seen in UV data from the OM instrument of . X-ray data from the PN detector yield only a marginally consistent value of log(f) = -3.28 ± 0.07 Hz, while division of the PN light curve into soft and hard bands yields a value of log(f) = -3.40 ± 0.09 Hz (from the hard band with smaller error and better χ^2_ red). This break frequency is more pronounced in the hard band.

Finally, we can mention a possible presence of L_3 in the OM data while it is not found in the PN data. This would suggest only an optical origin which is in agreement with the disc as a source (). But because of low resolution of the lower PDS limit, this conclusion still rests on weak grounds.



 §.§ Time delays



If the variability characterized by the break frequency L_1 is generated by the central corona, the signal must be seen first in X-rays and later reprocessed by the geometrically thin disc (). A time lag, where reprocessed optical/UV radiation lags behind the X-ray is expected. If we assume that such irradiation generates the reprocessed signal with maximal distance from the source to be the outermost boundary of the disc, the light travel time is the expected maximal time lag. The typical radii of discs in CVs are ∼ 10^10 cm, corresponding to a travel time of 0.3 s. With the precision of our  data we are not able to measure such a small time delay between OM and PN data. Therefore, any discussion about reprocessing time lag is unfortunately not possible.

However, we found some other delays in the  light curves. First is the X-ray lagging behind the UV signal between L_1 and L_2 frequency with a specific characteristics, i.e. X-ray data lagging behind the UV at the rising branch of a flare, while the decrease is simultaneous. Such profile is not new in binary light curves. It was observed by <cit.> in the X-ray binary Cyg X-1, where the hard band lagged behind the soft on the rising part of a flare, while the decrease was simultaneous.

Furthermore, we showed that the soft X-ray component lags behind the hard component, yielding an even stronger lag of soft X-rays behind the UV signal. The lag of hard X-rays behind the UV is not significant as in other band combinations, but still possible. Finally, we did not find any credible evidence of a lag of entire flares, i.e. the peaks seem to be simultaneous at least any lags are smaller than the bin size.

Because of the uncertain hard X-ray - UV behaviour, we examine three scenarios, i.e. 1) entire X-ray lagging behind the UV, 2) only soft X-ray lagging behind simultaneous UV and hard X-ray (neglecting the small difference), 3) and hard X-ray lagging behind UV, with soft X-ray lagging behind hard X-ray, where all three scenarios have simultaneous peaks.

How to explain the first solution? Let us assume that a turbulent element is formed somewhere in the disc. It needs some time to propagate towards the disc surface and to start to evaporate. Once the evaporation is feeding the corona, the locally enhanced plasma density generates enhanced X-ray emission. Once the turbulent element in the disc dissipates it stops to feed the corona and the response is an immediate decrease of X-ray emission, simultaneously with the turbulent eddy dissipation and local UV radiation. Therefore, the rise delay can be attributed to the turbulent eddy penetration towards the evaporation conditions. Such penetration can also be understood as an expansion which approaches the eddy to the disc surface. The idea of an expanding and cooling sphere, the so-called fireball model, was proposed by <cit.> and used by <cit.> to explain a typical color delay in CVs flickering. This mechanism does not agree with the interpretation by <cit.> of L_1 that the variability is generated by the corona, but it has also problems to explain the detected soft lag behind the hard X-ray.

The second solution is the simultaneous hard X-ray and UV. Two different regions would be responsible for these three radiation bands. A similar behaviour of UV and hard X-rays suggests a region where the hard X-rays are generated by free-free transitions of the hot ionised plasma and the UV is the reprocessing of X-rays by the disc or synchrotron radiation of the same plasma. The soft X-rays with the delay should be generated in a different region of the corona after propagation of the flare generating inhomogeneity/turbulent eddy to different temperature conditions.

Finally, the third scenario can be a combination of the two previous, i.e. an expanding eddy evaporating at the disc-corona transition, with the evaporated matter propagating into a lower temperature region, while the X-rays are reprocessed into UV during the overall process explaining the common decline.



 §.§ Spectral analysis



The  spectra are well fitted with a 2 temperature VAPEC or a mkcflow model. Both fits are statistically equivalent suggesting that our interpretation must include a cooler and a hotter source with several points requiring further discussion. For the VAPEC model, the hotter component would have a higher Fe abundance which for a partial absorber, applied only to the hot component, would not be necessary.

The first remark is based on a partial absorber that may be associated with the boundary layer between disc and the white dwarf (Fig. <ref>). The soft component may be interpreted as coming from the corona. In such case the detected fast variability can originate in the boundary layer instead of the corona, because the former is the source of 90% of the emission. At first view, this would disagree with the original assumption based on <cit.> modelling, that the observed  L_1 break frequency is generated in the corona.

However, <cit.> mention the boundary layer as the source of the cooling flow radiation[Note that the mkcflow model is more of an empiric description of boundary layer emission than a rigorous model and might not be suitable for all boundary layers of accretion discs.]. If this localisation is a necessary condition for the cooling flow, a modification is suggested, i.e. X-rays are not generated in the corona nor in the disc-white dwarf boundary layer, but in the boundary layer between the corona and the white dwarf (Fig. <ref>). Such boundary layer is expected, because of different tangential velocities of the corona and the central star. In this case the fast variability can still be generated in the corona by accretion inhomogeneities summed all the way down to the white dwarf surface, where it modulates the X-ray radiation of the boundary layer. Actually this is the principle of the multiplicative accretion process where the inhomogeneities are generated at different radii and transported inwards which modulates the final mass accretion rate (, , ).

The second remark is based on the mass accretion rate derived from the mkcflow model, of the order of 10^-12 M_⊙ y^-1. Such value is acceptable for a dwarf nova during quiescence () but hardly probable for a nova-like in the high state (see Fig. <ref>). Furthermore, if the hot component is associated with the boundary layer it should have a temperature of the order of 10^6 - 10^7 K. <cit.> studied MV Lyr in its high state using  UV data and derived a boundary layer temperature of 10^5 K which is strongly inconsistent with our finding of 10^7 K.

However, in a quiescent dwarf nova the inner disc is evaporating and forming a geometrically thick corona (). All the matter is accreted via this structure, and a UV - X-ray delay is observed. <cit.> proposed a sandwiched model for MV Lyr, where a geometrically thick corona is surrounding a geometrically thin disc up to a certain radius, but the geometrically thin disc reaches down to the white dwarf surface. In such a case, only a small fraction of the matter is evaporated and accreted to the central compact object via the corona. This idea is in agreement with our non-detection of a UV - X-ray delay and with the very low mass accretion rate of the order of 10^-12 M_⊙ y^-1. In such case the disc-white dwarf boundary layer can have the temperature of 10^5 K derived by <cit.>, but the corona-white dwarf boundary layer is much hotter with the temperature found in this paper. The latter is in agreement with the proposed hot boundary layer merged with an accretion-dominated flow (equivalent to corona in our paper) by <cit.>.

Therefore, both the cool and hot X-ray sources may either be generated in the corona where the hotter component extends closer/down to the white dwarf surface where it would be partially absorbed, or it is associated with the corona-white dwarf boundary layer. In both solutions the mass accretion fluctuations generating the  L_1 break frequency may be localised in the corona as suggested by <cit.>.

Unfortunately, we cannot distinguish between the case with and without the partial absorber as both models reproduce the  data equally well.

In order to compare our observation with the  data we used WebPIMMS[The WebPIMMS tool is available at https://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3pimms/w3pimms.pl] and adopted the (very) simple black-body model of <cit.>. We estimated the 0.3-10 keV band luminosities of all three  observations to be 5.87 × 10^31 erg s^-1, 4.65 × 10^31 erg s^-1, and 1.45 < × 10^30 erg s^-1, being the  observed luminosity in high state at least a factor ≃ 8.3 larger than what obtained by using the  data alone and the 2T VAPEC model (see Table <ref>).

The results from the spectral analysis of the  data by <cit.> and <cit.> are inconsistent with the results presented here.
In particular, the models that fit the  data are not consistent with the high-state  spectrum. Since the  spectra were acquired in a much shorter exposure time than the  ones and because of the lower sensitivity, the signal to noise is considerable lower. Especially at high energies, the  spectra pose little constraint, and we have some doubt that the source characteristics have really changed during the high-state  spectrum and the new  observation. The low signal to noise could introduce random features that are not accounted for in the response
that can lead to inconsistent results when assuming the average response. Note that for , only 'canned' response files are used while for , a
customized response file is produced based on the most recent calibration, the location of the source on the CCD detector etc. We thus doubt that the source characteristics were really different during the times of observations and give more confidence to the results derived from the EPIC spectra.



 §.§ Comparison with past interpretations


Part of the X-ray flux is probably generated in a wind from the system (, ). This may seem like a contradicting scenario to the one presented in this paper, but in fact it is not. Our corona interpretation is based on a model where inefficient cooling evaporates the corona following <cit.>. <cit.> used archival  spectra and fitted them with a combination of cooling flow model with an additional black body component. The latter were not consistent with any optically thick boundary layer (suggested by ). The authors deduced that the boundary layer must be hot with inefficient cooling and an evaporated corona is expected. However, following <cit.> the evaporation yields two solutions, i.e. the material is partially accreted on to the white dwarf via a corona, but partially lost via a wind. Therefore, the existence of the corona is directly connected to a wind. Furthermore, X-ray observations of UX UMa (), an eclipsing VY Scl system, show an eclipsed hard (kT ∼ 5 keV) component and an uneclipsed soft component. In order to produce eclipses, the hard component must thus be generated close enough to the white dwarf, while the soft component is generated further away from the system to escape eclipses. Moreover, the authors identified the boundary layer as the source of hard X-rays.

MV Lyr is a low inclination system without eclipses. If the UV is radiated by the disc, the non detection of the orbital modulation seen in Fig. <ref> is natural. If the hard X-rays are generated in a corona or corona-white dwarf boundary layer, the secondary is not big enough to generate any modulation, but an extended corona in the form of a wind can reach far enough away to be partially occulted. This explains well the presence/absence of orbital modulations of hard/soft X-rays in Fig. <ref>.

<cit.> folded the X-ray light curve of MV Lyr on the orbital period. They did the same for the whole energy interval and using two different intervals for soft (0.3 - 1.0 keV) and hard (2.5 - 7.0 keV) light curves. The authors did not detect significant changes in the light curve shape given the statistical errors of the low count rates, and concluded no significant energy dependence in orbital variation. Our modulation finding in Figs. <ref> and <ref> suggests the opposite, which is highlighted by our phased and binned soft and hard[Any possible out of error patterns in hard PN data can be cause by the division of the light curve into two bands which does not exactly separate the spectral components. Therefore, some residual soft X-ray counts are still present in the hard light curve.]  data in Fig. <ref> using an orbital period of 3.19 h by <cit.>.




Finally, <cit.> pointed out a dilemma in radiation from nova like systems, i.e. the optical and UV accretion rates and luminosities resembling those of dwarf novae in outburst, while their X-ray analysis during the same brightness state with resulting accretion rates and luminosities resemble those in quiescent dwarf novae. The sandwiched model offers a connection between these antagonistic aspects (Fig. <ref>), i.e. the existence of both, the standard disc with an optically thick boundary layer (resembling dwarf novae in outburst, proposed by ), together with a corona with an optically thin hot boundary layer (resembling quiescent dwarf novae, proposed by ).



 §.§ Proposed model


Finally a complex model of the multicomponent PDS shape in MV Lyr, that appears the most plausible one to us is illustrated in Fig. <ref>. <cit.> proposed that the highest break frequency L_1 detected in Kepler data is generated by a geometrically thick corona, while <cit.> found agreement with L_3 and the whole geometrically thin disc (the closer to the white dwarf, the more flares are produced) with enhanced activity of the outer disc rim causing the break frequency L_4.




However, despite the fact that the model by <cit.> is based on a geometrically thin disc, they found the same solution with disc parameters consistent with the geometrically thick disc solution by <cit.> for the frequency L_1. <cit.> used a statistical method developed by <cit.> based on various simplifications, while <cit.> used a proper physical method developed by <cit.>. Therefore, the geometrically thick nature of the disc solution is more relevant. Besides a coincidence, this can yield an important suggestion, that the flare statistics of a geometrically thin disc is similar to the geometrically thick disc, i.e. the thin disc behaves in a similar way as the corona, and the angular momentum redistribution/gradient, turbulence dimension scales (see  for details) are similar.

The only problematic break frequency is L_2. The models of <cit.> yield L_3 and L_4 with α = 0.1 - 0.4 and an outer disc radius of 0.5 and 0.9 times the primary Roche lobe. L_1 was simulated with a small outer disc radius of 10^10 cm with α = 0.9 as found by <cit.>. But there is a so far ignored solution for L_2 in Fig. 4 of <cit.> as well, i.e. an outer disc radius equal to the coronal radius and α = 0.3 - 0.4. These values of α are consistent with solutions for L_3 and L_4. In RU Peg the low X-ray break frequency was detected also in UV (observed by ) and is interpreted as fluctuations of the inner disc transported to the boundary layer or corona (), i.e. every fluctuation of mass transfer at the inner disc modulates the UV light and X-ray, because the inner disc is feeding the corona or boundary layer, and every inhomogeneity causes inhomogeneous evaporation or feed of the boundary layer. A similar phenomenon can be applied to MV Lyr, i.e. the inhomogeneous mass accretion rate in the geometrically thin disc below the corona modulates the UV and causes inhomogeneous evaporation and subsequent X-ray radiation. The necessary condition is the detection of L_2 in UV and X-ray, which seems to be fulfilled, and this speculation would agree with the detected X-ray lag behind the UV (eddy propagation mechanism described in Section <ref>). But the inner disc fluctuations as turbulences were already accounted for in modelling the L_3 frequency. Therefore, the inner disc variability should be somehow enhanced to generate an additional component L_2. Perhaps the evaporated inhomogeneity is first radiating in X-rays and subsequently they are reprocessed by the geometrically thin disc, which causes the enhanced optical variability of the inner disc below the corona. Perhaps the overall X-ray vs. UV behaviour is a combination of two scenarios; 1) X-ray variability generated by accretion fluctuations in the corona, subsequently reprocessed into optical by the disc, 2) inner thin disc fluctuations radiating in UV and evaporating to the corona where it radiates in X-rays generating additional reprocessed optical signal. Therefore, two processes can be at work leading to the detected complicated band behaviour.

Moreover, a solution where the geometrically thin disc below the corona has enhanced activity because of disc-corona interaction is attractive as well. Different densities and different physical conditions of the disc and corona can yield different tangential velocities of both, and such a velocity gradient is a necessary condition for the Kelvin-Helmholtz instability. Therefore, such transition region can be highly turbulent. A similar interpretation was proposed by <cit.> in the case of V1504 Cyg and V344 Lyr, where a break frequency in the optical PDS is seen during outburst, but absent during quiescence. The disc is developed down to the white dwarf surface in the former case, while it is absent in the latter. The disc-corona interaction during the outburst is proposed to be more turbulent and generating the additional PDS component.



§ SUMMARY



In this paper we analyse new X-ray and UV data of the nova-like system MV Lyr taken by . The results can be summarized as follows.

(i) The X-ray spectra are well fitted with a 2 temperature VAPEC or a mkcflow model. The VAPEC would require different Fe abundance for each temperature component.

(ii) A partial absorber model, applied only to the hotter component, tests the possibility that the hotter component is generated in the boundary layer, while the cooler component originates from the inner geometrically thick corona. The Fe abundance would be the same for both components, giving strong preference to the partial absorber model.

(iii) The derived mass accretion rate and independent boundary layer (between disc and white dwarf) temperature measurements by <cit.> suggest that the hotter component is probably generated in the inner hot corona or a corona-white dwarf boundary layer. This supports the sandwiched model proposed by <cit.>.

(iv) We confirm the presence of the expected break frequency (log(f) = -3.06 ± 0.02 Hz) around the optical signal at log(f) = -3.01 ± 0.06 Hz detected in  data (, ). An UV equivalent from  data at log(f) = -3.08 ± 0.03 Hz is detected as well.

(v) The second optical break frequency from  data of log(f) = -3.39 ± 0.04 is present in OM (log(f) = -3.34 ± 0.04 Hz) and PN data (log(f) = -3.40 ± 0.09 Hz), while it is more pronounced in the harder PN band after light curve decomposition of the X-ray data to soft and hard band based on the two spectral components.

(vi) Soft X-rays are lagging behind both the UV and hard X-rays. The two latter are (almost) simultaneous. Without this decomposition, the entire X-ray flares lag behind the UV, but with probably simultaneous peaks, similar to what is observed by <cit.> in Cyg X-1. The lags are most pronounced during the rising branch of the flares with an almost simultaneous decline.

(vii) We confirm the sandwiched model suggested by <cit.>, i.e. a central hot geometrically thick optically thin disc is surrounding a geometrically thin and optically thick disc up to a distance of approximately 10^10 cm from the central white dwarf. The fastest variability around log(f) = -3 Hz is generated by accretion fluctuations in the corona and radiated as X-rays. The optical signal is generated by reprocessing of these X-rays by the geometrically thin disc.

(viii) The sandwiched model explains the contradictory resemblance of VY Scl systems to dwarf novae, i.e. the corona-white dwarf boundary layer resembles dwarf novae in quiescence, while the underlying geometrically thin disc-white dwarf boundary layer resembles dwarf novae in outburst.

(ix) A complex model is summarized based on additional work of <cit.>, where the other three break frequencies detected in  data come from the outer disc rim, the whole accretion disc and possibly (just a speculation) from the inner disc region below the corona.



§ ACKNOWLEDGEMENT


AD was supported by the Slovak grant VEGA 1/0335/16 and by the ERDF - Research and Development Operational Programme under the project "University Scientific Park Campus MTF STU - CAMBO" ITMS: 26220220179. AAN acknowledges the support by the INFN project TAsP. Furthermore, we acknowledge with thanks the variable star observations from the AAVSO International Database contributed by observers worldwide and used in this research. We thank also Koji Mukai and Polina Zemko for advice and help in spectral modelling, and the anonymous reviewer for valuable comments.

mn2e




