{eq:diffusion_energy}
E_{\rm D} \equiv E_0 = \frac{\langle \phi_0 | \hat{H} | \Psi_{\rm T}
\rangle}{\langle \phi_0 | \Psi_{\rm T} \rangle} & = & \frac{\int f({\bf R})
E_{\rm L}({\bf R}) \, d{\bf R}}{\int f({\bf R}) \, d{\bf R}} \\ & \simeq &
\frac{1}{M} \sum_{i=1}^M E_{\rm L}({\bf R}_i) \;.\end{aligned}$$ Some example DMC data are shown in figure \[fig:graphit\_silane\].

![DMC data for a silane (SiH$_4$) molecule, with the ions represented by pseudopotentials. The upper panel shows the fluctuations in the population of configurations arising from the branching process used to simulate equation (\[eq:short time GB\]). The reference energy, $E_{\rm T}$, is altered during the run to control the population. Specifically, the reference energy is set to return the population to the target population (128,000 configurations) on the same time-scale as the autocorrelation period of the energy data [@foulkes_2001]. The total energy is shown in the lower panel as a function of the move number. The black line shows the instantaneous value of the local energy averaged over the current population of configurations, the red line is the reference energy $E_{\rm T}$ and the green line is the best estimate of the DMC energy as the simulation progresses. The configurations at move number zero are from the output of a VMC simulation, and the energy decays rapidly from its initial VMC value of about -6.250 a.u. and reaches a plateau with a DMC energy of about -6.305 a.u. The data up to move 1000 are deemed to form the equilibration phase, and are discarded. []{data-label="fig:graphit_silane"}](graphit_silane){width="100.00000%"}

