Since the start of ALMA observatory operation, new and important chemistry of infrared cold core was revealed. Molecular transitions at millimeter range are being used to identify and to characterize these sources. We have investigated the 231 GHz ALMA archive observations of the infrared dark cloud region C9, focusing on the brighter source that we called as IRDC-C9 Main. We report the existence of two sub-structures on the continuum map of this source: a compact bright spot with high chemistry diversity that we labelled as core, and a weaker and extended one, that we labelled as tail. In the core, we have identified lines of the molecules OCS(19-18), $^{13}$CS(5-4) and CH$_{3}$CH$_{2}$CN, several lines of CH$_{3}$CHO and the k-ladder emission of $^{13}$CH$_{3}$CN.We report two different temperature regions: while the rotation diagram of CH$_{3}$CHO indicates a temperature of 25 K, the rotation diagram of $^{13}$CH$_{3}$CN indicates a warmer phase at temperature of $\sim450$K. In the tail, only the OCS(19-18) and $^{13}$CS(5-4) lines were detected. We used the $Nautilus$ and the \textsc{Radex} codes to estimate the column densities and the abundances. The existence of hot gas in the core of IRDC-C9 Main suggests the presence of a protostar, which is not present in the tail.

We studied the correlation between brightness and polarization variations in 3C~279 at different wavelengths, over time intervals long enough to cover the time lags due to opacity effects. We used these correlations together with VLBI images to constrain the radio and high energy source position.We made 7 mm radio continuum and $R$-band polarimetric observations of 3C~279 between 2009 and 2014. The radio observations were performed at the Itapetinga Radio Observatory, while the polarimetric data were obtained at Pico dos Dias Observatory, both in Brazil. We compared our observations with the $\gamma$-ray Fermi/LAT and $R$-band SMARTS light curves. We found a good correlation between 7~mm and $R$-band light curves, with a delay of $170 \pm 30$ days in radio, but no correlation with the $\gamma$ rays. However, a group of several $\gamma$-ray flares in April 2011 could be associated with the start of the 7 mm strong activity observed at the end of 2011.We also detected an increase in $R$-band polarization degree and rotation of the polarization angle simultaneous with these flares. Contemporaneous VLBI images at the same radio frequency show two new strong components close to the core, ejected in directions very different from that of the jet.The good correlation between radio and $R$-band variability suggests that their origin is synchrotron radiation. The lack of correlation with $\gamma$-rays produced by the Inverse Compton process on some occasions could be due to the lack of low energy photons in the jet direction or to absorption of the high energy photons by the broad line region clouds. The variability of the polarization parameters during flares can be easily explained by the combination of the jet polarization parameters and those of newly formed jet components.

High resolution ALMA observations of the recent (2.52 yrs old) shell of Nova V5668 Sgr (2015) show a highly structured ionised gas distribution with small (10$^{15}$ cm) clumps. These are the smallest structures ever observed in the remnant of a stellar thermonuclear explosion. No extended contiguous emission could be found above the 2.5 $\sigma$ level in our data, while the peak hydrogen densities in the clumps reach 10$^6$ cm$^{-3}$. The millimetre continuum image suggests that large scale structures previously distinguished in other recent nova shells may result from the distribution of bright unresolved condensations.

The bipolar nebula Menzel 3 (Mz 3) was observed as part of the \textit{Herschel} Planetary Nebula Survey (\textit{HerPlaNS}), which used the PACS and SPIRE instruments aboard the \textit{Herschel Space Observatory} to study a sample of planetary nebulae (PNe). In this paper, one of the series describing \textit{HerPlaNS} results, we report the detection of H I recombination lines (HRLs) in the spectrum of Mz 3. Inspection of the spectrum reveals the presence of 12 HRLs in the 55 to 680 $\mu$m range covered by the PACS and SPIRE instruments (H11$\alpha$ to H21$\alpha$ and H14$\beta$). The presence of HRLs in this range is unusual for PNe and has not been reported in Mz 3 before. Our analysis indicates that the HRLs we observed are enhanced by laser effect occurring in the core of Mz 3. Our arguments for this are: (i) the available Mz 3 optical to submillimetre HRL $\alpha$ line intensity ratios are not well reproduced by the spontaneous emission of optically thin ionized gas, as would be typical for nebular gas in PNe; (ii) the compact core of Mz 3 is responsible for a large fraction of the Herschel HRLs emission; (iii) the line intensity ratios for Mz 3 are very similar to those in the core emission of the well known star MWC 349A, where laser effect is responsible for the enhancement of HRLs in the Herschel wavelength range; (iv) the physical characteristics relevant to cause laser effect in the core of MWC 349A are very similar to those in the core of Mz 3.

