zeta Ori RGS model fitting (5T apec)*vwindtab*tbabs

~/xspec/2015/zpup_rgs/

model specification has abundances tied among the model components (5 bvapecs and vwindtab); velocities are tied too

temps are frozen and norms are free

Sigma_* is a free parameter in vwindtab
but nH is frozen at the literature value in tbabs (.028)

Free params are (potentially): 

norms 1 - 5
velocity
C, N, O, Fe (He too in principle)
Sigma_*

7jun16

allowing all the parameters to be free (saved log from 27may had Fe frozen at 1 still and inadvertently v = 500 km/s frozen)

OK, v goes to about 700

but when I thaw Fe, though the fit statistic improves significantly, the CNO abundances all go down significantly. 

Let's step back and reset Fe to 1.0 and freeze it for now

First model fit: 

5T_nonsolar.dat

T = .15, .3, .6, 1.2, 2.4 keV

norm 1 - 5

5.34e-3
4.12e-3
1.09e-3
4.01e-4
0.0

CNO
1.14
1.06
1.01

vel = 703
Sigma = 2.41e-3

chi2 = 10116.10 
N = 2865

Next model - solar

5T_solar.dat

T = .15, .3, .6, 1.2, 2.4 keV

norm 1 - 5

5.41e-3
3.93e-3
1.08e-3
3.99e-4
8.20e-6

CNO = 1, 1, 1

vel = 698
Sigma = 2.21e-3

chi2 = 10177.97 
N = 2865

Next, try same two schemes, but with a different set of 5Ts: 0.13, 0.26, 0.52, 1.04, 2.08 

*do solar first from now on*

5T_solar_8jun16a.dat

T = 0.13, 0.26, 0.52, 1.04, 2.08 keV

norm 1 - 5

5.40e-3
5.39e-3
1.61e-3
5.53e-4
0

CNO = 1, 1, 1

vel = 702
Sigma = 3.36e-3

chi2 = 9925.78 
N = 2865

5T_nonsolar_8jun16a.dat

T = 0.13, 0.26, 0.52, 1.04, 2.08 keV

norm 1 - 5

5.82e-3
5.65e-3
1.47e-3
5.69e-4
0.0

CNO
0.925
0.987
0.916

vel = 703
Sigma = 3.36e-3

chi2 = 9854.40 
N = 2865

Now, we'll add a free redshift to mimic mild line asymmetry (expect blue shift)
otherwise the same as the 8jun16a model

5T_nonsolar_10jun16a.dat

T = 0.13, 0.26, 0.52, 1.04, 2.08 keV

norm 1 - 5

5.72e-3
5.15e-3
1.37e-3
5.62e-4
0.0

CNO
0.904
0.956
0.899

redshift = -4.85e-4

vel = 692
Sigma = 2.69e-3

chi2 = 8750.16 
N = 2865

(30jun16):

5T_nonsolar_30jun16a.dat

T = 0.12, 0.192, 0.307, 0.492, 0.786 keV (I'd like to add a 6th one: 1.26, but for now will not)
note that this is 0.12 * 1.6^i

norm 1 - 5

5.22e-3
2.47e-3
4.16e-3
7.57e-5
9.76e-4

CNO
0.861
0.946
0.917

redshift = -4.93e-4

vel = 687
Sigma = 2.38e-3

chi2 = 8486.89 
N = 2865

6T_nonsolar_30jun16b.dat
(mod file with date 30jun16b)

T = 0.12, 0.192, 0.307, 0.492, 0.786, 1.26 keV
same 0.12 * 1.6^i only with i up to 5 rather than 4

norm 1 - 6

5.11e-3
2.53e-3
3.93e-3
1.29e-4
9.21e-4
4.91e-5

CNO
0.853
0.934
0.906

redshift = -4.95e-4

vel = 686
Sigma = 2.17e-3

chi2 = 8482.29 
N = 2865
  
6T_nonsolar_5jul16a.dat

T = 0.11, 0.187, 0.318, 0.540, 0.919, 1.56 keV
same 0.11 * 1.7^i 

norm 1 - 6

4.23e-3
3.89e-3
2.89e-3
5.95e-4
6.68e-4
0.0 

CNO
0.835
0.920
0.896

redshift = -4.96e-4

vel = 687
Sigma = 1.67e-3

chi2 = 8648.49 
N = 2865

Note: chi2 worse; 6th norm to zero
So, try the closer (1.6) spacing but start with T = 0.11

6T_nonsolar_5jul16b.dat

T = 0.11, 0.176, 0.282, 0.451, 0.721, 1.15 keV
same 0.11 * 1.6^i 

norm 1 - 6

7.08e-3
2.01e-3
5.21e-3
2.20e-4
1.01e-3
1.68e-4 

CNO
0.698
0.876
0.870

redshift = -4.97e-4

vel = 687
Sigma = 3.02e-3

chi2 = 8260.15 
N = 2865