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

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

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

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

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

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

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

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

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

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

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

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

X-shooting Herbig Ae/Be stars: Accretion probed by near-infrared He I emission

The Herbig Ae/Be stars are intermediate mass pre-main sequence stars that bridge the gap between the low mass T Tauri stars and the Massive Young Stellar Objects. In this mass range, the acting star forming mechanism switches from magnetically controlled accretion to an as yet unknown mechanism, but which is likely to be direct disk accretion onto the star. We observed a large sample of Herbig Ae/Be stars with X-shooter to address this issue from a multi-wavelength perspective. It is the largest such study to date, not only because of the number of objects involved, but also because of the large wavelength coverage from the blue to the near-infrared. This allows many accretion diagnostics to be studied simultaneously. By correlating the various properties with mass, temperature and age, we aim to determine where and whether the magnetically controlled mass accretion mechanism halts and the proposed direct disk accretion takes over. Here, we will give an overview of the background, present some observations and discuss our initial results. We will introduce a new accretion diagnostic for the research of Herbig Ae/Be stars, the HeI 1.083 µm line. © 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

CSI 2264: Simultaneous Optical and Infrared Light Curves of Young Disk-bearing Stars in NGC 2264 with CoRoT and Spitzer—Evidence for Multiple Origins of Variability

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30 day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" with discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with Hα emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk. Based on data from the Spitzer and CoRoT missions. The CoRoT space mission was developed and is operated by the French space agency CNES, with participation of ESA's RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, and Spain.

Hα emission fluxes and lithium abundances of low-mass stars in the young open cluster IC 4665

As part of a long term effort to understand pre-main sequence Li burning, we have obtained high resolution spectroscopic observations of 14 late type stars (G0-M1) in the young open cluster IC 4665. Most of the stars have Hα filled-in and Li I absorption, as expected for their young age. From the equivalent widths of Hα emission excess (obtained using the spectral subtraction technique) and the the Li i λ6708 feature, we have derived Hα emission fluxes and photospheric Li abundances. The mean Li abundance of IC 4665 solar-type stars is log N(Li) = 3.1; the same as in other young clusters (α Per, Pleiades) and T Tauri stars. Our results support the conclusions from previous works that PMS Li depletion is very small for masses ∼ 1 M_⨀ . Among the IC 4665 late-G and early K-type stars, there is a spread in Li abundances of about one order of magnitude. The Li-poor IC 4665 members have low Hα excess and vsini≤10. Hence, the Li-activity-rotation connection which has been clearly established in the Pleiades also seems to hold in IC 4665. One M-type IC 4665 star that we have observed does not show Li, implying a very efficient Li depletion as observed in α Per stars of the same spectral type. The level of chromospheric activity and Li depletion among the low-mass stars of IC 4665 is similar to that in the Pleiades. In fact, we note that the Li abundance distributions in several young clusters (α Per, Pleiades, IC 2391, IC 4665) and in post T Tauri stars are strikingly similar. This result suggests that Hα emission and Li abundance not well correlated with age for low-mass stars between 20 and 100 Myr old. We argue that a finer age indicator, the ''LL-clock'', would be the luminosity at which the transition between efficient Li depletion and preservation takes place for fully convective objects. The LL-clock could allow in the near future to derive the relative ages of young open clusters, and clarify the study of PMS evolution of cool stars.

Pre-Main-Sequence Accretion in the Low Metallicity Galactic Star-forming Region Sh 2-284

We present optical spectra of pre-main-sequence (PMS) candidates around the Ha region taken with the Southern African Large Telescope in the low metallicity (Z) Galactic region Sh 2-284, which includes the open cluster Dolidze 25 with an atypical low metallicity of Z similar to 1/5 Z(circle dot). It has been suggested on the basis of both theory and observations that PMS mass-accretion rates, (M) over dot(acc), are a function of Z. We present the first sample of spectroscopic estimates of mass-accretion rates for PMS stars in any low-Z star-forming region. Our data set was enlarged with literature data of H alpha emission in intermediate-resolution R-band spectroscopy. Our total sample includes 24 objects spanning a mass range between 1 and 2 M-circle dot and with a median age of approximately 3.5 Myr. The vast majority (21 out of 24) show evidence for a circumstellar disk on the basis of Two Micron All Sky Survey and Spitzer infrared photometry. We find (M) over dot(acc) in the 1-2 M-circle dot interval to depend quasi-quadratically on stellarmass, with (M) over dot(acc) proportional to M-*(2.4 +/- 0.35), and inversely with stellar age, with (M) over dot(acc) proportional to t(*)(-0.7 +/- 0.4). Furthermore, we compare our spectroscopic (M) over dot(acc) measurements with solar Z Galactic PMS stars in the same mass range, but, surprisingly find no evidence for a systematic change in (M) over dot(acc) with Z. We show that literature accretion-rate studies are influenced by detection limits, and we suggest that (M) over dot(acc) may be controlled by factors other than Z(*), M-*, and age.

