abstract: The interplay between string theory and cosmology is very promising. Since string theory will yield a quantum theory of space-time and unify all forces of nature, it has the potential of addressing many of the conceptual problems of the current models of early Universe cosmology. In turn, cosmology is the most obvious testing ground in the effort to construct non-perturbative string theory, and can provide the crucial connection between theory and experiment/observation [^1].
author: Robert H. Brandenberger,
Department of Physics, Brown University, Providence, RI 02912, USA.
date: 2024-10-03
title: The Promise of String Cosmology [^2]

**Early Universe cosmology** has made spectacular progress over the past two decades, driven both by a wealth of new observational data and by theoretical breakthroughs. In particular, a scenario for the very early Universe, the **inflationary Universe** scenario , has emerged which has led to solutions of some of the deep mysteries of standard big bang cosmology, and which has made some quite specific predictions for the geometry of the Universe and for the spectrum of cosmic microwave background (CMB) anisotropies which have been confirmed by the most recent observations .

In spite of its success, the inflationary Universe scenario suffers from some deep conceptual problems which have prevented it from achieving the status of an honest physical theory. For example, scalar-field-driven inflationary cosmology is incomplete in the same sense as the standard big-bang cosmology was: in the context of inflationary cosmology an initial singularity is inevitable, thus making it impossible to formulate a consistent initial value problem. More seriously, the successful predictions of inflationary cosmology for the spectrum of the CMB depend - in the current models of inflation in which the exponential expansion of space is driven by the potential energy of a scalar matter field - on extrapolating the physical theory into a regime where its foundations (namely classical general relativity and weakly coupled quantum field theory) are known to break down [^3]. The Achilles heel of inflationary cosmology is the cosmological constant problem. How do we know that the unknown mechanism which cancels the vacuum energy today does not also cancel the transient vacuum energy which is required to drive inflation? Quite apart from these conceptual problems, a convincing realization of inflation in the context of our current models of particle physics based on quantum field theory is lacking.

Recent observations are also beginning to challenge some of the premises of our current cosmological models. These observations point to the existence of a dark energy component in the Universe which makes up the bulk of the energy density but which does not cluster . Furthermore, recent measurements and numerical simulations show a mismatch between the predicted and observed galaxy halo structures which may imply that the dark matter is not cold (see e.g. and references therein).

**String (or M-) theory** is the best candidate for a unified theory of all forces and as such for a quantum theory of space-time and matter. Thus, it will resolve the conflict between quantum field theory and general relativity, and therefore it has the potential of addressing many of the open issues in early Universe cosmology. If successful in this challenge, string theory will have made an important contact with experiment/observation, and will no longer be able to be criticized for lacking a connection to data and for thus being pure mathematics rather than a physical science. The successful determination of the entropy of an extremal black hole in the context of string theory is a first major success in establishing a link between string theory and cosmology.

Until recently there has not been much work on the interface between string theory and cosmology. However, it has been clear since the late 1980's that string theory has the potential to solve some of the problems of standard cosmology. For example, it was argued , in the context of perturbative string theory, and assuming that all spatial dimensions are toroidally compactified, that as a consequence of t-duality the physics at very small radii is equivalent to that at very large radii, hence eliminating the cosmological singularities. Moreover, string winding modes may yield a mechanism which allows at most three spatial dimensions to become large, thus providing a dynamical resolution of the potential embarrassment that superstring theory is consistent only in nine (or ten in the context of M-theory) spatial dimensions. This example demonstrates that with string theory, questions about the Universe can now be addressed scientifically which before the advent of string theory were exclusively in the realm of philosophy.

In the past five years string theory has undergone a "second revolution". With the discovery of **D-branes** it has become clear that there are many more fundamental degrees of freedom than just strings. Dualities between the previously known five consistent string theories were discovered which indicate that there is a common underlying theory, **M-theory**. The different string theories (and also 11-d supergravity) correspond to certain corners of moduli space of M-theory. A further new important tool is the **AdS-CFT correspondence** which relates classical gravity on an anti-de-Sitter background to a conformal field theory living on the boundary. Recently, much attention has been focused on **brane world scenarios**, scenarios based on a combination of the above ideas in which the matter fields of our Universe are confined to a three-brane and only gravity lives in the higher-dimensional bulk (see e.g. ). In general, the spatial dimensions transverse to the branes must be very small in order to avoid too large deviations from Newton's gravitational laws. An interesting suggestion first put forward without a direct connection to string theory [^4] is that two of the extra dimensions are large (mm scale). Since the four dimensional effective gravitational constant is related to the fundamental higher dimensional one via the size of the extra dimensions, the fundamental gravitational constant can be reduced to the TeV scale, provided the large extra dimensions have a size of the order of 1mm. This provides the potential for solving the mass hierarchy problem of elementary particle physics. By constructing models in which gravity remains confined to the brane, it is possible to make the extra dimensions large (in fact even infinite ).

On the other hand, there is still no non-perturbative formulation of M-theory; not even its fundamental physical degrees of freedom are known. In at least one approach to M-theory , the underlying theory is postulated to be a quantum mechanical matrix model. If this is true, then some of the non-commutativity of the matrix model will be reflected in non-commutativity of the space-time coordinates at a microscopic level. Returning to the more traditional description of M-theory, another serious problem results from the fact that the moduli space of vacuum states corresponding to any of the string theories is huge, and many of the predictions for particle physics as well as for cosmology depend on which vacuum state is chosen. A key challenge to string cosmology is to find a mechanism which distinguishes the vacuum state corresponding to our Universe. Since the moduli correspond to massless fields, new problems for cosmology appear: why don't we see the massless particles associated with the moduli fields, and if this question is answered in the usual way by invoking some non-perturbative mass generation for the moduli, then why don't the moduli over-close the Universe ?

However, even without first having to solve these difficult problems, there are interesting avenues for string cosmology. A particularly interesting issue is to determine whether there are string-specific ways to obtain a period of inflationary expansion in the early Universe. If the answer to this question is affirmative, then it is of utmost importance to try to find observational signatures which could distinguish these string-specific inflationary models from scalar-field-driven inflation. The answer to this question will inevitably require looking beyond the power spectra of density fluctuations and of CMB anisotropies, since any set of power spectra can be obtained by suitable choice of a toy model scalar field potential. It may, however, also turn out that string cosmology does not yield anything that looks like the usual inflationary scenario, but that the resulting theory nevertheless provides an alternative to inflation in solving the problems of standard big bang cosmology such as the horizon and flatness problems. This possibility is realized in the pre-big-bang approach to string cosmology based on dilaton gravity as the low-energy effective theory consistent with the symmetries of string theory. In this scenario the background in the Einstein frame is in fact contracting super-exponentially in the pre-big-bang phase, leaving an effective Hubble constant whose absolute value is increasing in time. In pre-big-bang cosmology the spectrum of scalar metric fluctuations stemming from the dilaton gravity sector is not scale-invariant as it would be in inflationary cosmology . A more radical alternative to inflation which can solve the homogeneity and flatness problems is the varying-speed-of-light scenario in which the speed of light in the early stages of the evolution of the Universe is postulated to be many orders of magnitude larger than at the present time .

Another direct observational issue which string cosmology should address concerns the origin of dark energy. Since string theory leads to a large number of scalar moduli fields, it is of great interest to explore the potential of these fields for explaining the origin of dark energy. Possibly related to this issue is the question of the smallness of the observed cosmological constant. Finally, even without knowing the exact vacuum state of string theory it should be possible to study whether string theory will indeed wash out the singularities of standard and inflationary cosmology.

Interest in string cosmology is growing. From July 24 - August 4, 2000, about 45 physicist working on string theory, non-perturbative field theory and cosmology gathered at the University of British Columbia (UBC) in Vancouver for a workshop devoted to string cosmology. This workshop was co-sponsored by the Pacific Institute for Mathematical Sciences (PIMS), the Canadian Institute for Advanced Research (CIAR) and by the APCTP. The invited participants included B. Greene (Columbia Univ.), N. Kaloper (Stanford), L. Kofman (CITA, Univ. of Toronto), D. Lowe (Brown Univ.), B. Ovrut (Univ. of Pennsylvania), S. Ramgoolam (Brown Univ.), S. Sin (Hanyang Univ.), D. Son (Columbia Univ.), P. Steinhardt (Princeton Univ.), H. Verlinde (Princeton Univ.), G. Veneziano (CERN) and A. Zhitnitsky (UBC) [^5]. Many different approaches to string cosmology were represented at this workshop. In the following discussion of various approaches to string cosmology I will limit myself to those which were discussed extensively at the workshop.

