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.

THREE-DIMENSIONAL NUMERICAL SIMULATIONS OF MAGNETIZED WINDS OF SOLAR-LIKE STARS

By means of self-consistent three-dimensional magnetohydrodynamics (MHD) numerical simulations, we analyze magnetized solar-like stellar winds and their dependence on the plasma-beta parameter (the ratio between thermal and magnetic energy densities). This is the first study to perform such analysis solving the fully ideal three-dimensional MHD equations. We adopt in our simulations a heating parameter described by gamma, which is responsible for the thermal acceleration of the wind. We analyze winds with polar magnetic field intensities ranging from 1 to 20 G. We show that the wind structure presents characteristics that are similar to the solar coronal wind. The steady-state magnetic field topology for all cases is similar, presenting a configuration of helmet streamer-type, with zones of closed field lines and open field lines coexisting. Higher magnetic field intensities lead to faster and hotter winds. For the maximum magnetic intensity simulated of 20 G and solar coronal base density, the wind velocity reaches values of similar to 1000 km s(-1) at r similar to 20r(0) and a maximum temperature of similar to 6 x 10(6) K at r similar to 6r(0). The increase of the field intensity generates a larger ""dead zone"" in the wind, i.e., the closed loops that inhibit matter to escape from latitudes lower than similar to 45 degrees extend farther away from the star. The Lorentz force leads naturally to a latitude-dependent wind. We show that by increasing the density and maintaining B(0) = 20 G the system recover back to slower and cooler winds. For a fixed gamma, we show that the key parameter in determining the wind velocity profile is the beta-parameter at the coronal base. Therefore, there is a group of magnetized flows that would present the same terminal velocity despite its thermal and magnetic energy densities, as long as the plasma-beta parameter is the same. This degeneracy, however, can be removed if we compare other physical parameters of the wind, such as the mass-loss rate. We analyze the influence of gamma in our results and we show that it is also important in determining the wind structure.

Metallicity effects on the modified wind momentum of CSPN

Recent investigations on the central stars of planetary nebulae (CSPN) indicate that the masses based on model atmospheres can be much larger than the masses derived from theoretical mass-luminosity relations. Also, the dispersion in the relation between the modified wind momentum and the luminosity depends on the mass spread of the CSPN, and is larger than observed in massive hot stars. Since the wind characteristics probably depend on the metallicity, we analyze the effects on the modified wind momentum by considering the dispersion in this quantity caused by the stellar metallicity. Our CSPN masses are based on a relation between the core mass and the nebular abundances. We conclude that these masses agree with the known mass distribution both for CSPN and white dwarfs, and that the spread in the modified wind momentum can be explained by the observed metallicity variations.

Modelling the line variations from the wind-wind shock emissions of WR 30a

The study of Wolf-Rayet stars plays an important role in evolutionary theories of massive stars. Among these objects, similar to 20 per cent are known to be in binary systems and can therefore be used for the mass determination of these stars. Most of these systems are not spatially resolved and spectral lines can be used to constrain the orbital parameters. However, part of the emission may originate in the interaction zone between the stellar winds, modifying the line profiles and thus challenging us to use different models to interpret them. In this work, we analysed the He II lambda 4686 angstrom + C IV lambda 4658 angstrom blended lines of WR 30a (WO4+O5) assuming that part of the emission originate in the wind-wind interaction zone. In fact, this line presents a quiescent base profile, attributed to the WO wind, and a superposed excess, which varies with the orbital phase along the 4.6-d period. Under these assumptions, we were able to fit the excess spectral line profile and central velocity for all phases, except for the longest wavelengths, where a spectral line with constant velocity seems to be present. The fit parameters provide the eccentricity and inclination of the binary orbit, from which it is possible to constrain the stellar masses.

Revisiting 2D numerical models for the 19th century outbursts of eta Carinae

