T_{b}(v) \simeq \tau_{\rm tot} \cdot T_{\rm ex}
\label{eqn_Tmb_thin}\end{aligned}$$ indicating that the $T_{\rm ex}$ and $\tau_{\rm tot}$ can be determined arbitrarily. Therefore, we conduct non-LTE analysis in Section \[subsec::non-lte\] to derive the physical conditions.

![image](non_lte_model_bary.pdf){width="\linewidth"}

$N_2H^+$ RADEX Non-LTE Modeling {#subsec::non-lte}
-------------------------------

Using RADEX [@radex], a non-LTE radiative transfer code, we construct spectra models in $\rm N_2H^+$ 3–2 and 1–0. The synthetic spectra are constructed with the equation $$\begin{aligned}
&&T_{b}(v) = \Psi \left( \frac{\sum J(T_{\rm ex}^i)\,\tau_i(v)}{\sum \tau_i(v)} - J(T_{\rm bg}) \right) \left( 1 - e^{-\sum \tau_i(v)} \right)
\nonumber \\
\label{Tb_non-lte}\end{aligned}$$ where $T_{\rm ex}^i$ and $\tau_i(v)$ represent the excitation temperature and optical depth for all hyperfine components in the 3–2 or 1–0 transitions, and $\Psi$ is the beam filling factor. In our models, the beam filling factors are assumed to be unity. We also construct an intensity ratio model by dividing the integrated spectra model in $\rm N_2H^+$ 3–2 with that in $\rm N_2H^+$ 1–0.

The constructed models can be represented by a three-dimensional grid, where the three axes are $\rm H_2$ density ($n_{\rm H_2}$) ranging from $10^4$ to $10^9$ $\rm cm^{-3}$, kinetic temperature ($T_{\rm kin}$) ranging from $8$ to $60$ K, and the ratio of $\rm N_2H^+$ column density to linewidth ($N(\rm{N_2H^+})$/$\Delta v$) ranging from $5\times10^{7}$ to $5\times10^{8}$ $\rm s\ cm^{-3}$. The step sizes of the grid are 1
T_{b}(v) \simeq \tau_{\rm tot } \cdot T_{\rm ex } 
 \label{eqn_Tmb_thin}\end{aligned}$$ indicating that the $ T_{\rm ex}$ and $ \tau_{\rm tot}$ can be determined arbitrarily. consequently, we impart non - LTE analysis in Section \[subsec::non - lte\ ] to derive the forcible conditions. 

! [ image](non_lte_model_bary.pdf){width="\linewidth " } 

 $ N_2H^+$ RADEX Non - LTE Modeling { # subsec::non - lte } 
 ------------------------------- 

 use RADEX [ @radex ], a non - LTE radiative transfer code, we reconstruct spectra model in $ \rm N_2H^+$ 3–2 and 1–0. The man-made spectra are constructed with the equality $ $ \begin{aligned } 
 & & T_{b}(v) = \Psi \left (\frac{\sum J(T_{\rm ex}^i)\,\tau_i(v)}{\sum \tau_i(v) } - J(T_{\rm bg }) \right) \left (1 - e^{-\sum \tau_i(v) } \right) 
 \nonumber \\ 
 \label{Tb_non - lte}\end{aligned}$$ where $ T_{\rm ex}^i$ and $ \tau_i(v)$ represent the excitation temperature and optical depth for all hyperfine component in the 3–2 or 1–0 transitions, and $ \Psi$ is the beam filling factor. In our models, the radio beam filling factors are assumed to be oneness. We also construct an intensity proportion model by separate the integrated spectra model in $ \rm N_2H^+$ 3–2 with that in $ \rm N_2H^+$ 1–0. 

 The constructed model can be represented by a three - dimensional grid, where the three axes are $ \rm H_2 $ density ($ n_{\rm H_2}$) ranging from $ 10 ^ 4 $ to $ 10 ^ 9 $ $ \rm cm^{-3}$, kinetic temperature ($ T_{\rm kin}$) range from $ 8 $ to $ 60 $ K, and the ratio of $ \rm N_2H^+$ column density to linewidth ($ N(\rm{N_2H^+})$/$\Delta v$) ranging from $ 5\times10^{7}$ to $ 5\times10^{8}$ $ \rm s\ cm^{-3}$. The step sizes of the grid are 1
T_{b}(v) \simeq \tau_{\rm tot} \cdot T_{\ro ex}
\label{eqn_Tmb_jhun}\end{anigned}$$ indicatkng that the $T_{\rm ex}$ and $\tau_{\rl rot}$ cqn be determined arbitfarily. Thvrefore, ww coiduct non-LTE analysis in Sectikk \[subvxc::non-lte\] to derlve the phyvical conditiots.

![kmcge](non_lte_model_bary.pdf){width="\linewidth"}

$N_2R^+$ RADEX Nln-LTE Modeling {#subxqc::noh-lte}
-------------------------------

Using RADEX [@radex], a non-LTE radjative uransfer code, we cpnstruct spectra models in $\rm J_2H^+$ 3–2 and 1–0. The synthftic spectrq arq constructed with the tqbation $$\begih{aligned}
&&T_{b}(v) = \Psi \left( \frac{\sum J(G_{\rm er}^i)\,\tau_i(v)}{\sum \jcy_i(v)} - J(T_{\rm bg}) \rijht) \lest( 1 - e^{-\sum \tam_p(v)} \righd)
\nonumbrr \\
\label{Tb_non-lbe}\end{eligbed}$$ where $T_{\rm ex}^i$ and $\tau_i(v)$ represent tre excitadikn temperature ane ipticdl dapth dor alm iypsrfine cojponents ih the 3–2 or 1–0 rransitions, and $\Psi$ if the beam filming fwceor. In our models, the beam filling factmrs are assumed to be unitt. We also construct aj intensiey ratio model by dividing the integrated spectra modem in $\ri B_2H^+$ 3–2 with that in $\rm N_2H^+$ 1–0.

The constructed models cwh ne represented bn a three-dimensionsl gtyd, where the jhree arss are $\rm H_2$ density ($j_{\rm H_2}$) rwngint from $10^4$ tj $10^9$ $\rk cm^{-3}$, kinetic temperature ($T_{\rn kin}$) ranginy feom $8$ to $60$ K, and the ratio of $\ro N_2H^+$ colukn density to linewidth ($N(\rm{H_2H^+})$/$\Delta v$) rwnging frko $5\times10^{7}$ to $5\times10^{8}$ $\rm s\ cm^{-3}$. The step sizes of the drid are 1
T_{b}(v) \simeq \tau_{\rm tot} \cdot T_{\rm ex} that $T_{\rm ex}$ $\tau_{\rm tot}$ can conduct analysis in Section to derive the conditions. ![image](non_lte_model_bary.pdf){width="\linewidth"} $N_2H^+$ RADEX Non-LTE Modeling ------------------------------- Using RADEX [@radex], a non-LTE radiative transfer code, we construct spectra models $\rm N_2H^+$ 3–2 and 1–0. The synthetic spectra are constructed with the equation &&T_{b}(v) \Psi \frac{\sum ex}^i)\,\tau_i(v)}{\sum \tau_i(v)} - J(T_{\rm bg}) \right) \left( 1 - e^{-\sum \tau_i(v)} \right) \nonumber \\ \label{Tb_non-lte}\end{aligned}$$ where ex}^i$ and $\tau_i(v)$ represent the excitation temperature and depth for all hyperfine in the 3–2 or 1–0 and is the filling In models, the beam factors are assumed to be unity. We also construct an intensity ratio model by dividing the integrated model in 3–2 with in N_2H^+$ The constructed models represented by a three-dimensional grid, where are $\rm H_2$ density ($n_{\rm H_2}$) ranging from to $10^9$ cm^{-3}$, kinetic temperature ($T_{\rm kin}$) ranging $8$ to $60$ K, and the ratio of N_2H^+$ column density to linewidth ($N(\rm{N_2H^+})$/$\Delta v$) ranging from $5\times10^{7}$ to $5\times10^{8}$ $\rm s\ cm^{-3}$. sizes of the grid 1
T_{b}(v) \simeq \tau_{\rm tot} \cdot T_{\rm ex}
\Label{eqn_TmB_thin}\End{AliGnEd}$$ inDicaTing that the $T_{\rm EX}$ and $\Tau_{\rm tot}$ can be determineD arbiTrARily. tHeReforE, we condUCt NON-LTe aNaLysIs IN SEctioN \[suBsec::non-Lte\] to derivE thE pHysical condiTIoNs.

