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.

The periodic variations of a white-light flare observed with ULTRACAM

High time resolution observations of a white-light flare on the active star EQ PegB show evidence of intensity variations with a period of ≈10 s. The period drifts to longer values during the decay phase of the flare. If the oscillation is interpreted as an impulsively-excited, standing-acoustic wave in a flare loop, the period implies a loop length of ≈3.4 Mm and ≈6.8 Mm for the case of the fundamental mode and the second harmonic, respectively. However, the small loop lengths imply a very high modulation depth making the acoustic interpretation unlikely. A more realistic interpretation may be that of a fast-MHD wave, with the modulation of the emission being due to the magnetic field. Alternatively, the variations could be due to a series of reconnection events. The periodic signature may then arise as a result of the lateral separation of individual flare loops or current sheets with oscillatory dynamics (i.e., periodic reconnection).

M Dwarf Flare Continuum Variations on One-second Timescales: Calibrating and Modeling of ULTRACAM Flare Color Indices

We present a large data set of high-cadence dMe flare light curves obtained with custom continuum filters on the triple-beam, high-speed camera system ULTRACAM. The measurements provide constraints for models of the near-ultraviolet (NUV) and optical continuum spectral evolution on timescales of ≈1 s. We provide a robust interpretation of the flare emission in the ULTRACAM filters using simultaneously obtained low-resolution spectra during two moderate-sized flares in the dM4.5e star YZ CMi. By avoiding the spectral complexity within the broadband Johnson filters, the ULTRACAM filters are shown to characterize bona fide continuum emission in the NUV, blue, and red wavelength regimes. The NUV/blue flux ratio in flares is equivalent to a Balmer jump ratio, and the blue/red flux ratio provides an estimate for the color temperature of the optical continuum emission. We present a new “color-color” relationship for these continuum flux ratios at the peaks of the flares. Using the RADYN and RH codes, we interpret the ULTRACAM filter emission using the dominant emission processes from a radiative-hydrodynamic flare model with a high nonthermal electron beam flux, which explains a hot, T ≈ 104 K, color temperature at blue-to-red optical wavelengths and a small Balmer jump ratio as observed in moderate-sized and large flares alike. We also discuss the high time resolution, high signal-to-noise continuum color variations observed in YZ CMi during a giant flare, which increased the NUV flux from this star by over a factor of 100. Based on observations obtained with the Apache Point Observatory 3.5 m telescope, which is owned and operated by the Astrophysical Research Consortium, based on observations made with the William Herschel Telescope operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofsica de Canarias, and observations, and based on observations made with the ESO Telescopes at the La Silla Paranal Observatory under programme ID 085.D-0501(A).

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.

The Decaying Long-period Oscillation of a Stellar Megaflare

We analyze and interpret the oscillatory signal in the decay phase of the U-band light curve of a stellar megaflare observed on 2009 January 16 on the dM4.5e star YZ CMi. The oscillation is well approximated by an exponentially decaying harmonic function. The period of the oscillation is found to be 32 minutes, the decay time about 46 minutes, and the relative amplitude 15%. As this observational signature is typical of the longitudinal oscillations observed in solar flares at extreme ultraviolet and radio wavelengths, associated with standing slow magnetoacoustic waves, we suggest that this megaflare may be of a similar nature. In this scenario, macroscopic variations of the plasma parameters in the oscillations modulate the ejection of non-thermal electrons. The phase speed of the longitudinal (slow magnetoacoustic) waves in the flaring loop or arcade, the tube speed, of about 230 km s-1 would require a loop length of about 200 Mm. Other mechanisms, such as standing kink oscillations, are also considered.

Searching the Sky with CONFIGR-STARS

