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.

Towards the main sequence: detailed analysis of weak line and post-T Tauri stars

A detailed study was performed for a sample of low-mass pre-main-sequence (PMS) stars, previously identified as weak-line T Tauri stars, which are compared to members of the Tucanae and Horologium Associations. Aiming to verify if there is any pattern of abundances when comparing the young stars at different phases, we selected objects in the range from 1 to 100 Myr, which covers most of PMS evolution. High-resolution optical spectra were acquired at European Southern Observatory and Observatorio do Pico dos Dias. The stellar fundamental parameters effective temperature and gravity were calculated by excitation and ionization equilibria of iron absorption lines. Chemical abundances were obtained via equivalent width calculations and spectral synthesis for 44 per cent of the sample, which shows metallicities within 0.5 dex solar. A classification was developed based on equivalent width of Li I 6708 angstrom and Ha lines and spectral types of the studied stars. This classification allowed a separation of the sample into categories that correspond to different evolutive stages in the PMS. The position of these stars in the Hertzsprung-Russell diagram was also inspected in order to estimate their ages and masses. Among the studied objects, it was verified that our sample actually contains seven weak-line T Tauri stars, three are Classical T Tauri, 12 are Fe/Ge PMS stars and 21 are post-T Tauri or young main-sequence stars. An estimation of circumstellar luminosity was obtained using a disc model to reproduce the observed spectral energy distribution. Most of the stars show low levels of circumstellar emission, corresponding to less than 30 per cent of the total emission.

SIMULATIONS OF WINDS OF WEAK-LINED T TAURI STARS. II. THE EFFECTS OF A TILTED MAGNETOSPHERE AND PLANETARY INTERACTIONS

Based on our previous work, we investigate here the effects on the wind and magnetospheric structures of weak-lined T Tauri stars due to a misalignment between the axis of rotation of the star and its magnetic dipole moment vector. In such a configuration, the system loses the axisymmetry presented in the aligned case, requiring a fully three-dimensional (3D) approach. We perform 3D numerical magnetohydrodynamic simulations of stellar winds and study the effects caused by different model parameters, namely the misalignment angle theta(t), the stellar period of rotation, the plasma-beta, and the heating index.. Our simulations take into account the interplay between the wind and the stellar magnetic field during the time evolution. The system reaches a periodic behavior with the same rotational period of the star. We show that the magnetic field lines present an oscillatory pattern. Furthermore, we obtain that by increasing theta(t), the wind velocity increases, especially in the case of strong magnetic field and relatively rapid stellar rotation. Our 3D, time-dependent wind models allow us to study the interaction of a magnetized wind with a magnetized extrasolar planet. Such interaction gives rise to reconnection, generating electrons that propagate along the planet`s magnetic field lines and produce electron cyclotron radiation at radio wavelengths. The power released in the interaction depends on the planet`s magnetic field intensity, its orbital radius, and on the stellar wind local characteristics. We find that a close-in Jupiter-like planet orbiting at 0.05 AU presents a radio power that is similar to 5 orders of magnitude larger than the one observed in Jupiter, which suggests that the stellar wind from a young star has the potential to generate strong planetary radio emission that could be detected in the near future with LOFAR. This radio power varies according to the phase of rotation of the star. For three selected simulations, we find a variation of the radio power of a factor 1.3-3.7, depending on theta(t). Moreover, we extend the investigation done in Vidotto et al. and analyze whether winds from misaligned stellar magnetospheres could cause a significant effect on planetary migration. Compared to the aligned case, we show that the timescale tau(w) for an appreciable radial motion of the planet is shorter for larger misalignment angles. While for the aligned case tau(w) similar or equal to 100 Myr, for a stellar magnetosphere tilted by theta(t) = 30 degrees, tau(w) ranges from similar to 40 to 70 Myr for a planet located at a radius of 0.05 AU. Further reduction on tau(w) might occur for even larger misalignment angles and/or different wind parameters.

SIMULATIONS OF WINDS OF WEAK-LINED T TAURI STARS: THE MAGNETIC FIELD GEOMETRY AND THE INFLUENCE OF THE WIND ON GIANT PLANET MIGRATION

By means of numerical simulations, we investigate magnetized stellar winds of pre-main-sequence stars. In particular, we analyze under which circumstances these stars will present elongated magnetic features (e.g., helmet streamers, slingshot prominences, etc). We focus on weak-lined T Tauri stars, as the presence of the tenuous accretion disk is not expected to have strong influence on the structure of the stellar wind. We show that the plasma-beta parameter (the ratio of thermal to magnetic energy densities) is a decisive factor in defining the magnetic configuration of the stellar wind. Using initial parameters within the observed range for these stars, we show that the coronal magnetic field configuration can vary between a dipole-like configuration and a configuration with strong collimated polar lines and closed streamers at the equator (multicomponent configuration for the magnetic field). We show that elongated magnetic features will only be present if the plasma-beta parameter at the coronal base is beta(0) << 1. Using our self-consistent three-dimensional magnetohydrodynamics model, we estimate for these stellar winds the timescale of planet migration due to drag forces exerted by the stellar wind on a hot-Jupiter. In contrast to the findings of Lovelace et al., who estimated such timescales using the Weber and Davis model, our model suggests that the stellar wind of these multicomponent coronae are not expected to have significant influence on hot-Jupiters migration. Further simulations are necessary to investigate this result under more intense surface magnetic field strengths (similar to 2-3 kG) and higher coronal base densities, as well as in a tilted stellar magnetosphere.

