The SARS algorithm: detrending CoRoT light curves with Sysrem using simultaneous external parameters

Surveys for exoplanetary transits are usually limited not by photon noise but rather by the amount of red noise in their data. In particular, although the CoRoT space-based survey data are being carefully scrutinized, significant new sources of systematic noises are still being discovered. Recently, a magnitude-dependant systematic effect was discovered in the CoRoT data by Mazeh et al. and a phenomenological correction was proposed. Here we tie the observed effect to a particular type of effect, and in the process generalize the popular Sysrem algorithm to include external parameters in a simultaneous solution with the unknown effects. We show that a post-processing scheme based on this algorithm performs well and indeed allows for the detection of new transit-like signals that were not previously detected.

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

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

Alternative diagnostic diagrams and the `forgotten` population of weak line galaxies in the SDSS

A numerous population of weak line galaxies (WLGs) is often left out of statistical studies on emission-line galaxies (ELGs) due to the absence of an adequate classification scheme, since classical diagnostic diagrams, such as [O iii]/H beta versus [N ii]/H alpha (the BPT diagram), require the measurement of at least four emission lines. This paper aims to remedy this situation by transposing the usual divisory lines between star-forming (SF) galaxies and active galactic nuclei (AGN) hosts and between Seyferts and LINERs to diagrams that are more economical in terms of line quality requirements. By doing this, we rescue from the classification limbo a substantial number of sources and modify the global census of ELGs. More specifically, (1) we use the Sloan Digital Sky Survey Data Release 7 to constitute a suitable sample of 280 000 ELGs, one-third of which are WLGs. (2) Galaxies with strong emission lines are classified using the widely applied criteria of Kewley et al., Kauffmann et al. and Stasinska et al. to distinguish SF galaxies and AGN hosts and Kewley et al. to distinguish Seyferts from LINERs. (3) We transpose these classification schemes to alternative diagrams keeping [N ii]/H alpha as a horizontal axis, but replacing H beta by a stronger line (H alpha or [O ii]), or substituting the ionization-level sensitive [O iii]/H beta ratio with the equivalent width of H alpha (W(H alpha)). Optimized equations for the transposed divisory lines are provided. (4) We show that nothing significant is lost in the translation, but that the new diagrams allow one to classify up to 50 per cent more ELGs. (5) Introducing WLGs in the census of galaxies in the local Universe increases the proportion of metal-rich SF galaxies and especially LINERs. In the course of this analysis, we were led to make the following points. (i) The Kewley et al. BPT line for galaxy classification is generally ill-used. (ii) Replacing [O iii]/H beta by W(H alpha) in the classification introduces a change in the philosophy of the distinction between LINERs and Seyferts, but not in its results. Because the W(H alpha) versus [N ii]/H alpha diagram can be applied to the largest sample of ELGs without loss of discriminating power between Seyferts and LINERs, we recommend its use in further studies. (iii) The dichotomy between Seyferts and LINERs is washed out by WLGs in the BPT plane, but it subsists in other diagnostic diagrams. This suggests that the right wing in the BPT diagram is indeed populated by at least two classes, tentatively identified with bona fide AGN and `retired` galaxies that have stopped forming stars and are ionized by their old stellar populations.

A transiting giant planet with a temperature between 250 K and 430 K

Of the over 400 known(1) exoplanets, there are about 70 planets that transit their central star, a situation that permits the derivation of their basic parameters and facilitates investigations of their atmospheres. Some short-period planets(2), including the first terrestrial exoplanet(3,4) (CoRoT-7b), have been discovered using a space mission(5) designed to find smaller and more distant planets than can be seen from the ground. Here we report transit observations of CoRoT-9b, which orbits with a period of 95.274 days on a low eccentricity of 0.11 +/- 0.04 around a solar-like star. Its periastron distance of 0.36 astronomical units is by far the largest of all transiting planets, yielding a `temperate` photospheric temperature estimated to be between 250 and 430 K. Unlike previously known transiting planets, the present size of CoRoT-9b should not have been affected by tidal heat dissipation processes. Indeed, the planet is found to be well described by standard evolution models(6) with an inferred interior composition consistent with that of Jupiter and Saturn.

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.

A mechanism for the formation of oxygen and iron bimodal radial distribution in the disc of our Galaxy

Recently, it has been proposed that there are two type Ia supernova progenitors: short-lived and long-lived. On the basis of this idea, we develop a theory of a unified mechanism for the formation of the bimodal radial distribution of iron and oxygen in the Galactic disc. The underlying cause for the formation of the fine structure of the radial abundance pattern is the influence of the spiral arms, specifically the combined effect of the corotation resonance and turbulent diffusion. From our modelling, we conclude that in order to explain the bimodal radial distributions simultaneously for oxygen and iron and to obtain approximately equal total iron output from different types of supernovae, the mean ejected iron mass per supernova event should be the same as quoted in the literature if the maximum mass of stars, which eject heavy elements, is 50 M(circle dot). For the upper mass limit of 70 M(circle dot), the production of iron by a type II supernova explosion should increase by about 1.5 times.

Dynamic picture of the inner asteroid belt: implications for the density, size and taxonomic distributions of real objects

In this paper, we construct a dynamic portrait of the inner asteroidal belt. We use information about the distribution of test particles, which were initially placed on a perfectly rectangular grid of initial conditions, after 4.2 Myr of gravitational interactions with the Sun and five planets, from Mars to Neptune. Using the spectral analysis method introduced by Michtchenko et al., the asteroidal behaviour is illustrated in detail on the dynamical, averaged and frequency maps. On the averaged and frequency maps, we superpose information on the proper elements and proper frequencies of real objects, extracted from the data base, AstDyS, constructed by Milani and Knezevic. A comparison of the maps with the distribution of real objects allows us to detect possible dynamical mechanisms acting in the domain under study; these mechanisms are related to mean-motion and secular resonances. We note that the two- and three-body mean-motion resonances and the secular resonances (strong linear and weaker non-linear) have an important role in the diffusive transportation of the objects. Their long-lasting action, overlaid with the Yarkovsky effect, may explain many observed features of the density, size and taxonomic distributions of the asteroids.

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

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

