2 aug 11 Ne X Ly alpha at 12.134 There's some blending with an iron line on the red side of the line, so try fitting on 11.98 to 12.18 ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 1.87112 +/- 0.277221 3 1 windprof u0 0.684375 +/- 9.78171E-02 4 1 windprof h 0.0 frozen 5 1 windprof tau0star 0.0 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 0 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.07210E-04 +/- 1.60268E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.49909E-03 +/- 5.27175E-04 ________________________________________________________________________ Chi-Squared = 49.27 using 37 PHA bins. Reduced chi-squared = 1.493 for 33 degrees of freedom Null hypothesis probability = 3.406711e-02 *Note formally adequate fit and so XSPEC will allow "error" calculations XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 1.63797 2.29524 (-0.233156,0.424121) XSPEC12>error 1. 3 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0.99 3 0.613383 0 (-0.0709921,-0.684375) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000299137 0.000315417 (-8.03737e-06,8.24241e-06) XSPEC12>error 1. 19 Parameter Confidence Range (1) 19 0.0012022 0.00179632 (-0.000297488,0.000296632) XSPEC12>iplot PLT> wdata 12.134_nonporous.dat anisotropic porosity free hinf: Zero porosity ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 1.82478 +/- 1.92372 3 1 windprof u0 0.675423 +/- 0.293251 4 1 windprof h 6.84283E-10 +/- 0.911501 5 1 windprof tau0star 0.0 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.07058E-04 +/- 1.77392E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.50739E-03 +/- 6.38332E-04 ________________________________________________________________________ Chi-Squared = 49.31 using 37 PHA bins. Reduced chi-squared = 1.541 for 32 degrees of freedom Null hypothesis probability = 2.597460e-02 XSPEC12>error 1. 4 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 4 0 0.132268 (-1.90269e-11,0.132268) next, anisotropic porosity with fixed porosity lengths: aniso, hinf = 0.5 ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 3.75107 +/- 1.57711 3 1 windprof u0 0.743440 +/- 0.257291 4 1 windprof h 0.500000 frozen 5 1 windprof tau0star 0.0 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.16540E-04 +/- 1.62384E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.17641E-03 +/- 5.95619E-04 ________________________________________________________________________ Chi-Squared = 55.86 using 37 PHA bins. Reduced chi-squared = 1.693 for 33 degrees of freedom Null hypothesis probability = 7.722148e-03 XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 3.14342 4.4571 (-0.607757,0.705923) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.659833 0.84408 (-0.0837463,0.100501) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000308473 0.000324684 (-8.05721e-06,8.15339e-06) XSPEC12>error 1. 19 Parameter Confidence Range (1) 19 0.000863853 0.001485 (-0.000312663,0.000308479) XSPEC12>iplot PLT> wdata 12.134_aniso_h0.5.dat aniso, hinf = 1 ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 5.71220 +/- 2.64109 3 1 windprof u0 0.707600 +/- 0.182731 4 1 windprof h 1.00000 frozen 5 1 windprof tau0star 0.0 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.21921E-04 +/- 1.69854E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 9.96628E-04 +/- 6.38320E-04 ________________________________________________________________________ Chi-Squared = 64.79 using 37 PHA bins. Reduced chi-squared = 1.963 for 33 degrees of freedom Null hypothesis probability = 7.798381e-04 XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 4.70745 6.93991 (-1.00477,1.2277) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.640283 0.784756 (-0.0673163,0.0771566) XSPEC12>error 1. 17 Parameter Confidence Range (1) ***Warning: Identical values of the parameter give different values of the statistic. Please check your result for the low end of the confidence range. 17 0.00031344 0.000330181 (-8.48072e-06,8.26006e-06) XSPEC12>steppar 17 3.08e-4 3.35e-4 20 Chi-Squared Delta norm Chi-Squared 17 67.512 2.7245 0 0.000308 66.979 2.1918 1 0.00030935 66.518 1.7305 2 0.0003107 66.119 1.3316 3 0.00031205 65.787 0.9994 4 0.0003134 65.479 0.69187 5 0.00031475 65.232 0.44426 6 0.0003161 65.043 0.25548 7 0.00031745 64.918 0.13031 8 0.0003188 64.821 0.033556 9 0.00032015 64.78 -0.0073866 10 0.0003215 64.811 0.023646 11 0.00032285 64.869 0.081521 12 0.0003242 64.98 0.19298 13 0.00032555 65.162 0.37422 14 0.0003269 65.37 0.58288 15 0.00032825 65.635 0.84815 16 0.0003296 65.966 1.1783 17 0.00033095 66.324 1.5368 18 0.0003323 66.74 1.9531 19 0.00033365 67.209 2.4216 20 0.000335 XSPEC12>error 1. 19 Parameter Confidence Range (1) 19 0.00067822 0.00132136 (-0.000330625,0.000312516) XSPEC12>iplot PLT> wdata 12.134_aniso_h1.dat aniso, hinf = 5 ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 63.9887 +/- 20.7665 3 1 windprof u0 0.746901 +/- 9.71545E-02 4 1 windprof h 5.00000 frozen 5 1 windprof tau0star 0.0 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 1 frozen 9 1 windprof anisotro 1 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.42448E-04 +/- 1.79824E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.30215E-04 +/- 7.20274E-04 ________________________________________________________________________ Chi-Squared = 101.31 using 37 PHA bins. Reduced chi-squared = 3.0699 for 33 degrees of freedom Null hypothesis probability = 7.267806e-09 XSPEC12>iplot PLT> wdata 12.134_aniso_h5.dat Now, isotropic porosity hinf free (zero is preferred): ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 1.85276 +/- 2.64557 3 1 windprof u0 0.683048 +/- 0.134091 4 1 windprof h 9.72570E-14 +/- 9.70878 5 1 windprof tau0star 0.0 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.07038E-04 +/- 1.64883E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.50685E-03 +/- 6.07640E-04 ________________________________________________________________________ Chi-Squared = 49.28 using 37 PHA bins. Reduced chi-squared = 1.540 for 32 degrees of freedom Null hypothesis probability = 2.612004e-02 XSPEC12>error 1. 4 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0 4 0 2.44348 (-6.59837e-15,2.44348) porosity length is relatively unconstrained iso, hinf = 0.5 ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 2.17415 +/- 0.362881 3 1 windprof u0 0.673853 +/- 8.04448E-02 4 1 windprof h 0.500000 frozen 5 1 windprof tau0star 0.0 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.07207E-04 +/- 1.59808E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.48845E-03 +/- 5.27506E-04 ________________________________________________________________________ Chi-Squared = 49.37 using 37 PHA bins. Reduced chi-squared = 1.496 for 33 degrees of freedom Null hypothesis probability = 3.336359e-02 XSPEC12>iplot PLT> wdata 12.134_iso_h0.5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 1.87738 2.60698 (-0.296743,0.43286) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.611098 0.772819 (-0.0627542,0.0989661) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000299393 0.000315625 (-8.122e-06,8.11023e-06) XSPEC12>error 1. 19 Parameter Confidence Range (1) 19 0.00118384 0.00177628 (-0.000299385,0.00029305) iso, hinf = 1 ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 2.56286 +/- 0.490936 3 1 windprof u0 0.668222 +/- 7.22546E-02 4 1 windprof h 1.00000 frozen 5 1 windprof tau0star 0.0 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.07903E-04 +/- 1.59186E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.45330E-03 +/- 5.26491E-04 ________________________________________________________________________ Chi-Squared = 49.51 using 37 PHA bins. Reduced chi-squared = 1.500 for 33 degrees of freedom Null hypothesis probability = 3.243285e-02 XSPEC12>iplot PLT> wdata 12.134_iso_h1.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 2.176 3.12868 (-0.386864,0.565814) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.610248 0.746661 (-0.0579625,0.078451) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000299734 0.000315899 (-8.06278e-06,8.10306e-06) XSPEC12>error 1. 19 Parameter Confidence Range (1) 19 0.00116219 0.00175383 (-0.000297692,0.000293944) iso, hinf = 5 ======================================================================== Model windprof<1> + powerlaw<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 windprof q 0.0 frozen 2 1 windprof taustar 14.0913 +/- 7.58494 3 1 windprof u0 0.651741 +/- 6.55242E-02 4 1 windprof h 5.00000 frozen 5 1 windprof tau0star 0.0 frozen 6 1 windprof beta 1.00000 frozen 7 1 windprof betaSob 0.0 frozen 8 1 windprof numerica 0 frozen 9 1 windprof anisotro 0 frozen 10 1 windprof rosselan 1 frozen 11 1 windprof expansio 0 frozen 12 1 windprof thick 0 frozen 13 1 windprof waveleng A 12.1340 frozen 14 1 windprof shift mA 0.0 frozen 15 1 windprof velocity km/s 2250.00 frozen 16 1 windprof verbose 0 frozen 17 1 windprof norm 3.15026E-04 +/- 1.55939E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 1.13128E-03 +/- 5.73272E-04 ________________________________________________________________________ Chi-Squared = 53.75 using 37 PHA bins. Reduced chi-squared = 1.629 for 33 degrees of freedom Null hypothesis probability = 1.269853e-02 XSPEC12>iplot PLT> wdata 12.134_iso_h5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 9.79761 20.7905 (-4.23216,6.76072) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.607436 0.70205 (-0.0441168,0.0504974) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000307028 0.000322909 (-7.94938e-06,7.93083e-06) XSPEC12>error 1. 19 Parameter Confidence Range (1) 19 0.00081882 0.00143602 (-0.000313668,0.000303535) confidence limit in 2D param space XSPEC12>log 12.134_iso_hinf_vs_taustar_35by35.log Logging to file:12.134_iso_hinf_vs_taustar_35by35.log XSPEC12>steppar 4 0. 5. 35 2 0. 11. 35