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.

Stability and dynamics of terrestrial planets in binary star systems: application on Alpha Centauri

Binary stars are frequent in the universe, with about 50% of the known main sequence stars being located at a multiple star system (Abt, 1979). Even though, they are universally thought as second rate sites for the location of exo-planets and the habitable zone, due to the difficulties of detection and high perturbation that could prevent planet formation and long term stability. In this work we show that planets in binary star systems can have regular orbits and remain on the habitable zone. We introduce a stability criterium based on the solution of the restricted three body problem and apply it to describe the short period planar and three-dimentional stability zones of S-type orbits around each star of the Alpha Centauri system. We develop as well a semi-analytical secular model to study the long term dynamics of fictional planets in the habitable zone of those stars and we verify that planets on the habitable zone would be in regular orbits with any eccentricity and with inclination to the binary orbital plane up until 35 degrees. We show as well that the short period oscillations on the semi-major axis is 100 times greater than the Earth's, but at all the time the planet would still be found inside the Habitable zone.

Brown dwarfs: At last filling the gap between stars and planets

Until the mid-1990s a person could not point to any celestial object and say with assurance that “here is a brown dwarf.” Now dozens are known, and the study of brown dwarfs has come of age, touching upon major issues in astrophysics, including the nature of dark matter, the properties of substellar objects, and the origin of binary stars and planetary systems.

Library of high-resolution UES echelle spectra of F, G, K and M field dwarf stars

We present a library of Utrecht echelle spectrograph (UES) observations of a sample of F, G, K and M field dwarf stars covering the spectral range from 4800 Å to 10600 Å with a resolution of 55000. These spectra include some of the spectral lines most widely used as optical and near-infrared indicators of chromospheric activity such as Hβ, Mg I b triplet, Na I D_1, D_2, He I D_3, Hα, and Ca II IRT lines, as well as a large number of photospheric lines which can also be affected by chromospheric activity. The spectra have been compiled with the aim of providing a set of standards observed at high-resolution to be used in the application of the spectral subtraction technique to obtain the active-chromosphere contribution to these lines in chromospherically active single and binary stars. This library can also be used for spectral classification purposes. A digital version with all the spectra is available via ftp and the World Wide Web (WWW) in both ASCII and FITS formats.

The supergiant B[e] star LHA 115-S 18 - binary and/or luminous blue variable?

