Recent high-contrast imaging surveys, looking for planets in young, nearby systems showed evidence of a small number of giant planets at relatively large separation beyond typically 20 au where those surveys are the most sensitive. Access to smaller physical separations between 5 and 20 au is the next step for future planet imagers on 10 m telescopes and ELTs in order to bridge the gap with indirect techniques (radial velocity, transit, astrometry with Gaia). In that context, we recently proposed a new algorithm, Keplerian-Stacker, combining multiple observations acquired at different epochs and taking into account the orbital motion of a potential planet present in the images to boost the ultimate detection limit. We showed that this algorithm is able to find planets in time series of simulated images of SPHERE even when a planet remains undetected at one epoch. Here, we validate the K-Stacker algorithm performances on real SPHERE datasets, to demonstrate its resilience to instrumental speckles and the gain offered in terms of true detection. This will motivate future dedicated multi-epoch observation campaigns in high-contrast imaging to search for planets in emitted and reflected light. Results. We show that K-Stacker achieves high success rate when the SNR of the planet in the stacked image reaches 7. The improvement of the SNR ratio goes as the square root of the total exposure time. During the blind test and the redetection of HD 95086 b, and betaPic b, we highlight the ability of K-Stacker to find orbital solutions consistent with the ones derived by the state of the art MCMC orbital fitting techniques, confirming that in addition to the detection gain, K-Stacker offers the opportunity to characterize the most probable orbital solutions of the exoplanets recovered at low signal to noise.

Using three fiducial Nbody+SPH simulations, we follow the merging of two disk galaxies with a hot gaseous halo component each, and examine whether the merger remnant can be a spiral galaxy. The stellar progenitor disks are destroyed by violent relaxation during the merging and most of their stars form a classical bulge, while the remaining form a thick disk and its bar. A new stellar disk forms subsequently and gradually in the remnant from the gas accreted mainly from the halo. It is vertically thin and well extended in its equatorial plane. A bar starts forming before the disk is fully in place, contrary to what is assumed in idealised simulations of isolated bar-forming galaxies. It has morphological features such as ansae and boxy/peanut bulges. Stars of different ages populate different parts of the box/peanut. A disky pseudobulge forms also, so that by the end of the simulation, all three types of bulges coexist. The oldest stars are found in the classical bulge, followed by those of the thick disk, then by those in the thin disk. The youngest stars are in the spiral arms and the disky pseudobulge. The disk surface density profiles are of type II (exponential with downbending), and the circular velocity curves are flat and show that the disks are submaximum in these examples: two clearly so and one near-borderline between maximum and submaximum. On average, only roughly between 10 and 20% of the stellar mass is in the classical bulge of the final models, i.e. much less than in previous simulations.

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.

Context: The solution of the nonlocal thermodynamical equilibrium (non-LTE) radiative transfer equation usually relies on stationary iterative methods, which may falsely converge in some cases. Furthermore, these methods are often unable to handle large-scale systems, such as molecular spectra emerging from, for example, cool stellar atmospheres. Aims: Our objective is to develop a new method, which aims to circumvent these problems, using nonstationary numerical techniques and taking advantage of parallel computers. Methods: The technique we develop may be seen as a generalization of the coupled escape probability method. It solves the statistical equilibrium equations in all layers of a discretized model simultaneously. The numerical scheme adopted is based on the generalized minimum residual method. Result:. The code has already been applied to the special case of the water spectrum in a red supergiant stellar atmosphere. This demonstrates the fast convergence of this method, and opens the way to a wide variety of astrophysical problems.

The structures of the outer atmospheres of red giants are very complex. The notion of large optically thick molecular spheres around the stars (MOLspheres) has been invoked in order to explain e.g. spectro-interferometric observations. However, high-resolution spectra in the mid-IR do not easily fit into this picture. They rule out any large sphere of water vapour in LTE surrounding red giants. Our aim here is to investigate high-resolution, mid-infrared spectra for a range of red giants, from early-K to mid M. We have recorded 12 microns spectra of 10 well-studied bright red giants, with TEXES on the IRTF. We find that all giants in our study cooler than 4300 K, spanning a range of effective temperatures, show water absorption lines stronger than expected. The strengths of the lines vary smoothly with spectral type. We identify several spectral features in the wavelength region that undoubtedly are formed in the photosphere. From a study of water-line ratios of the stars, we find that the excitation temperatures, in the line-forming regions, are several hundred Kelvin lower than expected from a classical photospheric model. This could either be due to an actually lower temperature structure in the outer regions of the photospheres caused by, for example, extra cooling, or due to non-LTE level populations, affecting the source function and line opacities. We have demonstrated that these diagnostically interesting water lines are a general feature of red giants across spectral types, and we argue for a general explanation of their formation rather than explanations requiring specific properties. Since the water lines are neither weak (filled in by emission) nor appear in emission, as predicted by LTE MOLsphere models in their simplest forms, the evidence for the existence of such large optically-thick, molecular spheres enshrouding the stars is weakened. (abbreviated)

