Numerical radiative transfer with state-of-the-art iterative methods made easy

This article presents an on-line tool (rttools.irap.omp.eu) and its accompanying software ressources for the numerical solution of basic radiation transfer out of local thermodynamic equilibrium (LTE). State-of-the-art stationary iterative methods such as Accelerated $\Lambda$-Iteration and Gauss-Seidel schemes, using a short characteristics-based formal solver are used. We also comment on typical numerical experiments associated to the basic non-LTE radiation problem. These ressources are intended for the largest use and benefit, in support to more classical radiation transfer lectures usually given at the Master level.Comment: 8 pages, 5 figures, accepted for Eur. J. Phys. - see also (and use!) http://rttools.irap.omp.e

A search for CO+ in planetary nebulae

We have carried out a systematic search for the molecular ion CO+ in a sample of 8 protoplanetary and planetary nebulae in order to determine the origin of the unexpectedly strong HCO+ emission previously detected in these sources. An understanding of the HCO+ chemistry may provide direct clues to the physical and chemical evolution of planetary nebulae. We find that the integrated intensity of the CO+ line may be correlated with that of HCO+, suggesting that the reaction of CO+ with molecular hydrogen may be an important formation route for HCO+ in these planetary nebulae.Comment: 6 pages, 4 figures, accepted for publication in MNRA

Dynamical localization of matter wave solitons in managed barrier potentials

The bright matter wave soliton propagation through a barrier with a rapidly oscillating position is investigated. The averaged over rapid oscillations Gross-Pitaevskii (GP) equation is derived. It is shown that the soliton is dynamically trapped by the effective double-barrier. The analytical predictions for the soliton effective dynamics is confirmed by the numerical simulations of the full GP equation.Comment: 10 pages, 6 figure

VLTI/AMBER spectro-interferometric imaging of VX Sgr's inhomogenous outer atmosphere

Aims. We aim to explore the photosphere of the very cool late-type star VX Sgr and in particular the existence and characterization of molecular layers above the continuum forming photosphere. Methods. We obtained interferometric observations with the VLTI/AMBER interferometer using the fringe tracker FINITO in the spectral domain 1.45-2.50 micron with a spectral resolution of about 35 and baselines ranging from 15 to 88 meters.We perform independent image reconstruction for different wavelength bins and fit the interferometric data with a geometrical toy model.We also compare the data to 1D dynamical models of Miras atmosphere and to 3D hydrodynamical simulations of red supergiant (RSG) and asymptotic giant branch (AGB) stars. Results. Reconstructed images and visibilities show a strong wavelength dependence. The H-band images display two bright spots whose positions are confirmed by the geometrical toy model. The inhomogeneities are qualitatively predicted by 3D simulations. At about 2,00 micron and in the region 2,35 - 2,50 micron, the photosphere appears extended and the radius is larger than in the H band. In this spectral region, the geometrical toy model locates a third bright spot outside the photosphere that can be a feature of the molecular layers. The wavelength dependence of the visibility can be qualitatively explained by 1D dynamical models of Mira atmospheres. The best-fitting photospheric models show a good match with the observed visibilities and give a photospheric diameter of theta = 8,82+-0,50 mas. The H2O molecule seems to be the dominant absorber in the molecular layers. Conclusions. We show that the atmosphere of VX Sgr rather resembles Mira/AGB star model atmospheres than RSG model atmospheres. In particular, we see molecular (water) layers that are typical for Mira stars.Comment: 9 Pages, Accepted for publication on Astronomy & Astrophysics, two references update

Late-Type Red Supergiants: Too Cool for the Magellanic Clouds?

