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.

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.

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.

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.

Verification of the role of the low level jets in Amazon squall lines

In this study we present a climatology of the Amazon squall lines (ASLs), between the years 2000 and 2008, using satellite imagery and European Centre for Medium-Range Weather Forecasts (ECMWF) reanalyses. The ASLs we are interested in are typically formed along the northern coast of Brazil and sometimes propagate for long distances inland. Results show that, on average, an ASL occurs every 2 days. ASLs are more frequent between April and June and less frequent between October and November. The years of 2005 and 2006 showed 25% more cases than the other years. This might be related to an increase of the Atlantic sea surface temperature. Of the total number of ASL cases, 54% propagated less than 170 km, 26% propagated between 170 and 400 km, and 20% propagated more than 400 km. We also studied the occurrence of low level jets (LLJs) associated with the coastal ASLs. Although LLJs are always present in the environment before the formation of the ASL and even on days without ASL cases, important differences were found, mainly related to the LLJ depths. (C) 2010 Elsevier B.V. All rights reserved.

GHASP: an H alpha kinematic survey of spiral and irregular galaxies - IX. The near-infrared, stellar and baryonic Tully-Fisher relations

We studied, for the first time, the near-infrared, stellar and baryonic Tully-Fisher relations for a sample of field galaxies taken from a homogeneous Fabry-Perot sample of galaxies [the Gassendi HAlpha survey of SPirals (GHASP) survey]. The main advantage of GHASP over other samples is that the maximum rotational velocities were estimated from 2D velocity fields, avoiding assumptions about the inclination and position angle of the galaxies. By combining these data with 2MASS photometry, optical colours, HI masses and different mass-to-light ratio estimators, we found a slope of 4.48 +/- 0.38 and 3.64 +/- 0.28 for the stellar and baryonic Tully-Fisher relation, respectively. We found that these values do not change significantly when different mass-to-light ratio recipes were used. We also point out, for the first time, that the rising rotation curves as well as asymmetric rotation curves show a larger dispersion in the Tully-Fisher relation than the flat ones or the symmetric ones. Using the baryonic mass and the optical radius of galaxies, we found that the surface baryonic mass density is almost constant for all the galaxies of this sample. In this study we also emphasize the presence of a break in the NIR Tully-Fisher relation at M(H,K) similar to -20 and we confirm that late-type galaxies present higher total-to-baryonic mass ratios than early-type spirals, suggesting that supernova feedback is actually an important issue in late-type spirals. Due to the well-defined sample selection criteria and the homogeneity of the data analysis, the Tully-Fisher relation for GHASP galaxies can be used as a reference for the study of this relation in other environments and at higher redshifts.

Can retired galaxies mimic active galaxies? Clues from the Sloan Digital Sky Survey

