Roche tomography of cataclysmic variables - V. A high-latitude star-spot on RU Pegasi

We present Roche tomograms of the secondary star in the dwarf nova system RU Pegasi derived from blue and red arm ISIS data taken on the 4.2-m William Herschel Telescope. We have applied the entropy landscape technique to determine the system parameters and obtained component masses of M1 = 1.06 Msun, M2 = 0.96 Msun, an orbital inclination angle of i = 43 degrees, and an optimal systemic velocity of gamma = 7 km/s. These are in good agreement with previously published values. Our Roche tomograms of the secondary star show prominent irradiation of the inner Lagrangian point due to illumination by the disc and/or bright spot, which may have been enhanced as RU Peg was in outburst at the time of our observations.We find that this irradiation pattern is axi-symmetric and confined to regions of the star which have a direct view of the accretion regions. This is in contrast to previous attempts to map RU Peg which suggested that the irradiation pattern was non-symmetric and extended beyond the terminator. We also detect additional inhomogeneities in the surface distribution of stellar atomic absorption that we ascribe to the presence of a large star-spot. This spot is centred at a latitude of about 82 degrees and covers approximately 4 per cent of the total surface area of the secondary. In keeping with the high latitude spots mapped on the cataclysmic variables AE Aqr and BV Cen, the spot on RU Peg also appears slightly shifted towards the trailing hemisphere of the star. Finally, we speculate that early mapping attempts which indicated non-symmetric irradiation patterns which extended beyond the terminator of CV donors could possibly be explained by a superposition of symmetric heating and a large spot.

Rapid Keck Spectroscopy of Cataclysmic Variables

We present an analysis of Rapid Keck Spectroscopy of the CVs AM Her (polar) and SS Cyg (dwarf nova). We decompose the spectra into constant and variable components and identify different types of variability in AM Her with different characteristic timescales. The variable flickering component of the accretion disc flux and the observational characteristics of a small flare in SS Cyg are isolated.

A Time-Series Analysis of the H alpha Emission Line in V3885 Sagitarii

Flickering is a phenomenon related to mass accretion observed among many classes of astrophysical objects. In this paper we present a study of flickering emission lines and the continuum of the cataclysmic variable V3885 Sgr. The flickering behavior was first analyzed through statistical analysis and the power spectra of lightcurves. Autocorrelation techniques were then employed to estimate the flickering timescale of flares. A cross-correlation study between the line and its underlying continuum variability is presented. The cross-correlation between the photometric and spectroscopic data is also discussed. Periodograms, calculated using emission-line data, show a behavior that is similar to those obtained from photometric datasets found in the literature, with a plateau at lower frequencies and a power-law at higher frequencies. The power-law index is consistent with stochastic events. The cross-correlation study indicates the presence of a correlation between the variability on Ha and its underlying continuum. Flickering timescales derived from the photometric data were estimated to be 25 min for two lightcurves and 10 min for one of them. The average timescales of the line flickering is 40 min, while for its underlying continuum it drops to 20 min.

Roche tomography of the secondary stars in CVs

The secondary stars in cataclysmic variables (CVs) are key to our understanding of the origin, evolution and behaviour of this class of interacting binary. In seeking a fuller understanding of these objects, the challenge for observers is to obtain images of the secondary star. This goal can be achieved through Roche tomography, an indirect imaging technique that can be used to map the Roche-lobe-filling secondary. The review begins with a description of the basic principles that underpin Roche tomography, including methods for determining the system parameters. Finally, we conclude with a look at the main scientific highlights to date, including the first unambiguous detection of starspots on AE Aqr B, and consider the future prospects of this technique.

Roche Tomography of the Secondary Stars in CVS

The secondary stars in cataclysmic variables (CVs) are key to our understanding of the origin evolution and behaviour of this class of interacting binary. In seeking a fuller understanding of these objects the challenge for observers is to obtain images of the secondary star. This goal can be achieved through Roche tomography an indirect imaging technique that can be used to map the Roche-lobe-filling secondary star. The review begins with a description of the basic principles that underpin Roche tomography including methods for determining the binary system parameters. Noise propagation onto Roche tomograms is also covered. Finally the review concludes with a look at the main scientific highlights to date and the future prospects for Roche tomography

