Hough Transform Approach to Identification of Flare Stars in Multi-exposure Plate Images

ACM Computing Classification System (1998): I.2.8 , I.2.10, I.5.1, J.2.

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.

Searching for s-Process-Enhanced Metal-Poor Stars

We present preliminary results for the estimation of barium [Ba/Fe], and strontium [Sr/Fe], abundances ratios using medium-resolution spectra (1-2 angstrom). We established a calibration between the abundance ratios and line indices for Ba and Sr, using multiple regression and artificial neural network techniques. A comparison between the two techniques (showing the advantage of the latter), as well as a discussion of future work, is presented.

Searching for the first stars with the Gaia mission

Aims. We construct a theoretical model to predict the number of orphan afterglows (OA) from gamma-ray bursts (GRBs) triggered by primordial metal-free (Pop III) stars expected to be observed by the Gaia mission. In particular, we consider primordial metal-free stars that were affected by radiation from other stars (Pop III. 2) as a possible target. Methods. We use a semi-analytical approach that includes all relevant feedback effects to construct cosmic star formation history and its connection with the cumulative number of GRBs. The OA events are generated using the Monte Carlo method, and realistic simulations of Gaia's scanning law are performed to derive the observation probability expectation. Results. We show that Gaia can observe up to 2.28 +/- 0.88 off-axis afterglows and 2.78 +/- 1.41 on-axis during the five-year nominal mission. This implies that a nonnegligible percentage of afterglows that may be observed by Gaia (similar to 10%) could have Pop III stars as progenitors.

Candidate counterparts to the soft gamma-ray flare in the direction od LSI+61303

A short duration burst reminiscent of a soft gamma-ray repeater/anomalous X-ray pulsar behaviour was detected in the direction of LS I +61 303 by the Swift satellite. While the association with this well known gamma-ray binary is likely, a different origin cannot be excluded. Aims. We explore the error box of this unexpected flaring event and establish the radio, near-infrared and X-ray sources in our search for any peculiar alternative counterpart. Methods. We carried out a combined analysis of archive Very Large Array radio data of LS I +61 303 sensitive to both compact and extended emission. We also reanalysed previous near infrared observations with the 3.5 m telescope of the Centro Astronómico Hispano Alemán and X-ray observations with the Chandra satellite. Results. Our deep radio maps of the LS I +61 303 environment represent a significant advancement on previous work and 16 compact radio sources in the LS I +61 303 vicinity are detected. For some detections, we also identify near infrared and X-ray counterparts. Extended emission features in the field are also detected and confirmed. The possible connection of some of these sources with the observed flaring event is considered. Based on these data, we are unable to claim a clear association between the Swift-BAT flare and any of the sources reported here. However, this study represents the most sophisticated attempt to determine possible alternative counterparts other than LS I +61 303.

Searching for swept-up hydrogen and helium in the late-time spectra of 11 nearby Type Ia supernovae

The direct detection of a stellar system that explodes as a Type Ia supernova (SN Ia) has not yet been successful. Various indirect methods have been used to investigate SN Ia progenitor systems but none have produced conclusive results. A prediction of single-degenerate models is that H- (or He-) rich material from the envelope of the companion star should be swept up by the SN ejecta in the explosion. Seven SNe Ia have been analysed to date looking for signs of H-rich material in their late-time spectra and none were detected. We present results from new late-time spectra of 11 SNe Ia obtained at the Very Large Telescope using XShooter and FORS2. We present the tentative detection of Hα emission for SN 2013ct, corresponding to ∼0.007 M_⊙ of stripped/ablated companion star material (under the assumptions of the spectral modelling). This mass is significantly lower than expected for single-degenerate scenarios, suggesting that >0.1 M_⊙ of H-rich is present but not observed. We do not detect Hα emission in the other 10 SNe Ia. This brings the total sample of normal SNe Ia with non-detections (<0.001–0.058 M_⊙) of H-rich material to 17 events. The simplest explanation for these non-detections is that these objects did not result from the explosion of a CO white dwarf accreting matter from a H-rich companion star via Roche lobe overflow or symbiotic channels. However, further spectral modelling is needed to confirm this. We also find no evidence of He-emission features, but models with He-rich companion stars are not available to place mass limits.

Abbondanze chimiche delle Blue Straggler Stars dell'ammasso aperto NGC 188