The recent discovery of a roughly simultaneous periodic variability in the light curves of the BL Lac object PG 1553+113 at several electromagnetic bands represents the first case of such odd behavior reported in the literature. Motivated by this, we analyzed 15 GHz interferometric maps of the parsec-scale radio jet of PG 1553+113 to verify the presence of a possible counterpart of this periodic variability. We used the Cross-Entropy statistical technique to obtain the structural parameters of the Gaussian components present in the radio maps of this source. We kinematically identified seven jet components formed coincidentally with flare-like features seen in the $\gamma$-ray light curve. From the derived jet component positions in the sky plane and their kinematics (ejection epochs, proper motions, and sky position angles), we modeled their temporal changes in terms of a relativistic jet that is steadily precessing in time. Our results indicate a precession period in the observer's reference frame of $2.24\pm0.03$ years, compatible with the periodicity detected in the light curves of PG 1553+113. However, the maxima of the jet Doppler boosting factor are systematically delayed relative to the peaks of the main $\gamma$-ray flares. We propose two scenarios that could explain this delay, both based on the existence of a supermassive black hole binary system in PG 1553+113. We estimated the characteristics of this putative binary system that also would be responsible for driving the inferred jet precession.

Aims: In this work, we analyse the multiwavelength brightness variations and flaring activity of FSRQ PKS1510-089, aiming to constrain the position of the emission sources. Methods: We report 7 mm (43 GHz) radio and R-band polarimetric observations of PKS1510-089. The radio observations were performed at the Itapetinga Radio Observatory, while the polarimetric data were obtained at the Pico dos Dias Observatory. The 7 mm observations cover the period between 2011 and 2013, while the optical polarimetric observations were made between 2009 and 2012. Results: At 7 mm, we detected a correlation between four radio and $\gamma$-ray flares with a delay of about 54 days between them; the higher frequency counterpart occurred first. Using optical polarimetry, we detected a large variation in polarization angle (PA) within two days associated with the beginning of a $\gamma$-ray flare. Complementing our data with other data obtained in the literature, we show that PA presented rotations associated with the occurrence of flares. Conclusions: Our results can be explained by a shock-in-jet model, in which a new component is formed in the compact core producing an optical and/or $\gamma$-ray flare, propagates along the jet, and after some time becomes optically thin and is detected as a flare at radio frequencies. The variability in the polarimetric parameters can also be reproduced; we can explain large variation in both PA and polarization degree (PD), in only one of them, or in neither, depending on the differences in PA and PD between the jet and the new component.

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

We report observations of $\eta$ Carinae obtained with ALMA in the continuum of 100, 230, 280 and 660 GHz in 2012 November, with a resolution that varied from 2.88" to 0.45" for the lower and higher frequencies respectively. The source is not resolved, even at the highest frequency; its spectrum is characteristic of thermal bremsstrahlung of a compact source, but different from the spectrum of optically thin wind. The recombination lines H42$\alpha$, He42$\alpha$, H40$\alpha$, He40$\alpha$, H50$\beta$, H28$\alpha$, He28$\alpha$, H21$\alpha$ and He21$\alpha$ were also detected and their intensities reveal non local thermodynamic equilibrium (NLTE) effects. We found that the line profiles could only be fit by an expanding shell of dense and ionized gas, which produces a slow shock in the surroundings of $\eta$ Carinae. Combined with fittings to the continuum, we were able to constrain the shell size, radius, density, temperature and velocity. The detection of the He recombination lines is compatible with the high temperature gas and requires a high energy ionizing photon flux, which must be provided by the companion star. The mass loss rate and wind velocity, necessary to explain the formation of the shell, are compatible with a LBV eruption. The position, velocity and physical parameters of the shell coincide with those of the Weigelt blobs. The dynamics found for the expanding shell corresponds to matter ejected by $\eta$ Carinae in 1941, in an event similar to that which formed the Little Homunculus ; for that reason we called the new ejecta the "Baby Homunculus".

We present 23 interferometric images of parsec-scale jet of the quasar PKS 1741--03 at 15, 24 and 43 GHz spanning about 13 yr. We model the images as a superposition of discrete two--dimensional elliptical Gaussian components, with parameters determined by the cross--entropy technique. All the images present a spatially unresolved component (core) and usually two or three components receding from it. The same components were found in simultaneous 24 and 43 GHz maps, showing the robustness of our model-fitting. The core-shift opacity effect between these frequencies is weak. We have identified seven components moving along straight lines at constant apparent superluminal speeds ($3.5\le\beta_\mathrm{obs}\le 6.1$), with different sky position angles ($-186\deg\le \eta\le -125\deg$). The core flux density tracks quite well the fluctuations seen in the historical single-dish light curve at 14.5 GHz, with no measurable delay. The total flux density from the moving jet components is delayed $\sim$2 yr in relation to the core light curve, roughly the same as the lag between the ejection epoch and the maximum flux density in the light curves of the jet components. We propose that there are three non-exclusive mechanisms for producing these delays. From the kinematics of the most robust jet components and the core brightness temperature, we determined the bulk Lorentz factor ($4.8\le\gamma\le 24.5$) and the jet viewing angle ($0.\deg 35 \le \theta \le 4.\deg 2$); these values agree with previous estimates from the spectral energy distribution of PKS 1741--03 and its radio variability.

