9May12 FITTING POROUS MODELS TO Fe XVII 15.014 both RGS1 and RGS2 XSPEC12>ignore 1-2:0.-14.80 980 channels (1-980) ignored in spectrum # 1 1029 channels (1-1029) ignored in spectrum # 2 XSPEC12>ignore 1-2:15.07-** 2386 channels (1007-3392) ignored in spectrum # 1 1982 channels (1056-3037) ignored in spectrum # 2 - non-porous fit (as before) ======================================================================== 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.76508 +/- 9.11001E-02 3 1 windprof u0 0.635302 +/- 2.25144E-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 15.0140 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 6.39349E-04 +/- 1.07805E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 3.13333E-03 +/- 2.33553E-04 ________________________________________________________________________ Chi-Squared = 129.54 using 52 PHA bins. Reduced chi-squared = 2.6988 for 48 degrees of freedom Null hypothesis probability = 2.017168e-09 NOTE: since the fit is poor, we can't calculate confidence limits on model parameters XSPEC12>iplot PLT> wdata fexvii_15.014_rgs12.dat - ANISOTROPIC POROSITY - hinf free ======================================================================== 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.76601 +/- 0.560402 3 1 windprof u0 0.635287 +/- 7.10472E-02 4 1 windprof h 4.11262E-08 +/- 0.306618 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 15.0140 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 6.39384E-04 +/- 1.08346E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 3.13245E-03 +/- 2.49384E-04 ________________________________________________________________________ Chi-Squared = 129.55 using 52 PHA bins. Reduced chi-squared = 2.7564 for 47 degrees of freedom Null hypothesis probability = 1.191920e-09 *same as non-porous XSPEC12>iplot PLT> wdata fexvii_15.014_rgs12_hfree.dat - h_inf = 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.28902 +/- 0.365358 3 1 windprof u0 0.694203 +/- 5.10818E-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 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 15.0140 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 6.47026E-04 +/- 1.07134E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.90620E-03 +/- 2.44120E-04 ________________________________________________________________________ Chi-Squared = 180.40 using 52 PHA bins. Reduced chi-squared = 3.7582 for 48 degrees of freedom Null hypothesis probability = 3.238082e-17 XSPEC12>iplot PLT> wdata fexvii_15.014_rgs12_h0.5.dat - h_inf = 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 4.66468 +/- 0.614634 3 1 windprof u0 0.668726 +/- 4.17636E-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 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 15.0140 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 6.51500E-04 +/- 1.08464E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.79367E-03 +/- 2.52311E-04 ________________________________________________________________________ Chi-Squared = 249.60 using 52 PHA bins. Reduced chi-squared = 5.2001 for 48 degrees of freedom Null hypothesis probability = 4.866367e-29 XSPEC12>iplot PLT> wdata fexvii_15.014_rgs12_h1.dat - h_inf = 2 ======================================================================== 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 8.23614 +/- 1.24225 3 1 windprof u0 0.637105 +/- 3.27472E-02 4 1 windprof h 2.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 15.0140 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 6.56305E-04 +/- 1.10197E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.66643E-03 +/- 2.60609E-04 ________________________________________________________________________ Chi-Squared = 370.11 using 52 PHA bins. Reduced chi-squared = 7.7106 for 48 degrees of freedom Null hypothesis probability = 2.657454e-51 XSPEC12>iplot PLT> wdata fexvii_15.014_rgs12_h2.dat - h_inf = 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 35.9550 +/- 4.47372 3 1 windprof u0 0.647950 +/- 2.55014E-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 15.0140 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 6.65923E-04 +/- 1.12506E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.37274E-03 +/- 2.72765E-04 ________________________________________________________________________ Chi-Squared = 577.23 using 52 PHA bins. Reduced chi-squared = 12.026 for 48 degrees of freedom Null hypothesis probability = 7.359872e-92 XSPEC12>iplot PLT> wdata fexvii_15.014_rgs12_h5.dat (Whoops; change all the file names to contain 'aniso') Can't make a confidence limit contour plot because even the best-fit is too poor - ISOTROPIC POROSITY - h_inf free ======================================================================== 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.74641 +/- 0.764722 3 1 windprof u0 0.635113 +/- 2.25243E-02 4 1 windprof h 9.67032E-13 +/- 3.22854 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 15.0140 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 6.39354E-04 +/- 1.08069E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 3.13549E-03 +/- 2.44842E-04 ________________________________________________________________________ Chi-Squared = 129.61 using 52 PHA bins. Reduced chi-squared = 2.7577 for 47 degrees of freedom Null hypothesis probability = 1.167824e-09 XSPEC12>iplot PLT> wdata fexvii_15.014_rgs12_iso_hfree.dat - h_inf = 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.01868 +/- 0.115890 3 1 windprof u0 0.630975 +/- 2.03525E-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 15.0140 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 6.39454E-04 +/- 1.07723E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 3.11449E-03 +/- 2.33752E-04 ________________________________________________________________________ Chi-Squared = 131.44 using 52 PHA bins. Reduced chi-squared = 2.7382 for 48 degrees of freedom Null hypothesis probability = 1.084596e-09 XSPEC12>iplot PLT> wdata fexvii_15.014_rgs12_iso_h0.5.dat - h_inf = 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.33447 +/- 0.150490 3 1 windprof u0 0.628612 +/- 1.91225E-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 15.0140 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 6.39535E-04 +/- 1.07636E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 3.09522E-03 +/- 2.34090E-04 ________________________________________________________________________ Chi-Squared = 133.62 using 52 PHA bins. Reduced chi-squared = 2.7838 for 48 degrees of freedom Null hypothesis probability = 5.265192e-10 XSPEC12>iplot PLT> wdata fexvii_15.014_rgs12_iso_h1.dat - h_inf = 2 ======================================================================== 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.23128 +/- 0.268633 3 1 windprof u0 0.626717 +/- 1.79565E-02 4 1 windprof h 2.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 15.0140 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 6.40029E-04 +/- 1.07426E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 3.04725E-03 +/- 2.35093E-04 ________________________________________________________________________ Chi-Squared = 139.31 using 52 PHA bins. Reduced chi-squared = 2.9023 for 48 degrees of freedom Null hypothesis probability = 7.845721e-11 XSPEC12>iplot PLT> wdata fexvii_15.014_rgs12_iso_h2.dat - h_inf = 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 9.95949 +/- 1.59997 3 1 windprof u0 0.619942 +/- 1.77602E-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 15.0140 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 6.44393E-04 +/- 1.06770E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.85197E-03 +/- 2.43867E-04 ________________________________________________________________________ Chi-Squared = 175.56 using 52 PHA bins. Reduced chi-squared = 3.6576 for 48 degrees of freedom Null hypothesis probability = 1.959389e-16 XSPEC12>iplot PLT> wdata fexvii_15.014_rgs12_iso_h5.dat