The physical structures of the outer atmospheres of red giants are not known. They are certainly complex and a range of recent observations are showing that we need to embrace to non-classical atmosphere models to interpret these regions. This region's properties is of importance, not the least, for the understanding of the mass-loss mechanism for these stars, which is not still understood. Here, we present observational constraints of the outer regions of red giants, based on mid-IR, high spectral resolution spectra. We also discuss possible non-LTE effects and highlight a new non-LTE code that will be used to analyse the spectra of these atmospheric layers. We conclude by mentioning our new SOFIA/EXES observations of red giants at 6 microns, where the vibration-rotation lines of water vapour can be detected and spectrally resolved for the first time.

We aim at addressing the questions related to galaxy mass assembly through major and minor wet merging processes in the redshift range 1<z<2. A consequent fraction of Milky Way like galaxies are thought to have undergone an unstable clumpy phase at this early stage. Using the adaptive mesh refinement code RAMSES, with a recent physically-motivated implementation of stellar feedback, we build the Merging and Isolated high-Redshift Adaptive mesh refinement Galaxies (MIRAGE) sample. It is composed of 20 mergers and 3 isolated idealized disks simulations with global physical properties in accordance with the 1<z<2 mass complete sample MASSIV. The numerical hydrodynamical resolution reaches 7 parsecs in the smallest Eulerian cells. Our simulations include: star formation, metal line cooling, metallicity advection, and a recent implementation of stellar feedback which encompasses OB-type stars radiative pressure, photo-ionization heating, and supernovae. The initial conditions are set to match the z~2 observations, thanks to a new public code DICE. The numerical resolution allows us to follow the formation and evolution of giant clumps formed in-situ from Jeans instabilities triggered by high initial gas fraction. The star formation history of isolated disks shows stochastic star formation rate, which proceeds from the complex behavior of the giant clumps. Our minor and major gas-rich merger simulations do not trigger starbursts, suggesting a saturation of the star formation in a turbulent and clumpy interstellar medium fed by substantial accretion from the circum-galactic medium. Our simulations are close to the normal regime of the disk-like star formation on a Schmidt-Kennicutt diagram. The mass-size relation and its rate of evolution matches observations, suggesting that the inside-out growth mechanisms of the stellar disk do not necessarily require to be achieved through a cold accretion.

The interpretation of water lines in red supergiant stellar atmospheres has been much debated over the past decade. The introduction of the so-called MOLspheres to account for near-infrared "extra" absorption has been controversial. We propose that non-LTE effects should be taken into account before considering any extra-photospheric contribution. After a brief introduction on the radiative transfer treatment and the inadequacy of classical treatments in the case of large-scale systems such as molecules, we present a new code, based on preconditioned Krylov subspace methods. Preliminary results suggest that NLTE effects lead to deeper water bands, as well as extra cooling.

In August 2002, the near-Earth asteroid 2002 NY40, made its closest approach to the Earth. This provided an opportunity to study a near-Earth asteroid with a variety of instruments. Several of the telescopes at the Maui Space Surveillance System were trained at the asteroid and collected adaptive optics images, photometry and spectroscopy. Analysis of the imagery reveals the asteroid is triangular shaped with significant self-shadowing. The photometry reveals a 20-hour period and the spectroscopy shows that the asteroid is a Q-type.

We study the effect of a central mass concentration (CMC) on the secular evolution of a barred disc galaxy. Unlike previous studies, we use fully self-consistent 3D $N$-body simulations with live haloes, which are known to be important for bar evolution. The CMC is introduced gradually, to avoid transients. In all cases where the mass of the CMC is of the order of, or more than, a few per cent of the mass of the disc, the strength of the bar decreases noticeably. The amount of this decrease depends strongly on the bar type. For the same CMC, bars with exponential surface-density profile, which formed in a disk-dominated galaxy (MD-type bars), can be totally destroyed, while strong bars with a flat surface-density profile, whose evolution is largely due to the halo (MH-type bars), witness only a decrease of their strength. This decrease occurs simultaneously from both the innermost and outermost parts of the bar. The CMC has a stronger effect on the Fourier components of higher azimuthal wave number $m$, leading to fatter and/or less rectangular bars. Furthermore, the CMC changes the side-on outline from peanut-shaped to boxy or, for massive CMCs, to elliptical. Similarly, side-on initially boxy outlines can be destroyed. The CMC also influences the velocity dispersion profiles. Most of the decrease of the bar strength occurs while the mass of the CMC increases and it is accompanied by an increase of the pattern speed. In all our simulations, the mass of the CMC necessary in order to destroy the bar is at least several per cent of the mass of the disc. This argues that observed super-massive black holes are not likely to destroy pre-existing bars.

