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.

Near-infrared integral field spectroscopy of the Homunculus nebula around eta Carinae using Gemini/CIRPASS

This work presents the first integral field spectroscopy of the Homunculus nebula around eta Carinae in the near-infrared spectral region (J band). We confirmed the presence of a hole on the polar region of each lobe, as indicated by previous near-IR long-slit spectra and mid-IR images. The holes can be described as a cylinder of height (i.e. the thickness of the lobe) and diameter of 6.5 and 6.0 x 10(16) cm, respectively. We also mapped the blue-shifted component of He I lambda 10830 seen towards the NW lobe. Contrary to previous works, we suggested that this blue-shifted component is not related to the Paddle but it is indeed in the equatorial disc. We confirmed the claim of N. Smith and showed that the spatial extent of the Little Homunculus matches remarkably well the radio continuum emission at 3 cm, indicating that the Little Homunculus can be regarded as a small H II region. Therefore, we used the optically thin 1.3 mm radio flux to derive a lower limit for the number of Lyman-continuum photons of the central source in eta Car. In the context of a binary system, and assuming that the ionizing flux comes entirely from the hot companion star, the lower limit for its spectral type and luminosity class ranges from O5.5 III to O7 I. Moreover, we showed that the radio peak at 1.7 arcsec NW from the central star is in the same line-of-sight of the `Sr-filament` but they are obviously spatially separated, while the blue-shifted component of He I lambda 10830 may be related to the radio peak and can be explained by the ultraviolet radiation from the companion star.

HI aperture synthesis and optical observations of the pair of galaxies NGC 6907 and 6908

NGC 6908, an S0 galaxy situated in the direction of NGC 6907, was only recently recognized as a distinct galaxy, instead of only a part of NGC 6907. We present 21-cm radio synthesis observations obtained with the Giant Metrewave Radio Telescope (GMRT) and optical images and spectroscopy obtained with the Gemini-North telescope of this pair of interacting galaxies. From the radio observations, we obtained the velocity field and the H I column density map of the whole region containing the NGC 6907/8 pair, and by means of the Gemini multi-object spectroscopy we obtained high-quality photometric images and 5 angstrom resolution spectra sampling the two galaxies. By comparing the rotation curve of NGC 6907 obtained from the two opposite sides around the main kinematic axis, we were able to distinguish the normal rotational velocity field from the velocity components produced by the interaction between the two galaxies. Taking into account the rotational velocity of NGC 6907 and the velocity derived from the absorption lines for NGC 6908, we verified that the relative velocity between these systems is lower than 60 km s(-1). The emission lines observed in the direction of NGC 6908, not typical of S0 galaxies, have the same velocity expected for the NGC 6907 rotation curve. Some emission lines are superimposed on a broader absorption profile, which suggests that they were not formed in NGC 6908. Finally, the H I profile exhibits details of the interaction, showing three components: one for NGC 6908, another for the excited gas in the NGC 6907 disc and a last one for the gas with higher relative velocities left behind NGC 6908 by dynamical friction, used to estimate the time when the interaction started in (3.4 +/- 0.6) x 10(7) yr ago.

Relativistic dispersion effects on SS433 jets

We here investigate the dispersion properties of radiation in the SS433 relativistic jets. We assume that the jet is composed of cold electron-proton plasma immersed in a predominantly parallel magnetic field to the jet axis. We find that for the mildly relativistic source SS433 (for which similar or equal to 79 degrees), the bulk velocity is too small (v similar or equal to 0.26c) to produce significant changes in the dispersion properties of the medium. Nonetheless, in the rarefied outer regions of the jets, where radio emission dominates, even a weak magnetic field has some influence on the dispersion properties and there appear two different electromagnetic branches that are slightly sensitive to the bulk relativistic motion. In the inner, X-ray region, the magnetic field is much stronger, but in this region the high electron density preserves the isotropic character of the local plasma and no branch separation occurs. In the region of the jet where the IR and optical emission dominates, the cold plasma may be also considered isotropic, i.e., neither the magnetic field nor the bulk velocity is able to affect the propagation of the radiation. Finally, we find that the Doppler line displacement in SS433 is affected by plasma dispersion only in a narrow frequency range in the far IR. As a consequence, although the shift (z) modulation due to precession of the SS433 jets is well described by previous work, it has to be corrected by plasma dispersion effects in the far-IR range.

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.

Gas streaming motions towards the nucleus of M81