Study on period search methods of variable stars: Application to ASAS and CRTS databases

Study on variable stars is an important topic of modern astrophysics. After the invention of powerful telescopes and high resolving powered CCD’s, the variable star data is accumulating in the order of peta-bytes. The huge amount of data need lot of automated methods as well as human experts. This thesis is devoted to the data analysis on variable star’s astronomical time series data and hence belong to the inter-disciplinary topic, Astrostatistics. For an observer on earth, stars that have a change in apparent brightness over time are called variable stars. The variation in brightness may be regular (periodic), quasi periodic (semi-periodic) or irregular manner (aperiodic) and are caused by various reasons. In some cases, the variation is due to some internal thermo-nuclear processes, which are generally known as intrinsic vari- ables and in some other cases, it is due to some external processes, like eclipse or rotation, which are known as extrinsic variables. Intrinsic variables can be further grouped into pulsating variables, eruptive variables and flare stars. Extrinsic variables are grouped into eclipsing binary stars and chromospheri- cal stars. Pulsating variables can again classified into Cepheid, RR Lyrae, RV Tauri, Delta Scuti, Mira etc. The eruptive or cataclysmic variables are novae, supernovae, etc., which rarely occurs and are not periodic phenomena. Most of the other variations are periodic in nature. Variable stars can be observed through many ways such as photometry, spectrophotometry and spectroscopy. The sequence of photometric observa- xiv tions on variable stars produces time series data, which contains time, magni- tude and error. The plot between variable star’s apparent magnitude and time are known as light curve. If the time series data is folded on a period, the plot between apparent magnitude and phase is known as phased light curve. The unique shape of phased light curve is a characteristic of each type of variable star. One way to identify the type of variable star and to classify them is by visually looking at the phased light curve by an expert. For last several years, automated algorithms are used to classify a group of variable stars, with the help of computers. Research on variable stars can be divided into different stages like observa- tion, data reduction, data analysis, modeling and classification. The modeling on variable stars helps to determine the short-term and long-term behaviour and to construct theoretical models (for eg:- Wilson-Devinney model for eclips- ing binaries) and to derive stellar properties like mass, radius, luminosity, tem- perature, internal and external structure, chemical composition and evolution. The classification requires the determination of the basic parameters like pe- riod, amplitude and phase and also some other derived parameters. Out of these, period is the most important parameter since the wrong periods can lead to sparse light curves and misleading information. Time series analysis is a method of applying mathematical and statistical tests to data, to quantify the variation, understand the nature of time-varying phenomena, to gain physical understanding of the system and to predict future behavior of the system. Astronomical time series usually suffer from unevenly spaced time instants, varying error conditions and possibility of big gaps. This is due to daily varying daylight and the weather conditions for ground based observations and observations from space may suffer from the impact of cosmic ray particles. Many large scale astronomical surveys such as MACHO, OGLE, EROS, xv ROTSE, PLANET, Hipparcos, MISAO, NSVS, ASAS, Pan-STARRS, Ke- pler,ESA, Gaia, LSST, CRTS provide variable star’s time series data, even though their primary intention is not variable star observation. Center for Astrostatistics, Pennsylvania State University is established to help the astro- nomical community with the aid of statistical tools for harvesting and analysing archival data. Most of these surveys releases the data to the public for further analysis. There exist many period search algorithms through astronomical time se- ries analysis, which can be classified into parametric (assume some underlying distribution for data) and non-parametric (do not assume any statistical model like Gaussian etc.,) methods. Many of the parametric methods are based on variations of discrete Fourier transforms like Generalised Lomb-Scargle peri- odogram (GLSP) by Zechmeister(2009), Significant Spectrum (SigSpec) by Reegen(2007) etc. Non-parametric methods include Phase Dispersion Minimi- sation (PDM) by Stellingwerf(1978) and Cubic spline method by Akerlof(1994) etc. Even though most of the methods can be brought under automation, any of the method stated above could not fully recover the true periods. The wrong detection of period can be due to several reasons such as power leakage to other frequencies which is due to finite total interval, finite sampling interval and finite amount of data. Another problem is aliasing, which is due to the influence of regular sampling. Also spurious periods appear due to long gaps and power flow to harmonic frequencies is an inherent problem of Fourier methods. Hence obtaining the exact period of variable star from it’s time series data is still a difficult problem, in case of huge databases, when subjected to automation. As Matthew Templeton, AAVSO, states “Variable star data analysis is not always straightforward; large-scale, automated analysis design is non-trivial”. Derekas et al. 2007, Deb et.al. 2010 states “The processing of xvi huge amount of data in these databases is quite challenging, even when looking at seemingly small issues such as period determination and classification”. It will be beneficial for the variable star astronomical community, if basic parameters, such as period, amplitude and phase are obtained more accurately, when huge time series databases are subjected to automation. In the present thesis work, the theories of four popular period search methods are studied, the strength and weakness of these methods are evaluated by applying it on two survey databases and finally a modified form of cubic spline method is intro- duced to confirm the exact period of variable star. For the classification of new variable stars discovered and entering them in the “General Catalogue of Vari- able Stars” or other databases like “Variable Star Index“, the characteristics of the variability has to be quantified in term of variable star parameters.