A possible stellar metallic enhancement in post-T Tauri stars by a planetesimal bombardment

The photospheres of stars hosting planets have larger metallicity than stars lacking planets. This could be the result of a metallic star contamination produced by the bombarding of hydrogen-deficient solid bodies. In the present work we study the possibility of an earlier metal enrichment of the photospheres by means of impacting planetesimals during the first 20-30 Myr. Here we explore this contamination process by simulating the interactions of an inward migrating planet with a disc of planetesimal interior to its orbit. The results show the percentage of planetesimals that fall on the star. We identified the dependence of the planet's eccentricity (e(p)) and time-scale of migration (tau) on the rate of infalling planetesimals. For very fast migrations (tau= 10(2) and 10(3) yr) there is no capture in mean motion resonances, independently of the value of e(p). Then, due to the planet's migration the planetesimals suffer close approaches with the planet and more than 80 per cent of them are ejected from the system. For slow migrations (tau= 10(5)and 10(6) yr) the percentage of collisions with the planet decreases with the increase of the planet's eccentricity. For e(p) = 0 and 0.1 most of the planetesimals were captured in the 2:1 resonance and more than 65 per cent of them collided with the star. Whereas migration of a Jupiter mass planet to very short pericentric distances requires unrealistic high disc masses, these requirements are much smaller for smaller migrating planets. Our simulations for a slowly migrating 0.1 M-Jupiter planet, even demanding a possible primitive disc three times more massive than a primitive solar nebula, produces maximum [Fe/H] enrichments of the order of 0.18 dex. These calculations open possibilities to explain hot Jupiter exoplanet metallicities.

Discovery of "isolated" co-moving T Tauri stars in Cepheus

Context. During the course of a large spectroscopic survey of X-ray active late-type stars in the solar neighbourhood, we discovered four lithium-rich stars packed within just a few degrees on the sky. Although located in a sky area rich in CO molecular regions and dark clouds, the Cepheus-Cassiopeia complex, these very young stars are projected several degrees away from clouds in front of an area void of interstellar matter. As such, they are very good "isolated" T Tauri star candidates. Aims. We present optical observations of these stars conducted with 1-2 m class telescopes. We acquired high-resolution optical spectra as well as photometric data allowing us to investigate in detail their nature and physical parameters with the aim of testing the "runaway" and "in-situ" formation scenarios. Their kinematical properties are also analyzed to investigate their possible connection to already known stellar kinematic groups. Methods. We use the cross-correlation technique and other tools developed by us to derive accurate radial and rotational velocities and perform an automatic spectral classification. The spectral subtraction technique is used to infer chromospheric activity level in the Hα line core and clean the spectra of photospheric lines before measuring the equivalent width of the lithium absorption line. Results. Both physical (lithium content, chromospheric, and coronal activities) and kinematical indicators show that all stars are very young, with ages probably in the range 10-30 Myr. In particular, the spectral energy distribution of TYC4496-780-1 displays a strong near-and far-infrared excess, typical of T Tauri stars still surrounded by an accretion disc. They also share the same Galactic motion, proving that they form a homogeneous moving group of stars with the same origin. Conclusions. The most plausible explanation of how these "isolated" T Tauri stars formed is the "in-situ" model, although accurate distances are needed to clarify their connection with the Cepheus-Cassiopeia complex. The discovery of this loose association of "isolated" T Tauri stars can help to shed light on atypical formation processes of stars and planets in low-mass clouds.

Post T Tauri stars in the solar neighborhood: isolated or members of young associations and moving groups

Post T Tauri stars (PTTS) are late-type stars in the age range between 10 and 100 Myr filling the gap between T Tauri (TTs) and zero-age: main sequence phases. This period of evolution remains ambiguous and until now different studies of young stars have failed to find the numbers of PTTS that are expected. In the last years, some PTTS have been identified among the X-ray detected pre-main sequence stars in some star-forming regions. More recently, additional PTTS have been identified in young associations and moving groups (β Pic, TW Hya, Tucana/Horologium, and the AB Dor). However, many isolated PTTS still remain undiscovered. In this contribution, we compiled the PTTS previously identified in the literature, and identified new candidates using the information provided by the high resolution spectra obtained during our surveys of late-type stars possible members to young moving groups, FGK stars in the solar neighborhood, and RasTyc sample. To identify PTTS we applied an age-oriented definition using relative age indicators (Li abundance, chromospheric and coronal emission and the kinematics) as well as color-magnitude diagrams and pre-main sequence isochrones.

New Insights from Aperiodic Variability of Young Stars

