/2\\
    d\ln Z_{\ch{H}}/d\ln T&=3/2+\phi_{\rm dis}/(2kT)\\
    d\ln Z_{\ch{H+}}/d\ln T&=3/2+(\phi_{\rm dis}+2\phi_{\rm ion})/(2kT)\\
    d\ln Z_{\ch{e-}}/d\ln T&=3/2\\
    d^{2}\ln Z_{\ch{H2}}/d\ln T^{2}&=0\\
    d^{2}\ln Z_{\ch{H}}/d\ln T^{2}&=-\phi_{\rm dis}/(2kT)\\
    d^{2}\ln Z_{\ch{H+}}/d\ln T^{2}&=-(\phi_{\rm dis}+2\phi_{\rm ion})/(2kT)\\
    d^{2}\ln Z_{\ch{e-}}/d\ln T^{2}&=0\end{aligned}$$ This would not be the case if one consider the vibrational and rotational degree of freedom of as in [@tomida2013; @tomida2015].

The only derivatives that we need to be calculated are $\left(\frac{\partial\ln n_{i}}{\partial\ln T}\right)_{\rho}$ and $\left(\frac{\partial\ln n_{i}}{\partial\ln\rho}\right)_{T}$. When all species present, i.e., $(\rho,T)\in$, let us re-write Equation \[eqn:nconserv\]-\[eqn:qion\] in log form while using \[eqn:econserv\], $$\begin{aligned}
\label{eqn:dlnndlnx}
\begin{pmatrix}
2n_{\ch{H2}}  &  n_{\ch{H}}  &  n_{\ch{n+}}  \\
-1  &  2  &  0  \\
0  &  -1  &  2
\end{pmatrix}
\begin{pmatrix}
d\ln n_{\ch{H2}}\\
d\ln n_{\ch{H}}\\
d\ln n_{\ch{H+}}
\end{pmatrix}=
\begin{pmatrix}
n_{\rm Htot}d\ln n_{\rm Htot}\\
d\ln q_{\rm dis}\\
d\ln q_{\
/2\\ 
     d\ln Z_{\ch{H}}/d\ln T&=3/2+\phi_{\rm dis}/(2kT)\\ 
     d\ln Z_{\ch{H+}}/d\ln T&=3/2+(\phi_{\rm dis}+2\phi_{\rm ion})/(2kT)\\ 
     d\ln Z_{\ch{e-}}/d\ln T&=3/2\\ 
     d^{2}\ln Z_{\ch{H2}}/d\ln T^{2}&=0\\ 
     d^{2}\ln Z_{\ch{H}}/d\ln T^{2}&=-\phi_{\rm dis}/(2kT)\\ 
     d^{2}\ln Z_{\ch{H+}}/d\ln T^{2}&=-(\phi_{\rm dis}+2\phi_{\rm ion})/(2kT)\\ 
     d^{2}\ln Z_{\ch{e-}}/d\ln T^{2}&=0\end{aligned}$$ This would not be the case if one consider the vibrational and rotational degree of freedom of equally in [ @tomida2013; @tomida2015 ]. 

 The lone derivatives that we need to be account are $ \left(\frac{\partial\ln n_{i}}{\partial\ln T}\right)_{\rho}$ and $ \left(\frac{\partial\ln n_{i}}{\partial\ln\rho}\right)_{T}$. When all coinage present, i.e., $ (\rho, T)\in$, let us re - write Equation \[eqn: nconserv\]-\[eqn: qion\ ] in log phase while using \[eqn: econserv\ ], $ $ \begin{aligned } 
 \label{eqn: dlnndlnx } 
 \begin{pmatrix } 
 2n_{\ch{H2 } }   &   n_{\ch{H } }   &   n_{\ch{n+ } }   \\ 
 -1   &   2   &   0   \\ 
 0   &   -1   &   2 
 \end{pmatrix } 
 \begin{pmatrix } 
 d\ln n_{\ch{H2}}\\ 
 d\ln n_{\ch{H}}\\ 
 d\ln n_{\ch{H+ } } 
 \end{pmatrix}= 
 \begin{pmatrix } 
 n_{\rm Htot}d\ln n_{\rm Htot}\\ 
 d\ln q_{\rm dis}\\ 
 d\ln q_{\
/2\\
    d\ln Z_{\ch{H}}/d\ln T&=3/2+\phi_{\rm dis}/(2kU)\\
    d\ln Z_{\ch{H+}}/d\ln J&=3/2+(\pyi_{\rm dms}+2\phi_{\rm ion})/(2kT)\\
    d\ln Z_{\ch{e-}}/d\ln T&=3/2\\
    d^{2}\ln Z_{\ch{H2}}/d\pn T^{2}&=0\\
    d^{2}\ln Z_{\ch{H}}/d\ln T^{2}&=-\phi_{\rm dir}/(2kT)\\
    d^{2}\lj Z_{\ch{H+}}/d\lb T^{2}&=-(\pii_{\rm dis}+2\phi_{\rm ioi})/(2iT)\\
    d^{2}\ln Z_{\ch{s-}}/f\ln C^{2}&=0\eid{aligned}$$ This eould not te the case if ove consider the vibrational and rotatyonal drggee of freedom of sf in [@nonida2013; @tomida2015].

The only derivativss that we need to br calculated are $\left(\frac{\pwrtiwl\ln n_{i}}{\partial\ln T}\gight)_{\rho}$ ane $\lesr(\frac{\partial\un n_{i}}{\partial\ln\rho}\right)_{J}$. When all species present, i.e., $(\rhu,T)\in$, ket us re-wtnre Fxuation \[eqn:iconsegv\]-\[eqn:qion\] in log form while ising \[eqn:econscrv\], $$\bxgin{qligned}
\label{eqn:dlnndlix}
\begin{pmatrix}
2n_{\ch{H2}}  &  n_{\ch{H}}  &  n_{\eh{n+}}  \\
-1  &  2  &  0  \\
0  &  -1  &  2
\wnd{pmdtrif}
\begkb{pmxtrjx}
v\ln n_{\ch{H2}}\\
d\pn i_{\ch{H}}\\
d\ln n_{\ch{G+}}
\end{pmatrix}=
\vegin{pmatrix}
n_{\rm Htou}d\lg n_{\rm Htot}\\
d\ln s_{\rm dif}\\
d\jn q_{\
/2\\ d\ln Z_{\ch{H}}/d\ln T&=3/2+\phi_{\rm dis}/(2kT)\\ d\ln Z_{\ch{H+}}/d\ln ion})/(2kT)\\ Z_{\ch{e-}}/d\ln T&=3/2\\ Z_{\ch{H2}}/d\ln T^{2}&=0\\ d^{2}\ln T^{2}&=-(\phi_{\rm ion})/(2kT)\\ d^{2}\ln Z_{\ch{e-}}/d\ln This would not the case if one consider the and rotational degree of freedom of as in [@tomida2013; @tomida2015]. The only derivatives we need to be calculated are $\left(\frac{\partial\ln n_{i}}{\partial\ln T}\right)_{\rho}$ and $\left(\frac{\partial\ln n_{i}}{\partial\ln\rho}\right)_{T}$. When species i.e., let re-write Equation \[eqn:nconserv\]-\[eqn:qion\] in log form while using \[eqn:econserv\], $$\begin{aligned} \label{eqn:dlnndlnx} \begin{pmatrix} 2n_{\ch{H2}} & n_{\ch{H}} & \\ -1 & 2 & 0 \\ 0 -1 & 2 \end{pmatrix} d\ln n_{\ch{H2}}\\ d\ln n_{\ch{H}}\\ d\ln \end{pmatrix}= n_{\rm Htot}d\ln Htot}\\ q_{\rm d\ln q_{\
/2\\
    d\ln Z_{\ch{H}}/d\ln T&=3/2+\phi_{\rm dis}/(2kT)\\
    d\ln Z_{\Ch{H+}}/d\ln T&=3/2+(\phi_{\Rm dis}+2\Phi_{\Rm iOn})/(2KT)\\
    d\lN Z_{\ch{E-}}/d\ln T&=3/2\\
    d^{2}\ln Z_{\ch{H2}}/d\lN t^{2}&=0\\
    d^{2}\ln z_{\ch{H}}/d\ln T^{2}&=-\phi_{\rm dis}/(2kT)\\
    d^{2}\ln Z_{\Ch{H+}}/d\lN T^{2}&=-(\PHi_{\rm DIs}+2\Phi_{\rm Ion})/(2kT)\\
    d^{2}\lN z_{\cH{E-}}/D\ln t^{2}&=0\eNd{AliGnED}$$ THis woUld Not be thE case if one ConSiDer the vibratIOnAl and rotatIonAl degree of frEedOm of as In [@TomIDa2013; @tomIda2015].

the onLy deriVAtives That we neeD tO Be calcULated arE $\LEfT(\fraC{\partial\ln n_{i}}{\partiAL\lN t}\right)_{\rho}$ and $\leFt(\frac{\PaRTiAL\Ln n_{I}}{\paRtial\ln\rho}\RiGht)_{T}$. WHEn all spECiES PResENt, i.e., $(\rho,T)\in$, let Us re-write EqUAtiOn \[eqn:nCoNseRV\]-\[eqn:qiOn\] in lOg FOrm While using \[eQn:ecOnserv\], $$\begIn{aligNEd}
\label{EQn:dlnndLnx}
\begIn{pMatRix}
2n_{\CH{H2}}  &  N_{\cH{H}}  &  n_{\Ch{N+}}  \\
-1  &  2  &  0  \\
0  &  -1  &  2
\End{PMaTriX}
\BegIn{pmatriX}
d\Ln N_{\ch{H2}}\\
d\Ln n_{\cH{h}}\\
D\LN n_{\ch{h+}}
\enD{pmaTrix}=
\bEgin{pmatrix}
n_{\rM HtOt}d\lN N_{\rm htot}\\
d\Ln q_{\rm Dis}\\
d\Ln Q_{\
/2\\
    d\ln Z_{\ch{H}}/d \ln T&=3/2 +\phi _{\ rmdi s}/( 2kT) \\
    d\ln Z_ { \ch{ H+}}/d\ln T&=3/2+(\phi _{\rm  d i s}+2 \ ph i_{\r m ion}) / (2 k T )\\ 
  d\ ln Z_ {\ch{ e-} }/d\lnT&=3/2\\
  d^ {2}\ln Z_{\c h {H 2}}/d\ln T ^{2 }&=0\\
    d ^{2 }\ln Z _{ \ch { H}}/d \ln  T^{2 }&=-\p h i_{\rm  dis}/(2k T) \ \
   d ^{2}\ln Z _{ \ch{ H+}}/d\ln T^{2}&= - (\ p hi_{\rm dis}+2 \phi_{ \r m  i o n })/ (2k T)\\
    d ^{ 2}\ln Z_{\ch{ e -} } / d \ln T^{2}&=0\end{ aligned}$$T his  would  n otb e thecaseif one  consider t he v ibrationa l andr otation a l degre e of f ree dom  ofa sin  [@ to m ida 2 01 3;@ tom ida2015] .
T he on ly d e r i v ativ esthat  we n eed to be cal cul ated are  $\le ft(\f rac{ \p artia l\ln n _{i}} {\ partial\ln T}\r ight )_{\rho}$  an d$\l ef t(\fr a c{\par tia l\l n n_{i} }{\part i al\ ln \ r h o} \right)_{T}$. When  a l l  s pecies p resent ,  i .e . , $(\rho ,T )\i n$,l e t usre-w r it e Equati on \[e q n: nc onserv\ ]- \[eqn: qi on\ ] i n log form  while  using \ [eqn: e conserv\], $$\ b egin{aligned} \l a b el { eqn: dln ndlnx}
\beg in{p m atri x}
2 n _{ \ch { H2}} &  n _{ \ ch { H}}  &  n_{\ch{n+}}  \\
-1 &  2   &  0  \\
0 &  -1  &2 \end{pma trix } 
\ b egin{pmatrix}d\lnn_{\ch{H2} } \\
d\lnn_{\c h{H}}\\d\ln n_{\ c h {H+}}
\e nd{ pma tri x}= \ be gin{pmatrix}n _ {\rm  H tot}d\l n n _{\rm H tot }\\ 
d\ lnq_ {\rm dis} \\
d\lnq_ {\ 
/2\\
 _ _d\ln Z_{\ch{H}}/d\ln T&=3/2+\phi_{\rm dis}/(2kT)\\
_ __d\ln Z_{\ch{H+}}/d\ln_T&=3/2+(\phi_{\rm_dis}+2\phi_{\rm ion})/(2kT)\\
 _  d\ln_Z_{\ch{e-}}/d\ln T&=3/2\\
  _ d^{2}\ln Z_{\ch{H2}}/d\ln_T^{2}&=0\\
_   d^{2}\ln Z_{\ch{H}}/d\ln T^{2}&=-\phi_{\rm dis}/(2kT)\\
    d^{2}\ln Z_{\ch{H+}}/d\ln T^{2}&=-(\phi_{\rm dis}+2\phi_{\rm_ion})/(2kT)\\
_ __d^{2}\ln_Z_{\ch{e-}}/d\ln T^{2}&=0\end{aligned}$$ This would not_be the case if one_consider the_vibrational and rotational degree of freedom of as_in_[@tomida2013; @tomida2015].

The only_derivatives that we need to be calculated are $\left(\frac{\partial\ln_n_{i}}{\partial\ln T}\right)_{\rho}$ and $\left(\frac{\partial\ln n_{i}}{\partial\ln\rho}\right)_{T}$. When_all species present,_i.e.,_$(\rho,T)\in$,_let us re-write Equation_\[eqn:nconserv\]-\[eqn:qion\] in log form while using_\[eqn:econserv\], $$\begin{aligned}
\label{eqn:dlnndlnx}
\begin{pmatrix}
2n_{\ch{H2}}  &  n_{\ch{H}}_ &  n_{\ch{n+}}  \\
-1 _&  2  & _0  \\
0  &_ -1_ &  2
\end{pmatrix}
\begin{pmatrix}
d\ln n_{\ch{H2}}\\
d\ln_n_{\ch{H}}\\
d\ln n_{\ch{H+}}
\end{pmatrix}=
\begin{pmatrix}
n_{\rm Htot}d\ln_n_{\rm Htot}\\
d\ln_q_{\rm dis}\\
d\ln q_{\