![image](non_Lte_Model_bary.pdf){WidTh="\lineWiDth"}

$n_2h^+$ RADEx NoN-LTE MOdelinG {#Subsec::Non-lte}
-------------------------------

UsiNg raDEX [@raDEx], a non-Lte RaDiatIve transfer code, we COnSTruct spectra moDels in $\Rm n_2h^+$ 3–2 aND 1–0. the SynThetic specTrA are cONstructED wITH The EQuation $$\begin{aLigned}
&&T_{b}(v) = \PsI \LefT( \frac{\sUm j(T_{\rM Ex}^i)\,\tau_I(v)}{\sum \TaU_I(v)} - J(t_{\rm bg}) \right) \lEft( 1 - e^{-\Sum \tau_i(v)} \rIght)
\noNUmber \\
\laBEl{Tb_non-Lte}\end{AliGneD}$$ wheRE $T_{\Rm Ex}^i$ AnD $\Tau_I(V)$ rEprESenT the exciTaTiOn temPeraTURE And oPtiCal dEpth fOr all hyperfinE coMponENts In the 3–2 Or 1–0 traNsitIoNs, and $\psi$ is tHe beaM fIlling factor. In oUr moDels, the beAm fIlLinG fActorS Are assUmeD to Be unity. we also cONstRuCT AN iNtensity ratio model By DIViDing the iNtegraTEd SpECtra modeL iN $\rm n_2H^+$ 3–2 wiTH That iN $\rm N_2h^+$ 1–0.

thE construCted moDElS cAn be repReSented By A thRee-DimenSIonaL grid, wHere the tHree aXEs are $\rm H_2$ densitY ($N_{\rm H_2}$) ranging frOM $10^4$ tO $10^9$ $\RM cM^{-3}$, KineTic Temperature ($t_{\rm kIN}$) ranGing FRoM $8$ to $60$ k, And thE ratiO oF $\Rm n_2h^+$ column density to linEwIdth ($N(\rM{N_2H^+})$/$\DeLta v$) ranging frOm $5\times10^{7}$ to $5\tIMES10^{8}$ $\rm s\ cm^{-3}$. ThE steP SiZEs of the grid are 1
T_{b}(v) \simeq \tau_{\rmtot} \cdot  T_{\ rmex} 
\ labe l{eq n_Tmb_thin}\en d {ali gned}$$ indicating tha t the  $ T _{\r m  e x}$ a nd $\ta u _{ \ r m t ot }$  ca nb edeter min ed arbi trarily. T her ef ore, we cond u ct  non-LTE a nal ysis in Sect ion  \[sub se c:: n on-lt e\]  to d erivet he phy sical con di t ions.![image ] ( no n_lt e_model_bary.pdf) { wi d th="\linewidth "}

$N _2 H ^+ $ RAD EXNon-LTE Mo de ling{ #subsec : :n o n - lte } 
------------ ----------- - --- ----
Us ing RADEX[@rad ex ] , a  non-LTE ra diat ive trans fer co d e, we c o nstruct  spect ramod elsi n$\ rmN_ 2 H^+ $  3 –2a nd1–0. The  s yn theti c sp e c t r a ar e c onst ructe d with the eq uat ion$ $\b egin{ align ed}&& T_{b} (v) =\Psi\l eft( \frac{\sum  J(T _{\rm ex} ^i) \, \ta u_ i(v)} { \sum \ tau _i( v)} - J (T_{\rm bg} )\ r i gh t) \left( 1 - e^{- \s u m  \ tau_i(v) } \rig h t) 
\ n onumber\\ 
\l abel { T b_non -lte } \e nd{align ed}$$w he re  $T_{\r mex}^i$  a nd$\t au_i( v )$ r eprese nt the e xcita t ion temperatur e  and opticald ep t h  f o r al l h yperfine co mpon e ntsin t h e3–2 or 1– 0 tra ns i ti o ns, and $\Psi$ is t he  beamfilli ng factor. In  our model s , the beam  fil l in g  factors are a ssume d to be un i ty. We a lso c onstruct  an inten s i ty ratio  mo del  by  di v i di ng the integr a t ed s pe ctra mo del  in $\r m N _2H ^+$  3– 2with that  in $\rm  N _2 H^ +$  1– 0.

T h e constr uc ted  m ode ls ca n  be re prese nted  b ya  th ree-dim e ns i o nalgr id , wh ere  t he th reea xes  are $\ rm H_2$ d ens i ty ( $n _{ \rm H_2 }$) ranging f ro m $10^4$ t o$10 ^9$ $\ r m  cm^{-3} $, kinetic temperature( $T_{\rm  ki n}$)rang ing from$8$  to $6 0$K , andthe ra tio o f$\r m N_2H^ + $  c olu mn  density t o lin ewidt h($N( \rm{N_2 H^+})$/$\Delta v$) ran ging from $5\ tim es10 ^ { 7} $ t o  $ 5 \ti me s 10^ { 8 }$ $\rm s\ cm^{ -3}$. Thest e psizes of t h e g ri d are 1 
T_{b}(v) \simeq_\tau_{\rm tot}_\cdot T_{\rm ex}
\label{eqn_Tmb_thin}\end{aligned}$$ indicating_that the_$T_{\rm_ex}$ and_$\tau_{\rm_tot}$ can be_determined arbitrarily. Therefore,_we conduct non-LTE analysis_in Section \[subsec::non-lte\]_to_derive the physical conditions.

![image](non_lte_model_bary.pdf){width="\linewidth"}

$N_2H^+$ RADEX Non-LTE Modeling {#subsec::non-lte}
-------------------------------

Using RADEX [@radex], a non-LTE radiative transfer_code,_we construct_spectra_models_in $\rm N_2H^+$ 3–2 and_1–0. The synthetic spectra are_constructed with_the equation $$\begin{aligned}
&&T_{b}(v) = \Psi \left( \frac{\sum J(T_{\rm_ex}^i)\,\tau_i(v)}{\sum_\tau_i(v)} - J(T_{\rm_bg}) \right) \left( 1 - e^{-\sum \tau_i(v)} \right)
\nonumber \\
\label{Tb_non-lte}\end{aligned}$$_where $T_{\rm ex}^i$ and $\tau_i(v)$ represent_the excitation temperature_and_optical_depth for all hyperfine_components in the 3–2 or 1–0_transitions, and $\Psi$ is the beam_filling factor. In our models, the beam_filling factors are assumed to be_unity. We also construct an_intensity ratio_model by dividing the integrated_spectra model in_$\rm N_2H^+$_3–2 with that_in $\rm N_2H^+$ 1–0.

The constructed models_can be represented_by a three-dimensional grid, where the_three_axes are $\rm_H_2$_density_($n_{\rm H_2}$)_ranging from $10^4$_to_$10^9$ $\rm_cm^{-3}$,_kinetic temperature ($T_{\rm kin}$) ranging from_$8$_to $60$ K, and the ratio of_$\rm N_2H^+$ column density_to_linewidth ($N(\rm{N_2H^+})$/$\Delta v$) ranging_from $5\times10^{7}$ to $5\times10^{8}$ $\rm_s\ cm^{-3}$. The step sizes of_the grid_are 1