 boundaries negligible, and determine the bound states similarly as in Section \[sec:com\]. To combine reasonable numerical efficiency with a good accuracy, we choose a segment length of $\Delta=2$$\mathrm{nm}$ (the last segment of each well or barrier is $\Delta\leq2$$\mathrm{nm}$).

Various models are available for including nonparabolicity [@1987PhRvB..35.7770N]; here, we use an energy dependent effective mass $m_{j}^{\ast}\left(  E\right)  =m_{j}^{\ast}\left[  1+\left(  E-V_{j}\right)
/E_{g,j}\right]  $ (with band gap energy $E_{g,j}$ at position $z_{j}$), which can straightforwardly be implemented into the transfer matrices. The Poisson equation is given by [@2000JAP....87.7931C; @2008SeScT..23l5040L]$$-\partial_{z}\left[  \epsilon\left(  z\right)  \partial_{z}\varphi\left(  z\right)
\right]  =e\left[  N\left(  z\right)  -\sum_{n}n_{\mathrm{2D},n}\left|
\psi_{n}\left(  z\right)  \right|  ^{2}\right] ,\label{poisson}%$$ leading to an additional potential -e$\varphi$ in (\[effm\]). Here, $\epsilon\left(  z\right)  $ is the permittivity, $e$ is the elementary charge, $N\left(  z\right)  $ is the doping concentration, and $n_{\mathrm{2D}%
,n}$ is the electron sheet density of level $n$ with wave function $\psi
_{n}\left(  z\right)  $. While for an operating QCL, $n_{\mathrm{2D},n}$ can only be exactly determined by detailed carrier transport simulations [@2007JAP...101h6109J], this is prohibitive for design optimizations of experimental QCL structures over an extended parameter range. Thus, for solving the Schrödinger-Poisson system, simpler and much faster approaches are commonly adopted, such as applying Fermi-Dirac statistics [@2000JAP....87.7931C; @2008SeScT..23l50
boundaries negligible, and determine the bound states similarly as in Section \[sec: com\ ]. To compound fair numerical efficiency with a good accuracy, we choose a section length of $ \Delta=2$$\mathrm{nm}$ (the final segment of each well or barrier is $ \Delta\leq2$$\mathrm{nm}$). 

 assorted models are available for including nonparabolicity [ @1987PhRvB.. 35.7770N ]; here, we practice an department of energy dependent effective mass $ m_{j}^{\ast}\left (  E\right)   = m_{j}^{\ast}\left [   1+\left (  east - V_{j}\right) 
 /E_{g, j}\right ]   $ (with band gap department of energy $ E_{g, j}$ at position $ z_{j}$), which can straightforwardly be follow through into the transfer matrices. The Poisson equality is given by [ @2000JAP.... 87.7931C; @2008SeScT.. 23l5040L]$$-\partial_{z}\left [   \epsilon\left (  z\right)   \partial_{z}\varphi\left (  z\right) 
 \right ]   = e\left [   N\left (  z\right)   -\sum_{n}n_{\mathrm{2D},n}\left| 
 \psi_{n}\left (  z\right)   \right|   ^{2}\right ], \label{poisson}%$$ leading to an additional likely -e$\varphi$ in (\[effm\ ]). Here, $ \epsilon\left (  z\right)   $ is the permittivity, $ e$ is the elementary charge, $ N\left (  z\right)   $ is the doping concentration, and $ n_{\mathrm{2D}% 
, n}$ is the electron sheet density of degree $ n$ with wave function $ \psi 
 _ { n}\left (  z\right)   $. While for an operating QCL, $ n_{\mathrm{2D},n}$ can only be exactly determined by detailed mailman transport simulations [ @2007JAP... 101h6109J ], this is prohibitive for design optimizations of experimental QCL structures over an extended parameter range. Thus, for solving the Schrödinger - Poisson system, simpler and much faster approaches are commonly adopted, such as apply Fermi - Dirac statistics [ @2000JAP.... 87.7931C; @2008SeScT.. 23l50
 bokndaries negligible, and aetermine the bound stetes sijilarly xs in Section \[sec:com\]. To combmne eeasobable numerical efficidncy with a good qccucacy, we choose a segment lengtg of $\Beota=2$$\mathrm{nm}$ (thg last segmett of each weln ur barrier is $\Delta\leq2$$\mathrm{nm}$).

Various models age available fjr imsludjng nonparabolicity [@1987PhRvB..35.7770N]; here, we use an energy depencent effective mass $m_{j}^{\ast}\lfft(  E\right)  =m_{j}^{\ast}\left[  1+\left(  E-V_{j}\tjghe)
/W_{g,j}\right]  $ (wkth band gap energy $E_{g,n}$ at position $z_{j}$), which can strakghtfprwardly bg jmonemented invo the transfer mabgices. Tve Poisxon equation ix gmven by [@2000JAP....87.7931C; @2008SeScT..23l5040L]$$-\partiel_{z}\left[  \epsilon\left(  z\right)  \pcrtial_{z}\varphi\left(  z\rugyt)
\rigvt]  =a\lefg[  N\uefu(  v\rifht)  -\skm_{n}i_{\mathrm{2D},n}\lert|
\psi_{n}\left(  z\right)  \right|  ^{2}\rignt] ,\oabel{poisson}%$$ meadind eo an additional potential -e$\varphi$ in (\[exfm\]). Here, $\epsilon\left(  z\rigyt)  $ is the permittivlty, $e$ is ehe elementary charge, $N\left(  z\right)  $ is the dopitg coicdntxqtion, qnf $n_{\mathrm{2D}%
,n}$ is the electron sheet density of jsvtl $k$ with wave funcbion $\psi
_{n}\left(  z\ribhh)  $. While for an operacjnf QCL, $n_{\mathrm{2D},n}$ caj only fe exqctly dettrminrd by detailed carrier tranwport simulanionw [@2007JAP...101h6109J], this is prlhibitive fur dgsign pptimizations of experioentzl QCL struftures ovsf an extended pafamvter range. Thus, for solving thq Schrödinjer-Ponsson syrtem, simplqr and mucj fasbar approaches are fommouly agopted, sucj as applying Fermi-Dirac statisvmcs [@2000JAP....87.7931C; @2008SeScJ..23l50
boundaries negligible, and determine the bound states in \[sec:com\]. To reasonable numerical efficiency choose segment length of (the last segment each well or barrier is $\Delta\leq2$$\mathrm{nm}$). models are available for including nonparabolicity [@1987PhRvB..35.7770N]; here, we use an energy dependent mass $m_{j}^{\ast}\left( E\right) =m_{j}^{\ast}\left[ 1+\left( E-V_{j}\right) /E_{g,j}\right] $ (with band gap energy $E_{g,j}$ position which straightforwardly implemented into the transfer matrices. The Poisson equation is given by [@2000JAP....87.7931C; @2008SeScT..23l5040L]$$-\partial_{z}\left[ \epsilon\left( z\right) \partial_{z}\varphi\left( \right] =e\left[ N\left( z\right) -\sum_{n}n_{\mathrm{2D},n}\left| \psi_{n}\left( z\right) \right| ,\label{poisson}%$$ leading to an potential -e$\varphi$ in (\[effm\]). Here, z\right) is the $e$ the charge, $N\left( z\right) is the doping concentration, and $n_{\mathrm{2D}% ,n}$ is the electron sheet density of level $n$ with wave $\psi _{n}\left( While for operating $n_{\mathrm{2D},n}$ only be exactly detailed carrier transport simulations [@2007JAP...101h6109J], this design optimizations of experimental QCL structures over an parameter range. for solving the Schrödinger-Poisson system, simpler much faster approaches are commonly adopted, such as Fermi-Dirac statistics [@2000JAP....87.7931C; @2008SeScT..23l50
 boundaries negligible, and deTermine the Bound StaTes SiMilaRly aS in Section \[sec:cOM\]. To cOmbine reasonable numeriCal efFiCIencY WiTh a goOd accurACy, WE ChoOsE a SegMeNT lEngth Of $\DElta=2$$\matHrm{nm}$ (the laSt sEgMent of each weLL oR barrier is $\delTa\leq2$$\mathrm{nM}$).

VaRious mOdEls ARe avaIlaBle foR incluDIng nonParabolicItY [@1987phRvB..35.7770N]; HEre, we usE AN eNergY dependent effectiVE mASs $m_{j}^{\ast}\left(  E\riGht)  =m_{j}^{\aSt}\LEfT[  1+\LEft(  e-V_{j}\Right)
/E_{g,j}\riGhT]  $ (with BAnd gap eNErGY $e_{G,j}$ aT Position $z_{j}$), whiCh can straigHTfoRwardlY bE imPLementEd intO tHE trAnsfer matriCes. THe Poisson EquatiON is giveN By [@2000JAP....87.7931C; @2008SEScT..23l5040L]$$-\ParTiaL_{z}\leFT[  \ePsIloN\lEFt(  z\RIgHt)  \pARtiAl_{z}\varphI\lEfT(  z\rigHt)
\riGHT]  =E\Left[  n\leFt(  z\rIght)  -\sUm_{n}n_{\mathrm{2D},n}\lEft|
\Psi_{n}\LEft(  Z\righT)  \righT|  ^{2}\rigHt] ,\Label{PoissoN}%$$ leadInG to an additional PoteNtial -e$\varPhi$ In (\[EffM\]). HEre, $\epSIlon\leFt(  z\RigHt)  $ is the PermittIVitY, $e$ IS THe Elementary charge, $N\lEfT(  Z\RiGht)  $ is the Doping COnCeNTration, aNd $N_{\maThrm{2d}%
,N}$ Is the ElecTRoN sheet deNsity oF LeVeL $n$ with wAvE functIoN $\psI
_{n}\lEft(  z\rIGht)  $. WHile foR an operaTing Qcl, $n_{\mathrm{2D},n}$ can oNLy be exactly deTErMINeD By deTaiLed carrier tRansPOrt sImulATiOns [@2007jaP...101h6109J], tHis is PrOHiBItive for design optimIzAtions Of expErimental QCL sTructures oVER An extendEd paRAmETer range. Thus, foR solvIng the SchrÖDinger-PoIsson System, siMpler and mUCH faster aPprOacHes Are COMmOnly adopted, suCH As apPlYing FerMi-DIrac staTisTicS [@2000JAp....87.7931C; @2008SESCT..23l50
 boundaries negligible, an d determin e the  bo und  s tate s si milarly as inS ecti on \[sec:com\]. To com binere a sona b le  nume rical e f fi c i enc ywi thag oo d acc ura cy, wechoose a s egm en t length of$ \D elta=2$$\m ath rm{nm}$ (the  la st seg me nto f eac h w ell o r barr i er is$\Delta\l eq 2 $$\mat h rm{nm}$ ) . 
Vari ous models are av a il a ble for includ ing no np a ra b o lic ity  [@1987PhR vB ..35. 7 770N];h er e , weu se an energydependent e f fec tive m as s $ m _{j}^{ \ast} \l e ft(   E\right) =m_ {j}^{\ast }\left [   1+\le f t(  E-V _{j}\r igh t)/E_{ g ,j }\ rig ht ]   $ (w ith ban d gap en er gy  $E_{ g,j} $ a t  pos iti on $ z_{j} $), which can  st raig h tfo rward ly be  imp le mente d into  thetr ansfer matrices . Th e Poisson  eq ua tio nis gi v en by[@2 000 JAP.... 87.7931 C ; @ 20 0 8 S eS cT..23l5040L]$$-\p ar t i al _{z}\lef t[  \e p si lo n \left( z\ rig ht)\ parti al_{ z }\ varphi\l eft( z \r ig ht)
\ri gh t]  =e \l eft [ N\lef t (  z \right )  -\sum _{n}n _ {\mathrm{2D},n } \left|
\psi_{ n }\ l e ft (   z\ rig ht)  \right |  ^ { 2}\r ight ]  , \la b el{po isson }% $ $l eading to an additi on al pot entia l -e$\varphi$  in (\[eff m \ ] ). Here,  $\e p si l on\left(  z\ri ght) $ is thep ermittiv ity,$e$ is t he elemen t a ry charg e,$N\ lef t(z \r ight)  $ is t h e  dop in g conce ntr ation,and  $n _{\ mat hr m{2D}%
,n }$ is th eel ec tr onsheet densityof  le ve l $ n$ wi t h wave  func tion  $ \p s i
_ {n}\lef t (z \rig ht ) $.Whi le  foran o p era ting QC L, $n_{\m ath r m{2D }, n} $ can o nly be exactl ydetermined  b y d etaile d carriertransport simulations [ @ 2007JAP ... 101h6 109J ], this i s p rohibi tiv e  for d esignoptim iz ati o n s ofe x pe rim en tal QCL st r u ctu res o ve r an  extend ed parameter range .  Th us, for solvi ngtheS c hr ödi n ge r -Po is s ons y stem, simpler a nd much fa st e rapproaches are  c ommonly  adopte d, su c h as ap plying Fe rmi-Dirac  s tati s t ics  [@2000JAP ....87.7 931C; @20 0 8SeSc T .. 23l50 
 boundaries_negligible, and_determine the bound states_similarly as_in_Section \[sec:com\]._To_combine reasonable numerical_efficiency with a_good accuracy, we choose_a segment length_of_$\Delta=2$$\mathrm{nm}$ (the last segment of each well or barrier is $\Delta\leq2$$\mathrm{nm}$).

Various models are available_for_including nonparabolicity_[@1987PhRvB..35.7770N];_here,_we use an energy dependent_effective mass $m_{j}^{\ast}\left(  E\right)_ =m_{j}^{\ast}\left[_ 1+\left(  E-V_{j}\right)
/E_{g,j}\right]  $ (with band_gap_energy $E_{g,j}$ at_position $z_{j}$), which can straightforwardly be implemented into the_transfer matrices. The Poisson equation is_given by [@2000JAP....87.7931C;_@2008SeScT..23l5040L]$$-\partial_{z}\left[__\epsilon\left(  z\right) _\partial_{z}\varphi\left(  z\right)
\right]  =e\left[ _N\left(  z\right)  -\sum_{n}n_{\mathrm{2D},n}\left|
\psi_{n}\left( _z\right)  \right|  ^{2}\right] ,\label{poisson}%$$ leading_to an additional potential -e$\varphi$ in_(\[effm\]). Here, $\epsilon\left(  z\right)_ $_is the permittivity, $e$ is_the elementary charge,_$N\left( _z\right)  $_is the doping concentration, and $n_{\mathrm{2D}%
,n}$_is the electron_sheet density of level $n$ with_wave_function $\psi
_{n}\left( _z\right)__$. While_for an operating_QCL,_$n_{\mathrm{2D},n}$ can_only_be exactly determined by detailed carrier_transport_simulations [@2007JAP...101h6109J], this is prohibitive for design_optimizations of experimental QCL_structures_over an extended parameter_range. Thus, for solving the_Schrödinger-Poisson system, simpler and much faster_approaches are_commonly adopted,_such as applying Fermi-Dirac statistics [@2000JAP....87.7931C; @2008SeScT..23l50