Based on the conclusions of one of the most conservative approach to string cosmology, namely **Pre-Big-Bang Cosmology** , it is clear that we should expect important differences in the evolution of the early Universe in string cosmology compared to scalar-field-driven inflationary cosmology. Pre-big-bang cosmology is a low energy description of physics which takes into account all of the low energy modes in string theory, i.e. not only the graviton but also the dilaton and the anti-symmetric tensor field (the latter, however, does not play an important role for the background dynamics at the level of the homogeneous equations of motion [^6]). The symmetries and equations of motion for dilaton gravity give rise to a scenario of cosmology in which the Universe starts out in a state near the perturbative string vacuum, goes through a dilaton-driven phase of contraction during which the Hubble radius $H^{-1}(t)$ ($H$ being the expansion rate of the Universe and $t$ denoting time) shrinks faster than the physical length of fixed comoving scales, thus simultaneously solving the horizon problem of standard cosmology and giving rise to a mechanism for structure formation similar as in the usual inflationary cosmologies: quantum vacuum fluctuations existing during the contraction phase on sub-Hubble lengths get frozen in when the Hubble radius becomes smaller than the wavelength of the fluctuation, they are amplified super-adiabatically while the wavelength is larger than the Hubble radius, and then re-enter the Hubble radius in the late Universe as classical density fluctuations. As the Universe contracts, the curvature increases. Eventually, in a way which unfortunately cannot be described by the equations of dilaton gravity, the Universe makes a transition to an expanding Friedmann cosmology related to the contracting phase via a duality transformation. Note that in pre-big-bang cosmology, both the background evolution and the spectrum of induced fluctuations is different from that of scalar-field-driven inflationary models. The density fluctuations from the dilaton sector are not scale invariant. In fact, they are completely unimportant on large length scales. The other fields present in pre-big-bang cosmology (such as the axion) will, however, generate fluctuations which can be scale-invariant . Since these fluctuations start out as isocurvature fluctuations, they lead to different predictions for observables such as the spectrum of CMB anisotropies than the scalar metric fluctuations of inflationary cosmology. The lesson we learn is that it is quite likely that the predictions of string cosmology will differ enough from those of standard inflationary cosmology to enable observations to probe the physics of the Planck scale.

Pre-big-bang cosmology suffers from the graceful exit problem: in the context of the equations of motion of dilaton gravity, the contracting branch of the evolution terminates in a singularity, and the expanding branch emerges from an equivalent singularity. Since the singularity occurs in the high curvature regime, it is clear that the effective dilaton gravity action no longer describes the stringy physics in a correct way close to the singularity. A challenge for string cosmology is to understand if true stringy physics can solve the graceful exit problem and lead to a smooth transition from the contracting dilaton-dominated phase to the expanding phase with fixed value of the dilaton. Although there have been several attempts to go beyond dilaton gravity in order to solve the graceful exit problem (see e.g. for a recent overview), there is not yet a convincing approach. An important open problem for string cosmology is to develop a truly stringy, non-perturbative version of the pre-big-bang scenario, or to find out how the scenario must be modified in order to become a truly stringy cosmological model.

Another promising approach to string cosmology is based on Horava-Witten theory and is sometimes called **heterotic M-theory cosmology** . The theory is based on one particular ray in the moduli space of M-theory, namely the compactification of 11-dimensional supergravity on an orbifold $S^1 / Z_2$. This leads to two distinguished co-dimension one planes, the orbifold fixed planes. As usual, it is assumed that, in addition, six spatial dimensions are compactified on a special type of Calabi-Yau manifold, and that there is a set of $E_8$ gauge supermultiplets on each orbifold plane. The reason for adopting this particular orbifold compactification is that it allows for chiral fermions . The resulting theory is dual to the strongly coupled $E_8 \times E_8$ heterotic string theory. Matching at tree level to the phenomenological gravitational and grand unified gauge couplings, one finds that the orbifold must be larger than the Calabi-Yau radius. This suggests that there is a regime in which the dynamics of the Universe is effectively five-dimensional. The resulting five-dimensional effective theory is a gauged $N = 1$ supergravity coupled to four-dimensional boundary gauge and matter supermultiplets (see e.g. for a recent survey). Initially, only the effective four-dimensional equations were studied. However, it has recently been realized that the dynamics in the orbifold direction is quite nontrivial and may lead to new effects of great interest to cosmology. In particular, the equations of motion in five dimensions admit as a static solution a pair of three-branes parallel to the orbifold fixed planes. It is tempting to try to obtain inflation from the nontrivial dynamics involving the orbifold direction. Note, however, that as in pre-big-bang cosmology, many truly stringy effects are frozen out in the effective action-based analysis of heterotic M-theory cosmology.

A rather different approach to string cosmology is based on proceeding in analogy to the usual starting point of big bang cosmology, which is the assumption that the Universe starts out with Planck size and emerges from a hot initial state in which all degrees of freedom were highly excited. It seems quite natural to generalize this starting point to string cosmology and to assume that the Universe starts out with all spatial dimensions of string scale and all fundamental degrees of freedom excited according to their thermal distribution. In particular, this means that the winding modes of strings and branes are highly excited. This approach is called the **brane gas approach** to string cosmology. In this approach, all spatial dimensions start out compact and the question for cosmology is not why certain dimensions are compactified, but why a fixed number of dimensions can dynamically decompactify (in the sense of their radii becoming larger than the present Hubble radius of the Universe). There have been some interesting recent developments along the lines of this approach to string cosmology. In particular, assuming that all spatial dimensions are compactified on a torus, it has been shown that whereas at early times (small sizes) the higher dimensional branes dominate the evolution, it is the fundamental strings which have the most important effect at late times . Thus, the mechanism of , which only allows three spatial dimensions to become large, survives in the modern view of string theory. In addition, the duality under inversion of the radius of the torus is maintained, thus assuring that no physical cosmological singularities will arise. It is important to generalize these considerations to more realistic background spaces.

The approach to string cosmology which at the present time is attracting most attention is the class of **brane-world scenarios** (see e.g. for work along these lines presented at the APCTP/CIAR/PIMS string cosmology workshop). The starting point of this approach is the assumption that our matter fields are confined to a three-brane (four space-time dimensions) embedded in a higher-dimensional (usually five-dimensional) bulk space-time. An important realization is that it is possible to confine gravity on one of the branes by having the four-dimensional metric depend on the fifth coordinate (the orbifold coordinate in the language of heterotic M-theory cosmology). In many cases, the bulk metric is AdS, and this gives rise to a tantalizing connection with the AdS-CFT correspondence. However, most work along these lines is not directly justified from string theory. There appear to be an unlimited number of versions of the basic scenario. As was stressed at the APCTP/CIAR/PIMS string cosmology workshop, an urgent task for string cosmology is to derive the basic rules which a brane-world scenario must obey if it is to emerge from string theory. Is the brane an orbifold fixed plane (as it is in heterotic M-theory cosmology) or is it a D-brane? How many branes are there in the scenario? Are there matter fields in the bulk? At the present time, most work on brane-world cosmology is based on studying the background gravitational equations of motion of the brane in the presence of the bulk, making use of the general relativistic Israel matching conditions. It has been realized that the equations for the dynamics of our brane in the brane-world scenario may be quite different from the usual four-dimensional general relativistic equations. This provides stringent constraints on brane world cosmologies. However, it also provides the chance to obtain new effects, or to provide new sources of inflation. In order to make contact with the wealth of recent cosmological observations on the distribution of light and CMB light, it is crucial to extend the theory of cosmological perturbations, the basis for making such a contact, from the usual four-dimensional setting to the brane-world scenarios. Initial work along these lines was presented at the workshop (see and references therein). A limitation which the brane-world scenarios share with some of the other approaches to string cosmology is that most of the analysis is based on a truncation of the dynamics to low energy modes, thus eliminating many truly stringy effects. It is, however, possible to obtain interesting new effects for cosmology. For example, if our Universe is a three-brane moving in a bulk black hole background, it is possible to find a realization of the varying light scenario for solving the problems of standard cosmology.

Ultimately, however, a complete cosmological model should be based on the non-perturbative formulation of M-theory. Starting from the matrix theory approach to M-theory (see e.g. ), it seems likely that the resulting space-time theory will be non-commutative at the string scale (see e.g. for work along these lines presented at the workshop). One limitation of the matrix theory approach to M-theory is that the theory is at the present time only defined in light-cone coordinates, which makes the connection to a space-time formulation very challenging.

A generic problem which appears in all approaches to string cosmology is the moduli problem. What freezes the value of the dilaton and what sets the scale of the extra dimensions? Usually it is assumed that some non-perturbative physics is responsible for the answers. If an effective potential provides the required mechanism (see e.g. for some recent work presented in Vancouver), then why does the value of the potential at its minimum not lead to another source for the cosmological constant? Thus, the moduli problem and the cosmological constant problem appear to be related. There are some new approaches to solving the cosmological constant problem based on the brane-world scenario, e.g. a model by Verlinde in the context of a theory with two branes, the physical brane and the Planck brane, which assumes that the cosmological constant on the Planck brane is protected by supersymmetry, and that it flows to a small cosmological constant on the physical brane via renormalization group flow.

The field of string cosmology is still in its infancy. Although many interesting approaches to connecting string theory and cosmology exist, most of them suffer from our ignorance about the true non-perturbative formulation of string theory. The promises of the field, however, are great, as witnessed by the excitement expressed at the APCTP/CIAR/PIMS string cosmology workshop

**Acknowledgements**

I wish to thank the APCTP, CIAR and PIMS for their generous support of the *String Cosmology Workshop* held at UBC in Vancouver as part of the *Frontiers of Mathematical Physics* workshop series jointly organized by PIMS and APCTP. I also thank R. MacKenzie for useful comments on the draft. The author's research was supported in part by the U.S. Department of Energy under Contract DE-FG02-91ER40688, TASK A.

[^1]: For a more technical overview of the existing approaches to connecting string theory and cosmology, the reader is referred in .

[^2]: This article was commissioned in conjunction of the Workshop on String Cosmology held from July 24 - Aug. 4 2000 at PIMS, Univ. of British Columbia and is published in APCTP Bulletin \# 6, (2001). Brown preprint BROWN-HET-1243.

[^3]: This is the so-called trans-Planckian problem for inflationary cosmology .

[^4]: See for a realization of this idea in the context of string theory.

[^5]: See the web site http://kepler.physics.ubc.ca/p̃fs99 for details about the workshop.

[^6]: As was mentioned at the workshop, a non-vanishing tensor field decreases the set of inhomogeneous initial conditions which can lead to the pre-big-bang evolution.


