Dense molecular globulettes and the dust arc towards the runaway O star AE Aur (HD 34078)

Some runaway stars are known to display IR arc-like structures around them, resulting from their interaction with surrounding interstellar material. The properties of these features as well as the processes involved in their formation are still poorly understood. We aim at understanding the physical mechanisms that shapes the dust arc observed near the runaway O star AEAur (HD34078). We obtained and analyzed a high spatial resolution map of the CO(1-0) emission that is centered on HD34078, and that combines data from both the IRAM interferometer and 30m single-dish antenna. The line of sight towards HD34078 intersects the outer part of one of the detected globulettes, which accounts for both the properties of diffuse UV light observed in the field and the numerous molecular absorption lines detected in HD34078's spectra, including those from highly excited H2 . Their modeled distance from the star is compatible with the fact that they lie on the 3D paraboloid which fits the arc detected in the 24 {\mu}m Spitzer image. Four other compact CO globulettes are detected in the mapped area. These globulettes have a high density and linewidth, and are strongly pressure-confined or transient. The good spatial correlation between the CO globulettes and the IR arc suggests that they result from the interaction of the radiation and wind emitted by HD 34078 with the ambient gas. However, the details of this interaction remain unclear. A wind mass loss rate significantly larger than the value inferred from UV lines is favored by the large IR arc size, but does not easily explain the low velocity of the CO globulettes. The effect of radiation pressure on dust grains also meets several issues in explaining the observations. Further observational and theoretical work is needed to fully elucidate the processes shaping the gas and dust in bow shocks around runaway O stars. (Abridged)Comment: Accepted for publication in Astronomy & Astrophysic

Interstellar H2 toward HD 37903

We present an analysis of interstellar H2 toward HD 37903, which is a hot, B 1.5 V star located in the NGC 2023 reflection nebula. Meyer et al. (2001) have used a rich spectrum of vibrationally excited H2 observed by the HST to calculate a model of the interstellar cloud toward HD 37903. We extend Mayer's analysis by including the v"=0 vibrational level observed by the FUSE satellite. The T01 temperature should not be interpreted as a rotational temperature, but rather as a temperature of thermal equilibrium between the ortho and para H2. The ortho to para H2 ratio is lower for collisionally populated levels than for the levels populated by fluorescence. The PDR model of the cloud located in front of HD 37903 points to a gas temperature Tkin=110-377 K, hydrogen density nH=1874-544 cm^-3 and the star-cloud distance of 0.45 pc

Low-ionization iron-rich Broad Absorption-Line Quasar SDSS J1652+2650: Physical conditions in the ejected gas from excited FeII and metastable HeI

We present high-resolution VLT/UVES spectroscopy and a detailed analysis of the unique Broad Absorption-Line system towards the quasar SDSS J165252.67+265001.96. This system exhibits low-ionization metal absorption lines from the ground states and excited energy levels of Fe II and Mn II, and the meta-stable 2^3S excited state of He I. The extended kinematics of the absorber encompasses three main clumps with velocity offsets of -5680, -4550, and -1770 km s$^{-1}$ from the quasar emission redshift, $z=0.3509\pm0.0003$, derived from [O II] emission. Each clump shows moderate partial covering of the background continuum source, $C_f \approx [0.53; 0.24; 0.81]$. We discuss the excitation mechanisms at play in the gas, which we use to constrain the distance of the clouds from the Active Galactic Nucleus (AGN) as well as the density, temperature, and typical sizes of the clouds. The number density is found to be $n_{\rm H} \sim 10^4\rm cm^{-3}$ and the temperature $T_e \sim 10^4\rm\,K$, with longitudinal cloudlet sizes of $\gtrsim0.01$ pc. Cloudy photo-ionization modelling of He I$^{*}$, which is also produced at the interface between the neutral and ionized phases, assuming the number densities derived from Fe II, constrains the ionization parameter to be $\log U \sim -3$. This corresponds to distances of a few 100 pc from the AGN. We discuss these results in the more general context of associated absorption-line systems and propose a connection between FeLoBALs and the recently-identified molecular-rich intrinsic absorbers. Studies of significant samples of FeLoBALs, even though rare per se, will soon be possible thanks to large dedicated surveys paired with high-resolution spectroscopic follow-ups.Comment: Accepted for publication in MNRAS, 27 pages, 21 Figure