We report 4 years of observations of 3C273 at 7 mm obtained with the Itapetinga Radiotelescope, in Brazil, between 2009 and 2013. We detected a flare in 2010 March, when the flux density increased by 50% and reached 35 Jy. After the flare, the flux density started to decrease and reached values lower than 10 Jy. We suggest that the 7 mm flare is the radio counterpart of the $\gamma$-ray flare observed by Fermi/LAT in 2009 September, in which the flux density at high energies reached a factor of fifty of its average value. A delay of 170 days between the radio and $\gamma$-ray flares was revealed using the Discrete Correlation Function (DCF) that can be interpreted in the context of a shock model, in which each flare corresponds to the formation of a compact superluminal component that expands and becomes optically thin at radio frequencies at latter epochs. The difference in flare intensity between frequencies and at a different times, is explained as a consequence of an increase in the Doppler factor $\delta$, as predicted by the 16 year precession model proposed by Abraham & Romero, which has a large effect on boosting at high frequencies while does not affect too much the observed optically thick radio emission. We discuss other observable effects of the variation in $\delta$, as the increase in the formation rate of superluminal components, the variations in the time delay between flares and the periodic behaviour of the radio light curve that we found compatible with changes in the Doppler factor.

We present the result of 6 years monitoring of SgrA*, radio source associated to the supermassive black hole at the centre of the Milky Way. Single dish observations were performed with the Itapetinga radio telescope at 7 mm, and the contribution of the SgrA complex that surrounds SgrA* was subtracted and used as instantaneous calibrator. The observations were alternated every 10 min with those of the HII region SrgB2, which was also used as a calibrator. The reliability of the detections was tested comparing them with simultaneous observations using interferometric techniques. During the observing period we detected a continuous increase in the SgrA* flux density starting in 2008, as well as variability in timescales of days and strong intraday fluctuations. We investigated if the continuous increase in flux density is compatible with free-free emission from the tail of the disrupted compact cloud that is falling towards SgrA* and concluded that the increase is most probably intrinsic to SgrA*. Statistical analysis of the light curve using Stellingwerf and Structure Function methods revealed the existence of two minima, 156 +/- 10 and 220 +/- 10 days. The same statistical tests applied to a simulated light curve constructed from two quadratic sinusoidal functions superimposed to random variability reproduced very well the results obtained with the real light curve, if the periods were 57 and 156 days. Moreover, when a daily sampling was used in the simulated light curve, it was possible to reproduce the 2.3 GHz structure function obtained by Falcke in 1999, which revealed the 57 days period, while the 106 periodicity found by Zhao et al in 2001 could be a resonance of this period.

BL Lacertae is the prototype of the BL Lac class of active galactic nuclei, exhibiting intensive activity on parsec (pc) scales, such as intense core variability and multiple ejections of jet components. In particular, in previous works the existence of precession motions in the pc-scale jet of BL Lacertae has been suggested. In this work we revisit this issue, investigating temporal changes of the observed right ascension and declination offsets of the jet knots in terms of our relativistic jet-precession model. The seven free parameters of our precession model were optimized via a heuristic cross-entropy method, comparing the projected precession helix with the positions of the jet components on the plane of the sky and imposing constraints on their maximum and minimum superluminal velocities. Our optimized best model is compatible with a jet having a bulk velocity of 0.9824c, which is precessing with a period of about 12.1 yr in the observer's reference frame and changing its orientation in relation to the line of sight between 4 and 5 degrees, approximately. Assuming that the jet precession has its origin in a supermassive binary black hole system, we show that the 2.3-yr periodic variation in the structural position angle of the very-long-baseline interferometry (VLBI) core of BL Lacertae reported by Stirling et al. is compatible with a nutation phenomenon if the secondary black hole has a mass higher than about six times that of the primary black hole.

We report a study of the stellar content of the Near-infrared cluster [DBS2003]\,157 embedded in the extended H\,{\sc ii} region GAL\,331.31-00.34, which is associated with the IRAS source 16085-5138. $JHK$ photometry was carried out in order to identify potential ionizing candidates, and the follow-up NIR spectroscopy allowed the spectral classification of some sources, including two O-type stars. A combination of NIR photometry and spectroscopy data was used to obtain the distance of these two stars, with the method of spectroscopic parallax: IRS\,298 (O6\,{\sc V}, $3.35 \pm 0.61$\,kpc) and IRS\,339 (O9\,{\sc V}, $3.24 \pm 0.56$\,kpc). Adopting the average distance of $3.29 \pm 0.58$\,kpc and comparing the Lyman continuum luminosity of these stars with that required to account for the radio continuum flux of the H\,{\sc ii} region, we conclude that these two stars are the ionizing sources of GAL\,331.31-00.34. Young stellar objects (YSOs) were searched by using our NIR photometry and MIR data from the GLIMPSE survey. The analysis of NIR and MIR colour-colour diagrams resulted in 47 YSO candidates. The GLIMPSE counterpart of IRAS\,16085-5138, which presents IRAS colour indices compatible with an ultra-compact H\,{\sc ii} region, has been identified. The analysis of its spectral energy distribution between 2 and $100\,\mu$m revealed that this source shows a spectral index $\alpha = 3.6$ between 2 and $25\,\mu$m, which is typical of a YSO immersed in a protostellar envelope. Lower limits to the bolometric luminosity and the mass of the embedded protostar have been estimated as $L=7.7\times10^3L_{\sun}$ and $M = 10\,M_{\sun}$, respectively, which corresponds to a B0--B1\,{\sc V} ZAMS star.