author: Serena Viti
bibliography: references.bib
date: Received; accepted
title: Molecular transitions as probes of the physical conditions of extragalactic environments

We present a method to interpret molecular observations and molecular line ratios in nearby extragalactic regions. Ab initio grids of time dependent chemical models, varying in gas density, temperature, cosmic ray ionization rate, and radiation field, are used as inputs into RADEX calculations. Tables of abundances, column densities, theoretical line intensities, and line ratios for some of the most used dense gas tracers are provided. The degree of correlation as well as degeneracy inherent in molecular ratios is discussed. Comparisons of the theoretical intensities with example observations are also provided. We find that, within the parameters space explored, chemical abundances can be constrained by a well-defined set of gas density, gas temperature, and cosmic ray ionization rates for the species we investigate here. However, line intensities, and more importantly line ratios, from different chemical models can be very similar, thereby leading to a clear degeneracy. We also find that the gas subjected to a galactic cosmic ray ionization rate will not necessarily have reached steady state in 1 million years. The species most affected by time dependency effects are HCN and CS, which are both high density tracers. We use our ab initio method to fit an example set of data from two galaxies, i.e. M82 and NGC 253. We find that (i) molecular line ratios can be easily matched even with erroneous individual line intensities, (ii) no set of species can be matched by a one-component interstellar medium (ISM), and (iii) a species may be a good tracer of an energetic process but only under specific density and temperature conditions. We provide tables of chemical abundances and line intensities ratios for some of the most commonly observed extragalactic tracers of dense gas for a grid of models. We show that by taking the chemistry behind each species and the individual line intensities into consideration, many degeneracies that arise by just using molecular line ratios can be avoided. Finally we show that using a species or a ratio as a tracer of an individual energetic process, such as cosmic rays and UV, ought to be done with caution.

# Introduction

Following the first detection of a molecular species in external galaxies in the 1970s, molecular line emission studies have been routinely performed to determine the physical characteristics of nearby galaxies . Galaxies span a variety of shapes, sizes, and physical conditions. The physical parameters, such as gas densities, ultraviolet (UV) fields, cosmic ray ionization rates, and dust properties, which determine the appearance of our own Galaxy may therefore be very different in other galaxies. Nevertheless external galaxies, in general, contain the variety of regions and sources similar to those that we can identify in the Milky Way; hence molecules could be ideal probes of the dense gas that is the reservoir for star formation and could help in describing the process of star formation itself as well as the influence of newly formed stars on their environments .

Because of their range of critical densities across the different molecular species and across the transitions of the same molecule, molecular observations are an ideal tool to trace a wide range of gas densities in the interstellar medium and the dependencies of chemical reactions on the energy available to the system, i.e the kinetic temperature of the gas. However, what makes observations of molecules so adept is also what complicates their interpretation. In particular, even with the highest spatial resolution available with the Atacama Large Millimeter/Submillimeter Array (ALMA), the beam usually encompasses emissions from many types of regions and the molecular emission is an ensemble of multi-phase gas in which the spatial and temporal effects are diluted in the beam. Concomitant detection of multiple molecules in a spatially unresolved galaxy therefore does not necessarily mean that they are emitted from the same gas . While this is not surprising, many multi-species studies of nearby galaxies still need to assume coexistence of molecular emission when interpreting their observations . This is partly because we are forced to use multiple transitions and molecules to overcome the well-known degeneracy between density and temperature for linear molecules, which implies that single transitions of these molecules cannot be easily used to determine quantitatively the gas density or temperature. There have been several theoretical studies on the sensitivity of chemistry to the range of physical parameters that determine a galaxy . Most of these sensitivity studies concentrated on the modelling of photon-dominated regions (PDR) or X-ray dominated regions (XDR), rather than the modelling of dense molecular gas. Observationally, many studies have shown that molecular line ratios, such as HCO$^{+}$/HCN or HCN/CO, differ across different types of galaxies especially between Active Galactic Nuclei (AGN)-dominated galaxies and starburst-dominated galaxies , although some studies have found that enhanced HCN/HCO$^+$ are not unique to AGN environments but can also be found in systems dominated by star formation . However the derived abundance ratios for an individual galaxy also highly differ across studies depending on the transitions observed, available resolution, and method used in deriving the column densities. An understanding of the chemistry behind each molecule and its dependencies on the density and temperature of the gas is therefore essential and will lead to a more coordinated approach to understanding the nature of galaxies using molecular line emissions.

In this paper we use, as a test bed, a relatively small grid of chemical models to present a method for the interpretation of molecular observations and molecular line ratios in nearby extragalactic regions, in particular, starburst and AGN dominated galaxies as these are the most targeted types of galaxies for dense gas studies; similar studies can be repeated for Ultra Luminous Infrared Galaxies (ULIRG), dwarf galaxies, etc., where, however, much of the gas may be in ionic or atomic form . We analyse the theoretical abundances and column densities of some of the most common tracers for a parameter space in gas density, temperature, cosmic ray ionization rate, and radiation field. The theoretical abundances are then used as inputs to the radiative transfer code RADEX , which in turn provides the theoretical intensities as well as common line ratios for the parameter space investigated by the chemical models. We investigate the degree of correlation and degeneracy inherent in each ratio and show how one can use theoretically derived abundances and line intensities for comparison with observations by using published observational data as examples.

The aim of this study is twofold. First, we seek to aid the observer in determining the most likely scenario that fits the data when multiple solutions are present. Second, we aim to provide online tables that one can use directly to compare with their observed intensities and intensity ratios for the most common used transitions; the tables are simply meant to be an example subset but include some of the most observed molecular species and transitions. We present the chemical models and discuss their degeneracies in Section 2. In Section 3 we first calculate theoretical line intensities and line ratios for four of the most observed gas tracers: CO, HCN, HCO$^+$, and CS. While the focus of this paper is purely on the characterization of the dense gas, CO (whose critical density is low and hence ubiquitously traces all gas above $\sim$ 100 cm$^{-3}$) is included because it is still routinely used to interpret molecular ratios or determine the star formation rate density , especially if high-J CO observations are available . Secondly, for completeness, we also extend our abundance and line intensities analyses to further species, not as commonly observed in external galaxies as those above, but often used as tracers of particular energetic events (e.g. shocks). In Section 4 we make use of our calculations to partly reinterpret observations in two very well studied galaxies to give examples of how this methodology can be used.

All the tables referred to in the paper are online and can be found at https://uclchem.github.io/viti2017.pdf.

# Chemical models: Theoretical molecular abundances

We used the `UCL_CHEM`[^1] gas-grain time dependent chemical code to run a grid of chemical models covering a range of gas densities (10$^4$–10$^6$ cm$^{-3}$), temperatures (50–200 K), radiation fields (1–500 Draine), cosmic ray ionization rates (1–500 standard galactic cosmic ray ionization field), and visual extinctions (1–50 mags); we used the cosmic ray ionization flux to simulate an enhancement in X-ray flux. As previously noted , this approximation has its limitations in that the X-ray flux heats the gas more efficiently than cosmic rays. However, the chemistry arising from these two fluxes should be similar.

While there are several more parameters that can be explored when performing chemical modelling, for the purpose of this study we chose those that affect most the line intensities of dense gas tracers in starburst and AGN-dominated galaxies, as predicted by Large Velocity Gradient (LVG) models. For example, gas density, temperature, and visual extinction (directly correlated to column density) are all input parameters to RADEX and other LVG models. Cosmic ray ionization and radiation fluxes are known to differ from the Galactic values in starburst and AGN-dominated galaxies, and hence are also considered. Moreover they directly affect the temperature of the gas and dust. The choice of the lower and upper limits for the gas densities and temperatures reflect the best fits found for the dense molecular gas in the most commonly studied starburst and AGN galaxies; they are by no means meant to represent the only possible values. The range of cosmic ray ionization rates and radiation fields are subsets of the wider range explored in previous theoretical studies. As already mentioned , our grid of models is not meant to be exhaustive but simply represent a plausible range of physical conditions. We also ran some models simulating the presence of shocks. For these models we adopted a representative shock velocity of 40 km s$^{-1}$, which depending on the pre-shock densities, leads to a maximum temperature for the shock ranging from 2200 K to 4000 K. The shocks in UCL_CHEM are treated in a parametric form as in . Below we briefly describe UCL_CHEM. Further details of this version of UCL_CHEM, coupled with the shock module, can be found in , and .

UCL_CHEM is an ab initio model and computes the evolution as a function of time of chemical abundances of the gas and on the ices starting from a diffuse and atomic gas (at a density of 10 cm$^{-3}$). The parameter regulating the visual extinction has two effects on the models. First, this parameter affects the photochemistry since the radiation field is very efficiently attenuated for A$_V$ $>$ 4-5 mags. However, and more importantly, this paramater is effectively a measure of the size of the emission region, since in order to derive column densities from the fractional abundances calculated by the models we multiply the latter by the hydrogen column density as measured at 1 magnitude times the visual extinction. This is the so-called "on the spot" approximation . As the chemical model is single point, we estimate the theoretical column densities, a posteriori.

