... | ... | @@ -31,6 +31,7 @@ The routine that computes the free-free gaunt factors has been extended to inclu |
|
|
|
|
|
The helium ionization balance at low photon and particle densities, and at high particle densities, has always been exact, and this continues to be the case. There was a problem in the helium ion for high radiation densities, in versions 87 and before. The code used three pseudo levels to represent the levels between 7 and 1000, for H, He, and He+. This seemed to work well for the atoms for the cases of high densities, but testing has shown that it did not represent the physics of the high radiation density limit well. The problem is that the pseudo-levels had very large statistical weights, they represented line energies in the far infrared, and had A’s appropriate for lower levels. As a result they had very large induced rates when the photon occupation numbers were large, and this affected populations of lower levels. As a result the atom became too ionized - as much as a factor of two for He+. The following test illustrates this problem:
|
|
|
|
|
|
|
|
|
title helium ionization in high photon density limit
|
|
|
|
|
|
print departure coef
|
... | ... | @@ -54,7 +55,7 @@ In versions 88 and later no pseudo levels are used for any hydrogen or helium at |
|
|
|
|
|
## Heavy elements
|
|
|
|
|
|
The atomic data base, the organization of aspects of the code dealing with storing heavy element ionization, and all the associated routines, have been totally rewritten. The lightest 30 elements are now included. Photoionization data are from Verner et al. (1996), recombination data partially from Verner and Ferland (1996), and roughly 104 lines of the heavy elements have been added (Verner, Verner, and Ferland 1996).
|
|
|
The atomic database, the organization of aspects of the code dealing with storing heavy element ionization, and all the associated routines, have been totally rewritten. The lightest 30 elements are now included. Photoionization data are from Verner et al. (1996), recombination data partially from Verner and Ferland (1996), and roughly 104 lines of the heavy elements have been added (Verner, Verner, and Ferland 1996).
|
|
|
|
|
|
The number of resonance lines has increased by more than an order of magnitude. All resonance lines listed by Verner, Verner, and Ferland (1996) are included. As a result of these many additional lines the cooling function tends to be larger and smoother.
|
|
|
|
... | ... | @@ -65,21 +66,21 @@ The default solar mixture has been changed to Grevesse and Noel (1993). The big |
|
|
|
|
|
## Free-free, line heating and cooling
|
|
|
|
|
|
These are counted in a different but equivalent manner. Now the difference between cooling and heating is used, since this is more numerically stable at high radiation densities. This difference has no physical affect on the predictions, but the printed contributors to the total heating and cooling do appear different.
|
|
|
These are counted in a different but equivalent manner. Now the difference between cooling and heating is used, since this is more numerically stable at high radiation densities. This difference has no physical effect on the predictions, but the printed contributors to the total heating and cooling do appear different.
|
|
|
|
|
|
|
|
|
## Excited state photoionization cross sections
|
|
|
## Excited-state photoionization cross-sections
|
|
|
|
|
|
OP data are now used. For the excited state of Mg+ this is nearly ten times smaller than old screened hydrogenic values. This affects the intensity of Mg II l2798 in some BLR calculations.
|
|
|
|
|
|
|
|
|
### The O+ photoionization cross section
|
|
|
|
|
|
Figure 17 This figure shows the O+ photoionization cross sections now used, compared to the Reilman and Manson values, and used in version 84 and before.
|
|
|
Figure 17 This figure shows the O+ photoionization cross-sections now used, compared to the Reilman and Manson values, and used in version 84 and before.
|
|
|
|
|
|
The Reilman and Manson photoionization cross sections, used before version 87, show a jump in the photoionization cross section at the 2s - 2p edge, and low values above that threshold extending up to the valence electron threshold. Opacity Project cross sections are used in the current version of the code, and these do not show the 2s edge (the OP calculations find that the 2s and 2p electrons are highly correlated). The cross section remains large up to the valence threshold. The difference approaches a factor of two, and this affects high ionization parameter clouds since O+ is the dominant opacity for some energies.
|
|
|
The Reilman and Manson photoionization cross-sections, used before version 87, show a jump in the photoionization cross section at the 2s - 2p edge, and low values above that threshold extending up to the valence electron threshold. Opacity Project cross sections are used in the current version of the code, and these do not show the 2s edge (the OP calculations find that the 2s and 2p electrons are highly correlated). The cross-section remains large up to the valence threshold. The different approaches a factor of two, and this affects high ionization parameter clouds since O+ is the dominant opacity for some energies.
|
|
|
|
|
|
Verner, et al. (1996) comment on all other cases where the photoionization cross sections have changed. There are generally atoms and first ions where Opacity Project data are now available.
|
|
|
Verner, et al. (1996) comment on all other cases where the photoionization cross-sections have changed. There are generally atoms and first ions where Opacity Project data are now available.
|
|
|
|
|
|
-------
|
|
|
|
... | ... | |