The KX method for producing K-band flux-limited samples of quasars

The longstanding question of the extent to which the quasar population is affected by dust extinction, within host galaxies or galaxies along the line of sight, remains open. More generally, the spectral energy distributions of quasars vary significantly and flux-limited samples defined at different wavelengths include different quasars. Surveys employing flux measurements at widely separated wavelengths are necessary to characterise fully the spectral properties of the quasar population. The availability of panoramic near-infrared detectors on large telescopes provides the opportunity to undertake surveys capable of establishing the importance of extinction by dust on the observed population of quasars. We introduce an efficient method for selecting K-band, flux-limited samples of quasars, termed ``KX'' by analogy with the UVX method. This method exploits the difference between the power-law nature of quasar spectra and the convex spectra of stars: quasars are relatively brighter than stars at both short wavelengths (the UVX method) and long wavelengths (the KX method). We consider the feasibility of undertaking a large-area KX survey for damped Ly-alpha galaxies and gravitational lenses using the planned UKIRT wide-field near-infrared camera.Comment: 6 pages, 3 figures, to appear in MNRA

CO emission and variable CH and CH+ absorption towards HD34078: evidence for a nascent bow shock ?

The runaway star HD34078, initially selected to investigate small scale structure in a foreground diffuse cloud has been shown to be surrounded by highly excited H2. We first search for an association between the foreground cloud and HD34078. Second, we extend previous investigations of temporal absorption line variations (CH, CH+, H2) in order to better characterize them. We have mapped the CO(2-1) emission at 12 arcsec resolution around HD34078's position, using the 30 m IRAM antenna. The follow-up of CH and CH+ absorption lines has been extended over 5 more years. In parallel, CH absorption towards the reddened star Zeta Per have been monitored to check the homogeneity of our measurements. Three more FUSE spectra have been obtained to search for N(H2) variations. CO observations show a pronounced maximum near HD34078's position, clearly indicating that the star and diffuse cloud are associated. The optical spectra confirm the reality of strong, rapid and correlated CH and CH+ fluctuations. On the other hand, N(H2, J=0) has varied by less than 5 % over 4 years. We also discard N(CH) variations towards Zeta Per at scales less than 20 AU. Observational constraints from this work and from 24 micron dust emission appear to be consistent with H2 excitation but inconsistent with steady-state bow shock models and rather suggest that the shell of compressed gas surrounding HD34078, is seen at an early stage of the interaction. The CH and CH+ time variations as well as their large abundances are likely due to chemical structure in the shocked gas layer located at the stellar wind/ambient cloud interface. Finally, the lack of variations for both N(H2, J=0) towards HD34078 and N(CH) towards Zeta Per suggests that quiescent molecular gas is not subject to pronounced small-scale structure.Comment: 19 pages, 15 figures, accepted for publication in A&

Metallicities, dust and molecular content of a QSO-Damped Lyman-{\alpha} system reaching log N (H i) = 22: An analog to GRB-DLAs