My Thesis work is aimed at searching for evidence of CO-depletion in the atmosphere of these objects. To this end, I am analyzing high-resolution spectra recently acquired with the spectrograph HDS at the Subaru Telescope. The global sample is composed of 5 BSSs with hot WD companions, 3 additional binary BSSs with no a significant far-UV excess, and one "normal" stars (along the RGB star) needed for comparison. The analysis will provide us with the surface chemical abundances (especially of C and O) of the target stars. The final goal is to verify the expected {chemical signature} and put constraints on its characteristic time scale (for instance, a positive detection in the 5 BSSs with hot companions and a non-detection in the other 3 binary BSSs would imply that CO-depletion is a transient phenomenon, lasting approx 300 Myr only). The analysis will also provide us with the rotational velocities of each target, thus allowing to investigate their kinematical properties and shed new light on their evolutionary path.

A SEARCH FOR FIELD HORIZONTAL-BRANCH STARS NEAR THE SOUTH GALACTIC POLE

Searches for field horizontal-branch (FHB) stars in the halo of the Galaxy in the past have been carried out by several techniques, such as objective-prism surveys and visual or infrared photometric surveys. By choosing adequate color criteria, it is possible to improve the efficiency of identifying bona fide FHB stars among the other objects that exhibit similar characteristics, such as main-sequence A-stars, blue stragglers, subdwarfs, etc. In this work, we report the results of a spectroscopic survey carried out near the south Galactic pole intended to validate FHB stars originally selected from the HK objective-prism survey of Beers and colleagues, based on near-infrared color indices. A comparison between the stellar spectra obtained in this survey with theoretical stellar atmosphere models allows us to determine T(eff), log g, and [Fe/H] for 13 stars in the sample. Stellar temperatures were calculated from measured (B-V)(o), when this measurement was available (16 stars). The color index criteria adopted in this work are shown to correctly classify 30% of the sample as FHB, 25% as non-FHB (main-sequence stars and subdwarfes), whereas 40% could not be distinguished between FHB and main-sequence stars. We compare the efficacy of different color criteria in the literature intended to select FHB stars, and discuss the use of the Mg II 4481 line to estimate the metallicity.

Searching for molecular markers in head and neck squamous cell carcinomas (HNSCC) by statistical and bioinformatic analysis of larynx-derived SAGE libraries

Background: Head and neck squamous cell carcinoma (HNSCC) is one of the most common malignancies in humans. The average 5-year survival rate is one of the lowest among aggressive cancers, showing no significant improvement in recent years. When detected early, HNSCC has a good prognosis, but most patients present metastatic disease at the time of diagnosis, which significantly reduces survival rate. Despite extensive research, no molecular markers are currently available for diagnostic or prognostic purposes. Methods: Aiming to identify differentially-expressed genes involved in laryngeal squamous cell carcinoma (LSCC) development and progression, we generated individual Serial Analysis of Gene Expression (SAGE) libraries from a metastatic and non-metastatic larynx carcinoma, as well as from a normal larynx mucosa sample. Approximately 54,000 unique tags were sequenced in three libraries. Results: Statistical data analysis identified a subset of 1,216 differentially expressed tags between tumor and normal libraries, and 894 differentially expressed tags between metastatic and non-metastatic carcinomas. Three genes displaying differential regulation, one down-regulated (KRT31) and two up-regulated (BST2, MFAP2), as well as one with a non-significant differential expression pattern (GNA15) in our SAGE data were selected for real-time polymerase chain reaction (PCR) in a set of HNSCC samples. Consistent with our statistical analysis, quantitative PCR confirmed the upregulation of BST2 and MFAP2 and the downregulation of KRT31 when samples of HNSCC were compared to tumor-free surgical margins. As expected, GNA15 presented a non-significant differential expression pattern when tumor samples were compared to normal tissues. Conclusion: To the best of our knowledge, this is the first study reporting SAGE data in head and neck squamous cell tumors. Statistical analysis was effective in identifying differentially expressed genes reportedly involved in cancer development. The differential expression of a subset of genes was confirmed in additional larynx carcinoma samples and in carcinomas from a distinct head and neck subsite. This result suggests the existence of potential common biomarkers for prognosis and targeted-therapy development in this heterogeneous type of tumor.

ESPECTRA: Searching the Bipolar Spectrum in Eating Disorder patients

