Normalized to: Delahaye, F.
[1]
oai:arXiv.org:1804.03365 [pdf] - 1663482
A Quantitative Comparison of Opacities Calculated Using the Distorted-
Wave and $\boldsymbol{R}$-Matrix Methods
Delahaye, F.;
Badnell, N. R.;
Ballance, C. P.;
Palmeri, P.;
Preval, S.;
Quinet, P.;
Ramsbottom, C.;
Smyth, R. T.;
Turkington, M.;
Zeippen, C. J.
Submitted: 2018-04-10
The present debate on the reliability of astrophysical opacities has reached
a new climax with the recent measurements of Fe opacities on the Z-machine at
the Sandia National Laboratory \citep{Bailey2015}. To understand the
differences between theoretical results, on the one hand, and experiments on
the other, as well as the differences among the various theoretical results,
detailed comparisons are needed. Many ingredients are involved in the
calculation of opacities; deconstructing the whole process and comparing the
differences at each step are necessary to quantify their importance and impact
on the final results. We present here such a comparison using the two main
approaches to calculate the required atomic data, the $R$-Matrix and
distorted-wave methods, as well as sets of configurations and coupling schemes
to quantify the effects on the opacities for the $Fe\ XVII$ and $Ni\ XIV$ ions.
[2]
oai:arXiv.org:1511.07260 [pdf] - 1358992
IPOPv2 online service for the generation of opacity tables
Submitted: 2015-11-23
In the framework of the present phase -- IPOPv2 -- of the international
Opacity Project (OP), a new web service has been implemented based on the
latest release of the OP opacities. The user may construct online opacity
tables to be conveniently included in stellar evolution codes in the format
most commonly adopted by stellar physicists, namely the OPAL format. This
facility encourages the use and comparison of both the OPAL and OP data sets in
applications. The present service allows for the calculation of multi-element
mixtures containing the 17 species (H, He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar,
Ca, Cr, Mn, Fe and Ni) considered by the OP, and underpins the latest release
of OP opacities. This new service provides tables of Rosseland mean opacites
using OP atomic data. We provide an alternative to the OPAL opacity services
allowing direct comparison as well as study of the effect of uncertainties in
stellar modeling due to mean opacities.
[3]
oai:arXiv.org:1401.3122 [pdf] - 1550245
Ipopv2: Photoionization of Ni XIV -- a test case
Submitted: 2014-01-14
Several years ago, M. Asplund and coauthors (2004) proposed a revision of the
Solar composition. The use of this new prescription for Solar abundances in
standard stellar models generated a strong disagreement between the predictions
and the observations of Solar observables. Many claimed that the Standard Solar
Model (SSM) was faulty, and more specifically the opacities used in such
models. As a result, activities around the stellar opacities were boosted. New
experiments (J. Bailey at Sandia on Z-Pinch, The OPAC consortium at LULI200) to
measure directly absorbtion coefficients have been realized or are underway.
Several theoretical groups (CEA-OPAS, Los Alamos Nat. Lab., CEA-SCORCG, The
Opacity Project - The Iron Project (IPOPv2)) have started new sets of
calculations using different approaches and codes. While the new results seem
to confirm the good quality of the opacities used in SSM, it remains important
to improve and complement the data currently available. We present recent
results in the case of the photoionization cross sections for Ni XIV (Ni13+ )
from IPOPv2 and possible implications on stellar modelling.
[4]
oai:arXiv.org:1312.3885 [pdf] - 1202240
The chemical composition of the Sun from helioseismic and solar neutrino
data
Submitted: 2013-12-13
We perform a quantitative analysis of the solar composition problem by using
a statistical approach that allows us to combine the information provided by
helioseimic and solar neutrino data in an effective way. We include in our
analysis the helioseismic determinations of the surface helium abundance and of
the depth of the convective envelope, the measurements of the $^7{\rm Be}$ and
$^8{\rm B}$ neutrino fluxes, the sound speed profile inferred from helioseismic
frequencies. We provide all the ingredients to describe how these quantities
depend on the solar surface composition and to evaluate the (correlated)
uncertainties in solar model predictions. We include errors sources that are
not traditionally considered such as those from inversion of helioseismic data.
We, then, apply the proposed approach to infer the chemical composition of the
Sun. We show that the opacity profile of the Sun is well constrained by the
solar observational properties. In the context of a two parameter analysis in
which elements are grouped as volatiles (i.e. C, N, O and Ne) and refractories
(i.e Mg, Si, S, Fe), the optimal composition is found by increasing the the
abundance of volatiles by $\left( 45\pm 4\right)\%$ and that of refractories by
$\left( 19\pm 3\right)\%$ with respect to the values provided by AGSS09. This
corresponds to the abundances $\varepsilon_{\rm O}=8.85\pm 0.01$ and
$\varepsilon_{\rm Fe}=7.52\pm0.01$. As an additional result of our analysis, we
show that the observational data prefer values for the input parameters of the
standard solar models (radiative opacities, gravitational settling rate, the
astrophysical factors $S_{34}$ and $S_{17}$) that differ at the $\sim 1\sigma$
level from those presently adopted.
[5]
oai:arXiv.org:1307.1010 [pdf] - 1172466
Simulations of protostellar collapse using multigroup radiation
hydrodynamics. II. The second collapse
Submitted: 2013-07-03
Star formation begins with the gravitational collapse of a dense core inside
a molecular cloud. As the collapse progresses, the centre of the core begins to
heat up as it becomes optically thick. The temperature and density in the
centre eventually reach high enough values where fusion reactions can ignite;
the protostar is born. This sequence of events entail many physical processes,
of which radiative transfer is of paramount importance. Many simulations of
protostellar collapse make use of a grey treatment of radiative transfer
coupled to the hydrodynamics. However, interstellar gas and dust opacities
present large variations as a function of frequency. In this paper, we
follow-up on a previous paper on the collapse and formation of Larson's first
core using multigroup radiation hydrodynamics (Paper I) by extending the
calculations to the second phase of the collapse and the formation of Larson's
second core. We have made the use of a non-ideal gas equation of state as well
as an extensive set of spectral opacities in a spherically symmetric fully
implicit Godunov code to model all the phases of the collapse of a 0.1, 1 and
10 solar mass cloud cores. We find that, for a same central density, there are
only small differences between the grey and multigroup simulations. The first
core accretion shock remains supercritical while the shock at the second core
border is found to be strongly subcritical with all the accreted energy being
transfered to the core. The size of the first core was found to vary somewhat
in the different simulations (more unstable clouds form smaller first cores)
while the size, mass and temperature of the second cores are independent of
initial cloud mass, size and temperature. Our simulations support the idea of a
standard (universal) initial second core size of 0.003 AU and mass 0.0014 solar
masses.
[6]
oai:arXiv.org:1201.6245 [pdf] - 469945
Comparison of Fe and Ni opacity calculations for a better understanding
of pulsating stellar envelopes
Gilles, D.;
Turck-Chièze, S.;
Loisel, G.;
Piau, L.;
Ducret, J. -E.;
Poirier, M.;
Blenski, T.;
Thais, F.;
Blancard, C.;
Cossé, P.;
Faussurier, G.;
Gilleron, F.;
Pain, J. -C.;
Porcherot, Q.;
Guzik, J. A.;
Kilcrease, D. P.;
Magee, N. H.;
Harris, J.;
Busquet, M.;
Delahaye, F.;
Zeippen, C. J.;
Bastiani-Ceccotti, S.
Submitted: 2012-01-30
Opacity is an important ingredient of the evolution of stars. The calculation
of opacity coefficients is complicated by the fact that the plasma contains
partially ionized heavy ions that contribute to opacity dominated by H and He.
Up to now, the astrophysical community has greatly benefited from the work of
the contributions of Los Alamos [1], Livermore [2] and the Opacity Project (OP)
[3]. However unexplained differences of up to 50% in the radiative forces and
Rosseland mean values for Fe have been noticed for conditions corresponding to
stellar envelopes. Such uncertainty has a real impact on the understanding of
pulsating stellar envelopes, on the excitation of modes, and on the
identification of the mode frequencies. Temperature and density conditions
equivalent to those found in stars can now be produced in laboratory
experiments for various atomic species. Recently the photo-absorption spectra
of nickel and iron plasmas have been measured during the LULI 2010 campaign,
for temperatures between 15 and 40 eV and densities of ~3 mg/cm3. A large
theoretical collaboration, the "OPAC", has been formed to prepare these
experiments. We present here the set of opacity calculations performed by eight
different groups for conditions relevant to the LULI 2010 experiment and to
astrophysical stellar envelope conditions.
[7]
oai:arXiv.org:1201.4693 [pdf] - 467147
On plasma radiative properties in stellar conditions
Submitted: 2012-01-23
The knowledge of stellar evolution is evolving quickly thanks to an increased
number of opportunities to scrutinize the stellar internal plasma properties by
stellar seismology and by 1D and 3D simulations. These new tools help us to
introduce the internal dynamical phenomena in stellar modeling. A proper
inclusion of these processes supposes a real confidence in the microscopic
physics used, partly checked by solar or stellar acoustic modes. In the present
paper we first recall which fundamental physics has been recently verified by
helioseismology. Then we recall that opacity is an important ingredient of the
secular evolution of stars and we point out why it is necessary to measure
absorption coefficients and degrees of ionization in the laboratory for some
well identified astrophysical conditions. We examine two specific experimental
conditions which are accessible to large laser facilities and are suitable to
solve some interesting questions of the stellar community: are the solar
internal radiative interactions properly estimated and what is the proper role
of the opacity in the excitation of the non radial modes in the envelop of the
$\beta$ Cephei and the Be stars ? At the end of the paper we point out the
difficulties of the experimental approach that we need to overcome.
[8]
oai:arXiv.org:1201.4692 [pdf] - 1093173
Iron and Nickel spectral opacity calculations in conditions relevant for
pulsating stellar envelopes and experiments
Gilles, D.;
Turck-Chièze, S.;
Busquet, M.;
Thais, F.;
Loisel, G.;
Piau, L.;
Ducret, J. E.;
Blenski, T.;
Blancard, C.;
Cossé, P.;
Faussurier, G.;
Gilleron, F.;
Pain, J. C.;
Porcherot, Q.;
Guzik, J. A.;
Kilcrease, D. P.;
Magee, N. H.;
Harris, J.;
Bastiani-Ceccotti, S.;
Delahaye, F.;
Zeippen, C. J.
Submitted: 2012-01-23
Seismology of stars is strongly developing. To address this question we have
formed an international collaboration OPAC to perform specific experimental
measurements, compare opacity calculations and improve the opacity calculations
in the stellar codes [1]. We consider the following opacity codes: SCO,
CASSANDRA, STA, OPAS, LEDCOP, OP, SCO-RCG. Their comparison has shown large
differences for Fe and Ni in equivalent conditions of envelopes of type II
supernova precursors, temperatures between 15 and 40 eV and densities of a few
mg/cm3 [2, 3, 4]. LEDCOP, OPAS, SCO-RCG structure codes and STA give similar
results and differ from OP ones for the lower temperatures and for spectral
interval values [3]. In this work we discuss the role of Configuration
Interaction (CI) and the influence of the number of used configurations. We
present and include in the opacity code comparisons new HULLAC-v9 calculations
[5, 6] that include full CI. To illustrate the importance of this effect we
compare different CI approximations (modes) available in HULLAC-v9 [7]. These
results are compared to previous predictions and to experimental data.
Differences with OP results are discussed.
[9]
oai:arXiv.org:1110.1873 [pdf] - 423198
How to Make an Atomic Blog in Your Own Kitchen. Summary of the Workshop:
Uncertainties in Atomic Data and How They Propagate in Chemical Abundances
Luridiana, Valentina;
García-Rojas, Jorge;
Aggarwal, Kanti;
Bautista, Manuel;
Bergemann, Maria;
Delahaye, Franck;
del Zanna, Giulio;
Ferland, Gary;
Lind, Karin;
Manchado, Arturo;
Mendoza, Claudio;
Delgado, Adal Mesa;
Díaz, Manuel Núñez;
Shaw, Richard A.;
Wesson, Roger
Submitted: 2011-10-09
This workshop brought together scientists (including atomic physicists,
theoretical astrophysicists and astronomers) concerned with the completeness
and accuracy of atomic data for astrophysical applications. The topics covered
in the workshop included the evaluation of uncertainties in atomic data, the
propagation of such uncertainties in chemical abundances, and the feedback
between observations and calculations. On a different level, we also discussed
communication issues such as how to ensure that atomic data are correctly
understood and used, and which forum is the best one for a fluid interaction
between all communities involved in the production and use of atomic data. This
paper reports on the discussions held during the workshop and introduces
AstroAtom, a blog created as a platform for timely and open discussions on the
needs and concerns over atomic data, and their effects on astronomical
research.
The complete proceedings will be published on
http://astroatom.wordpress.com/.
[10]
oai:arXiv.org:1101.4955 [pdf] - 1051690
A numerical model for multigroup radiation hydrodynamics
Submitted: 2011-01-25
We present in this paper a multigroup model for radiation hydrodynamics to
account for variations of the gas opacity as a function of frequency. The
entropy closure model (M1) is applied to multigroup radiation transfer in a
radiation hydrodynamics code. In difference from the previous grey model, we
are able to reproduce the crucial effects of frequency-variable gas opacities,
a situation omnipresent in physics and astrophysics. We also account for the
energy exchange between neighbouring groups which is important in flows with
strong velocity divergence. These terms were computed using a finite volume
method in the frequency domain. The radiative transfer aspect of the method was
first tested separately for global consistency (reversion to grey model) and
against a well established kinetic model through Marshak wave tests with
frequency dependent opacities. Very good agreement between the multigroup M1
and kinetic models was observed in all tests. The successful coupling of the
multigroup radiative transfer to the hydrodynamics was then confirmed through a
second series of tests. Finally, the model was linked to a database of
opacities for a Xe gas in order to simulate realistic multigroup radiative
shocks in Xe. The differences with the previous grey models are discussed.
[11]
oai:arXiv.org:1101.1170 [pdf] - 1042920
Radiative properties of stellar plasmas and open challenges
Turck-Chièze, S.;
Loisel, G.;
Gilles, D.;
Piau, L.;
Blancard, C.;
Blenski, T.;
Busquet, M.;
Caillaud, T.;
Cossé, P.;
Delahaye, F.;
Faussurier, G.;
Fariaut, J.;
Gilleron, F.;
Guzik, J. A.;
Harris, J.;
Kilcrease, D. P.;
Magee, N. H.;
Pain, J. C.;
Porcherot, Q.;
Poirier, M.;
Soullier, G.;
Zeippen, C. J.;
Bastiani-Ceccotti, S.;
Reverdin, C.;
Silvert, V.;
Thais, F.;
Villette, B.
Submitted: 2011-01-06
The lifetime of solar-like stars, the envelope structure of more massive
stars, and stellar acoustic frequencies largely depend on the radiative
properties of the stellar plasma. Up to now, these complex quantities have been
estimated only theoretically. The development of the powerful tools of helio-
and astero- seismology has made it possible to gain insights on the interiors
of stars. Consequently, increased emphasis is now placed on knowledge of the
monochromatic opacity coefficients. Here we review how these radiative
properties play a role, and where they are most important. We then concentrate
specifically on the envelopes of $\beta$ Cephei variable stars. We discuss the
dispersion of eight different theoretical estimates of the monochromatic
opacity spectrum and the challenges we need to face to check these calculations
experimentally.
[12]
oai:arXiv.org:1101.1182 [pdf] - 1042921
Theoretical and experimental activities on opacities for a good
interpretation of seismic stellar probes
Turck-Chièze, S;
Loisel, G;
Gilles, D;
Thais, F;
Bastiani, S;
Blancard, C;
Busquet, M;
Caillaud, T;
Cosse, P;
Blenski, T;
Delahaye, F;
Ducret, J E;
Faussurier, G;
Gilleron, F;
Guzik, J;
Harris, J W;
Kilcrease, D P;
Magee, N H;
Piau, L;
Pain, J C;
Poirier, M;
Porcherot, Q;
Reverdin, C;
Silvert, V;
Villette, B;
Zeippen, C
Submitted: 2011-01-06
Opacity calculations are basic ingredients of stellar modelling. They play a
crucial role in the interpretation of acoustic modes detected by SoHO, COROT
and KEPLER. In this review we present our activities on both theoretical and
experimental sides. We show new calculations of opacity spectra and comparisons
between eight groups who produce opacity spectra calculations in the domain
where experiments are scheduled. Real differences are noticed with real
astrophysical consequences when one extends helioseismology to cluster studies
of different compositions. Two cases are considered presently: (1) the solar
radiative zone and (2) the beta Cephei envelops. We describe how our
experiments are performed and new preliminary results on nickel obtained in the
campaign 2010 at LULI 2000 at Polytechnique.
[13]
oai:arXiv.org:1005.0423 [pdf] - 159257
Helioseismic Constraints on the Solar Ne/O Ratio and Heavy Element
Abundances
Submitted: 2010-05-03
We examine the constraints imposed by helioseismic data on the solar heavy
element abundances. In prior work we argued that the measured depth of the
surface convection zone R_CZ and the surface helium abundance Y_surf were good
metallicity indicators which placed separable constraints on light metals
(CNONe) and the heavier species with good relative meteoritic abundances. The
resulting interiors-based abundance scale was higher than some published
studies based on 3D model atmospheres at a highly significant level. In this
paper we explore the usage of the solar sound speed in the radiative interior
as an additional diagnostic, and find that it is sensitive to changes in the
Ne/O ratio even for models constructed to have the same R_CZ and Y_surf. Three
distinct helioseismic tests (opacity in the radiative core, ionization in the
convection zone, and the core mean molecular weight) yield consistent results.
Our preferred O, Ne and Fe abundances are 8.86 +/-0.04, 8.15 +/-0.17 and 7.50
+/-0.05 respectively. They are consistent with the midrange of recently
published 3D atmospheric abundances measurements. The values for O, Ne and Fe
which combine interiors and atmospheric inferences are 8.83 +/-0.04, 8.08
+/-0.09 and 7.49 +/-0.04 respectively.
[14]
oai:arXiv.org:0907.1082 [pdf] - 26069
A Photometric Metallicity Estimate of the Virgo Stellar Overdensity
Submitted: 2009-07-06, last modified: 2009-11-05
We determine photometric metal abundance estimates for individual
main-sequence stars in the Virgo Overdensity (VOD), which covers almost 1000
deg^2 on the sky, based on a calibration of the metallicity sensitivity of
stellar isochrones in the gri filter passbands using field stars with
well-determined spectroscopic metal abundances. Despite the low precision of
the method for individual stars, we derive [Fe/H] = -2.0 +/-0.1 (internal)
+/-0.5 (systematic) for the metal abundance of the VOD from photometric
measurements of 0.7 million stars in the Northern Galactic hemisphere with
heliocentric distances from ~10 kpc to ~20 kpc. The metallicity of the VOD is
indistinguishable, within Delta [Fe/H] < 0.2, from that of field halo stars
covering the same distance range. This initial application suggests that the
SDSS gri passbands can be used to probe the properties of main-sequence stars
beyond ~10 kpc, complementing studies of nearby stars from more
metallicity-sensitive color indices that involve the u passband.
[15]
oai:arXiv.org:0910.4731 [pdf] - 1018081
Metallicity Mapping with gri Photometry: The Virgo Overdensity and the
Halos of the Galaxy
Beers, Timothy C.;
An, Deokkeun;
Johnson, Jennifer A.;
Pinsonneault, Marc H.;
Terndrup, Donald M.;
Delahaye, Franck;
Lee, Young Sun;
Masseron, Thomas;
Carollo, Daniela;
Yanny, Brian
Submitted: 2009-10-25
We describe the methodology required for estimation of photometric estimates
of metallicity based on the SDSS gri passbands, which can be used to probe the
properties of main-sequence stars beyond ~ 10 kpc, complementing studies of
nearby stars from more metallicity-sensitive color indices that involve the u
passband. As a first application of this approach, we determine photometric
metal abundance estimates for individual main-sequence stars in the Virgo
Overdensity, which covers almost 1000 square degrees on the sky, based on a
calibration of the metallicity sensitivity of stellar isochrones in the gri
filter passbands using field stars with well-determined spectroscopic metal
abundances. Despite the low precision of the method for individual stars,
internal errors of in [Fe/H] ~ +/- 0.1 dex can be achieved for bulk stellar
populations. The global metal abundance of the Virgo Overdensity determined in
this way is <[Fe/H]> = -2.0 +/- 0.1 (internal) +/- 0.5 (systematic), from
photometric measurements of 0.7 million stars with heliocentric distances from
~ 10 kpc to ~ 20 kpc. A preliminary metallicity map, based on results for 2.9
million stars in the northern SDSS DR-7 footprint, exhibits a shift to lower
metallicities as one proceeds from the inner- to the outer-halo population,
consistent with recent interpretation of the kinematics of local samples of
stars with spectroscopically available metallicity estimates and full space
motions.
[16]
oai:arXiv.org:0905.3743 [pdf] - 24561
Galactic Globular and Open Clusters in the Sloan Digital Sky Survey. II.
Test of Theoretical Stellar Isochrones
Submitted: 2009-05-22
We perform an extensive test of theoretical stellar models for main-sequence
stars in ugriz, using cluster fiducial sequences obtained in the previous paper
of this series. We generate a set of isochrones using the Yale Rotating
Evolutionary Code (YREC) with updated input physics, and derive magnitudes and
colors in ugriz from MARCS model atmospheres. These models match cluster main
sequences over a wide range of metallicity within the errors of the adopted
cluster parameters. However, we find a large discrepancy of model colors at the
lower main sequence (Teff < ~4500 K) for clusters at and above solar
metallicity. We also reach similar conclusions using the theoretical isochrones
of Girardi et al. and Dotter et al., but our new models are generally in better
agreement with the data. Using our theoretical isochrones, we also derive
main-sequence fitting distances and turn-off ages for five key globular
clusters, and demonstrate the ability to derive these quantities from
photometric data in the Sloan Digital Sky Survey. In particular, we exploit
multiple color indices (g - r, g - i, and g - z) in the parameter estimation,
which allows us to evaluate internal systematic errors. Our distance estimates,
with an error of sigma(m - M) = 0.03-0.11 mag for individual clusters, are
consistent with Hipparcos-based subdwarf fitting distances derived in the
Johnson-Cousins or Stromgren photometric systems.
[17]
oai:arXiv.org:0704.1583 [pdf] - 351
OPserver: interactive online-computations of opacities and radiative
accelerations
Mendoza, C.;
Seaton, M. J.;
Buerger, P.;
Bellorin, A.;
Melendez, M.;
Gonzalez, J.;
Rodriguez, L. S.;
Delahaye, F.;
Palacios, E.;
Pradhan, A. K.;
Zeippen, C. J.
Submitted: 2007-04-12
Codes to compute mean opacities and radiative accelerations for arbitrary
chemical mixtures using the Opacity Project recently revised data have been
restructured in a client--server architecture and transcribed as a subroutine
library. This implementation increases efficiency in stellar modelling where
element stratification due to diffusion processes is depth dependent, and thus
requires repeated fast opacity reestimates. Three user modes are provided to
fit different computing environments, namely a web browser, a local workstation
and a distributed grid.
[18]
oai:arXiv.org:astro-ph/0606077 [pdf] - 82517
The Solar Heavy Element Abundances: II. Constraints from Stellar
Atmospheres
Submitted: 2006-06-04
Estimates of the bulk metal abundance of the Sun derived from the latest
generation of model atmospheres are significantly lower than the earlier
standard values. In Paper I we demonstrated that a low solar metallicity is
inconsistent with helioseismology if the quoted errors in the atmospheres
models (of order 0.05 dex) are correct. In this paper we undertake a critical
analysis of the solar metallicity and its uncertainty from a model atmospheres
perspective, focusing on CNO. We argue that the non-LTE corrections for
abundances derived from atomic features are overestimated in the recent
abundance studies, while systematic errors in the absolute abundances are
underestimated. If we adopt the internal consistency between different
indicators as a measure of goodness of fit, we obtain intermediate abundances
[C/H] = 8.44 +/- 0.06, [N/H] = 7.96 +/- 0.10 and [O/H] = 8.75 +/- 0.08. The
errors are too large to conclude that there is a solar abundance problem, and
permit both the high and low scales. However, the center-to-limb continuum flux
variations predicted in the simulations appear to be inconsistent with solar
data, which would favor the traditional thermal structure and lead to high CNO
abundances of (8.52, 7.96, 8.80) close to the seismic scale. We argue that
further empirical tests of non-LTE corrections and the thermal structure are
required for precise absolute abundances. The implications for beryllium
depletion and possible sources of error in the numerical simulations are
discussed.
[19]
oai:arXiv.org:astro-ph/0603482 [pdf] - 80691
Electron impact excitation of Helium-like ions up to n=4 levels
including radiation damping
Submitted: 2006-03-17
Helium-like ions provide the most important X-ray spectral diagnostics in
high temperature fusion and astrophysical plasmas. We previously presented
computed collision strengths for O~VII including relativistic fine structure,
levels up to the $n=4$ complex and radiation damping of autoionizing
resonances. We have extended this work to other He-like ions (N, Ne, Mg, Al,
Si, S, Ca). The calculations are carried out using the Breit-Pauli R-matrix
(BPRM) method with a 31-level eigenfunction expansion. Collision strengths for
the principal lines important in X-ray plasma diagnostics, w, x, y and z,
corresponding to the 4 transitions to the ground level 1s^2(^1S_0) <-
1s2p(^1P^o_1), 1s2p(^3P^o_2), 1s2p(^3P^o_1), 1s2s(^3S_1), are explicitly shown.
We find the effect of radiation damping to be significant for the forbidden
transitions in heavier He-like ions, which should affect the diagnostic line
ratios. We extrapolated the collision strengths to their values at infinite
energy using the Burgess-Tully extrapolation technique. This is required to
calculate the Maxwellian average collision strengths at high temperature. We
show that the coupling between dipole allowed and inter-combination transitions
affects increasingly the effective collision strengths for the n ^1S_0 - n'
^3P_1 transition as the charge of the ion increases. This clearly affects the
treatment of the extrapolation toward the infinite energy point of the
collision strength. This work is carried out as part of the Iron Project-RmaX
Network.
[20]
oai:arXiv.org:astro-ph/0511779 [pdf] - 78137
The Solar Heavy Element Abundances: I. Constraints from Stellar
Interiors
Submitted: 2005-11-28
The latest solar atmosphere models include non-LTE corrections and 3D
hydrodynamic convection simulations. These models predict a significant
reduction in the solar metal abundance, which leads to a serious conflict
between helioseismic data and the predictions of solar interiors models. We
demonstrate that the helioseismic constraints on the surface convection zone
depth and helium abundance combined with stellar interiors models can be used
to define the goodness of fit for a given chemical composition. After a
detailed examination of the errors in the theoretical models we conclude that
models constructed with the older solar abundances are consistent (<2 \sigma)
with the seismic data. Models constructed with the proposed new low abundance
scale are strongly disfavored, disagreeing at the 15 \sigma level. We then use
the sensitivity of the seismic properties to abundance changes to invert the
problem and infer a seismic solar heavy element abundance mix with two
components: meteoritic abundances, and the light metals CNONe. Seismic
degeneracies between the best solutions for the elements arise for changes in
the relative CNONe abundances and their effects are quantified. We obtain
Fe/H=7.50+/-0.045+/-0.003(CNNe) and O/H=8.86+/-0.041+/-0.025(CNNe) for the
relative CNNe in the GS98 mixture. The inferred solar oxygen abundance disagree
with the abundance inferred from the 3D hydro models. Changes in the Ne
abundance can mimic changes in O for the purposes of scalar constraints.Models
constructed with low oxygen and high neon are inconsistent with the solar sound
speed profile. The implications for the solar abundance scale are discussed.
[21]
oai:arXiv.org:astro-ph/0410744 [pdf] - 68593
Up-dated opacities from the Opacity Project
Submitted: 2004-10-29, last modified: 2004-12-18
Using the code AUTOSTRUCTURE, extensive calculations of inner-shell atomic
data have been made for the chemical elements He, C, N, O, Ne, Na, Mg, Al, Si,
S, Ar, Ca, Cr, Mn, Fe and Ni. The results are used to obtain up-dated opacities
from the Opacity Project, OP. A number of other improvements on earlier work
have also been included.
Rosseland-mean opacities from OP are compared with those from OPAL.
Differences of 5 to 10% occur. OP gives the `Z-bump', at log(T)=5.2, to be
shifted to slightly higher temperatures. The opacities from OP, as functions of
temperature and density, are smoother than those from OPAL. Extensive tests
show that the numerical accuracy of the OP opacities following integration over
frequency mesh and interpolation on temperature-density mesh is better than 1%.
Prior to a number of recent investigations which have indicated a need for a
downward revision in the solar abundances of oxygen and other elements, there
was good agreement between properties of the sun deduced from helioseismology
and from stellar evolution models calculated using OPAL opacities. The
revisions in abundances destroy that agreement. In a recent paper Bahcall et al
argue that the agreement would be restored if opacities for the regions of the
sun between 0.7 and 0.4 solar radii were larger than those given by OPAL by
about 10%. In the region concerned, the present results from OP do not differ
from those of OPAL by more than 2.5%.
[22]
oai:arXiv.org:astro-ph/0411799 [pdf] - 69398
Comparison of radiative accelerations obtained with atomic data from OP
and OPAL
Submitted: 2004-11-30
Here, we compare radiative accelerations (g_rad) derived from the new Opacity
Project (OP) data with those computed from OPAL and some previous data from OP.
For the case where we have full data from OPAL, the differences in the
Rosseland mean opacities between OPAL and the new OP data are within 12% and
less than 30% between new OP and previous OP data (OP1 at CDS). The radiative
accelerations g_rad differ at up to the 17% level when compared to OPAL and up
to the 38% level when compared to OP1. The comparison with OP1 on a larger
(rho-T) space gives a difference of up to 40% for g_rad(C). And it increases
for heavier elements. The differences increase for heavier elements reaching
60% for Si and 65% for S and Fe. We also constructed four representative
stellar models in order to compare the new OP accelerations with prior
published results that used OPAL data. The Rosseland means overall agree better
than 10% for all of our cases. For the accelerations, the comparisons with
published values yield larger differences in general. The published OPAL
accelerations for carbon are even larger relative to OP compare to what our
direct comparisons would indicate. Potential reasons for this puzzling behavior
are discussed. In light of the significant differences in the inferred
acceleration rates, theoretical errors should be taken into account when
comparing models with observations. The implications for stellar evolution are
briefly discussed. The sensitivity of g_rad to the atomic physics may provide a
useful test of different opacity sources.
[23]
oai:arXiv.org:astro-ph/0406472 [pdf] - 65629
Resolution and accuracy of resonances in R-matrix cross sections
Submitted: 2004-06-21
We investigate the effect of resonances in photoionization and recombination
cross sections computed using the R-matrix method. Photoionization and
recombination rates derived from high-resolution cross sections for oxygen ions
are compared with earlier works with less resolution and accuracy, such as in
the widely used Opacity Project data. We find significant differences in
photoionization rates for O II metastable states, averaged over Planck
functions corresponding to ionizing radiation fields, with respect to the
intrinsic accuracy of the calculations and improved resolution. Furthermore,
for highly charged ions other physical effects are also important.
Recombination rate coefficients, averaged over a Maxwellian distribution, are
extremely sensitive to the position and resolution of near-threshold
resonances, and radiation damping, in (e + O VII) --> O VI + hnu. Surprisingly
however, the effect on the monochromatic and the mean Rosseland and Planck
bound-free opacities is relatively small, but may be potentially significant.
[24]
oai:arXiv.org:astro-ph/0302238 [pdf] - 54865
X-ray absorption via K-alpha resonance complexes in oxygen ions
Submitted: 2003-02-12
The K-alpha resonance complexes in oxygen ions O I - O VI are theoretically
computed and resonance oscillator strengths and wavelengths are presented. The
highly resolved photoionization cross sections, with relativistic fine
structure, are computed in the coupled channel approximation using the
Breit-Pauli R-matrix method. A number of strong K-alpha resonances are found to
be appreciable, with resonance oscillator strengths f_r > 0.1. The K-alpha
resonance wavelengths of O I-O VI lie in a relatively narrow wavelength range
22 - 23.5 A, and the X-ray opacity in this region should therefore be
significantly affected by K --> L transitions in oxygen. The results should be
useful in the interpretation of soft X-ray spectra observed from Chandra and
XMM-Newton.
[25]
oai:arXiv.org:astro-ph/0204116 [pdf] - 48661
X-Ray Resonance Opacity of Oxygen and Iron in AGN MCG--6-30-15
Submitted: 2002-04-08, last modified: 2002-04-09
Theoretical differential oscillator strengths related to monochromatic
opacity due to K- and L-shell absorption from oxygen and iron ions are directly
compared with the Chandra and Xmm-Newton spectra of Seyfert~1 galaxy
MCG--6-30-15. We compute the highly resolved continuum with resonances due to O
I - VI and Fe XVI. It is found that the KLn (n --> infty) resonance series
limits in O VI, and the lowest LMM resonance in Fe XVI, lie at the prominent
17.5 A break in the observed spectra. We also calculate and identify, for the
first time, the observed gap in spectral flux at 22-23 A in the Chandra spectra
due to K --> L resonant absorption features from all O-ions O I - VI,
potentially leading to abundance determination of several or all ionization
stages; similar signature gaps may be predicted due to other elements. The
precise atomic parameters are computed in the relativistic close coupling
approximation using the Breit-Pauli R-matrix method. The new X-ray opacities
might possibly distinguish between models of a dusty warm absorber, and/or
gravitational redshift and broadening due to a massive black hole in
MCG--6-30-15.
[26]
oai:arXiv.org:astro-ph/0201535 [pdf] - 47479
Electron impact excitation of helium-like oxygen up to n = 4 levels
including radiation damping
Submitted: 2002-01-31
The primary X-ray diagnostic lines in He-like ions are mainly excited by
electron impact from the ground level to the n = 2 levels, but at high
temperatures n > 2 levels are also excited. In order to describe the atomic
processes more completely collision strengths are computed for OVII including
for the first time all of the following: (i) relativistic fine structure, (ii)
levels up to the n = 4, and (iii) radiation damping of autoionizing resonances.
The calculations are carried out using the Breit-Pauli R-matrix (BPRM) method
with a 31-level eigenfunction expansion. Resonance structures in collision
strengths are delineated in detail up to the n = 4 thresholds. For highly
charged He-like ions radiation damping of autoionizing resonances is known to
be significant. We investigate this effect in detail and find that while
resonances are discernibly damped radiatively as the series limit n --> infty
is approached from below, the overall effect on effective cross sections and
rate coefficients is found to be very small. Collision strengths for the
principal lines important in X-ray plasma diagnostics, w,x,y and z,
corresponding to the 4 transitions to the ground level 1s^2 (^1S_0) <-- 1s2p
(^1P^o_1), 1s2p (^3P^o_2), 1s2p (^3P^o_1), 1s2s (^3S_1), are explicitly shown.
It is found that the effective collision strength of the forbidden z-line is up
to a factor of 4 higher at T < 10^6 K than previous values. This is likely to
be of considerable importance in the diagnostics of photoionized astrophysical
plasmas. Significant differences are also found with previous works for several
other transitions. This work is carried out as part of the Iron Project-RmaX
Network.
[27]
oai:arXiv.org:astro-ph/0002230 [pdf] - 34573
Atomic data from the Iron Project.XLIII. Transition probabilities for Fe
V
Submitted: 2000-02-10
An extensive set of dipole-allowed, intercombination, and forbidden
transition probabilities for Fe V is presented. The Breit-Pauli R-matrix (BPRM)
method is used to calculate 1.46 x 10^6 oscillator strengths for the allowed
and intercombination E1 transitions among 3,865 fine-structure levels dominated
by configuration complexes with n <= 10 and l <= 9. These data are complemented
by an atomic structure configuration interaction (CI) calculation using the
SUPERSTRUCTURE program for 362 relativistic quadrupole (E2) and magnetic dipole
(M1) transitions among 65 low-lying levels dominated by the 3d^4 and 3d^ 4s
configurations. Procedures have been developed for the identification of the
large number of fine-structure levels and transitions obtained through the BPRM
calculations. The target ion Fe VI is represented by an eigenfunction expansion
of 19 fine-structure levels of 3d^3 and a set of correlation configurations. Fe
V bound levels are obtained with angular and spin symmetries SL\pi and J\pi of
the (e + Fe VI) system such that 2S+1 = 5,3,1, L <= 10, J <= 8 of even and odd
parities. The completeness of the calculated dataset is verified in terms of
all possible bound levels belonging to relevant LS terms and transitions in
correspondence with the LS terms. The fine-structure averaged relativistic
values are compared with previous Opacity Project LS coupling data and other
works. The 362 forbidden transition probabilities considerably extend the
available data for the E2 and M1 transtions, and are in good agreement with
those computed by Garstang for the 3d^4 transitions.