Context. The mechanism by which supergiant (sg)B[e] stars support cool, dense dusty discs/tori and their physical relationship with other evolved, massive stars such as luminous blue variables is uncertain. Aims. In order to investigate both issues we have analysed the long term behaviour of the canonical sgB[e] star LHA 115-S 18. Methods. We employed the OGLE II-IV lightcurve to search for (a-)periodic variability and supplemented these data with new and historic spectroscopy. Results. In contrast to historical expectations for sgB[e] stars, S18 is both photometrically and spectroscopically highly variable. The lightcurve is characterised by rapid aperiodic ` aring' throughout the 16 years of observations. Changes in the high excitation emission line component of the spectrum imply evolution in the stellar temperature - as expected for luminous blue variables - although somewhat surprisingly, spectroscopic and photometric variability appears not to be correlated. Characterised by emission in low excitation metallic species, the cool circumstellar torus appears largely unaffected by this behaviour. Finally, in conjunction with intense, highly variable He ii emission, X-ray emission implies the presence of an unseen binary companion. Conclusions. S18 provides observational support for the putative physical association of (a subset of) sgB[e] stars and luminous blue variables. Given the nature of the circumstellar environment of S18 and that luminous blue variables have been suggested as SN progenitors, it is tempting to draw a parallel to the progenitors of SN1987A and SN2009ip. Moreover the likely binary nature of S18 strengthens the possibility that the dusty discs/tori that characterise sgB[e] stars are the result of binary-driven mass-loss; consequently such stars may provide a window on the short lived phase of mass-transfer in massive compact binaries.

The Rotational Evolution and Magnetospheric Emission of the Magnetic Early B-type Stars

How do the magnetic fields of massive stars evolve over time? Are their gyrochronological ages consistent with ages inferred from evolutionary tracks? Why do most stars predicted to host Centrifugal Magnetospheres (CMs) display no H$\alpha$ emission? Does plasma escape from CMs via centrifugal breakout events, or by a steady-state leakage mechanism? This thesis investigates these questions via a population study with a sample of 51 magnetic early B-type stars. The longitudinal magnetic field \bz~was measured from Least Squares Deconvolution profiles extracted from high-resolution spectropolarimetric data. New rotational periods $P_{\rm rot}$ were determined for 15 stars from \bz, leaving only 3 stars for which $P_{\rm rot}$ is unknown. Projected rotational velocities \vsini~were measured from multiple spectral lines. Effective temperatures and surface gravities were measured via ionization balances and line profile fitting of H Balmer lines. Fundamental physical parameters, \bz, \vsini, and $P_{\rm rot}$ were then used to determine radii, masses, ages, dipole oblique rotator model, stellar wind, magnetospheric, and spindown parameters using a Monte Carlo approach that self-consistently calculates all parameters while accounting for all available constraints on stellar properties. Dipole magnetic field strengths $B_{\rm d}$ follow a log-normal distribution similar to that of Ap stars, and decline over time in a fashion consistent with the expected conservation of fossil magnetic flux. $P_{\rm rot}$ increases with fractional main sequence age, mass, and $B_{\rm d}$, as expected from magnetospheric braking. However, comparison of evolutionary track ages to maximum spindown ages $t_{\rm S,max}$ shows that initial rotation fractions may be far below critical for stars with $M_*>10 M_\odot$. Computing $t_{\rm S,max}$ with different mass-loss prescriptions indicates that the mass-loss rates of B-type stars are likely much lower than expected from extrapolation from O-type stars. Stars with H$\alpha$ in emission and absorption occupy distinct regions in the updated rotation-magnetic confinement diagram: H$\alpha$-bright stars are found to be younger, more rapidly rotating, and more strongly magnetized than the general population. Emission strength is sensitive both to the volume of the CM and to the mass-loss rate, favouring leakage over centrifugal breakout.

Exploring the high frequency emission of radio loud X-ray binaries

We report millimetre-wave continuum observations of the X-ray binaries Cygnus X-3, SS 433, LSI+61 303, Cygnus X-1 and GRS 1915+105. The observations were carried out with the IRAM 30 m-antenna at 250 GHz (1.25 mm) from 1998 March 14 to March 20. These millimetre measurements are complemented with centimetre observations from the Ryle Telescope, at 15 GHz (2.0 cm) and from the Green Bank Interferometer at 2.25 and 8.3 GHz (13 and 3.6 cm). Both Cygnus X-3 and SS 433 underwent moderate flaring events during our observations, whose main spectral evolution properties are described and interpreted. A significant spectral steepening was observed in both sources during the flare decay, that is likely to be caused by adiabatic expansion, inverse Compton and synchrotron losses. Finally, we also report 250 GHz upper limits for three additional undetected X-ray binary stars: LSI+65 010, LSI+61 235 and X Per.

Dynamically tuned interferometers for the observation of gravitational waves from coalescing compact binaries

We propose a new method of operating laser interferometric gravitational-wave detectors when observing chirps of gravitational radiation from coalescing compact binary stars. This technique consists of the use of narrow-band dual recycling to increase the signal but with the tuning frequency of the detector arranged to follow the frequency of a chirp. We consider the response of such an instrument to chirps, including the effect of inevitable errors in tracking. Different possible tuning strategies are discussed. Both the final signal-to-noise ratio and timing accuracy are evaluated and are shown to be significantly improved by the use of dynamic tuning. This should allow an accurate and reliable measurement of Hubble's constant.

The Besançon Galaxy Model renewed. I. Constraints on the Galactic thin disc evolution from Tycho data

Context. The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allow to test hypotheses on the star formation history, star evolution, as well as chemical and dynamical evolution of the Galaxy. Aims. The new version of the Besanc¸on Galaxy Model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the Initial Mass Function (IMF) and Star Formation Rate (SFR) of the Galactic disc. Methods. We present a new strategy for the generation of thin disc stars which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. We keep in this new scheme the local dynamical self-consistency as in Bienayme et al (1987). We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. Results. For the first time, a full sky comparison is performed between BGM and data. This strategy allows to constrain the IMF slope at high masses which is found to be close to 3.0, excluding a shallower slope such as Salpeter"s one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 M pc−3 implying that there is no compelling evidence for significant amount of dark matter in the disc. While the model is fitted to Tycho2 data, a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. Conclusions. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

The Besançon Galaxy Model renewed. I. Constraints on the Galactic thin disc evolution from Tycho data

Context. The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allow to test hypotheses on the star formation history, star evolution, as well as chemical and dynamical evolution of the Galaxy. Aims. The new version of the Besanc¸on Galaxy Model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the Initial Mass Function (IMF) and Star Formation Rate (SFR) of the Galactic disc. Methods. We present a new strategy for the generation of thin disc stars which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. We keep in this new scheme the local dynamical self-consistency as in Bienayme et al (1987). We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. Results. For the first time, a full sky comparison is performed between BGM and data. This strategy allows to constrain the IMF slope at high masses which is found to be close to 3.0, excluding a shallower slope such as Salpeter"s one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 M pc−3 implying that there is no compelling evidence for significant amount of dark matter in the disc. While the model is fitted to Tycho2 data, a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. Conclusions. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

Sobre a natureza dos processos físicos de aquecimento da atmosfera de estrelas do tipo solar

In the present work, we have studied the nature of the physical processes of the coronal heating, considering as basis significant samples of single and binary evolved stars, that have been achieved with the ROSAT satellite. In a total of 191 simple stars were studied, classified in the literature as giants with spectral type F, G and K. The results were compared with those obtained from 106 evolved stars of spectral type F, G and K, which belong to the spectroscopic binary systems. Accurate measurements on rotation and information about binarity were obtained from De Medeiros s catalog. We have analysed the behavior of the coronal activity in function of diverse stellar parameters. With the purpose to better clarify the profile of the stars evolution, the HR diagram was built for the two samples of stars, the single and the binary ones. The evolved traces added in the diagram were obtained from the Toulouse-Geneve code, Nascimento et al. (2000). The stars were segregated in this diagram not only in range of rotational speed but also in range of X-ray flux. Our analysis shows clearly that the single stars and the binary ones have coronal activity controlled by physical process independent on the rotation. Non magnetic processes seem to be strongly influencing the coronal heating. For the binary stars, we have also studied the behavior of the coronal emission as a function of orbital parameters, such as period and eccentricity, in which it was revealed the existence of a discontinuity in the emission of X-rays around an orbital period of 100 days. The study helped to conclude that circular orbits of the binary stars are presented as a necessary property for the existence of a higher level ofX-rays emission, suggesting that the effect of the gravitational tide has an important role in the coronal activity level. When applied the Kolmogorov-Smirnov test (KS test ) for the Vsini and FX parameters to the samples of single and binary stars, we could evidence very relevant aspects for the understanding of the mechanisms inherent to the coronal activity. For the Vsini parameter, the differences between the single stars and the binary ones for rotation over 6.3 km/s were really remarkable. We believe, therefore, that the existence of gravitational tide is, at least, one of the factors that most contribute for this behavior. About the X-rays flux, the KS test showed that the behavior of the single and the binary stars, regarding the coronal activity, comes from the same origin

Planet formation in a triple stellar system: implications of the third star's orbital inclination

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Researches on the evolution of the stellar systems ...

v. 1. On the universality of the law of gravitation and on the orbits and general characteristics of binary stars.--v. 2. The capture theory of cosmical evolution, founded on dynamical principles and illustrated by phenomena observed in the spiral nebulae, the planetary system, the double and multiple stars and clusters and the star-clouds of the Milky way.

Análise wavelet em curvas de luz estelares de sistemas binários da missão espacial CoRoT

Binary systems are key environments to study the fundamental properties of stars. In this work, we analyze 99 binary systems identified by the CoRoT space mission. From the study of the phase diagrams of these systems, our sample is divided into three groups: those whose systems are characterized by the variability relative to the binary eclipses; those presenting strong modulations probably due to the presence of stellar spots on the surface of star; and those whose systems have variability associated with the expansion and contraction of the surface layers. For eclipsing binary stars, phase diagrams are used to estimate the classification in regard to their morphology, based on the study of equipotential surfaces. In this context, to determine the rotation period, and to identify the presence of active regions, and to investigate if the star exhibits or not differential rotation and study stellar pulsation, we apply the wavelet procedure. The wavelet transform has been used as a powerful tool in the treatment of a large number of problems in astrophysics. Through the wavelet transform, one can perform an analysis in time-frequency light curves rich in details that contribute significantly to the study of phenomena associated with the rotation, the magnetic activity and stellar pulsations. In this work, we apply Morlet wavelet (6th order), which offers high time and frequency resolution and obtain local (energy distribution of the signal) and global (time integration of local map) wavelet power spectra. Using the wavelet analysis, we identify thirteen systems with periodicities related to the rotational modulation, besides the beating pattern signature in the local wavelet map of five pulsating stars over the entire time span.

Análise wavelet em curvas de luz estelares de sistemas binários da missão espacial CoRoT

Binary systems are key environments to study the fundamental properties of stars. In this work, we analyze 99 binary systems identified by the CoRoT space mission. From the study of the phase diagrams of these systems, our sample is divided into three groups: those whose systems are characterized by the variability relative to the binary eclipses; those presenting strong modulations probably due to the presence of stellar spots on the surface of star; and those whose systems have variability associated with the expansion and contraction of the surface layers. For eclipsing binary stars, phase diagrams are used to estimate the classification in regard to their morphology, based on the study of equipotential surfaces. In this context, to determine the rotation period, and to identify the presence of active regions, and to investigate if the star exhibits or not differential rotation and study stellar pulsation, we apply the wavelet procedure. The wavelet transform has been used as a powerful tool in the treatment of a large number of problems in astrophysics. Through the wavelet transform, one can perform an analysis in time-frequency light curves rich in details that contribute significantly to the study of phenomena associated with the rotation, the magnetic activity and stellar pulsations. In this work, we apply Morlet wavelet (6th order), which offers high time and frequency resolution and obtain local (energy distribution of the signal) and global (time integration of local map) wavelet power spectra. Using the wavelet analysis, we identify thirteen systems with periodicities related to the rotational modulation, besides the beating pattern signature in the local wavelet map of five pulsating stars over the entire time span.