(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 8.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     56430,       1274]
NotebookOptionsPosition[     55543,       1237]
NotebookOutlinePosition[     55970,       1254]
CellTagsIndexPosition[     55927,       1251]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Initialization", "Section",
 CellChangeTimes->{{3.57938348280345*^9, 3.579383494443326*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"yr", " ", "=", " ", 
   RowBox[{"3600", "*", "24", "*", "365.25"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pc", " ", "=", " ", 
   RowBox[{"3.086", " ", 
    RowBox[{"10", "^", "18"}]}]}], ";"}], "\n", 
 RowBox[{
  RowBox[{"k", ":=", 
   RowBox[{"1.38", " ", 
    RowBox[{"10", "^", 
     RowBox[{"-", "16"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mH", ":=", 
   RowBox[{"1.67", " ", 
    RowBox[{"10", "^", 
     RowBox[{"-", "24"}]}]}]}], ";"}], "\n", 
 RowBox[{"Tste", ":=", 
  RowBox[{"{", 
   RowBox[{
   "51200", ",", "48700", ",", "47400", ",", "46100", ",", "44800", ",", 
    "43600", ",", "42300", ",", "41000", ",", "39700", ",", "38400", ",", 
    "37200", ",", "35900", ",", "34600", ",", " ", "33300", ",", " ", 
    "32000"}], "}"}]}], "\n", 
 RowBox[{"logQ", ":=", 
  RowBox[{"{", 
   RowBox[{
   "49.87", ",", "49.70", ",", "49.61", ",", "49.53", ",", "49.43", ",", 
    "49.34", ",", "49.23", ",", "49.12", ",", "49.00", ",", "48.87", ",", 
    "48.72", ",", "48.56", ",", "48.38", ",", " ", "48.18", ",", " ", 
    "47.90"}], "}"}]}], "\[IndentingNewLine]", 
 RowBox[{"Qint", ":=", 
  RowBox[{"Interpolation", "[", 
   RowBox[{"Transpose", "[", 
    RowBox[{"N", "@", 
     RowBox[{"{", 
      RowBox[{"Tste", ",", 
       RowBox[{"(", 
        RowBox[{"10", "^", "logQ"}], ")"}]}], "}"}]}], "]"}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{"\[Alpha]Bcoeff", "[", "T_", "]"}], ":=", 
  RowBox[{"2.59", " ", 
   RowBox[{"10", "^", 
    RowBox[{"-", "13"}]}], " ", "*", " ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"T", "/", " ", 
      RowBox[{"10", "^", "4"}]}], ")"}], "^", 
    RowBox[{"(", 
     RowBox[{"-", ".7"}], ")"}]}]}]}], "\n", 
 RowBox[{
  RowBox[{"csform", "[", "T_", "]"}], ":=", 
  RowBox[{"Sqrt", "[", 
   RowBox[{"k", "*", 
    RowBox[{"T", "/", 
     RowBox[{"(", "mH", ")"}]}]}], "]"}]}], "\n", 
 RowBox[{
  RowBox[{"rs", "[", 
   RowBox[{"Q0_", ",", "nH_", ",", "\[Alpha]B_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"3", " ", 
     RowBox[{"Q0", "/", 
      RowBox[{"(", 
       RowBox[{"4", "\[Pi]", " ", 
        RowBox[{
         RowBox[{"(", "nH", ")"}], "^", "2"}], " ", 
        RowBox[{"(", "\[Alpha]B", ")"}]}], ")"}]}]}], ")"}], "^", 
   RowBox[{"(", 
    RowBox[{"1", "/", "3"}], ")"}]}]}], "\n", 
 RowBox[{
  RowBox[{"trec", "[", 
   RowBox[{"nH_", ",", "\[Alpha]B_"}], "]"}], ":=", 
  RowBox[{"1", "/", 
   RowBox[{"(", 
    RowBox[{"\[Alpha]B", "  ", "*", "nH"}], ")"}]}]}], "\n", 
 RowBox[{
  RowBox[{"rpos", "[", 
   RowBox[{"t_", ",", " ", "Q0_", ",", "nH_", ",", "\[Alpha]B_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"rs", "[", 
    RowBox[{"Q0", ",", "nH", ",", "\[Alpha]B"}], "]"}], "*", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"1", "-", 
      RowBox[{"Exp", "[", 
       RowBox[{
        RowBox[{"-", "\[Alpha]B"}], " ", "*", " ", "nH", "*", " ", "t"}], 
       "]"}]}], ")"}], "^", 
    RowBox[{"(", 
     RowBox[{"1", "/", "3"}], ")"}]}]}]}], "\n", 
 RowBox[{
  RowBox[{"ts", "[", 
   RowBox[{"csb_", ",", "Q0b_", ",", "nHb_", ",", "\[Alpha]Bb_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"rs", "[", 
    RowBox[{"Q0b", ",", "nHb", ",", "\[Alpha]Bb"}], "]"}], "/", 
   "csb"}]}], "\n", 
 RowBox[{
  RowBox[{"rd", "[", 
   RowBox[{"csa_", ",", "Q0a_", ",", "nHa_", ",", "\[Alpha]Ba_", ",", "ta_"}],
    "]"}], ":=", 
  RowBox[{
   RowBox[{"rs", "[", 
    RowBox[{"Q0a", ",", "nHa", ",", "\[Alpha]Ba"}], "]"}], "*", 
   RowBox[{"Power", "[", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"7", "/", 
        RowBox[{"Sqrt", "[", "12", "]"}]}], "*", 
       RowBox[{"ta", "/", 
        RowBox[{"ts", "[", 
         RowBox[{"csa", ",", "Q0a", ",", "nHa", ",", "\[Alpha]Ba"}], 
         "]"}]}]}], ")"}], ",", 
     RowBox[{"(", 
      RowBox[{"4", "/", "7"}], ")"}]}], "]"}]}]}]}], "Input",
 InitializationCell->True,
 CellChangeTimes->{3.579335569498302*^9},
 CellID->319438733]
}, Open  ]],