In the past few decades detailed observations of radio and X-rays emission from massive binary systems revealed a whole new physics present in such systems. Both thermal and non-thermal components of this emission indicate that most of the radiation at these bands originates in shocks. OB and WR stars present supersonic and massive winds that, when colliding, emit largely due to the free-free radiation. The non-thermal radio and X-ray emissions are due to synchrotron and inverse compton processes, respectively. In this case, magnetic fields are expected to play an important role on the emission distribution. In the past few years the modeling of the free-free and synchrotron emissions from massive binary systems have been based on purely hydrodynamical simulations, and ad hoc assumptions regarding the distribution of magnetic energy and the field geometry. In this work we provide the first full MHD numerical simulations of wind-wind collision in massive binary systems. We study the free-free emission characterizing its dependence on the stellar and orbital parameters. We also study self-consistently the evolution of the magnetic field at the shock region, obtaining also the synchrotron energy distribution integrated along different lines of sight. We show that the magnetic field in the shocks is larger than that obtained when the proportionality between $B$ and the plasma density is assumed. Also, we show that the role of the synchrotron emission relative to the total radio emission has been underestimated.

We present a new technique for obtaining model fittings to VLBI 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 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. We show that our technique is capable of recovering the parameters of the sources with a similar accuracy to that obtained from the traditional AIPS task IMFIT when the image is relatively simple (e.g., few components). For more complex 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 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.

Several galaxy clusters are known to present multiple and misaligned pairs of cavities seen in X-rays, as well as twisted kiloparsec-scale jets at radio wavelengths. It suggests that the AGN precessing jets play a role in the formation of the misaligned bubbles. Also, X-ray spectra reveal that typically these systems are also able to supress cooling flows, predicted theoretically. The absence of cooling flows in galaxy clusters has been a mistery for many years since numerical simulations and analytical studies suggest that AGN jets are highly energetic, but are unable to redistribute it at all directions. We performed 3D hydrodynamical simulations of the interaction between a precessing AGN jet and the warm intracluster medium plasma, which dynamics is coupled to a NFW dark matter gravitational potential. Radiative cooling has been taken into account and the cooling flow problem was studied. We found that precession is responsible for multiple pairs of bubbles, as observed. The misaligned bubbles rise up to scales of tens of kiloparsecs, where the thermal energy released by the jets are redistributed. After $\sim 150$ Myrs, the temperature of the gas within the cavities is kept of order of $\sim 10^7$ K, while the denser plasma of the intracluster medium at the central regions reaches $T \sim 10^5$ K. The existence of multiple bubbles, at diferent directions, result in an integrated temperature along the line of sight much larger than the simulations of non-precessing jets. This result is in agreement with the observations. The simulations reveal that the cooling flows cessed $\sim 50 - 70$ Myr after the AGN jets are started.

We present a map of the spiral structure of the Galaxy, as traced by molecular CS emission associated with IRAS sources which are believed to be compact HII regions. The CS line velocities are used to determine the kinematic distances of the sources, in order to investigate their distribution in the galactic plane. This allows us to use 870 objects to trace the arms, a number larger than that of previous studies based on classical HII regions. The distance ambiguity of the kinematic distances, when it exists, is solved by different procedures, including the latitude distribution and an analysis of the longitude-velocity diagram. The well defined spiral arms are seen to be confined inside the co-rotation radius, as is often the case in spiral galaxies. We identify a square-shaped sub-structure in the CS map with that predicted by stellar orbits at the 4:1 resonance (4 epicycle oscillations in one turn around the galactic center). The sub-structure is found at the expected radius, based on the known pattern rotation speed and epicycle frequency curve. An inner arm presents an end with strong inward curvature and intense star formation that we tentatively associate with the region where this arm surrounds the extremity of the bar, as seen in many barred galaxies. Finally, a new arm with concave curvature is found in the Sagitta to Cepheus region of the sky.

In this paper we present multiband optical polarimetric observations of the VHE blazar PKS 2155-304 made simultaneously with a H.E.S.S./Fermi high-energy campaign in 2008, when the source was found to be in a low state. The intense daily coverage of the dataset allowed us to study in detail the temporal evolution of the emission and we found that the particle acceleration timescales 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 polarised 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 organised 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.

The Perseus galaxy cluster is known to present multiple and misaligned pairs of cavities seen in X-rays, as well as twisted kiloparsec-scale jets at radio wavelengths; both morphologies suggest that the AGN jet is subject to precession. In this work we performed 3D hydrodynamical simulations of the interaction between a precessing AGN jet and the warm intracluster medium plasma, which dynamics is coupled to a NFW dark matter gravitational potential. The AGN jet inflates cavities that become buoyantly unstable and rise up out of the cluster core. We found that under certain circumstances precession can originate multiple pairs of bubbles. For the physical conditions in the Perseus cluster, multiple pairs of bubbles are obtained for a jet precession opening angle > 40 degrees acting for at least three precession periods, reproducing well both radio and X-ray maps. Based on such conditions, assuming that the Bardeen-Peterson effect is dominant, we studied the evolution of the precession opening angle of this system. We were able to constrain the ratio between the accretion disc and black hole angular momenta as 0.7 - 1.4. We were also able to constrain the present precession angle to 30 - 40 degrees, as well as the approximate age of the inflated bubbles to 100 - 150 Myrs.

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 with the axis of the binary orbit. We also explained the short quasi-periodic oscillations by assuming nutation of the rotation axis, with amplitude of about 5 degrees and 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.

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 optimisation. 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. Aiming to recover these parameters, we applied the cross-entropy method to our precession model, varying exhaustively the quantities associated to 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 1%-level. Even for a non-precessing jet, our optimization method could point out successfully the lack of precession.

