Aiming for a new and more comprehensive DIB catalog between 4,000 and 9,000{\AA}, we revisited the Atlas Catalog based on the observations of HD 183143 and HD 204827 (Hobbs et al. 2008 and 2009). Twenty-five medium-to-highly reddened sight lines were selected, sampling a variety of spectral types of the background star and the interstellar environments. The median SNR of these spectra is ~ 1,300 around 6,400{\AA}. Compared to the Atlas Catalog, 22 new DIBs were found, and the boundaries of 27 (sets) of DIBs were adjusted, resulting in an updated catalog containing 559 DIBs that we refer to as the Apache Point Observatory Catalog of Optical Diffuse Interstellar Bands. Measurements were then made based on this catalog. We found our survey most sensitive between 5,500 and 7,000{\AA} due largely to the local SNR of the spectra, the relative absence of interfering stellar lines, and the weakness of telluric residuals. For our data sample, the number of DIBs detected in a given sight line is mostly dependent on EB-V and less on the spectral type of the background star. Some dependence on fH2 is observed, but less well-determined due to the limited size of the data sample. The variations of the wavelengths of each DIB in different sight lines are generally larger than those of the known interstellar lines CH+, CH, and K I. Those variations could be due to the inherent error in the measurement, or to differences in the velocity components among sight lines.

We study the behavior of eight diffuse interstellar bands (DIBs) in different interstellar environments, as characterized by the fraction of hydrogen in molecular form [$f$(H$_2$)], with comparisons to the corresponding behavior of various known atomic and molecular species. The equivalent widths of the five "normal" DIBs ($\lambda\lambda$5780.5, 5797.1, 6196.0, 6283.8, and 6613.6), normalized to $E(B-V)$, show a "Lambda-shaped" behavior: they increase at low $f$(H$_2$), peak at $f$(H$_2$) ~ 0.3, and then decrease. The similarly normalized column densities of Ca, Ca$^+$, Ti$^+$, and CH$^+$ also decline for $f$(H$_2$) > 0.3. In contrast, the normalized column densities of Na, K, CH, CN, and CO increase monotonically with $f$(H$_2$), and the trends exhibited by the three C$_2$ DIBs ($\lambda\lambda$4726.8, 4963.9, and 4984.8) lie between those two general behaviors. These trends with $f$(H$_2$) are accompanied by cosmic scatter, the dispersion at any given $f$(H$_2$) being significantly larger than the individual errors of measurement. The Lambda-shaped trends suggest the balance between creation and destruction of the DIB carriers differs dramatically between diffuse atomic and diffuse molecular clouds; additional processes besides ionization and shielding are needed to explain those observed trends. Except for several special cases, the highest $W$(5780)/$W$(5797) ratios, characterizing the so-called "sigma-zeta effect", occur only at $f$(H$_2$) < 0.2. We propose a sequence of DIBs based on trends in their pair-wise strength ratios with increasing $f$(H$_2$). In order of increasing environmental density, we find the $\lambda$6283.8 and $\lambda$5780.5 DIBs, the $\lambda$6196.0 DIB, the $\lambda$6613.6 DIB, the $\lambda$5797.1 DIB, and the C$_2$ DIBs.

