5 aug 11 reanalyzing the O VIII Ly-alpha line, systematically looking at iso- and aniso- porosity with specific, fixed values of hinf using the wavelength range and continuum level from the old fitting (Cohen et al. 2010) MEG only 18.80 to 19.13 n = 2, norm = 2.40e-3 ======================================================================== 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.99760 +/- 0.722425 3 1 windprof u0 0.822835 +/- 0.653010 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 18.9690 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.70272E-04 +/- 5.20279E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 150.88 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 18.969_nonporous.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 2.45747 3.53594 (-0.549721,0.528756) XSPEC12>error 1. 3 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0.99 3 0.629705 0 (-0.24317,-0.872875) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000342678 0.000399366 (-2.74808e-05,2.92077e-05) ANISOTROPIC POROSITY hinf free: consistent with hinf = 0; recover non-porous fit ======================================================================== 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.04974 +/- 0.849218 3 1 windprof u0 0.819458 +/- 0.607427 4 1 windprof h 1.61356E-02 +/- 0.281450 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 18.9690 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.70388E-04 +/- 5.20233E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 150.83 using 130 PHA bins and 126 degrees of freedom. 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 4.25698 +/- 1.51623 3 1 windprof u0 0.697418 +/- 0.186954 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 18.9690 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.70778E-04 +/- 5.19887E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 152.91 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 18.969_aniso_h0.5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 3.10276 5.44395 (-1.03872,1.30247) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.575194 0.821954 (-0.120931,0.125829) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000344387 0.000398222 (-3.06504e-05,2.31845e-05) 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.22083 +/- 2.36525 3 1 windprof u0 0.638609 +/- 0.132151 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 18.9690 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.62058E-04 +/- 5.30241E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 155.06 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 18.969_aniso_h1.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 3.6095 7.39066 (-1.60986,2.17131) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.544971 0.746257 (-0.0934609,0.107826) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000342982 0.000399782 (-1.908e-05,3.77198e-05) 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 21.0993 +/- 14.6303 3 1 windprof u0 0.575679 +/- 8.23278E-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 18.9690 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.56780E-04 +/- 5.35072E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 164.98 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 18.969_aniso_h5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 10.5963 37.0188 (-10.5002,15.9223) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.506897 0.648037 (-0.0688006,0.07234) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.00034225 0.000399127 (-1.8847e-05,3.80301e-05) ISOTROPIC POROSITY hinf free: hinf = 0.21 is best fit, but the constraints on this parameter are minimal ======================================================================== 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.33199 +/- 2.74632 3 1 windprof u0 0.786469 +/- 0.631625 4 1 windprof h 0.208017 +/- 1.96678 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 18.9690 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.69714E-04 +/- 5.19826E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 150.82 using 130 PHA bins and 126 degrees of freedom. XSPEC12>iplot PLT> wdata 18.969_iso_hfree.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 12.9871, 13.5161 and delta stat 0.908698, 1.10312 but latest trial 13.2038 gives 1.15282 Suggest that you check this result using the steppar command. 2 2.49045 13.2516 (-1.06984,9.6913) XSPEC12>steppar 2 0. 12. 30 C-Statistic Delta taustar C-Statistic 2 196.86 45.903 0 0 172.98 22.022 1 0.4 163.3 12.339 2 0.8 158.05 7.0876 3 1.2 155.03 4.0645 4 1.6 153.22 2.2627 5 2 152.07 1.107 6 2.4 151.04 0.077471 7 2.8 150.85 -0.11463 8 3.2 150.84 -0.12027 9 3.6 150.98 0.021258 10 4 151.06 0.10241 11 4.4 151.13 0.17081 12 4.8 151.1 0.13505 13 5.2 151.08 0.12193 14 5.6 151.1 0.14084 15 6 151.1 0.14368 16 6.4 151.11 0.15418 17 6.8 151.13 0.17142 18 7.2 151.16 0.19444 19 7.6 151.18 0.22239 20 8 151.22 0.25457 21 8.4 151.25 0.29039 22 8.8 151.29 0.32932 23 9.2 151.33 0.37092 24 9.6 151.38 0.41479 25 10 151.42 0.46064 26 10.4 151.47 0.50818 27 10.8 151.52 0.55714 28 11.2 151.57 0.60733 29 11.6 151.62 0.65855 30 12 XSPEC12> error 1. 3 Parameter Confidence Range (1) ***Warning: Parameter pegged at hard limit: 0.99 3 0.549914 0 (-0.20631,-0.756224) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000343666 0.000397991 (-3.03819e-05,2.39439e-05) XSPEC12>log 18.969_iso_hinf_vs_taustar_35by35.log Logging to file:18.969_iso_hinf_vs_taustar_35by35.log XSPEC12>steppar 4 0. 5. 35 2 0. 11. 35 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 3.85478 +/- 1.13767 3 1 windprof u0 0.756437 +/- 0.195913 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 18.9690 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.69599E-04 +/- 5.19323E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 150.88 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 18.969_iso_h0.5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) Apparent non-monotonicity in statistic space detected. Current bracket values 3.75743, 4.89511 and delta stat 0, 1.18524 but latest trial 4.80243 gives 1.5605 Suggest that you check this result using the steppar command. 2 2.75188 4.32627 (-1.00555,0.568835) XSPEC12>steppar 2 2.5 4.5 20 C-Statistic Delta taustar C-Statistic 2 152.36 1.3326 0 2.5 152.22 1.1982 1 2.6 152.04 1.0179 2 2.7 151.88 0.85889 3 2.8 151.74 0.71638 4 2.9 151.61 0.58774 5 3 151.5 0.47087 6 3.1 151.39 0.36292 7 3.2 151.29 0.26242 8 3.3 151.2 0.17217 9 3.4 151.12 0.094602 10 3.5 151.06 0.029781 11 3.6 151 -0.021064 12 3.7 150.94 -0.090725 13 3.8 150.93 -0.097163 14 3.9 150.93 -0.093907 15 4 150.97 -0.052275 16 4.1 151.08 0.057753 17 4.2 151.09 0.061991 18 4.3 151.29 0.25982 19 4.4 151.33 0.30026 20 4.5 XSPEC12> error 1. 3 Parameter Confidence Range (1) 3 0.588354 0.890607 (-0.153658,0.148595) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000342494 0.000398818 (-1.98127e-05,3.65115e-05) 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 4.42406 +/- 1.59133 3 1 windprof u0 0.691722 +/- 0.124236 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 18.9690 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.70734E-04 +/- 5.19009E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 151.04 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 18.969_iso_h1.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 2.95006 6.16047 (-1.47428,1.73613) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.55975 0.818238 (-0.13196,0.126529) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000342458 0.00039871 (-1.97602e-05,3.64921e-05) 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 13.2261 +/- 8.88659 3 1 windprof u0 0.589924 +/- 6.84432E-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 18.9690 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.59699E-04 +/- 5.32385E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 152.10 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 18.969_iso_h5.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 7.4338 21.9702 (-5.79367,8.74274) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.525704 0.649594 (-0.0642375,0.0596525) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000342618 0.000399195 (-1.92197e-05,3.73573e-05) 14may12 CHECKING THE H_INF = 2 CASE - ANISO ======================================================================== 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 7.82684 +/- 4.39243 3 1 windprof u0 0.602749 +/- 0.102531 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 18.9690 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.61707E-04 +/- 5.30194E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 158.30 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 18.969_aniso_h2.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 4.90061 11.7472 (-2.91228,3.93434) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.522172 0.684713 (-0.0800732,0.0824682) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000342808 0.000399636 (-1.89038e-05,3.79246e-05) - ISOTROPIC ======================================================================== 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.85220 +/- 2.65529 3 1 windprof u0 0.640474 +/- 8.92852E-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 18.9690 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.62140E-04 +/- 5.29639E-05 18 2 powerlaw PhoIndex 2.00000 frozen 19 2 powerlaw norm 2.40000E-03 frozen ________________________________________________________________________ C-statistic = 151.11 using 130 PHA bins and 127 degrees of freedom. XSPEC12>iplot PLT> wdata 18.969_iso_h2.dat XSPEC12>error 1. 2 Parameter Confidence Range (1) 2 3.65983 8.73475 (-2.19237,2.88255) XSPEC12>error 1. 3 Parameter Confidence Range (1) 3 0.544115 0.729632 (-0.0963593,0.0891572) XSPEC12>error 1. 17 Parameter Confidence Range (1) 17 0.000342555 0.00039886 (-1.95851e-05,3.67198e-05)