The code is run in two phases in time (Phase I and Phase II). The initial density in Phase I is $\sim$ 10 cm$^{-3}$. The gas is allowed to collapse and to reach a high density gas by means of a free-fall collapse. The collapse is not meant to represent the formation of protostars, but it is simply a way to compute the chemistry of high density gas in a self-consistent way starting from a diffuse atomic gas, i.e. without assuming the chemical composition at its final density. Of course, it is unknown how the gas reached its high density and we could have assumed any other number of density functions with time. For simplicity, we adopt free-fall collapse. The temperature during this phase is kept constant at 10 K, and the cosmic ray ionization rates and radiation fields are at their standard Galactic values of $\zeta_o$ = 5$\times$`<!-- -->`{=html}10$^{-17}$ s$^{-1}$ and 1 Draine, or 2.74$\times$`<!-- -->`{=html}10$^{-3}$ erg/s/cm$^{2}$, respectively. In this way the chemistry of the high density gas is computed in a self-consistent way starting from a diffuse atomic gas, i.e without assuming the chemical composition at its final density. During Phase I atoms and molecules are allowed to freeze onto the dust grains and react with each other, forming icy mantles. In Phase II, for the non-shock models, UCL_CHEM computes the chemical evolution of the gas and the icy mantles after either an assumed burst of star formation or AGN activity has occurred. The temperature of the gas increases from 10 K to a value set by the user and sublimation from the icy mantles occurs. The chemical evolution of the gas is then followed for 10$^7$ years. In the shock models, Phase II considers a plane-parallel C-Shock propagating with a velocity, $V_s$, through the ambient medium. Details of the code and its use for the determination of abundances in external galaxies can be found in . A detailed description of the shock module in UCL_CHEM can be found in and in . In both phases of the UCL_CHEM model, the gas phase chemical network is based on the UMIST13 database with updates from the KIDA database . The surface reactions included in this model are assumed to be mainly hydrogenation reactions, allowing chemical saturation when possible.

Tables 1 and 2 list the grid of models with the resultant abundances (Table 1), column densities, and ratios (Table 2) for four of the most observed molecules in galaxies: CO, HCO$^+$, HCN, and CS. Both tables show the results at the final time step (10$^7$ years). However for the shock models, the shock occurs at the beginning of Phase II and the models are stopped once the dissipation length has been reached; hence the Tables show the chemistry at the last time step, representing the post-shocked gas (at $\sim$ 10$^5$ years). For the shock models the $A_V$ is the initial visual extinction before the passage of the shock. In Table 3 we tabulate the column densities for selected models at an earlier time of 1 Myr (see Section 2.2).

Tables 4, 5, and 6 list the column densities for further species, often used as tracers of shocks (e.g. SiO, CH$_3$OH, and SO) or surface grain chemistry (e.g. HNCO). HNC has been added because several studies use the HNC/HNC ratio as a tracer of dense PDR and XDRs environments . We note that HNCO and CH$_3$OH do not survive until chemical equilibrium is reached and they are only listed for 10$^6$ years for all non-shock models. Table 6 lists the column densities for CH$_3$OH and HNCO for Models 65-67 at two different times: during the peak temperature of the shock and once the shock has passed. For the latter time the column densities were calculated assuming an average visual extinction of 10 mags.

## Chemical 'degeneracy'

From Tables 1–6 we find a chemical degeneracy worth mentioning in relation to this study. Models differing by one or more parameters can still lead to column densities similar enough to be within the observational error in the column density estimation. For example, the HCN column density remains the same when cosmic ray ionization rate is enhanced from 1 to 10 $\times$ the Galactic value (cf. Models 1 and 9); however this does not mean HCN is not sensitive to changes in the cosmic ray ionization rates as higher cosmic rate affect the HCN chemistry. In most cases this degeneracy can be dealt with if one observes multiple molecular species at the same spatial resolution. For example, Models 3 and 6 lead to almost identical CO, and similar HCO$^+$ and HCN column densities but they differ by $\sim$ one order of magnitude in CS. It is key to know a priori which species one needs to observe to differentiate among models by referring to Tables such as those presented here. However a more interesting degeneracy is caused when one obtains the same line intensities for different combinations of \[n$_H$,T$_{gas}$, $N(X)$\]. We discuss this in Section 3.

## The importance of time dependency

One chemical effect that is difficult to discern observationally is the time dependency of the abundance of each species. We do not discuss the shock models in this section, as chemical equilibrium is never reached during the passage of the shock. As shown by several studies in the past, not all chemistry is equally affected by time. Hence, while for some of the models we ran, equilibrium is reached well before 1 million years, for others it is not reached until much later; this is important because it is not obvious whether the large scale gas we observe, especially in starburst galaxies, will have remained subjected to the same energetic and dynamical conditions for longer than 1 million years because, for example, the flux of cosmic ray will vary across episodes of starbursts. Time dependencies effects have been discussed before , in the context of specific galaxies as well . In Tables 3 and 5 we tabulate the column densities for the models where chemical equilibrium has not been reached yet at 1 million years, and hence where the column densities differ from Tables 2 and 4. In Tables 3 and 5 there are no models with $\zeta$ higher than the standard value; it is clear therefore that a high cosmic ray ionization rate leads to steady state more quickly, as also found by previous work . On the other hand, an enhanced radiation field (at least up to 500 the interstellar galactic field) does not help the chemistry in reaching steady state at the visual extinctions considered here, since the UV does not penetrate. Interestingly, of the four most observed species (CO, HCN, HCO$^+$, and CS), the species most affected by time dependency is HCN followed by CS: both of these species are routinely used as high density gas tracers.

Some species are always short-lived; in the parameter space explored, these are CH$_3$OH and HNCO. In all the non-shock models these two species are always destroyed by the time chemical equilibrium is reached. In models not listed in Table 5, these species do not even survive for a 1 Myr; hence an observed high abundance of these two species is either a signature of a young gas or, more likely, a dynamically excited gas. In fact, we find that for some of the models in which HNCO and/or CH$_3$OH are very low before or at 10$^6$ years, their abundance is well enhanced at earlier times.

# Radiative transfer calculations: Intensities and intensities ratios

We estimated the theoretical integrated intensity for the first 10 rotational lines of CO, HCN, HCO$^+$, and CS, and for selected lines of SiO, HNC, SO, CH$_3$OH, and HNCO, where the line selection is observationally motivated. The intensities were calculated using the radiative transfer code RADEX from the grid of models in Tables 2-6. The collisional rates used for the RADEX calculations were taken from the LAMBDA database and were computed with H$_2$ as collisional partner, where available. We used a background temperature of 2.7 K. The approach we take in this paper differs from the usual use of RADEX or any LVG modelling in external galaxies in that we do not produce a grid in order to fit a set of observations. Rather, we used the input from our chemical modelling to fix the three-parameter set \[n$_H$, T$_{gas}$, $N(X)$\]. Hence the only two free parameters in RADEX are the linewidth and geometry. We take two representative linewidths of 50 and 100 kms$^{-1}$, but as the results do not affect our main conclusions we only present the 100 kms$^{-1}$ models. Although RADEX can be used with different geometries, for the purpose of these calculations we used spherical geometry. While with Tables 2-6 such calculations can be repeated for more geometries, used all three geometries available to characterize the dense gas in the AGN dominated nucleus of NGC 1068 but did not find significant differences in their fitting.

Also, estimating theoretical intensities for the shock models using RADEX falls well short of a proper radiative transfer calculation since the gas temperatures and densities vary quickly with the passage of the shock, and so do the abundances; hence one must be aware that these line intensities represent a gas that has been shocked (by a single or multiple shock episodes) but that has cooled back down to its 'equilibrium' temperature (in this case 100 K). Clearly a large variation of the parameter space can be performed when modelling shocked gas, which goes beyond the scope of this paper.

The results of our calculations are tabulated in Tables 7 and 8, for the most commonly observed species, namely CO, HCO$^+$, HCN and CS; in Table 7 we list chemical abundances as obtained at 10$^7$ years or at the last time step in the case of the shock models and in Table 8 we list chemical abundances as obtained at 10$^6$ years for the models in which chemical equilibrium was not reached at this time. In Tables 10–16 we tabulate selected line intensities for SiO, SO, HNC, and HNCO. Because of the large number of CH$_3$OH transitions observed as $k$ bands, we chose not to tabulate its line intensities (which can however be easily derived from the column densities tables.) For some models, at 10$^6$ years, the RADEX calculation of the HCN does not converge; what these models have in common, at 1 million years, is that the HCN column density is very high and the gas is either very dense (10$^6$ cm$^{-3}$) and/or very hot (200 K). The issue of non-convergence is most likely due to high-opacity effects, and in these cases Monte Carlo techniques combined with an accelerated lambda convergence method have been shown to perform better .

What it is often used in molecular studies in external galaxies are molecular line ratios. In Table 9 we tabulate the HCN/HCO$^+$ ratios for each transition; we deliberately approximate to its integer value as it is seldom that one can claim a better accuracy from observational values. It is clear even from a quick glance that the often used ratios, for example the HCN(4-3)/HCO$^+$ (4-3) or the HCN(1-0)/HCO$^+$(1-0), are the same for several models within the accuracy assumed. While this is unavoidable, by approaching the interpretation of the observations with an a priori method such as presented here, one may, at least, be able to constrain this degeneracy. As a practical example, let us consider the HCN/HCO$^+$ line ratios, which we plot in Figure 1 for four different transitions. These ratios are routinely used as tracers of AGNs because they are found to be significantly higher in AGN when compared to those found in starbursts . Model 1 has the same HCN/HCO$^+$(1-0) and similar HCN/HCO$^+$(2-1) ratio to Model 7 but all the other ratios are very different. These two models differ substantially only in gas density, which is a physical characteristic that the lower transition ratios could not have traced. Figure 1 shows how a relatively standard grid of chemical models can lead to large ranges for the HCN(1-0)/HCO$^+$ ratio of key transitions. This large range of ratios is consistent with the fact that, as pointed out in , among others, differences in the intensity ratios may arise due to different densities, temperatures, radiation fields, and evolutionary state of activity. However, as long as more than two transitions per molecules are observed, degeneracies can indeed be broken, at least within the parameter space explored here.