Absorption spectra toward Herschel 36 for the A^1Pi <-- X^1Sigma transitions of CH+ in the J=1 excited rotational level and the A^2Delta <-- X^2Pi transition of CH in the J=3/2 excited fine structure level have been analyzed. These excited levels are above their ground levels by 40.1 K and ~25.7 K and indicate high radiative temperatures of the environment, 14.6 K and 6.7 K, respectively. The effect of the high radiative temperature is more spectacular in some diffuse interstellar bands (DIBs) observed toward Her 36; remarkable extended tails toward red (ETR) were observed. We interpret these ETRs as due to a small decrease of rotational constants upon excitation of excited electronic states. Along with radiative pumping of a great many high-J rotational levels, this causes the ETRs. In order to study this effect quantitatively, we have developed a model calculation in which the effects of collision and radiation are treated simultaneously. The simplest case of linear molecules is considered. It has been found that the ETR is reproduced if the fraction of the variation of the rotational constant, beta = (B'-B)/B, is sufficiently high (3-5%) and the radiative temperature is high (T_r > 50 K). Although modeling for general molecules is beyond the scope of this paper, the results indicate that the prototypical DIBs at 5780.5, 5797.1, and 6613.6 A which show the pronounced ETRs are due to polar molecules sensitive to the radiative excitation. The requirement of high beta favors relatively small molecules with 3-6 heavy atoms. DIBs at 5849.8, 6196.0, and 6379.3 A which do not show the pronounced ETRs are likely due to non-polar molecules or large polar molecules with small beta.

Anomalously broad diffuse interstellar bands (DIBs) at 5780.5, 5797.1, 6196.0, and 6613.6 A are found in absorption along the line of sight to Herschel 36, the star illuminating the bright Hourglass region of the H II region Messier 8. Interstellar absorption from excited CH+ in the J=1 level and from excited CH in the J=3/2 level is also seen. To our knowledge, neither those excited molecular lines nor such strongly extended DIBs have previously been seen in absorption from interstellar gas. These unusual features appear to arise in a small region near Herschel 36 which contains most of the neutral interstellar material in the sight line. The CH+ and CH in that region are radiatively excited by strong far-IR radiation from the adjacent infrared source Her 36 SE. Similarly, the broadening of the DIBs toward Herschel 36 may be due to radiative pumping of closely spaced high-J rotational levels of relatively small, polar carrier molecules. If this picture of excited rotational states for the DIB carriers is correct and applicable to most DIBs, the 2.7 degree cosmic microwave background may set the minimum widths (about 0.35 A) of known DIBs, with molecular processes and/or local radiation fields producing the larger widths found for the broader DIBs. Despite the intense local UV radiation field within the cluster NGC 6530, no previously undetected DIBs stronger than 10 mA in equivalent width are found in the optical spectrum of Herschel 36, suggesting that neither dissociation nor ionization of the carriers of the known DIBs by this intense field creates new carriers with easily detectable DIB-like features. Possibly related profile anomalies for several other DIBs are noted.

In this paper we analyse a set of 33 optical spectra, which were acquired with the ARCES echelle spectrograph (R = 38,000) on the 3.5-m telescope at the Apache Point Observatory. We examine the H{\alpha} profile in each of these observations in order to determine the geometry and velocity structure of the previously discovered bipolar jet, which originates from the secondary star of HD 44179 located at the centre of the Red Rectangle nebula. Using a 3D geometric model we are able to determine the orbital coverage during which the jet occults the primary star. During the occultation, part of the H{\alpha} line profile appears in absorption. The velocity structure of the jet was determined by modelling the absorption line profile using the Sobolev approximation for each orbital phase during which we have observations. The results indicate the presence of a wide angle jet, likely responsible for observed biconical structure of the outer nebula. Furthermore, we were able to determine a likely velocity structure and rule out several others. We find that the jet is comprised of low-density, high-velocity, central region and a higher-density, lower-velocity, conical shell.

In this paper we examine the profiles of the complex Na I D-lines in the Red Rectangle. The spectra were acquired with the ARCES \'{e}chelle spectrograph R = 38,000 on the 3.5-m telescope at the Apache Point Observatory. Additional spectra taken with STIS were acquired from the Hubble Legacy Archive (HLA) and were used to independently confirm the spatial origin of the spectral features. The profile of a single D-line consists of double-peaked emission, red-shifted absorption and blue-shifted absorption. We find that the double-peaked emission originates from the bipolar outflow, the red-shifted absorption feature is due to the photospheric line, and the blue-shifted absorption arises from the bipolar outflow as seen against the photosphere of the luminous post-AGB component in HD 44179. In order to better understand the Na I D-line profile, we examined the periodically variable asymmetric photospheric absorption lines. The asymmetric lines are interpreted as a signature of slow self-accretion following enhanced mass-loss around periastron. An empirical model was constructed to remove the photospheric component from the Na I D-line profile in order to study the nebular emission and absorption of sodium along the line-of-sight to the primary. This paper also discusses the different origins of the single-peaked emission, the double-peaked emission and the blue-shifted and red-shifted absorption components.

We establish correlations between equivalent widths of eight diffuse interstellar bands (DIBs), and examine their correlations with atomic hydrogen, molecular hydrogen, and EB-V . The DIBs are centered at \lambda\lambda 5780.5, 6204.5, 6283.8, 6196.0, 6613.6, 5705.1, 5797.1, and 5487.7, in decreasing order of Pearson\^as correlation coefficient with N(H) (here defined as the column density of neutral hydrogen), ranging from 0.96 to 0.82. We find the equivalent width of \lambda 5780.5 is better correlated with column densities of H than with E(B-V) or H2, confirming earlier results based on smaller datasets. We show the same is true for six of the seven other DIBs presented here. Despite this similarity, the eight strong DIBs chosen are not well enough correlated with each other to suggest they come from the same carrier. We further conclude that these eight DIBs are more likely to be associated with H than with H2, and hence are not preferentially located in the densest, most UV shielded parts of interstellar clouds. We suggest they arise from different molecules found in diffuse H regions with very little H (molecular fraction f<0.01). Of the 133 stars with available data in our study, there are three with significantly weaker \lambda 5780.5 than our mean H-5780.5 relationship, all of which are in regions of high radiation fields, as previously noted by Herbig. The correlations will be useful in deriving interstellar parameters when direct methods are not available. For instance, with care, the value of N(H) can be derived from W{\lambda}(5780.5).

We utilize high-resolution (R~60,000), high S/N (~100) spectroscopy of 17 cool Pleiades dwarfs to examine the confounding star-to-star scatter in the 6707 Li I line strengths in this young cluster. Our Pleiads, selected for their small projected rotational velocity and modest chromospheric emission, evince substantial scatter in the linestrengths of 6707 Li I feature that is absent in the 7699 K I resonance line. The Li I scatter is not correlated with that in the high-excitation 7774 O I feature, and the magnitude of the former is greater than the latter despite the larger temperature sensitivity of the O I feature. These results suggest that systematic errors in linestrength measurements due to blending, color (or color-based T_eff) errors, or line formation effects related to an overlying chromosphere are not the principal source of Li I scatter in our stars. There do exist analytic spot models that can produce the observed Li scatter without introducing scatter in the K I line strengths or the color-magnitude diagram. However, these models predict factor of >3 differences in abundances derived from the subordinate 6104 and resonance 6707 Li I features; we find no difference in the abundances determined from these two features. These analytic spot models also predict CN line strengths significantly larger than we observe in our spectra. The simplest explanation of the Li, K, CN, and photometric data is that there must be a real abundance component to the Pleiades Li dispersion. We suggest that this real abundance component is the manifestation of relic differences in erstwhile pre-main-sequence Li burning caused by effects of surface activity on stellar structure. We discuss observational predictions of these effects.

Echelle spectra of HD 183143 [B7Iae, E(B-V) = 1.27] were obtained on three nights, at a resolving power R = 38,000 and with a signal-to-noise ratio ~1000 at 6400 A in the final, combined spectrum. A catalog is presented of 414 diffuse interstellar bands (DIBs) measured between 3900 and 8100 A in this spectrum. The central wavelengths, the widths (FWHM), and the equivalent widths of nearly all of the bands are tabulated, along with the minimum uncertainties in the latter. Among the 414 bands, 135 (or 33%) were not reported in four previous, modern surveys of the DIBs in the spectra of various stars, including HD 183143. The principal result of this study is that the great majority of the bands in the catalog are very weak and fairly narrow. Typical equivalent widths amount to a few mA, and the bandwidths (FWHM) are most often near 0.7 A. No preferred wavenumber spacings among the 414 bands are identified which could provide clues to the identities of the large molecules thought to cause the DIBs. At generally comparable detection limits in both spectra, the population of DIBs observed toward HD 183143 is systematically redder, broader, and stronger than that seen toward HD 204827 (Paper II). In addition, interstellar lines of C2 molecules have not been detected toward HD 183143, while a very high value of N(C2)/E(B-V) is observed toward HD 204827. Therefore, either the abundances of the large molecules presumed to give rise to the DIBs, or the physical conditions in the absorbing clouds, or both, must differ significantly between the two cases.

The ionization rate of interstellar material by cosmic rays has been a major source of controversy, with different estimates varying by three orders of magnitude. Observational constraints of this rate have all depended on analyzing the chemistry of various molecules that are produced following cosmic-ray ionization, and in many cases these analyses contain significant uncertainties. Even in the simplest case (H3+) the derived ionization rate depends on an (uncertain) estimate of the absorption path length. In this paper we examine the feasibility of inferring the cosmic-ray ionization rate using the 10830 A absorption line of metastable helium. Observations through the diffuse clouds toward HD 183143 are presented, but yield only an upper limit on the metastable helium column density. A thorough investigation of He+ chemistry reveals that only a small fraction of He+ will recombine into the triplet state and populate the metastable level. In addition, excitation to the triplet manifold of helium by secondary electrons must be accounted for as it is the dominant mechanism which produces He* in some environments. Incorporating these various formation and destruction pathways, we derive new equations for the steady state abundance of metastable helium. Using these equations in concert with our observations, we find zeta_He < 1.2*10^-15 s^-1, an upper limit about 5 times larger than the ionization rate previously inferred for this sight line using H3+. While observations of interstellar He* are extremely difficult at present...

The proto-planetary Red Rectangle nebula is powered by HD 44179, a spectroscopic binary (P = 318 d), in which a luminous post-AGB component is the primary source of both luminosity and current mass loss. Here, we present the results of a seven-year, eight-orbit spectroscopic monitoring program of HD 44179, designed to uncover new information about the source of the Lyman/far-ultraviolet continuum in the system as well as the driving mechanism for the bipolar outflow producing the current nebula. Our observations of the H-alpha line profile around the orbital phase of superior conjunction reveal the secondary component to be the origin of the fast (max. v~560$ km s$^{-1}$) bipolar outflow in the Red Rectangle. The variation of total H-alpha flux from the central H II region with orbital phase also identifies the secondary or its surroundings as the source of the far-ultraviolet ionizing radiation in the system. The estimated mass of the secondary (~0.94 M$\sun$) and the speed of the outflow suggest that this component is a main sequence star and not a white dwarf, as previously suggested. We identify the source of the Lyman/far-ultraviolet continuum in the system as the hot, inner region (T$_{max} \ge 17,000$ K) of an accretion disk surrounding the secondary, fed by Roche lobe overflow from the post-AGB primary at a rate of about $2 - 5\times10^{-5}$ M$\sun$ yr$^{-1}$. The total luminosity of the accretion disk around the secondary is estimated to be at least 300 L$\sun$, about 5% of the luminosity of the entire system. (abridged)

New optical spectra of Ca II and Na I toward HD 72127AB provide additional evidence for both spatial and temporal variations in the complex interstellar absorption along the two sight lines; archival UV spectra yield information on the abundances, depletions, and physical conditions in the gas toward HD 72127A. Similarities in the strengths of various tracers of interstellar material in the two lines of sight suggest that the total hydrogen column densities [N(H) ~ 2.5 x 10^{20} cm^{-2}] and the depletions and ionization in the main components at low LSR velocities also are similar. Toward HD 72127A, the main components are relatively cool (T < 900 K), but with depletions resembling those found in warm, diffuse disc clouds; the generally weaker components at higher velocities have much milder depletions, more like those found in halo clouds. Several trace neutral species -- Ca I, Cr I, and Fe I -- are much stronger toward HD 72127B, however. The column density of Cr I, for example, is about 30 times the value determined toward zeta Oph (the only previous detection of that species in the ISM). Dielectronic recombination in warmer gas (T > 5000 K) may be largely responsible for the enhanced abundances of those trace neutral species toward HD 72127B. If the main components toward HD 72127AB are associated with material in the Vela SNR, the differences in abundances and physical conditions occur on scales of about 1100 AU.

Echelle spectra of the double-lined spectroscopic binary HD 204827 were obtained on five nights, at a resolving power R = 38,000 and with a S/N ratio = 750 near 6000 A in the final, combined spectrum. The stars show E(B-V) = 1.11 and spectral types near O9.5V and B0.5III. A catalog is presented of 380 diffuse interstellar bands (DIBs) measured between 3900 and 8100 A in the stars' spectrum. The central wavelengths, the widths (FWHM), and the equivalent widths of nearly all of the bands are tabulated, along with the minimum uncertainties in the latter. The reliable removal of very weak stellar lines from the catalog, and of some stellar lines from the less severe blends with DIBs, is made generally easy by the highly variable radial velocities of both stars. The principal result of this investigation is that the great majority of the bands in the catalog are very weak and relatively narrow. Typical equivalent widths amount to a few mA, and the band widths (FWHM) are most often near 0.55 A. Therefore, most of these DIBs can be detected only in spectra obtained at a resolving power and a S/N ratio at least comparable to those used here. In addition, the anomalous interstellar reddening and the very high value of the ratio N(C2)/E(B-V) seen toward HD 204827 indicate that the physical conditions in one or more of the several interstellar clouds seen in this direction differ significantly from those found toward the prototypical DIB target HD 183143, for example. Probably primarily for these reasons, 113 of the 380 bands, or 30%, were not detected in four previous, modern surveys of the DIBs seen in the spectra of stars other than HD 204827. Both the tabulated data for the DIBs and a plot of our spectrum of HD 204827 are available online at http://dibdata.org

Oxygen abundances have been derived from the near-IR, high-excitation Lambda 7774 O I triplet in high-resolution, high signal-to-noise spectra of 45 Hyades dwarfs using standard one dimensional, plane-parallel LTE models. Effective temperatures of the stellar sample range from 4319-6301 K, and the derived relative O abundances as a function of T_eff evince a trichotomous morphology. At T_eff > 6100 K, there is evidence of an increase in the O abundances with increasing T_eff, consistent with non-LTE (NLTE) predictions. At intermediate T_eff (5450 < T_eff < 6100 K), the O abundances are flat, and star-to-star values are in good agreement, having a mean value of [O/H] = +0.25 +/- 0.02; however, systematic errors at the ~0.10 dex level might exist. The O abundances for stars with T_eff < 5450 K show a striking increase with decreasing T_eff, in stark contrast to expectations and canonical NLTE calculations. The cool Hyades triplet results are compared to those recently reported for dwarfs in the Pleiades cluster and the UMa moving group; qualitative differences between the trends observed in these stellar aggregates point to a possible age-related diminution of triplet abundance trends in cool open cluster dwarfs. Correlations with age-related phenomena, i.e., chromospheric activity and photospheric spots, faculae, and/or plages, are investigated. No correlation with Ca II H+K chromospheric activity indicators is observed. Multi-component LTE ``toy'' models have been constructed in order to simulate photospheric temperature inhomogeneities that could arise from the presence of starspots, and we demonstrate that photospheric spots are a plausible source of the triplet trends among the cool dwarfs.

We present high-resolution (FWHM ~ 0.3-1.5 km/s) spectra of interstellar Ca I absorption toward 30 Galactic stars. Comparisons of the column densities of Ca I, Ca II, K I, and other species -- for individual components identified in the line profiles and also when integrated over entire lines of sight -- yield information on relative electron densities and depletions. There is no obvious relationship between the ratio N(Ca I)/N(Ca II) [equal to n_e/(Gamma/alpha_r) for photoionization equilibrium] and the fraction of hydrogen in molecular form f(H2) (often taken to be indicative of the local density n_H). For a smaller sample of sightlines for which the thermal pressure (n_H T) and local density can be estimated via analysis of the C I fine-structure excitation, the average electron density inferred from C, Na, and K (assuming photoionization equilibrium) seems to be independent of n_H and n_H T. While the n_e obtained from the ratio N(Ca I)/N(Ca II) is often significantly higher than the values derived from other elements, the patterns of relative n_e derived from different elements show both similarities and differences for different lines of sight -- suggesting that additional processes besides photoionization and radiative recombination commonly and significantly affect the ionization balance of heavy elements in diffuse IS clouds. Such additional processes may also contribute to the (apparently) larger than expected fractional ionizations (n_e/n_H) found for some lines of sight with independent determinations of n_H. In general, inclusion of ``grain-assisted'' recombination does reduce the inferred n_e, but it does not reconcile the n_e estimated from different elements. The depletion of calcium may have a much weaker dependence on density than was suggested by earlier comparisons with CH and CN.

As part of an extensive survey of diffuse interstellar bands (DIBs), we have obtained optical spectra of the moderately reddened B5V star HD 62542, which is known to have an unusual UV extinction curve of the type usually identified with dark clouds. The typically strongest of the commonly catalogued DIBs covered by the spectra -- those at 5780, 5797, 6270, 6284, and 6614 A -- are essentially absent in this line of sight, in marked contrast with other lines of sight of similar reddening. We compare the HD 62542 line of sight with others exhibiting a range of extinction properties and molecular abundances and interpret the weakness of the DIBs as an extreme case of deficient DIB formation in a dense cloud whose more diffuse outer layers have been stripped away. We comment on the challenges these observations pose for identifying the carriers of the diffuse bands.

Using the new high resolution (~ 8 km/s) echelle spectrograph on the 3.5-m telescope at the Apache Point Observatory, we have begun a high sensitivity survey of the diffuse interstellar bands in a large sample of reddened stars. Now that we are two years into this long-term survey, our sample includes over 20 reddened stars which show at least one of the DIBs that had been suggested to be caused by C7-, based on the gas phase measurement of the C7- spectrum by J. P. Maier's group. The high quality astronomical data from this larger sample of stars, along with the spectroscopic constants from the new laboratory work recently reported by Maier's group, have enabled us to examine more carefully the agreement between C7- and the DIBs. We find that none of the C7- bands matches the DIBs in wavelength or expected profile. One of the DIBs (lambda5748) attributed to C7- is actually a stellar line. The two strongest DIBs attributed to C7- (lambda 6270 and lambda4964) do not vary together in strength, indicating that they do not share the same carrier. On the whole, we find no evidence supporting the hypothesis that C7- is a carrier of the diffuse interstellar bands.

We present the results of exploratory HST observations to detect the emission lines of ions in the Beta Pic disk with the spectrograph slit placed off the star. Possible emission lines from FeII have been detected at $\sim 10^{-14}$erg cm-2 s-1 at 0.5 arcsec from Beta Pic, which would suggest a total FeII density of > 2.$10^{-2}$ cm-3 at 10 AU. This detection needs confirmation. If real, it requires a large production rate of gas and dust equivalent to the total disruption of ten bodies 30 kilometers in radius per year.

Spectra obtained with the Hubble Space Telescope GHRS are combined with high-resolution optical spectra and UV spectra from Copernicus to study the abundances and physical conditions in the diffuse interstellar clouds seen along the line of sight to the star 23 Ori. Strong low-velocity (SLV) absorption, due to cool, moderately dense neutral gas and representing about 92% of the total N(H I), is seen for various neutral and singly ionized species at +20 km/s < v < +26 km/s. Weaker low-velocity absorption, probably largely due to warmer neutral gas, is seen primarily for various singly ionized species at 0 km/s < v < +30 km/s. Absorption from a number of singly and doubly ionized species, perhaps due to a radiative shock, is seen at -108 km/s < v < -83 km/s. Weak absorption components from ionized gas are seen at intermediate velocities (-43 km/s < v < -4 km/s). The large range in n_e (from 0.04 cm^{-3} to 0.95 cm^{-3}) derived independently from nine pairs of neutral and singly ionized species in the SLV gas suggests that additional processes besides simple photoionization and radiative recombination affect the ionization balance. The adopted average SLV electron density, n_e = 0.15+/-0.05 cm^{-3}, implies a relatively large n_e/n_H ~ 0.01, and thus some ionization of hydrogen in these predominantly neutral components. Comparisons of the SLV depletions and n_H with those found for the strong ``component B'' (v ~ -14 km/s) blend toward zeta Oph hint at a possible relationship between depletion and local density for relatively cold interstellar clouds. Calcium appears to be more severely depleted in warm, low density gas than has generally been assumed. An appendix summarizes the most reliable oscillator strengths currently available for a number of the interstellar absorption lines analyzed in this work.

The Goddard High Resolution Spectrograph (GHRS) of the Hubble Space Telescope (HST) has been used to observe the boron 2500 A region of BD-13 3442. At a metallicity of [Fe/H]=-3.00, this is the most metal-poor star ever observed for B. Nearly 26 hours of exposure time resulted in a detection. Spectrum synthesis using the latest Kurucz model atmospheres yields an LTE boron abundance of log e(B)=0.01 +/- 0.20. This value is consistent with the linear relation of slope ~1.0 between log e(B[LTE]) and [Fe/H] found for 10 halo and disk stars by Duncan et al. 1997. Using the NLTE correction of Kiselman and Carlsson (1996), the NLTE boron abundance is log e(B)=0.93 +/- 0.20. This is also consistent with the NLTE relation determined by Duncan et al. (1997) where the slope of log e(B[NLTE]) vs. [Fe/H] is ~0.7. These data support a model in which most production of B and Be comes from the spallation of energetic C and O nuclei onto protons and He nuclei, probably in the vicinity of massive supernovae in star-forming regions, rather than the spallation of cosmic ray protons and alpha particles onto CNO nuclei in the general interstellar medium.