Nearly all young stars are variable, with the variability traditionally divided into two classes: periodic variables and aperiodic or "irregular" variables. Periodic variables have been studied extensively, typically using periodograms, while aperiodic variables have received much less attention due to a lack of standard statistical tools. However, aperiodic variability can serve as a powerful probe of young star accretion physics and inner circumstellar disk structure. For my dissertation, I analyzed data from a large-scale, long-term survey of the nearby North America Nebula complex, using Palomar Transient Factory photometric time series collected on a nightly or every few night cadence over several years. This survey is the most thorough exploration of variability in a sample of thousands of young stars over time baselines of days to years, revealing a rich array of lightcurve shapes, amplitudes, and timescales.

I have constrained the timescale distribution of all young variables, periodic and aperiodic, on timescales from less than a day to ~100 days. I have shown that the distribution of timescales for aperiodic variables peaks at a few days, with relatively few (~15%) sources dominated by variability on tens of days or longer. My constraints on aperiodic timescale distributions are based on two new tools, magnitude- vs. time-difference (Δm-Δt) plots and peak-finding plots, for describing aperiodic lightcurves; this thesis provides simulations of their performance and presents recommendations on how to apply them to aperiodic signals in other time series data sets. In addition, I have measured the error introduced into colors or SEDs from combining photometry of variable sources taken at different epochs. These are the first quantitative results to be presented on the distributions in amplitude and time scale for young aperiodic variables, particularly those varying on timescales of weeks to months.

Roche tomography of cataclysmic variables - III. Star-spots on AE Aqr

We present a Roche tomography reconstruction of the secondary star in the cataclysmic variable AE Aqr. The tomogram reveals several surface inhomogeneities that are due to the presence of large, cool star-spots. In addition to a number of lower latitude spots, the maps also show the presence of a large, high-latitude spot similar to that seen in Doppler images of rapidly rotating isolated stars, and a relative paucity of spots at a latitude of 40 degrees. In total, we estimate that some 18 per cent of the Northern hemisphere of AE Aqr is spotted.

The masses of the cataclysmic variables AC Cancri and V363 Aurigae

We present time-resolved spectroscopy and photometry of the double-lined eclipsing cataclysmic variables AC Cnc and V363 Aur (=Lanning 10). There is evidence of irradiation on the inner hemisphere of the secondary star in both systems, which we correct for using a model that reproduces the observations remarkably well. We find the radial velocity of the secondary star in AC Cnc to be K-R=176+/-3 km s(-1) and its rotational velocity to be v sin i=135+/-3 km s(-1). From these parameters we obtain masses of M-1=0.76+/-0.03 M-circle dot for the white-dwarf primary and M-2=0.77+/-0.05 M-circle dot for the K2+/-1 V secondary star, giving a mass ratio of q=1.02+/-0.04. We measure the radial and rotational velocities of the G7+/-2V secondary star in V363 Aur to be K-R=168+/-5 km s(-1) and v sin i=143+/-5 km s(-1), respectively. The component masses of V363 Aur are M-1=0.90+/-0.06M(circle dot) and M-2=1.06+/-0.11 M-circle dot giving a mass ratio of q=1.17+/-0.07. The mass ratios for AC Cnc and V363 Aur fall within the theoretical limits for dynamically and thermally stable mass transfer. Both systems are similar to the SW Sex stars, exhibiting single-peaked emission lines with transient absorption features, high-velocity S-wave components and phase-offsets in their radial-velocity curves. The Balmer lines in V363 Aur show a rapid increase in flux around phase 0 followed by a rapid decrease, which we attribute to the eclipse of an optically thick region at the centre of the disc. This model could also account for the behaviour of other SW Sex stars where the Balmer lines show only a shallow eclipse compared to the continuum.

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 cataclysmic variables - I. Artefacts and techniques

Roche tomography is a technique used for imaging the Roche-lobe-filling secondary stars in cataclysmic variables (CVs). In order to interpret Roche tomograms correctly, one must determine whether features in the reconstruction are real, or the result of statistical or systematic errors. We explore the effects of systematic errors using reconstructions of simulated data sets, and show that systematic errors result in characteristic distortions of the final reconstructions that can be identified and corrected. In addition, we present a new method of estimating statistical errors on tomographic reconstructions using a Monte Carlo bootstrapping algorithm, and show this method to be much more reliable than Monte Carlo methods which 'jiggle' the data points in accordance with the size of their error bars.

The POINT-AGAPE survey - I. The variable stars in M31

For the purposes of identifying microlensing events, the POINT-AGAPE collaboration has been monitoring the Andromeda galaxy (M31) for three seasons (1999-2001) with the Wide Field Camera on the Isaac Newton Telescope. In each season, data are taken for one hour per night for roughly 60 nights during the six months that M31 is visible. The two 33 x 33 arcmin(2) fields of view straddle the central bulge, northwards and southwards. We have calculated the locations, periods and brightness of 35 414 variable stars in M31 as a by-product of the microlensing search. The variables are classified according to their period and brightness. Rough correspondences with classical types of variable star (such as Population I and II Cepheids, Miras and semiregular long-period variables) are established. The spatial distribution of Population I Cepheids is clearly associated with the spiral arms, while the central concentration of the Miras and long-period variables varies noticeably, the brighter and the shorter period Miras being much more centrally concentrated.