The approach we outline here is not proposed as a substitute for the more canonical use of LVG modelling but as a complementary analysis. In particular one could use tables, such as those presented in this study, to determine a priori whether similar line intensities arise from very different physical conditions and viceversa. This may be of particular importance in light of the results from who find that LVG models struggle to recover any parameter better than to within half a dex and they find no evidence of systemic offsets. While these Tables are not exhaustive we hope that they can be used to track the origin of multiple solutions.

# Characterizing the dense gas in nearby galaxies

In this section we provide examples of how one can use theoretically determined molecular intensities and ratios to interpret observations. We refer to two commonly observed nearby galaxies as our test beds and we compare an arbitrary selection of single-dish observations for one of the galaxies, and interferometric observations for the second galaxy. What is presented here is not a full model of any of these two galaxies. In order to do that a much more comprehensive set of observations, together with a larger grid of models, ought to be used.

## M82

M82 (at a distance of $\sim$ 3.50 Mpc) is one of the most studied nearby starburst galaxies. Among the many works, we chose, as an example, the observations presented in who performed a line survey of the north-east (NE) lobe of M82 with the single dish IRAM 30 m telescope. The spatial resolution obtained ranged from 158 to 333 kpc.

Of the molecular transitions they observed, we first compared CS (3-2), CS (5-4), and the HCO$^+$(3-2). We directly compared the integrated intensities (Table A.1 in by making the crude assumption that the emission region fills the beam. In fact, according to , and based on previous interferometric observations, the emission region should be within a factor of 1–1.5 of the beam size, hence our assumption is justified. However, substructures likely still exist within the beam and hence the observed line temperatures may be in fact smaller than the actual temperatures. The CS (3-2), observed to have an integrated line intensity of $\sim$ 11.5 K km s$^{-1}$, is reasonably well matched by Models 8, 24, and to a lesser extent, 32. The first two models simulate very dense gas (10$^6$ cm$^{-3}$) and only differ in the strength of the radiation field, while the latter simulates a lower density but cosmic ray dominated gas. From this transition and single dish observations alone, we cannot distinguish between the two types of gas. On the other hand, for the CS (5-4) transition with an observed integrated line intensity of $\sim$ 4.8 K km s$^{-1}$, the closest (but not a good fit) match is achieved by Models 49, 7 and 23; these all have a very high density gas but differ from those matching the lower $J$ transition either in visual extinction or temperature. The HCO$^+$ (3-2) transition, observed to have an intensity of $\sim$ 4.2 K kms$^{-1}$, is instead best matched by Model 44, simulating a low density (10$^4$ cm$^{-3}$), high temperature (200 K), and a low visual extinction gas, as one would expect as ions do not survive in very optically thick gas unless the ionization rates are enhanced. In fact, Model 44 best matches the HCO$^+$ (3-2) transition at earlier times.

Of course, molecular ratios of lines are routinely used. This however makes the assumption that the two transitions arise from the same gas and that the individual column densities are both correct. For example, the CS(5-4)/CS(3-2) ratio is best matched by Models 13, and then 6 and 22; none of these models match either of the individual line intensities.

Finally, comparing the CS LTE column density derived by with our best matching models we find that all models matching the CS (3-2) can reproduce it, but, of the models that best match the CS (5-4) transition, Model 49, with the highest temperature and visual extinction, matches best the column density. The column density of HCO$^+$ is very well matched by the model that fits best its line intensity.

also observed SO(4$_3$-3$_2$) and several methanol transitions. For the former its line intensity (0.94 K km s$^{-1}$) is not well matched by any of our models, although its column density is close to that we derive for Model 33, highlighting the danger of just relying on column densities comparison alone (e.g. as in ). The observed column density for methanol, $\sim$ 10$^{14}$ cm$^{-2}$, is well matched by several models, of which all but one have a gas density of 10$^5$ cm$^{-3}$; interestingly all of the shock models over-predict its column density.

In summary, while a handful of single dish transitions are clearly not enough to characterize the multi-gas components of the nucleus of a galaxy, our ab initio method is able to certainly exclude ranges of densities, temperatures, and other physical parameters that do not affect the observed gas. It is clear that single dish observations cannot be interpreted by assuming different species and transitions as unique tracers of a particular environment. However this qualitative comparison already seems to indicate that while the CS (3-2) transition traces very dense gas at $\sim$ 100 K, it is likely that the CS (5-4) traces a hotter component and the HCO$^+$ traces a hot, possibly young, gas, at a lower density. Comparisons with interferometric observations will offer stringent tests to our method, as outlined by the example below.

## NGC 253

Again, as for the case of M82, here we simply perform an exercise to see how a priori modelling can help constrain some of the characteristics of the dense gas in another commonly studied galaxy. spatially resolved several molecular species in the central kiloparsec of NGC 253 (at a distance of $\sim$ 350 pc) with ALMA at a resolution of $\sim$ 50 pc. They identified 27 molecular lines in 10 individual locations. performed a detailed analysis highlighting the rich chemistry in the central zone of this galaxy. They reported the integrated line intensities, derived column densities, and LVG analyses. Their Figure 5 shows an HCN/HCO$^+$ ratio as a function of position going from $\sim$ 0.8 to 1.6. Looking at Table 9 we find that this range is hard to match with any model and in fact those that are close to their best fit LVG model for the inner nuclear disk (positions 5–7), i.e at a gas density of n$_H$ =10$^{4.5-5.0}$ cm$^{-3}$ and a gas temperature of T = 60-120 K, all lead to a much higher HCN(1-0)/HCO$^+$1-0 ratio. This is consistent with conclusions that the main isotopologue HCN/HCO$^+$ line ratios cannot alone be used as a unique diagnostic for starburst or AGN dominated galaxies, in this case probably owing to their high opacities. On the other hand, if we look at the individual theoretical intensities from Tables 7 and 8, we find that a good match for the HCN(1-0) transition can be obtained with Models 46 and 61 (10$^5$-10$^6$ cm$^{-3}$, 50-200K) but only at early times, i.e the gas does not seem to be in equilibrium. Interestingly, the gas traced by this transition does not seem to be affected by enhanced cosmic ray ionization rates. The HCO$^+$ (1-0) transition in the same nuclear region can be matched by Models 27, 39, and 54, differing in physical characteristics but all sharing the commonality that the cosmic ray ionization rate is higher than standard; more importantly the gas is at equilibrium, implying that the gas that HCO$^+$ is tracing is not the same as that traced by HCN(1-0), which is another possible reason why their ratio is a poor indicator in this occasion.

also observes typical shock tracers: SiO and HNCO. Making the assumption that SiO traces shocks they interpret the SiO/HCN line ratio as an indicator of gas that is experiencing shocks. Yet, their SiO abundance, $\sim$ 10$^{-9}$, which depending on the visual extinction gives column densities of 2–8$\times$`<!-- -->`{=html}10$^{13}$ cm$^{-2}$, is not matched by our shock models but rather by models that have an enhanced cosmic ray ionization rate (Models 33, 36, 39, see Table 10). On the other hand, HNCO can be matched by several models, including a shock model (Model 67). Hence, in this particular case, it seems that SiO can be used as a tracer of enhanced cosmic ray ionization gas.

Speculating further on the physical characteristics of each galaxy is beyond the scope of this paper but it is clear that a comparison of several transitions and species with large grids of chemical models may lead to a more complete picture than simply comparing molecular line ratios, especially if both column densities and individual line intensities, with appropriate corrections based on the filling factor of each observation, are taken into consideration.

# Conclusions

We present a grid of chemical models and provide the abundances, abundance ratios, and column densities of the some of the most common tracers for a large parameter space in gas density, temperature, cosmic ray ionization rate, and radiation field (Tables 1–6). The theoretical abundances are then used as inputs to a radiative transfer code to derive the theoretical intensities for the parameter space investigated by the chemical models (Tables 7, 8, 10–16). We finally provide examples of comparisons of observations from two galaxies with our theoretical line intensities. Our main conclusions, which are only valid within the parameter space of chemical models that we ran here, are as follows: (i) Time dependent effects on chemical abundances are only important for models in which the cosmic ray ionization rate is close to the galactic ionization rate. However the most affected species are, in fact, those most used to trace dense gas (HCN and CS). (ii) Line intensities ratios, such as HCN(1-0)/HCO$^+$(1-0), are not unique across the grid, leading to clear degeneracies in any RADEX solution. (iii) Individual species can be unique tracers of a particular energetic process, but only within specific physical conditions. (iv) A brief analysis of molecular observations from two galaxies have shown how only relying on molecular line ratios can lead to erroneous conclusions on the physical conditions of the emitting gas and how therefore individual intensities, as well as column densities or fractional abundances, need to be taken into consideration.

[^1]: `This code has now being made public `` and can be downloaded at https://uclchem.github.io/`


author: S. N. Shore, G. M. Wahlgren
date: Received —; accepted —
title: The O I\] 1641Å line as a probe of symbiotic star winds