Cell[CellGroupData[{

Cell["R-type Expansion", "Section",
 CellChangeTimes->{{3.57938348280345*^9, 3.579383512726465*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"Graphics", "[", 
    RowBox[{
     RowBox[{"{", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"Thick", ",", "Dashed", ",", 
         RowBox[{"Circle", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", "0"}], "}"}], ",", 
           RowBox[{
            RowBox[{"1", "/", "pc"}], "*", 
            RowBox[{"rs", "[", 
             RowBox[{
              RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
              RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}]}]}], 
          "]"}]}], "}"}], ",", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{"Red", ",", 
         RowBox[{"Disk", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", "0"}], "}"}], ",", 
           RowBox[{
            RowBox[{"1", "/", "pc"}], "*", 
            RowBox[{"rpos", "[", 
             RowBox[{
              RowBox[{"t", "*", "yr"}], ",", 
              RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
              RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}]}]}], 
          "]"}]}], "}"}]}], "}"}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"\"\<x [pc]\>\"", ",", "\"\<y [pc]\>\""}], "}"}]}], ",", 
     RowBox[{"LabelStyle", "\[Rule]", "Bold"}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "Teff", ",", "32000", ",", "\"\<Stellar eff. temperature [K]\>\""}], 
      "}"}], ",", "32000", ",", "51000"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"Te", ",", "10000", ",", "\"\<Electron temperature [K]\>\""}], 
      "}"}], ",", "5000", ",", "50000"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "nH", ",", "1", ",", 
       "\"\<Neutral density  [\!\(\*SuperscriptBox[\(cm\), \(-3\)]\)]\>\""}], 
      "}"}], ",", "1", ",", "100"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"t", ",", 
       RowBox[{".01", "*", 
        RowBox[{
         RowBox[{"trec", "[", 
          RowBox[{"nH", ",", 
           RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", "yr"}]}],
        ",", " ", "\"\<Time [yr]\>\""}], "}"}], ",", 
     RowBox[{".01", "*", 
      RowBox[{
       RowBox[{"trec", "[", 
        RowBox[{"nH", ",", 
         RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", "yr"}]}], 
     ",", 
     RowBox[{"10", "*", 
      RowBox[{
       RowBox[{"trec", "[", 
        RowBox[{"nH", ",", 
         RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", "yr"}]}], 
     ",", 
     RowBox[{"AppearanceElements", "\[Rule]", "All"}]}], "}"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"Dynamic", "[", "\[IndentingNewLine]", 
    RowBox[{"LogLinearPlot", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"rpos", "[", 
        RowBox[{
         RowBox[{"tt", "*", "yr"}], ",", 
         RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
         RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", "pc"}], 
      ",", 
      RowBox[{"{", 
       RowBox[{"tt", ",", 
        RowBox[{".01", "*", 
         RowBox[{
          RowBox[{"trec", "[", 
           RowBox[{"nH", ",", 
            RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
          "yr"}]}], ",", 
        RowBox[{"10", "*", 
         RowBox[{
          RowBox[{"trec", "[", 
           RowBox[{"nH", ",", 
            RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
          "yr"}]}]}], "}"}], ",", 
      RowBox[{"AxesLabel", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{"\"\<t [yr]\>\"", ",", "\"\<R [pc]\>\""}], "}"}]}], ",", 
      RowBox[{"LabelStyle", "\[Rule]", "Bold"}], ",", 
      RowBox[{"ImageSize", "\[Rule]", "300"}], ",", 
      RowBox[{"AspectRatio", "\[Rule]", ".5"}], ",", 
      RowBox[{"Epilog", "\[Rule]", 
       RowBox[{"{", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"Point", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"Log", "@", "t"}], ",", 
            RowBox[{
             RowBox[{"rpos", "[", 
              RowBox[{
               RowBox[{"t", "*", "yr"}], ",", 
               RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
               RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
             "pc"}]}], "}"}], "]"}], ",", "\[IndentingNewLine]", 
         RowBox[{"{", 
          RowBox[{"Dashed", ",", " ", 
           RowBox[{"Line", "[", 
            RowBox[{"{", "\[IndentingNewLine]", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"Log", "@", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"trec", "[", 
                    RowBox[{"nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                   "yr"}], ")"}]}], ",", 
                RowBox[{
                 RowBox[{"rs", "[", 
                  RowBox[{
                   RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                   RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                 "pc"}]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"Log", "@", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"trec", "[", 
                    RowBox[{"nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                   "yr"}], ")"}]}], ",", 
                RowBox[{
                 RowBox[{
                  RowBox[{"rs", "[", 
                   RowBox[{
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                  "pc"}], "*", ".1"}]}], "}"}]}], "\[IndentingNewLine]", 
             "}"}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", 
         RowBox[{"{", 
          RowBox[{"Text", "[", 
           RowBox[{"\"\< \!\(\*SubscriptBox[\(t\), \(rec\)]\)\>\"", ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"Log", "[", 
               RowBox[{"1.3", "*", 
                RowBox[{
                 RowBox[{"trec", "[", 
                  RowBox[{"nH", ",", 
                   RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                 "yr"}]}], "]"}], ",", 
              RowBox[{
               RowBox[{
                RowBox[{"rs", "[", 
                 RowBox[{
                  RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                  RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                "pc"}], " ", "*", ".75"}]}], "}"}]}], "]"}], 
          "\[IndentingNewLine]", "}"}]}], "}"}]}]}], "]"}], "]"}], 
   "\[IndentingNewLine]", ",", "\[IndentingNewLine]", 
   RowBox[{"ControlPlacement", "\[Rule]", "Left"}]}], "]"}]], "Input",
 CellOpen->False,
 CellChangeTimes->{{3.5793320316026087`*^9, 3.579332037502488*^9}, {
  3.579332410150779*^9, 3.579332410320136*^9}, {3.579332464588944*^9, 
  3.579332466665523*^9}}],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`nH$$ = 
    1, $CellContext`t$$ = (
     1.9390824192495932`/$CellContext`nH$$) $CellContext`Te$$^0.7, \
$CellContext`Te$$ = 10000, $CellContext`Teff$$ = 32000, Typeset`show$$ = True,
     Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", 
    Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = 
    "\"untitled\"", Typeset`specs$$ = {{{
       Hold[$CellContext`Teff$$], 32000, "Stellar eff. temperature [K]"}, 
      32000, 51000}, {{
       Hold[$CellContext`Te$$], 10000, "Electron temperature [K]"}, 5000, 
      50000}, {{
       Hold[$CellContext`nH$$], 1, 
       "Neutral density  [\!\(\*SuperscriptBox[\(cm\), \(-3\)]\)]"}, 1, 
      100}, {{
       Hold[$CellContext`t$$], (
        1.9390824192495932`/$CellContext`nH$$) $CellContext`Te$$^0.7, 
       "Time [yr]"}, 
      Dynamic[0.01 ($CellContext`trec[$CellContext`nH$$, 
          $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr)], 
      Dynamic[10 ($CellContext`trec[$CellContext`nH$$, 
          $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr)]}, {
      Hold[
       Dynamic[
        LogLinearPlot[$CellContext`rpos[$CellContext`tt $CellContext`yr, 
           $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
           $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`pc, \
{$CellContext`tt, 0.01 ($CellContext`trec[$CellContext`nH$$, 
             $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr),
           10 ($CellContext`trec[$CellContext`nH$$, 
             $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr)}\
, AxesLabel -> {"t [yr]", "R [pc]"}, LabelStyle -> Bold, ImageSize -> 300, 
         AspectRatio -> 0.5, Epilog -> {
           Point[{
             
             Log[$CellContext`t$$], $CellContext`rpos[$CellContext`t$$ \
$CellContext`yr, 
               $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
               $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`\
pc}], {Dashed, 
            Line[{{
               Log[$CellContext`trec[$CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], $CellContext`rs[
                 $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                 $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`\
pc}, {
               Log[$CellContext`trec[$CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], ($CellContext`rs[
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`pc) 0.1}}]}, {
            Text[" \!\(\*SubscriptBox[\(t\), \(rec\)]\)", {
              Log[1.3 ($CellContext`trec[$CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr)], ($CellContext`rs[
                 $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                 $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`\
pc) 0.75}]}}]]], Manipulate`Dump`ThisIsNotAControl}}, Typeset`size$$ = {
    310., {149., 154.}}, Typeset`update$$ = 0, Typeset`initDone$$, 
    Typeset`skipInitDone$$ = True, $CellContext`Teff$13054$$ = 
    0, $CellContext`Te$13055$$ = 0, $CellContext`nH$13056$$ = 
    0, $CellContext`t$13057$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, 
      "Variables" :> {$CellContext`nH$$ = 
        1, $CellContext`t$$ = (
          1.9390824192495932`/$CellContext`nH$$) $CellContext`Te$$^0.7, \
$CellContext`Te$$ = 10000, $CellContext`Teff$$ = 32000}, 
      "ControllerVariables" :> {
        Hold[$CellContext`Teff$$, $CellContext`Teff$13054$$, 0], 
        Hold[$CellContext`Te$$, $CellContext`Te$13055$$, 0], 
        Hold[$CellContext`nH$$, $CellContext`nH$13056$$, 0], 
        Hold[$CellContext`t$$, $CellContext`t$13057$$, 0]}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> Graphics[{{Thick, Dashed, 
          Circle[{0, 0}, (1/$CellContext`pc) $CellContext`rs[
             $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
             $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]}, {Red, 
          
          Disk[{0, 
           0}, (1/$CellContext`pc) $CellContext`rpos[$CellContext`t$$ \
$CellContext`yr, 
             $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
             $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]}}, Frame -> 
        True, FrameLabel -> {"x [pc]", "y [pc]"}, LabelStyle -> Bold], 
      "Specifications" :> {{{$CellContext`Teff$$, 32000, 
          "Stellar eff. temperature [K]"}, 32000, 
         51000}, {{$CellContext`Te$$, 10000, "Electron temperature [K]"}, 
         5000, 50000}, {{$CellContext`nH$$, 1, 
          "Neutral density  [\!\(\*SuperscriptBox[\(cm\), \(-3\)]\)]"}, 1, 
         100}, {{$CellContext`t$$, (
           1.9390824192495932`/$CellContext`nH$$) $CellContext`Te$$^0.7, 
          "Time [yr]"}, 
         Dynamic[0.01 ($CellContext`trec[$CellContext`nH$$, 
             $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr)]\
, 
         Dynamic[10 ($CellContext`trec[$CellContext`nH$$, 
             $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr)]\
, AppearanceElements -> All}, 
        Dynamic[
         LogLinearPlot[$CellContext`rpos[$CellContext`tt $CellContext`yr, 
            $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
            $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`pc, \
{$CellContext`tt, 0.01 ($CellContext`trec[$CellContext`nH$$, 
              $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr)\
, 10 ($CellContext`trec[$CellContext`nH$$, 
              $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr)\
}, AxesLabel -> {"t [yr]", "R [pc]"}, LabelStyle -> Bold, ImageSize -> 300, 
          AspectRatio -> 0.5, Epilog -> {
            Point[{
              
              Log[$CellContext`t$$], $CellContext`rpos[$CellContext`t$$ \
$CellContext`yr, 
                $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`\
pc}], {Dashed, 
             Line[{{
                Log[$CellContext`trec[$CellContext`nH$$, 
                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], $CellContext`rs[
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`pc}, {
                Log[$CellContext`trec[$CellContext`nH$$, 
                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], ($CellContext`rs[
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`pc) 0.1}}]}, {
             Text[" \!\(\*SubscriptBox[\(t\), \(rec\)]\)", {
               Log[1.3 ($CellContext`trec[$CellContext`nH$$, 
                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr)], ($CellContext`rs[
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`pc) 0.75}]}}]]}, "Options" :> {ControlPlacement -> Left}, 
      "DefaultOptions" :> {}],
     ImageSizeCache->{686., {200., 206.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{3.579335607083054*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["D-type Expansion", "Section",
 CellChangeTimes->{{3.57938348280345*^9, 3.5793835240006638`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Manipulate", "[", 
  RowBox[{
   RowBox[{"Graphics", "[", 
    RowBox[{
     RowBox[{"{", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"Red", ",", 
         RowBox[{"Disk", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", "0"}], "}"}], ",", 
           RowBox[{
            RowBox[{"1", "/", "pc"}], "*", 
            RowBox[{"rd", "[", 
             RowBox[{
              RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
              RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
              RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}], ",", 
              RowBox[{"t", "*", "yr"}]}], "]"}]}]}], "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"Thick", ",", "Dashed", ",", 
         RowBox[{"Circle", "[", 
          RowBox[{
           RowBox[{"{", 
            RowBox[{"0", ",", "0"}], "}"}], ",", 
           RowBox[{
            RowBox[{"1", "/", "pc"}], "*", 
            RowBox[{"rs", "[", 
             RowBox[{
              RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
              RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}]}]}], "]"}], 
         ",", 
         RowBox[{"Text", "[", "\[IndentingNewLine]", 
          RowBox[{
           RowBox[{"StringForm", "[", 
            RowBox[{
            "\"\<Initial \!\(\*SubscriptBox[\(R\), \(s\)]\) = `` pc\>\"", ",", 
             RowBox[{
              RowBox[{"rs", "[", 
               RowBox[{
                RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
              "pc"}]}], "]"}], ",", 
           RowBox[{"{", 
            RowBox[{"0", ",", "0"}], "}"}]}], "]"}]}], "}"}]}], "}"}], ",", 
     RowBox[{"Frame", "\[Rule]", "True"}], ",", 
     RowBox[{"FrameLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"\"\<x [pc]\>\"", ",", "\"\<y [pc]\>\""}], "}"}]}], ",", 
     "\[IndentingNewLine]", 
     RowBox[{"LabelStyle", "\[Rule]", "Bold"}]}], "]"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "Teff", ",", "32000", ",", "\"\<Stellar eff. temperature [K]\>\""}], 
      "}"}], ",", "32000", ",", "51000"}], "}"}], ",", "\[IndentingNewLine]", 
   
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"Te", ",", "10000", ",", "\"\<Electron temperature [K]\>\""}], 
      "}"}], ",", "5000", ",", "50000"}], "}"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "nH", ",", "1", ",", 
       "\"\<Neutral density  [\!\(\*SuperscriptBox[\(cm\), \(-3\)]\)]\>\""}], 
      "}"}], ",", "1", ",", "100"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"t", ",", 
       RowBox[{
        RowBox[{
         RowBox[{"Sqrt", "[", "12", "]"}], "/", "7"}], "*", 
        RowBox[{
         RowBox[{"ts", "[", 
          RowBox[{
           RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
           RowBox[{"Qint", "[", "Teff", "]"}], ",", " ", "nH", ",", " ", 
           RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", "yr"}]}],
        ",", "\"\<Time [yr]\>\""}], "}"}], ",", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"trec", "[", 
         RowBox[{"nH", ",", 
          RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", "yr"}], 
       "*", ".9"}], ")"}], ",", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"5", "/", "yr"}], "*", 
       RowBox[{"ts", "[", 
        RowBox[{
         RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
         RowBox[{"Qint", "[", "Teff", "]"}], ",", " ", "nH", ",", " ", 
         RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}]}], ")"}], ",", 
     RowBox[{"AppearanceElements", "\[Rule]", "All"}]}], "}"}], ",", 
   "\[IndentingNewLine]", 
   RowBox[{"Dynamic", "[", "\[IndentingNewLine]", 
    RowBox[{"LogLinearPlot", "[", 
     RowBox[{
      RowBox[{
       RowBox[{"1", "/", "pc"}], "*", 
       RowBox[{"rd", "[", 
        RowBox[{
         RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
         RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
         RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}], ",", 
         RowBox[{"tt", "*", "yr"}]}], "]"}]}], ",", 
      RowBox[{"{", 
       RowBox[{"tt", ",", 
        RowBox[{
         RowBox[{
          RowBox[{"trec", "[", 
           RowBox[{"nH", ",", 
            RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", "yr"}], 
         "*", ".9"}], ",", 
        RowBox[{
         RowBox[{"5", "/", "yr"}], "*", 
         RowBox[{"ts", "[", 
          RowBox[{
           RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
           RowBox[{"Qint", "[", "Teff", "]"}], ",", " ", "nH", ",", " ", 
           RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}]}]}], "}"}], 
      ",", "\[IndentingNewLine]", 
      RowBox[{"AxesLabel", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{"\"\<t [yr]\>\"", ",", "\"\<R [pc]\>\""}], "}"}]}], ",", 
      RowBox[{"LabelStyle", "\[Rule]", "Bold"}], ",", 
      RowBox[{"ImageSize", "\[Rule]", "350"}], ",", 
      RowBox[{"AspectRatio", "\[Rule]", ".5"}], ",", "\[IndentingNewLine]", 
      RowBox[{"Epilog", "\[Rule]", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{"Point", "[", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"Log", "@", "t"}], ",", 
            RowBox[{
             RowBox[{"rd", "[", 
              RowBox[{
               RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
               RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
               RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}], ",", 
               RowBox[{"t", "*", "yr"}]}], "]"}], "/", "pc"}]}], "}"}], "]"}],
          ",", "\[IndentingNewLine]", 
         RowBox[{"{", 
          RowBox[{"Dashed", ",", " ", 
           RowBox[{"Line", "[", 
            RowBox[{"{", "\[IndentingNewLine]", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"Log", "@", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"ts", "[", 
                    RowBox[{
                    RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", " ", "nH", ",", 
                    " ", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                   "yr"}], ")"}]}], ",", 
                RowBox[{
                 RowBox[{
                  RowBox[{"rs", "[", 
                   RowBox[{
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                  "pc"}], "*", ".01"}]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"Log", "@", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"ts", "[", 
                    RowBox[{
                    RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", " ", "nH", ",", 
                    " ", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                   "yr"}], ")"}]}], ",", 
                RowBox[{
                 RowBox[{
                  RowBox[{"rs", "[", 
                   RowBox[{
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                  "pc"}], "*", "2"}]}], "}"}]}], "}"}], "\[IndentingNewLine]",
             "]"}]}], "}"}], ",", "\[IndentingNewLine]", 
         RowBox[{"{", 
          RowBox[{"Text", "[", 
           RowBox[{"\"\< \!\(\*SubscriptBox[\(t\), \(s\)]\)\>\"", ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"Log", "[", 
               RowBox[{
                RowBox[{
                 RowBox[{"ts", "[", 
                  RowBox[{
                   RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                   RowBox[{"Qint", "[", "Teff", "]"}], ",", " ", "nH", ",", 
                   " ", 
                   RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                 "yr"}], "*", ".9"}], "]"}], ",", 
              RowBox[{
               RowBox[{
                RowBox[{"rs", "[", 
                 RowBox[{
                  RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                  RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                "pc"}], " ", "*", "1.75"}]}], "}"}]}], "]"}], "}"}], 
         "\[IndentingNewLine]", ",", "\[IndentingNewLine]", 
         RowBox[{"{", 
          RowBox[{"Dashed", ",", " ", 
           RowBox[{"Line", "[", 
            RowBox[{"{", "\[IndentingNewLine]", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"Log", "@", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"trec", "[", 
                    RowBox[{"nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                   "yr"}], ")"}]}], ",", 
                RowBox[{
                 RowBox[{
                  RowBox[{"rd", "[", 
                   RowBox[{
                    RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}], ",", 
                    RowBox[{"trec", "[", 
                    RowBox[{"nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}]}], 
                   "]"}], "/", "pc"}], "*", "1.3"}]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"Log", "@", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{"trec", "[", 
                    RowBox[{"nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                   "yr"}], ")"}]}], ",", 
                RowBox[{
                 RowBox[{
                  RowBox[{"rd", "[", 
                   RowBox[{
                    RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}], ",", 
                    RowBox[{"trec", "[", 
                    RowBox[{"nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}]}], 
                   "]"}], "/", "pc"}], "*", ".1"}]}], "}"}]}], 
             "\[IndentingNewLine]", "}"}], "]"}]}], "}"}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{"{", 
          RowBox[{"Dashed", ",", " ", 
           RowBox[{"Line", "[", 
            RowBox[{"{", "\[IndentingNewLine]", 
             RowBox[{
              RowBox[{"{", 
               RowBox[{
                RowBox[{"Log", "@", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"Sqrt", "[", "12", "]"}], "/", "7"}], "*", 
                   RowBox[{
                    RowBox[{"ts", "[", 
                    RowBox[{
                    RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", " ", "nH", ",", 
                    " ", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                    "yr"}]}], ")"}]}], ",", 
                RowBox[{
                 RowBox[{"rd", "[", 
                  RowBox[{
                   RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                   RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                   RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}], ",", 
                   RowBox[{
                    RowBox[{
                    RowBox[{"Sqrt", "[", "12", "]"}], "/", "7"}], "*", 
                    RowBox[{"ts", "[", 
                    RowBox[{
                    RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", " ", "nH", ",", 
                    " ", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}]}]}], 
                  "]"}], "/", "pc"}]}], "}"}], ",", 
              RowBox[{"{", 
               RowBox[{
                RowBox[{"Log", "@", 
                 RowBox[{"(", 
                  RowBox[{
                   RowBox[{
                    RowBox[{"Sqrt", "[", "12", "]"}], "/", "7"}], "*", 
                   RowBox[{
                    RowBox[{"ts", "[", 
                    RowBox[{
                    RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", " ", "nH", ",", 
                    " ", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                    "yr"}]}], ")"}]}], ",", 
                RowBox[{
                 RowBox[{
                  RowBox[{"rd", "[", 
                   RowBox[{
                    RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}], ",", 
                    RowBox[{
                    RowBox[{
                    RowBox[{"Sqrt", "[", "12", "]"}], "/", "7"}], "*", 
                    RowBox[{"ts", "[", 
                    RowBox[{
                    RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                    RowBox[{"Qint", "[", "Teff", "]"}], ",", " ", "nH", ",", 
                    " ", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}]}]}], 
                   "]"}], "/", "pc"}], "*", ".1"}]}], "}"}]}], "}"}], 
            "\[IndentingNewLine]", "]"}]}], "}"}], ",", "\[IndentingNewLine]", 
         RowBox[{"{", 
          RowBox[{"Text", "[", 
           RowBox[{"\"\< \!\(\*SubscriptBox[\(t\), \(rec\)]\)\>\"", ",", 
            RowBox[{"{", 
             RowBox[{
              RowBox[{"Log", "[", 
               RowBox[{"1.2", "*", 
                RowBox[{
                 RowBox[{"trec", "[", 
                  RowBox[{"nH", ",", 
                   RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}], "/", 
                 "yr"}]}], "]"}], ",", 
              RowBox[{"1.8", "*", 
               RowBox[{
                RowBox[{"rd", "[", 
                 RowBox[{
                  RowBox[{"csform", "[", "Te", "]"}], ",", " ", 
                  RowBox[{"Qint", "[", "Teff", "]"}], ",", "nH", ",", 
                  RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}], ",", 
                  RowBox[{"trec", "[", 
                   RowBox[{"nH", ",", 
                    RowBox[{"\[Alpha]Bcoeff", "[", "Te", "]"}]}], "]"}]}], 
                 "]"}], "/", "pc"}]}]}], "}"}]}], "]"}], 
          "\[IndentingNewLine]", "}"}]}], "\[IndentingNewLine]", "}"}]}]}], 
     "]"}], "]"}], "\[IndentingNewLine]", ",", 
   RowBox[{"ControlPlacement", "\[Rule]", "Left"}]}], " ", "]"}]], "Input",
 CellOpen->False,
 CellChangeTimes->{{3.579332512508893*^9, 3.579332544828951*^9}, {
  3.5793356415693493`*^9, 3.579335641943233*^9}},
 CellID->1581038721],

Cell[BoxData[
 TagBox[
  StyleBox[
   DynamicModuleBox[{$CellContext`nH$$ = 
    1, $CellContext`t$$ = (
     1.957392187082916*^-9 $CellContext`Te$$^
      Rational[-1, 2]) (($CellContext`nH$$^(-2) $CellContext`Te$$^0.7) 
      InterpolatingFunction[{{32000., 51200.}}, {
        4, 7, 0, {15}, {4}, 0, 0, 0, 0, Automatic}, {{32000., 33300., 34600., 
        35900., 37200., 38400., 39700., 41000., 42300., 43600., 44800., 
        46100., 47400., 48700., 51200.}}, {
        Developer`PackedArrayForm, {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 
         13, 14, 15}, {7.943282347242788*^47, 1.5135612484362073`*^48, 
         2.3988329190195045`*^48, 3.6307805477010326`*^48, 
         5.248074602497712*^48, 7.413102413009131*^48, 1.*^49, 
         1.3182567385563993`*^49, 1.698243652461732*^49, 
         2.1877616239495697`*^49, 2.691534803926914*^49, 
         3.388441561392034*^49, 4.073802778041122*^49, 5.011872336272756*^49, 
         7.413102413009131*^49}}, {Automatic}][$CellContext`Teff$$])^
     Rational[1, 3], $CellContext`Te$$ = 10000, $CellContext`Teff$$ = 32000, 
    Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, 
    Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 
    1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{{
       Hold[$CellContext`Teff$$], 32000, "Stellar eff. temperature [K]"}, 
      32000, 51000}, {{
       Hold[$CellContext`Te$$], 10000, "Electron temperature [K]"}, 5000, 
      50000}, {{
       Hold[$CellContext`nH$$], 1, 
       "Neutral density  [\!\(\*SuperscriptBox[\(cm\), \(-3\)]\)]"}, 1, 
      100}, {{
       Hold[$CellContext`t$$], (
        1.957392187082916*^-9 $CellContext`Te$$^
         Rational[-1, 2]) (($CellContext`nH$$^(-2) $CellContext`Te$$^0.7) 
         InterpolatingFunction[{{32000., 51200.}}, {
           4, 7, 0, {15}, {4}, 0, 0, 0, 0, Automatic}, {{32000., 33300., 
           34600., 35900., 37200., 38400., 39700., 41000., 42300., 43600., 
           44800., 46100., 47400., 48700., 51200.}}, {
           Developer`PackedArrayForm, {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 
            12, 13, 14, 15}, {7.943282347242788*^47, 1.5135612484362073`*^48, 
            2.3988329190195045`*^48, 3.6307805477010326`*^48, 
            5.248074602497712*^48, 7.413102413009131*^48, 1.*^49, 
            1.3182567385563993`*^49, 1.698243652461732*^49, 
            2.1877616239495697`*^49, 2.691534803926914*^49, 
            3.388441561392034*^49, 4.073802778041122*^49, 
            5.011872336272756*^49, 7.413102413009131*^49}}, {
           Automatic}][$CellContext`Teff$$])^Rational[1, 3], "Time [yr]"}, 
      Dynamic[($CellContext`trec[$CellContext`nH$$, 
          $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr) 
       0.9], 
      Dynamic[(5/$CellContext`yr) $CellContext`ts[
         $CellContext`csform[$CellContext`Te$$], 
         $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
         $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]}, {
      Hold[
       Dynamic[
        LogLinearPlot[(1/$CellContext`pc) $CellContext`rd[
           $CellContext`csform[$CellContext`Te$$], 
           $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
           $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], $CellContext`tt \
$CellContext`yr], {$CellContext`tt, ($CellContext`trec[$CellContext`nH$$, 
             $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr) 
          0.9, (5/$CellContext`yr) $CellContext`ts[
            $CellContext`csform[$CellContext`Te$$], 
            $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
            $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]}, 
         AxesLabel -> {"t [yr]", "R [pc]"}, LabelStyle -> Bold, ImageSize -> 
         350, AspectRatio -> 0.5, Epilog -> {
           Point[{
             Log[$CellContext`t$$], $CellContext`rd[
               $CellContext`csform[$CellContext`Te$$], 
               $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
               $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], \
$CellContext`t$$ $CellContext`yr]/$CellContext`pc}], {Dashed, 
            Line[{{
               Log[$CellContext`ts[
                  $CellContext`csform[$CellContext`Te$$], 
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], ($CellContext`rs[
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`pc) 0.01}, {
               Log[$CellContext`ts[
                  $CellContext`csform[$CellContext`Te$$], 
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], ($CellContext`rs[
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`pc) 2}}]}, {
            Text[" \!\(\*SubscriptBox[\(t\), \(s\)]\)", {
              Log[($CellContext`ts[
                  $CellContext`csform[$CellContext`Te$$], 
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr) 0.9], ($CellContext`rs[
                 $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                 $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`\
pc) 1.75}]}, {Dashed, 
            Line[{{
               Log[$CellContext`trec[$CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], ($CellContext`rd[
                  $CellContext`csform[$CellContext`Te$$], 
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], 
                  $CellContext`trec[$CellContext`nH$$, 
                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]/$\
CellContext`pc) 1.3}, {
               Log[$CellContext`trec[$CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], ($CellContext`rd[
                  $CellContext`csform[$CellContext`Te$$], 
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], 
                  $CellContext`trec[$CellContext`nH$$, 
                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]/$\
CellContext`pc) 0.1}}]}, {Dashed, 
            Line[{{
               Log[(Sqrt[12]/7) ($CellContext`ts[
                   $CellContext`csform[$CellContext`Te$$], 
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr)], $CellContext`rd[
                 $CellContext`csform[$CellContext`Te$$], 
                 $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                 $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], (Sqrt[12]/
                  7) $CellContext`ts[
                   $CellContext`csform[$CellContext`Te$$], 
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]/$\
CellContext`pc}, {
               Log[(Sqrt[12]/7) ($CellContext`ts[
                   $CellContext`csform[$CellContext`Te$$], 
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr)], ($CellContext`rd[
                  $CellContext`csform[$CellContext`Te$$], 
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], (Sqrt[12]/
                   7) $CellContext`ts[
                    $CellContext`csform[$CellContext`Te$$], 
                    $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                    $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]/$\
CellContext`pc) 0.1}}]}, {
            Text[" \!\(\*SubscriptBox[\(t\), \(rec\)]\)", {
              Log[1.2 ($CellContext`trec[$CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr)], 1.8 ($CellContext`rd[
                 $CellContext`csform[$CellContext`Te$$], 
                 $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                 $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], 
                 $CellContext`trec[$CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]/$\
CellContext`pc)}]}}]]], Manipulate`Dump`ThisIsNotAControl}}, 
    Typeset`size$$ = {323., {156., 160.}}, Typeset`update$$ = 0, 
    Typeset`initDone$$, Typeset`skipInitDone$$ = 
    True, $CellContext`Teff$13554$$ = 0, $CellContext`Te$13555$$ = 
    0, $CellContext`nH$13556$$ = 0, $CellContext`t$13557$$ = 0}, 
    DynamicBox[Manipulate`ManipulateBoxes[
     1, StandardForm, 
      "Variables" :> {$CellContext`nH$$ = 
        1, $CellContext`t$$ = (
          1.957392187082916*^-9 $CellContext`Te$$^
           Rational[-1, 2]) (($CellContext`nH$$^(-2) $CellContext`Te$$^0.7) 
           InterpolatingFunction[{{32000., 51200.}}, {
             4, 7, 0, {15}, {4}, 0, 0, 0, 0, Automatic}, {{32000., 33300., 
             34600., 35900., 37200., 38400., 39700., 41000., 42300., 43600., 
             44800., 46100., 47400., 48700., 51200.}}, {
             Developer`PackedArrayForm, {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11,
               12, 13, 14, 15}, {7.943282347242788*^47, 
              1.5135612484362073`*^48, 2.3988329190195045`*^48, 
              3.6307805477010326`*^48, 5.248074602497712*^48, 
              7.413102413009131*^48, 1.*^49, 1.3182567385563993`*^49, 
              1.698243652461732*^49, 2.1877616239495697`*^49, 
              2.691534803926914*^49, 3.388441561392034*^49, 
              4.073802778041122*^49, 5.011872336272756*^49, 
              7.413102413009131*^49}}, {Automatic}][$CellContext`Teff$$])^
          Rational[1, 3], $CellContext`Te$$ = 10000, $CellContext`Teff$$ = 
        32000}, "ControllerVariables" :> {
        Hold[$CellContext`Teff$$, $CellContext`Teff$13554$$, 0], 
        Hold[$CellContext`Te$$, $CellContext`Te$13555$$, 0], 
        Hold[$CellContext`nH$$, $CellContext`nH$13556$$, 0], 
        Hold[$CellContext`t$$, $CellContext`t$13557$$, 0]}, 
      "OtherVariables" :> {
       Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, 
        Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, 
        Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$,
         Typeset`skipInitDone$$}, "Body" :> Graphics[{{Red, 
          Disk[{0, 0}, (1/$CellContext`pc) $CellContext`rd[
             $CellContext`csform[$CellContext`Te$$], 
             $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
             $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], $CellContext`t$$ \
$CellContext`yr]]}, {Thick, Dashed, 
          Circle[{0, 0}, (1/$CellContext`pc) $CellContext`rs[
             $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
             $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]], 
          Text[
           StringForm[
           "Initial \!\(\*SubscriptBox[\(R\), \(s\)]\) = `` pc", \
$CellContext`rs[
              $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
              $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`pc]\
, {0, 0}]}}, Frame -> True, FrameLabel -> {"x [pc]", "y [pc]"}, LabelStyle -> 
        Bold], "Specifications" :> {{{$CellContext`Teff$$, 32000, 
          "Stellar eff. temperature [K]"}, 32000, 
         51000}, {{$CellContext`Te$$, 10000, "Electron temperature [K]"}, 
         5000, 50000}, {{$CellContext`nH$$, 1, 
          "Neutral density  [\!\(\*SuperscriptBox[\(cm\), \(-3\)]\)]"}, 1, 
         100}, {{$CellContext`t$$, (
           1.957392187082916*^-9 $CellContext`Te$$^
            Rational[-1, 2]) (($CellContext`nH$$^(-2) $CellContext`Te$$^0.7) 
            InterpolatingFunction[{{32000., 51200.}}, {
              4, 7, 0, {15}, {4}, 0, 0, 0, 0, Automatic}, {{32000., 33300., 
              34600., 35900., 37200., 38400., 39700., 41000., 42300., 43600., 
              44800., 46100., 47400., 48700., 51200.}}, {
              Developer`PackedArrayForm, {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 
               11, 12, 13, 14, 15}, {7.943282347242788*^47, 
               1.5135612484362073`*^48, 2.3988329190195045`*^48, 
               3.6307805477010326`*^48, 5.248074602497712*^48, 
               7.413102413009131*^48, 1.*^49, 1.3182567385563993`*^49, 
               1.698243652461732*^49, 2.1877616239495697`*^49, 
               2.691534803926914*^49, 3.388441561392034*^49, 
               4.073802778041122*^49, 5.011872336272756*^49, 
               7.413102413009131*^49}}, {Automatic}][$CellContext`Teff$$])^
           Rational[1, 3], "Time [yr]"}, 
         Dynamic[($CellContext`trec[$CellContext`nH$$, 
             $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr) 
          0.9], 
         Dynamic[(5/$CellContext`yr) $CellContext`ts[
            $CellContext`csform[$CellContext`Te$$], 
            $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
            $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]], 
         AppearanceElements -> All}, 
        Dynamic[
         LogLinearPlot[(1/$CellContext`pc) $CellContext`rd[
            $CellContext`csform[$CellContext`Te$$], 
            $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
            $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], $CellContext`tt \
$CellContext`yr], {$CellContext`tt, ($CellContext`trec[$CellContext`nH$$, 
              $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$CellContext`yr)
             0.9, (5/$CellContext`yr) $CellContext`ts[
             $CellContext`csform[$CellContext`Te$$], 
             $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
             $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]}, 
          AxesLabel -> {"t [yr]", "R [pc]"}, LabelStyle -> Bold, ImageSize -> 
          350, AspectRatio -> 0.5, Epilog -> {
            Point[{
              Log[$CellContext`t$$], $CellContext`rd[
                $CellContext`csform[$CellContext`Te$$], 
                $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], \
$CellContext`t$$ $CellContext`yr]/$CellContext`pc}], {Dashed, 
             Line[{{
                Log[$CellContext`ts[
                   $CellContext`csform[$CellContext`Te$$], 
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], ($CellContext`rs[
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`pc) 0.01}, {
                Log[$CellContext`ts[
                   $CellContext`csform[$CellContext`Te$$], 
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], ($CellContext`rs[
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`pc) 2}}]}, {
             Text[" \!\(\*SubscriptBox[\(t\), \(s\)]\)", {
               Log[($CellContext`ts[
                   $CellContext`csform[$CellContext`Te$$], 
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr) 0.9], ($CellContext`rs[
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`pc) 1.75}]}, {Dashed, 
             Line[{{
                Log[$CellContext`trec[$CellContext`nH$$, 
                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], ($CellContext`rd[
                   $CellContext`csform[$CellContext`Te$$], 
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], 
                   $CellContext`trec[$CellContext`nH$$, 
                    $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]/$\
CellContext`pc) 1.3}, {
                Log[$CellContext`trec[$CellContext`nH$$, 
                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr], ($CellContext`rd[
                   $CellContext`csform[$CellContext`Te$$], 
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], 
                   $CellContext`trec[$CellContext`nH$$, 
                    $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]/$\
CellContext`pc) 0.1}}]}, {Dashed, 
             Line[{{
                Log[(Sqrt[12]/7) ($CellContext`ts[
                    $CellContext`csform[$CellContext`Te$$], 
                    $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                    $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr)], $CellContext`rd[
                  $CellContext`csform[$CellContext`Te$$], 
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], (Sqrt[12]/
                   7) $CellContext`ts[
                    $CellContext`csform[$CellContext`Te$$], 
                    $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                    $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]/$\
CellContext`pc}, {
                Log[(Sqrt[12]/7) ($CellContext`ts[
                    $CellContext`csform[$CellContext`Te$$], 
                    $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                    $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr)], ($CellContext`rd[
                   $CellContext`csform[$CellContext`Te$$], 
                   $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, \

                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], (Sqrt[12]/
                    7) $CellContext`ts[
                    $CellContext`csform[$CellContext`Te$$], 
                    $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                    $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]/$\
CellContext`pc) 0.1}}]}, {
             Text[" \!\(\*SubscriptBox[\(t\), \(rec\)]\)", {
               Log[1.2 ($CellContext`trec[$CellContext`nH$$, 
                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]/$\
CellContext`yr)], 1.8 ($CellContext`rd[
                  $CellContext`csform[$CellContext`Te$$], 
                  $CellContext`Qint[$CellContext`Teff$$], $CellContext`nH$$, 
                  $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$], 
                  $CellContext`trec[$CellContext`nH$$, 
                   $CellContext`\[Alpha]Bcoeff[$CellContext`Te$$]]]/$\
CellContext`pc)}]}}]]}, "Options" :> {ControlPlacement -> Left}, 
      "DefaultOptions" :> {}],
     ImageSizeCache->{712., {200., 206.}},
     SingleEvaluation->True],
    Deinitialization:>None,
    DynamicModuleValues:>{},
    SynchronousInitialization->True,
    UnsavedVariables:>{Typeset`initDone$$},
    UntrackedVariables:>{Typeset`size$$}], "Manipulate",
   Deployed->True,
   StripOnInput->False],
  Manipulate`InterpretManipulate[1]]], "Output",
 CellChangeTimes->{{3.5793356169117327`*^9, 3.579335643264846*^9}}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{740, 652},
WindowMargins->{{20, Automatic}, {Automatic, 44}},
FrontEndVersion->"8.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (February 23, \
2011)",
StyleDefinitions->FrontEnd`FileName[{"Creative"}, "PastelColor.nb", 
  CharacterEncoding -> "UTF-8"]
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 98, 1, 72, "Section"],
Cell[680, 25, 3936, 120, 278, "Input",
 InitializationCell->True,
 CellID->319438733]
}, Open  ]],
Cell[CellGroupData[{
Cell[4653, 150, 100, 1, 72, "Section"],
Cell[CellGroupData[{
Cell[4778, 155, 7191, 187, 19, "Input",
 CellOpen->False],
Cell[11972, 344, 8115, 155, 433, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[20136, 505, 102, 1, 72, "Section"],
Cell[CellGroupData[{
Cell[20263, 510, 15285, 364, 19, "Input",
 CellOpen->False,
 CellID->1581038721],
Cell[35551, 876, 19964, 357, 433, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)