The corotation gap in the Galactic H i distribution

We used the H i data from the LAB Survey to map the ring-shaped gap in H i density that lies slightly outside the solar circle. Adopting R(0) = 7.5 kpc, we find an average gap radius of 8.3 kpc and an average gap width of 0.8 kpc. The characteristics of the H i gap correspond closely to the expected ones, as predicted by theory and by numerical simulations of the gas flow near the corotation resonance.

Chemically tagging the Hyades stream: does it partly originate from the Hyades cluster?

The Hyades stream has long been thought to be a dispersed vestige of the Hyades cluster. However, recent analyses of the parallax distribution, of the mass function, and of the action-space distribution of stream stars have shown it to be rather composed of orbits trapped at a resonance of a density disturbance. This resonant scenario should leave a clearly different signature in the element abundances of stream stars than the dispersed cluster scenario, since the Hyades cluster is chemically homogeneous. Here, we study the metallicity as well as the element abundances of Li, Na, Mg, Fe, Zr, Ba, La, Ce, Nd and Eu for a random sample of stars belonging to the Hyades stream, and compare them with those of stars from the Hyades cluster. From this analysis: (i) we independently confirm that the Hyades stream cannot be solely composed of stars originating in the Hyades cluster; (ii) we show that some stars (namely 2/21) from the Hyades stream nevertheless have abundances compatible with an origin in the cluster; (iii) we emphasize that the use of Li as a chemical tag of the cluster origin of main-sequence stars is very efficient in the range 5500 K <= T(eff) <= 6200 K, since the Li sequence in the Hyades cluster is very tight, while at the same time spanning a large abundance range; (iv) we show that, while this evaporated population has a metallicity excess of similar to 0.2 dex with respect to the local thin-disc population, identical to that of the Hyades cluster, the remainder of the Hyades stream population has still a metallicity excess of similar to 0.06-0.15 dex, consistent with an origin in the inner Galaxy and (v) we show that the Hyades stream can be interpreted as an inner 4:1 resonance of the spiral pattern: this then also reproduces an orbital family compatible with the Sirius stream, and places the origin of the Hyades stream up to 1 kpc inwards from the solar radius, which might explain the observed metallicity excess of the stream population.

Hydrodynamical simulations of Galactic fountains - II. Evolution of multiple fountains

The ejection of gas out of the disc in late-type galaxies is related to star formation and is mainly due to the explosion of Type II supernovae (SN II). In a previous paper, we considered the evolution of a single Galactic fountain, that is, a fountain powered by a single SN cluster. Using three-dimensional hydrodynamical simulations, we studied in detail the fountain flow and its dependence with several factors, such as the Galactic rotation, the distance to the Galactic centre and the presence of a hot gaseous halo. As a natural followup, this paper investigates the dynamical evolution of multiple generations of fountains generated by similar to 100 OB associations. We have considered the observed size-frequency distribution of young stellar clusters within the Galaxy in order to appropriately fuel the multiple fountains in our simulations. Most of the results of the previous paper have been confirmed, like for example the formation of intermediate velocity clouds above the disc by the multiple fountains. Also, this work confirms the localized nature of the fountain flows: the freshly ejected metals tend to fall back close to the same Galactocentric region where they are delivered. Therefore, the fountains do not change significantly the radial profile of the disc chemical abundance. The multiple fountain simulations also allowed us to consistently calculate the feedback of the star formation on the halo gas. We found that the hot gas gains about 10 per cent of all the SN II energy produced in the disc. Thus, the SN feedback more than compensate for the halo radiative losses and allow for a quasi steady-state disc-halo circulation to exist. Finally, we have also considered the possibility of mass infall from the intergalactic medium and its interaction with the clouds that are formed by the fountains. Though our simulations are not suitable to reproduce the slow rotational pattern that is typically observed in the haloes around the disc galaxies, they indicate that the presence of an external gas infall may help to slow down the rotation of the gas in the clouds and thus the amount of angular momentum that they transfer to the coronal gas, as previously suggested in the literature.

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

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

Constraining the orbital orientation of eta Carinae from H Paschen lines

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

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.

Analysis of total visual and CCD V-broadband observations of Comet C/1995 O1 (Hale-Bopp): 1995-2001

We developed a general method for determination of water production rates from groundbased visual observations and applied it to Comet Hale-Bopp. Our main objective is to extend the method to include total visual magnitude observations obtained with CCD detector and V filter in the analysis of total visual magnitudes. We compare the CCD V-broadband careful observations of Liller [Liller, W. Pre-perihelion CCD photometry of Comet 1995 01 (Hale-Bopp). Planet. Space Sci. 45, 1505-1513, 1997; Liller, W. CCD photometry of Comet C/1995 O1 (Hale-Bopp): 1995-2000. Int. Comet Quart. 23(3), 93-97, 2001] with the total visual magnitude observations from experienced international observers found in the International Comet Quarterly (ICQ) archive. A data set of similar to 400 CCD observations covering about the same 6 years time span of the similar to 12,000 ICQ selected total visual magnitude observations were used in the analysis. A least-square method applied to the water production rates, yields power laws as a function of the heliocentric distances for the pre- and post-perihelion phases. The average dimension of the nucleus as well as its effective active area is determined and compared with values published in the literature. (C) 2009 COSPAR. Published by Elsevier Ltd. All rights reserved.