SyNAPSE program of the Defense Advanced Projects Research Agency (HRL Laboratories LLC, subcontract #801881-BS under DARPA prime contract HR0011-09-C-0001); CELEST, a National Science Foundation Science of Learning Center (SBE-0354378)

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.

A search for line intensity enhancements in the far-UV spectra of active late-type stars arising from opacity

Context. Radiative transfer calculations have predicted intensity enhancements for optically thick emission lines, as opposed to the normal intensity reductions, for astrophysical plasmas under certain conditions. In particular, the results are predicted to be dependent both on the geometry of the emitting plasma and the orientation of the observer. Hence in principle the detection of intensity enhancement may provide a way of determining the geometry of an unresolved astronomical source.
Aims. To investigate such enhancements we have analysed a sample of active late-type stars observed in the far ultraviolet spectral region.
Methods. Emission lines of O vi in the FUSE satellite spectra of ϵ Eri, II Peg and Prox Cen were searched for intensity enhancements due to opacity.
Results. We have found strong evidence for line intensity enhancements due to opacity during active or flare-like activity for all three stars. The O vi 1032/1038 line intensity ratios, predicted to have a value of 2.0 in the optically thin case, are found to be up to ~30% larger during several orbital phases.
Conclusions. Our measurements, combined with radiative transfer models, allow us to constrain both the geometry of the O vi emitting regions in our stellar sources and the orientation of the observer. A spherical emitting plasma can be ruled out, as this would lead to no intensity enhancement. In addition, the theory tells us that the line-of-sight to the plasma must be close to perpendicular to its surface, as observations at small angles to the surface lead to either no intensity enhancement or the usual line intensity decrease over the optically thin value. For the future, we outline a laboratory experiment, that could be undertaken with current facilities, which would provide an unequivocal test of predictions of line intensity enhancement due to opacity, in particular the dependence on plasma geometry.

New identifications for blue objects towards the Galactic center: post-AGE stars, Be/disk stars and others

As part of a programme to investigate spatial variations in the Galactic chemical composition, we have been searching for normal B-type stars and A-type supergiants near the Galactic center. During this search we have found eleven peculiar stars, and in some cases performed detailed abundance analyses of them which suggest that they may be at a post-AGB evolutionary stage.

Finding exoplanets orbiting young active stars - I. Technique

Stellar activity, such as starspots, can induce radial velocity (RV) variations that can mask or even mimic the RV signature of orbiting exoplanets. For this reason RV exoplanet surveys have been unsuccessful when searching for planets around young, active stars and are therefore failing to explore an important regime which can help to reveal how planets form and migrate. This paper describes a new technique to remove spot signatures from the stellar line-profiles of moderately rotating, active stars (v sin i ranging from 10 to 50 km s(-1)). By doing so it allows planetary RV signals to be uncovered. We used simulated models of a G5V type star with differing dark spots on its surface along with archive data of the known active star HD 49933 to validate our method. The results showed that starspots could be effectively cleaned from the line-profiles so that the stellar RV jitter was reduced by more than 80 per cent. Applying this procedure to the same models and HD 49933 data, but with fake planets injected, enabled the effective removal of starspots so that Jupiter mass planets on short orbital periods were successfully recovered. These results show that this approach can be useful in the search for hot-Jupiter planets that orbit around young, active stars with a v sin i of similar to 10-50 km/s.

Evolution of egg size and fertilisation efficiency in sea stars: large eggs are not fertilised more readily than small eggs in Patiriella (Echinodermata: Asteroidea)

Fertilisation kinetics theory suggests that, when sperm are limiting, the larger eggs of broadcast-spawning marine organisms ought to be fertilised more frequently than smaller eggs, because they provide a bigger target for searching sperm. Whilst this effect has been demonstrated within species, it is not known if this pattern holds among species. We tested whether a large difference in egg size between congeneric seastars with contrasting planktotrophic and lecithotrophic modes of development results in differences in the likelihood of eggs being fertilised in sperm-limiting situations. Measurement of egg sizes and sperm swimming speeds led to the prediction that the sperm–egg collision rate constant for Patiriella calcar (420-µm-diameter egg) should be nine times greater than for P. regularis (140-µm-diameter egg). Although the eggs of P. calcar should be fertilised at greater rates in low sperm concentrations, they were not. When gametes were allowed to mix for 10 s, the hypothesis that P. calcar eggs required less sperm than P. regularis to ensure 50% of eggs were fertilised was rejected. When gametes were mixed for 5 min, P. regularis eggs were more frequently fertilised, but the difference was not statistically significant. We conclude there must be a difference between these species in the likelihood that when a sperm finds a conspecific egg it can successfully fertilise. This apparent uncoupling of egg size and likelihood of fertilisation suggests that fertilisation is not a major constraint on the evolution of egg size in these seastars.

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.