Showing 201–245 of 245 results for author: Mathis, S

Accretion of Saturn's mid-sized moons during the viscous spreading of young massive rings: solving the paradox of silicate-poor rings versus silicate-rich moons

Authors: Sébastien Charnoz, Aurélien Crida, Julie C. Castillo-Rogez, Valery Lainey, Luke Dones, Özgür Karatekin, Gabriel Tobie, Stephane Mathis, Christophe Le Poncin-Lafitte, Julien Salmon

Abstract: The origin of Saturn's inner mid-sized moons (Mimas, Enceladus, Tethys, Dione and Rhea) and Saturn's rings is debated. Charnoz et al. (2010) introduced the idea that the smallest inner moons could form from the spreading of the rings' edge while Salmon et al. (2010) showed that the rings could have been initially massive, and so was the ring's progenitor itself. One may wonder if the mid-sized moo… ▽ More The origin of Saturn's inner mid-sized moons (Mimas, Enceladus, Tethys, Dione and Rhea) and Saturn's rings is debated. Charnoz et al. (2010) introduced the idea that the smallest inner moons could form from the spreading of the rings' edge while Salmon et al. (2010) showed that the rings could have been initially massive, and so was the ring's progenitor itself. One may wonder if the mid-sized moons may have formed also from the debris of a massive ring progenitor, as also suggested in Canup (2010). However, the process driving mid-sized moons accretion from the icy debris disks has not been investigated in details. In particular, this process does not seem able to explain the varying silicate contents of the mid-sized moons (from 6% to 57% in mass). Here, we explore the formation of large objects from a massive ice-rich ring (a few times Rhea's mass) and describe the fundamental properties and implications of this new process. Using a hybrid computer model, we show that accretion within massive icy rings can form all mid-sized moons from Mimas to Rhea. However in order to explain their current locations, intense dissipation within Saturn (with Qp<2000) would be required. Our results are consistent with a satellite origin tied to the rings formation at least 2.5 Gy ago, both compatible with either a formation concurrent to Saturn or during the Late Heavy Bombardment. Tidal heating related to high-eccentricity post-accretional episodes may induce early geological activity. If some massive irregular chunks of silicates were initially present within the rings, they would be present today inside the satellites' cores which would have accreted icy shells while being tidally expelled from the rings (via a heterogeneous accretion process)while those still present in the rings are interpreted as today Saturn's rings' propellers and ring-moons (like Pan or Daphnis). △ Less

Submitted 15 September, 2011; originally announced September 2011.

Comments: 41 pages, 14 figures. Accepted in ICARUS

The dynamics of the radiative zone of the Sun

Authors: S Turck-Chieze, S Couvidat, V Duez, S Mathis, J Marques, A Palacios, L Piau

Abstract: Helioseismology puts strong constraints on the internal sound speed and on the rotation profile in the radiative zone. Young stars of solar type are more active and faster rotators than the Sun. So we begin to build models which include different rotation histories and compare the results with all the solar observations. The profiles of the rotation we get have interesting consequence for the intr… ▽ More Helioseismology puts strong constraints on the internal sound speed and on the rotation profile in the radiative zone. Young stars of solar type are more active and faster rotators than the Sun. So we begin to build models which include different rotation histories and compare the results with all the solar observations. The profiles of the rotation we get have interesting consequence for the introduction of magnetic field in the radiative zone. We discuss also the impact of mass loss deduced from measured flux of young stars. We deduce from these comparisons some quantitative effect of the dynamical processes (rotation, magnetic field and mass loss) of these early stages on the present sound speed and density. We show finally how we can improve our present knowledge of the radiative zone with PICARD and GOLFNG. △ Less

Submitted 6 January, 2011; originally announced January 2011.

Comments: 4 pages, 2 figures, contribution to SOHO24 in press

Journal ref: Journal of Physics: Conf. Series 2011

New insights on the solar core

Authors: R. A. Garcia, D. Salabert, J. Ballot, A. Eff-Darwich, R. Garrido, A. Jimenez, S. Mathis, S. Mathur, A. Moya, P. L. Palle, C. Regulo, K. Sato, J. C. Suarez, S. Turck-Chieze

Abstract: Since the detection of the asymptotic properties of the dipole gravity modes in the Sun, the quest to find individual gravity modes has continued. An extensive and deeper analysis of 14 years of continuous GOLF/SoHO observational data, unveils the presence of a pattern of peaks that could be interpreted as individual dipole gravity modes in the frequency range between 60 and 140 microHz, with ampl… ▽ More Since the detection of the asymptotic properties of the dipole gravity modes in the Sun, the quest to find individual gravity modes has continued. An extensive and deeper analysis of 14 years of continuous GOLF/SoHO observational data, unveils the presence of a pattern of peaks that could be interpreted as individual dipole gravity modes in the frequency range between 60 and 140 microHz, with amplitudes compatible with the latest theoretical predictions. By collapsing the power spectrum we have obtained a quite constant splitting for these patterns in comparison to regions where no g modes were expected. Moreover, the same technique applied to simultaneous VIRGO/SoHO data unveils some common signals between the power spectra of both instruments. Thus, we are able to identify and characterize individual g modes with their central frequencies, amplitudes and splittings allowing to do seismic inversions of the rotation profile inside the solar core. These results open a new light on the physics and dynamics of the solar deep core. △ Less

Submitted 2 December, 2010; originally announced December 2010.

Comments: Proceedings of the SoHO-24 meeting. 6 pages and 7 figures

Journal ref: Journal of Physics: Conference Series, 2011, Volume 271

A first family of stable configurations to model magnetic fields in stellar radiation zones

Authors: V. Duez, J. Braithwaite, S. Mathis

Abstract: We conduct 3D magneto-hydrodynamic (MHD) simulations in order to test the stability of the magnetic equilibrium configuration described by Duez & Mathis (2010). This analytically-derived configuration describes the lowest energy state for a given helicity in a stellar radiation zone. The necessity of taking into account the non force-free property of the large-scale, global field is here emphasize… ▽ More We conduct 3D magneto-hydrodynamic (MHD) simulations in order to test the stability of the magnetic equilibrium configuration described by Duez & Mathis (2010). This analytically-derived configuration describes the lowest energy state for a given helicity in a stellar radiation zone. The necessity of taking into account the non force-free property of the large-scale, global field is here emphasized. We then show that this configuration is stable. It therefore provides a useful model to initialize the magnetic topology in upcoming MHD simulations and stellar evolution codes taking into account magneto-rotational transport processes. △ Less

Submitted 28 September, 2010; originally announced September 2010.

Comments: 2 pages, 1 figure, to be published in the proceedings of IAU 272 "Active OB stars: structure, evolution, mass loss and critical limit", Paris 19-23 July 2010

On the Stability of Non Force-Free Magnetic Equilibria in Stars

Authors: V. Duez, J. Braithwaite, S. Mathis

Abstract: The existence of stable magnetic configurations in white dwarfs, neutron stars and various non-convective stellar {regions} is now well recognized. It has recently been shown numerically that various families of equilibria, including axisymmetric mixed poloidal-toroidal configurations, are stable. Here we test the stability of an analytically-derived non force-free magnetic equilibrium, using thre… ▽ More The existence of stable magnetic configurations in white dwarfs, neutron stars and various non-convective stellar {regions} is now well recognized. It has recently been shown numerically that various families of equilibria, including axisymmetric mixed poloidal-toroidal configurations, are stable. Here we test the stability of an analytically-derived non force-free magnetic equilibrium, using three-dimensional magnetohydrodynamic simulations: the mixed configuration is compared with the dynamical evolution of its purely poloidal and purely toroidal components, both known to be unstable. The mixed equilibrium shows no sign of instability under white noise perturbations. {This configuration therefore provides a good description of magnetic equilibrium topology inside non-convective stellar objects and will be useful to initialize magneto-rotational transport in stellar evolution codes. △ Less

Submitted 27 September, 2010; originally announced September 2010.

Comments: Accepted for publication in ApJL

Effects of rotational mixing on the asteroseismic properties of solar-type stars

Authors: P. Eggenberger, G. Meynet, A. Maeder, A. Miglio, J. Montalban, F. Carrier, S. Mathis, C. Charbonnel, S. Talon

Abstract: The influence of rotational mixing on the evolution and asteroseismic properties of solar-type stars is studied. Rotational mixing changes the global properties of a solar-type star with a significant increase of the effective temperature resulting in a shift of the evolutionary track to the blue part of the HR diagram. These differences are related to changes of the chemical composition, because… ▽ More The influence of rotational mixing on the evolution and asteroseismic properties of solar-type stars is studied. Rotational mixing changes the global properties of a solar-type star with a significant increase of the effective temperature resulting in a shift of the evolutionary track to the blue part of the HR diagram. These differences are related to changes of the chemical composition, because rotational mixing counteracts the effects of atomic diffusion leading to larger helium surface abundances for rotating models than for non-rotating ones. Higher values of the large frequency separation are then found for rotating models than for non-rotating ones at the same evolutionary stage, because the increase of the effective temperature leads to a smaller radius and hence to an increase of the stellar mean density. Rotational mixing also has a considerable impact on the structure and chemical composition of the central stellar layers by bringing fresh hydrogen fuel to the core, thereby enhancing the main-sequence lifetime. The increase of the central hydrogen abundance together with the change of the chemical profiles in the central layers result in a significant increase of the values of the small frequency separations and of the ratio of the small to large separations for models including shellular rotation. This increase is clearly seen for models with the same age sharing the same initial parameters except for the inclusion of rotation as well as for models with the same global stellar parameters and in particular the same location in the HR diagram. By computing rotating models of solar-type stars including the effects of a dynamo that possibly occurs in the radiative zone, we find that the efficiency of rotational mixing is strongly reduced when the effects of magnetic fields are taken into account, in contrast to what happens in massive stars. △ Less

Submitted 23 September, 2010; originally announced September 2010.

Comments: 11 pages, 15 figures, accepted for publication in A&A

The MiMeS Project: Overview and Current Status

Authors: G. A. Wade, E. Alecian, D. A. Bohlender, J. -C. Bouret, D. Cohen, V. Duez, M. Gagné, J. Grunhut, H. F. Henrichs, N. Hill, O. Kochukhov, S. Mathis, C. Neiner, M. Oksala, S. Owocki, V. Petit, M. Shultz, T. Rivinius, R. Townsend, J. S. Vink, the MiMeS Collaboration

Abstract: The Magnetism in Massive Stars (MiMeS) Project is a consensus collaboration among many of the foremost international researchers of the physics of hot, massive stars, with the basic aim of understanding the origin, evolution and impact of magnetic fields in these objects. At the time of writing, MiMeS Large Programs have acquired over 950 high-resolution polarised spectra of about 150 individual s… ▽ More The Magnetism in Massive Stars (MiMeS) Project is a consensus collaboration among many of the foremost international researchers of the physics of hot, massive stars, with the basic aim of understanding the origin, evolution and impact of magnetic fields in these objects. At the time of writing, MiMeS Large Programs have acquired over 950 high-resolution polarised spectra of about 150 individual stars with spectral types from B5-O4, discovering new magnetic fields in a dozen hot, massive stars. The quality of this spectral and magnetic matériel is very high, and the Collaboration is keen to connect with colleagues capable of exploiting the data in new or unforeseen ways. In this paper we review the structure of the MiMeS observing programs and report the status of observations, data modeling and development of related theory. △ Less

Submitted 18 September, 2010; originally announced September 2010.

Comments: Proceedings of IAUS272: Active OB stars

Towards a detection of individual g modes in the Sun

Authors: R. A. Garcia, J. Ballot, A. Eff-Darwich, R. Garrido, A. Jimenez, S. Mathis, S. Mathur, A. Moya, P. L. Palle, C. Regulo, D. Salabert, J. C. Suarez, S. Turck-Chieze

Abstract: Since the detection of the asymptotic properties of the dipole gravity modes in the Sun, the quest to find the individual gravity modes has continued. A deeper analysis of the GOLF/SoHO data unveils the presence of a pattern of peaks that could be interpreted as individual dipole gravity modes. The computed collapsed spectrum -around these candidate modes- uncovers the presence of a quasi constant… ▽ More Since the detection of the asymptotic properties of the dipole gravity modes in the Sun, the quest to find the individual gravity modes has continued. A deeper analysis of the GOLF/SoHO data unveils the presence of a pattern of peaks that could be interpreted as individual dipole gravity modes. The computed collapsed spectrum -around these candidate modes- uncovers the presence of a quasi constant frequency splitting, in contrast with regions where no g modes are expected in which the collapsogram gives random results. Besides, the same technique applied to VIRGO/SoHO unveils some common signals between both power spectra. Thus, we can identify and characterize the modes, for example, with their central frequency and splittings. This would open the path towards new investigations to better constrain the solar core. △ Less

Submitted 26 July, 2010; originally announced July 2010.

Comments: Proceedings of the HELAS IV meeting. 4 pages, 5 Figures

Asteroseismology of Solar-type Stars with Kepler I: Data Analysis

Authors: C. Karoff, W. J. Chaplin, T. Appourchaux, Y. Elsworth, R. A. Garcia, G. Houdek, T. S. Metcalfe, J. Molenda-Zakowicz, M. J. P. F. G. Monteiro, M. J. Thompson, J. Christensen-Dalsgaard, R. L. Gilliland, H. Kjeldsen, S. Basu, T. R. Bedding, T. L. Campante, P. Eggenberger, S. T. Fletcher, P. Gaulme, R. Handberg, S. Hekker, M. Martic, S. Mathur, B. Mosser, C. Regulo , et al. (24 additional authors not shown)

Abstract: We report on the first asteroseismic analysis of solar-type stars observed by Kepler. Observations of three G-type stars, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation can clearly be distinguished in each star. We discuss the appearance of the oscillation spectra,… ▽ More We report on the first asteroseismic analysis of solar-type stars observed by Kepler. Observations of three G-type stars, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation can clearly be distinguished in each star. We discuss the appearance of the oscillation spectra, including the presence of a possible signature of faculae, and the presence of mixed modes in one of the three stars. △ Less

Submitted 19 July, 2010; v1 submitted 4 May, 2010; originally announced May 2010.

Comments: 5 pages, 4 figure, submitted to Astronomische Nachrichten

CP and related phenomena in the context of Stellar Evolution

Authors: J. Braithwaite, T. Akgün, E. Alecian, A. F. Kholtygin, J. D. Landstreet, S. Mathis, G. Michaud, J. Portnoy, G. Alecian, V. D. Bychkov, L. V. Bychkova, N. Drake, S. N. Fabrika, A. Reisenegger, R. Steinitz, M. Vick

Abstract: We review the interaction in intermediate and high mass stars between their evolution and magnetic and chemical properties. We describe the theory of Ap-star `fossil' fields, before touching on the expected secular diffusive processes which give rise to evolution of the field. We then present recent results from a spectropolarimetric survey of Herbig Ae/Be stars, showing that magnetic fields of… ▽ More We review the interaction in intermediate and high mass stars between their evolution and magnetic and chemical properties. We describe the theory of Ap-star `fossil' fields, before touching on the expected secular diffusive processes which give rise to evolution of the field. We then present recent results from a spectropolarimetric survey of Herbig Ae/Be stars, showing that magnetic fields of the kind seen on the main-sequence already exist during the pre-main sequence phase, in agreement with fossil field theory, and that the origin of the slow rotation of Ap/Bp stars also lies early in the pre-main sequence evolution; we also present results confirming a lack of stars with fields below a few hundred gauss. We then seek which macroscopic motions compete with atomic diffusion in determining the surface abundances of AmFm stars. While turbulent transport and mass loss, in competition with atomic diffusion, are both able to explain observed surface abundances, the interior abundance distribution is different enough to potentially lead to a test using asterosismology. Finally we review progress on the turbulence-driving and mixing processes in stellar radiative zones. △ Less

Submitted 22 February, 2010; originally announced February 2010.

Comments: Proceedings of IAU GA in Rio, JD4 on Ap stars; 10 pages, 7 figures

The asteroseismic potential of Kepler: first results for solar-type stars

Authors: W. J. Chaplin, T. Appourchaux, Y. Elsworth, R. A. Garcia, G. Houdek, C. Karoff, T. S. Metcalfe, J. Molenda-Zakowicz, M. J. P. F. G. Monteiro, M. J. Thompson, T. M. Brown, J. Christensen-Dalsgaard, R. L. Gilliland, H. Kjeldsen, W. J. Borucki, D. Koch, J. M. Jenkins, J. Ballot, S. Basu, M. Bazot, T. R. Bedding, O. Benomar, A. Bonanno, I. M. Brandao, H. Bruntt , et al. (83 additional authors not shown)

Abstract: We present preliminary asteroseismic results from Kepler on three G-type stars. The observations, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation may be clearly distinguished in each star. We discuss the appearance of the oscillation spectra, use the frequencies a… ▽ More We present preliminary asteroseismic results from Kepler on three G-type stars. The observations, made at one-minute cadence during the first 33.5d of science operations, reveal high signal-to-noise solar-like oscillation spectra in all three stars: About 20 modes of oscillation may be clearly distinguished in each star. We discuss the appearance of the oscillation spectra, use the frequencies and frequency separations to provide first results on the radii, masses and ages of the stars, and comment in the light of these results on prospects for inference on other solar-type stars that Kepler will observe. △ Less

Submitted 18 January, 2010; v1 submitted 4 January, 2010; originally announced January 2010.

Comments: 16 pages, 3 figures, 1 table; now accepted for publication in ApJ Letters

Detecting individual gravity modes in the Sun: Chimera or reality?

Authors: Rafael A. Garcia, J. Ballot, A. Eff-Darwich, R. Garrido, A. Jimenez, S. Mathis, S. Mathur, A. Moya, P. L. Palle, C. Regulo, D. Salabert, K. Sato, J. C. Suarez, S. Turck-Chieze

Abstract: Over the past 15 years, our knowledge of the interior of the Sun has tremendously progressed by the use of helioseismic measurements. However, to go further in our understanding of the solar core, we need to measure gravity (g) modes. Thanks to the high quality of the Doppler-velocity signal measured by GOLF/SoHO, it has been possible to unveil the signature of the asymptotic properties of the s… ▽ More Over the past 15 years, our knowledge of the interior of the Sun has tremendously progressed by the use of helioseismic measurements. However, to go further in our understanding of the solar core, we need to measure gravity (g) modes. Thanks to the high quality of the Doppler-velocity signal measured by GOLF/SoHO, it has been possible to unveil the signature of the asymptotic properties of the solar g modes, thus obtaining a hint of the rotation rate in the core. However, the quest for the detection of individual g modes is not yet over. In this work, we apply the latest theoretical developments to guide our research using GOLF velocity time series. In contrary to what was thought till now, we are maybe starting to identify individual low-frequency g modes... △ Less

Submitted 30 November, 2009; originally announced November 2009.

Comments: Highlight of Astronomy (HoA) proceedings of the JD-11, IAU 2009. 2 pages, 1 figure

Detecting individual gravity modes in the Sun

Authors: R. A. Garcia, J. Ballot, A. Eff-Darwich, R. Garrido, A. Jimenez, S. Mathis, A. Moya, P. L. Palle, C. Regulo, D. Salabert, J. C. Suarez, S. Turck-Chieze

Abstract: Many questions are still open regarding the structure and the dynamics of the solar core. By constraining more this region in the solar evolution models, we can reduce the incertitudes on some physical processes and on momentum transport mechanisms. A first big step was made with the detection of the signature of the dipole-gravity modes in the Sun, giving a hint of a faster rotation rate inside… ▽ More Many questions are still open regarding the structure and the dynamics of the solar core. By constraining more this region in the solar evolution models, we can reduce the incertitudes on some physical processes and on momentum transport mechanisms. A first big step was made with the detection of the signature of the dipole-gravity modes in the Sun, giving a hint of a faster rotation rate inside the core. A deeper analysis of the GOLF/SoHO data unveils the presence of a pattern of peaks that could be interpreted as dipole gravity modes. In that case, those modes can be characterized, thus bringing better constraints on the rotation of the core as well as some structural parameters such as the density at these very deep layers of the Sun interior. △ Less

Submitted 30 November, 2009; originally announced November 2009.

Comments: Proc. HELAS Workshop on `New insights into the Sun', eds M. Cunha, M. Thompson

Effect of a fossil magnetic field on the structure of a young Sun

Authors: Vincent Duez, Stéphane Mathis, Sylvaine Turck-Chièze

Abstract: We study the impact of a fossil magnetic field on the physical quantities which describe the structure of a young Sun of 500 Myr. We consider for the first time a non force-free field composed of a mixture of poloidal and toroidal magnetic fields and we propose a specific configuration to illustrate our purpose. In the present paper, we estimate the relative role of the different terms which app… ▽ More We study the impact of a fossil magnetic field on the physical quantities which describe the structure of a young Sun of 500 Myr. We consider for the first time a non force-free field composed of a mixture of poloidal and toroidal magnetic fields and we propose a specific configuration to illustrate our purpose. In the present paper, we estimate the relative role of the different terms which appear in the modified stellar structure equations. We note that the Lorentz tension plays a non negligible role in addition to the magnetic pressure. This is interesting because most of the previous stellar evolution codes ignored that term and the geometry of the field. The solar structure perturbations are, as already known, small and consequently we have been able to estimate each term semi-analytically. We develop a general treatment to calculate the global modification of the structure and of the energetic balance. We estimate also the gravitational multipolar moments associated with the presence of a fossil large-scale magnetic field in radiative zone. The values given for the young Sun help the future implementation in stellar evolution codes. This work can be repeated for any other field configuration and prepares the achievement of a solar MHD model where we will follow the transport of such field on secular timescales and the associated transport of momentum and chemicals. The described method will be applied at the present Sun and the results will be compared with the coming balloon or space measurements. △ Less

Submitted 4 November, 2009; originally announced November 2009.

Comments: 11 pages, accepted for publication in MNRAS

Mode excitation by turbulent convection in rotating stars. I. Effect of uniform rotation

Authors: K. Belkacem, S. Mathis, M. J. Goupil, R. Samadi

Abstract: We focus on the influence of the Coriolis acceleration on the stochastic excitation of oscillation modes in convective regions of rotating stars. Our aim is to estimate the asymmetry between excitation rates of prograde and retrograde modes. We extend the formalism derived for obtaining stellar $p$- and $g$-mode amplitudes (Samadi & Goupil 2001, Belkacem et al. 2008) to include the effect of the… ▽ More We focus on the influence of the Coriolis acceleration on the stochastic excitation of oscillation modes in convective regions of rotating stars. Our aim is to estimate the asymmetry between excitation rates of prograde and retrograde modes. We extend the formalism derived for obtaining stellar $p$- and $g$-mode amplitudes (Samadi & Goupil 2001, Belkacem et al. 2008) to include the effect of the Coriolis acceleration. We then study the special case of uniform rotation for slowly rotating stars by performing a perturbative analysis. This allows us to consider the cases of the Sun and the CoRoT target HD 49933. We find that, in the subsonic regime, the influence of rotation as a direct contribution to mode driving is negligible in front of the Reynolds stress contribution. In slow rotators, the indirect effect of the modification of the eigenfunctions on mode excitation is investigated by performing a perturbative analysis of the excitation rates. It turns out that the excitation of solar $p$ modes is affected by rotation with excitation rates asymmetries between prograde and retrograde modes of the order of several percents. Solar low-order $g$ modes are also affected by uniform rotation and their excitation rates asymmetries are found to reach up to 10 %. The CoRoT target HD 49933 is rotating faster than the Sun ($Ω/ Ω_\odot \approx 8$) and we show that the resulting excitation rates asymmetry is about 10 % for the excitation rates of $p$ modes. We have then demonstrated that $p$ and $g$ mode excitation rates are modified by uniform rotation through the Coriolis acceleration. Study of the effect of differential rotation is dedicated to a forthcoming paper. △ Less

Submitted 8 October, 2009; v1 submitted 8 September, 2009; originally announced September 2009.

Comments: 9 pages, 4 figures, accepted in A&A

Large-Scale Magnetic Fields in Stellar Interiors I. Equilibrium Configurations to Model Fossil Fields

Authors: V. Duez, S. Mathis

Abstract: Context. The understanding of fossil fields origin, topology and stability is one of the corner stones of the stellar magnetism theory. On one hand, since they survive over secular time-scales, they may modify the structure and the evolution of their host stars. On the other hand, they must have a complex stable structure since it has been demonstrated by Tayler and collaborators that simplest pur… ▽ More Context. The understanding of fossil fields origin, topology and stability is one of the corner stones of the stellar magnetism theory. On one hand, since they survive over secular time-scales, they may modify the structure and the evolution of their host stars. On the other hand, they must have a complex stable structure since it has been demonstrated by Tayler and collaborators that simplest purely poloidal or toroidal fields are unstable on dynamical time-scales. In this context, the only stable configuration which has been found today is the one resulting of a numerical simulation by Braithwaite and collaborators who have studied the evolution of an initial stochastic magnetic field, which is found to relax on a mixed stable configuration (poloidal and toroidal) that seems to be in equilibrium and then diffuses. Aims. In this work, we thus go on the track of such type of equilibrium field in a semi-analytical way. Methods. In this first article, we study the barotropic magnetohydrostatic equilibrium states; the problem reduces to a Grad-Shafranov-like equation with arbitrary functions. Those latters are constrained by deriving the lowest-energy equilibrium states for given invariants of the considered axisymmetric problem and in particular for a given helicity which is known to be one of the main actor of such problems. Then, we obtain the generalization of the force-free Taylor's relaxation states obtained in laboratory experiments (in spheromaks) that become non force-free in the self-gravitating stellar case. The case of general baroclinic equilibrium states will be studied in Paper II. Results. Those theoretical results are applied to realistic stellar cases, namely to the solar radiative core and to the envelope of an Ap star, and discussed. In both cases we assume that the field is initially confined in the stellar radiation zone. △ Less

Submitted 15 March, 2010; v1 submitted 9 April, 2009; originally announced April 2009.

Comments: Accepted for publication by A&A on the 6th of March, 2010.

Diagnoses to unravel secular hydrodynamical processes in rotating main sequence stars

Authors: T. Decressin, S. Mathis, A. Palacios, L. Siess, S. Talon, C. Charbonnel, J. -P. Zahn

Abstract: (Abridged) We present a detailed analysis of the main physical processes responsible for the transport of angular momentum and chemical species in the radiative regions of rotating stars. We focus on cases where meridional circulation and shear-induced turbulence only are included in the simulations. Our analysis is based on a 2-D representation of the secular hydrodynamics, which is treated usi… ▽ More (Abridged) We present a detailed analysis of the main physical processes responsible for the transport of angular momentum and chemical species in the radiative regions of rotating stars. We focus on cases where meridional circulation and shear-induced turbulence only are included in the simulations. Our analysis is based on a 2-D representation of the secular hydrodynamics, which is treated using expansions in spherical harmonics. We present a full reconstruction of the meridional circulation and of the associated fluctuations of temperature and mean molecular weight along with diagnosis for the transport of angular momentum, heat and chemicals. In the present paper these tools are used to validate the analysis of two main sequence stellar models of 1.5 and 20 Msun for which the hydrodynamics has been previously extensively studied in the literature. We obtain a clear visualization and a precise estimation of the different terms entering the angular momentum and heat transport equations in radiative zones. This enables us to corroborate the main results obtained over the past decade by Zahn, Maeder, and collaborators concerning the secular hydrodynamics of such objects. We focus on the meridional circulation driven by angular momentum losses and structural readjustements. We confirm quantitatively for the first time through detailed computations and separation of the various components that the advection of entropy by this circulation is very well balanced by the barotropic effects and the thermal relaxation during most of the main sequence evolution. This enables us to derive simplifications for the thermal relaxation on this phase. The meridional currents in turn advect heat and generate temperature fluctuations that induce differential rotation through thermal wind thus closing the transport loop. △ Less

Submitted 1 December, 2008; originally announced December 2008.

Comments: 16 pages, 18 figures. Accepted for publication in A&A

Is it necessary to go beyond the ponctual mass approximation for tidal perturber in close systems?

Authors: Stephane Mathis, Christophe Le Poncin-lafitte

Abstract: With the discoveries of very close star-planet systems with planet orbiting sometimes at several star radius but also with well-known situations in our solar system where natural satellites are very close to their parent planet the validity of the ponctual mass approximation for the tidal perturber (respectively the parent star or planet when we study the close planet or natural satellite dynami… ▽ More With the discoveries of very close star-planet systems with planet orbiting sometimes at several star radius but also with well-known situations in our solar system where natural satellites are very close to their parent planet the validity of the ponctual mass approximation for the tidal perturber (respectively the parent star or planet when we study the close planet or natural satellite dynamics) has to be examined. In this short paper, we consider this problematic using results coming from a complete formalism that allows to treat the tidal interaction between extended bodies. We focus on its application to a simplified configuration. △ Less

Submitted 13 October, 2008; originally announced October 2008.

Tidal dynamics of extended bodies in planetary systems and multiple stars

Authors: Stephane Mathis, Christophe Le Poncin-Lafitte

Abstract: With the discovery during the past decade of a large number of extrasolar planets orbiting their parent stars at a distance lower than 0.1 astronomical unit (and the launch and the preparation of dedicated space missions such as CoRoT and KEPLER), with the position of inner natural satellites around giant planets in our Solar System and with the existence of very closed but separated binary star… ▽ More With the discovery during the past decade of a large number of extrasolar planets orbiting their parent stars at a distance lower than 0.1 astronomical unit (and the launch and the preparation of dedicated space missions such as CoRoT and KEPLER), with the position of inner natural satellites around giant planets in our Solar System and with the existence of very closed but separated binary stars, tidal interaction has to be carefully studied. In particular, a question arises about the validity of usual approximations used in the modelling of this interaction. The purpose of this paper is to examine the step beyond the ponctual approximation for the tidal perturber. To achieve this aim, the gravitational interaction between two extended bodies and more precisely the interaction between mass multipole moments of their gravitational fields and the associated tidal phenomena are studied. Use of Cartesian Symmetric Trace Free (STF) tensors, of their relation with spherical harmonics and of the Kaula's transform enables to derive analytically the tidal and mutual interaction potentials as well as the associated disturbing functions in extended bodies systems. The tidal and mutual interaction potentials of two extended bodies are derived. Next, the external gravitational potential of such tidally disturbed extended body is obtained, using the Love's number theory, as well as the associated disturbing function. Finally, the dynamical evolution equations for such systems are given in their more general form without any linearization. The dynamical equations for the gravitational and tidal interactions between extended bodies and associated dynamics are derived in a form where they could be directly implemented to perform coherent numerical simulations of planetary systems or multiple stars tidal evolution. △ Less

Submitted 12 December, 2008; v1 submitted 12 November, 2007; originally announced November 2007.

Comments: accepted for publication in Astronomy and Astrophysics

On magnetic instabilities and dynamo action in stellar radiation zones

Authors: J. -P. Zahn, A. S. Brun, S. Mathis

Abstract: We examine the MHD instabilities arising in the radiation zone of a differentially rotating star, in which a poloidal field of fossil origin is sheared into a toroidal field. We focus on the non-axisymmetric instability that affects the toroidal magnetic field in a rotating star, which was first studied by Pitts and Tayler in the non-dissipative limit. According to Spruit, it could also drive a… ▽ More We examine the MHD instabilities arising in the radiation zone of a differentially rotating star, in which a poloidal field of fossil origin is sheared into a toroidal field. We focus on the non-axisymmetric instability that affects the toroidal magnetic field in a rotating star, which was first studied by Pitts and Tayler in the non-dissipative limit. According to Spruit, it could also drive a dynamo. The Pitts & Tayler instability is manifestly present in our simulations, with its conspicuous m=1 dependence in azimuth. But its analytic treatment used so far is too simplified to be applied to the real stellar situation. Although the instability generated field reaches an energy comparable to that of the mean poloidal field, that field seems unaffected by the instability: it undergoes Ohmic decline, and is neither eroded nor regenerated by the instability. The toroidal field is produced by shearing the poloidal field and it draws its energy from the differential rotation. The small scale motions behave as Alfven waves; they cause negligible eddy-diffusivity and contribute little to the net transport of angular momentum. In our simulations we observe no sign of dynamo action, of either mean field or fluctuation type, up to a magnetic Reynolds number of 10^5. However the Pitts & Tayler instability is sustained as long as the differential rotation acting on the poloidal field is able to generate a toroidal field of sufficient strength. △ Less

Submitted 23 July, 2007; originally announced July 2007.

Comments: 11 pages, 6 figures, Accepted in Astronomy & Astrophysics

Transport and mixing by internal waves in stellar interiors: effect of the Coriolis force

Authors: S. Mathis, J. -P. Zahn

Abstract: We briefly recall the physical background of the transport of angular momentum and the mixing of chemicals inside stellar radiation zones and its importance for stellar evolution. Then, we describe its present modeling, its successes and its weaknesses. Next, we introduce the new theoretical developments that allow us to treat in a self-consistent way the effect of the Coriolis force on the low-… ▽ More We briefly recall the physical background of the transport of angular momentum and the mixing of chemicals inside stellar radiation zones and its importance for stellar evolution. Then, we describe its present modeling, its successes and its weaknesses. Next, we introduce the new theoretical developments that allow us to treat in a self-consistent way the effect of the Coriolis force on the low-frequencies internal waves and its consequences for the transport processes. This research is aimed at improving the modeling of stellar interiors in the perspective of future astero and helioseismology missions such as COROT and GOLF-NG. △ Less

Submitted 16 June, 2007; originally announced June 2007.

Comments: 4 pages, SF2A-2005: Semaine de l'Astrophysique Francaise, meeting held in Strasbourg, France, June 27 - July 1, 2005, Edited by F. Casoli, T. Contini, J.M. Hameury and L. Pagani. Published by EdP-Sciences, Conference Series, 2005, p. 319

Advances in secular magnetohydrodynamics of stellar interiors dedicated to asteroseismic spatial missions

Authors: S. Mathis, P. Eggenberger, T. Decressin, A. Palacios, L. Siess, C. Charbonnel, S. Turck-Chieze, J. -P. Zahn

Abstract: With the first light of COROT, the preparation of KEPLER and the future helioseismology spatial projects such as GOLF-NG, a coherent picture of the evolution of rotating stars from their birth to their death is needed. We describe here the modelling of the macroscopic transport of angular momentum and matter in stellar interiors that we have undertaken to reach this goal. First, we recall in det… ▽ More With the first light of COROT, the preparation of KEPLER and the future helioseismology spatial projects such as GOLF-NG, a coherent picture of the evolution of rotating stars from their birth to their death is needed. We describe here the modelling of the macroscopic transport of angular momentum and matter in stellar interiors that we have undertaken to reach this goal. First, we recall in detail the dynamical processes that are driving these mechanisms in rotating stars and the theoretical advances we have achieved. Then, we present our new results of numerical simulations which allow us to follow in 2D the secular hydrodynamics of rotating stars, assuming that anisotropic turbulence enforces a shellular rotation law. Finally, we show how this work is leading to a dynamical vision of the Hertzsprung-Russel diagram with the support of asteroseismology and helioseismology, seismic observables giving constraints on the modelling of the internal transport and mixing processes. In conclusion, we present the different processes that should be studied in the next future to improve our description of stellar radiation zones. △ Less

Submitted 6 March, 2007; originally announced March 2007.

Comments: 14 pages, 3 figures, Proceeding of the Joint HELAS and CoRoT/ESTA Workshop (20-23 November 2006, CAUP, Porto - Portugal)

The Frequency of Mid-Infrared Excess Sources in Galactic Surveys

Authors: B. Uzpen, H. A. Kobulnicky, A. J. Monson, M. J. Pierce, D. P. Clemens, D. E. Backman, M. R. Meade, B. L. Babler, R. Indebetouw, B. A. Whitney, C. Watson, M. G. Wolfire, R. A. Benjamin, S. Bracker, T. M. Bania, M. Cohen, C. J. Cyganowski, K. E. Devine, F. Heitsch, J. M. Jackson, J. S. Mathis, E. P. Mercer, M. S. Povich, J. Rho, T. P. Robitaille , et al. (5 additional authors not shown)

Abstract: We have identified 230 Tycho-2 Spectral Catalog stars that exhibit 8 micron mid-infrared extraphotospheric excesses in the MidCourse Space Experiment (MSX) and Spitzer Space Telescope Galactic Legacy MidPlane Survey Extraordinaire (GLIMPSE) surveys. Of these, 183 are either OB stars earlier than B8 in which the excess plausibly arises from a thermal bremsstrahlung component or evolved stars in w… ▽ More We have identified 230 Tycho-2 Spectral Catalog stars that exhibit 8 micron mid-infrared extraphotospheric excesses in the MidCourse Space Experiment (MSX) and Spitzer Space Telescope Galactic Legacy MidPlane Survey Extraordinaire (GLIMPSE) surveys. Of these, 183 are either OB stars earlier than B8 in which the excess plausibly arises from a thermal bremsstrahlung component or evolved stars in which the excess may be explained by an atmospheric dust component. The remaining 47 stars have spectral classifications B8 or later and appear to be main sequence or late pre-main-sequence objects harboring circumstellar disks. Six of the 47 stars exhibit multiple signatures characteristic of pre-main-sequence circumstellar disks, including emission lines, near-infrared K-band excesses, and X-ray emission. Approximately one-third of the remaining 41 sources have emission lines suggesting relative youth. Of the 25 GLIMPSE stars with SST data at >24 microns, 20 also show an excess at 24 microns. Three additional objects have 24 micron upper limits consistent with possible excesses, and two objects have photospheric measurements at 24 microns. Six MSX sources had a measurement at wavelengths >8 microns. We modeled the excesses in 26 stars having two or more measurements in excess of the expected photospheres as single-component blackbodies. We determine probable disk temperatures and fractional infrared luminosities in the range 191 < T < 787 and 3.9x10^-4 < L_IR/L_* < 2.7x10^-1. We estimate a lower limit on the fraction of Tycho-2 Spectral Catalog main-sequence stars having mid-IR, but not near-IR, excesses to be 1.0+-0.3%. △ Less

Submitted 8 December, 2006; originally announced December 2006.

Comments: Accepted to ApJ

Journal ref: Astrophys.J.658:1264-1288,2007

Absolute diffuse calibration of IRAC through mid-infrared and radio study of HII regions

Authors: Martin Cohen, Anne J. Green, Marilyn R. Meade, Brian Babler, Remy Indebetouw, Barbara A. Whitney, Christer Watson, Mark Wolfire, Mike J. Wolff, John S. Mathis, Edward B. Churchwell, .

Abstract: We investigate the diffuse absolute calibration of the InfraRed Array Camera on the Spitzer Space Telescope at 8.0microns using a sample of 43 HII regions with a wide range of morphologies near GLON=312deg. For each region we carefully measure sky-subtracted,point-source- subtracted, areally-integrated IRAC 8.0-micron fluxes and compare these with Midcourse Space eXperiment (MSX) 8.3-micron imag… ▽ More We investigate the diffuse absolute calibration of the InfraRed Array Camera on the Spitzer Space Telescope at 8.0microns using a sample of 43 HII regions with a wide range of morphologies near GLON=312deg. For each region we carefully measure sky-subtracted,point-source- subtracted, areally-integrated IRAC 8.0-micron fluxes and compare these with Midcourse Space eXperiment (MSX) 8.3-micron images at two different spatial resolutions, and with radio continuum maps. We determine an accurate median ratio of IRAC 8.0-micron/MSX\8.3-micron fluxes, of 1.55+/-0.15. From robust spectral energy distributions of these regions we conclude that the present 8.0-micron diffuse calibration of the SST is 36% too high compared with the MSX validated calibration, perhaps due to scattered light inside the camera. This is an independent confirmation of the result derived for the diffuse calibration of IRAC by the Spitzer Science Center (SSC). From regression analyses we find that 843-MHz radio fluxes of HII regions and mid-infrared (MIR) fluxes are linearly related for MSX at 8.3-microns and Spitzer at 8.0 microns, confirming the earlier MSX result by Cohen & Green. The median ratio of MIR/843-MHz diffuse continuum fluxes is 600 times smaller in nonthermal than thermal regions, making it a sharp discriminant. The ratios are largely independent of morphology up to a size of ~24 arcsec. We provide homogeneous radio and MIR morphologies for all sources. MIR morphology is not uniquely related to radio structure. Compact regions may have MIR filaments and/or diffuse haloes, perhaps infrared counter- parts to weakly ionized radio haloes found around compact HII regions. We offer two IRAC colour-colour plots as quantitative diagnostics of diffuse HII regions. △ Less

Submitted 19 October, 2006; originally announced October 2006.

Comments: 29 pages, LaTeX (aastex), incl. 31 PostScript (ps,eps) figures and 5 tables. Accepted by MNRAS (main journal). Replaced an unused file and added this URL for people wishing to download a version with high-resolution images: http://www.astro.wisc.edu/sirtf/martin.hii.accepted.pdf

Journal ref: Mon.Not.Roy.Astron.Soc.374:979-998,2007

The magnetism of the solar interior for a complete MHD solar vision

Authors: S. Turck-Chieze, T. Appourchaux, J. Ballot, G. Berthomieu, P. Boumier, A. S. Brun, A. Cacciani, J. Christensen-Dalsgaard, T. Corbard, S. Couvidat, A. Eff-Darwich, B. Dintrans, E. Fossat, R. A. Garcia, B. Gelly, L. Gizon, D. Gough, A. Jimenez, S. Jimenez-Reyes, A. Kosovishev, P. Lambert, I. Lopes, M. Martic, S. Mathis, N. Meunier , et al. (15 additional authors not shown)

Abstract: The solar magnetism is no more considered as a purely superficial phenomenon. The SoHO community has shown that the length of the solar cycle depends on the transition region between radiation and convection. Nevertheless, the internal solar (stellar) magnetism stays poorly known. Starting in 2008, the American instrument HMI/SDO and the European microsatellite PICARD will enrich our view of the… ▽ More The solar magnetism is no more considered as a purely superficial phenomenon. The SoHO community has shown that the length of the solar cycle depends on the transition region between radiation and convection. Nevertheless, the internal solar (stellar) magnetism stays poorly known. Starting in 2008, the American instrument HMI/SDO and the European microsatellite PICARD will enrich our view of the Sun-Earth relationship. Thus obtaining a complete MHD solar picture is a clear objective for the next decades and it requires complementary observations of the dynamics of the radiative zone. For that ambitious goal, space prototypes are being developed to improve gravity mode detection. The Sun is unique to progress on the topology of deep internal magnetic fields and to understand the complex mechanisms which provoke photospheric and coronal magnetic changes and possible longer cycles important for human life. We propose the following roadmap in Europe to contribute to this "impressive" revolution in Astronomy and in our Sun-Earth relationship: SoHO (1995-2007), PICARD (2008-2010), DynaMICS (2009-2017) in parallel to SDO (2008-2017) then a world-class mission located at the L1 orbit or above the solar poles. △ Less

Submitted 31 October, 2005; originally announced October 2005.

Comments: 10 pages, 8 figures, will appear in the proceedings of 2005 ESLAB symposium 19-21 April 2005

First GLIMPSE Results on the Stellar Structure of the Galaxy

Authors: Robert A. Benjamin, E. Churchwell, B. L. Babler, R. Indebetouw, M. R. Meade, B. A. Whitney, C. Watson, M. G. Wolfire, M. J. Wolff, R. Ignace, T. M. Bania, S. Bracker, D. P. Clemens, L. Chomiuk, M. Cohen, J. M. Dickey, J. M. Jackson, H. A. Kobulnicky, E. P. Mercer, J. S. Mathis, S. R. Stolovy, B. Uzpen

Abstract: The GLIMPSE (Galactic Legacy Mid-Plane Survey Extraordinaire) Point Source Catalog of ~ 30 million mid-infrared sources towards the inner Galaxy, 10 < |l| < 65 degrees and |b| < 1 degree, was used to determine the distribution of stars in Galactic longitude, latitude, and apparent magnitude. The counts versus longitude can be approximated by the modified Bessel function N=N_0*(l/l_0)*K_1(l/l_0),… ▽ More The GLIMPSE (Galactic Legacy Mid-Plane Survey Extraordinaire) Point Source Catalog of ~ 30 million mid-infrared sources towards the inner Galaxy, 10 < |l| < 65 degrees and |b| < 1 degree, was used to determine the distribution of stars in Galactic longitude, latitude, and apparent magnitude. The counts versus longitude can be approximated by the modified Bessel function N=N_0*(l/l_0)*K_1(l/l_0), where l_0 is insensitive to limiting magnitude, band choice, and side of Galactic center: l_0= 17-30 degrees with a best fit value in the the 4.5 micron band of l_0=24 +/- 4 degrees. Modeling the source distribution as an exponential disk yields a radial scale length of H= 3.9 +/- 0.6 kpc. There is a pronounced north-south asymmetry in source counts for |l| < 30 degrees, with ~ 25% more stars in the north. For l=10-30 degrees, there is a strong enhancement of stars of m= 11.5-13.5 mag. A linear bar passing through the Galactic center with half-length R_bar=4.4 +/- 0.5 kpc, tilted by phi=44 +/- 10 degrees to the Sun-Galactic Center line, provides the simplest interpretation of this data. We examine the possibility that enhanced source counts at l=26-28 degrees, 31.5-34 degrees, and 306-309 degrees are related to Galactic spiral structure. Total source counts are depressed in regions where the counts of red objects (m_K-m_[8.0] >3) peak. In these areas, the counts are reduced by extinction due to molecular gas and/or high diffuse backgrounds associated with star formation. △ Less

Submitted 15 August, 2005; originally announced August 2005.

Comments: Accepted for publication in Astrophysical Journal Letters. Version with full res figures and info on GLIMPSE Point Source Catalog at http://www.astro.wisc.edu/glimpse

Journal ref: Astrophys.J.630:L149-L152,2005

Transport and mixing in the radiation zones of rotating stars II. Axisymmetric magnetic field

Authors: S. Mathis, J. -P. Zahn

Abstract: The purpose of this paper is to improve the modeling of the mixing of chemical elements that occurs in stellar radiation zones. In addition to the classical rotational mixing considered in our previous paper, which results of the combined action of the thermally-driven meridional circulation and of the turbulence generated by the shear of differential rotation, we include here the effect of an a… ▽ More The purpose of this paper is to improve the modeling of the mixing of chemical elements that occurs in stellar radiation zones. In addition to the classical rotational mixing considered in our previous paper, which results of the combined action of the thermally-driven meridional circulation and of the turbulence generated by the shear of differential rotation, we include here the effect of an axisymmetric magnetic field in a self-consistent way. We treat the advection of the field by the meridional circulation, its Ohmic diffusion, and the production of its toroidal component through the shear of differential rotation. The Lorentz force is assumed not to exceed the centrifugal force; it acts on the baroclinic balance and therefore on the meridional flow, and it has a strong impact on the transport of angular momentum. All variables and governing equations are expanded in spherical or spherical vectorial functions, to arbitrary order: this yields a system of partial differential equations in time and in the radial coordinate, which is ready to be implemented in a stellar structure code. △ Less

Submitted 6 June, 2005; originally announced June 2005.

Comments: 28 pages, 1 figure, accepted for publication in Astronomy and Astrophysics

FUSE Measurements of Far Ultraviolet Extinction. II. Magellanic Cloud Sight Lines

Authors: Stefan I. B. Cartledge, Geoffrey C. Clayton, Karl D. Gordon, Brian L. Rachford, B. T. Draine, P. G. Martin, John S. Mathis, K. A. Misselt, Ulysses J. Sofia, D. C. B. Whittet, Michael J. Wolff

Abstract: We present an extinction analysis of 9 paths through the LMC and SMC based on FUSE observations. To date, just two LMC sight lines have probed dust grain composition and size distributions in the Clouds using spectra including wavelengths as short as 950 A. We supplement these with results from 4 regions distinguished by their IR through UV extinction curves and grouped as LMCAvg, LMC2, SMC bar… ▽ More We present an extinction analysis of 9 paths through the LMC and SMC based on FUSE observations. To date, just two LMC sight lines have probed dust grain composition and size distributions in the Clouds using spectra including wavelengths as short as 950 A. We supplement these with results from 4 regions distinguished by their IR through UV extinction curves and grouped as LMCAvg, LMC2, SMC bar and SMC wing. Despite the distinct characters of Milky Way and Magellanic Cloud extinction, our results are generally analogous to those found for Galactic curves in that the FUSE portions of each extinction curve are described reasonably well by FM curves fitted only to longer wavelength data and lack any dramatic new extinction features, and any deviations from the CCM formalism continue into FUV wavelengths. An MEM analysis of these curves suggests that LMCAvg and SMC wing sight lines require more silicon and/or carbon in dust than current abundance measurements would allow, while the requirements for LMC2 and SMC bar sight lines do not fully tax the available reservoirs. An intermediate product of this analysis is the measurement of new H_2 abundances in the Magellanic Clouds. Collectively considering Cloud sight lines that possess significant H_2 column densities, E(B-V)/N(HI) ratios are reduced by significant factors relative to the Galactic mean, whereas the corresponding E(B-V)/N(H_2) values more closely resemble their Galactic counterpart. These trends reflect the fact that among these sight lines f(H_2)-values are lower than those common in the Milky Way for paths with similar degrees of reddening. △ Less

Submitted 19 May, 2005; originally announced May 2005.

Comments: 40 pages, 8 figures, accepted for publication in ApJ

Identification of Main Sequence Stars with Mid-Infrared Excesses Using GLIMPSE: Beta-Pictoris Analogs?

Authors: B. Uzpen, H. A. Kobulnicky, K. A. G. Olsen, D. P. Clemens, T. L. Laurance, M. R. Meade, B. L. Babler, R. Indebetouw, B. A. Whitney, C. Watson, M. G. Wolfire, M. J. Wolff, R. A. Benjamin, T. M. Bania, M. Cohen, K. E. Devine, J. M. Dickey, F. Heitsch, J. M. Jackson, A. P. Marston, J. S. Mathis, E. P. Mercer, J. R. Stauffer, S. R. Stolovy, D. E. Backman , et al. (1 additional authors not shown)

Abstract: Spitzer IRAC 3.6-8 micron photometry obtained as part of the GLIMPSE survey has revealed mid-infrared excesses for 33 field stars with known spectral types in a 1.2 sq. degree field centered on the southern Galactic HII region RCW49. These stars comprise a subset of 184 stars with known spectral classification, most of which were pre-selected to have unusually red IR colors. We propose that the… ▽ More Spitzer IRAC 3.6-8 micron photometry obtained as part of the GLIMPSE survey has revealed mid-infrared excesses for 33 field stars with known spectral types in a 1.2 sq. degree field centered on the southern Galactic HII region RCW49. These stars comprise a subset of 184 stars with known spectral classification, most of which were pre-selected to have unusually red IR colors. We propose that the mid-IR excesses are caused by circumstellar dust disks that are either very late remnants of stellar formation or debris disks generated by planet formation. Of these 33 stars, 29 appear to be main-sequence stars based on optical spectral classifications. Five of the 29 main-sequence stars are O or B stars with excesses that can be plausibly explained by thermal bremsstrahlung emission, and four are post main-sequence stars. The lone O star is an O4V((f)) at a spectrophotometric distance of 3233+ 540- 535 pc and may be the earliest member of the Westerlund 2 cluster. Of the remaining 24 main-sequence stars, 18 have SEDs that are consistent with hot dusty debris disks, a possible signature of planet formation. Modeling the excesses as blackbodies demonstrates that the blackbody components have fractional bolometric disk-to-star luminosity ratios, LIR/L*, ranging from 10^-3 to 10^-2 with temperatures ranging from 220 to 820 K. The inferred temperatures are more consistent with asteroid belts rather than the cooler temperatures expected for Kuiper belts. Mid-IR excesses are found in all spectral types from late B to early K. △ Less

Submitted 20 April, 2005; originally announced April 2005.

Journal ref: Astrophys.J. 629 (2005) 512-525

The Effects of Clumping on Derived Abundances in H II Regions

Authors: John S. Mathis, Kenneth Wood

Abstract: We have compared Monte Carlo photoionization models of H II regions with a uniform density distribution with models with the same central stars and chemical compositions but with 3-D hierarchical clumps. We compare the abundances of He, N, O, Ne, and S obtained from emission line strengths and [O III] and [N II] temperatures to those in our models. We consider stellar temperatures in the range 3… ▽ More We have compared Monte Carlo photoionization models of H II regions with a uniform density distribution with models with the same central stars and chemical compositions but with 3-D hierarchical clumps. We compare the abundances of He, N, O, Ne, and S obtained from emission line strengths and [O III] and [N II] temperatures to those in our models. We consider stellar temperatures in the range 37.5 -- 45kK and ionizing luminosities from 10^{48} to 10^{51} photons/s. Clumped models have different ionic abundances than uniform. For hot stars, He^0/He^+ is 2 -- 3%, much larger than with uniform models. This amount of He I is independent of metallicity and so impacts the determination of the primordial abundance of He. The total abundances of O, Ne, and S obtained by the usual methods of analysis, using T([OIII) for high stages of ionization and T([NII]) for low, are about as accurate for clumped models as for uniform and within about 20% of the true values. If T([OIII]) is used for analyzing all ions, the derived (O/H) is 40 to 60% too large for cool stars but is good for hot stars. Uniform models have similar errors, so the clumping does not change the accuracy of abundance analysis. The physical causes of the ionic abundance errors are present in real nebulae. In clumped models, helium ionizing radiation from zones of high ionization (low He^0 and low UV opacity) can penetrate nearby regions near the edge of the ionized zone. This effect allows He^0 to absorb more stellar photons than in uniform or radially symmetrical geometries. In turn, these absorptions compete with O+, etc., for those energetic stellar photons. △ Less

Submitted 26 March, 2005; originally announced March 2005.

Comments: 11 pages, 10 figures. Accepted by MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc. 360 (2005) 227-235

G313.3+00.3: A New Planetary Nebula discovered by the Australia Telescope Compact Array and the Spitzer Space Telescope

Authors: Martin Cohen, Anne J. Green, Mallory S. E. Roberts, Marilyn R. Meade, Brian Babler, Remy Indebetouw, Barbara A. Whitney, Christer Watson, Mark Wolfire, Mike J. Wolff, John S. Mathis, Edward B. Churchwell

Abstract: We present a new planetary nebula, first identified in images from the Australia Telescope Compact Array, although not recognized at that time. Recent observations with the Spitzer Space Telescope during the GLIMPSE Legacy program have rediscovered the object. The high-resolution radio and infrared images enable the identification of the central star or its wind, the recognition of the radio emi… ▽ More We present a new planetary nebula, first identified in images from the Australia Telescope Compact Array, although not recognized at that time. Recent observations with the Spitzer Space Telescope during the GLIMPSE Legacy program have rediscovered the object. The high-resolution radio and infrared images enable the identification of the central star or its wind, the recognition of the radio emission as thermal, and the probable presence of polycylic aromatic hydrocarbons in and around the source. These lead to the conclusion that G313.3+00.3 is a planetary nebula. This object is of particular interest because it was discovered solely through radio and mid-infrared imaging, without any optical (or near-infrared) confirmation, and acts as a proof of concept for the discovery of many more highly extinguished planetary nebulae. G313.3+00.3 is well-resolved by both the instruments with which it was identified, and suffers extreme reddening due to its location in the Scutum-Crux spiral arm. △ Less

Submitted 15 March, 2005; originally announced March 2005.

Comments: 18 pages, LaTeX (aastex), incl. 8 PostScript (eps) figures and 1 table. Accepted by ApJ (Part 1)

Journal ref: Astrophys.J. 627 (2005) 446-453

FUSE Measurements of Far Ultraviolet Extinction. I. Galactic Sight Lines

Authors: U. J. Sofia, M. J. Wolff, B. Rachford, K. D. Gordon, G. C. Clayton, S. I. B. Cartledge, P. G. Martin, B. T. Draine, J. S. Mathis, T. P. Snow, D. C. B. Whittet

Abstract: We present extinction curves that include data down to far ultraviolet wavelengths (FUV; 1050 - 1200 A) for nine Galactic sight lines. The FUV extinction was measured using data from the Far Ultraviolet Spectroscopic Explorer. The sight lines were chosen for their unusual extinction properties in the infrared through the ultraviolet; that they probe a wide range of dust environments is evidenced… ▽ More We present extinction curves that include data down to far ultraviolet wavelengths (FUV; 1050 - 1200 A) for nine Galactic sight lines. The FUV extinction was measured using data from the Far Ultraviolet Spectroscopic Explorer. The sight lines were chosen for their unusual extinction properties in the infrared through the ultraviolet; that they probe a wide range of dust environments is evidenced by the large spread in their measured ratios of total-to-selective extinction, R_V = 2.43 - 3.81. We find that extrapolation of the Fitzpatrick & Massa relationship from the ultraviolet appears to be a good predictor of the FUV extinction behavior. We find that predictions of the FUV extinction based upon the Cardelli, Clayton & Mathis (CCM) dependence on R_V give mixed results. For the seven extinction curves well represented by CCM in the infrared through ultraviolet, the FUV extinction is well predicted in three sight lines, over-predicted in two sight lines, and under-predicted in 2 sight lines. A Maximum Entropy Method analysis using a simple three component grain model shows that seven of the nine sight lines in the study require a larger fraction of grain materials to be in dust when FUV extinction is included in the models. Most of the added grain material is in the form of small (radii < 200 A) grains. △ Less

Submitted 24 February, 2005; originally announced February 2005.

Comments: Accepted for publication in the Astrophysical Journal. 31 pages with 7 figures

Journal ref: Astrophys.J. 625 (2005) 167-180

Discovery of a New Low-Latitude Milky Way Globular Cluster using GLIMPSE

Authors: Chip Kobulnicky, A. J. Monson, B. A Buckalew, J. M. Darnel, B. Uzpen, B. A. Whitney, R. Indebetouw, B. L. Babler, M. R. Meade, C. Watson, E. Churchwell, M. J. Wolff, M. G. Wolfire, D. P. Clemens, R. Shah, T. M. Bania, R. A. Benjamin, M. Cohen, K. E. Devine, J. M. Dickey, F. Heitsch, J. M. Jackson, A. P. Marston, J. S. Mathis, E. P. Mercer , et al. (6 additional authors not shown)

Abstract: Spitzer Space Telescope imaging as part of the Galactic Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) reveals a previously unidentified low-latitude rich star cluster near l=31.3 degrees, b=-0.1 degrees. Near-infrared JHK' photometry from the Wyoming Infrared Observatory indicates an extinction of A_V ~ 15+/-3 mag for cluster members. Analysis of 13CO features along the same sightline suggest… ▽ More Spitzer Space Telescope imaging as part of the Galactic Legacy Mid-Plane Survey Extraordinaire (GLIMPSE) reveals a previously unidentified low-latitude rich star cluster near l=31.3 degrees, b=-0.1 degrees. Near-infrared JHK' photometry from the Wyoming Infrared Observatory indicates an extinction of A_V ~ 15+/-3 mag for cluster members. Analysis of 13CO features along the same sightline suggests a probable kinematic distance of 3.1 -- 5.2 kpc. The new cluster has an angular diameter of ~1-2 pc, a total magnitude m_{K_0}=2.1 corrected for extinction, and a luminosity of M_K ~ -10.3 at 3.1 kpc. In contrast to young massive Galactic clusters with ages <100 Myr, the new cluster has no significant radio emission. Comparison to theoretical K-band luminosity functions indicates an age of at least several Gyr and a mass of at least 10^5 solar masses. Unlike known old open clusters, this new cluster lies in the inner Galaxy at R_{GC} ~ 6.1 kpc. We designate this object ``GLIMPSE-C01'' and classify it as a Milky Way globular cluster passing through the Galactic disk. We also identify a region of star formation and fan-shaped outflows from young stellar objects in the same field as the cluster. The cluster's passage through the Galactic molecular layer may have triggered this star formation activity. △ Less

Submitted 16 October, 2004; originally announced October 2004.

Comments: Accepted for Publication in AJ Report with full resolution images available at http://physics.uwyo.edu/~chip/Glob

Journal ref: Astron.J. 129 (2005) 239-250

Transport and mixing in the radiation zones of rotating stars: I-Hydrodynamical processes

Authors: S. Mathis, J. -P. Zahn

Abstract: The purpose of this paper is to improve the modelization of the rotational mixing which occurs in stellar radiation zones, through the combined action of the thermally driven meridional circulation and of the turbulence generated by the shear of differential rotation. The turbulence is assumed to be anisotropic, due to the stratification, with stronger transport in the horizontal directions than… ▽ More The purpose of this paper is to improve the modelization of the rotational mixing which occurs in stellar radiation zones, through the combined action of the thermally driven meridional circulation and of the turbulence generated by the shear of differential rotation. The turbulence is assumed to be anisotropic, due to the stratification, with stronger transport in the horizontal directions than in the vertical. The main difference with the former treatments by Zahn (1992) and Maeder & Zahn (1998) is that we expand here the departures from spherical symmetry to higher order, and include explicitly the differential rotation in latitude, to first order. This allows us to treat simultaneously the bulk of a radiation zone and its tachocline(s). Moreover, we take fully into account the non-stationarity of the problem, which will enable us to tackle the rapid phases of evolution. The system of partial differential equations, which govern the transport of angular momentum, heat and chemical elements, is written in a form which makes it ready to implement in a stellar evolution code. Here the effect of a magnetic field is deliberately ignored; it will be included in forthcoming papers. △ Less

Submitted 18 June, 2004; originally announced June 2004.

Comments: 16 pages, no figures, accepted for publication in A&A

Journal ref: Astron.Astrophys. 425 (2004) 229-242

The Wavelength Dependence of Interstellar Extinction from 1.25 to 8.0 microns Using GLIMPSE Data

Authors: R. Indebetouw, J. S. Mathis, B. L. Babler, M. R. Meade, C. Watson, B. A. Whitney, M. J. Wolff, M. G. Wolfire, T. M. Bania, R. A. Benjamin, D. P. Clemens, Martin Cohen, J. M. Dickey, J. M. Jackson, H. A. Kobulnicky, A. P. Marston, E. P. Mercer, J. R. Stauffer, S. R. Stolovy, E. Churchwell

Abstract: We determine and tabulate A(lambda)/A(K), the wavelength dependence of interstellar extinction, in the Galactic plane for 1.25um<lambda<8.0um along two lines of sight: l~42 and l~284. The first is a relatively quiescent and unremarkable region; the second contains the giant HII region RCW 49 as well as a ``field'' region unrelated to the cluster and nebulosity. Areas near these Galactic longitud… ▽ More We determine and tabulate A(lambda)/A(K), the wavelength dependence of interstellar extinction, in the Galactic plane for 1.25um<lambda<8.0um along two lines of sight: l~42 and l~284. The first is a relatively quiescent and unremarkable region; the second contains the giant HII region RCW 49 as well as a ``field'' region unrelated to the cluster and nebulosity. Areas near these Galactic longitudes were imaged at J, H, and K bands by 2MASS and at 3--8um by Spitzer for the GLIMPSE Legacy program. We measure the mean values of the color excess ratios (A(lambda)-A(K))/(A(J)-A(K)) directly from the color distributions of observed stars. The extinction ratio between two of the filters, e.g. A(J)/A(K), is required to calculate A(lambda)/A(K) from those measured ratios. We use the apparent JHK magnitudes of giant stars along our two sightlines, and fit the reddening as a function of magnitude (distance) to determine A(J)/kpc, A(K)/kpc, and A(J)/A(K). Our values of A(lambda)/A(K) show a flattening across the 3--8um wavelength range, roughly consistent with the Lutz 1996 extinction measurements derived for the sightline toward the Galactic center. △ Less

Submitted 28 October, 2004; v1 submitted 17 June, 2004; originally announced June 2004.

Comments: 20pp, 6fig, accepted version: description of the method has been revised, results are unchanged, discussion of implications slightly expanded

Journal ref: Astrophys.J.619:931-938,2005

A GLIMPSE of Star Formation in the Giant H II Region RCW 49

Authors: B. A. Whitney, R. Indebetouw, B. L. Babler, M. R. Meade, C. Watson, M. J. Wolff, M. G. Wolfire, D. P. Clemens, T. M. Bania, R. A. Benjamin, M. Cohen, K. E. Devine, J. M. Dickey, F. Heitsch, J. M. Jackson, H. A. Kobulnicky, A. P. Marston, J. S. Mathis, E. P. Mercer, J. R. Stauffer, S. R. Stolovy, E. Churchwell

Abstract: GLIMPSE imaging using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope indicates that star formation is ongoing in the RCW 49 giant H II region. A photometric comparison of the sources in RCW 49 to a similar area to its north finds that at least 300 stars brighter than 13th magnitude in band [3.6] have infrared excesses inconsistent with reddening due to foreground extinction. The… ▽ More GLIMPSE imaging using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope indicates that star formation is ongoing in the RCW 49 giant H II region. A photometric comparison of the sources in RCW 49 to a similar area to its north finds that at least 300 stars brighter than 13th magnitude in band [3.6] have infrared excesses inconsistent with reddening due to foreground extinction. These are likely young stellar objects (YSOs) more massive than 2.5 M_sun suggesting that thousands more low-mass stars are forming in this cloud. Some of the YSOs are massive (B stars) and therefore very young, suggesting that a new generation of star formation is occurring, possibly triggered by stellar winds and shocks generated by the older (2-3 Myr) central massive cluster. The Spitzer IRAC camera has proven to be ideally suited for distinguishing young stars from field stars, and the GLIMPSE survey of the Galactic Plane will likely find thousands of new star formation regions. △ Less

Submitted 3 June, 2004; originally announced June 2004.

Comments: Accepted by ApJS, Spitzer edition. 9 pages, 6 figures. Due to space limitations on astro-ph, figures 1-4 are submitted as separate .jpg files at http://arxiv.org/ps/astro-ph/0406100. Full-resolution figures and text can be downloaded from http://www.astro.wisc.edu/glimpse/glimpsepubs.html

Journal ref: Astrophys.J.Suppl. 154 (2004) 315-321

Dust heating by the interstellar radiation field in models of turbulent molecular clouds

Authors: Thomas Bethell, Ellen Zweibel, Fabian Heitsch, J. S. Mathis

Abstract: We have calculated the radiation field, dust grain temperatures, and far infrared emissivity of numerical models of turbulent molecular clouds. When compared to a uniform cloud of the same mean optical depth, most of the volume inside the turbulent cloud is brighter, but most of the mass is darker. There is little mean attenuation from center to edge, and clumping causes the radiation field to b… ▽ More We have calculated the radiation field, dust grain temperatures, and far infrared emissivity of numerical models of turbulent molecular clouds. When compared to a uniform cloud of the same mean optical depth, most of the volume inside the turbulent cloud is brighter, but most of the mass is darker. There is little mean attenuation from center to edge, and clumping causes the radiation field to be somewhat bluer. There is also a large dispersion, typically by a few orders of magnitude, of all quantities relative to their means. However, despite the scatter, the 850 micron emission maps are well correlated with surface density. The fraction of mass as a function of intensity can be reproduced by a simple hierarchical model of density structure. △ Less

Submitted 2 April, 2004; originally announced April 2004.

Comments: 32 pages, 14 figures, submitted to ApJ

Journal ref: Astrophys.J. 610 (2004) 801-812

On shear-induced turbulence in rotating stars

Authors: S. Mathis, A. Palacios, J. -P. Zahn

Abstract: We review various prescriptions which have been proposed for the turbulent transport of matter and angular momentum in differentially rotating stellar radiation zones. A new prescription is presented for the horizontal transport associated with the anisotropic shear turbulence which is produced by the differential rotation in latitude; this `beta-viscosity' is drawn from torque measurements in t… ▽ More We review various prescriptions which have been proposed for the turbulent transport of matter and angular momentum in differentially rotating stellar radiation zones. A new prescription is presented for the horizontal transport associated with the anisotropic shear turbulence which is produced by the differential rotation in latitude; this `beta-viscosity' is drawn from torque measurements in the classical Couette-Taylor experiment (Richard & Zahn 1999). Its implementation in a stellar evolution code leads to enhanced mixing, as illustrated by models of a rotating main-sequence star of 1.5 solar mass. △ Less

Submitted 24 March, 2004; originally announced March 2004.

Comments: 5 pages, 3 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: Astron.Astrophys. 425 (2004) 243-247

Monte Carlo Photoionization Simulations of Diffuse Ionized Gas

Authors: Kenneth Wood, John S. Mathis

Abstract: We demonstrate that the observed increase of some nebular line ratios with height above the midplane in the diffuse ionized gas (DIG) in the Milky Way and other galaxies is a natural consequence of the progressive hardening of the radiation field far from the midplane ionizing sources. To obtain increasing temperatures and line ratios away from the midplane, our photoionization simulations of a… ▽ More We demonstrate that the observed increase of some nebular line ratios with height above the midplane in the diffuse ionized gas (DIG) in the Milky Way and other galaxies is a natural consequence of the progressive hardening of the radiation field far from the midplane ionizing sources. To obtain increasing temperatures and line ratios away from the midplane, our photoionization simulations of a multi-component interstellar medium do not require as much additional heating (over and above that from photoionization) as previous studies that employed one dimensional, spherically averaged models. Radiation leaking into the DIG from density bounded H II regions is generally harder in the H-ionizing continuum and has its He-ionizing photons suppressed compared to the ionizing source of the H II region. In line with other recent investigations, we find that such leaky H II region models can provide elevated temperatures and line ratios, and a lower He+ fraction in the DIG. For a composite model representing the relative spectral types of O stars in the solar neighbourhood, we find that additional heating less than 1.0E-26 n_e ergs/s/cm^3 can reproduce the observed elevated line ratios in the DIG. This additional heating is considerably less than previous estimates due to the natural hardening of the radiation field reaching large heights in our simulations. △ Less

Submitted 17 February, 2004; originally announced February 2004.

Comments: Accepted for publication in MNRAS

GLIMPSE: I. A SIRTF Legacy Project to Map the Inner Galaxy

Authors: Robert A. Benjamin, E. Churchwell, Brian L. Babler, T. M. Bania, Dan P. Clemens, Martin Cohen, John M. Dickey, Remy Indebetouw, James M. Jackson, Henry A. Kobulnicky, Alex Lazarian, A. P. Marston, John S. Mathis, Marilyn R. Meade, Sara Seager, S. R. Stolovy, C. Watson, Barbara A. Whitney, Michael J. Wolff, Mark G. Wolfire

Abstract: GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire), a SIRTF Legacy Science Program, will be a fully sampled, confusion-limited infrared survey of the inner two-thirds of the Galactic disk with a pixel resolution of \~1.2" using the Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 microns. The survey will cover Galactic latitudes |b| <1 degree and longitudes |l|=10 to 65 deg… ▽ More GLIMPSE (Galactic Legacy Infrared Mid-Plane Survey Extraordinaire), a SIRTF Legacy Science Program, will be a fully sampled, confusion-limited infrared survey of the inner two-thirds of the Galactic disk with a pixel resolution of \~1.2" using the Infrared Array Camera (IRAC) at 3.6, 4.5, 5.8, and 8.0 microns. The survey will cover Galactic latitudes |b| <1 degree and longitudes |l|=10 to 65 degrees (both sides of the Galactic center). The survey area contains the outer ends of the Galactic bar, the Galactic molecular ring, and the inner spiral arms. The GLIMPSE team will process these data to produce a point source catalog, a point source data archive, and a set of mosaicked images. We summarize our observing strategy, give details of our data products, and summarize some of the principal science questions that will be addressed using GLIMPSE data. Up-to-date documentation, survey progress, and information on complementary datasets are available on the GLIMPSE web site: www.astro.wisc.edu/glimpse. △ Less

Submitted 12 June, 2003; originally announced June 2003.

Comments: Description of GLIMPSE, a SIRTF Legacy project (Aug 2003 PASP, in press). Paper with full res.color figures at http://www.astro.wisc.edu/glimpse/glimpsepubs.html

Journal ref: Publ.Astron.Soc.Pac.115:953-964,2003

Can Reflection from Grains Diagnose the Albedo?

Authors: John S. Mathis, Barbara A. Whitney, Kenneth Wood

Abstract: By radiation transfer models with a realistic power spectra of the projected density distributions, we show that the optical properties of grains are poorly constrained by observations of reflection nebulae. The ISM is known to be hierarchically clumped from a variety of observations (molecules, H I, far-infrared). Our models assume the albedo and phase parameter of the dust, the radial optical… ▽ More By radiation transfer models with a realistic power spectra of the projected density distributions, we show that the optical properties of grains are poorly constrained by observations of reflection nebulae. The ISM is known to be hierarchically clumped from a variety of observations (molecules, H I, far-infrared). Our models assume the albedo and phase parameter of the dust, the radial optical depth of the sphere averaged over all directions, and random distributions of the dust within the sphere. The outputs are the stellar extinction, optical depth, and flux of scattered light as seen from various viewing angles. Observations provide the extinction and scattered flux from a particular direction. Hierarchical geometry has a large effect on the flux of scattered light emerging from a nebula for a particular extinction of the exciting star. There is a very large spread in both scattered fluxes and extinctions for any distribution of dust. Consequently, an observed stellar extinction and scattered flux can be fitted by a wide range of albedos. With hierarchical geometry it is not completely safe to determine even relative optical constants from multiwavelength observations of the same reflection nebula. The geometry effectively changes with wavelength as the opacity of the clumps varies. Limits on the implications of observing the same object in various wavelengths are discussed briefly. Henry (2002) uses a recipe to determine the scattered flux from a star with a given extinction. It is claimed to be independent of the geometry. It provides considerably more scattering than our models, probably leading to an underestimate of the grain albedos from the UV Diffuse Galactic Light. △ Less

Submitted 2 April, 2002; originally announced April 2002.

Comments: 27 pages, including 7 figures. Accepted by ApJ

ISO spectroscopy of compact HII regions in the Galaxy. II Ionization and elemental abundances

Authors: N. L. Martin-Hernandez, E. Peeters, C. Morisset, A. G. G. M. Tielens, P. Cox, P. R. Roelfsema, J. -P. Baluteau, D. Schaerer, J. S. Mathis, F. Damour, E. Churchwell, M. F. Kessler

Abstract: Based on the ISO spectral catalogue of compact HII regions by Peeters et al. (2001), we present a first analysis of the hydrogen recombination and atomic fine-structure lines originated in the ionized gas. The sample consists of 34 HII regions located at galactocentric distances between Rgal = 0 and 15 kpc. The SWS HI recombination lines between 2 and 8 mum are used to estimate the extinction la… ▽ More Based on the ISO spectral catalogue of compact HII regions by Peeters et al. (2001), we present a first analysis of the hydrogen recombination and atomic fine-structure lines originated in the ionized gas. The sample consists of 34 HII regions located at galactocentric distances between Rgal = 0 and 15 kpc. The SWS HI recombination lines between 2 and 8 mum are used to estimate the extinction law at these wavelengths for 14 HII regions. An extinction in the K band between 0 and $\sim$ 3 mag. has been derived. The fine-structure lines of N, O, Ne, S and Ar are detected in most of the sources. Most of these elements are observed in two different ionization stages probing a range in ionization potential up to 41 eV. The ISO data, by itself or combined with radio data taken from the literature, is used to derive the elemental abundances relative to hydrogen. The present data thus allow us to describe for each source its elemental abundance, its state of ionization and to constrain the properties of the ionizing star(s). △ Less

Submitted 30 October, 2001; originally announced October 2001.

Comments: Accepted in Astronomy and Astrophysics, 22 pages, 20 figures, 9 tables

Journal ref: Astron.Astrophys.381:606-627,2002

The Warm Ionized Medium in the Milky Way and Other Galaxies

Authors: John S. Mathis

Abstract: Observations of the "Warm Ionized Medium" (or, equivalently, the "Diffuse Ionized Gas") of the local ISM, the Perseus arm in the Milky Way, and also in several other galaxies show strong [NII]6563 (~H-alpha in some cases) and [SII]6717/[NII]6583 = 0.6 - 0.7 in all locations and objects. Other line ratios (e.g., [O III]5007/H-beta) vary considerably. Simple photoionization models reproduce the ob… ▽ More Observations of the "Warm Ionized Medium" (or, equivalently, the "Diffuse Ionized Gas") of the local ISM, the Perseus arm in the Milky Way, and also in several other galaxies show strong [NII]6563 (~H-alpha in some cases) and [SII]6717/[NII]6583 = 0.6 - 0.7 in all locations and objects. Other line ratios (e.g., [O III]5007/H-beta) vary considerably. Simple photoionization models reproduce the observed spectra, providing extra heating beyond that supplied by photoionization is assumed (Reynolds, Haffner, & Tufte 1999). With observed gas-phase abundances (not solar), the line ratios in the local arm at b = 0 deg are fitted with no extra heating and (S/H) = 13 ppm (solar is 20 ppm). Local gas observed at b = -35 deg requires extra heating of about gamma = 0.75, where gamma is the extra heating in units of 10^{-25} erg H^{-1} s^{-1}. In the Perseus arm, there are similar results, with a domposition consistent with the Galactic abundance gradient. The requirements for NGC 891 are similar to the Perseus arm: little or no extra heating at |z| = 1 kpc and gamma 3 at 2 kpc. In NGC 891 there is also an increase of 5007/H-alpha with |z| that can only come about if most of the ionizing radiation is supplied by stars with T~50000 K. Either their radiation must propagate from the plane to high |z| through very little intervening matter, or else the stars are located at high |z|. The total power requirement of the extra heating is <15% of the photoionization power. [O~II]3727/H-beta can serve as a useful diagnostic of extra heating, but [S~III] 9065,9531/H-alpha is not useful in this regard. △ Less

Submitted 10 October, 2000; originally announced October 2000.

Comments: 32 pages, including 2 figures. To appear in November 20 ApJ

Dependence of Gas Phase Abundances in the ISM on Column Density

Authors: B. P. Wakker, J. S. Mathis

Abstract: Sightlines through high- and intermediate-velocity clouds allow measurements of ionic gas phase abundances, A, at very low values of HI column density, N(HI). Present observations cover over 4 orders of magnitude in N(HI). Remarkably, for several ions we find that the A vs N(HI) relation is the same at high and low column density and that the abundances have a relatively low dispersion (factors… ▽ More Sightlines through high- and intermediate-velocity clouds allow measurements of ionic gas phase abundances, A, at very low values of HI column density, N(HI). Present observations cover over 4 orders of magnitude in N(HI). Remarkably, for several ions we find that the A vs N(HI) relation is the same at high and low column density and that the abundances have a relatively low dispersion (factors of 2-3) at any particular N(HI). Halo gas tends to have slightly higher values of A than disk gas at the same N(HI), suggesting that part of the dispersion may be attributed to the environment. We note that the dispersion is largest for NaI; using NaI as a predictor of N(HI) can lead to large errors. Important implications of the low dispersions regarding the physical nature of the ISM are: (a) because of clumping, over sufficiently long pathlengths N(HI) is a reasonable measure of the_local_ density of_most_ of the H atoms along the sight line; (b) the destruction of grains does not mainly take place in catastrophic events such as strong shocks, but is a continuous function of the mean density; (c) the cycling of the ions becoming attached to grains and being detached must be rapid, and the two rates must be roughly equal under a wide variety of conditions; (d) in gas that has a low average density the attachment should occur within denser concentrations. △ Less

Submitted 2 October, 2000; originally announced October 2000.

Calculating Cross Sections of Composite Interstellar Grains

Authors: Nikolai V. Voshchinnikov, John S. Mathis

Abstract: Interstellar grains may be composite collections of particles of distinct materials, including voids, agglomerated together. We determine the various optical cross sections of such composite grains, given the optical properties of each constituent, using an approximate model of the composite grain. We assume it consists of many concentric spherical layers of the various materials, each with a sp… ▽ More Interstellar grains may be composite collections of particles of distinct materials, including voids, agglomerated together. We determine the various optical cross sections of such composite grains, given the optical properties of each constituent, using an approximate model of the composite grain. We assume it consists of many concentric spherical layers of the various materials, each with a specified volume fraction. In such a case the usual Mie theory can be generalized and the extinction, scattering, and other cross sections determined exactly. We find that the ordering of the materials in the layering makes some difference to the derived cross sections, but averaging over the various permutations of the order of the materials provides rapid convergence as the number of shells (each of which is filled by all of the materials proportionately to their volume fractions) is increased. Three shells, each with one layer of a particular constituent material, give a very satisfactory estimate of the average cross section produced by larger numbers of shells. We give the formulae for the Rayleigh limit (small size parameter) for multi-layered spheres and use it to propose an ``Effective Medium Theory'' (EMT), in which an average optical constant is taken to represent the ensemble of materials. Multi-layered models are used to compare the accuracies of several EMTs already in the literature. △ Less

Submitted 21 August, 1999; originally announced August 1999.

Comments: 29 pages, 6 figures, accepted for publication in the Astrophysical Journal (part 1, scheduled in Vol. 526, #1, Nov. 20)