We present here new results of two-dimensional hydrodynamical simulations of the eruptive events of the 1840s (the great) and the 1890s (the minor) eruptions suffered by the massive star eta Carinae (Car). The two bipolar nebulae commonly known as the Homunculus and the little Homunculus (LH) were formed from the interaction of these eruptive events with the underlying stellar wind. We assume here an interacting, non-spherical multiple-phase wind scenario to explain the shape and the kinematics of both Homunculi, but adopt a more realistic parametrization of the phases of the wind. During the 1890s eruptive event, the outflow speed decreased for a short period of time. This fact suggests that the LH is formed when the eruption ends, from the impact of the post-outburst eta Car wind (that follows the 1890s event) with the eruptive flow (rather than by the collision of the eruptive flow with the pre-outburst wind, as claimed in previous models; Gonzalez et al.). Our simulations reproduce quite well the shape and the observed expansion speed of the large Homunculus. The LH (which is embedded within the large Homunculus) becomes Rayleigh-Taylor unstable and develop filamentary structures that resemble the spatial features observed in the polar caps. In addition, we find that the interior cavity between the two Homunculi is partially filled by material that is expelled during the decades following the great eruption. This result may be connected with the observed double-shell structure in the polar lobes of the eta Car nebula. Finally, as in previous work, we find the formation of tenuous, equatorial, high-speed features that seem to be related to the observed equatorial skirt of eta Car.

Precession and nutation in the eta Carinae binary system: evidence from the X-ray light curve

It is believed that eta Carinae is actually a massive binary system, with the wind-wind interaction responsible for the strong X-ray emission. Although the overall shape of the X-ray light curve can be explained by the high eccentricity of the binary orbit, other features like the asymmetry near periastron passage and the short quasi-periodic oscillations seen at those epochs have not yet been accounted for. In this paper we explain these features assuming that the rotation axis of eta Carinae is not perpendicular to the orbital plane of the binary system. As a consequence, the companion star will face eta Carinae on the orbital plane at different latitudes for different orbital phases and, since both the mass-loss rate and the wind velocity are latitude dependent, they would produce the observed asymmetries in the X-ray flux. We were able to reproduce the main features of the X-ray light curve assuming that the rotation axis of eta Carinae forms an angle of 29 degrees +/- 4 degrees with the axis of the binary orbit. We also explained the short quasi-periodic oscillations by assuming nutation of the rotation axis, with an amplitude of about 5 degrees and a period of about 22 days. The nutation parameters, as well as the precession of the apsis, with a period of about 274 years, are consistent with what is expected from the torques induced by the companion star.

Constraining the orbital orientation of eta Carinae from H Paschen lines

During the past decade, several observational and theoretical works have provided evidence of the binary nature of eta Carinae. Nevertheless, there is still no direct determination of the orbital parameters, and the different current models give contradictory results. The orbit is, in general, assumed to coincide with the Homunculus equator although the observations are not conclusive. Among all systems, eta Car has the advantage that it is possible to observe both the direct emission of line transitions in the central source and its reflection by the Homunculus, which is dependent on the orbital inclination. In this work, we studied the orbital phase-dependent hydrogen Paschen spectra reflected by the south-east lobe of the Homunculus to constrain the orbital parameters of eta Car and determine its inclination with respect to the Homunculus axis. Assuming that the emission excess originates in the wind-wind shock region, we were able to model the latitude dependence of the spectral line profiles. For the first time, we were able to estimate the orbital inclination of eta Car with respect to the observer and to the Homunculus axis. The best fit occurs for an orbital inclination to the line of sight of i similar to 60 degrees +/- 10 degrees, and i* similar to 35 degrees +/- 10 degrees with respect to the Homunculus axis, indicating that the angular momenta of the central object and the orbit are not aligned. We were also able to fix the phase angle of conjunction as similar to -40 degrees, showing that periastron passage occurs shortly after conjunction.

ON THE NATURE OF THE PROTOTYPE LUMINOUS BLUE VARIABLE AG CARINAE. II. WITNESSING A MASSIVE STAR EVOLVING CLOSE TO THE EDDINGTON AND BISTABILITY LIMITS