We present the elemental abundance and H2 content measurements of a Damped Lyman-{\alpha} (DLA) system with an extremely large H i column density, log N(H i) (cm-2) = 22.0+/-0.10, at zabs = 3.287 towards the QSO SDSS J 081634+144612. We measure column densities of H2, C i, C i^*, Zn ii, Fe ii, Cr ii, Ni ii and Si ii from a high signal-to-noise and high spectral resolution VLT-UVES spectrum. The overall metallicity of the system is [Zn/H] = -1.10 +/- 0.10 relative to solar. Two molecular hydrogen absorption components are seen at z = 3.28667 and 3.28742 (a velocity separation of \approx 52 km s-1) in rotational levels up to J = 3. We derive a total H2 column density of log N(H2) (cm-2) = 18.66 and a mean molecular fraction of f = 2N(H2)/[2N(H2) + N(H i)] = 10-3.04+/-0.37, typical of known H2-bearing DLA systems. From the observed abundance ratios we conclude that dust is present in the Interstellar Medium (ISM) of this galaxy, with a enhanced abundance in the H2-bearing clouds. However, the total amount of dust along the line of sight is not large and does not produce any significant reddening of the background QSO. The physical conditions in the H2-bearing clouds are constrained directly from the column densities of H2 in different rotational levels, C i and C i^* . The kinetic temperature is found to be T = 75 K and the particle density lies in the range nH = 50-80 cm-3 . The neutral hydrogen column density of this DLA is similar to the mean H i column density of DLAs observed at the redshift of {\gamma}-ray bursts (GRBs). We explore the relationship between GRB-DLAs and high column density end of QSO-DLAs finding that the properties (metallicity and depletion) of DLAs with log N(H i) > 21.5 in the two populations do not appear to be significantly different

Interstellar absorptions towards the LMC: Small-scale density variations in Milky Way disc gas

Observations show that the ISM contains sub-structure on scales less than 1 pc, detected in the form of spatial and temporal variations in column densities or optical depth. Despite the number of detections, the nature and ubiquity of the small-scale structure in the ISM is not yet fully understood. We use UV absorption data mainly from FUSE and partly from STIS of six LMC stars (Sk-67{\deg}111, LH54-425, Sk-67{\deg}107, Sk-67{\deg}106, Sk-67{\deg}104, and Sk-67{\deg}101), all located within 5 arcmin of each other, and analyse the physical properties of the Galactic disc gas in front of the LMC on sub-pc scales. We analyse absorption lines of a number of ions within the UV spectral range. Most importantly, interstellar molecular hydrogen, neutral oxygen, and fine-structure levels of neutral carbon have been used in order to study changes in the density and the physical properties of the Galactic disc gas over small angular scales. While most species do not show any significant variation in their column densities, we find an enhancement of almost 2 dex for H_2 from Sk-67{\deg}111 to Sk-67{\deg}101, accompanied by only a small variation in the OI column density. Based on the formation-dissociation equilibrium, we trace these variations to the actual density variations in the molecular gas. On the smallest spatial scale of < 0.08 pc, between Sk-67{\deg}107 and LH54-425, we find a gas density variation of a factor of 1.8. The line of sight towards LH54-425 does not follow the relatively smooth change seen from Sk-67{\deg}101 to Sk-67{\deg}111, suggesting that sub-structure might exist on a smaller spatial scale than the linear extent of our sight-lines. Our observations suggest that the detected H_2 in these six lines of sight is not necessarily physically connected, but that we are sampling molecular cloudlets with pathlengths < 0.1-1.8 pc and possibly different densities.Comment: 17 pages, 8 figures. Accepted for publication in A&

Radiative equilibrium in Monte Carlo radiative transfer using frequency distribution adjustment

The Monte Carlo method is a powerful tool for performing radiative equilibrium calculations, even in complex geometries. The main drawback of the standard Monte Carlo radiative equilibrium methods is that they require iteration, which makes them numerically very demanding. Bjorkman & Wood recently proposed a frequency distribution adjustment scheme, which allows radiative equilibrium Monte Carlo calculations to be performed without iteration, by choosing the frequency of each re-emitted photon such that it corrects for the incorrect spectrum of the previously re-emitted photons. Although the method appears to yield correct results, we argue that its theoretical basis is not completely transparent, and that it is not completely clear whether this technique is an exact rigorous method, or whether it is just a good and convenient approximation. We critically study the general problem of how an already sampled distribution can be adjusted to a new distribution by adding data points sampled from an adjustment distribution. We show that this adjustment is not always possible, and that it depends on the shape of the original and desired distributions, as well as on the relative number of data points that can be added. Applying this theorem to radiative equilibrium Monte Carlo calculations, we provide a firm theoretical basis for the frequency distribution adjustment method of Bjorkman & Wood, and we demonstrate that this method provides the correct frequency distribution through the additional requirement of radiative equilibrium. We discuss the advantages and limitations of this approach, and show that it can easily be combined with the presence of additional heating sources and the concept of photon weighting. However, the method may fail if small dust grains are included... (abridged)Comment: 17 pages, 2 figures, accepted for publication in New Astronom