During the past decade several observational and theoretical works provided evidences 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 Carinae 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 are dependent on the orbital inclination. In this work, we studied the orbital phase dependent hydrogen Paschen spectra reflected by the SE lobe of the Homunculus to constrain the orbital parameters of etaCar and determine its inclination with respect to the Homunculus axis. Assuming that the emission excess is originated 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 etaCar 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 ~ 60 degrees, and i* ~ 35 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 ~ -40 degrees, showing that the periastron passage occurs shortly after conjunction.

K-band spectra of young stellar candidates in four southern hemisphere clusters have been obtained with the near-infrared spectrograph GNIRS in Gemini South. The clusters are associated with IRAS sources that have colours characteristic of ultracompact HII regions. Spectral types were obtained by comparison of the observed spectra with those of a NIR library; the results include the spectral classification of nine massive stars and seven objects confirmed as background late-type stars. Two of the studied sources have K-band spectra compatible with those characteristic of very hot stars, as inferred from the presence of Civ, Niii, and Nv emission lines at 2.078 micron, 2.116 micron, and 2.100 micron respectively. One of them, I16177 IRS1, has a K-band spectrum similar to that of Cyg OB2 7, an O3If* supergiant star. The nebular K-band spectrum of the associated UC Hii region shows the s-process [Kriii] and [Seiv] high excitation emission lines, previously identified only in planetary nebula. One young stellar object (YSO) was found in each cluster, associated with either the main IRAS source or a nearby resolved MSX component, confirming the results obtained from previous NIR photometric surveys. The distances to the stars were derived from their spectral types and previously determined JHK magnitudes; they agree well with the values obtained from the kinematic method, except in the case of IRAS15408-5356, for which the spectroscopic distance is about a factor two smaller than the kinematic value.