In this paper we measure the two-body relaxation time from the angular deflection of test particles launched in a rigid configuration of field particles. We find that centrally concentrated configurations have relaxation times that can be shorter than those of the corresponding homogeneous distributions by an order of magnitude or more. For homogeneous distributions we confirm that the relaxation time is proportional to the number of particles. On the other hand centrally concentrated configurations have a much shallower dependence, particularly for small values of the softening. The relaxation time increases with the inter-particle velocities and with softening. The latter dependence is not very strong, of the order of a factor of two when the softening is increased by an order of magnitude. Finally we show that relaxation times are the same on GRAPE-3 and GRAPE-4, dedicated computer boards with limited and high precision respectively.

In N-body simulations the force calculated between particles representing a given mass distribution is usually softened, to diminish the effect of graininess. In this paper we study the effect of such a smoothing, with the aim of finding an optimal value of the softening parameter. As already shown by Merritt (1996), for too small a softening the estimates of the forces will be too noisy, while for too large a softening the force estimates are systematically misrepresented. In between there is an optimal softening, for which the forces in the configuration approach best the true forces. The value of this optimal softening depends both on the mass distribution and on the number of particles used to represent it. For higher number of particles the optimal softening is smaller. More concentrated mass distributions necessitate smaller softening, but the softened forces are never as good an approximation of the true forces as for not centrally concentrated configurations. We give good estimates of the optimal softening for homogeneous spheres, Plummer spheres, and Dehnen spheres. We also give a rough estimate of this quantity for other mass distributions, based on the harmonic mean distance to the $k$th neighbour ($k$ = 1, .., 12), the mean being taken over all particles in the configuration. Comparing homogeneous Ferrers ellipsoids of different shapes we show that the axial ratios do not influence the value of the optimal softening. Finally we compare two different types of softening, a spline softening (Hernquist & Katz 1989) and a generalisation of the standard Plummer softening to higher values of the exponent. We find that the spline softening fares roughly as well as the higher powers of the power-law softening and both give a better representation of the forces than the standard Plummer softening.

We report on the results of the \hethet\ experiment at the underground accelerator facility LUNA (Gran Sasso). For the first time the lowest projectile energies utilized for the cross section measurement correspond to energies below the center of the solar Gamow peak ($E_{\rm 0}$=22 keV). The data provide no evidence for the existence of a hypothetical resonance in the energy range investigated. Although no extrapolation is needed anymore (except for energies at the low-energy tail of the Gamow peak), the data must be corrected for the effects of electron screening, clearly observed the first time for the \hethet\ reaction. The effects are however larger than expected and not understood, leading presently to the largest uncertainty on the quoted $S_{\rm b}(E_{\rm 0})$ value for bare nuclides (=5.40 MeV b).

The performance and accuracy of a GRAPE-3 system for collisionless N-body simulations is discussed. After a description of the hardware configurations available to us at Marseille, and the usefulness of on-line analysis, we concentrate on the performance and accuracy of the direct summation and tree code software. For the former we discuss the sources of round-off errors. The standard Barnes-Hut tree code has been modified for GRAPE-3, by dividing particles into blocks and performing the tree traversal for the whole block. The forces are then calculated by direct summation over the whole interaction list. The performance of the tree code depends on the number of particles in the block, the optimum number depending on the speed of the front end and the number of boards. The time per step scales with the number of particles, and decreases only weakly as the tolerance increases. We discuss the effect of the front end and compare performances with those of general purpose computers. The accuracy of both direct summation and the tree code is discussed as function of the number of particles and the softening. The increased role of direct summation in the force calculation enhances the accuracy compared to the standard tree code. Finally, we follow the evolution of an isolated barred galaxy using different hardware and software in order to assess the reliability of our results. We find excellent agreement between the pattern speed of the bar in direct summation simulations run on the high precision GRAPE-4 machines and on our GRAPE-3 system. The agreement with the tree code is very good for tolerance values smaller than about 1.0. We conclude that GRAPE-3 systems are well suited for collisionless simulations, in particular those of galaxies.