A far UV study of interstellar gas towards HD34078: high excitation H2 and small scale structure - Based on observations performed by the FUSE mission and at the CFHT telescope

To investigate the presence of small scale structure in the spatial distribution of H2 molecules we have undertaken repeated FUSE UV observations of the runaway O9.5V star, HD34078. In this paper we present five spectra obtained between January 2000 and October 2002. These observations reveal an unexpectedly large amount of highly excited H2. Column densities for H2 levels from (v = 0, J = 0) up to (v = 0, J = 11) and for several v = 1 and v = 2 levels are determined. These results are interpreted in the frame of a model involving essentially two components: i) a foreground cloud (unaffected by HD34078) responsible for the H2 (J = 0, 1), CI, CH, CH+ and CO absorptions; ii) a dense layer of gas (n = 10E4 cm-3) close to the O star and strongly illuminated by its UV flux which accounts for the presence of highly excited H2. Our model successfully reproduces the H2 excitation, the CI fine-structure level populations as well as the CH, CH+ and CO column densities. We also examine the time variability of H2 absorption lines tracing each of these two components. From the stability of the J = 0, 1 and 2 damped H2 profiles we infer a 3 sigma upper limit on column density variations Delta(N(H2))/N(H2) of 5% over scales ranging from 5 to 50 AU. This result clearly rules out any pronounced ubiquitous small scale "density" structure of the kind apparently seen in HI. The lines from highly excited gas are also quite stable (equivalent to Delta(N)/N <= 30%) indicating i) that the ambient gas through which HD34078 is moving is relatively uniform and ii) that the gas flow along the shocked layer is not subject to marked instabilitie

Molecular hydrogen in high-redshift Damped Lyman-alpha systems: The VLT/UVES database

We present the current status of ongoing searches for molecular hydrogen in high-redshift (1.8 < zabs <= 4.2) Damped Lyman-alpha systems (DLAs) capitalising on observations performed with the ESO Very Large Telescope (VLT) Ultraviolet and Visual Echelle Spectrograph (UVES). We identify 77 DLAs/strong sub-DLAs, with log N(HI) >= 20 and z_abs > 1.8, which have data that include redshifted H2 Lyman and/or Werner-band absorption lines. This sample of HI, H2 and metal line measurements, performed in an homogeneous manner, is more than twice as large as our previous sample (Ledoux et al. 2003) and considers every system in which searches for H2 could be completed so far, including all non-detections. H2 is detected in thirteen of the systems with molecular fractions of values between f=5x10^-7 and f=0.1, where f=2N(H2)/(2N(H2)+N(HI)). Upper limits are measured for the remaining 64 systems with detection limits of typically log N(H2)=14.3, corresponding to log f<-5. We find that about 35% of the DLAs with metallicities relative to solar [X/H]>=-1.3 (i.e., 1/20th solar), with X = Zn, S or Si, have molecular fractions log f>-4.5, while H2 is detected -- regardless of the molecular fraction -- in 50% of them. In contrast, only about 4% of the [X/H]-4.5. We show that the presence of H2 does not strongly depend on the total neutral hydrogen column density, although the probability of finding log f>-4.5 is higher for log N(HI)>=20.8 than below this limit (19% and 7% respectively). The overall H2 detection rate in log N(HI)>=20 DLAs is found to be about 16% (10% considering only log f>-4.5 detections) after correction for a slight bias towards large N(HI). [truncated]Comment: 11 pages, 1 table, 10 figures. Accepted for publication in A&