We have identified seven red supergiants (RSGs) in the Large Magellanic Cloud (LMC) and four RSGs in the Small Magellanic Cloud (SMC), all of which have spectral types that are considerably later than the average type observed in their parent galaxy. Using moderate-resolution optical spectrophotometry and the MARCS stellar atmosphere models, we determine their physical properties and place them on the H-R diagram for comparison with the predictions of current stellar evolutionary tracks. The radial velocities of these stars suggest that they are likely all members of the Clouds rather than foreground dwarfs or halo giants. Their locations in the H-R diagram also show us that those stars are cooler than the current evolutionary tracks allow, appearing to the right of the Hayashi limit, a region in which stars are no longer in hydrodynamic equilibrium. These stars exhibit considerable variability in their V magnitudes, and three of these stars also show changes in their effective temperatures (and spectral types) on the time-scales of months. One of these stars, [M2002] SMC 055188, was caught in an M4.5 I state, as late as that seen in HV 11423 at its recent extreme: considerable later, and cooler, than any other supergiant in the SMC. In addition, we find evidence of variable extinction due to circumstellar dust and changes in the stars' luminosities, also consistent with our recent findings for HV 11423 - when these stars are hotter they are also dustier and more luminous. We suggest that these stars have unusual properties because they are in an unstable (and short-lived) evolutionary phase.Comment: 26 pages, 6 figures; submitted to the Astrophysical Journa

Modeling active electrolocation in weakly electric fish

In this paper, we provide a mathematical model for the electrolocation in weakly electric fishes. We first investigate the forward complex conductivity problem and derive the approximate boundary conditions on the skin of the fish. Then we provide a dipole approximation for small targets away from the fish. Based on this approximation, we obtain a non-iterative location search algorithm using multi-frequency measurements. We present numerical experiments to illustrate the performance and the stability of the proposed multi-frequency location search algorithm. Finally, in the case of disk- and ellipse-shaped targets, we provide a method to reconstruct separately the conductivity, the permittivity, and the size of the targets from multi-frequency measurements.Comment: 37 pages, 11 figure

It's a wonderful tail: the mass loss history of Mira

Recent observations of the Mira AB binary system have revealed a surrounding arc-like structure and a stream of material stretching 2 degrees away in opposition to the arc. The alignment of the proper motion vector and the arc-like structure shows the structures to be a bow shock and accompanying tail. We have successfully hydrodynamically modelled the bow shock and tail as the interaction between the asymptotic giant branch (AGB) wind launched from Mira A and the surrounding interstellar medium. Our simulations show that the wake behind the bow shock is turbulent: this forms periodic density variations in the tail similar to those observed. We investigate the possiblity of mass-loss variations, but find that these have limited effect on the tail structure. The tail is estimated to be approximately 450,000 years old, and is moving with a velocity close to that of Mira itself. We suggest that the duration of the high mass-loss phase on the AGB may have been underestimated. Finally, both the tail curvature and the rebrightening at large distance can be qualitatively understood if Mira recently entered the Local Bubble. This is estimated to have occured 17 pc downstream from its current location.Comment: 12 pages, 3 colour figures, accepted by ApJ Part II (Letters

Molecular Hydrogen Emission Lines in Far Ultraviolet Spectroscopic Explorer Observations of Mira B

We present new Far Ultraviolet Spectroscopic Explorer (FUSE) observations of Mira A's wind-accreting companion star, Mira B. We find that the strongest lines in the FUSE spectrum are H2 lines fluoresced by H I Lyman-alpha. A previously analyzed Hubble Space Telescope (HST) spectrum also shows numerous Lyman-alpha fluoresced H2 lines. The HST lines are all Lyman band lines, while the FUSE H2 lines are mostly Werner band lines, many of them never before identified in an astrophysical spectrum. We combine the FUSE and HST data to refine estimates of the physical properties of the emitting H2 gas. We find that the emission can be reproduced by an H2 layer with a temperature and column density of T=3900 K and log N(H2)=17.1, respectively. Another similarity between the HST and FUSE data, besides the prevalence of H2 emission, is the surprising weakness of the continuum and high temperature emission lines, suggesting that accretion onto Mira B has weakened dramatically. The UV fluxes observed by HST on 1999 August 2 were previously reported to be over an order of magnitude lower than those observed by HST and the International Ultraviolet Explorer (IUE) from 1979--1995. Analysis of the FUSE data reveals that Mira B was still in a similarly low state on 2001 November 22.Comment: 23 pages, 6 figures; AASTEX v5.0 plus EPSF extensions in mkfig.sty; accepted by Ap