We show that the significantly different effective temperatures (T(eff)) achieved by the luminous blue variable AG Carinae during the consecutive visual minima of 1985-1990 (T(eff) similar or equal to 22,800 K) and 2000-2001 (T(eff) similar or equal to 17,000 K) place the star on different sides of the bistability limit, which occurs in line-driven stellar winds around T(eff) similar to 21,000 K. Decisive evidence is provided by huge changes in the optical depth of the Lyman continuum in the inner wind as T(eff) changes during the S Dor cycle. These changes cause different Fe ionization structures in the inner wind. The bistability mechanism is also related to the different wind parameters during visual minima: the wind terminal velocity was 2-3 times higher and the mass-loss rate roughly two times smaller in 1985-1990 than in 2000-2003. We obtain a projected rotational velocity of 220 +/- 50 km s(-1) during 1985-1990 which, combined with the high luminosity (L(star) = 1.5 x 10(6) L(circle dot)), puts AG Car extremely close to the Eddington limit modified by rotation (Omega Gamma limit): for an inclination angle of 90 degrees, Gamma(Omega) greater than or similar to 1.0 for M(circle dot) less than or similar to 60. Based on evolutionary models and mass budget, we obtain an initial mass of similar to 100 M(circle dot) and a current mass of similar to 60-70 M(circle dot) for AG Car. Therefore, AG Car is close to, if not at, the Omega Gamma limit during visual minimum. Assuming M = 70 M(circle dot), we find that Gamma(Omega) decreases from 0.93 to 0.72 as AG Car expands toward visual maximum, suggesting that the star is not above the Eddington limit during maximum phases.

Wolf-Rayet optically thick winds with Alfven waves

The Wolf-Rayet (WR) stars are hot luminous objects which are suffering an extreme mass loss via a continuous stellar wind. The high values of mass loss rates and high terminal velocities of the WR stellar winds constitute a challenge to the theories of radiation driven winds. Several authors incorporated magnetic forces to the line driven mechanism in order to explain these characteristics of the wind. Observations indicate that the WR stellar winds may reach, at the photosphere, velocities of the order of the terminal values, which means that an important part of the wind acceleration occurs at the optically thick region. The aim of this study is to analyze a model in which the wind in a WR star begins to be accelerated in the optically thick part of the wind. We used as initial conditions stellar parameters taken from the literature and solved the energy, mass and momentum equations. We demonstrate that the acceleration only by radiative forces is prevented by the general behavior of the opacities. Combining radiative forces plus a flux of Alfven waves, we found in the simulations a fast drop in the wind density profile which strongly reduces the extension of the optically thick region and the wind becomes optically thin too close its base. The understanding how the WR wind initiate is still an open issue. (C) 2010 COSPAR. Published by Elsevier Ltd. All rights reserved.

Alfven waves as a driving mechanism in stellar winds

Alfven waves have been invoked as an important mechanism of particle acceleration in stellar winds of cool stars. After their identification in the solar wind they started to be studied in winds of stars located in different regions of the FIR diagram. We discuss here some characteristics of these waves and we present a direct application in the acceleration of late-type stellar winds. (C) 2009 COSPAR. Published by Elsevier Ltd. All rights reserved.

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.

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.

Tomographic Simulations of accretion disks in cataclysmic variables - Flickering and wind

Both continuum and emission line flickering are phenomena directly associated with the mass-accretion process. In this work we simulated accretion-disk Doppler maps, including the effects of winds and flickering flares. Synthetic flickering Doppler maps were calculated and the effect of the flickering parameters on the maps was explored. Jets and winds occur in many astrophysical objects where accretion disks are present. Jets are generally absent among the cataclysmic variables (CVs), but there is evidence of mass loss by wind in many objects. CVs are ideal objects to study accretion disks, and consequently to study the wind associated with these disks. We also present simulations of accretion disks, including the presence of a wind with orbital phase resolution. Synthetic Ha line profiles in the optical region were obtained and their corresponding Doppler maps were calculated. The effect of the wind simulation parameters on the wind line profiles was also explored. From this study we verified that optically thick lines and/or emission by diffuse material into the primary Roche lobe are necessary to generate single peaked line profiles, often seen in CVs. The future accounting of these effects is suggested for interpreting Doppler tomography reconstructions.

On the mass distribution of neutron stars

The distribution of masses for neutron stars is analysed using the Bayesian statistical inference, evaluating the likelihood of the proposed Gaussian peaks by using 54 measured points obtained in a variety of systems. The results strongly suggest the existence of a bimodal distribution of the masses, with the first peak around 1.37 M(circle dot) and a much wider second peak at 1.73 M(circle dot). The results support earlier views related to the different evolutionary histories of the members for the first two peaks, which produces a natural separation (even if no attempt to `label` the systems has been made here). They also accommodate the recent findings of similar to M(circle dot) masses quite naturally. Finally, we explore the existence of a subgroup around 1.25 M(circle dot), finding weak, if any, evidence for it. This recently claimed low-mass subgroup, possibly related to the O-Mg-Ne core collapse events, has a monotonically decreasing likelihood and does not stand out clearly from the rest of the sample.