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.

Some Topics in Relation to José Celestino Mutis and Medicine

José Celestino Mutis showed different interest areas of study and worked in multiple knowledge fields: mathematics, mineralogy, astronomy, zoology and botanic, the one in which is best known. However, he was a medical doctor, whit profound an up to date knowledge. He studied at Sevilla University, and in the Royal College of Cadiz. He arrived to the “Nuevo Reino de Granada” as personal medical doctor of the the viceroy Pedro Messía de la Cerda and during 47 years, period of time he lived there, he always exercised his profession. We showed a special interest in the study of cinchona and the only work he published was El Arcano de la quina. He advised in public health many viceroys and he contributed in different fields but particularly in the prevention of smallpox by inoculation. He played a decisive role in the reopening of the medical school of the Colegio Mayor de Nuestra Señora del Rosario. Some of his achievements justifies the appellative “Father of the Medicine in Colombia”.

Families of periodic horseshoe orbits in the restricted three-body problem

We compute families of symmetric periodic horseshoe orbits in the restricted three-body problem. Both the planar and three-dimensional cases are considered and several families are found.We describe how these families are organized as well as the behavior along and among the families of parameters such as the Jacobi constant or the eccentricity. We also determine the stability properties of individual orbits along the families. Interestingly, we find stable horseshoe-shaped orbit up to the quite high inclination of 17◦

Cooled filters 16-42µm for FIR astronomy

The design and manufacture of dielectric-film interference filters for cooled FIR stronmy is described. The bands are 16.5-21.5µm, 17.5-19.5µm, 19.5-21.5µm and 27µm cut on. The films are PbTe/CdSe and the substrates are CdTe (some 1/2 mm thick), without absorption in the region: KRS-6 films are used for antireflection.

The Maunder minimum (1645-1715) was indeed a grand minimum: a reassessment of multiple datasets

Aims. Although the time of the Maunder minimum (1645–1715) is widely known as a period of extremely low solar activity, it is still being debated whether solar activity during that period might have been moderate or even higher than the current solar cycle (number 24). We have revisited all existing evidence and datasets, both direct and indirect, to assess the level of solar activity during the Maunder minimum. Methods. We discuss the East Asian naked-eye sunspot observations, the telescopic solar observations, the fraction of sunspot active days, the latitudinal extent of sunspot positions, auroral sightings at high latitudes, cosmogenic radionuclide data as well as solar eclipse observations for that period. We also consider peculiar features of the Sun (very strong hemispheric asymmetry of the sunspot location, unusual differential rotation and the lack of the K-corona) that imply a special mode of solar activity during the Maunder minimum. Results. The level of solar activity during the Maunder minimum is reassessed on the basis of all available datasets. Conclusions. We conclude that solar activity was indeed at an exceptionally low level during the Maunder minimum. Although the exact level is still unclear, it was definitely lower than during the Dalton minimum of around 1800 and significantly below that of the current solar cycle #24. Claims of a moderate-to-high level of solar activity during the Maunder minimum are rejected with a high confidence level.

A METHOD FOR THE STUDY OF ACCRETION DISK EMISSION IN CATACLYSMIC VARIABLES. I. THE MODEL

We have developed a spectrum synthesis method for modeling the ultraviolet (UV) emission from the accretion disk from cataclysmic variables (CVs). The disk is separated into concentric rings, with an internal structure from the Wade & Hubeny disk-atmosphere models. For each ring, a wind atmosphere is calculated in the comoving frame with a vertical velocity structure obtained from a solution of the Euler equation. Using simple assumptions, regarding rotation and the wind streamlines, these one-dimensional models are combined into a single 2.5-dimensional model for which we compute synthetic spectra. We find that the resulting line and continuum behavior as a function of the orbital inclination is consistent with the observations, and verify that the accretion rate affects the wind temperature, leading to corresponding trends in the intensity of UV lines. In general, we also find that the primary mass has a strong effect on the P Cygni absorption profiles, the synthetic emission line profiles are strongly sensitive to the wind temperature structure, and an increase in the mass-loss rate enhances the resonance line intensities. Synthetic spectra were compared with UV data for two high orbital inclination nova-like CVs-RW Tri and V347 Pup. We needed to include disk regions with arbitrary enhanced mass loss to reproduce reasonably well widths and line profiles. This fact and a lack of flux in some high ionization lines may be the signature of the presence of density-enhanced regions in the wind, or alternatively, may result from inadequacies in some of our simplifying assumptions.

