 (for example [@Kimura04]), the deformation is usually constrained by using the rotational invariant value $D\equiv Tr(QQ)/Tr^2(Q)$, where the matrix $Q$ is calculated by quadrupole operators as $Q_{\sigma\rho}=\langle \sum_i \hat\sigma_i
\hat\rho_i \rangle$ ($\hat\sigma=\hat{x},\hat{y},\hat{z}$ and $\hat\rho=\hat{x},\hat{y},\hat{z}$) [@Dote97]. Here $D$ is approximately related to the quadrupole deformation parameter $\beta$ as $D(\beta)= (5\beta^2/2\pi+1)/3$. In the present work, we used the modified quadrupole matrix $Q'_{\sigma\rho}\equiv Q_{\sigma\rho}-A\delta_{\sigma\rho}$ ($A$ is the mass number) instead of the original $Q_{\sigma\rho}$ and imposed the constraint on the $D'\equiv Tr(Q'Q')/Tr^2(Q')$. This is useful for He isotopes to obtain basis wave functions with various configurations on mesh points of the two-dimensional parameters, $\beta$ and $\langle \hat N^{\rm ho} \rangle$. The energy variation with the constraint values $N_{\rm const}$ and $\beta_{\rm const}$ was performed with respect to the parity-projected AMD wave function by minimizing the energy defined as, $$E\equiv \langle \hat{H}\rangle+
V^{N}(N_{\rm const}-\langle \hat N^{\rm ho} \rangle)^2
+V^{\beta}(D(\beta_{\rm const})-D')^2.\label{eq:const}$$ Here the artificial potentials are introduced to satisfy the condition of the constraints. With a given set of constraint values $(N_{\rm const},\beta_{\rm const})$ the optimum wave function $\Phi^\pm_{\rm AMD}(N_{\rm const},\beta_{\rm const})$ was obtained. Finally, we superposed the spin-parity eigen states projected from the obtained wave functions, $$|^8{\rm He}(J^\pm_n)\rangle = 
\sum_{N_{\rm const},\beta_{\rm const}}c^{J\pm}_n(N_{\rm const},\beta_{\rm const})
|P^{J}_{MK}\Phi^\pm_{\rm AMD}(N_{\rm const
(for example [ @Kimura04 ]), the deformation is usually constrained by use the rotational changeless value $ D\equiv Tr(QQ)/Tr^2(Q)$, where the matrix $ Q$ is calculated by quadrupole operators as $ Q_{\sigma\rho}=\langle \sum_i \hat\sigma_i 
 \hat\rho_i \rangle$ ($ \hat\sigma=\hat{x},\hat{y},\hat{z}$ and $ \hat\rho=\hat{x},\hat{y},\hat{z}$) [ @Dote97 ]. Here $ D$ is approximately related to the quadrupole distortion argument $ \beta$ as $ D(\beta)= (5\beta^2/2\pi+1)/3$. In the present work, we used the modified quadrupole matrix $ Q'_{\sigma\rho}\equiv Q_{\sigma\rho}-A\delta_{\sigma\rho}$ ($ A$ is the batch number) instead of the original $ Q_{\sigma\rho}$ and enforce the constraint on the $ D'\equiv Tr(Q'Q')/Tr^2(Q')$. This is useful for He isotopes to prevail basis wave functions with respective configuration on mesh points of the two - dimensional parameters, $ \beta$ and $ \langle \hat N^{\rm ho } \rangle$. The energy variation with the restraint values $ N_{\rm const}$ and $ \beta_{\rm const}$ was performed with respect to the parity bit - projected AMD wave function by minimizing the energy defined as, $ $ E\equiv \langle \hat{H}\rangle+ 
 V^{N}(N_{\rm const}-\langle \hat N^{\rm ho } \rangle)^2 
 + V^{\beta}(D(\beta_{\rm const})-D')^2.\label{eq: const}$$ Here the artificial potential are introduced to satisfy the condition of the constraints. With a given hardening of constraint values $ (N_{\rm const},\beta_{\rm const})$ the optimum wave function $ \Phi^\pm_{\rm AMD}(N_{\rm const},\beta_{\rm const})$ was obtained. Finally, we superposed the spin - parity eigen states projected from the obtained wave functions, $ $ |^8{\rm He}(J^\pm_n)\rangle = 
 \sum_{N_{\rm const},\beta_{\rm const}}c^{J\pm}_n(N_{\rm const},\beta_{\rm const }) 
 |P^{J}_{MK}\Phi^\pm_{\rm AMD}(N_{\rm const
 (fog example [@Kimura04]), the defurmation is usucoly coistrainsd by uskng the rotational invariant vqlue $E\equiv Tr(QQ)/Tr^2(Q)$, where tfe matrix $Q$ is caoculeted by quadrupole operabjrs zd $Q_{\sngna\rho}=\langle \suk_i \hat\sigmd_i
\hat\rho_i \rangne$ ($\hct\sigma=\hat{x},\hat{y},\hat{z}$ and $\hat\rho=\hat{x},\hae{y},\hat{z}$) [@Cohe97]. Here $D$ is akproxpmwtelg related to the quadrupole deformztion perameter $\beta$ ax $D(\beta)= (5\beta^2/2\pi+1)/3$. In the presfnt aork, we used the mldified quaerupjoe matrix $Q'_{\skgma\rho}\equpr Q_{\sigma\rho}-Z\delta_{\sigma\rho}$ ($A$ is the mass nuober) nnstead of jkw oghginal $Q_{\signa\rho}$ and imposed the convtraint on the $D'\equiv Tr(X'Q')/Te^2(Q')$. This is useful for He isotopes to objain basis wcve functions with vaeiius cmnficuraguonr oh keah poijts of the twk-dimensionao parameters, $\beta$ amd $\oangle \hat N^{\rj ho} \rwndle$. The energy variation with the constgainf values $N_{\rm const}$ and $\veta_{\rm const}$ was perflrmed witr respect to the parity-projected AMD wave functiot by jkniniding rhf energy defined as, $$E\equiv \langle \hat{H}\rangle+
V^{G}(H_{\rk bonst}-\langle \hat N^{\vm ho} \rangle)^2
+V^{\beta}(D(\neha_{\ti const})-D')^2.\label{gq:const}$$ Hede the artificial ootentiwls aee introdtced to satisfy the condition od the constrcinrs. With a given sec of constrannt vakues $(M_{\rm const},\beta_{\rm const})$ tke optjmum wave fknction $\Pgk^\pm_{\rm AMD}(N_{\rm conrt},\bvta_{\rk const})$ was obtained. Finaljy, we supxrposzd the soin-psrity qigen statfs projected from the obhaineb wava functiond, $$|^8{\rm He}(J^\pm_n)\rangle = 
\sum_{N_{\rm const},\uxta_{\rm const}}c^{J\km}_n(T_{\rm const},\beca_{\rm cpnst})
|P^{J}_{MK}\Phi^\pm_{\wm AMD}(N_{\rm consj
(for example [@Kimura04]), the deformation is usually using rotational invariant $D\equiv Tr(QQ)/Tr^2(Q)$, where by operators as $Q_{\sigma\rho}=\langle \hat\sigma_i \hat\rho_i \rangle$ and $\hat\rho=\hat{x},\hat{y},\hat{z}$) [@Dote97]. Here $D$ is related to the quadrupole deformation parameter $\beta$ as $D(\beta)= (5\beta^2/2\pi+1)/3$. In the present we used the modified quadrupole matrix $Q'_{\sigma\rho}\equiv Q_{\sigma\rho}-A\delta_{\sigma\rho}$ ($A$ is the mass number) of original and the constraint on the $D'\equiv Tr(Q'Q')/Tr^2(Q')$. This is useful for He isotopes to obtain basis wave with various configurations on mesh points of the parameters, $\beta$ and $\langle N^{\rm ho} \rangle$. The energy with constraint values const}$ $\beta_{\rm was performed with to the parity-projected AMD wave function by minimizing the energy defined as, $$E\equiv \langle \hat{H}\rangle+ V^{N}(N_{\rm const}-\langle N^{\rm ho} const})-D')^2.\label{eq:const}$$ Here artificial are to satisfy the the constraints. With a given set $(N_{\rm const},\beta_{\rm const})$ the optimum wave function $\Phi^\pm_{\rm const},\beta_{\rm const})$ obtained. Finally, we superposed the spin-parity states projected from the obtained wave functions, $$|^8{\rm = \sum_{N_{\rm const},\beta_{\rm const}}c^{J\pm}_n(N_{\rm const},\beta_{\rm const}) |P^{J}_{MK}\Phi^\pm_{\rm AMD}(N_{\rm const
 (for example [@Kimura04]), the deformAtion is usuAlly cOnsTraInEd by UsinG the rotational INvarIant value $D\equiv Tr(QQ)/Tr^2(Q)$, Where ThE MatrIX $Q$ Is calCulated BY qUADruPoLe OpeRaTOrS as $Q_{\sIgmA\rho}=\lanGle \sum_i \hat\SigMa_I
\hat\rho_i \rangLE$ ($\hAt\sigma=\hat{X},\haT{y},\hat{z}$ and $\hat\Rho=\Hat{x},\haT{y},\Hat{Z}$) [@dote97]. HEre $d$ is apProximATely reLated to thE qUAdrupoLE deformATIoN parAmeter $\beta$ as $D(\beta)= (5\BEtA^2/2\Pi+1)/3$. In the present Work, we UsED tHE ModIfiEd quadrupoLe MatriX $q'_{\sigma\rHO}\eQUIV Q_{\sIGma\rho}-A\delta_{\sIgma\rho}$ ($A$ is tHE maSs numbEr) InsTEad of tHe oriGiNAl $Q_{\Sigma\rho}$ and ImpoSed the conStrainT On the $D'\eQUiv Tr(Q'Q')/tr^2(Q')$. ThiS is UseFul fOR HE iSotOpES to OBtAin BAsiS wave funCtIoNs witH varIOUS ConfIguRatiOns on Mesh points of tHe tWo-diMEnsIonal ParamEterS, $\bEta$ anD $\langlE \hat N^{\Rm Ho} \rangle$. The enerGy vaRiation wiTh tHe ConStRaint VAlues $N_{\Rm cOnsT}$ and $\betA_{\rm consT}$ Was PeRFORmEd with respect to the PaRITy-ProjecteD AMD waVE fUnCTion by miNiMizIng tHE EnergY defINeD as, $$E\equiV \langlE \HaT{H}\Rangle+
V^{n}(N_{\Rm consT}-\lAngLe \hAt N^{\rm HO} \ranGle)^2
+V^{\beTa}(D(\beta_{\rM consT})-d')^2.\label{eq:const}$$ HERe the artificiAL pOTEnTIals Are Introduced tO satISfy tHe coNDiTioN Of the ConstRaINtS. with a given set of consTrAint vaLues $(N_{\Rm const},\beta_{\rm Const})$ the opTIMUm wave fuNctiON $\PHI^\pm_{\rm AMD}(N_{\rm conSt},\betA_{\rm const})$ waS Obtained. finalLy, we supeRposed the SPIn-parity EigEn sTatEs pROJeCted from the obTAIned WaVe functIonS, $$|^8{\rm He}(J^\pM_n)\rAngLe = 
\sUm_{N_{\Rm Const},\beta_{\Rm const}}c^{j\pM}_n(n_{\rM cOnsT},\beta_{\RM const})
|P^{J}_{mK}\phi^\Pm_{\Rm AmD}(N_{\rm COnst
 (for example [@Kimura04]) , the defo rmati onisus uall y co nstrained by u s ingthe rotational invaria nt va lu e  $D\ e qu iv Tr (QQ)/Tr ^ 2( Q ) $,wh er e t he ma trix$Q$  is cal culated by  qu ad rupole opera t or s as $Q_{\ sig ma\rho}=\lan gle  \sum_ i\ha t \sigm a_i 
\hat \rho_i \rangl e$ ($\hat \s i gma=\h a t{x},\h a t {y },\h at{z}$ and $\hat\ r ho = \hat{x},\hat{y },\hat {z } $) [ @Do te9 7]. Here $ D$  is a p proxima t el y r ela t ed to the qua drupole def o rma tion p ar ame t er $\b eta$as $D( \beta)= (5\ beta ^2/2\pi+1 )/3$.I n the p r esent w ork, w e u sed  the mo di fie dq uad r up ole mat rix $Q'_ {\ si gma\r ho}\ e q u i v Q_ {\s igma \rho} -A\delta_{\si gma \rho } $ ( $A$ i s the  mas snumbe r) ins teadof  the original $ Q_{\ sigma\rho }$an d i mp osedt he con str ain t on th e $D'\e q uiv  T r ( Q 'Q ')/Tr^2(Q')$. This  i s us eful for  He is o to pe s  to obta in  ba sisw a ve fu ncti o ns  with va riousc on fi guratio ns  on me sh  po int s oft he t wo-dim ensional  para m eters, $\beta$ and $\langle\ ha t N^ { \rmho}  \rangle$.Thee nerg y va r ia tio n  with  theco n st r aint values $N_{\rm  c onst}$  and$\beta_{\rm c onst}$ was p e rformedwith re s pect to the pa rity- projectedA MD wavefunct ion by m inimizing t he energ y d efi ned  as , $$ E\equiv \lang l e  \ha t{ H}\rang le+ 
V^{N}( N_{ \rm  co nst }- \langle \ hat N^{\ rm  h o}  \ ran gle)^ 2 
+V^{\be ta }(D (\ bet a_{\r m  const })-D' )^2. \l ab e l{e q:const } $$ H ereth earti fic ia l pot enti a lsare int roduced t o s a tisf yth e condi tion of the c on straints.Wi tha give n set of c onstraint values $(N_{\ r m const },\ beta_ {\rm  const})$  th e opti mum wave f unctio n $\P hi ^\p m _ {\rmA M D} (N_ {\ rm const}, \ b eta _{\rm  c onst })$ was  obtained. Finally ,  we  superposed t hespin - p ar ity ei g enst a tes p rojected from t he obtaine dw av e function s , $ $| ^8{\rmHe}(J^\ pm_n) \ rangle= 
\sum_{ N_{\rm co ns t},\ b e ta_ {\rm const }}c^{J\p m}_n(N_{\ r m con s t} ,\bet a_{ \rm co ns t}) 
|P^{ J}_{MK } \Ph i^\pm _{\rmAM D}(N_{ \rm c on st
 (for_example [@Kimura04]),_the deformation is usually_constrained by_using_the rotational_invariant_value $D\equiv Tr(QQ)/Tr^2(Q)$,_where the matrix_$Q$ is calculated by_quadrupole operators as_$Q_{\sigma\rho}=\langle_\sum_i \hat\sigma_i
\hat\rho_i \rangle$ ($\hat\sigma=\hat{x},\hat{y},\hat{z}$ and $\hat\rho=\hat{x},\hat{y},\hat{z}$) [@Dote97]. Here $D$ is approximately related to the_quadrupole_deformation parameter_$\beta$_as_$D(\beta)= (5\beta^2/2\pi+1)/3$. In the present_work, we used the modified_quadrupole matrix_$Q'_{\sigma\rho}\equiv Q_{\sigma\rho}-A\delta_{\sigma\rho}$ ($A$ is the mass number) instead_of_the original $Q_{\sigma\rho}$_and imposed the constraint on the $D'\equiv Tr(Q'Q')/Tr^2(Q')$. This_is useful for He isotopes to_obtain basis wave_functions_with_various configurations on mesh_points of the two-dimensional parameters, $\beta$_and $\langle \hat N^{\rm ho} \rangle$._The energy variation with the constraint values_$N_{\rm const}$ and $\beta_{\rm const}$ was_performed with respect to the_parity-projected AMD_wave function by minimizing the_energy defined as,_$$E\equiv \langle_\hat{H}\rangle+
V^{N}(N_{\rm const}-\langle \hat_N^{\rm ho} \rangle)^2
+V^{\beta}(D(\beta_{\rm const})-D')^2.\label{eq:const}$$ Here the_artificial potentials are_introduced to satisfy the condition of_the_constraints. With a_given_set_of constraint_values $(N_{\rm const},\beta_{\rm_const})$_the optimum_wave_function $\Phi^\pm_{\rm AMD}(N_{\rm const},\beta_{\rm const})$ was_obtained._Finally, we superposed the spin-parity eigen states_projected from the obtained_wave_functions, $$|^8{\rm He}(J^\pm_n)\rangle =_
\sum_{N_{\rm const},\beta_{\rm const}}c^{J\pm}_n(N_{\rm const},\beta_{\rm const})
|P^{J}_{MK}\Phi^\pm_{\rm_AMD}(N_{\rm const