Background: Bipolar Disorder (BD) is a chronic, recurrent and highly prevalent illness. Despite the need for correct diagnosis to allow proper treatment, studies have shown that reaching a diagnosis can take up to ten years due to the lack of recognition of the broader presentations of BD. Frequent comorbidities with other psychiatric disorders are a major cause of misdiagnosis and warrant thorough evaluation. Methods/Design: ESPECTRA (Occurrence of Bipolar Spectrum Disorders in Eating Disorder Patients) is a single-site cross-sectional study involving a comparison group, designed to evaluate the prevalence of bipolar spectrum in an eating disorder sample. Women aged 18-45 years will be evaluated using the SCID-P and Zurich criteria for diagnosis and the HAM-D, YOUNG, SCI-MOODS, HCL-32, BIS-11, BSQ, WHOQoL and EAS instruments for rating symptoms and measuring clinical correlates. Discussion: The classificatory systems in psychiatry are based on categorical models that have been criticized for simplifying the diagnosis and leading to an increase in comorbidities. Some dimensional approaches have been proposed aimed at improving the validity and reliability of psychiatric disorder assessments, especially in conditions with high rates of comorbidity such as BD and Eating Disorder (ED). The Bipolar Spectrum (BS) remains under-recognized in clinical practice and its definition is not well established in current diagnostic guidelines. Broader evaluation of psychiatric disorders combining categorical and dimensional views could contribute to a more realistic understanding of comorbidities and help toward establishing a prognosis.

On the mass determination of super-Earths orbiting active stars: the CoRoT-7 system

Context. Star activity makes the mass determination of CoRoT-7b and CoRoT 7c uncertain. Investigators of the CoRoT team proposed several solutions, but all but one of them are larger than the initial determinations of 4.8 +/- 0.8 M(Earth) for CoRoT-7b and 8.4 +/- 0.9 M(Earth) for CoRoT 7c. Aims. This investigation uses the excellent HARPS radial velocity measurements of CoRoT-7 to redetermine the planet masses and to explore techniques for determining mass and orbital elements of planets discovered around active stars when the relative variation in the radial velocity due to the star activity cannot be considered as just noise and can exceed the variation due to the planets. Methods. The main technique used here is a self-consistent version of the high-pass filter used by Queloz et al. (2009, A&A, 506, 303) in the first mass determination of CoRoT-7b and CoRoT-7c. The results are compared to those given by two alternative techniques: (1) the approach proposed by Hatzes et al. (2010, A&A, 520, A93) using only those nights in which two or three observations were done; (2) a pure Fourier analysis. In all cases, the eccentricities are taken equal to zero as indicated by the study of the tidal evolution of the system. The periods are also kept fixed at the values given by Queloz et al. Only the observations done in the time interval BJD 2 454 847-873 are used because they include many nights with multiple observations; otherwise, it is not possible to separate the effects of the rotation fourth harmonic (5.91 d = P(rot)/4) from the alias of the orbital period of CoRoT-7b (0.853585 d). Results. The results of the various approaches are combined to give planet mass values 8.0 +/- 1.2 M(Earth) for CoRoT-7b and 13.6 +/- 1.4 M(Earth) for CoRoT 7c. An estimation of the variation of the radial velocity of the star due to its activity is also given. Conclusions. The results obtained with three different approaches agree to give higher masses than those in previous determinations. From the existing internal structure models they indicate that CoRoT-7b is a much denser super-Earth. The bulk density is 11 +/- 3.5 g cm(-3), so CoRoT-7b may be rocky with a large iron core.

Abundances in the Galactic bulge: results from planetary nebulae and giant stars

Context. Our understanding of the chemical evolution (CE) of the Galactic bulge requires the determination of abundances in large samples of giant stars and planetary nebulae (PNe). Studies based on high resolution spectroscopy of giant stars in several fields of the Galactic bulge obtained with very large telescopes have allowed important progress. Aims. We discuss PNe abundances in the Galactic bulge and compare these results with those presented in the literature for giant stars. Methods. We present the largest, high-quality data-set available for PNe in the direction of the Galactic bulge (inner-disk/bulge). For comparison purposes, we also consider a sample of PNe in the Large Magellanic Cloud (LMC). We derive the element abundances in a consistent way for all the PNe studied. By comparing the abundances for the bulge, inner-disk, and LMC, we identify elements that have not been modified during the evolution of the PN progenitor and can be used to trace the bulge chemical enrichment history. We then compare the PN abundances with abundances of bulge field giant. Results. At the metallicity of the bulge, we find that the abundances of O and Ne are close to the values for the interstellar medium at the time of the PN progenitor formation, and hence these elements can be used as tracers of the bulge CE, in the same way as S and Ar, which are not expected to be affected by nucleosynthetic processes during the evolution of the PN progenitors. The PN oxygen abundance distribution is shifted to lower values by 0.3 dex with respect to the distribution given by giants. A similar shift appears to occur for Ne and S. We discuss possible reasons for this PNe-giant discrepancy and conclude that this is probably due to systematic errors in the abundance derivations in either giants or PNe (or both). We issue an important warning concerning the use of absolute abundances in CE studies.