Can retired galaxies mimic active galaxies? Clues from the Sloan Digital Sky Survey

The classification of galaxies as star forming or active is generally done in the ([O III]/H beta, [N II]/H alpha) plane. The Sloan Digital Sky Survey (SDSS) has revealed that, in this plane, the distribution of galaxies looks like the two wings of a seagull. Galaxies in the right wing are referred to as Seyfert/LINERs, leading to the idea that non-stellar activity in galaxies is a very common phenomenon. Here, we argue that a large fraction of the systems in the right wing could actually be galaxies which stopped forming stars. The ionization in these `retired` galaxies would be produced by hot post-asymptotic giant branch stars and white dwarfs. Our argumentation is based on a stellar population analysis of the galaxies via our STARLIGHT code and on photoionization models using the Lyman continuum radiation predicted for this population. The proportion of LINER galaxies that can be explained in such a way is, however, uncertain. We further show how observational selection effects account for the shape of the right wing. Our study suggests that nuclear activity may not be as common as thought. If retired galaxies do explain a large part of the seagull`s right wing, some of the work concerning nuclear activity in galaxies, as inferred from SDSS data, will have to be revised.

Discovery in IC10 of the farthest known symbiotic star

We report the discovery of the first known symbiotic star in IC10, a starburst galaxy belonging to the Local Group, at a distance of similar to 750 kpc. The symbiotic star was identified during a survey of emission-line objects. It shines at V = 24.62 +/- 0.04, V - R(C) = 2.77 +/- 0.05 and R(C) - I(C) = 2.39 +/- 0.02, and suffers from E(B-V) = 0.85 +/- 0.05 reddening. The spectrum of the cool component well matches that of solar neighbourhood M8III giants. The observed emission lines belong to Balmer series, [S II], [N II] and [O III]. They suggest a low electronic density, negligible optical depth effects and 35 000 < T(eff) < 90 000 K for the ionizing source. The spectrum of the new symbiotic star in IC10 is an almost perfect copy of that of Hen 2-147, a well-known Galactic symbiotic star and Mira.

Contribution to the ground-based follow-up of the Gaia space mission

The Gaia Space Mission [Mignard, F., 2005. The three-dimensional universe with Gaia. ESA/SP-576; Perryman, M., 2005. The three-dimensional universe with Gaia. ESA/SP-576] will observe several transient events as supernovae, microlensing, gamma ray bursts and new Solar System objects. The satellite, due to its scanning law, will detect these events but will not be able to monitor them. So, to take these events into consideration and to perform further studies it is necessary to follow them with Earth-based observations. These observations could be efficiently done by a ground-based network of well-equipped telescopes scattered in both hemispheres. Here we focus our attention at the new Solar System objects to be discovered and observed by the Gaia satellite [Mignard, F., 2002. Observations of Solar System objects by Gaia I. Detection of NEOS. Astron. Astrophys. 393, 727] mainly asteroids, NEOs and comets. A dedicated ground-based network of telescopes as proposed by Thuillot [2005. The three-dimensional universe with Gaia. ESA/SP-576] will allow to monitor those events, to avoid losing them and to perform a quick characterization of some physical properties which will be important for the identification of these objects in further measurements by Gaia. We present in this paper, the beginning of the organization of a Latin-American ground-based network of telescopes and observers joining several institutions in Argentina, Bolivia, Brazil and other Latin-American countries aiming to contribute to the follow-up of Gaia science alerts for Solar System objects. (C) 2008 Elsevier Ltd. All rights reserved.

Transition redshift: new kinematic constraints from supernovae

The transition redshift (deceleration/acceleration) is discussed by expanding the deceleration parameter to first order around its present value. A detailed study is carried out by considering two different parametrizations, q = q(0) + q(1)z and q = q(0) + q(1)z(1 + z)(-1), and the associated free parameters (q(0), q(1)) are constrained by three different supernovae (SNe) samples. A previous analysis by Riess et al. using the first expansion is slightly improved and confirmed in light of their recent data (Gold07 sample). However, by fitting the model with the Supernova Legacy Survey (SNLS) type Ia sample, we find that the best fit to the redshift transition is z(t) = 0.61, instead of z(t) = 0.46 as derived by the High-z Supernovae Search (HZSNS) team. This result based in the SNLS sample is also in good agreement with the sample of Davis et al., z(t) = 0.60(-0.11)(+0.28) (1 sigma). Such results are in line with some independent analyses and accommodate more easily the concordance flat model (Lambda CDM). For both parametrizations, the three SNe Ia samples considered favour recent acceleration and past deceleration with a high degree of statistical confidence level. All the kinematic results presented here depend neither on the validity of general relativity nor on the matter-energy contents of the Universe.

The embedded cluster DBSB 48 in the nebula Hoffleit 18: Comparison with Trumpler 14

We derive fundamental parameters of the embedded cluster DBSB 48 in the southern nebula Hoffleit 18 and the very young open cluster Trumpler 14, by means of deep JHK(s) infrared photometry. We build colour-magnitude and colour-colour diagrams to derive reddening and age, based on main sequence and pre-main sequence distributions. Radial stellar density profiles are used to study cluster structure and guide photometric diagram extractions. Field-star decontamination is applied to uncover the intrinsic cluster sequences in the diagrams. Ages are inferred from K-excess fractions. A prominent pre-main sequence population is present in DBSB 48, and the K-excess fraction f(K) = 55 +/- 6% gives an age of 1.1 +/- 0.5 Myr. A mean reddening of A(Ks) = 0.9 +/- 0.03 was found, corresponding to A(v) = 8.2 +/- 0.3. The cluster CMD is consistent with the far kinematic distance of 5 kpc for Hoffleit 18. For Trumpler 14 we derived similar parameters as in previous studies in the optical, in particular an age of 1.7 +/- 0.7 Myr. The fraction of stars with infrared excess in Trumpler 14 is f(K) = 28 +/- 4%. Despite the young ages, both clusters are described by a King profile with core radii R-core = 0.46 +/- 0.05 pc and R-core = 0.35 +/- 0.04 pc, respectively, for DBSB 48 and Trumpler 14. Such cores are smaller than those of typical open clusters. Small cores are probably related to the cluster formation and/or parent molecular cloud fragmentation. In DBSB 48, the magnitude extent of the upper main sequence is Delta K-s approximate to 2 mag, while in Trumpler 14 it is Delta K-s approximate to 5 mag, consistent with the estimated ages. (c) 2008 Elsevier B.V. All rights reserved.

Non-LTE Monte Carlo radiative transfer. II. Nonisothermal solutions for viscous Keplerian disks

We discuss the basic hydrodynamics that determines the density structure of the disks around hot stars. Observational evidence supports the idea that these disks are Keplerian (rotationally supported) gaseous disks. A popular scenario in the literature, which naturally leads to the formation of Keplerian disks, is the viscous decretion model. According to this scenario, the disks are hydrostatically supported in the vertical direction, while the radial structure is governed by the viscous transport. This suggests that the temperature is one primary factor that governs the disk density structure. In a previous study we demonstrated, using three-dimensional non-LTE Monte Carlo simulations, that viscous Keplerian disks can be highly nonisothermal. In this paper we build on our previous work and solve the full problem of the steady state nonisothermal viscous diffusion and vertical hydrostatic equilibrium. We find that the self-consistent solution departs significantly from the analytic isothermal density, with potentially large effects on the emergent spectrum. This implies that nonisothermal disk models must be used for a detailed modeling of Be star disks.

Hydrodynamical simulations of Galactic fountains - I. Evolution of single fountains

The ejection of the gas out of the disc in late-type galaxies is related to star formation and is due mainly to Type II supernovae. In this paper, we studied in detail the development of the Galactic fountains in order to understand their dynamical evolution and their influence on the redistribution of the freshly delivered metals over the disc. To this aim, we performed a number of 3D hydrodynamical radiative cooling simulations of the gas in the Milky Way where the whole Galaxy structure, the Galactic differential rotation and the supernova explosions generated by a single OB association are considered. A typical fountain powered by 100 Type II supernovae may eject material up to similar to 2 kpc which than collapses back mostly in the form of dense, cold clouds and filaments. The majority of the gas lifted up by the fountains falls back on the disc remaining within a radial distance Delta R = 0.5 kpc from the place where the fountain originated. This localized circulation of disc gas does not influence the radial chemical gradients on large scale, as required by the chemical models of the Milky Way which reproduce the metallicity distribution without invoking large fluxes of metals. Simulations of multiple fountains fuelled by Type II supernovae of different OB associations will be presented in a companion paper.