SDSS J150722.30+523039.8: a cataclysmic variable formed directly from a detached white dwarf/brown dwarf binary?

We present high-speed, three-colour photometry of the eclipsing cataclysmic variable SDSS J150722.30+523039.8 (hereafter SDSS J1507). This system has an orbital period of 66.61 min, placing it below the observed `period minimum' for cataclysmic variables. We determine the system parameters via a parametrized model of the eclipse fitted to the observed lightcurve by ?2 minimization. We obtain a mass ratio of q = 0.0623 +/- 0.0007 and an orbital inclination . The primary mass is Mw = 0.90 +/- 0.01Msolar. The secondary mass and radius are found to be Mr = 0.056 +/- 0.001Msolar and Rr = 0.096 +/- 0.001Rsolar, respectively. We find a distance to the system of 160 +/- 10pc. The secondary star in SDSS J1507 has a mass substantially below the hydrogen burning limit, making it the second confirmed substellar donor in a cataclysmic variable. The very short orbital period of SDSS J1507 is readily explained if the secondary star is nuclearly evolved, or if SDSS J1507 formed directly from a detached white dwarf/brown dwarf binary. Given the lack of any visible contribution from the secondary star, the very low secondary mass and the low HeI ?6678/Ha emission-line ratio, we argue that SDSS J1507 probably formed directly from a detached white dwarf/brown dwarf binary. If confirmed, SDSS J1507 will be the first such system identified. The implications for binary star evolution, the brown dwarf desert and the common envelope phase are discussed.

The geometry and ionization structure of the wind in the eclipsing nova-like variables RW Tri and UX UMa

The UV spectra of nova-like variables are dominated by emission from the accretion disk, modified by scattering in a wind emanating from the disk. Here, we model the spectra of RW Tri and UX UMa, the only two eclipsing nova-like variables which have been observed with the Hubble Space Telescope in the far-ultraviolet, in an attempt to constrain the geometry and the ionization structure of their winds. Using our Monte Carlo radiative transfer code, we computed spectra for simply parameterized axisymmetric biconical outflow models and were able to find plausible models for both systems. These reproduce the primary UV resonance lines-N v, Si iv, and C iv-in the observed spectra in and out of eclipse. The distribution of these ions in the wind models is similar in both cases as is the extent of the primary scattering regions in which these lines are formed. The inferred mass-loss rates are 6%-8% of the mass accretion rates for the systems. We discuss the implication of our point models for our understanding of accretion disk winds in cataclysmic variables. © 2010. The American Astronomical Society. All rights reserved.

A J-band detection of the sub-stellar mass donor in SDSS J1433+1011

We present time-resolved J-band spectroscopy of the short-period cataclysmic variable SDSS J143317.78+101123.3. We detect absorption lines from the sub-stellar donor star in this system, which contributes 38 +/- 5 per cent to the J-band light. From the relative strengths of the absorption lines in the J band, we estimate the spectral type of the donor star to be L2 +/- 1. These data are the first spectroscopic detection of a donor with a confirmed sub-stellar mass in a cataclysmic variable, and the spectral type is consistent with that expected from semi-empirical evolutionary models.

Using skew mapping, we have been able to derive an estimate for the radial velocity of the donor of K-d = 520 +/- 60 km/s. This value is consistent with, though much less precise than, predictions from mass determinations found via photometric fitting of the eclipse light curves.

The Unusually Luminous Extragalactic Nova SN 2010U

We present observations of the unusual optical transient SN 2010U, including spectra taken 1.03 days to 15.3 days after maximum light that identify it as a fast and luminous Fe II type nova. Our multi-band light curve traces the fast decline (t 2 = 3.5 ± 0.3 days) from maximum light (MV = -10.2 ± 0.1 mag), placing SN 2010U in the top 0.5% of the most luminous novae ever observed. We find typical ejecta velocities of ≈1100 km s-1 and that SN 2010U shares many spectral and photometric characteristics with two other fast and luminous Fe II type novae, including Nova LMC 1991 and M31N-2007-11d. For the extreme luminosity of this nova, the maximum magnitude versus rate of decline relationship indicates a massive white dwarf (WD) progenitor with a low pre-outburst accretion rate. However, this prediction is in conflict with emerging theories of nova populations, which predict that luminous novae from massive WDs should preferentially exhibit an alternate spectral type (He/N) near maximum light.

Self-consistent modelling of line-driven hot-star winds with Monte Carlo radiation hydrodynamics

Radiative pressure exerted by line interactions is a prominent driver of outflows in astrophysical systems, being at work in the outflows emerging from hot stars or from the accretion discs of cataclysmic variables, massive young stars and active galactic nuclei. In this work, a new radiation hydrodynamical approach to model line-driven hot-star winds is presented. By coupling a Monte Carlo radiative transfer scheme with a finite volume fluid dynamical method, line-driven mass outflows may be modelled self-consistently, benefiting from the advantages of Monte Carlo techniques in treating multiline effects, such as multiple scatterings, and in dealing with arbitrary multidimensional configurations. In this work, we introduce our approach in detail by highlighting the key numerical techniques and verifying their operation in a number of simplified applications, specifically in a series of self-consistent, one-dimensional, Sobolev-type, hot-star wind calculations. The utility and accuracy of our approach are demonstrated by comparing the obtained results with the predictions of various formulations of the so-called CAK theory and by confronting the calculations with modern sophisticated techniques of predicting the wind structure. Using these calculations, we also point out some useful diagnostic capabilities our approach provides. Finally, we discuss some of the current limitations of our method, some possible extensions and potential future applications.

Simulations of angular momentum evolution in wind-fed CV precursors

It is thought that the secondary stars in cataclysmic variables (CVs) may undergo a period of mass loss in the form of a wind during the evolution of the system (Mullan et al. 1992). This wind is thought to magnetically brake the secondary star with a time-scale ~ 10^8 yr (e.g. van Paradijs 1986). When the secondary’s spin has been brought close to synchronism with the orbit it is possible for tidal torques to lock the secondary in synchronous rotation.

A two-armed pattern in flickering maps of the nova-like variable UU aquarii

We report the analysis of a uniform sample of 31 light curves of the nova-like variable UU Aqr with eclipse-mapping techniques. The data were combined to derive eclipse maps of the average steady-light component, the long-term brightness changes, and the low- and high-frequency flickering components. The long-term variability responsible for the ""low-brightness`` and ""high-brightness`` states is explained in terms of the response of a viscous disk to changes of 20%-50% in the mass transfer rate from the donor star. Low- and high-frequency flickering maps are dominated by emission from two asymmetric arcs reminiscent of those seen in the outbursting dwarf nova IP Peg, and they are similarly interpreted as manifestations of a tidally induced spiral shock wave in the outer regions of a large accretion disk. The asymmetric arcs are also seen in the map of the steady light aside from the broad brightness distribution of a roughly steady-state disk. The arcs account for 25% of the steady-light flux and are a long-lasting feature in the accretion disk of UU Aqr. We infer an opening angle of 10 degrees +/- 3 degrees for the spiral arcs. The results suggest that the flickering in UU Aqr is caused by turbulence generated after the collision of disk gas with the density-enhanced spiral wave in the accretion disk.

Estudo sistemático de estrelas variáveis na era dos grandes surveys

The interest in the systematic analysis of astronomical time series data, as well as development in astronomical instrumentation and automation over the past two decades has given rise to several questions of how to analyze and synthesize the growing amount of data. These data have led to many discoveries in the areas of modern astronomy asteroseismology, exoplanets and stellar evolution. However, treatment methods and data analysis have failed to follow the development of the instruments themselves, although much effort has been done. In present thesis, we propose new methods of data analysis and two catalogs of the variable stars that allowed the study of rotational modulation and stellar variability. Were analyzed the photometric databases fromtwo distinctmissions: CoRoT (Convection Rotation and planetary Transits) and WFCAM (Wide Field Camera). Furthermore the present work describes several methods for the analysis of photometric data besides propose and refine selection techniques of data using indices of variability. Preliminary results show that variability indices have an efficiency greater than the indices most often used in the literature. An efficient selection of variable stars is essential to improve the efficiency of all subsequent steps. Fromthese analyses were obtained two catalogs; first, fromtheWFCAMdatabase we achieve a catalog with 319 variable stars observed in the photometric bands Y ZJHK. These stars show periods ranging between ∼ 0, 2 to ∼ 560 days whose the variability signatures present RR-Lyrae, Cepheids , LPVs, cataclysmic variables, among many others. Second, from the CoRoT database we selected 4, 206 stars with typical signatures of rotationalmodulation, using a supervised process. These stars show periods ranging between ∼ 0, 33 to ∼ 92 days, amplitude variability between ∼ 0, 001 to ∼ 0, 5 mag, color index (J - H) between ∼ 0, 0 to ∼ 1, 4 mag and spectral type CoRoT FGKM. The WFCAM variable stars catalog is being used to compose a database of light curves to be used as template in an automatic classifier for variable stars observed by the project VVV (Visible and Infrared Survey Telescope for Astronomy) moreover it are a fundamental start point to study different scientific cases. For example, a set of 12 young stars who are in a star formation region and the study of RR Lyrae-whose properties are not well established in the infrared. Based on CoRoT results we were able to show, for the first time, the rotational modulation evolution for an wide homogeneous sample of field stars. The results are inagreement with those expected by the stellar evolution theory. Furthermore, we identified 4 solar-type stars ( with color indices, spectral type, luminosity class and rotation period close to the Sun) besides 400 M-giant stars that we have a special interest to forthcoming studies. From the solar-type stars we can describe the future and past of the Sun while properties of M-stars are not well known. Our results allow concluded that there is a high dependence of the color-period diagram with the reddening in which increase the uncertainties of the age-period realized by previous works using CoRoT data. This thesis provides a large data-set for different scientific works, such as; magnetic activity, cataclysmic variables, brown dwarfs, RR-Lyrae, solar analogous, giant stars, among others. For instance, these data will allow us to study the relationship of magnetic activitywith stellar evolution. Besides these aspects, this thesis presents an improved classification for a significant number of stars in the CoRoT database and introduces a new set of tools that can be used to improve the entire process of the photometric databases analysis

A possible origin of viscosity in Keplerian accretion disks due to secondary perturbation: Turbulent transport without magnetic fields

The origin of hydrodynamic turbulence in rotating shear flow is a long standing puzzle. Resolving it is especially important in astrophysics when the flow's angular momentum profile is Keplerian which forms an accretion disk having negligible molecular viscosity. Hence, any viscosity in such systems must be due to turbulence, arguably governed by magnetorotational instability, especially when temperature T greater than or similar to 10(5). However, such disks around quiescent cataclysmic variables, protoplanetary and star-forming disks, and the outer regions of disks in active galactic nuclei are practically neutral in charge because of their low temperature, and thus are not expected to be coupled with magnetic fields enough to generate any transport due to the magnetorotational instability. This flow is similar to plane Couette flow including the Coriolis force, at least locally. What drives their turbulence and then transport, when such flows do not exhibit any unstable mode under linear hydrodynamic perturbation? We demonstrate that the three-dimensional secondary disturbance to the primarily perturbed flow that triggers elliptical instability may generate significant turbulent viscosity in the range 0.0001 less than or similar to nu(t) less than or similar to 0.1, which can explain transport in accretion flows.