-
Radial Distributions of Coronal Electron Temperatures: specificities of the DYN model
Authors:
Joseph F. Lemaire,
Athanassios C. Katsiyannis
Abstract:
This paper is a follow up of the article where Lemaire and Stegen (2016) introduced their DYN method to calculate coronal temperature profiles from given radial distributions of the coronal and solar wind (SW) electron densities. Several such temperature profiles are calculated and presented corresponding to a set of given empirical density models derived from eclipse observations and in-situ meas…
▽ More
This paper is a follow up of the article where Lemaire and Stegen (2016) introduced their DYN method to calculate coronal temperature profiles from given radial distributions of the coronal and solar wind (SW) electron densities. Several such temperature profiles are calculated and presented corresponding to a set of given empirical density models derived from eclipse observations and in-situ measurements of the electron density and bulk velocity at 1 AU. The DYN temperature profiles obtained for the equatorial and polar regions of the corona challenge the results deduced since 1958 from singular hydrodynamical models of the SW. In these models - where the expansion velocity transits through a singular saddle point - the maximum coronal temperature is predicted to be located at the base of the corona, while in all DYN models the altitude of the maximum temperature is found at significantly higher altitudes in the mid-corona. Furthermore, the maximum of the DYN-estimated temperatures is found at much higher altitudes over the polar regions and coronal holes, than over the equator. However, at low altitudes, in the inner corona, the DYN temperatures are always smaller at high latitudes, than at low equatorial latitudes. This appears well in agreement with existing coronal hole observations. These findings have serious implications on the open questions: what is the actual source of the coronal heating, and where is the maximum energy deposited within the solar corona?
△ Less
Submitted 25 February, 2021; v1 submitted 18 February, 2020;
originally announced February 2020.
-
X-ray induced deuterium enrichment on N-rich organics in protoplanetary disks: an experimental investigation using synchrotron light
Authors:
Lisseth Gavilan,
Laurent Remusat,
Mathieu Roskosz,
Horia Popescu,
Nicolas Jaouen,
Christophe Sandt,
Cornelia Jager,
Thomas Henning,
Alexandre Simionovici,
Jean Louis Lemaire,
Denis Mangin,
Nathalie Carrasco
Abstract:
The deuterium enrichment of organics in the interstellar medium, protoplanetary disks and meteorites has been proposed to be the result of ionizing radiation. The goal of this study is to quantify the effects of soft X-rays (0.1 - 2 keV), a component of stellar radiation fields illuminating protoplanetary disks, on the refractory organics present in the disks. We prepared tholins, nitrogen-rich co…
▽ More
The deuterium enrichment of organics in the interstellar medium, protoplanetary disks and meteorites has been proposed to be the result of ionizing radiation. The goal of this study is to quantify the effects of soft X-rays (0.1 - 2 keV), a component of stellar radiation fields illuminating protoplanetary disks, on the refractory organics present in the disks. We prepared tholins, nitrogen-rich complex organics, via plasma deposition and used synchrotron radiation to simulate X-ray fluences in protoplanetary disks. Controlled irradiation experiments at 0.5 and 1.3 keV were performed at the SEXTANTS beam line of the SOLEIL synchrotron, and were followed by ex-situ infrared, Raman and isotopic diagnostics. Infrared spectroscopy revealed the loss of singly-bonded groups (N-H, C-H and R-N$\equiv$C) and the formation of sp$^3$ carbon defects. Raman analysis revealed the introduction of defects and structural amorphization. Finally, tholins were measured via secondary ion mass spectrometry (SIMS), revealing that significant D-enrichment is induced by X-ray irradiation. Our results are compared to previous experimental studies involving the thermal degradation and electron irradiation of organics. The penetration depth of soft X-rays in $μ$m-sized tholins leads to volume rather than surface modifications: lower energy X-rays (0.5 keV) induce a larger D-enrichment than 1.3 keV X-rays, reaching a plateau for doses larger than 5 $\times$ 10$^{27}$ eV cm$^{-3}$. Our work provides experimental evidence of a new non-thermal pathway to deuterium fractionation of organic matter.
△ Less
Submitted 5 May, 2017;
originally announced May 2017.
-
Hard X-ray irradiation of cosmic silicate analogs: structural evolution and astrophysical implications
Authors:
L. Gavilan,
C. Jäger,
A. Simionovici,
J. L. Lemaire,
T. Sabri,
E. Foy,
S. Yagoubi,
T. Henning,
D. Salomon,
G. Martinez-Criado
Abstract:
Protoplanetary disks, interstellar clouds, and active galactic nuclei, contain X-ray dominated regions. X-rays interact with the dust and gas present in such environments. While a few laboratory X-ray irradiation experiments have been performed on ices, X-ray irradiation experiments on bare cosmic dust analogs have been scarce up to now. Our goal is to study the effects of hard X-rays on cosmic du…
▽ More
Protoplanetary disks, interstellar clouds, and active galactic nuclei, contain X-ray dominated regions. X-rays interact with the dust and gas present in such environments. While a few laboratory X-ray irradiation experiments have been performed on ices, X-ray irradiation experiments on bare cosmic dust analogs have been scarce up to now. Our goal is to study the effects of hard X-rays on cosmic dust analogs via in-situ X-ray diffraction. By using a hard X-ray synchrotron nanobeam, we seek to simulate cumulative X-ray exposure on dust grains during their lifetime in these astrophysical environments, and provide an upper limit on the effect of hard X-rays on dust grain structure.
We prepared enstatite nanograins, analogs to cosmic silicates, via the melting-quenching technique. These amorphous grains were then annealed to obtain polycrystalline grains. These were characterized via scanning electron microscopy and high-resolution transmission electron microscopy before irradiation. Powder samples were prepared in X-ray transparent substrates and were irradiated with hard X-rays nanobeams (29.4 keV) provided by beamline ID16B of the European Synchrotron Radiation Facility. X-ray diffraction images were recorded in transmission mode. We detected the amorphization of polycrystalline silicates embedded in an organic matrix after an accumulated X-ray exposure of 6.4 x 10$^{27}$ eV cm$^{-2}$. Pure crystalline silicate grains (without resin) did not exhibit amorphization. None of the amorphous silicate samples (pure and embedded in resin) underwent crystallization. We analyzed the evolution of the polycrystalline sample embedded in an organic matrix as a function of X-ray exposure. Loss of diffraction peak intensity, peak broadening, and the disappearance of discrete spots and arcs, revealed the amorphization of the resin embedded (originally polycrystalline) silicate sample.
△ Less
Submitted 20 January, 2016;
originally announced January 2016.
-
Formation of hydroxylamine on dust grains via ammonia oxidation
Authors:
Jiao He,
Gianfranco Vidali,
Jean-Louis Lemaire,
Robin T. Garrod
Abstract:
The quest to detect prebiotic molecules in space, notably amino acids, requires an understanding of the chemistry involving nitrogen atoms. Hydroxylamine (NH$_2$OH) is considered a precursor to the amino acid glycine. Although not yet detected, NH$_2$OH is considered a likely target of detection with ALMA. We report on an experimental investigation of the formation of hydroxylamine on an amorphous…
▽ More
The quest to detect prebiotic molecules in space, notably amino acids, requires an understanding of the chemistry involving nitrogen atoms. Hydroxylamine (NH$_2$OH) is considered a precursor to the amino acid glycine. Although not yet detected, NH$_2$OH is considered a likely target of detection with ALMA. We report on an experimental investigation of the formation of hydroxylamine on an amorphous silicate surface via the oxidation of ammonia. The experimental data are then fed into a simulation of the formation of NH$_2$OH in dense cloud conditions. On ices at 14 K and with a modest activation energy barrier, NH$_2$OH is found to be formed with an abundance that never falls below a factor 10 with respect to NH$_3$. Suggestions of conditions for future observations are provided.
△ Less
Submitted 19 January, 2015;
originally announced January 2015.
-
Collisionless transport equations derived from a kinetic exospheric solar wind model with kappa velocity distribution functions
Authors:
Gabriel Voitcu,
Viviane Pierrard,
Herve Lamy,
Joseph Lemaire,
Marius Echim
Abstract:
In this paper we discuss the collisionless transport equations, continuity, momentum and energy conservation, derived from a kinetic exospheric model of the solar wind based on a kappa velocity distribution function of the electrons. The model is stationary and is based on a non-monotonic potential energy for the protons. The present study is carried out for an exobase located at 1.5 solar radii a…
▽ More
In this paper we discuss the collisionless transport equations, continuity, momentum and energy conservation, derived from a kinetic exospheric model of the solar wind based on a kappa velocity distribution function of the electrons. The model is stationary and is based on a non-monotonic potential energy for the protons. The present study is carried out for an exobase located at 1.5 solar radii and for two different values of the kappa index. The maximum radial distance considered is equal to one astronomical unit. The moments of the velocity distribution function computed with the kinetic exospheric model for both electrons and protons are introduced into the mass continuity equation, momentum conservation equation and energy conservation equation. The relative importance of various terms in the macroscopic transport equations for each component species are analyzed and discussed. The results obtained show that the kinetic description based on kappa velocity distribution functions satisfies rigorously the transport equations that give a macroscopic description of the solar wind plasma.
△ Less
Submitted 14 October, 2014;
originally announced October 2014.
-
A review on solar wind modeling: kinetic and fluid aspects
Authors:
Marius Echim,
Joseph Lemaire,
Oystein Lie-Svendsen
Abstract:
We review the main advantages and limitations of the kinetic exospheric and fluid models of the solar wind (SW). We discuss the hydrostatic model imagined by Chapman, the first supersonic hydrodynamic models published by Parker and the first generation subsonic kinetic model proposed by Chamberlain. It is shown that a correct estimation of the electric field as in the second generation kinetic exo…
▽ More
We review the main advantages and limitations of the kinetic exospheric and fluid models of the solar wind (SW). We discuss the hydrostatic model imagined by Chapman, the first supersonic hydrodynamic models published by Parker and the first generation subsonic kinetic model proposed by Chamberlain. It is shown that a correct estimation of the electric field as in the second generation kinetic exospheric models developed by Lemaire and Scherer, provides a supersonic expansion of the corona, reconciling the hydrodynamic and the kinetic approach. The third generation kinetic exospheric models considers kappa velocity distribution function (VDF) instead of a Maxwellian at the exobase and in addition they treat a non-monotonic variation of the electric potential with the radial distance; the fourth generation exospheric models include Coulomb collisions based on the Fokker--Planck collision term. Multi-fluid models of the solar wind provide a coarse grained description and reproduce with success the spatio-temporal variation of SW macroscopic properties (density, bulk velocity). The main categories of multi-fluid SW models are reviewed: the 5-moment, or Euler, models, originally proposed by Parker to describe the supersonic SW expansion; the 8-moment and 16-moment fluid models, the gyrotropic approach with improved collision terms as well as the gyrotropic models based on observed VDFs. The outstanding problem of collisions, including the long range Coulomb encounters, is also discussed, both in the kinetic and multi-fluid context. Although for decades the two approaches have been seen as opposed, in this paper we emphasize their complementary. The review of the kinetic and fluid models of the solar wind contributes also to a better evaluation of the open questions still existent in SW modeling and suggests possible future developments.
△ Less
Submitted 4 June, 2013;
originally announced June 2013.
-
What happens when the geomagnetic field reverses?
Authors:
Joseph F. Lemaire,
Fred Singer
Abstract:
During geomagnetic field reversals the radiation belt high-energy proton populations become depleted. Their energy spectra become softer, with the trapped particles of highest energies being lost first, and eventually recovering after a field reversal. The radiation belts rebuild in a dynamical way with the energy spectra flattening on the average during the course of many millennia, but without e…
▽ More
During geomagnetic field reversals the radiation belt high-energy proton populations become depleted. Their energy spectra become softer, with the trapped particles of highest energies being lost first, and eventually recovering after a field reversal. The radiation belts rebuild in a dynamical way with the energy spectra flattening on the average during the course of many millennia, but without ever reaching complete steady state equilibrium between successive geomagnetic storm events determined by southward turnings of the IMF orientation. Considering that the entry of galactic cosmic rays and the solar energetic particles with energies above a given threshold are strongly controlled by the intensity of the northward component of the interplanetary magnetic field, we speculate that at earlier epochs when the geomagnetic dipole was reversed, the entry of these energetic particles into the geomagnetic field was facilitated when the interplanetary magnetic field was directed northward. Unlike in other complementary work where intensive numerical simulations have been used, our demonstration is based on a simple analytical extension of Stormer's theory. The access of GCR and SEP beyond geomagnetic cut-off latitudes is enhanced during epochs when the Earth's magnetic dipole is reduced, as already demonstrated earlier.
△ Less
Submitted 21 August, 2012;
originally announced August 2012.
-
The effect of a southward interplanetary magnetic field on Störmer's allowed regions
Authors:
Joseph F. Lemaire
Abstract:
The motion of a charged particle in a magnetic dipole has first been studied by Stormer. The different applications of Stormer's theory to aurorae, cosmic rays and Van Allen radiation belt particles are recalled in an historical perspective. In this paper, we expand the Stormer theory in order to take into account the effects produced by an additional uniform and stationary interplanetary magnetic…
▽ More
The motion of a charged particle in a magnetic dipole has first been studied by Stormer. The different applications of Stormer's theory to aurorae, cosmic rays and Van Allen radiation belt particles are recalled in an historical perspective. In this paper, we expand the Stormer theory in order to take into account the effects produced by an additional uniform and stationary interplanetary magnetic field (IMF) whose orientation is parallel or antiparallel to the magnetic moment of the dipole. A new expression is derived for the Stormer potential taking into account the additional IMF component. It is shown how Stormer's allowed and forbidden zones are influenced by the implementation of a northward or a southward IMF, and how a southward turning of the IMF orientation makes it easier for Solar Energetic Particle and Galactic Cosmic Rays to enter into the inner part of the geomagnetic field along interconnected magnetic field lines.
△ Less
Submitted 21 July, 2012;
originally announced July 2012.
-
Coronal temperature profiles obtained from kinetic models and from coronal brightness measurements obtained during solar eclipses
Authors:
V. Pierrard,
K. Borremans,
J. F. Lemaire
Abstract:
Coronal density, temperature and heat flux distributions for the equatorial and polar corona have been deduced by Lemaire [2012] from Saito's model of averaged coronal white light (WL) brightness and polarization observations. They are compared with those determined from a kinetic collisionless/exospheric model of the solar corona. This comparison indicates rather similar distributions at large ra…
▽ More
Coronal density, temperature and heat flux distributions for the equatorial and polar corona have been deduced by Lemaire [2012] from Saito's model of averaged coronal white light (WL) brightness and polarization observations. They are compared with those determined from a kinetic collisionless/exospheric model of the solar corona. This comparison indicates rather similar distributions at large radial distances (> 7 Rs) in the collisionless region. However, rather important differences are found close to the Sun in the acceleration region of the solar wind. The exospheric heat flux is directed away from the Sun, while that inferred from all WL coronal observations is in the opposite direction, i.e., conducting heat from the inner corona toward the chromosphere. This could indicate that the source of coronal heating rate extends up into the inner corona where it maximizes at r > 1.5 Rs well above the transition region.
△ Less
Submitted 20 July, 2012;
originally announced July 2012.
-
Convective Instability Of The Solar Corona: Why The Solar Wind Blows
Authors:
Joseph Lemaire
Abstract:
Chapman's (1957) conductive model of the solar corona is characterized by a temperature varying as r**(-2/7) with heliocentric distance r. The density distribution in this non-isothermal hydrostatic model has a minimum value at 123 RS, and increases with r above that altitude. It is shown that this hydrostatic model becomes convectively unstable above r = 35 RS, where the temperature lapse rate be…
▽ More
Chapman's (1957) conductive model of the solar corona is characterized by a temperature varying as r**(-2/7) with heliocentric distance r. The density distribution in this non-isothermal hydrostatic model has a minimum value at 123 RS, and increases with r above that altitude. It is shown that this hydrostatic model becomes convectively unstable above r = 35 RS, where the temperature lapse rate becomes superadiabatic. Beyond this radial distance heat conduction fails to be efficient enough to keep the temperature gradient smaller than the adiabatic lapse rate. We report the results obtained by Lemaire (1968) who showed that an additional mechanism is then required to transport the energy flux away from the Sun into interplanetary space. He pointed out that this additional mechanism is advection: i.e. the stationary hydrodynamic expansion of the corona. In other words the corona is unable to stay in hydrostatic equilibrium. The hydrodynamic solar wind expansion is thus a physical consequence of the too steep (superadiabatic) temperature gradient beyond the peak of coronal temperature that can be determined from white light brightness distributions observed during solar eclipses. The thermodynamic argument for the existence of a continuous solar wind expansion which is presented here, complements Parker's classical argument based on boundary conditions imposed to the solutions of the hydrodynamic equations for the coronal expansion: i.e. the inability of the mechanical forces to hold the corona in hydrostatic equilibrium. The thermodynamic argument presented here is based on the energy transport equation. It relies on the temperature distribution which becomes super-adiabatic above a certain altitude in the inner corona.
△ Less
Submitted 24 May, 2012;
originally announced May 2012.
-
Half a Century of Kinetic Solar Wind Models
Authors:
Joseph Lemaire
Abstract:
I outline the development of four generations of kinetic models, starting with Chamberlain's solar breeze exospheric model. It is shown why this first kinetic model did not give apposite supersonic evaporation velocities, like early hydrodynamic models of the solar wind. When a self-consistent polarization electric potential distribution is used in the coronal plasma, instead of the Pannekoek-Ross…
▽ More
I outline the development of four generations of kinetic models, starting with Chamberlain's solar breeze exospheric model. It is shown why this first kinetic model did not give apposite supersonic evaporation velocities, like early hydrodynamic models of the solar wind. When a self-consistent polarization electric potential distribution is used in the coronal plasma, instead of the Pannekoek-Rosseland's one, supersonic bulk velocities are readily obtained in the second generation of kinetic models. It is outlined how the third and fourth generations of these models have improved the agreement with observations of slow and fast speed solar wind streams.
△ Less
Submitted 24 May, 2012;
originally announced May 2012.
-
A 3D view of the outflow in the Orion Molecular Cloud 1 (OMC-1)
Authors:
H. D. Nissen,
N. J. Cunningham,
M. Gustafsson,
J. Bally,
J. -L. Lemaire,
C. Favre,
D. Field
Abstract:
The fast outflow emerging from a region associated with massive star formation in the Orion Molecular Cloud 1 (OMC-1), located behind the Orion Nebula, appears to have been set in motion by an explosive event. Here we study the structure and dynamics of outflows in OMC-1. We combine radial velocity and proper motion data for near-IR emission of molecular hydrogen to obtain the first 3-dimensional…
▽ More
The fast outflow emerging from a region associated with massive star formation in the Orion Molecular Cloud 1 (OMC-1), located behind the Orion Nebula, appears to have been set in motion by an explosive event. Here we study the structure and dynamics of outflows in OMC-1. We combine radial velocity and proper motion data for near-IR emission of molecular hydrogen to obtain the first 3-dimensional (3D) structure of the OMC-1 outflow. Our work illustrates a new diagnostic tool for studies of star formation that will be exploited in the near future with the advent of high spatial resolution spectro-imaging in particular with data from the Atacama Large Millimeter Array (ALMA). We use published radial and proper motion velocities obtained from the shock-excited vibrational emission in the H2 v=1-0 S(1) line at 2.122 $μ$m obtained with the GriF instrument on the Canada-France-Hawaii Telescope, the Apache Point Observatory, the Anglo-Australian Observatory and the Subaru Telescope. These data give the 3D velocity of ejecta yielding a 3D reconstruction of the outflows. This allows one to view the material from different vantage points in space giving considerable insight into the geometry. Our analysis indicates that the ejection occurred <720 years ago from a distorted ring-like structure of ~15" (6000 AU) in diameter centered on the proposed point of close encounter of the stars BN, source I and maybe also source n. We propose a simple model involving curvature of shock trajectories in magnetic fields through which the origin of the explosion and the centre defined by extrapolated proper motions of BN, I and n may be brought into spatial coincidence.
△ Less
Submitted 14 March, 2012;
originally announced March 2012.
-
Sticking coefficient of hydrogen and deuterium on silicates under interstellar conditions
Authors:
H. Chaabouni,
H. Bergeron,
S. Baouche,
F. Dulieu,
E. Matar,
E. Congiu,
L. Gavilan,
J. L. Lemaire
Abstract:
Sticking of H and D atoms on interstellar dust grains is the first step in molecular hydrogen formation, which is a key reaction in the InterStellar Medium (ISM). After studying the sticking coefficients of H2 and D2 molecules on amorphous silicate surfaces experimentally and theoretically, we extrapolate the results to the sticking coefficient of atoms and propose a formulae that gives the sticki…
▽ More
Sticking of H and D atoms on interstellar dust grains is the first step in molecular hydrogen formation, which is a key reaction in the InterStellar Medium (ISM). After studying the sticking coefficients of H2 and D2 molecules on amorphous silicate surfaces experimentally and theoretically, we extrapolate the results to the sticking coefficient of atoms and propose a formulae that gives the sticking coefficients of H and D on both silicates and icy dust grains. In our experiments, we used the King and Wells method for measuring the sticking coefficients of H2 and D2 molecules on a silicate surface held at 10 K. It consists of measuring with a QMS (quadrupole mass spectrometer) the signals of H2 and D2 molecules reflected by the surface during the exposure of the sample to the molecular beam at a temperature ranging from 20 K to 340 K. We tested the efficiency of a physical model, developed previously for sticking on water-ice surfaces. We applied this model to our experimental results for the sticking coefficients of H2 and D2 molecules on a silicate surface and estimated the sticking coefficient of atoms by a single measurement of atomic recombination and propose an extrapolation. Sticking of H, D, HD, H2, and D2 on silicates grains behaves the same as on icy dust grains. The sticking decreases with the gas temperature, and is dependent on the mass of the impactor. The sticking coefficient for both surfaces and impactors can be modeled by an analytical formulae S(T), which describes both the experiments and the thermal distribution expected in an astrophysical context. The parameters S0 and T0 are summarized in a table. Previous estimates for the sticking coefficient of H atoms are close to the new estimation; however, we find that, when isotopic effects are taken into account, the sticking coefficient variations can be as much as a factor of 2 at T=100 K.
△ Less
Submitted 11 February, 2012;
originally announced February 2012.
-
Determination of coronal temperatures from electron density profiles
Authors:
J. F. Lemaire
Abstract:
The most popular method for determining coronal temperatures is the scale-height-method (shm). It is based on electron density profiles inferred from White Light (WL) brightness measurements of the corona during solar eclipses. This method has been applied to several published coronal electron density models. The calculated temperature distributions reach a maximum at r > 1.3 RS, and therefore do…
▽ More
The most popular method for determining coronal temperatures is the scale-height-method (shm). It is based on electron density profiles inferred from White Light (WL) brightness measurements of the corona during solar eclipses. This method has been applied to several published coronal electron density models. The calculated temperature distributions reach a maximum at r > 1.3 RS, and therefore do not satisfy one of the conditions for applying the shm method. Another method is the hydrostatic equilibrium method (hst), which enables coronal temperature distributions to be determined, providing solutions to the hydrostatic equilibrium equation. The temperature maximas using the hst method are almost equal to those obtained using the shm method, but the temperature peak is always at significantly lower altitude when the hst-method is used than when the shm-method is used. A third and more recently developed method, dyn, can be used for the same published electron density profiles. The temperature distributions obtained using the dyn method are regular solutions of the hydrodynamic equations. They depend on the expansion velocity of the coronal plasma considered as a free input parameter in the calculations. At the base of the solar corona, where the coronal bulk velocity is small (subsonic), the dyn and hst methods give similar temperature values. Larger differences are found at higher altitudes where the expansion velocity approaches and exceeds the velocity of sound. This paper shows also the effects of the flow tube geometry and ion temperatures on the electron temperature distributions. The dyn method is a generalization of the hst method, where the electron temperature distribution tends to zero at infinite radial distances. It is a new diagnostic tool for determining from WL coronal observations the range of radial distances where the coronal heating rate is at its maximum.
△ Less
Submitted 18 January, 2012; v1 submitted 16 December, 2011;
originally announced December 2011.
-
Changes in the morphology of interstellar ice analogues after hydrogen atom exposure
Authors:
Mario Accolla,
Emanuele Congiu,
François Dulieu,
Giulio Manicò,
Henda Chaabouni,
Elie Matar,
Hakima Mokrane,
Jean Louis Lemaire,
Valerio Pirronello
Abstract:
The morphology of water ice in the interstellar medium is still an open question. Although accretion of gaseous water could not be the only possible origin of the observed icy mantles covering dust grains in cold molecular clouds, it is well known that water accreted from the gas phase on surfaces kept at 10 K forms ice films that exhibit a very high porosity. It is also known that in the dark clo…
▽ More
The morphology of water ice in the interstellar medium is still an open question. Although accretion of gaseous water could not be the only possible origin of the observed icy mantles covering dust grains in cold molecular clouds, it is well known that water accreted from the gas phase on surfaces kept at 10 K forms ice films that exhibit a very high porosity. It is also known that in the dark clouds H2 formation occurs on the icy surface of dust grains and that part of the energy (4.48 eV) released when adsorbed atoms react to form H2 is deposited in the ice. The experimental study described in the present work focuses on how relevant changes of the ice morphology result from atomic hydrogen exposure and subsequent recombination. Using the temperature-programmed desorption (TPD) technique and a method of inversion analysis of TPD spectra, we show that there is an exponential decrease in the porosity of the amorphous water ice sample following D-atom irradiation. This decrease is inversely proportional to the thickness of the ice and has a value of Phi_0 = 2 x 10^16 D-atoms/cm^2 per layer of H2O. We also use a model which confirms that the binding sites on the porous ice are destroyed regardless of their energy depth, and that the reduction of the porosity corresponds in fact to a reduction of the effective area. This reduction appears to be compatible with the fraction of D2 formation energy transferred to the porous ice network. Under interstellar conditions, this effect is likely to be efficient and, together with other compaction processes, provides a good argument to believe that interstellar ice is amorphous and non-porous.
△ Less
Submitted 21 March, 2011; v1 submitted 24 November, 2010;
originally announced November 2010.
-
VLT/NACO near-infrared imaging and spectroscopy of N88A in the SMC
Authors:
G. Testor,
J. L. Lemaire,
M. Heydari-Malayeri,
L. E. Kristensen,
S. Diana,
D. Field
Abstract:
We present near-infrared imaging and spectroscopic high spatial resolution observations of the SMC region N88 containing the bright, excited, extincted and compact H II region N88A of size ~ 1 pc. To investigate its stellar content and reddening, N88 was observed using spectroscopy and imagery in the JHKs- and L'-band at a spatial resolution of ~ 0.1-0.3", using the VLT UT4 equipped with the NAO…
▽ More
We present near-infrared imaging and spectroscopic high spatial resolution observations of the SMC region N88 containing the bright, excited, extincted and compact H II region N88A of size ~ 1 pc. To investigate its stellar content and reddening, N88 was observed using spectroscopy and imagery in the JHKs- and L'-band at a spatial resolution of ~ 0.1-0.3", using the VLT UT4 equipped with the NAOS adaptive optics system. In order to attempt to establish if the origin of the infra-red (IR) excess is due to bright nebulosity, circumstellar material and/or local dust, we used Ks vs J-K colour-magnitude (CM) and JHK colour-colour (CC) diagrams, as well as L' imagery.Our IR-data reveal in the N88 area an IR-excess fraction of geq 30 per cent of the detected stars,as well as an unprecedently detailed morphology of N88A. It consists of an embedded cluster of ~3.5" (~ 1 pc) in diameter, of at least thirteen resolved stars superposed with an unusual bright continuum centered on a very bright star. The four brightest stars in this cluster lie red-ward of H-K geq 0.45 mag, and could be classified as young stellar object (YSO) candidates. Four other probable YSO candidates are also detected in N88 along a south-north bow-shaped thin H2 filament at ~ 7" east of the young central bright star. At 0.2" east of this star, a heavily embedded core is detected that could be a massive class I protostar candidate. The 2.12 mu H2 image of N88A resembles a shell of diameter ~ 3" ~ 0.9 pc) centered on the bright star. The line ratios of H2 2-1 S(1) and 1-0 S(0) relative to 1-0 S(1), as well as the presence of high v lines, are indicative of photodissociation regions, rather than shocks.
△ Less
Submitted 29 November, 2009;
originally announced November 2009.
-
Experimental evidence of water formation on interstellar dust grains
Authors:
F. Dulieu,
L. Amiaud,
J-H. Fillion,
E. Matar,
A. Momeni,
V. Pirronello,
J. L. Lemaire
Abstract:
The synthesis of water is one necessary step in the origin and development of life. It is believed that pristine water is formed and grows on the surface of icy dust grains in dark interstellar clouds. Until now, there has been no experimental evidence whether this scenario is feasible or not. We present here the first experimental evidence of water synthesis under interstellar conditions. After…
▽ More
The synthesis of water is one necessary step in the origin and development of life. It is believed that pristine water is formed and grows on the surface of icy dust grains in dark interstellar clouds. Until now, there has been no experimental evidence whether this scenario is feasible or not. We present here the first experimental evidence of water synthesis under interstellar conditions. After D and O deposition on a water ice substrate (HO) held at 10 K, we observe production of HDO and DO. The water substrate itself has an active role in water formation, which appears to be more complicated than previously thought. Amorphous water ice layers are the matrices where complex organic prebiotic species may be synthesized. This experiment opens up the field of a little explored complex chemistry that could occur on interstellar dust grains, believed to be the site of key processes leading to the molecular diversity and complexity observed in our universe.
△ Less
Submitted 11 November, 2009;
originally announced November 2009.
-
Experimental evidence for water formation via ozone hydrogenation on dust grains at 10 K
Authors:
H. Mokrane,
H. Chaabouni,
M. Accolla,
E. Congiu,
F. Dulieu,
M. Chehrouri,
J. L. Lemaire
Abstract:
The formation of water molecules from the reaction between ozone (O3) and D-atoms is studied experimentally for the first time. Ozone is deposited on non-porous amorphous solid water ice (H2O), and D-atoms are then sent onto the sample held at 10 K. HDO molecules are detected during the desorption of the whole substrate where isotope mixing takes place, indicating that water synthesis has occurr…
▽ More
The formation of water molecules from the reaction between ozone (O3) and D-atoms is studied experimentally for the first time. Ozone is deposited on non-porous amorphous solid water ice (H2O), and D-atoms are then sent onto the sample held at 10 K. HDO molecules are detected during the desorption of the whole substrate where isotope mixing takes place, indicating that water synthesis has occurred. The efficiency of water formation via hydrogenation of ozone is of the same order of magnitude of that found for reactions involving O atoms or O2 molecules and exhibits no apparent activation barrier. These experiments validate the assumption made by models using ozone as one of the precursors of water formation via solid-state chemistry on interstellar dust grains.
△ Less
Submitted 5 October, 2009; v1 submitted 29 July, 2009;
originally announced July 2009.
-
Experimental evidence for water formation on interstellar dust grains by hydrogen and oxygen atoms
Authors:
F. Dulieu,
L. Amiaud,
E. Congiu,
J-H. Fillion,
E. Matar,
A. Momeni,
V. Pirronello,
J. L. Lemaire
Abstract:
Context. The synthesis of water is one necessary step in the origin and development of life. It is believed that pristine water is formed and grows on the surface of icy dust grains in dark interstellar clouds. Until now, there has been no experimental evidence whether this scenario is feasible or not on an astrophysically relevant template and by hydrogen and oxygen atom reactions.
Aims. We p…
▽ More
Context. The synthesis of water is one necessary step in the origin and development of life. It is believed that pristine water is formed and grows on the surface of icy dust grains in dark interstellar clouds. Until now, there has been no experimental evidence whether this scenario is feasible or not on an astrophysically relevant template and by hydrogen and oxygen atom reactions.
Aims. We present here the first experimental evidence of water synthesis by such a process on a realistic grain surface analogue in dense clouds, i.e., amorphous water ice.
Methods. Atomic beams of oxygen and deuterium are aimed at a porous water ice substrate (H2O) held at 10 K. Products are analyzed by the temperature-programmed desorption technique.
Results. We observe production of HDO and D2O, indicating that water is formed under conditions of the dense interstellar medium from hydrogen and oxygen atoms. This experiment opens up the field of a little explored complex chemistry that could occur on dust grains, believed to be the site where key processes lead to the molecular diversity and complexity observed in the Universe.
△ Less
Submitted 18 March, 2009;
originally announced March 2009.
-
Mobility of D atoms on porous amorphous water ice surfaces under interstellar conditions
Authors:
E. Matar,
E. Congiu,
F. Dulieu,
A. Momeni,
J. L. Lemaire
Abstract:
Aims. The mobility of H atoms on the surface of interstellar dust grains at low temperature is still a matter of debate. In dense clouds, the hydrogenation of adsorbed species (i.e., CO), as well as the subsequent deuteration of the accreted molecules depend on the mobility of H atoms on water ice. Astrochemical models widely assume that H atoms are mobile on the surface of dust grains even if c…
▽ More
Aims. The mobility of H atoms on the surface of interstellar dust grains at low temperature is still a matter of debate. In dense clouds, the hydrogenation of adsorbed species (i.e., CO), as well as the subsequent deuteration of the accreted molecules depend on the mobility of H atoms on water ice. Astrochemical models widely assume that H atoms are mobile on the surface of dust grains even if controversy still exists. We present here direct experimental evidence of the mobility of H atoms on porous water ice surfaces at 10 K. Methods. In a UHV chamber, O2 is deposited on a porous amorphous water ice substrate. Then D atoms are deposited onto the surface held at 10 K. Temperature-Programmed Desorption (TPD) is used and desorptions of O2 and D2 are simultaneously monitored. Results. We find that the amount of O2 that desorb during the TPD diminishes if we increase the deposition time of D atoms. O2 is thus destroyed by D atoms even though these molecules have previously diffused inside the pores of thick water ice. Our results can be easily interpreted if D is mobile at 10 K on the water ice surface. A simple rate equation model fits our experimental data and best fit curves were obtained for a D atoms diffusion barrier of 22(+-)2 meV. Therefore hydrogenation can take place efficiently on interstellar dust grains. These experimental results are in line with most calculations and validate the hypothesis used in several models.
△ Less
Submitted 13 October, 2008;
originally announced October 2008.
-
A method for detection of structure
Authors:
M. Gustafsson,
J. L. Lemaire,
D. Field
Abstract:
Context. In order to understand the evolution of molecular clouds it is important to identify the departures from self-similarity associated with the scales of self-gravity and the driving of turbulence.
Aims. A method is described based on structure functions for determining whether a region of gas, such as a molecular cloud, is fractal or contains structure with characteristic scale sizes.…
▽ More
Context. In order to understand the evolution of molecular clouds it is important to identify the departures from self-similarity associated with the scales of self-gravity and the driving of turbulence.
Aims. A method is described based on structure functions for determining whether a region of gas, such as a molecular cloud, is fractal or contains structure with characteristic scale sizes.
Methods. Using artificial data containing structure it is shown that derivatives of higher order structure functions provide a powerful way to detect the presence of characteristic scales should any be present and to estimate the size of such structures. The method is applied to observations of hot H2 in the Kleinman-Low nebula, north of the Trapezium stars in the Orion Molecular Cloud, including both brightness and velocity data. The method is compared with other techniques such as Fourier transform and histogram techniques.
Results. It is found that the density structure, represented by H2 emission brightness in the K-band (2-2.5micron), exhibits mean characteristic sizes of 110, 550, 1700 and 2700AU. The velocity data show the presence of structure at 140, 1500 and 3500AU. Compared with other techniques such as Fourier transform or histogram, the method appears both more sensitive to characteristic scales and easier to interpret.
△ Less
Submitted 7 June, 2006;
originally announced June 2006.
-
Oscillator Strengths and Predissociation Rates for Rydberg Transitions in 12C16O, 13C16O, and 13C18O Involving the E 1Pi, B 1Sigma+, and W 1Pi States
Authors:
M. Eidelsberg,
Y. Sheffer,
S. R. Federman,
J. L. Lemaire,
J. H. Fillion,
F. Rostas,
J. Ruiz
Abstract:
One of the processes controlling the interstellar CO abundance and the ratio of its isotopologues is photodissociation. Accurate oscillator strengths and predissociation rates for Rydberg transitions are needed for modeling this process. We present results on absorption from the E ^1Pi-X ^1Sigma^+ (1-0) and B ^1Sigma^+-X ^1Sigma^+ (6-0) bands at 1051 and 1002 Å, respectively, and the vibrational…
▽ More
One of the processes controlling the interstellar CO abundance and the ratio of its isotopologues is photodissociation. Accurate oscillator strengths and predissociation rates for Rydberg transitions are needed for modeling this process. We present results on absorption from the E ^1Pi-X ^1Sigma^+ (1-0) and B ^1Sigma^+-X ^1Sigma^+ (6-0) bands at 1051 and 1002 Å, respectively, and the vibrational progression W ^1Pi-X ^1Sigma^+ (v'-0) bands with v' = 0 to 3 at 972, 956, 941, and 925 Å, respectively. The corresponding spectra were acquired at the high resolution (R ~ 30,000) SU5 beam line at the Super ACO Synchrotron in Orsay, France. Spectra were obtained for the ^12C^16O, ^13C^16O, and ^13C^18O isotopologues. These represent the most complete set of measurements available. Comparison is made with earlier results, both empirical and theoretical. While earlier determinations of oscillator strengths based on absorption from synchrotron radiation tend to be somewhat smaller than ours, the suite of measurements from a variety of techniques agree for the most part considering the mutual uncertainties. For the bands studied here, their relative weakness, or their significant line widths arising from predissociation, minimizes potential problems from large optical depths at line center in absorption measurements. Predissociating line widths could generally be extracted from the spectra thanks to the profile simulations used in the analysis. In many cases, these simulations allowed us to consider e and f parity levels separately and to determine the dependence of the width on rotational quantum number, J. Our results are consistent with earlier determinations, especially the widths inferred from laser experiments.
△ Less
Submitted 6 May, 2006;
originally announced May 2006.
-
The nature of turbulence in OMC1 at the star forming scale: observations and simulations
Authors:
M. Gustafsson,
A. Brandenburg,
J. L. Lemaire,
D. Field
Abstract:
Aim: To study turbulence in the Orion Molecular Cloud (OMC1) by comparing observed and simulated characteristics of the gas motions.
Method: Using a dataset of vibrationally excited H2 emission in OMC1 containing radial velocity and brightness which covers scales from 70AU to 30000AU, we present the transversal structure functions and the scaling of the structure functions with their order. The…
▽ More
Aim: To study turbulence in the Orion Molecular Cloud (OMC1) by comparing observed and simulated characteristics of the gas motions.
Method: Using a dataset of vibrationally excited H2 emission in OMC1 containing radial velocity and brightness which covers scales from 70AU to 30000AU, we present the transversal structure functions and the scaling of the structure functions with their order. These are compared with the predictions of two-dimensional projections of simulations of supersonic hydrodynamic turbulence.
Results: The structure functions of OMC1 are not well represented by power laws, but show clear deviations below 2000AU. However, using the technique of extended self-similarity, power laws are recovered at scales down to 160AU. The scaling of the higher order structure functions with order deviates from the standard scaling for supersonic turbulence. This is explained as a selection effect of preferentially observing the shocked part of the gas and the scaling can be reproduced using line-of-sight integrated velocity data from subsets of supersonic turbulence simulations. These subsets select regions of strong flow convergence and high density associated with shock structure. Deviations of the structure functions in OMC1 from power laws cannot however be reproduced in simulations and remains an outstanding issue.
△ Less
Submitted 21 April, 2006; v1 submitted 8 December, 2005;
originally announced December 2005.
-
Observations of spatial and velocity structure in the Orion Molecular Cloud
Authors:
H. D. Nissen,
M. Gustafsson,
J. L. Lemaire,
Y. Clenet,
D. Rouan,
D. Field
Abstract:
Observations are reported of H2 IR emission in the S(1) v=1-0 line at 2.121 microns in the Orion Molecular Cloud, OMC1, using the GriF instrument on the Canada-France-Hawaii Telescope. GriF is a combination of adaptive optics and Fabry-Perot interferometry, yielding a spatial resolution of 0.15" to 0.18" and a velocity discrimination as high as 1 km/s. Thanks to the high spatial and velocity res…
▽ More
Observations are reported of H2 IR emission in the S(1) v=1-0 line at 2.121 microns in the Orion Molecular Cloud, OMC1, using the GriF instrument on the Canada-France-Hawaii Telescope. GriF is a combination of adaptive optics and Fabry-Perot interferometry, yielding a spatial resolution of 0.15" to 0.18" and a velocity discrimination as high as 1 km/s. Thanks to the high spatial and velocity resolution of the GriF data, 193 bright H2 emission regions can be identified in OMC1. The general characteristics of these features are described in terms of radial velocities, brightness and spatial displacement of maxima of velocity and brightness, the latter to yield the orientation of flows in the plane of the sky. Strong spatial correlation between velocity and bright H2 emission is found and serves to identify many features as shocks. Important results are: (i) velocities of the excited gas illustrate the presence of a zone to the south of BN-IRc2 and Peak 1, and the west of Peak 2, where there is a powerful blue-shifted outflow with an average velocity of -18 km/s. This is shown to be the NIR counterpart of an outflow identified in the radio from source I, a very young O-star. (ii) There is a band of weak velocity features (<5 km/s) in Peak 1 which may share a common origin through an explosive event, in the BN-IRc2 region, with the fast-moving fingers (or bullets) to the NW of OMC1. (iii) A proportion of the flows are likely to represent sites of low mass star formation and several regions show multiple outflows, probably indicative of multiple star formation within OMC1. The high spatial and velocity resolution of the GriF data show these and other features in more detail than has previously been possible.
△ Less
Submitted 28 February, 2006; v1 submitted 8 November, 2005;
originally announced November 2005.
-
Probing Turbulence with Infrared Observations in OMC1
Authors:
M. Gustafsson,
D. Field,
J. L. Lemaire,
F. P. Pijpers
Abstract:
A statistical analysis is presented of the turbulent velocity structure in the Orion Molecular Cloud at scales ranging from 70 AU to 30000 AU. Results are based on IR Fabry-Perot interferometric observations of shock and photon-excited H2 in the K-band S(1) v=1-0 line at 2.121micron and refer to the dynamical characteristics of warm perturbed gas. Observations establish that the Larson size-line…
▽ More
A statistical analysis is presented of the turbulent velocity structure in the Orion Molecular Cloud at scales ranging from 70 AU to 30000 AU. Results are based on IR Fabry-Perot interferometric observations of shock and photon-excited H2 in the K-band S(1) v=1-0 line at 2.121micron and refer to the dynamical characteristics of warm perturbed gas. Observations establish that the Larson size-linewidth relation is obeyed to the smallest scales studied here extending the range of validity of this relationship by nearly 2 orders of magnitude. The velocity probability distribution function (PDF) is constructed showing extended exponential wings, providing evidence of intermittency, further supported by the skewness and kurtosis of the velocity distribution. Variance and kurtosis of the PDF of velocity differences are constructed as a function of lag. The variance shows an approximate power law dependence on lag, with exponent significantly lower than the Kolmogorov value, and with deviations below 2000AU which are attributed to outflows and possibly disk structures associated with low mass star formation within OMC1. The kurtosis shows strong deviation from a gaussian velocity field, providing evidence of velocity correlations at small lags. Results agree accurately with semi-empirical simulations in Eggers & Wang (1998).
In addition, 170 individual H2 emitting clumps have been analysed with sizes between 500 and 2200 AU. These show considerable diversity with regard to PDFs and variance functions. Our analysis constitutes the first characterization of the turbulent velocity field at the scale of star formation and provide a dataset which models of star-forming regions should aim to reproduce.
△ Less
Submitted 13 December, 2005; v1 submitted 30 August, 2005;
originally announced August 2005.
-
Oscillator strengths for transitions to Rydberg levels in $^{12}C^{16}O$, $^{13}C^{16}O$ and $^{13}C^{18}O$ between 967 and 972 A
Authors:
M. Eidelsberg,
J. L. Lemaire,
J. H. Fillion,
F. Rostas,
S. Federman,
Y. Sheffer
Abstract:
Absorption oscillator strengths have been determined from high-resolution spectra in the 967-972 Åregion of three CO isotopomers for transitions to the Rydberg levels 4{\it p$π$}(0), 3{\it d$π$}(1) and 4{\it p$σ$}(0), as well as to the mixed {\it E(6)} level recently characterized by Eidelsberg et al. (2004). Synchrotron radiation from the Super-ACO electron storage ring at Orsay (LURE) was used…
▽ More
Absorption oscillator strengths have been determined from high-resolution spectra in the 967-972 Åregion of three CO isotopomers for transitions to the Rydberg levels 4{\it p$π$}(0), 3{\it d$π$}(1) and 4{\it p$σ$}(0), as well as to the mixed {\it E(6)} level recently characterized by Eidelsberg et al. (2004). Synchrotron radiation from the Super-ACO electron storage ring at Orsay (LURE) was used as a light source. Oscillator strengths were extracted from the recorded spectra by least-squares fitting of the experimental profiles with synthetic spectra taking into account the homogeneous and heterogeneous interactions of the four levels. Column densities were derived from fits to the 3{\it p$π$}(0) absorption band whose oscillator strength is well established. These are the first reported measurements for $^{13}$C$^{18}$O. For $^{12}$C$^{16}$O, our results are consistent with the larger values obtained in the most recent laboratory and astronomical studies.
△ Less
Submitted 31 May, 2004;
originally announced May 2004.
-
VLT/NACO infrared adaptive optics images of small scale structures in OMC1
Authors:
Francois Lacombe,
Eric Gendron,
Daniel Rouan,
Yann Clenet,
Jean-Louis Lemaire,
Anne-Marie Lagrange,
David Mouillet,
Gerard Rousset,
Claude Marlot,
Philippe Feautrier,
M. Gustafsson,
David Field,
Bernard Servan
Abstract:
Near-infrared observations of line emission from excited H2 and in the continuum are reported in the direction of the Orion molecular cloud OMC1, using the European Southern Observatory Very Large Telescope UT4, equipped with the NAOS adaptive optics system and the CONICA infrared array camera. Spatial resolution has been achieved at close to the diffraction limit of the telescope (0.08" - 0.12"…
▽ More
Near-infrared observations of line emission from excited H2 and in the continuum are reported in the direction of the Orion molecular cloud OMC1, using the European Southern Observatory Very Large Telescope UT4, equipped with the NAOS adaptive optics system and the CONICA infrared array camera. Spatial resolution has been achieved at close to the diffraction limit of the telescope (0.08" - 0.12") and images show a wealth of morphological detail. Structure is not fractal but shows two preferred scale sizes of 2.4" (1100 AU) and 1.2" (540 AU), where the larger scale may be associated with star formation.
△ Less
Submitted 3 December, 2003;
originally announced December 2003.
-
High Resolution HST-STIS Spectra of CI and CO in the Beta Pictoris Circumstellar Disk
Authors:
A. Roberge,
P. D. Feldman,
A. M. Lagrange,
A. Vidal-Madjar,
R. Ferlet,
A. Jolly,
J. L. Lemaire,
F. Rostas
Abstract:
High resolution FUV echelle spectra showing absorption features arising from CI and CO gas in the Beta Pictoris circumstellar (CS) disk were obtained on 1997 December 6 and 19 using the Space Telescope Imaging Spectrograph (STIS). An unsaturated spin-forbidden line of CI at 1613.376 A not previously seen in spectra of Beta Pictoris was detected, allowing for an improved determination of the colu…
▽ More
High resolution FUV echelle spectra showing absorption features arising from CI and CO gas in the Beta Pictoris circumstellar (CS) disk were obtained on 1997 December 6 and 19 using the Space Telescope Imaging Spectrograph (STIS). An unsaturated spin-forbidden line of CI at 1613.376 A not previously seen in spectra of Beta Pictoris was detected, allowing for an improved determination of the column density of CI at zero velocity relative to the star (the stable component), N = (2-4) x 10^{16} cm^{-2}. Variable components with multiple velocities, which are the signatures of infalling bodies in the Beta Pictoris CS disk, are observed in the CI 1561 A and 1657 A multiplets. Also seen for the first time were two lines arising from the metastable singlet D level of carbon, at 1931 A and 1463 A The results of analysis of the CO A-X (0-0), (1-0), and (2-0) bands are presented, including the bands arising from {13}^CO, with much better precision than has previously been possible, due to the very high resolution provided by the STIS echelle gratings. Only stable CO gas is observed, with a column density N(CO) = (6.3 +/- 0.3) x 10^{14} cm{-2}. An unusual ratio of the column densities of {12}^CO to {13}^CO is found (R = 15 +/- 2). The large difference between the column densities of CI and CO indicates that photodissociation of CO is not the primary source of CI gas in the disk, contrary to previous suggestion.
△ Less
Submitted 29 March, 2000;
originally announced March 2000.