The neutral oxygen resonance $\lambda$`<!-- -->`{=html}1302Å line can, if the optical depth is sufficiently high, de-excite by an intercombination transition at $\lambda$`<!-- -->`{=html}1641Å to a metastable state. This has been noted in a number of previous studies but never systematically investigated as a diagnostic of the neutral red giant wind in symbiotic stars and symbiotic-like recurrent novae. We used archival $IUE$ high resolution, and GHRS and STIS medium and high resolution, spectra to study a sample of symbiotic stars. The integrated fluxes were measured, where possible, for the O I $\lambda$`<!-- -->`{=html}1302Å and O I\] $\lambda$`<!-- -->`{=html}1641Å lines. The intercombination 1641Å line is detected in a substantial number of symbiotic stars with optical depths that give column densities comparable with direct eclipse measures (EG And) and the evolution of the recurrent nova RS Oph 1985 in outburst. In four systems (EG And, Z And, V1016 Cyg, and RR Tel), we find that the O I\] variations are strongly correlated with the optical light curve and outburst activity. This transition can also be important for the study of a wide variety of sources in which an ionization-bounded H II region is imbedded in an extensive neutral medium, including active galactic nuclei, and not only for evaluations of extinction.

# Introduction

Symbiotic stars present the unusual situation of a nearly neutral, stable environment centered on a cool giant star, in which a hot source, along with its surrounding ionized region, is imbedded. The radius of the H II region is determined only by the mass gainer's effective temperature and luminosity, which in turn depend only on the accretion rate from the wind (or in the cases where a disk forms, from the flux distribution of the surrounding disk along with that of the underlying star). Since these can be separated using multiwavelength observations, and the incident spectra are simple at ultraviolet (UV) wavelengths (a hot white dwarf and/or an accretion disk continuum and emission line continuum), it is possible to model the formation of the spectrum comparatively easily. This is mainly because the wind of the companion red giant is at nearly its terminal velocity (see Vogel 1991, Pereira et al 1999) and, even if structured by the orbital motion and hydrodynamic processes related to the accretion (e.g. Dumm et al. 2000; Walder et al. 2008) this happens on a length scale far larger than the gainer and its ionized zone.

In such an environment several radiative processes, not usually encountered under nebular conditions, are observable. Principal among these are fluorescence due to various scattering mechanisms. Accidental resonances account for much of the down-conversion of UV emission to emission in the optical and near infrared. Perhaps the best known are those Fe II and related ions that can be excited by UV resonance transitions of highly-ionized species, e.g. C IV and its coincidence with ground state multiplets of Fe II that de-excite through optical forbidden transitions (Johansson 1983, 1988). Raman scattering (e.g. Schmid 1989), a nearly coincident resonance process that produces broad, down-converted emission lines, is particularly spectacular in the symbiotics, the most notable lines being those of O VI $\lambda\lambda$`<!-- -->`{=html}6825, 7082Å that are produced by the near coincidence of the resonant O VI doublet $\lambda\lambda$`<!-- -->`{=html}1031, 1037Å and H Ly$\beta$. There is, in addition, a process whereby the UV resonance line of a *neutral* species can, by virtue of absorption in a surrounding neutral gas, produce both optical and UV emission lines through otherwise inaccessible forbidden transitions. This happens because the ionization potential of several neutral atoms, in particular oxygen, is slightly higher than that of hydrogen and can therefore form in the H II region along with those formed by recombination. Furthermore, resonant absorption by the neutral gas at energies significantly below the ionization limit can, if the optical depth is sufficiently large, lead to emission in alternate channels even in the resonant scattering case.

The O I spectrum is a case in point. The $\lambda$`<!-- -->`{=html}1302 resonance line is one of the strongest emission features in the spectrum of late-type symbiotics. It forms in the H II region around the degenerate gainer since O I has a slightly higher ionization potential than neutral hydrogen. In addition to the ground state, the O I $^3$S$\rm^o$ $-$ $^3$P ($\lambda\lambda$`<!-- -->`{=html}1302, 1304, 1305Å) multiplet 5 is connected to two long-lived states through emission at $\lambda$`<!-- -->`{=html}1641Å and $\lambda$`<!-- -->`{=html}2324Å, both spin forbidden (intercombination) transitions (see Fig. 1), and their associated decay channels to the ground state. One such decay channel, $\lambda$`<!-- -->`{=html}6300 Å, is well known from terrestrial auroral spectra. Also, the $\lambda$`<!-- -->`{=html}1641 line has been used as a proxy measure of solar activity variability and its effect on the atmosphere (e.g. Bowers et al. 1987, see below). These lines are also well known from planetary nebulae (e.g. Feibelman 1997) and have been discussed in the literature for studies of interstellar extinction in Seyfert galaxies (Grandi 1983) and the determination of the oxygen abundance in cool stars (\[O I\] $\lambda$`<!-- -->`{=html}6300, Nissen et al. 2002).

A difficulty presented by any neutral or singly-ionized resonance transitions is that the interstellar medium, possessing the same resonance transitions, is opaque along many lines of sight, especially for distances of several kiloparsecs that are typical of symbiotic and planetary nebular targets. This is exacerbated for cosmological distances where the intervening Ly$\alpha$ forest potentially contaminates the whole redshift range from that of the host galaxy to nearly the local standard of rest. These systems should, therefore, present sufficient line of sight optical depths to produce detectable O I\] emission.

In a study of the UV spectra of the recurrent nova RS Oph during its 1985 outburst, Shore & Aufdenberg (1991) noted the presence of a transient emission line on the red wing of He II $\lambda$`<!-- -->`{=html}1640 relatively early in the outburst and identified this as O I\] $\lambda$`<!-- -->`{=html}1641. This line was also identified by Aufdenberg (1993) in the $IUE$ spectrum of RR Tel. In a recent study of the 2006-2009 outburst of the S-type symbiotic star AG Dra, we discussed the variations of the optical spectra, concentrating on the optical Raman features (Shore, Wahlgren, Genovali, et al. 2010). This survey included an examination of archival material as well as optical high-resolution spectra. The absence of the \[O I\] $\lambda$`<!-- -->`{=html}6300 line was noted but it was suggested that it would be worthwhile checking the existing archive of high resolution UV data for O I\] $\lambda$`<!-- -->`{=html}1641. This symbiotic is a special case, having a radial velocity of $-$`<!-- -->`{=html}144 km s$^{-1}$; any resonance line originating from the star is well shifted in wavelength with respect to its ISM components. In this paper we report on our search of the archives for the presence of O I\] $\lambda$`<!-- -->`{=html}1641, as well as other O I lines, in the spectra of symbiotic stars. As mentioned above, the O I spectrum originates in the vicinity of the red-giant star. Observations of the O I\] line may prove to be useful diagnostics of the red-giant wind and its sources of excitation. Correlation of the UV lines with optical and near-infrared (near-IR) O I lines would therefore enable studies of symbiotic star properties and behavior in the absence of UV spectra.

# Observations

We retrieved all MAST archival spectra for symbiotic stars taken with the *International Ultraviolet Explorer* (*IUE*) satellite at high resolution (R $\approx$ 10000, large aperture) and *Hubble Space Telescope* ($HST$) GHRS and STIS medium resolution (G140M, G160M) spectra. No $HST$ echelle spectra exist for symbiotic stars. For AG Dra, these were supplemented with Telescopio Nazionale Galileo (TNG) high resolution optical spectra. Note was taken of literature sources presenting symbiotic spectra that included O I lines. Figure 2 presents a sample of O I\] $\lambda$`<!-- -->`{=html}1641 + He II $\lambda$`<!-- -->`{=html}1640 line profiles for six symbiotics from $HST$/STIS spectra. The high spectral resolution of the STIS instrument clearly shows that the O I\] $\lambda$`<!-- -->`{=html}1641 line can be located either within or outside of the He II $\lambda$`<!-- -->`{=html}1640 profile.

## Stellar sample

The available data sets are not of a homogeneous quality, since they were obtained for various scientific purposes, as well as serving as a calibration target. In addition, these data neither represent a thorough nor even sampling of light curves or eruptive events. Therefore, statistics and correlations are directed to the observed occurrence of O I\] $\lambda$`<!-- -->`{=html}1641 and its possible excitation mechanism. Stars for which the line was detected or suspected of being present in the available spectra are the following:

-   Z And: The line is visible in high resolution $IUE$ spectra, its variations are discussed in the following section. The system is both eclipsing (as for CH Cyg) and a very active variable with jet-like outflows having been detected.

-   EG And: Two GHRS spectra show the He II and O I\] lines. The forbidden line is strong and easily visible on the lower resolution spectra. There are five STIS pointings, see Crowley et al. (2008) for details. They note the presence of the line but do not study its variations relative to the orbital phase. The main point is that there is almost no variability in the STIS spectra while the He II line is strongly Fe-curtained. The line is present in all spectra taken outside of eclipse, its variation, based on the $IUE$ data, is discussed in the next section.

-   CH Cyg: The UV absorption Fe-curtain spectrum is among the strongest of any symbiotic star. Although cited by Hack & Selvelli (1982), the O I\] line is weak when present, blended with He II, and visible in two IUE high resolution spectra, SWP8940 (MJD 44365) and SWP10878 (MJD 44596).

-   CI Cyg: There is one GHRS spectrum, showing one of the strongest and most unusual O I\] lines; the intensity relative to He II is very high. There is one published study (Mikolajewska et al. 2006) that discusses the O I\], but it does not discuss the line formation. Possibly present in $IUE$ spectra during the mid-1990s.

-   V1016 Cyg: For the high resolution $IUE$ spectra, many are saturated at He II and partially mask any weaker emission lines, but the O I\] line is apparent in spectrum SWP05612 and possibility detected in other spectra. Only a single G140L STIS spectrum is available from the archives, the resolution of which is insufficient to detect the O I\]. The line variations are discussed in the next section.

-   V1329 Cyg: The O I\] line is possibly present at a level just above the approximate continuum noise level for several $IUE$ spectra, especially in SWP29816.