Trial wave functions {#sec:psi_trial}
====================

Trial wave functions are of central importance in VMC and DMC calculations because they introduce importance sampling and control both the statistical efficiency and accuracy obtained. The accuracy of a DMC calculation depends on the nodal surface of the trial wave function via the fixed-node approximation, while in V
{ eq: diffusion_energy } 
 E_{\rm D } \equiv E_0 = \frac{\langle \phi_0 | \hat{H } | \Psi_{\rm T } 
 \rangle}{\langle \phi_0 | \Psi_{\rm T } \rangle } & = & \frac{\int f({\bf R }) 
 E_{\rm L}({\bf R }) \, d{\bf R}}{\int f({\bf R }) \, d{\bf R } } \\ & \simeq & 
 \frac{1}{M } \sum_{i=1}^M E_{\rm L}({\bf R}_i) \;.\end{aligned}$$ Some example DMC data are shown in name \[fig: graphit\_silane\ ]. 

! [ DMC datum for a silane (SiH$_4 $) molecule, with the ions represent by pseudopotentials. The upper gore shows the fluctuations in the population of shape arise from the branching procedure used to simulate equality (\[eq: short time GB\ ]). The reference energy, $ E_{\rm T}$, is interpolate during the run to see the population. Specifically, the reference energy is sic to return the population to the target population (128,000 configurations) on the like time - scale as the autocorrelation time period of the energy datum [ @foulkes_2001 ]. The total department of energy is shown in the lower panel as a function of the motion number. The black line shows the instantaneous value of the local energy average over the current population of configurations, the red line is the reference energy $ E_{\rm T}$ and the green line is the best estimate of the DMC energy as the model progresses. The configurations at motion numeral zero are from the output of a VMC simulation, and the energy decay rapidly from its initial VMC value of about -6.250 a.u.   and reaches a plateau with a DMC energy of about -6.305 a.u. The datum up to move 1000 are deemed to form the equilibration phase, and are discarded. [ ] { data - label="fig: graphit_silane"}](graphit_silane){width="100.00000% " } 

 Trial wave functions { # sec: psi_trial } 
 = = = = = = = = = = = = = = = = = = = = 

 Trial wave functions are of central importance in VMC and DMC calculations because they introduce importance sampling and see both the statistical efficiency and accuracy obtained. The accuracy of a DMC calculation depend on the nodal surface of the trial wave function via the fixed - node estimate, while in V
{eq:dlffusion_energy}
E_{\rm D} \equiy E_0 = \frac{\langle \khu_0 | \hat{I} | \Psi_{\rj T}
\rangld}{\langle \phi_0 | \Psi_{\rm T} \rangle} & = & \frac{\unt f({\bf R})
E_{\rm L}({\bf R}) \, d{\bw R}}{\int f({\bv R}) \, d{\bf R}} \\ & \wimeq &
\frac{1}{M} \sum_{i=1}^M C_{\xm L}({\br R}_i) \;.\znv{aligned}$$ Some ewample DMC gata are shown iv yigure \[fig:graphit\_silane\].

![DMC data for w silanr (DiH$_4$) molecule, wyth uhe iona represented by pseudopotentials. Fhe upptr panel shows the fluctuations in the populwtioj of configurationd arising ftkm eye branching process used to simulzte equation (\[eq:short time GB\]). Thd refzrence enerty, $E_{\rl T}$, is alterxd durpng the run to contron the pppulation. Speclficanly, the reference energy is set to return jhe populadikn to the target populathon (128,000 conwugufatjois) kn the saje time-scame as the aytocorrelation peripd if the energy data [@sotlkes_2001]. The total energy is shown in the nowsr panel as a function if the move number. Thg black lige shows the instantaneous value of the local enesgy atefagtd over rhf current population of configurations, the req lonv is the referencc energy $E_{\rm T}$ and tje dreen line is the bzat estimate of the DLC enerdy as the simujatipn progresses. The configurarions at movv nunber zero are from the output of s VMC simulation, and the enexgy dedays rapidlj from ita initial VMC valje pf about -6.250 a.u. and reaches a plwteau witi a DKC enerey og aboue -6.305 a.u. The fata m[ to move 1000 are deeled tl xorm the eeuilibration phase, and are discecded. []{data-labek="fhg:ggaphit_silcne"}](graihit_silane){width="100.00000%"}

Erial wave funetions {#szc:psi_tfial}
====================

Trial save fuictions are jf central im[lrtance in VKC and DIC cqlcuoations cecause they imtroduce pmkortance sqmpling and controk bufh the statistiecl efficiency and aczurwcj outainqg. The accurawy ow a CMC cxlculation beoendx on the nodal surfawe or the trial wave fingtion via the fixqd-node approxomation, while in V
{eq:diffusion_energy} E_{\rm D} \equiv E_0 = \frac{\langle \hat{H} \Psi_{\rm T} \phi_0 | \Psi_{\rm \frac{\int R}) E_{\rm L}({\bf \, d{\bf R}}{\int R}) \, d{\bf R}} \\ & & \frac{1}{M} \sum_{i=1}^M E_{\rm L}({\bf R}_i) \;.\end{aligned}$$ Some example DMC data are shown figure \[fig:graphit\_silane\]. ![DMC data for a silane (SiH$_4$) molecule, with the ions represented pseudopotentials. upper shows fluctuations in the population of configurations arising from the branching process used to simulate equation (\[eq:short GB\]). The reference energy, $E_{\rm T}$, is altered the run to control population. Specifically, the reference energy set return the to target (128,000 configurations) on same time-scale as the autocorrelation period of the energy data [@foulkes_2001]. The total energy is shown in lower panel function of move The line shows the of the local energy averaged over of configurations, the red line is the reference $E_{\rm T}$ the green line is the best of the DMC energy as the simulation progresses. configurations at move number zero are from the output of a VMC simulation, and the rapidly from its initial value of about a.u. reaches plateau a DMC of about -6.305 a.u. The data up to move 1000 are to form the equilibration phase, and are discarded. []{data-label="fig:graphit_silane"}](graphit_silane){width="100.00000%"} Trial {#sec:psi_trial} Trial wave functions of central importance in and calculations because they introduce and both and obtained. accuracy of a DMC depends on the nodal surface the trial wave function in V
{eq:diffusion_energy}
E_{\rm D} \equiV E_0 = \frac{\langLe \phi_0 | \Hat{h} | \PsI_{\rM T}
\raNgle}{\Langle \phi_0 | \Psi_{\rm t} \RangLe} & = & \frac{\int f({\bf R})
E_{\rm L}({\bf R}) \, d{\bF R}}{\int F({\bF r}) \, d{\bf r}} \\ & \SiMeq &
\frAc{1}{M} \sum_{i=1}^m e_{\rM l}({\Bf R}_I) \;.\eNd{AliGnED}$$ SOme exAmpLe DMC daTa are shown In fIgUre \[fig:graphiT\_SiLane\].

![DMC datA foR a silane (SiH$_4$) mOleCule, wiTh The IOns rePreSenteD by pseUDopoteNtials. The UpPEr paneL Shows thE FLuCtuaTions in the populatIOn OF configurationS arisiNg FRoM THe bRanChing proceSs Used tO SimulatE EqUATIon (\[EQ:short time GB\]). THe reference ENerGy, $E_{\rm T}$, Is AltERed durIng thE rUN to Control the pOpulAtion. SpecIficalLY, the refERence enErgy is Set To rEturN ThE pOpuLaTIon TO tHe tARgeT populatIoN (128,000 cOnfigUratIONS) On thE saMe tiMe-scaLe as the autocoRreLatiON peRiod oF the eNergY dAta [@foUlkes_2001]. THe totAl Energy is shown in The lOwer panel As a FuNctIoN of thE Move nuMbeR. ThE black lIne showS The InSTANtAneous value of the loCaL ENeRgy averaGed oveR ThE cURrent popUlAtiOn of CONfiguRatiONs, The red liNe is thE ReFeRence enErGy $E_{\rm T}$ AnD thE grEen liNE is tHe best Estimate Of the dmC energy as the sIMulation progrESsES. thE ConfIguRations at moVe nuMBer zEro aRE fRom THe outPut of A Vmc sIMulation, and the energY dEcays rApidlY from its initiAl VMC value OF ABout -6.250 a.u. anD reaCHeS A plateau with a DmC eneRgy of about -6.305 A.U. The data Up to mOve 1000 are deEmed to forM THe equiliBraTioN phAse, AND aRe discarded. []{daTA-LabeL="fIg:graphIt_sIlane"}](grAphIt_sIlaNe){wIdTh="100.00000%"}

Trial waVe functiOnS {#sEc:PsI_trIal}
====================

TrIAl wave fuNcTioNs Are Of cenTRal impOrtanCe in vMc aND DMc calculATiONS becAuSe They IntRoDuce iMporTAncE sampliNg and contRol BOth tHe StAtisticAl efficiency aNd Accuracy obTaIneD. The acCURacy of a DmC calculation depends on tHE nodal sUrfAce of The tRial wave fUncTion viA thE Fixed-nOde appRoximAtIon, WHIle in v
{eq:diffusion_energy}
E_{\ rm D} \equ iv E_ 0 =  \f ra c{\l angl e \phi_0 | \ha t {H}| \Psi_{\rm T}
\rangle }{\la ng l e \p h i_ 0 | \ Psi_{\r m  T } \ra ng le } &  = &\frac {\i nt f({\ bf R})
E_{ \rm  L }({\bf R}) \ ,  d {\bf R}}{\ int  f({\bf R})\,d{\bfR} } \ \  & \s ime q &
\ frac{1 } {M} \s um_{i=1}^ ME _{\rmL }({\bfR } _i ) \; .\end{aligned}$$S om e  example DMC d ata ar es ho w n  in  fi gure \[fig :g raphi t \_silan e \] . ![D M C data for asilane (SiH $ _4$ ) mole cu le, with t he io ns rep resented by  pse udopotent ials.T he uppe r  panelshowsthe  fl uctu a ti on s i nt hep op ula t ion  of conf ig ur ation s ar i s i n g fr omthebranc hing processuse d to sim ulate  equa tion  ( \[eq: shorttimeGB \]). The refere nceenergy, $ E_{ \r m T }$ , isa ltered  du rin g the r un to c o ntr ol t h epopulation. Specif ic a l ly , the re ferenc e  e ne r gy is se ttoretu r n  thepopu l at ion to t he tar g et  p opulati on  (128, 00 0 c onf igura t ions ) on t he sametime- s cale as the au t ocorrelationp er i o do f th e e nergy data[@fo u lkes _200 1 ].  Th e  tota l ene rg y  i s  shown in the lower  p anel a s a f unction of th e move num b e r . The bl ackl in e  shows the ins tanta neous valu e  of thelocal  energyaveragedo v er the c urr ent  po pul a t io n of configur a t ions ,the red  li ne is t heref ere nce  e nergy $E_ {\rm T}$  a nd  t he  gr een l i ne is th ebes test imate of the  DMCener gy  a s  th e simul a ti o n  pro gr es ses.  Th econfi gura t ion s at mo ve number  ze r o ar efr om theoutput of a V MC  simulatio n,  an d thee n ergy dec ays rapidly from its in i tial VM C v alueof a bout -6.2 50a.u. a ndr eaches  a pla teauwi tha DMC e n e rg y o fabout -6.3 0 5  a. u. Th edata  up tomove 1000 are deem e d t o form the eq uil ibra t i on  ph a se ,  an da red i scarded. []{dat a-label="f ig : gr aphit_sila n e"} ]( graphit _silane ){wid t h="100. 00000%"}
Trial wa ve  fun c t ion s {#sec:ps i_trial} 
======== = ===== = == ===
Tri al wav efun ction s areo f c entra l impo rt ance i n VMC  a nd DMC c alculations because the y intr oduce  im portancesam p lin g and con trol  both thesta tis tical  ef f icien cy a n dacc u racyobta i ned. Thea cc ura c y  o f a DMC cal c u l ati on de pen d s on t he n odal surface of t h e trial wave f unct i o n v iat he f ix ed-node approx ima ti o n , whilein  V
{eq:diffusion_energy}
E_{\rm D}_\equiv E_0_= \frac{\langle \phi_0 |_\hat{H} |_\Psi_{\rm_T}
\rangle}{\langle \phi_0_|_\Psi_{\rm T} \rangle}_& = &_\frac{\int f({\bf R})
E_{\rm L}({\bf_R}) \, d{\bf_R}}{\int_f({\bf R}) \, d{\bf R}} \\ & \simeq &
\frac{1}{M} \sum_{i=1}^M E_{\rm L}({\bf R}_i) \;.\end{aligned}$$_Some_example DMC_data_are_shown in figure \[fig:graphit\_silane\].

![DMC data_for a silane (SiH$_4$) molecule,_with the_ions represented by pseudopotentials. The upper panel shows_the_fluctuations in the_population of configurations arising from the branching process used_to simulate equation (\[eq:short time GB\])._The reference energy,_$E_{\rm_T}$,_is altered during the_run to control the population. Specifically,_the reference energy is set to_return the population to the target population_(128,000 configurations) on the same time-scale_as the autocorrelation period of_the energy_data [@foulkes_2001]. The total energy_is shown in_the lower_panel as a_function of the move number. The_black line shows_the instantaneous value of the local_energy_averaged over the_current_population_of configurations,_the red line_is_the reference_energy_$E_{\rm T}$ and the green line_is_the best estimate of the DMC energy_as the simulation progresses._The_configurations at move number_zero are from the output_of a VMC simulation, and the_energy decays_rapidly from_its initial VMC value of about -6.250 a.u. and reaches a plateau_with a DMC energy of about_-6.305 a.u. The data_up to_move_1000 are deemed_to_form the_equilibration phase, and are discarded. []{data-label="fig:graphit_silane"}](graphit_silane){width="100.00000%"}

Trial wave_functions {#sec:psi_trial}
====================

Trial_wave functions are of central importance_in VMC and DMC_calculations_because they introduce importance sampling and_control both the statistical efficiency and_accuracy obtained. The accuracy of_a_DMC_calculation depends on the nodal_surface of the trial wave function_via the fixed-node_approximation, while in V