The classification of galaxies as star forming or active is generally done in the ([O III]/H beta, [N II]/H alpha) plane. The Sloan Digital Sky Survey (SDSS) has revealed that, in this plane, the distribution of galaxies looks like the two wings of a seagull. Galaxies in the right wing are referred to as Seyfert/LINERs, leading to the idea that non-stellar activity in galaxies is a very common phenomenon. Here, we argue that a large fraction of the systems in the right wing could actually be galaxies which stopped forming stars. The ionization in these `retired` galaxies would be produced by hot post-asymptotic giant branch stars and white dwarfs. Our argumentation is based on a stellar population analysis of the galaxies via our STARLIGHT code and on photoionization models using the Lyman continuum radiation predicted for this population. The proportion of LINER galaxies that can be explained in such a way is, however, uncertain. We further show how observational selection effects account for the shape of the right wing. Our study suggests that nuclear activity may not be as common as thought. If retired galaxies do explain a large part of the seagull`s right wing, some of the work concerning nuclear activity in galaxies, as inferred from SDSS data, will have to be revised.

Intragroup diffuse light in compact groups of galaxies - II. HCG 15, 35 and 51

This continuing study of intragroup light in compact groups of galaxies aims to establish new constraints to models of formation and evolution of galaxy groups, specially of compact groups, which are a key part in the evolution of larger structures, such as clusters. In this paper we present three additional groups (HCG 15, 35 and 51) using deep wide-field B- and R-band images observed with the LAICA camera at the 3.5-m telescope at the Calar Alto observatory (CAHA). This instrument provides us with very stable flat-fielding, a mandatory condition for reliably measuring intragroup diffuse light. The images were analysed with the OV_WAV package, a wavelet technique that allows us to uncover the intragroup component in an unprecedented way. We have detected that 19, 15 and 26 per cent of the total light of HCG 15, 35 and 51, respectively, are in the diffuse component, with colours that are compatible with old stellar populations and with mean surface brightness that can be its low as 28.4 B mag arcsec(-2). Dynamical masses, crossing times and mass-to-light ratios were recalculated using the new group parameters. Also tidal features were analysed using the wavelet technique.

XMM-Newton temperature maps for five intermediate redshift clusters of galaxies

We have analyzed XMM-Newton archive data for five clusters of galaxies (redshifts 0.223-0.313) covering a wide range of dynamical states, from relaxed objects to clusters undergoing several mergers. We present here temperature maps of the X-ray gas together with a preliminary interpretation of the formation history of these clusters. (c) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.

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.

Breaks in the radial oxygen abundance and corotation radius of three spiral galaxies

The correlation between the breaks in the metallicity distribution and the corotation radius of spiral galaxies has been already advocated in the past and is predicted by a chemodynamical model of our Galaxy that effectively introduces the role of spiral arms in the star formation rate. In this work, we present photometric and spectroscopic observations made with the Gemini Telescope for three of the best candidates of spiral galaxies to have the corotation inside the optical disc: IC 0167, NGC 1042 and NGC 6907. We observed the most intense and well-distributed H ii regions of these galaxies, deriving reliable galactocentric distances and oxygen abundances by applying different statistical methods. From these results, we confirm the presence of variations in the gradients of metallicity of these galaxies that are possibly correlated with the corotation resonance.

Morphological properties of superclusters of galaxies

We studied superclusters of galaxies in a volume-limited sample extracted from the Sloan Digital Sky Survey Data Release 7 and from mock catalogues based on a semi-analytical model of galaxy evolution in the Millennium Simulation. A density field method was applied to a sample of galaxies brighter than M(r) = -21+5 log h(100) to identify superclusters, taking into account selection and boundary effects. In order to evaluate the influence of the threshold density, we have chosen two thresholds: the first maximizes the number of objects (D1) and the second constrains the maximum supercluster size to similar to 120 h(-1) Mpc (D2). We have performed a morphological analysis, using Minkowski Functionals, based on a parameter, which increases monotonically from filaments to pancakes. An anticorrelation was found between supercluster richness (and total luminosity or size) and the morphological parameter, indicating that filamentary structures tend to be richer, larger and more luminous than pancakes in both observed and mock catalogues. We have also used the mock samples to compare supercluster morphologies identified in position and velocity spaces, concluding that our morphological classification is not biased by the peculiar velocities. Monte Carlo simulations designed to investigate the reliability of our results with respect to random fluctuations show that these results are robust. Our analysis indicates that filaments and pancakes present different luminosity and size distributions.

Nuclear and extended spectra of NGC 1068-II. Near-infrared stellar population synthesis

We performed stellar population synthesis on the nuclear and extended regions of NGC 1068 by means of near-infrared spectroscopy to disentangle their spectral energy distribution components. This is the first time that such a technique is applied to the whole 0.8-2.4 mu m wavelength interval in this galaxy. NGC 1068 is one of the nearest and probably the most studied Seyfert 2 galaxy, becoming an excellent laboratory to study the interaction between black holes, the jets that they can produce and the medium in which they propagate. Our main result is that traces of young stellar population are found at similar to 100 pc south of the nucleus. The contribution of a power-law continuum in the centre is about 25 per cent, which is expected if the light is scattered from a Seyfert 1 nucleus. We find peaks in the contribution of the featureless continuum about 100-150 pc from the nucleus on both sides. They might be associated with regions where the jet encounters dense clouds. Further support to this scenario is given by the peaks of hot dust distribution found around these same regions and the H(2) emission-line profile, leading us to propose that the peaks might be associated to regions where stars are being formed. Hot dust also has an important contribution to the nuclear region, reinforcing the idea of the presence of a dense, circumnuclear torus in this galaxy. Cold dust appears mostly in the south direction, which supports the view that the south-west emission is behind the plane of the galaxy and is extinguished very likely by dust in the plane. Intermediate-age stellar population contributes significantly to the continuum, especially in the inner 200 pc.

MAGNETIC FIELD EFFECTS NEAR THE LAUNCHING REGION OF ASTROPHYSICAL JETS

One of the fundamental properties of astrophysical magnetic fields is their ability to change topology through reconnection and in doing so, to release magnetic energy, sometimes violently. In this work, we review recent results on the role of magnetic reconnection and associated heating and particle acceleration in jet/accretion disk systems, namely young stellar objects (YSOs), microquasars, and active galactic nuclei (AGNs).