-   AG Dra: O I\] is strong, one STIS spectrum shows this perfectly. Weakly present in the red wing of He II $\lambda$`<!-- -->`{=html}1640 in the $IUE$ spectrum SWP25444. The O I\] line is, however, too weak to measure in the archival $IUE$ spectra and there is only one *STIS* spectrum (2003 Apr. 19, see Shore et al. 2010), for which the O I\] $\lambda$`<!-- -->`{=html}1641 flux is $8.13\times 10^{-14}$ erg s$^{-1}$cm$^{-2}$ and the (O I $\lambda$`<!-- -->`{=html}1302/O I\] $\lambda$`<!-- -->`{=html}1641) ratio is 10.2. Mikolajewska et al. (1995) propose a moderately small inclination and no eclipses.

-   RW Hya: Two GHRS and one STIS low resolution spectra show only a hint of O I\], mainly as a redward-extended wing on He II that could mask the unresolved line.

-   SY Mus: High resolution $IUE$ spectra (including SWP14236) show the O I\] line.

-   AG Peg: Two GHRS spectra give a possible indication of the line, but it is clearly weak. Possibly present in the $IUE$ spectrum SWP37420, the shortest exposure available.

-   RX Pup: There is one STIS spectrum but no detectable emission at O I\]. An $IUE$ spectrum (SWP14240) shows the line.

-   HM Sge: Only high resolution $IUE$ are available, and there is a possibility that these show the O I\] line.

-   RR Tel: The most studied of the symbiotic stars in this sample, spectra were taken with GHRS and STIS as a standard for calibration of wavelengths and comparisons between instruments. The O I\] line is very strong. The line is also strong in $IUE$ spectra.

-   KX TrA: Weakly present in the $IUE$ spectra SWP38741, SWP38742.

## Notes and correlations

For the stars AS210, S190, R Aqr, AE Ara, T CrB, BF Cyg, and AX Per the O I\] line was either absent or too weak to measure in the available data, in all cases only $IUE$ spectra exist. Based on the available spectra mentioned, two-thirds of the symbiotic stars have shown the presence of O I\] $\lambda$`<!-- -->`{=html}1641. Among the remaining classical symbiotics from Kenyon (1986), Allen (1984), and Belcyński et al (2000), Y CrA, V443 Her, BX Mon, AR Pav, CL Sco, HK Sco, CL Sco, and AS296 have only low resolution *IUE* archival spectra, and are therefore not useful for detecting O I\] $\lambda$`<!-- -->`{=html}1641.

Although not measured in the previous studies (especially the outburst analysis by Shore et al. (1996)), the O I\] line is detected in *IUE* high dispersion spectra of the 1985 outburst of the symbiotic-like recurrent nova RS Oph. In Fig. 3 we show three early spectra, from 10 days (SWP 25248), 25 days (SWP 25290), and 30 days (SWP 25328) after optical maximum. The O I\] line is clearly detected in only one of these, SWP 25290, which was obtained during the broad-line phase of the permitted lines and when the inferred neutral column density was about $10^{23}$cm$^{-2}$ based on the Fe-curtain absorption (derived from the narrow line components of C IV, see Shore et al. (1996)). Its radial velocity is consistent with that of the binary system, $v_{\rm rad} = -10$ km s$^{-1}$. To our knowledge, the line O I $\lambda$`<!-- -->`{=html}1.12895 $\mu$m has only been identified in symbiotic and related star spectra by Evans et al. (2007), where it was recorded in post outburst spectra of RS Oph.

We show four examples of long term variability of this line in Fig. 4 for EG And, Fig. 5 for V1016 Cyg, Fig.6 and Fig.7 for RR Tel, and Fig. 8 for Z And. The precipitous drop in the RR Tel flux corresponds to a step in the visual magnitude light curve, obtained from the AAVSO website. Similarly, the spike in the Z And O I\] flux corresponds to a spike in the visual magnitude light curve. The AAVSO visual magnitude light curve for V1016 Cyg shows a weak correlation with peaks in the O I\] flux in Fig. 5. The sparser EG And data in Fig. 4 is more difficult to correlate with its light curve; however, the amplitude of the O I\] flux is similar to the visual magnitude amplitude over a period that is similar to the 482 day period of its (orbit/pulsation). This latter is an eclipsing system in which the line of sight optical depth is so large, due to the absorption line and Rayleigh scattering opacities, as to obscure the entire region, including the O I\] line. For non-eclipsing systems, the optical depth appears never to be so large that the O I\], which is always optically thin, is not seen. As with the Raman feature, the line is formed close to the ionized-neutral interface in the red giant wind and for low inclination systems should always be visible.

Parimucha et al. (2002) have shown that the fluxes for resonance lines of ionized species (C III-C IV, N III-N V, etc.) in V1016 Cyg showed a well defined minimum at around MJD 47000. It is possible that the O I\] line, which shows an almost lightcurve, may reach minimum strength slightly earlier.

The secular development of the O I\] strength in RR Tel is the same as the slow decline of the $V$ magnitude, according to the AAVSO light curves, for the entire period of the $IUE$ measurements. In addition, since in this system the stellar radial velocity suffices to displace the O I $\lambda$`<!-- -->`{=html}1302Å line from within the interstellar absorption, it is possible to study the long term variation of the $\lambda$`<!-- -->`{=html}1302Å to $\lambda$`<!-- -->`{=html}1641Å line flux ratio, shown in Fig. 7. There is an apparently asymptotic trend with F($\lambda$`<!-- -->`{=html}1302)/F($\lambda$`<!-- -->`{=html}1641)$\approx$ 0.4 at late times in the IUE data set. This flux ratio agrees with the STIS observations from almost a decade later.

The integrated flux variations of the $\lambda$`<!-- -->`{=html}1641Å line for Z And (Fig. 8) shows a strong correlation with the long term optical variations (based on the AAVSO archive), especially the strong outburst between MJD 46000 and 47000. In the figure, the zero flux values are the conjunctions when it is heavily obscured (according to Friedjung et al. 2010) that were otherwise well exposed (not, as in several symbiotic stars, over- or underexposures that we have ignored in the analysis).

The *HST* data set includes extended wavelength coverage at high spectral resolution and signal-to-noise for the stars AG Dra, EG And, and RR Tel. This allows for a search of multiple lines from O I for the purpose of investigating emission line excitation mechanisms. STIS spectra for AG Dra extend from the vacuum UV to the red. Emission is observed for $\lambda\lambda$`<!-- -->`{=html}1302, 1304, 1305, 1355, 1358, 1641. The $\lambda$`<!-- -->`{=html}1304 line is nearly entirely removed by saturated absorption from the ISM feature Si II $\lambda$`<!-- -->`{=html}1304.370, while the O I $\lambda$`<!-- -->`{=html}1302 emission line is impinged upon by its saturated ISM counterpart (as an example we show the 1300Å region of AG Dra in Fig. 9). Several lines ($\lambda\lambda$`<!-- -->`{=html}2324, 8446, 9204, 9260) may have marginal detections, and no evidence exists for the presence of others ($\lambda\lambda$`<!-- -->`{=html}5577, 6300, 6363, 6391). The $\lambda$`<!-- -->`{=html}2324 line is suspect due to broadening of the 2324 Å feature (C II + O I).

For EG And the $HST$ data are limited to UV wavelengths, showing emission for $\lambda\lambda$`<!-- -->`{=html}1302,1304,1305,1355, and 1358. Emission in O I $\lambda$`<!-- -->`{=html}2324.738, if present, is blended with the emission feature C II $\lambda$`<!-- -->`{=html}2324.69, which is one of five lines that comprise the C II UV2 2s$^2$ 2p $^2$P$^o$ - 2s 2p$^2$ $^4$P multiplet, all of which are found in emission. O I $\lambda$`<!-- -->`{=html}1727.106 does not appear to be present. Emission is not found for lines at the longer wavelengths ($\lambda\lambda$`<!-- -->`{=html}2958, 2972) as the continuum flux of the cool star, the increased number of absorption lines, and the lower transition probabilities render these lines difficult to detect. The RR Tel data set clearly shows emission from O I $\lambda\lambda$`<!-- -->`{=html}1302, 1304, 1305, 1641, 5577, 6300, 6363. The detection of $\lambda$`<!-- -->`{=html}2324 line is complicated by blending with C II, with the lines of this multiplet not being in proportion to their relative $gA$ values. Strong emission at $\lambda$`<!-- -->`{=html}1025 is not evident (see below).

The O I\] $\lambda$`<!-- -->`{=html}1641 emission appears to be positively correlated with the O VI Raman emission at $\lambda\lambda$`<!-- -->`{=html}6825, 7082. Based on ORFEUS spectra, Schmid et al. (1999) compared emission from the O VI $\lambda\lambda$`<!-- -->`{=html}1032, 1038 Å resonance doublet with the Raman scattering features for six symbiotic stars. Of these six, five have UV spectra that can be searched for O I\]. We find that the strongest O VI Raman emitters (AG Dra, RR Tel) are also the strongest O I\] emitters, while three others (AG Peg, Z And, V1016 Cyg) present both weak O I\] and weak O VI Raman lines. The O VI FUV doublet has two contributors, the white dwarf wind and surrounding ionized cavity within the red giant wind, that are Raman scattered by Ly$\beta$ absorption by the H I in the wind. The O I\] line samples the same warm neutral H wind zone. Crowley et al. (2008), for instance, find that this lies about 10$^{13}$cm (a small fraction of the semimajor axis of the system) from the white dwarf and it is in the cooler zone that both oxygen and hydrogen are neutrals. Since the O I 1302Å line arises from recombination in the H II region, whose boundary is virtually the same for the two species, the 1641Å line probes the same region as the Raman feature. During outburst events, for instance for Z And around MJD 46700 (Fig. 7), the O I\] line notably and briefly strengthened and when the He II line was hidden by a sufficient optical depth in the Fe-curtain absorption the O I\] was also obscured.