The metal content of bulge field stars from FLAMES-GIRAFFE spectra - I. Stellar parameters and iron abundances

Aims. We determine the iron distribution function (IDF) for bulge field stars, in three different fields along the Galactic minor axis and at latitudes b = -4 degrees, b = -6 degrees, and b = -12 degrees. A fourth field including NGC 6553 is also included in the discussion. Methods. About 800 bulge field K giants were observed with the GIRAFFE spectrograph of FLAMES@VLT at spectral resolution R similar to 20 000. Several of them were observed again with UVES at R similar to 45 000 to insure the accuracy of the measurements. The LTE abundance analysis yielded stellar parameters and iron abundances that allowed us to construct an IDF for the bulge that, for the first time, is based on high-resolution spectroscopy for each individual star. Results. The IDF derived here is centered on solar metallicity, and extends from [Fe/H] similar to -1.5 to [Fe/H] similar to + 0.5. The distribution is asymmetric, with a sharper cutoff on the high-metallicity side, and it is narrower than previously measured. A variation in the mean metallicity along the bulge minor axis is clearly between b = -4 degrees and b = -6 degrees ([Fe/H] decreasing similar to by 0.6 dex per kpc). The field at b = -12 degrees. is consistent with the presence of a gradient, but its quantification is complicated by the higher disk/bulge fraction in this field. Conclusions. Our findings support a scenario in which both infall and outflow were important during the bulge formation, and then suggest the presence of a radial gradient, which poses some challenges to the scenario in which the bulge would result solely from the vertical heating of the bar.

Chemical similarities between Galactic bulge and local thick disk red giant stars

Context. The evolution of the Milky Way bulge and its relationship with the other Galactic populations is still poorly understood. The bulge has been suggested to be either a merger-driven classical bulge or the product of a dynamical instability of the inner disk. Aims. To probe the star formation history, the initial mass function and stellar nucleosynthesis of the bulge, we performed an elemental abundance analysis of bulge red giant stars. We also completed an identical study of local thin disk, thick disk and halo giants to establish the chemical differences and similarities between the various populations. Methods. High-resolution infrared spectra of 19 bulge giants and 49 comparison giants in the solar neighborhood were acquired with Gemini/Phoenix. All stars have similar stellar parameters but cover a broad range in metallicity. A standard 1D local thermodynamic equilibrium analysis yielded the abundances of C, N, O and Fe. A homogeneous and differential analysis of the bulge, halo, thin disk and thick disk stars ensured that systematic errors were minimized. Results. We confirm the well-established differences for [O/Fe] (at a given metallicity) between the local thin and thick disks. For the elements investigated, we find no chemical distinction between the bulge and the local thick disk, which is in contrast to previous studies relying on literature values for disk dwarf stars in the solar neighborhood. Conclusions. Our findings suggest that the bulge and local thick disk experienced similar, but not necessarily shared, chemical evolution histories. We argue that their formation timescales, star formation rates and initial mass functions were similar.

High spatial resolution spectroscopy of W51IRS 2E and IRS 2W: Two very massive young stars in early formation stages

We present K-band spectra of the near infrared counterparts to IRS 2E and IRS 2W which is associated with the ultracompact H II region W51d, both of them embedded sources in the Galactic compact H II region W51 IRS 2. The high spatial resolution observations were obtained with the laser guide star facility and Near-infrared Integral Field Spectrograph (NIFS) mounted at the Gemini-North observatory. The spectrum of the ionizing source of W51d shows the photospheric features N III ( 21155 angstrom) in emission and He II ( 21897 angstrom) in absorption which lead us to classify it as a young O3 type star. We detected CO overtone in emission at 23000 angstrom in the spectrum of IRS 2E, suggesting that it is a massive young object still surrounded by an accretion disk, probably transitioning from the hot core phase to an ultracompact H II region.