We present two-dimensional stellar and gaseous kinematics of the inner 120 x 250 pc2 of the LINER/Seyfert 1 galaxy M81, from optical spectra obtained with the Gemini Multi-Object Spectrograph (GMOS) integral field spectrograph on the Gemini-North telescope at a spatial resolution of approximate to 10 pc. The stellar velocity field shows circular rotation and, overall, is very similar to the published large-scale velocity field, but deviations are observed close to the minor axis which can be attributed to stellar motions possibly associated with a nuclear bar. The stellar velocity dispersion of the bulge is 162 +/- 15 km s-1, in good agreement with previous measurements and leading to a black hole mass of M(BH) = 5.5+3.6(-2.0) x 107 M(circle dot) based on the M(BH)-Sigma relationship. The gas kinematics is dominated by non-circular motions and the subtraction of the stellar velocity field reveals blueshifts of approximate to-100 km s-1 on the far side of the galaxy and a few redshifts on the near side. These characteristics can be interpreted in terms of streaming towards the centre if the gas is in the plane. On the basis of the observed velocities and geometry of the flow, we estimate a mass inflow rate in ionized gas of approximate to 4.0 x 10-3 M(circle dot) yr-1, which is of the order of the accretion rate necessary to power the LINER nucleus of M81. We have also applied the technique of principal component analysis (PCA) to our data, which reveals the presence of a rotating nuclear gas disc within approximate to 50 pc from the nucleus and a compact outflow, approximately perpendicular to the disc. The PCA combined with the observed gas velocity field shows that the nuclear disc is being fed by gas circulating in the galaxy plane. The presence of the outflow is supported by a compact jet seen in radio observations at a similar orientation, as well as by an enhancement of the [O i]/H alpha line ratio, probably resulting from shock excitation of the circumnuclear gas by the radio jet. With these observations we are thus resolving both the feeding - via the nuclear disc and observed gas inflow, and the feedback - via the outflow, around the low-luminosity active nucleus of M81.

Are the magnetic fields of millisecond pulsars similar to 10(8) G?

It is generally assumed that the magnetic fields of millisecond pulsars (MSPs) are similar to 10(8) G. We argue that this may not be true and the fields may be appreciably greater. We present six evidences for this: (1) The similar to 10(8)G field estimate is based on magnetic dipole emission losses which is shown to be questionable; (2) The MSPs in low mass X-ray binaries (LMXBs) are claimed to have < 10(11) G on the basis of a Rayleygh-Taylor instability accretion argument. We show that the accretion argument is questionable and the upper limit 10(11) G may be much higher; (3) Low magnetic field neutron stars have difficulty being produced in LMXBs; (4) MSPs may still be accreting indicating a much higher magnetic field; (5) The data that predict similar to 10(8) G for MSPs also predict ages on the order of, and greater than, ten billion years, which is much greater than normal pulsars. If the predicted ages are wrong, most likely the predicted similar to 10(8) G fields of MSPs are wrong; (6) When magnetic fields are measured directly with cyclotron lines in X-ray binaries, fields a parts per thousand << 10(8) G are indicated. Other scenarios should be investigated. One such scenario is the following. Over 85% of MSPs are confirmed members of a binary. It is possible that all MSPs are in large separation binaries having magnetic fields > 10(8) G with their magnetic dipole emission being balanced by low level accretion from their companions.

Particle acceleration and magnetic field structure in PKS 2155-304: optical polarimetric observations

In this paper, we present multiband optical polarimetric observations of the very-high energy blazar PKS 2155-304 made simultaneously with a HESS/Fermi high-energy campaign in 2008, when the source was found to be in a low state. The intense daily coverage of the data set allowed us to study in detail the temporal evolution of the emission, and we found that the particle acceleration time-scales are decoupled from the changes in the polarimetric properties of the source. We present a model in which the optical polarimetric emission originates at the polarized mm-wave core and propose an explanation for the lack of correlation between the photometric and polarimetric fluxes. The optical emission is consistent with an inhomogeneous synchrotron source in which the large-scale field is locally organized by a shock in which particle acceleration takes place. Finally, we use these optical polarimetric observations of PKS 2155-304 at a low state to propose an origin for the quiescent gamma-ray flux of the object, in an attempt to provide clues for the source of its recently established persistent TeV 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.

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.

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.

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.

A faraway quasar in the direction of the highest energy Auger event

The highest energy cosmic ray event reported by the Auger Observatory has an energy of 148 EeV. It does not correlate with any nearby (z<0.024) object capable of originating such a high energy event. Intrigued by the fact that the highest energy event ever recorded (by the Fly`s Eye collaboration) points to a faraway quasar with very high radio luminosity and large Faraday rotation measurement, we have searched for a similar source for the Auger event. We find that the Auger highest energy event points to a quasar with similar characteristics to the one correlated to the Fly`s Eye event. We also find the same kind of correlation for one of the highest energy AGASA events. We conclude that so far these types of quasars are the best source candidates for both Auger and Fly`s Eye highest energy events. We discuss a few exotic candidates that could reach us from gigaparsec distances.