# Discussion

Several fluorescence processes can lead to the population of levels that will ultimately lead to emission of O I\] 1641 Å. Coincidence of H Ly$\beta$ $\lambda$`<!-- -->`{=html}1025.722 with O I $\lambda$`<!-- -->`{=html}1025.762 will populate the O I 3d $^3$D$^o$ (97488 cm$^{-1}$) level from the ground level. The dominant decay chain (according to their Einstein transition probabilities, see Fig. 1) from the 3d $^3$D$^o$ level ($\lambda\lambda$`<!-- -->`{=html}11285 and 8446) will populate the 3s $^3$S$^o$ (76794 cm$^{-1}$) level, which subsequently decays through three channels ($\lambda\lambda$`<!-- -->`{=html}1302, 1641, 2324), two of which are commonly detected in symbiotic star spectra. In addition to the Ly$\beta$ O I pumping, McMurry & Jordan (2000) identified CO emission fluorescently-excited by O I UV 2 resonance line emission near $\lambda$`<!-- -->`{=html}1302 in the UV spectrum of $\alpha$ Tau.

A second pumping mechanism for the O I 76794 cm$^{-1}$ level is He II $\lambda$`<!-- -->`{=html}1640 for those stars which have a broad He II line. This is evident for RR Tel and EG And, and less so for RW Hya and AG Peg, in Fig. 2. Other possible pumping mechanisms that might lead to population enhancement of this O I level are: 1) H Ly$\epsilon$ $\lambda$`<!-- -->`{=html}937.803 Å coincident with O I $\lambda$`<!-- -->`{=html}937.841 to pump the O I 106765 cm$^{-1}$ level, which can decay to the 76794 cm $^{-1}$ level through the chain $\lambda\lambda$`<!-- -->`{=html}14110, 4368 Å, among others. 2) H Ly$_6$ $\lambda$`<!-- -->`{=html}930.748 and He II $\lambda$`<!-- -->`{=html}930.342 can pump the O I 8s $^3$S$^o$ 107497 cm$^{-1}$ level, which decays to the 76794 cm$^{-1}$ level through the chains $\lambda\lambda$`<!-- -->`{=html}12790, 4368 or $\lambda\lambda$`<!-- -->`{=html}5298, 8446. 3) C II $\lambda$`<!-- -->`{=html}2324.69 emission is coincident with O I $\lambda$`<!-- -->`{=html}2324.738 and can pump the O I 76794 cm$^{-1}$ level from the metastable O I level 33792 cm$^{-1}$. The five lines of C II multiplet 2 are seen in emission in a number of symbiotics and symbiotic novae (RR Tel). Direct excitation by the resonance line should, however, be more effective in symbiotics, as in the terrestrial case since the Ly lines are so optically thick and the illumination is from the companion, not *in situ* from the chromosphere (there will, of course, be a contribution from the spectrum of the late-type component but this is small compared to that from white dwarf environment).

Population of the O I 76794 cm$^{-1}$ level via electron recombination is possible through additional decay chains. Spectral observations at infrared wavelengths may offer a means of determining the dominant excitation mechanisms by detecting other emission lines. The number of lines from the O I spectrum that have been observed in astronomical targets, in particular symbiotic stars and novae, are few. Common UV lines detected include transitions at wavelengths $\lambda\lambda$`<!-- -->`{=html}1302, 1304, 1305, 1355, 1358, 1641. At optical wavelengths $\lambda\lambda$`<!-- -->`{=html}6300, 6363 are found in planetary nebula spectra with $\lambda$`<!-- -->`{=html}6300 commonly used for abundance analysis in cool stars. For near-IR wavelengths, detections, or suspicions of detections, have been mentioned for 8446 Å in AG Dra (Iijima et al. 1987), and $\lambda$`<!-- -->`{=html}11289 (Evans et al. 2007). Absorption lines at $\lambda\lambda$`<!-- -->`{=html}7771, 7773, 7774 are commonly used in abundance analysis in a variety of stars. There is also the curious appearance of an undiscussed weak emission line near 2.9 $\mu$m (Schild, Boyle & Schmid 1992) in spectra of several symbiotics. Conspicuous by their absence from discussion and published spectra are lines of large transition probability ($\lambda\lambda$`<!-- -->`{=html}4368, 9204, 9260 for example) and small transition probability ($\lambda\lambda$`<!-- -->`{=html}1727, 2324, 2958, 2972). A full accounting of O I lines for any target would be useful for determining the excitation conditions and better enable the physical modeling.

The importance of the O I\] $\lambda$`<!-- -->`{=html}1641 line for the symbiotics is as a possible tool for as long as ultraviolet spectroscopy is available. Oddly, although this line has been included in a number of identification lists at high resolution, it has never been exploited as a diagnostic for symbiotics or related systems. It has, however, been noted as a contributor to the energetics of AGN when the O I resonance line is sufficiently optically thick. Grandi (1983), in discussing reddening determinations for AGN using the resonance and Bowen fluorescence O I UV2 lines ($\lambda\lambda$`<!-- -->`{=html}1302, 1304 vs $\lambda$`<!-- -->`{=html}8446) noted that the line ratio O I\] $\lambda$`<!-- -->`{=html}1641 (UV146) to the resonance multiplet is often unusually large, given the branching ratio. This can be accounted for by a large enough optical depth to strongly self-absorb the ground state lines. With a Ly$\alpha$ optical depth as large as $10^6$ the reduction in $\lambda$`<!-- -->`{=html}1302 is sufficient to produce an integrated flux of only a factor of 2 larger than the forbidden transition. The inhomogeneous regions around the central engine often show such large opacities while still permitting observation of the nucleus along a given sight line. More recently, the chromospheric O I spectrum has been rediscussed for a few main-sequence and evolved F, G, and K stars by Koncewicz (2005) and Koncewicz & Jordan (2007). There is another mode to produce the O I\] $\lambda$`<!-- -->`{=html}1641 emission, the coincidence of Ly$\beta$ $\lambda$`<!-- -->`{=html}1025.72 and the O I resonance line (UV4) at $\lambda$`<!-- -->`{=html}1025.77 that pumps the $^3$D$^o$ 97488 cm$^{-1}$ level , which then decays through the $\lambda\lambda$`<!-- -->`{=html}11286, 8446, 1641, 6300 chain.

To date, however, most of the literature deals with the O I lines in the context of planetary – specifically, terrestrial – atmospheric structure and composition. Atomic oxygen forms in excited states by dissociative collisions between O$_2$ and electrons. These transitions have been used for studying the oxygen abundance and temperature structure of the troposphere in a number of papers, e.g. Meier & Conway (1985), and Conway et al. (1988). Doering & Gulcicek (1989) include the $\lambda\lambda$`<!-- -->`{=html}1355, 1358 lines. Since the O I\] transition is always optically thin and absent in the reflected solar spectrum this transition probes almost the entire terrestrial stratosphere and ionosphere. The branching ratio (Garstang 1961, Erdman & Zipf 1986) is O I\] $\lambda$`<!-- -->`{=html}1641/O I $\lambda$`<!-- -->`{=html}1302 = 5.1$\times$`<!-- -->`{=html}10$^{-6}$ with an uncertainty of 30% . The most recent compilation, Wiese et al. (1996), gives $A$($\lambda$`<!-- -->`{=html}1641)/$A$($\lambda$`<!-- -->`{=html}1302) = 5.4$\times$`<!-- -->`{=html}10$^{-6}$. Following Grandi (1983), based on escape probability formalism (Kwan & Krolik 1981), we can estimate the required column density in the resonance line. The observed branching ratio is $\approx$ 1 for all systems in which the O I\] 1641Å line is detected in our survey, the implied optical depth for $\lambda$`<!-- -->`{=html}1302 is $\tau_{\lambda1302} \approx 7\times 10^{-14}f v_{50}^{-1}
N_O \approx 2\times10^5$, where $f$ is the oscillator strength and $v_{50}$ is the wind velocity in km s$^{-1}$, that for a solar O/H ratio ($5\times 10^{-4}$, see Asplund et al. 2009) gives a column density N$_H$ $\approx$ 10$^{23}$ cm$^{-2}$ for the neutral absorption region. This is the same order of magnitude as the column density in absorption required to explain the narrow UV emission line variations during the early RS Oph outburst (Shore et al. 1996) and similar to that derived by Crowley et al. (2008) from eclipse spectra of EG And. Using the length scale from the photoionization modeling in Crowley et al. (2008), who obtain a standoff distance for the neutral region from the white dwarf in EG And of about 10$^{13}$ cm, gives a characteristic number density of about 10$^{10}$ cm$^{-3}$. For a wind velocity of 50 km s$^{-1}$ and using R$_{13}$ =(R/10$^{13}$ cm) gives an estimate of the mass loss rate for the red giant of $\approx$ 10$^{-6}$ R$^2$$_{13}$ M$_{\sun}$ yr$^{-1}$. This estimate is different between systems (there are several with lower branching ratios, others with higher, and in many cases the stellar radial velocity is not sufficient to displace the $\lambda$`<!-- -->`{=html}1302 line from within the interstellar absorption. The optical depth expected for the red-giant wind in the FUV O I doublet, combined with the opacity of the Ly$\beta$ transition, suggest that this is not a dominant mechanism in producing the $\lambda$`<!-- -->`{=html}1641 line and that the O I opacity suffices.