{tr}(-1)^F$ [@Witten:1982df].[^5] For the massive free theory $W(\Phi)=m\Phi^2/2$, one easily sees that $\Delta=1$ and reproduces the correct Witten index $\operatorname{tr}(-1)^F=1$. Thus assuming that the proportionality constant between the partition function defined by the functional integral and $\operatorname{tr}(-1)^F$ is independent of the interaction, we have $\Delta=\operatorname{tr}(-1)^F$. It is known that for a homogeneous superpotential $W(\Phi)=\lambda\Phi^n/n$, $\operatorname{tr}(-1)^F=n-1$. Therefore, the relation $\Delta=n-1$ provides a quantitative test of the simulation [@Kawai:2010yj].

A disadvantage of the present simulation algorithm is that it is not *a priori* clear how many solutions $(\Tilde{A}(p)_i,\Tilde{A}(p)^*_i)$ that Eqs.  and  have. Thus we cannot be completely sure whether we have found all the solutions or not. Another limitation of the algorithm is that it is applicable only to supersymmetric boundary conditions; so it cannot explore physics with the finite temperature, for example.

Simulation parameters
=====================

In this paper, we fix the coupling constant in Eq.  to be $$a\lambda=0.3.$$ This is the same choice of the coupling constant as Ref. [@Kawai:2010yj], provided that our parameter $a$ in Eq.  is identified with the lattice spacing in Ref. [@Kawai:2010yj]. The size of the momentum grid $N_0\times N_1$ is varied from $16\times16$ to $36\times36$. For each value of $N_0\times N_1$, we generated $1280$ configurations of the Gaussian random numbers $(\Tilde{N}(p),\Tilde{N}(p)^*)$. We then solved Eqs.  and  by using the Newton–Raphson (NR) method [@NR3]. The solution $(\Tilde{A}(p),\Tilde{A}(p)^*)$ depends on the initial guess in the NR method and, as the initial guess, we used $100$ random configurations $(\T
{ tr}(-1)^F$   [ @Witten:1982df].[^5 ] For the massive free theory $ W(\Phi)=m\Phi^2/2 $, one easily sees that $ \Delta=1 $ and reproduce the right Witten index   $ \operatorname{tr}(-1)^F=1$. Thus assume that the proportion constant between the partition function specify by the functional integral and   $ \operatorname{tr}(-1)^F$ is autonomous of the interaction, we take $ \Delta=\operatorname{tr}(-1)^F$. It is known that for a homogeneous superpotential $ W(\Phi)=\lambda\Phi^n / n$, $ \operatorname{tr}(-1)^F = n-1$. Therefore, the relation back   $ \Delta = n-1 $ provides a quantitative test of the model   [ @Kawai:2010yj ]. 

 A disadvantage of the present simulation algorithm is that it is not * a priori * clear how many solutions $ (\Tilde{A}(p)_i,\Tilde{A}(p)^*_i)$ that Eqs.    and    have. Thus we cannot be wholly sure whether we have found all the solutions or not. Another limit of the algorithm is that it is applicable only to supersymmetric boundary conditions; so it cannot research physics with the finite temperature, for exemplar. 

 Simulation parameters 
 = = = = = = = = = = = = = = = = = = = = = 

 In this newspaper, we fix the coupling constant in   Eq.    to be $ $ a\lambda=0.3.$$ This is the like choice of the coupling constant as   Ref.   [ @Kawai:2010yj ], provided that our parameter   $ a$ in   Eq.    is identified with the wicket space in   Ref.   [ @Kawai:2010yj ]. The size of the momentum grid $ N_0\times N_1 $ is varied from   $ 16\times16 $ to   $ 36\times36$. For each value of   $ N_0\times N_1 $, we beget $ 1280 $ configurations of the Gaussian random numbers   $ (\Tilde{N}(p),\Tilde{N}(p)^*)$. We then solved Eqs.    and    by using the Newton – Raphson (NR) method   [ @NR3 ]. The solution $ (\Tilde{A}(p),\Tilde{A}(p)^*)$ depends on the initial guess in the NR method and, as the initial guess, we used $ 100 $ random configuration   $ (\T
{tr}(-1)^F$ [@Aitten:1982df].[^5] For the massive free theory $W(\Pku)=m\Phi^2/2$, mne eaaily seer that $\Delta=1$ and reproduces vhe xorrext Witten index $\operatofname{tr}(-1)^F=1$. Nhus assuning rhat the pckportiokclity gonstcnv between the psrtition fgnction defineg cy the functional integral and $\operatowname{tr}(-1)^G$ ls independent of uhe intsgagtion, we have $\Delta=\operatorname{td}(-1)^F$. It iv known that gor a homogeneous superpotfntiwl $W(\Phi)=\lambda\Phi^n/n$, $\operatornane{tr}(-1)^S=b-1$. Therefore, ghe relation $\Delta=n-1$ probides a quantitative test of thd simblation [@Kawau:2010yh].

A fhsadvantage of tre present slkulatimn algotithm is that lt is nor *a priori* clear how kany solutions $(\Tilqe{A}(p)_i,\Tilda{A}(')^*_i)$ that Eqs.  and  have. Rhys we cantot cw cumpmevelg sure whxther we habe found alo the solutions or moe. Another limifation os the algorithm is that it is applicablt onlg to supersymmetric boubdary conditions; so ij cannot evplore physics with the finite temperature, for exdmple.

Akmuoabion pagameters
=====================

In this paper, we fix the coupling conffamt in Eq.  to be $$a\lcmbda=0.3.$$ This is thr dakg choice of thg couplnhg constant as Ref. [@Kawwi:2010yj], prjvidee that ouw patameter $a$ in Eq.  is identified with the lanticw spacing in Ref. [@Kawci:2010yj]. The sizz of tne mokentum grid $N_0\times N_1$ is varjed from $16\timfs16$ to $36\timea36$. For each value uf $N_0\nimev N_1$, we gtverated $1280$ configurwtions of the Gaussixn rsndom gumbers $(\Tilfe{N}(p),\Tlnde{N}(p)^*)$. We then solvfd Eqd.  dnd  by usijg the Newton–Raphson (NR) method [@NC3]. The solutiom $(\Dilge{A}(p),\Tildz{A}(p)^*)$ deiends on the inytial guess in the NR methud and, as nhe initiel guess, we tsed $100$ random wlnfigurationv $(\T
{tr}(-1)^F$ [@Witten:1982df].[^5] For the massive free theory easily that $\Delta=1$ reproduces the correct that proportionality constant between partition function defined the functional integral and $\operatorname{tr}(-1)^F$ is of the interaction, we have $\Delta=\operatorname{tr}(-1)^F$. It is known that for a homogeneous $W(\Phi)=\lambda\Phi^n/n$, $\operatorname{tr}(-1)^F=n-1$. Therefore, the relation $\Delta=n-1$ provides a quantitative test of the simulation A of present algorithm is that it is not *a priori* clear how many solutions $(\Tilde{A}(p)_i,\Tilde{A}(p)^*_i)$ that Eqs. and Thus we cannot be completely sure whether we found all the solutions not. Another limitation of the is it is only supersymmetric conditions; so it explore physics with the finite temperature, for example. Simulation parameters ===================== In this paper, we fix the constant in be $$a\lambda=0.3.$$ is same of the coupling Ref. [@Kawai:2010yj], provided that our parameter is identified with the lattice spacing in Ref. The size the momentum grid $N_0\times N_1$ is from $16\times16$ to $36\times36$. For each value of N_1$, we generated $1280$ configurations of the Gaussian random numbers $(\Tilde{N}(p),\Tilde{N}(p)^*)$. We then solved Eqs. using the Newton–Raphson (NR) [@NR3]. The solution depends the guess the NR and, as the initial guess, we used $100$ random configurations $(\T
{tr}(-1)^F$ [@Witten:1982df].[^5] For the massive fRee theory $W(\phi)=m\PHi^2/2$, oNe eAsIly sEes tHat $\Delta=1$ and repROducEs the correct Witten indeX $\operAtORnamE{Tr}(-1)^f=1$. Thus AssuminG ThAT The PrOpOrtIoNAlIty coNstAnt betwEen the partItiOn Function defiNEd By the functIonAl integral anD $\opEratorNaMe{tR}(-1)^f$ is inDepEndenT of the INteracTion, we havE $\DELta=\opeRAtornamE{TR}(-1)^F$. it is Known that for a homoGEnEOus superpotentIal $W(\PhI)=\lAMbDA\phi^N/n$, $\oPeratornamE{tR}(-1)^F=n-1$. ThERefore, tHE rELATioN $\delta=n-1$ provideS a quantitatIVe tEst of tHe SimULation [@kawai:2010Yj].

a DisAdvantage of The pResent simUlatioN AlgoritHM is that It is noT *a pRioRi* clEAr HoW maNy SOluTIoNs $(\TILde{a}(p)_i,\Tilde{a}(p)^*_I)$ tHat EqS.  and  HAVE. thus We cAnnoT be coMpletely sure wHetHer wE HavE founD all tHe soLuTions Or not. ANotheR lImitation of the aLgorIthm is thaT it Is AppLiCable ONly to sUpeRsyMmetric BoundarY ConDiTIONs; So it cannot explore pHySICs With the fInite tEMpErATure, for eXaMplE.

SimULAtion ParaMEtErs
=====================

In thiS paper, WE fIx The coupLiNg consTaNt iN Eq.  To be $$a\LAmbdA=0.3.$$ This iS the same ChoicE Of the coupling cONstant as Ref. [@KaWAi:2010YJ], PrOVideD thAt our parameTer $a$ IN Eq.  iS ideNTiFieD With tHe latTiCE sPAcing in Ref. [@Kawai:2010yj]. ThE sIze of tHe momEntum grid $N_0\timEs N_1$ is varieD FROm $16\times16$ tO $36\timES36$. FOR each value of $N_0\tImes N_1$, We generateD $1280$ ConfigurAtionS of the GaUssian ranDOM numbers $(\tilDe{N}(P),\TiLde{n}(P)^*)$. we Then solved Eqs.  AND  by uSiNg the NeWtoN–RaphsoN (NR) MetHod [@nR3]. THe Solution $(\TIlde{A}(p),\TiLdE{A}(P)^*)$ dEpEndS on thE Initial gUeSs iN tHe Nr methOD and, as The inItiaL gUeSS, we Used $100$ ranDOm CONfigUrAtIons $(\t
{tr}(-1)^F$ [@Witten:1982d f].[^5] Fo r the  ma ssi ve  fre e th eory $W(\Phi)= m \Phi ^2/2$, one easily sees  that  $ \ Delt a =1 $ and  reprod u ce s the  c or rec tW it ten i nde x $\ope ratorname{ tr} (- 1)^F=1$. Thu s  a ssuming th atthe proporti ona lity c on sta n t bet wee n the  parti t ion fu nction de fi n ed byt he func t i on al i ntegral and $\ope r at o rname{tr}(-1)^ F$ isin d ep e n den t o f the inte ra ction ,  we hav e  $ \ D e lta = \operatorname {tr}(-1)^F$ .  It  is kn ow n t h at for  a ho mo g ene ous superpo tent ial $W(\P hi)=\l a mbda\Ph i ^n/n$,$\oper ato rna me{t r }( -1 )^F =n - 1$. Th ere f ore , the re la ti on $\ Delt a = n - 1$ p rov ides  a qu antitative te stof t h e s imula tion[@Ka wa i:201 0yj].
A di sa dvantage of the  pre sent simu lat io n a lg orith m  is th atitis not*a prio r i*cl e a r  h ow many solutions$( \ T il de{A}(p) _i,\Ti l de {A } (p)^*_i) $tha t Eq s .   and   ha v e.  Thus we  canno t  b ecomplet el y sure  w het her  we h a ve f ound a ll the s oluti o ns or not. Ano t her limitatio n  o f th e  alg ori thm is that  iti s ap plic a bl e o n ly to  supe rs y mm e tric boundary condi ti ons; s o itcannot explor e physicsw i t h the fi nite te m perature, forexamp le.

Simul a tion par amete rs
===== ========= = = =====

I n t his  pa per , we  fix the coup l i ng c on stant i n E q.  tobe$$a \la mbd a= 0.3.$$ Th is is th esa me  c hoi ce of the coup li ngco nst ant a s  Ref.[@Kaw ai:2 01 0y j ],provide d  t h a t ou rpa rame ter  $ a$ in  Eq.  is  identi fied with  th e  lat ti ce  spacin g in Ref. [@K aw ai:2010yj] .The  sizeo f  the mom entum grid $N_0\times N _ 1$ is v ari ed fr om $ 16\times1 6$to $36 \ti m es36$.  For e ach v al ueo f  $N_0 \ t im esN_ 1$, we gen e r ate d $12 80 $ co nfigura tions of the Gauss i anrandom number s $ (\Ti l d e{ N}( p ), \ Til de { N}( p ) ^*)$. We then s olved Eqs.  a nd   by using the  N ewton–R aphson(NR)m ethod [ @NR3]. Th e solutio n$(\T i l de{ A}(p),\Til de{A}(p) ^*)$ depe n ds on th e ini tia l gues sinthe N R meth o d a nd, a s thein itialguess ,we used$100$ random configurat ions $ (\T
{tr}(-1)^F$ [@Witten:1982df].[^5] For_the massive_free theory $W(\Phi)=m\Phi^2/2$, one_easily sees_that_$\Delta=1$ and_reproduces_the correct Witten_index $\operatorname{tr}(-1)^F=1$. Thus assuming_that the proportionality constant_between the partition_function_defined by the functional integral and $\operatorname{tr}(-1)^F$ is independent of the interaction, we have $\Delta=\operatorname{tr}(-1)^F$._It_is known_that_for_a homogeneous superpotential $W(\Phi)=\lambda\Phi^n/n$, $\operatorname{tr}(-1)^F=n-1$._Therefore, the relation $\Delta=n-1$ provides a_quantitative test_of the simulation [@Kawai:2010yj].

A disadvantage of the present simulation_algorithm_is that it_is not *a priori* clear how many solutions $(\Tilde{A}(p)_i,\Tilde{A}(p)^*_i)$_that Eqs.  and  have. Thus we_cannot be completely_sure_whether_we have found all_the solutions or not. Another limitation_of the algorithm is that it_is applicable only to supersymmetric boundary conditions;_so it cannot explore physics with_the finite temperature, for example.

Simulation_parameters
=====================

In this_paper, we fix the coupling_constant in Eq.  to_be $$a\lambda=0.3.$$_This is the_same choice of the coupling constant_as Ref. [@Kawai:2010yj], provided that_our parameter $a$ in Eq.  is identified with_the_lattice spacing in Ref. [@Kawai:2010yj]._The_size_of the_momentum grid $N_0\times_N_1$_is varied_from $16\times16$_to $36\times36$. For each value of $N_0\times N_1$,_we_generated $1280$ configurations of the Gaussian random_numbers $(\Tilde{N}(p),\Tilde{N}(p)^*)$. We then solved_Eqs. _and  by using the_Newton–Raphson (NR) method [@NR3]. The solution_$(\Tilde{A}(p),\Tilde{A}(p)^*)$ depends on the initial guess_in the_NR method_and, as the initial guess, we used $100$ random configurations $(\T