CHANDRA AND SPITZER IMAGING OF THE INFRARED CLUSTER IN NGC 2071

We present results of a sensitive Chandra X-ray observation and Spitzer mid-infrared (mid-IR) observations of the IR cluster lying north of the NGC 2071 reflection nebula in the Orion B molecular cloud. We focus on the dense cluster core known as NGC 2071-IR, which contains at least nine IR sources within a 40 `` x 40 `` region. This region shows clear signs of active star formation including powerful molecular outflows, Herbig-Haro objects, and both OH and H(2)O masers. We use Spitzer Infrared Array Camera (IRAC) images to aid in X-ray source identification and to determine young stellar object (YSO) classes using mid-IR colors. Spitzer IRAC colors show that the luminous source IRS 1 is a class I protostar. IRS 1 is believed to be driving a powerful bipolar molecular outflow and may be an embedded B-type star or its progenitor. Its X-ray spectrum reveals a fluorescent Fe emission line at 6.4 keV, arising in cold material near the protostar. The line is present even in the absence of large flares, raising questions about the nature of the ionizing mechanism responsible for producing the 6.4 keV fluorescent line. Chandra also detects X-ray sources at or near the positions of IRS 2, IRS 3, IRS 4, and IRS 6 and a variable X-ray source coincident with the radio source VLA 1, located just 2 `` north of IRS 1. No IR data are yet available to determine a YSO classification for VLA 1, but its high X-ray absorption shows that it is even more deeply embedded than IRS 1, suggesting that it could be an even younger, less-evolved protostar.

Low-resolution spectroscopy and spectral energy distributions of selected sources towards σ Orionis

Aims. We study in detail nine sources in the direction of the young σ Orionis cluster, which is considered to be a unique site for studying stellar and substellar formation. The nine sources were selected because of their peculiar properties, such as extremely-red infrared colours or excessively strong Hα emission for their blue optical colours. Methods. We acquired high-quality, low-resolution spectroscopy (R ∼ 500) of the nine targets with ALFOSC at the Nordic Optical Telescope. We also re-analysed [24]-band photometry from MIPS/Spitzer and compiled the highest quality photometric dataset available at the ViJHK_s passbands and the four IRAC/Spitzer channels, for constructing accurate spectral energy distributions between 0.55 and 24 μm. Results. The nine targets were classified into: one Herbig Ae/Be star with a scattering edge-on disc; two G-type stars; one X-ray flaring, early-M, young star with chromospheric Hα emission; one very low-mass, accreting, young spectroscopic binary; two young objects at the brown-dwarf boundary with the characteristics of classical T Tauri stars; and two emission-line galaxies, one undergoing star formation, and another whose spectral energy distribution is dominated by an active galactic nucleus. We also discovered three infrared sources associated with overdensities in a cold cloud of the cluster centre. Conclusions. Low-resolution spectroscopy and spectral energy distributions are a vital tool for measuring the physical properties and evolution of young stars and candidates in the σ Orionis cluster.

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.

DYNAMICAL EVOLUTION OF VISCOUS DISKS AROUND Be STARS. I. PHOTOMETRY

Be stars possess gaseous circumstellar disks that modify in many ways the spectrum of the central B star. Furthermore, they exhibit variability at several timescales and for a large number of observables. Putting the pieces together of this dynamical behavior is not an easy task and requires a detailed understanding of the physical processes that control the temporal evolution of the observables. There is an increasing body of evidence that suggests that Be disks are well described by standard alpha-disk theory. This paper is the first of a series that aims at studying the possibility of inferring several disk and stellar parameters through the follow-up of various observables. Here we study the temporal evolution of the disk density for different dynamical scenarios, including the disk buildup as a result of a long and steady mass injection from the star, the disk dissipation that occurs after mass injection is turned off, as well as scenarios in which active periods are followed by periods of quiescence. For those scenarios, we investigate the temporal evolution of continuum photometric observables using a three-dimensional non-LTE radiative transfer code. We show that light curves for different wavelengths are specific of a mass loss history, inclination angle, and alpha viscosity parameter. The diagnostic potential of those light curves is also discussed.