Galactic outflows and the pollution of the galactic environment by supernovae

We here explore the effects of the SN explosions into the halo of star-forming galaxies like the Milky Way. Successive randomly distributed and clustered SNe explosions cause the formation of hot superbubbles that drive either fountains or galactic winds above the galactic disk, depending on the amount and concentration of energy that is injected by the SNe. In a galactic fountain, the ejected gas is re-captured by the gravitational potential and falls back onto the disk. From 3D non-equilibrium radiative cooling hydrodynamical simulations of these fountains, we find that they may reach altitudes up to about 5 kpc in the halo and thus allow for the formation of the so called intermediate-velocity-clouds (IVCs) which are often observed in the halos of disk galaxies. The high-velocity-clouds that are also observed but at higher altitudes (of up to 12 kpc) require another mechanism to explain their production. We argue that they could be formed either by the capture of gas from the intergalactic medium and/or by the action of magnetic fields that are carried to the halo with the gas in the fountains. Due to angular momentum losses to the halo, we find that the fountain material falls back to smaller radii and is not largely spread over the galactic disk. Instead, the SNe ejecta fall nearby the region where the fountain was produced, a result which is consistent with recent chemical models of the galaxy. The fall back material leads to the formation of new generations of molecular clouds and to supersonic turbulence feedback in the disk. (C) 2009 COSPAR. Published by Elsevier Ltd. All rights reserved.

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.

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.

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.

TOWARD UNDERSTANDING THE B[e] PHENOMENON. IV. MODELING OF IRAS 00470+6429

FS CMa type stars are a recently described group of objects with the B[e] phenomenon which exhibits strong emission-line spectra and strong IR excesses. In this paper, we report the first attempt for a detailed modeling of IRAS 00470+6429, for which we have the best set of observations. Our modeling is based on two key assumptions: the star has a main-sequence luminosity for its spectral type (B2) and the circumstellar (CS) envelope is bimodal, composed of a slowly outflowing disklike wind and a fast polar wind. Both outflows are assumed to be purely radial. We adopt a novel approach to describe the dust formation site in the wind that employs timescale arguments for grain condensation and a self-consistent solution for the dust destruction surface. With the above assumptions we were able to satisfactorily reproduce many observational properties of IRAS 00470+6429, including the Hi line profiles and the overall shape of the spectral energy distribution. Our adopted recipe for dust formation proved successful in reproducing the correct amount of dust formed in the CS envelope. Possible shortcomings of our model, as well as suggestions for future improvements, are discussed.

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.

The effect of a single supernova explosion on the cuspy density profile of a small-mass dark matter halo

Some observations of galaxies, and in particular dwarf galaxies, indicate a presence of cored density profiles in apparent contradiction with cusp profiles predicted by dark matter N-body simulations. We constructed an analytical model, using particle distribution functions (DFs), to show how a supernova (SN) explosion can transform a cusp density profile in a small-mass dark matter halo into a cored one. Considering the fact that an SN efficiently removes matter from the centre of the first haloes, we study the effect of mass removal through an SN perturbation in the DFs. We find that the transformation from a cusp into a cored profile occurs even for changes as small as 0.5 per cent of the total energy of the halo, which can be produced by the expulsion of matter caused by a single SN explosion.

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.

Cross-entropy optimizer: a new tool to study precession in astrophysical jets

Evidence of jet precession in many galactic and extragalactic sources has been reported in the literature. Much of this evidence is based on studies of the kinematics of the jet knots, which depends on the correct identification of the components to determine their respective proper motions and position angles on the plane of the sky. Identification problems related to fitting procedures, as well as observations poorly sampled in time, may influence the follow-up of the components in time, which consequently might contribute to a misinterpretation of the data. In order to deal with these limitations, we introduce a very powerful statistical tool to analyse jet precession: the cross-entropy method for continuous multi-extremal optimization. Only based on the raw data of the jet components (right ascension and declination offsets from the core), the cross-entropy method searches for the precession model parameters that better represent the data. In this work we present a large number of tests to validate this technique, using synthetic precessing jets built from a given set of precession parameters. With the aim of recovering these parameters, we applied the cross-entropy method to our precession model, varying exhaustively the quantities associated with the method. Our results have shown that even in the most challenging tests, the cross-entropy method was able to find the correct parameters within a 1 per cent level. Even for a non-precessing jet, our optimization method could point out successfully the lack of precession.

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.

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.

MODELING VERY LONG BASELINE INTERFEROMETRIC IMAGES WITH THE CROSS-ENTROPY GLOBAL OPTIMIZATION TECHNIQUE

We present a new technique for obtaining model fittings to very long baseline interferometric images of astrophysical jets. The method minimizes a performance function proportional to the sum of the squared difference between the model and observed images. The model image is constructed by summing N(s) elliptical Gaussian sources characterized by six parameters: two-dimensional peak position, peak intensity, eccentricity, amplitude, and orientation angle of the major axis. We present results for the fitting of two main benchmark jets: the first constructed from three individual Gaussian sources, the second formed by five Gaussian sources. Both jets were analyzed by our cross-entropy technique in finite and infinite signal-to-noise regimes, the background noise chosen to mimic that found in interferometric radio maps. Those images were constructed to simulate most of the conditions encountered in interferometric images of active galactic nuclei. We show that the cross-entropy technique is capable of recovering the parameters of the sources with a similar accuracy to that obtained from the very traditional Astronomical Image Processing System Package task IMFIT when the image is relatively simple (e. g., few components). For more complex interferometric maps, our method displays superior performance in recovering the parameters of the jet components. Our methodology is also able to show quantitatively the number of individual components present in an image. An additional application of the cross-entropy technique to a real image of a BL Lac object is shown and discussed. Our results indicate that our cross-entropy model-fitting technique must be used in situations involving the analysis of complex emission regions having more than three sources, even though it is substantially slower than current model-fitting tasks (at least 10,000 times slower for a single processor, depending on the number of sources to be optimized). As in the case of any model fitting performed in the image plane, caution is required in analyzing images constructed from a poorly sampled (u, v) plane.

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.

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.