A&A 415, 331-348 (2004)
DOI: 10.1051/0004-6361:20034002
I. Pagano 1 - J. L. Linsky2 - J. Valenti3 - D. K. Duncan4
1 - INAF, Catania Astrophysical Observatory, via Santa
Sofia 78, 95125 Catania, Italy
2 -
JILA, University of Colorado and NIST, Boulder, CO 80309-0440,
USA
3 -
Space Telescope Science Institute, 3700 San Martin Dr. Baltimore,
MD
21218, USA
4 -
Department of Astrophysical and Planetary Sciences, University of
Colorado, Boulder, CO 80309-0389, USA
Received 25 June 2003 / Accepted 9 October 2003
Abstract
We describe and analyze HST/STIS observations of the G2 V star
Centauri A (
Cen A, HD 128620), a star similar to the
Sun. The
high resolution echelle spectra obtained with the E140H and E230H gratings
cover the complete spectral range 1133-3150 Å with a resolution of 2.6 km s-1,
an absolute flux calibration accurate to
%, and an absolute wavelength
accuracy of 0.6-1.3 km s-1. We present here a study of the E140H spectrum
covering the 1140-1670 Å spectral range, which includes 671 emission lines
representing 37 different ions and the molecules CO and H2. For
Cen A and
the quiet and active Sun, we intercompare the redshifts, nonthermal line
widths, and parameters of two Gaussian representations of transition region
lines (e.g., Si IV, C IV), infer the electron density from the
O IV intersystem lines, and compare their differential emission measure
distributions.
One purpose of this study is to compare the
Cen A and solar UV spectra
to
determine how the atmosphere and heating processes in
Cen A differ
from the
Sun as a result of the small differences in gravity, age, and chemical
composition of the two stars. A second purpose is to provide an excellent high
resolution UV spectrum of a solar-like star that can serve as a proxy for the
Sun observed as a point source when comparing other stars to the Sun.
Key words: stars: individual: Cen A - stars: chromospheres
Our knowledge and understanding of phenomena related to magnetic activity in late-type stars is based largely on the analysis of observations of the Sun obtained with high spatial, spectral and temporal resolution. In particular, the different heating rates and emission measure distributions of stellar chromospheres and transition regions can be understood by comparing stellar UV spectra with corresponding solar spectra. However, as strange as this may at first appear, we lack a true "reference spectrum'' for the Sun observed as a star for such comparisons. In fact, the existing solar UV spectra provided by instruments on the Solar Maximum Mission (SMM) and the Solar and Heliospheric Observatory (SOHO) typically have moderate to high spectral resolution, but do not represent a full disk average, have uncertain wavelength and absolute flux calibrations, and consist of a stitching together of many small parts of the UV spectrum obtained at different times.
Table 1:
Ultraviolet spectral atlases of the Sun and Cen A.
One way to obtain a close approximation to a high resolution spectrum of the
whole Sun observed as a point source with excellent S/N, absolute flux
calibration, and wavelength accuracy is to observe a bright star with very
similar properties to the Sun. We have done this with the Space Telescope
Imaging Spectrograph (STIS) instrument on HST (Woodgate et al. 1998), obtaining a very
high S/N and high resolution (
)
spectrum of the star
Cen A, a nearby (d=1.34 pc) twin of the Sun with
the same
spectral type (G2 V). Although there are some small differences in effective
temperature and metal abundances between
Cen A and the Sun (see
below), this
STIS spectrum of
Cen A can be considered the best available
"reference
spectrum'' for the Sun viewed as a star, because it is a full disk average, has
excellent wavelength and flux calibration (Bohlin et al. 2001), and covers the
entire 1130-3100 Å UV range with high S/N and within a short period of
time.
Table 2:
Abundances of Cen A in log units.
Cen AB (G2 V + K1 V) is the binary system located closest to the Earth
(d=1.34 pc). It shows an eccentric orbit (e=0.519) with a period of almost 80 years (Pourbaix et al. 2002). Actually
Cen is a triple star system. The third member of the system,
Cen C or Proxima Cen, is a M5.5 Ve flare star (V=11.05) about 12 000 AU distant from
Cen and only d=1.29 pc from the Sun (Perryman et al. 1997). Thanks to the high apparent brightness (V=-0.01 and V=1.33 for the A and B component, respectively) and large parallax of the
Cen
stars, their surface abundances, other stellar properties, and astrometric
parameters are among the best known of any star except the Sun. Guenther & Demarque (2000),
Morel et al. (2000), and Pourbaix et al. (2002) have reviewed recent determinations of
the physical characteristics of
Cen AB. According to
Morel et al. (2000) and references therein,
Cen A has nearly the same
surface temperature of the Sun (
K), slightly lower
gravity than the Sun (
,
i.e. 0.76
), and a mass of
- which is probably an upper limit, given different estimates
reported in the literature starting from 1.08
(Guenther & Demarque 2000). The same authors give
a metal overabundance of
0.2 dex with respect
to the Sun, but similar Li and Be abundances to the Sun. In
Table 2 we list the
Cen A abundances
used in this paper, which were compiled from Feltzing & Gonzalez (2001) and
Morel et al. (2000).
The age of
Cen A is controversial:
Morel et al. (2000)
derive an age in the range 2.7-4.1 Gyr depending on the adopted convection
model, while Guenther & Demarque (2000) estimate an age in the range 6.8-7.6 Gyr. One could
argue that
Cen A is younger than the Sun on the basis that it is
formed of
metal enriched material, but the larger radius and lower gravity compared to
the Sun argue that the star is more evolved and somewhat older than the Sun,
even considering its somewhat larger mass.
A closer analog to the Sun is 18 Sco (V=5.50), but this star is too faint to get
high S/N high resolution UV spectra with STIS.
Cen has been extensively studied in the ultraviolet by IUE. Jordan et al. (1987)
used IUE data to create simple one-dimensional models of the atmospheric structure of the two stars. Hallam et al. (1991) have studied the rotational modulation of the most prominent lines in IUE spectra of
Cen A and found a rotation period of about 29 d. This is consistent with the
Boesgaard & Hagen (1974) estimate that the
Cen A rotation period is 10% larger than the solar one,
but is larger than the
22 d rotation period derived from the
km s-1rotational velocity measured by Saar & Osten (1997), assuming a radius of
1.2
and an orbital inclination of
79
.
Ayres et al. (1995) have studied the time variability of the most prominent UV lines of
Cen A
and B during about 11 years of observations. While a clear evidence of a solar-like activity cycle
was found for
Cen B, UV line fluxes from
Cen A do not give any clear
indication for an activity cycle.
In this paper we report on the Cen A spectrum recorded with the E140 grating
by HST/STIS between 1140-1670 Å, while the analysis of the E230H spectrum (1620-3150 Å) will be published in a forthcoming paper. Information on data acquisition and
reduction are
provided in Sect. 2, the spectral line identification and the
analysis of interesting lines are presented in Sect. 3. A
detailed comparison of our STIS
Cen A spectrum, with the SOHO/SUMER (Curdt et al. 2001) and the SMM/UVSP (Shine & Frank 2000) spectra of the
Sun is given in Sect. 5. Then, we derive the
Cen A transition
region electronic densities (Sect. 6), and emission measure distribution
(Sect. 7). In Sect. 8 we call the reader's attention on some
absorption features present in high exicitation lines, and give our conclusions in Sect. 9.
![]() |
Figure 1:
The E140H spectrum of ![]() |
Open with DEXTER |
The E140H spectrum of Cen A was acquired on 1999 Feb. 12 with 3 exposures of 4695 s each, centered at 1234, 1416, and 1598 Å, respectively.
The E140H mode ensures an average dispersion of
Å per pixel, which corresponds to a resolving power of 2.6 km s-1. The E140H grating
is used with the FUV-MAMA detector, which we operated in TIME-TAG mode. We
used
the
arcsec aperture.
The data were reduced using the STIS Science Team's IDL-based software, CALSTIS (version 6.6). CALSTIS performs a variety of functions including flat fielding, assignment of statistical errors, compensation for the Doppler shifts induced by the spacecraft's motion in orbit, conversion of counts to count rates, dark-rate image subtraction, and the removal of data from bad/hot pixels. Wavelength calibration was carried out assuming the post launch echelle dispersion coefficients and a dispersion coefficient correction for the Monthly MSM offsets released to the STIS Science Team on 1999 September (Lindler 1999b). The on-board Pt lamp spectra taken in association with the science observations were used to measure zero point adjustments. For the echelle observations, CALSTIS computes a wavelength offset for each spectral order. The adopted offset is the median of these offsets. As a check for the success of the algorithm used, we have verified that the measured offsets are all within one pixel of the median offset. As a further check on the accuracy of the wavelength scale, we measured the centroids of emission lines recorded in adjacent orders, and found that the results agree to within less than 1 pixel. The nominal absolute wavelength accuracy is 0.5-1 pixel (i.e. 0.6-1.3 km s-1) (Leitherer et al. 2001).
CALSTIS outputs a file containing wavelength, flux, and error vectors, which is
used in all subsequent processing. To remove the effects of scattered light
that are important near the Ly- line, we used the IDL ECHELLE_SCAT
routine (Lindler 1999a) in the STIS Science Team's software package.
This routine uses the first estimate of the spectrum and a scattering model of
the spectrograph to determine the intensity of the scattered light and to
estimate what the spectrum plus scattered light image should look like.
Comparison of this calculated spectrum with the observations yields
differences that indicate the errors in the first estimate of the spectrum.
This spectrum is then corrected and the process is iterated until acceptable
agreement is obtained between the prediction and the observed image.
After correction for scattered light, the spectrum was then analyzed using
software packages written in IDL. We used routines of the ICUR fitting
code, adapted to handle our
STIS data, which perform multi-Gaussian fits to the line profiles using
Bevington's (1969) CURFIT algorithm. To correct for instrumental broadening, we
convolved each proposed fit to an emission line profile with the instrumental
line spread function (LSF), which was assumed to be a Gaussian with the
nominal width ranging from
1.2 pixel at 1200 Å to
1 pixel at 1700 Å (Leitherer et al. 2001), as is appropriate for lines which are much
broader than the width of the LSF.
Table 3:
Lines detected in the STIS E140H spectrum of Cen A.
In Fig. 1 we show the E140H spectrum of Cen A. We have
measured a
total of 662 emission features of which 77 are due to blends of two or more
lines, 71 are due to unidentified transitions, and 514 are identified as due
to single emission lines. Taking into account the 157 lines identified in
blended features, we find a total of 671 emission lines in this spectrum. In
Table 3 we list all the ions that have been identified. Most of these
lines are due to Si I, Fe II, C I, which together
contribute 441 lines, but S I and Ni II are each represented by more than 30 lines.
Table 4 lists the line identifications, laboratory and measured wavelengths, radial velocity shifts corrected for the stellar radial velocity of -23.45 km s-1, computed using the orbital parameters and ephemeris given by Pourbaix et al. (2002), line full-widths at half-maxima (FWHM), and line fluxes. The laboratory wavelengths listed in Table 4 are from Sandlin et al. (1986), unless otherwise noted in the table. We used Gaussian fits to the line profiles to measure wavelengths, FWHM, and fluxes for single or blended emission lines which do not show central reversals. For the lines which have interstellar absorption components or central reversals, i.e. the most intense optically thick chromospheric lines of C I, O I, Si II, and C II, we instead integrated the flux contained in a suitable wavelength interval and tabulated the FWHM of the observed profile. In Table 4 these lines are indicated with "CR'' in the Notes column.
The strongest transition region lines show broad wings, and therefore do not have a Gaussian profile. For these lines, we list in Table 4 the line centroid, the FWHM of the observed profile and the flux integrated in a suitable wavelength interval. The analysis of these lines is reported in Sect. 3.3.
Absorption features due to the interstellar medium have been measured in a number of lines originating in transitions from the ground level. Such lines are indicated with "ISM'' in Table 4 (column Notes). They will be discussed in a separate paper, together with the derived properties of the interstellar medium along this line of sight.
Several intersystem lines are present in the spectrum of the Cen A, including the O IV UV 0.01 intercombination multiplet 2s22p2P0J-2s2p2 4PJ, that are
diagnostics of electron density (cf. Del Zanna et al. 2002; Brage et al. 1996, and
references therein), the N IV line at 1486 Å, and the
O III line at 1666 Å. We have used these lines to measure
densities in the
Cen A chromosphere and transition
region as discussed in Sect. 6.
The chromospheric Ly
emission line is altered greatly by the
superimposed narrow, weak deuterium (D I) interstellar absorption and by
very broad, saturated hydrogen (H I) interstellar, heliospheric, and
astrospheric absorption, and by geocoronal emission.
The Ly
line flux given in Table 4 was estimated by fitting a Gaussian to the
wings of the line profile, disregarding the central part of the line, which is strongly affected
by ISM absorption and geocoronal emission, and including a second Gaussian to account for
the Deuterium absorption. This is a very rough estimate of the Ly
flux. We refer to
the Linsky& Wood (1996) and Wood et al. (2001) papers for reliable estimates of the intrinsic Ly
in
Cen A.
The two O V lines that we have measured in the Cen A STIS
spectrum have radial velocities differing by about 1.8 km s-1, with the 1218 Å line less red-shifted than the line at 1371 Å. On the Sun, the 1371 Å line
has a Doppler shift of
5 km s-1 greater than the 1218 Å line,
but
Brekke (1993a) concluded that such a difference between the two lines can be
explained only by an error in the adopted laboratory wavelength of the
O V 1218 Å line, which is an intersystem line and thus difficult to
measure in the laboratory. However, if this were the case, adoption of the
wavelength 1218.325 Å suggested by Brekke (1993a) as the laboratory
wavelength, leads to a significant difference (2.8 km s-1) in the
opposite sense.
We suggest that the main reason of the slight wavelength disagreement, even on
the Sun, can be attributed to the difficulty in measuring the wavelength of the
O V 1218 Å line (see Fig. 2b)
in the sloping wing of the Ly
line. The
O V line at 1371 Å (see Fig. 2c)
shows a double peak with an apparent central
reversal. We know of no explanation for this effect as the line is unlikely
to be optically thick and thus self-reversed, and interstellar absorption is
also unlikely.
A blow-up of the region with a complex feature located near 1241.8 Å is
shown in Fig. 3. The feature is noisy, but its double-peak structure
is preserved even after smoothing with a boxcar average
of width as large as 13 pixels. We have therefore fitted the profile
with two Gaussians, and identified the two lines as
the S I 1241.9 Å and Fe XII 1242 Å lines.
Since the Fe XII line is formed at temperature
,
its
predicted
thermal width is
33 km s-1. We have frozen the line width of the
Fe XII line to its thermal width, and derived a flux of
erg s-1 cm-2. An a-posteriori check for the
accuracy of our measured flux
is given by the excellent agreement between the
emission measure derived by using this line at
and the emission
measure derived
at temperatures
and 6.3 from Chandra spectra (Raassen et al. 2003) (cf.
Sect. 7).
![]() |
Figure 2: Blow-up of the regions containing the S V 1199.134 Å line (panel a)), the O V 1218 & 1371 Å lines (panels b) and c)), and the Fe IV quadruplet between 1601.5 and 1606.5 Å (panel d)). Of this quadruplet we could measure only the 1602 Å line. Light-ink labels in panel d) indicate the positions of the missed Fe IV lines. The symbol * in panels a) and c) marks the absorption components due to the interstellar medium. In all of the panels the wavelength scale has been shifted according to the stellar radial velocity. |
Open with DEXTER |
![]() |
Figure 3: A double-peak structure identified as the S I at 1241.9 Å and the Fe XII 1242 Å lines. |
Open with DEXTER |
The weak emission feature observed in solar spectra
at 1356.88 Å was tentatively
attributed to the S III line at 1357.0 Å by Feldman et al. (1975). Its
laboratory wavelength makes this line slightly blue-shifted, in contrast with
the expectation (cf. Sect. 4), therefore we can argue that either
the identification is wrong or the laboratory wavelength given by
Feldman et al. (1975) is inaccurate.
We have measured the Fe IV line at 1602 Å that belongs to a multiplet
of four lines. A careful inspection of the spectrum shows slight flux
increments at the wavelengths corresponding to the 1603.181, and 1603.730 Å Fe IV lines, which, however, are below our detection limit as shown in
Fig. 2 (panel d), but we do not
find any appreciable emission feature corresponding to the fourth line of this
multiplet at 1606.333 Å. While all of the Fe IV lines have been
identified in the Kelly (1982) line database, no Fe IV lines have been
identified in the solar spectrum analyzed by Sandlin et al. (1986). We have inspected
the solar SMM/UVSP spectrum (Shine & Frank 2000) to look for Fe IV lines,
but even the strongest line measured in the Cen A STIS spectrum
at
1656 Å is missing in the solar spectrum, as shown in the left-bottom
panel
of Fig. 4.
![]() |
Figure 4:
Plots of interesting portions of the ![]() |
Open with DEXTER |
As shown by Wood et al. (1997), the strongest transition region emission lines of
Cen A have profiles
with broad wings. We find that broad wings are present in
the Si III
1206 Å, N V
1238 Å,
Si IV
1393 & 1402 Å, and C IV
1548 & 1502 Å line profiles. For these lines we used one narrow Gaussian component (NC) to fit the line core and one broad Gaussian component (BC) to fit the
broad wings (see Fig. 5). This bi-modal structure of the transition
region lines is typically observed for several RS CVn-type stars (i.e.,
Capella and HR 1099), main sequence type stars (i.e., AU Mic, Procyon,
Cen A, and
Cen B), and the giants 31 Com,
Cet,
Dra,
Gem, and AB Dor (Linsky et al. 1995; Pagano et al. 2000; Linsky & Wood 1994).
Wood et al. (1997) showed that the narrow components can be produced by turbulent
wave dissipation or Alfvén wave heating mechanisms, while the broad
components, that resemble the explosive events on the Sun, are diagnostics for
microflare heating. Analysis of SUMER data led Peter (2001) to propose
an alternative explanation for the broad Gaussians, which he calls the "tail
component'', seen in lines formed at temperatures between 50 000 and 300 000 K
in the chromospheric network. He argues that the tail component originates in
coronal funnels that magnetically connect the lower transition region with the
corona, and the broadening is by passing magneto-acoustic waves.
![]() |
Figure 5: The Si III, N V, Si IV, and C IV line profiles. The narrow and broad dashed lines indicate the narrow and broad Gaussian components, respectively, required to best fit the broad wings of these transition region lines. The vertical solid and dashed lines indicate the centroids of the narrow and broad Gaussians, respectively. |
Open with DEXTER |
Table 5 lists the parameters resulting from our multi-Gaussian
fits. Both
the narrow and broad components are redshifted with respect to the stellar
chromosphere, whose rest velocity is determined by the mean velocity of 80
selected Si I lines as discussed in Sect. 4. The narrow components
show larger redshifts as is seen in solar data (Peter 2001). This effect
was also noticed by Wood et al. (1997), who analyzed the Si IV 1393 Å line in
GHRS/HST spectra of
Cen A.
The broad and narrow Gaussian components have comparable intensity as
the flux-weighted mean ratio between the flux in the broad component and the
total flux is
.
This ratio is typical for the most active
stars studied by Wood et al. (1997), and it appears to be independent of the
activity level of the star.
Table 5:
Parameters derived from the multi-Gaussian fits to the transition
region emission lines of
Cen A which show broad wings. Flux is in units of 10-15.
The flux-weighted average of the FWHMs are
,
and
km s-1 for the narrow and broad components, respectively. By comparison,
explosive events on the Sun produce transition region lines as broad as
km s-1 (Dere et al. 1989).
Lines of the same ion generally form at nearly the same temperature in a collisional ionization equilibrium plasma. Therefore, for most ions we use all the measured lines to derive their mean Doppler shifts and nonthermal widths. The results, listed in Table 6, have been derived according to the following procedure. First, we have computed the standard deviation of the heliospheric velocities measured for all the unblended lines of each ion. Then, we have selected the lines whose velocity is different from the mean by less than 1 standard deviation in order to remove from the analysis lines that might be altered by unknown blends or have inaccurate wavelengths. With these selected lines we then computed the mean heliospheric velocity and standard deviation of the mean, as well as the mean FWHM. For some ions this procedure was not applied - e.g. in the case of ions for which less than 3 lines have been measured - as notated in the last column of Table 6.
Table 6:
Doppler shift and nonthermal velocities of chromospheric and
transition region lines measured in the STIS E140H spectrum of Cen A.
The most probable nonthermal
speeds ()
listed in Table 6 were
computed from the measured FWHM (in km s-1) by:
The Doppler shift and nonthermal
widths of the chromospheric and transition
region lines measured in the E140H spectrum of Cen A are plotted in
Figs. 6 and 7. In both figures polynomial fits to
the SUMER measurements of radial velocities and nonthermal widths in a
solar active region and in the quiet Sun, derived by Teriaca et al. (1999), are
represented with dotted and dashed curves, respectively. According to
Teriaca et al. (1999), in solar active regions the lines formed at temperatures
between
K and
K are red-shifted, with
a maximum red-shift about 15 km s-1 at
105 K (C IV). At
higher
temperatures the velocities decrease becoming blue-shifted (about -10 km s-1 at
K). However, in the quiet Sun the Doppler shift reaches a maximum
at a slightly higher temperature,
K (O IV,
N V), and then decreases to a blue-shift of about -2 km s-1 at
K (Ne VIII).
![]() |
Figure 6:
Doppler shifts of chromospheric and transition region lines of ![]() |
Open with DEXTER |
We performed a 2nd order polynomial fit to the Cen A line Doppler
shifts with
respect to the temperature of line formation, giving each data point a weight
equal to the square root of the inverse of its standard deviation. We find that
in the sampled temperature range (
K to
K) the redshift increases monotonically, but the data are not adequate to
infer the temperature of the turnover. Even though the lines of C II,
Si II, and Fe II are believed to be optically thick, these
lines do
not show Doppler shifts different from the optically thin lines formed at the
same formation temperature. The N V lines, especially the 1238 Å line, have broad wings, and their centroids show a smaller redshift than is
expected by the general distribution. In fact, the broad components are
generally marginally blueshifted relative to the narrow components, as
discussed in Sect. 3.3.
![]() |
Figure 7: Nonthermal velocities of chromospheric and transition region lines as function of the temperature of line formation. The solid line represents a third order polynomial fit to the data. Lines which are possibly affected by opacity (Fe II, Si II, and C II) are not included in the polynomial fit ( square data points). The dotted and dashed lines represent the nonthermal line widths in a solar active region and in the quiet Sun, respectively, as derived by Teriaca et al. (1999). |
Open with DEXTER |
For both active and quiet regions on the Sun, the distributions of nonthermal
line width versus line formation temperature, derived by including only those
lines that are not affected by opacity effects, show a peak at
K. To map the increase in turbulent velocity with line formation
temperature (and hence approximate height in the atmosphere), we fit the
Cen A data with a third order polynomial using the widths of the
optically
thin lines. In the sampled temperature interval, the turbulent velocity
distribution for
Cen A resembles the solar data, although slightly
larger
nonthermal line widths are measured for line formation temperatures greater
than
K. We note that the O III mean line width is probably
underestimated because the 1660 Å line may be blended, which can
alter its intensity and consequently the determination of its width.
Solar UV spectra with comparable spectral resolution to our HST/STIS spectrum
of Cen A are those observed by the UltraViolet Spectrometer and
Polarimeter
(UVSP) instrument which flew on the the Solar Maximum Mission (SMM) and by the
SUMER (Solar Ultraviolet Measurements of Emitted Radiation) spectrograph now
operating on SOHO (Solar and Heliospheric Observatory).
The UVSP/SMM spectrum was obtained during the minimum of solar cycle (1984) in the
range 1290-1772 Å with a
slit, oriented north-south
near solar disk center,
with spectral resolution of the order of 100 000. The atlas (2nd order) was
prepared
by Richard Shine and Zoe Frank of the Lockheed-Martin Space and Astrophysics
Lab.,
and was retrieved from the site
ftp://umbra.nascom.nasa.gov/pub/uv_atlases/. According to Shine (private communication), this spectrum was calibrated using
Rottman's quiet Sun data from rocket flights, which had accurate flux scales but
had low spectral resolution. For intercomparison purposes, the wavelength scale
of the solar spectrum was shifted by performing a cross correlation between
this spectrum and the
Cen A spectrum
in many selected wavelength
intervals.
The SUMER/SOHO spectrum is the FUV part of the spectrum that has been derived
from observations obtained in the range 670-1609 Å by Curdt et al. (2001).
These data were acquired with a dispersion of 41.2 mÅ/pixel (1st order) at 1500 Å, for an effective resolution of
8.2 km s-1. The
wavelengths are typically accurate to 10 mÅ, i.e. 2 to 5 km s-1. The
data represent the average radiance (mW sr-1m-2Å-1) for the
quiet Sun at disk center (April 20, 1997), a coronal hole (Oct. 12, 1996),
and a solar spot (Mar. 18, 1999). Hence, the quiet Sun and coronal hole SUMER spectra
were acquired during phases of minimum of the solar cycle.
To be comparable with the Cen A spectrum, we have computed the solar
irradiance at the
Cen distance from the quiet Sun radiance at disk
center by multiplying by
(cf.
Wilhelm et al. 1998) for the quiet Sun, sunspot, and coronal hole spectra. This
conversion does not take into account center-to-limb variations in the lines
and continuum.
The comparison between the STIS Cen A and the UVSP solar spectra can
be made in the 1192-1688 Å spectral range. In Fig. 4 we
plot interesting regions of the UVSP spectrum and the
Cen A spectrum
degraded to a resolution of 0.010 Å/pixel in order to be comparable to the
UVSP spectrum. The wavelength scale of the
Cen A spectrum was shifted
to compensate for the radial velocity of the system (-23.45 km s-1).
Since the UVSP data refer to the "mean intensity over the disk'', it is
possible to perform a radiometric comparison with the Cen A spectrum.
For the emission lines whose
integrated flux in the STIS spectrum exceed
erg s-1 cm-2,
we list in Table 7 line surface fluxes
and full widths at half maximum (FWHM) for both
Cen A
and the Sun. We find that the line widths for
the two stars are very similar for
most of the
chromospheric lines, whereas the transition region lines are typically broader
for
Cen A compared to the Sun. We show in Fig. 8 the FWHM
ratios versus the temperatures of line formation. A linear fit to these data
suggests that the two quantities are correlated with a correlation coefficient
of 0.83. Typically the
Cen A line surface
fluxes are slightly larger than
those of the
Sun (see Fig. 9), with a mean flux ratio
(Sun/
Cen A) of
,
but the Si II 1526 and 1533 Å,
the
He II 1640 Å and the Al II 1671 Å lines are stronger in the Sun than in
Cen A. The interstellar medium absorption in the Si II 1526 Å and Al II lines of
Cen A can partially explain the
high flux ratios. However, this is not the case for
the Si II 1533 Å
line. The common factor for the three
lines is the presence of central reversals.
In the Sun these lines form in the
chromosphere, where temperature increases with height. Line source
functions, however, first increase and then decrease with height
over the line formation region, due to non-LTE effects. A central
reversal occurs for an optically thick line when
the line core forms above the region where
the source function peaks. The "horns'' of the observed profile
form roughly where the source function peaks (Mauas et al. 1989). On the Sun we see
that the depth of the central reversal is a function of
position on the solar
surface. For example, the C I lines of the multiplet at 1560 and 1657 Å show profiles which have central reversals quite deep at the limb, and
less
pronounced both above the limb
and towards disk center (Roussel-Dupré 1983).
The UVSP spectrum was acquired near disk center, while the
Cen A spectrum is a
full disk average. As a consequence, we expect less pronounced central
reversals for the strongest chromospheric lines
in the solar UVSP spectrum than in
the
Cen A spectrum, as is observed.
Table 7:
Surface fluxes (in units of 103 erg cm-2 s-1), and FWHM (km s-1) of a selection of lines present in the spectra
of both Cen A and the Sun. Included lines
have fluxes in the
Cen A spectrum greater than
erg s-1 cm-2.
![]() |
Figure 8:
Ratios of the solar to ![]() |
Open with DEXTER |
We believe that the He II line, which is optically thin and not
self-reversed, is really weaker on Cen A
than on the Sun. Since
the He II line is extremely sensitive to the coronal activity, a flux
ratio of 1.8 suggests that the Sun is more active than
Cen A. This
conclusion is strengthened by the absence of a limb contribution in the solar
data since
the He II line is limb-brightened for the Sun and likely also for
Cen A.
![]() |
Figure 9:
Ratios of the solar to ![]() |
Open with DEXTER |
It should be mentioned here that some of the UV line flux
differences between the Cen A and the Sun may be due to
observing at different phases of the two stellar activity cycles. In fact transition
region lines on the Sun can vary up to factors 2-5 over a magnetic cycle. While
both UVSP and SUMER spectra were acquired during phases of solar minimum, we do not
know the phase of
Cen A activity cycle at which our observations
have been obtained. In fact, the long-time extended IUE data base of
Cen A
does not provide any hints of an activity cycle period (Ayres et al. 1995). The roughly
13 individual measurements of the C IV multiplet flux
obtained by IUE over a 13 year time period (see Fig. 11a in their paper)
have a mean value of
but a range from 1.5 to 3.6 in these units. The STIS flux for the C IV multiplet is 2.80
in the same units. This STIS flux is about 12% larger than the mean
IUE flux. Since the IUE data are low resolution (about 6 Å), the broad
line wings could be difficult to measure compared to the continuum and
nearby weak lines. Thus the C IV fluxes measured from IUE spectra are
likely somewhat low, and the C IV flux observed by STIS is probably very
close to the mean value observed by IUE. We therefore believe that
Cen A had average transition region fluxes when it was observed
by STIS.
In Fig. 10 we plot interesting regions of the Cen A spectrum,
and
the SUMER spectra of a coronal hole, a sunspot, and the quiet Sun,
respectively. The wavelength scale of the
Cen A spectrum was shifted
to remove the radial velocity of the star (-23.45 km s-1).
The SUMER spectrum
has the best photon statistics, therefore faint lines can be more easily seen
in the solar spectrum than in the
Cen STIS spectrum. On the other
hand, the STIS spectrum has better resolution, which can be useful in resolving
line blends and in studying line reversals due to optical thickness effects
better than with the solar spectrum. In Table 8 we summarize how
many lines we found in common between the two spectra. Many of the lines
present in the solar spectrum but not in
Cen A are located at
wavelengths below about 1500 Å, whereas many of the lines detected only in
the STIS spectrum are at wavelengths above 1500 Å. These differences probably
result from the different resolutions of the two data sets and the increasing
S/N of the SUMER data to shorter wavelengths. The emission lines in the
Cen A spectrum are much stronger than in the quiet Sun spectrum. This
is most likely because we are comparing the solar spectrum, which is an average
disk-center quiet Sun spectrum, with the
Cen A full disk irradiance
spectrum, which includes emission from the limb. For most of the lines, the
full disk irradiance is nearly a factor of two larger than the irradiance
derived from disk-center radiance data (cf., Wilhelm et al. 1999), but there is
no difference for the continuum. Because of this effect, the radiometric
comparison between the SUMER and STIS spectra is uncertain.
![]() |
Figure 10:
Plots of interesting portions of the ![]() ![]() |
Open with DEXTER |
Table 8: Measured features in the 1170-1610 Å common spectral range.
The ratios of lines emitted by the same ion can be sensitive to electron
density when the upper levels of the two transitions are depopulated in
different ways. However, misleading results can be obtained when using line
ratios that have a very small sensitivity at the inferred electron densities or
when temperature effects are not properly taken into account. For this reason
we have computed transition region densities for Cen A using both the
line
ratios method and the so-called L-functions method, as described by
Landi & Landini (1997). The main advantage of this method is that it gives an
overall view of all lines and clearly shows which lines (and not line ratios)
are more suitable in a particular density regime and when lines are at the
limit of their density sensitive regime. According to these authors, the contribution function for each line of a selected ion,
), can be
expressed as a product of two functions, one depending on electron density and
electron temperature, and the other on temperature alone:
The analysis of the Cen A electron densities have been carried out
with the
help of the CHIANTI database VERSION 4.0 (Dere et al. 1997; Young et al. 2003),
assuming the ionization equilibria described by Mazzotta et al. (1998) and the
Cen A photospheric abundances listed in Table 2.
The O IV intercombination multiplet near 1400 Å can be used as a
density diagnostic in the range
cm-3(Brage et al. 1996). The 5 lines of the multiplet are all measured in the
Cen A STIS spectrum. The 1399.780 and 1407.382 Å lines originate
from a
common upper level; their ratio is
,
consistent with the ratio
of their A-values (1.08), which is the expected value for a branching ratio in
the optical thin case (Jordan 1967). The O IV 1404.806 Å line is
blended with S IV 1404.808 Å, and possibly with another unknown line
(Del Zanna et al. 2002). The percentage of the blending attributable to the O IV line has been derived and discussed by many authors (cf. Del Zanna et al. 2002; Dufton et al. 1982; Brage et al. 1996). The analysis of the five O IV line ratios (1401/1399, 1401/1407, 1401/1404, 1404/1407, and 1404/1399) indicate that
is in the range 9.8-10.2, assuming that the effective
temperature of O IV formation is
and that the
O IV 1404 Å line accounts for about 70-80% of the blend
with the S IV 1404 Å line. This result does not change strongly with the assumed
temperature at which the O IV lines are formed. Different
estimates of
the O IV and S IV relative contributions to the 1404 Å blend
have been found to result in densities inconsistent with those obtained by
ratios not involving the 1404 Å line. This is true if one assumes that
either the O IV accounts for
92% of the blend, which is obtained
from the theoretical line ratio of the S IV 1404.808 Å and 1406.016 Å lines assuming the atomic calculations by Dufton et al. (1982), or that
the O IV line accounts for
50% of the blend, as derived by Del Zanna et al. (2002) in their analysis of a solar flare and of a GHRS spectrum of
Capella. On the other hand, in their analysis of FUSE and STIS data for the
dM1e star AU Mic Del Zanna et al. (2002) conclude that the O IV contribution to
the
blend is
80%, which is similar to our analysis.
In Fig. 11 (top panel)
the L-functions of the O IV lines, computed at
,
are plotted versus
.
Apart from
the 1397 Å line, which is very weak and results in large errors,
the other lines meet at
,
assuming that the O IV 1404 Å line accounts for 70% of the blend.
The L-functions of the S IV 1404, 1406, and 1407 Å lines, plotted in
Fig. 11 (middle panel),
show that in order to have
consistency among the lines, the S IV 1404 Å line must not exceed 10%
of the flux in the blend. Therefore, an unidentified line contributes 10-20% of the total flux to the 1404 Å blend. Moreover, the
L-functions of the
S IV lines clearly show that no reliable density measurements can be
derived from these lines.
![]() |
Figure 11:
The L-function curves (as defined by Landi & Landini 1998) plotted for
the O IV, S IV, and O V 1218 and 1371 Å lines
observed in the STIS spectrum of ![]() |
Open with DEXTER |
The observed O V 1218/1371 Å line ratio is
.
Using
the Chianti code, and
assuming the ionization equilibrium as described in
Mazzotta et al. (1998) and the
Cen A photospheric abundances listed in
Table 2, we find that:
Either an overestimation of the 1218 Å O V line flux or an
underestimation of the 1317 Å line flux can produce this higher than
expected flux ratio. The 1218 Å line would have to be less than half the
measured value, or the 1371 Å line would
have to be more than double the measured value, to be consistent with a density
at
.
As discussed in
Sect. 3.2 and shown in Fig. 2, the O V line at 1371 Å appears in the STIS spectrum with a double peak, indicating an apparent central reversal or
overlying absorption. We know of no explanation for this effect, but it could
be the cause of a slight flux understimation for this line. On the other hand,
in Sect. 7 we show that the density sensitive O V 1218 Å line does not match the differential emission measure distribution
determined from
the allowed lines at any density, unless we assume that its actual flux is from 20 to 30% less than measured. We think that such an error can be ascribed to
the difficulty in measuring the O V line in the sloping wing of the
Ly
line, as already discussed in Sect. 3.2, where we also
show a slight discrepancy in the radial velocities measured for the
O V 1218 and 1371 Å lines.
The physical properties of the transition region plasma can be determined by
means of
an emission measure distribution analysis (cf. Jordan & Brown 1981; Pagano et al. 2000; Dere & Mason 1981).
The frequency integrated flux Fji of an effectively thin emission line
between levels j and i of an atom, in units erg cm-2 s-1, can be
written as:
![]() |
(4) |
We have derived the
emission measure loci as a function of temperature for each emission line of
interest
using Eq. (6), the
functions computed with the
CHIANTI database VERSION 4.0 (Dere et al. 1997; Young et al. 2003),
the ionization equilibrium as in Mazzotta et al. (1998), and the
Cen A photospheric abundances listed in Table 2.
Following the procedure described in Pagano et al. (2000), we derive the
differential
emission measure distribution shown in Fig. 12. We used the
allowed lines of Si II, S II, C II, S III, Si III,
C III, S IV, O V, and Fe XII observed in the STIS spectrum, which are labelled in the last column of Table 3 with the
letters "emd''. We have also used the allowed lines C II 1036 & 1037 Å, N III 990 Å, S III 1077 Å, Si III 1108 Å,
S IV 1062 & 1063 Å, Ne V 1145 Å, and O VI 1037 Å observed in the FUSE spectrum of
Cen A (Redfield et al. 2002). When more
than
two lines for a given ion have been used, the errorbars in Fig. 12
represent the standard deviation of the emission measures computed for the
different lines. Otherwise, the errorbar is indicative of the uncertainty due
to the line flux measurement. The allowed C IV 1548 & 1502 Å,
Si IV 1393 & 1402 Å, N V 1238 & 1242 Å, and S VI 933 Å lines were not used to derive the differential emission measure
distribution,
because Del Zanna et al. (2002) showed that such lines from the Li and Na isoelectronic
sequences, which were commonly used in previous literature, produce erroneous
results in the determination of emission measures. In Fig. 12
emission measure loci of the C IV, Si IV, N V, and S VI lines are, in fact, anomalous with respect to the emission measure
distribution derived from the other ions.
![]() |
Figure 12:
The differential
emission measure distribution of ![]() ![]() |
Open with DEXTER |
In Fig. 12 we also plot for comparison the differential emission
measure
distribution for the quiet Sun (Landi & Landini 1998) and for solar active
regions
(Dere & Mason 1993).
There is close
agreement between the differential emission measure distributions of Cen A and the quiet
Sun in the
range 5.0-5.6. For temperatures below
,
the
emission measure is larger for
Cen A than for the quiet Sun. At
temperatures
higher than
we have only the
Fe XII 1242 Å line in the STIS spectrum,
therefore it is not possible
to constrain the
real slope of the emission measure distribution. There is, however, a reasonable
good agreement between the emission measure of Fe XII
and the values derived
at temperatures
and 6.3 from Chandra spectra (Raassen et al. 2003).
The spin-forbidden lines can be used to obtain information on the plasma
density by comparing their emission measure loci, computed for different values
of the density, with the emission measure distribution derived by using the
resonance lines. In fact, for collisionally de-exicited spin-forbidden lines
(i.e., when
and
cm-3 s-1 is the
total collision rate out of level j), the ratio
in
Eq. (3) is proportional to
.
Therefore, the emission
measure loci of spin-forbidden lines depend upon both electron density and
temperature. In Fig. 12 we have plotted the emission measures of
the intersystem O III 1666Å, S IV 1406 Å, N IV 1486 Å,
and O V 1218 Å lines, and the mean emission measure of the O IV lines at 1397, 1399, 1401, and 1407 Å. The O III and O IV ions
match the emission measure distribution derived from allowed lines for electron
density
,
and lie above at higher densities. The S IV,
N IV, and O V lines do not strictly match the emission measure
distribution at any density. A possible explanation for this behaviour is that
the fluxes of these lines are slightly overestimated (no more
than 50%). As shown in Fig. 13, the S IV 1406 Å and
N IV 1486 Å lines are very well detected in the STIS
Cen A spectrum,
and a flux overestimation could be caused only by unknown blends.
Alternatively, inaccurate atomic data could be the cause of the
observed discrepancy. For the O V 1218 Å line, a flux
from 20 to 30% less than measured would make
this line consistent with the
emission measure distribution derived from the allowed lines. As seen in
Sect. 6.2, the anomalous ratio
between 1218 and 1371 Å line
fluxes also suggests that the O V 1218 Å line flux was
overestimated,
possibly because of the sloping Ly
wings (cf. Sect. 3.2).
![]() |
Figure 13: The S IV] 1406 Å line (panel a)) and the N IV] 1486 Å line (panel b)) discussed in Sect. 7. |
Open with DEXTER |
The total power radiated per unit (surface) area from the stellar
atmosphere is:
With the same method described above, we have computed the power loss from
Solar
Active regions and from the Quiet Sun adopting their emission measure
distributions as in Fig. 12. For the temperature range
range log
we find a power loss of
erg s-1 cm-2 (
)
from Solar Active Regions, and of
erg s-1 cm-2 (
)
from the Quiet Sun. The
Cen A power loss in the same temperature range is therefore midway between those of the
Quiet Sun and the Solar Active
Regions, but closer to the former than to
the latter.
A number of lines in the Cen A spectrum contain interstellar
absorption
components (look for label "ISM'' in the Notes column of the
Table 3). The analysis of these components will be presented and
analyzed in a subsequent
paper.
We call here attention to the narrow
absorption features in the Si IV 1393 Å, C IV 1548 Å, and N V 1238 Å lines shown in Fig. 5. We think that these
features could be real. In fact, if they were an artifact of the
line spread function, then they would appear in all emission
lines, not just the high excitation lines. Also the possibility that these
features are the results of order
overlapping, in the case the wavelength
scales of two adjacent orders were not consistent, was examined and
disregarded.
In fact, not one of the above emission lines lies in the overlap region
of two adjacent orders. Although the central reversals in N V
and perhaps C IV could be noise, the reversal in the Si IV line appears to be
too deep and includes too many pixels to be just noise. The reality of the
Si IV feature suggests that the other features could be real. If the
reversals are real, then what is their cause? Self-absorption should
affect both lines in a doublet, but the weaker members of the doublet do
not show absorption features. Perhaps we are detecting narrow absorption
by some cool species above the hot regions. Since these possible
absorption features are too narrow to be thermal, we have no explanation
for them, and further observations and analysis are needed to verify
their reality and, if real, search for their cause.
We present our analysis of HST/STIS observations of Cen A and compare
its
spectrum with its near twin, the Sun:
(1) We present a high resolution (
) spectrum
of
Cen A obtained using the E140H mode of STIS that covers the
spectral range 1140-1670 Å with very high signal-to-noise. The spectrum has an absolute
flux calibration accurate to
5%, an absolute wavelength accuracy of 0.6-1.3 km s-1, and and is corrected for scattered light. To our knowlege this
is the best available ultraviolet spectrum of a solar-like star.
(2) As strange as this may at first appear, there is no available ultraviolet
reference spectrum of the Sun as a point source with the characteristics of the
Cen A spectrum that can be used to compare stellar spectra with the
Sun. Many
ultraviolet spectra of the Sun do exist, but they either have lower spectral
resolution, lack wavelength or flux accuracy, or do not include the
center-to-limb variation across the solar disk required to provide an accurate
spectrum of the Sun as a point source. Although
Cen A differs slightly
from
the Sun in effective temperature, gravity, and metal abundance, its spectrum
can serve as a representative solar spectrum for comparison with other stars.
(3) We compare the Cen A spectrum to the solar irradiance (the Sun
viewed as
a point source) derived from UVSP data for the "mean intensity over the disk''
by placing the Sun at the distance to
Cen A and shifting the
Cen A spectrum
by the star's radial velocity. The line widths of the two stars are similar for
chromospheric lines, but the transition region lines of
Cen A are
broader
than those of the Sun by roughly 20%. The line surface fluxes are typically larger on
Cen A, presumably due to
Cen A being somewhat metal rich.
However, the
He II 1640 Å line is stronger in the Sun, indicating that the solar
corona is more active.
(4) We also compare the Cen A spectrum to the solar irradiance derived
from
SUMER spectra of the disk center quiet Sun, assuming constant center-to-limb
radiance and shifting the
Cen A wavelength scale by the radial
velocity of
the star. A total of 671 emission lines are detected in the
Cen A spectrum
from 37 different ions and 2 molecules (CO and H2). In addition to the well
known chromospheric and transition region lines, we also identify lines of Al IV, Si VIII, S V, Ca VII, Fe IV, Fe V,
and Fe XII. A total of 172 emission lines observed in
Cen A are
not
seen in the SUMER spectrum.
(5) Broad wings are present in the strong resonance lines of C IV, N V, Si III, and Si IV, as are seen in solar observations of the chromospheric network. We fit the line profiles with two Gaussians: a narrow component ascribed to Alfvén waves in small magnetic loops, and a broad component ascribed to microflares or magneto-acoustic waves in large coronal funnels. Both components are redshifted with the narrow Gaussians having larger redshifts as is seen on the Sun. At line formation temperatures between 20 000 K and 200 000 K, there is a trend of increasing line redshift, similar to but with a somewhat lower magnitude than the quiet Sun. A similar trend of increasing nonthermal velocities with temperature is nearly identical to that which is observed in solar quiet and active regions.
(6) Using line ratios and L-functions, we infer that the O IV lines
are formed where the electron density is
.
The S IV and O V lines, however, do not provide reliable
values.
Values of
have been obtained for
the Sun and other solar type stars (cf. Cook et al. 1995).
It is hard to make any comparison with these results because they are strongly affected by
the adopted atomic calculation or by the choice of lines with a limited density sensitivity,
as the O IV 1400 Å line (see Del Zanna et al. 2002). Hence estimates of
from
different computations are often not consistent. We can, however, compare the electron
density we derive for
Cen A at
with the electron density derived by Del Zanna et al. (2002) for
Capella (G1 III + G8 III) and AU Mic (dM1e), because we use the same
computation methods. The comparison tell us that at
the electron density
is slightly less in
Cen A than in the more active Capella (
)
and AU Mic (
).
(7) The emission measure distribution of Cen A derived from emission
lines of
ions not in the Li and Na isoelectronic sequences is in close agreement with
that of
the quiet Sun in the temperature range
,
but lies somewhat
above the quiet Sun in the temperature range
.
This
could be explained by the higher metal abundance of
Cen A combined
with a
somewhat less active corona that provides less conductive heating to the upper
transition region. The estimated total radiative power loss from the transition
region (
)
is
erg s-1 cm-2,
corresponding to
.
Acknowledgements
This work is supported by NASA grant S-56500-D to NIST and the University of Colorado. We thanks Dr. Joseph B. Gurman, who kindly provided us information about the UVSP/SMM solar spectrum, and Dr. Richard Shine and Dr. Zoe Frank who made it available.