The study of Wolf-Rayet stars plays an important role in evolutionary theories of massive stars. Among these objects, ~ 20% 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 analyzed the HeII4686\AA + CIV4658\AA blended lines of WR30a (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 day 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.

Strong evidence for the presence of a warped Keplerian accretion disc in NGC4258 (M 106) has been inferred from the kinematics of water masers detected at sub-parsec scales. Assuming a power-law accretion disc and using constraints on the disc parameters derived from observational data, we have analyzed the relativistic Bardeen-Petterson effect driven by a Kerr black hole as the potential physical mechanism responsible for the disc warping. We found that the Bardeen-Petterson radius is comparable to or smaller than the inner radius of the maser disc (independent of the allowed value for the black hole spin parameter). Numerical simulations for a wide range of physical conditions have shown that the evolution of a misaligned disc due to the Bardeen-Petterson torques usually produces an inner flat disc and a warped transition region with a smooth gradient in the tilt and twist angles. Since this structure is similar to that seen in NGC 4258, we propose that the Bardeen-Petterson effect may be responsible for the disc warping in this galaxy. We estimated the time-scale necessary for the disc inside of the Bardeen-Petterson radius to align with the black hole's equator, as a function of the black hole spin. Our results show that the Bardeen-Petterson effect can align the disc within a few billion years in the case of NGC 4258. Finally, we show that if the observed curvature of the outer anomalous arms in the galactic disc of NGC 4258 is associated with the precession of its radio jet/counterjet, then the Bardeen-Petterson effect can provide the required precession period.

We modeled the HeII 4686 line profiles observed in the eta Carinae binary system close to the 2003.5 spectroscopic event, assuming that they were formed in the shocked gas that flows at both sides of the contact surface formed by wind-wind collision. We used a constant flow velocity and added turbulence in the form of a gaussian velocity distribution. We allowed emission from both the primary and secondary shocks but introduced infinite opacity at the contact surface, implying that only the side of the contact cone visible to the observer contributed to the line profile. Using the orbital parameters of the binary system derived from the 7 mm light curve during the last spectroscopic event (Paper II) we were able to reproduce the line profiles obtained with the HST at different epochs, as well as the line mean velocities obtained with ground based telescopes. A very important feature of our model is that the line profile depends on the inclination of the orbital plane; we found that to explain the latitude dependent mean velocity of the line, scattered into the line of sight by the Homunculus, the orbit inclination should be close to 90 degrees, meaning that it does not lie in the Homunculus equatorial plane, as usually assumed. This inclination, together with the relative position of the stars during the spectroscopic events, allowed us to explain most of the observational features, like the variation of the Purple Haze with the orbital phase, and to conciliate the X-ray absorption with the postulated shell effect used to explain the optical and UV light curves.

The Bardeen-Petterson general relativistic effect has been suggested as the mechanism responsible for precession in some accretion disk systems. Here we examine separately four mechanisms (tidally-induced, irradiation-induced, magnetically-induced and Bardeen-Petterson-induced) that can lead to warping and precession. We use a sample of eight X-ray binaries and four Active Galactic Nuclei (AGNs) that present signatures of warping and/or precession in their accretion disks to explore the viability of the different mechanisms. For the X-ray binaries SMC X-1 and 4U 1907+09 all four mechanisms provide precession periods compatible with those observed, while for Cyg X-1 and the active galaxies Arp 102B and NGC 1068, only two mechanisms are in agreement with the observations. The irradiation-driven instability seems incapable of producing the inferred precession of the active galaxies in our sample, and the tidally-induced precession can probably be ruled out in the case of Arp 102B. Perhaps the best case for a Bardeen-Petterson precession can be achieved for NGC 1068. Our results show that given the many observational uncertainties that still exist, it is extremely difficult to confirm unambiguously that the Bardeen-Petterson effect has been observed in any of the other sources of our sample.

One of the most intriguing spectral features of WR binary stars is the presence of time-dependent line profiles. Long term observations of several systems revealed the periodicity of this variability, synchronized with the orbital movement. Several partially successful models have been proposed to reproduce the observed data. The most promising assume that the origin of the emission is the wind-wind interaction zone. In this scenario, two high velocity and dense winds produce a strong shock layer, responsible for most of the X-rays observed from these systems. As the secondary star moves along its orbital path, the shock region of conical shape, changes its position with relation to the line of sight. As a consequence, the stream measured Doppler shift presents time variations resulting in position changes of the spectral line. In our work, we present an alternative model, introducing turbulence in the shock layer to account for the line broadening and opacity effects for the asymmetry in the line profiles. We showed that the gas turbulence avoids the need of an unnaturally large contact layer thickness to reproduce line broadening. Also, we demonstrated that if the emission from the opposing cone surface is absorbed, the result is a single peaked profile. This result fully satisfies the recent data obtained from massive binary systems, and can help on the determination of both winds and orbital parameters. We successfully applied this model to the Br22 system and determined its orbital parameters.

We report NIR broad and narrow band photometric observations in the direction of the IRAS17149-3916 source that reveal the presence of a young cluster of massive stars embedded in an HII region coincident with RCW121. These observations, together with published radio data, MSX and Spitzer images were used to determine some of the physical parameters of the region. We found 96 cluster member candidates in an area of about 1.5 x 2.0 square arcmin, 30% of them showing excess emission in the NIR. IRS 1, the strongest source in the cluster with an estimated spectral type of O5V-O6V ZAMS based on the color-magnitude diagram, is probably the main ionizing source of the HII region detected at radio wavelengths. Using the integrated Brgamma and the 5 GHz flux densities, we derived a mean visual extinction AV=5.49 magnitudes. From the observed size of the Brgamma extended emission, we calculated the emission measure E and the electron density ne, characteristic of compact HII regions.

We report the discovery of a small cluster of massive stars embedded in a NIR nebula in the direction of the IRAS15411-5352 point source, which is related to the alleged planetary nebula W16-185. The majority of the stars present large NIR excess characteristic of young stellar objects and have bright counterparts in the Spitzer IRAC images; the most luminous star (IRS1) is the NIR counterpart of the IRAS source. We found very strong unresolved Brgamma emission at the IRS1 position and more diluted and extended emission across the continuum nebula. From the sizes and electron volume densities we concluded that they represent ultra-compact and compact HII regions, respectively. Comparing the Brgamma emission with the 7 mm free-free emission, we estimated that the visual extinction ranges between 14 and 20 mag. We found that only one star (IRS1) can provide the number of UV photons necessary to ionize the nebula.

VLBA high spatial resolution observations of the disk structure of the active galactic nucleus NGC 1068 has recently revealed that the kinematics and geometry of this AGN is well characterized by an outer disk of H2O maser emission having a compact milliarcsecond (parsec) scale structure, which is encircling a thin rotating inner disk surrounding a ~10^7 M$_\sun$ compact mass, likely a black hole. A curious feature in this source is the occurrence of a misalignment between the inner and outer parts of the disk, with the galaxy's radio jet being orthogonal to the inner disk. We interpret this peculiar configuration as due to the Bardeen-Petterson effect, a general relativistic effect that warps an initially inclined (to the black hole equator) viscous disk, and drives the angular momentum vector of its inner part into alignment with the rotating black hole spin. We estimate the time-scale for both angular momenta to get aligned as a function the spin parameter of the Kerr black hole. We also reproduce the shape of the parsec and kiloparsec scale jets, assuming a model in which the jet is precessing with a period and aperture angle that decrease exponentially with time, as expected from the Bardeen-Petterson effect.

In this paper we use the 7 mm and 1.3 mm light curves obtained during the 2003.5 low excitation phase of the eta Carinae system to constrain the possible parameters of the binary orbit. To do that we assumed that the mm wave emission is produced in a dense disk surrounding the binary system; during the low excitation phase, which occurs close to periastron, the number of ionizing photons decreases, producing the dip in the radio emission. On the other hand, due to the large eccentricity, the density of the shock region at periastron is very high and the plasma is optically thick for free-free radiation at 7 mm, explaining the sharp peak that was observed at this frequency and lasted for about 10 days. From the shape and duration of the peak we were able to determine the orbital parameters of the binary system, independently of the stellar parameters, such as mass loss rates, wind velocities or temperature at the post-shock region.

In this paper we present observations of eta Carinae in the 1.3 mm and 7 mm radio continuum, during the 2003.5 low excitation phase. The expected minimum in the light curves was confirmed at both wavelengths and was probably due to a decrease in the number of UV photons available to ionize the gas surrounding the binary system. At 7 mm a very well defined peak was superimposed on the declining flux density. It presented maximum amplitude in 29 June 2003 and lasted for about 10 days. We show that its origin can be free-free emission from the gas at the shock formed by wind-wind collision, which is also responsible for the observed X-ray emission. Even though the shock strength is strongly enhanced as the two stars in the binary system approach each other, during periastron passage the X-ray emission is strongly absorbed and the 7 mm observations represent the only direct evidence of this event.

The exact nature of eta Carinae is still an open issue. In this paper we assume a binary system to explain the strong X-ray emission, but we also take into account that, near periastron and because of the highly eccentric orbit, the wind emerging from eta Carinae accumulates behind the shock and can mimic a shell-like ejection event. For this process to be effective, at periastron the secondary star should be located between eta Carinae and the observer, solving also the discrepancy between the orbital parameters derived from ground and space based observations. We show that, as the secondary moves in its orbit, the shell cools down and the number of available stellar ionizing photons is not enough to maintain the shell temperature at its equilibrium value of about 7500 K. The central part of the shell remains cold and under these conditions grain formation and growth can take place in timescales of hours. Using recently published HST data for the optical and UV continuum during the 2003 shell event, we estimated that a fraction of 0.01 of the interstellar dust to H ratio is enough to explain the observations. We also calculated the neutral gas column density intercepting the line of sight at each point of the orbit near periastron, and were able to reproduce the form and duration of the X-ray light curve without any change in the eta Carinae mass loss rate. This same column density can explain the observed H$\alpha$ light curve observed during the 2003 event.

We show that it is possible to explain the physical origin of jet precession in active galactic nuclei (AGNs) through the misalignment between the rotation axes of the accretion disk and of the Kerr black hole. We apply this scenario to quasars, Seyfert galaxies and also to the Galactic Center black hole Sgr A*, for which signatures of either jet or disk precession have been found. The formalism adopted is parameterized by the ratio of the precession period to the black hole mass and can be used to put constraints to the physical properties of the accretion disk as well as to the black hole spin in those systems.

We present optical light curves for the period 1996-2000, of two of the brightest known EGRET BL Lac objects: PKS 2005-489 and PKS 2155-304, the latter also one of the few known TeV sources. For both objects, quiescent epochs of low level of variability were followed by active periods, without any indication of periodicity. In PKS 2005-489, several variability events with duration of about 20 days were observed. In PKS 2155-304 fast drops and subsequent rises in luminosity occurred in time scales of days. We proposed an explanation in which a region moving along the relativistic jet is eclipsed by broad line region clouds or star clusters in the host galaxy. We compare our light curves with contemporaneous X-ray observations from All-Sky Monitor/RXTE. Correlations between optical and X-ray activity were not found in any of the sources at long time scales. However in PKS 2005-489 possible correlations were observed in 1997 and 1998 at short time scales, with optical variability preceding X-rays by 30 days in 1997 and succeeding them by about 10 days in 1998. The analysis of the SED, using the optical data presented here and BeppoSAX contemporaneous observations obtained from the literature, shows only small shifts in the synchrotron peak as the X-ray flux density changes.

We report the discovery of a young massive stellar cluster and infrared nebula in the direction of the CS molecular cloud associated to the IRAS point source 16132-5039. The analysis of the mid-infrared images from the more accurate MSX catalog, reveled that there are two independent components associated with the IRAS source. The integral of the spectral energy distribution for these components, between 8.28 $\mu$m and 100 $\mu$m, gave lower limits for the bolometric luminosity of the embedded objects of $8.7 \times 10^4 L_\odot$ and $9 \times 10^3 L_\odot$, which corresponds to ZAMS O8 and B0.5 stars, respectively. The number of Lyman continuum photons expected from the stars that lie along the reddening line for early-type stars is about 1.7 $ \times$ 10$^{49}$ s$^{-1}$, enough to produce the detected flux densities at 5 GHz. The NIR spectrum of the nebula increases with frequency, implying that free-free emission cannot be the main source of the extended luminosity, from which we conclude that the observed emission must be mainly dust scattered light. A comparison of the cluster described in this paper with the young stellar cluster associated with the IRAS source 16177-5018, which is located at the same distance and direction, shows that the mean visual absorption of the newly discovered cluster is about 10 magnitudes smaller and it contains less massive stars, suggesting that it was formed from a less massive molecular cloud.

Optical variability of 3C 120 is discussed in the framework of jet precession. Specifically, we assume that the observed long-term periodic variability is produced by the emission from an underlying jet with a time-dependent boosting factor driven by precession. The differences in the apparent velocities of the different superluminal components in the milliarcsecond jet can also be explained by the precession model as being related to changes in the viewing angle. The evolution of the jet components has been used to determine the parameters of the precession model, which also reproduce the helical structure seen at large scales. Among the possible mechanisms that could produce jet precession, we consider that 3C 120 harbours a super-massive black hole binary system in its nuclear region and that torques induced by misalignment between the accretion disc and the orbital plane of the secondary black hole are responsible for this precession; we estimated upper and lower limits for the black holes masses and their mean separation.

VLBI observations have shown that the parsec-jet of 3C 345 is formed by several components, ejected from the core with superluminal velocities and travelling along bent trajectories on the plane of the sky. We interpret the differences in velocity and position angle among the different features at formation time as the result of parsec-scale precession of the relativistic jet and calculate the aperture angle of the precession cone, the angle between the cone axis and the line of sight and the Lorentz factor associated with the jet bulk motion. We assumed a precession period of 10.1 yr, which is one of the B-band light curve long-term periods reported in the literature. We propose that boosting of the underlying jet emission, which is time-dependent due to precession, is responsible for this long-term optical variability. Jet precession with periods of several years can be produced in super-massive black hole binary systems, when the secondary black hole is in an orbit non-coplanar with the primary accretion disk, inducing torques in the inner parts of the disk. Assuming that this mechanism is responsible for the jet precession in 3C 345, we estimate upper and lower limits for the masses of the two black holes, as well as their mean separation. We found a correlation between the formation of jet components and the occurrence of strong optical flares, as well as a very strong anti-correlation between the intensity of these flares and the time required for the components to reach the maximum flux density at radio frequencies.

We present results of a microvariability polarization study in the violently variable quasar 3C279. We have resolved the polarization curves in the V band for this object down to timescales of minutes. We found two main components in the evolution of the degree of linear polarization, one consisting of a flicker with timescales of several tens of minutes and other component with far more significant variations on timescales of a few days. The linear polarization descended from $\sim 17$ % down to $\sim 8$ % in three nights. The polarization angle underwent a sudden change of more that 10 degrees in a few hours, perhaps indicating the injection of a new shock in the jet. The amplitude of the intranight flickering in the degree of polarization is at the level of $\sim 1$%. These are probably the best sampled polarization data ever obtained for this object. We also performed IR observations and we provide a follow-up of the evolution of this source at such energies after the main polarization outburst.

We report the discovery of a young massive stellar cluster embedded in an extended HII region, invisible at optical wavelengths where the extinction is $A_V \approx 28$ magnitudes, associated with the IRAS source 16177-5018. $J, H$ and nb$K$ imaging photometry combined with the $K_S$ 2MASS data show the presence of sources with infrared excess emission at 2.2 $\mu$m, concentrated in an area of about one square parsec around a massive young stellar object identified as the IRAS source. This object has a near-mid infrared spectral index betweem 2.2 and 25 $\mu$m $\alpha({\rm IR}) = d {\rm log}(\lambda F_\lambda)/d {\rm log} \lambda $ =4.78, characteristic of compact H II regions, with bolometric luminosity, inferred from the integrated near to far-infrared flux density of $2.8 \times 10^5 L_\odot$, which corresponds to a ZAMS star of about $42 M_\odot$. From the color-magnitude diagram we were able to classify the majority of the cluster members as reddened massive stars earlier than spectral type B5.

Recent IR observations indicate that many massive binary systems present dust formation episodes, in regions close to the stars, during the periastron passage. These systems are known to be high-energy sources, and it is believed that wind collisions are the origin of the emission. In this work we show that wind collisions not only increase the X-ray emission but also allow dust formation. As an application we study eta Car, which presents, near periastron, an increase in the X-ray emission followed by a sudden decrease that lasts for about a month. We reproduce this feature calculating the optical depth due to dust formation along the orbital period.

Massive binary star systems are not uncommon, and neither the supersonic collision of their winds. In the present work we study these shocks and the further consequences on wind structure. The post-shock gas is a warm and high-density environment, which allows dust to form and grow. We show that this growth is fast, of just a few hours. An application for etaCarinae shows that, probably, the decline of X-rays fluxes observed in its light curve is the consequence of its high absorption in periodic dust formation events, on the periastron passage.

The region centered in the NGC 6334 I(N) radio continuum source was surveyed in an extension of 6 arcmin in right ascension and 12 arcmin in declination, in the NH3(J,K) = (1,1) transition, using the Itapetinga radio telescope. The spectra show non-LTE behavior, and gradients of velocity and line-width were detected along the region. A detailed analysis of the spectra showed that the surveyed region is composed of at least three overlapped sources related to regions that are in different stages of star formation: NGC 6334 I, associated with an already known molecular bipolar outflow, NGC 6334 I(N)w, the brightest ammonia source, coincidental with the continuum source NGC 6334 I(N), and NGC 6334 I(N)e, weaker, more extended and probably less evolved than the others. The physical parameters of the last two sources were calculated in non-LTE conditions, assuming that their spectra are the superposition of the narrow line spectra produced by small dense clumps. The H2 density, NH3 column density, kinetic temperature, diameter and mass of the clumps were found to be very similar in the two regions, but the density of clumps is lower in the probably less evolved source NGC 6334 I(N)e. Differences between the physical parameters derived assuming LTE and non-LTE conditions are also discussed in this work.