 for different time steps and extrapolating to zero time step or by simply choosing a small enough time step. An example of time-step extrapolation is shown in figure \[fig:time\_step\_errors\].

-   Population control bias. The $f$ distribution is represented by a finite population of configurations which fluctuates due to branching. The population may be controlled in various ways, but this introduces a population control bias which is positive and falls off as the reciprocal of the population. In practice the population control bias is normally so small that it is difficult to detect [@umrigar_1993; @drummond_2004].

-   Finite size errors within periodic boundary conditions calculations. It is important to correct for finite size effects carefully, as mentioned in section \[sec:pbc\].

-   The pseudopotential approximation inevitably introduces errors. In DMC there is an additional error arising from the localisation [@mitas_1991] or semi-localisation [@casula_2006] of the non-local pseudopotential operator. The localisation error appears to be quite small in the cases for which it has been tested [@drummond_2006b].

![DMC energy against time step for a 64-electron ferromagnetic 2D hexagonal Wigner crystal at density parameter $r_s=50$ a.u. with a Slater-Jastrow wave function. The solid line is a linear fit to the data. []{data-label="fig:time_step_errors"}](ferro_rs50_N064_exp0002_E_v_dt){width=".7\textwidth"}

Practical methods for handling statistical errors in QMC results {#subsec:statistical errors}
----------------------------------------------------------------

Two main practical problems are encountered when dealing with errors in the QMC data: the data are serially correlated and the underlying probability distributions are non-Gaussian. The probability distribution of the local energies has $|E-E_0|^{-4}$ tails, where $E_0$ is a constant. These tails arise from singularities in the local energy such as the divergence at the nodal surface [@trail_2005_1; @trail_2005_2], as shown in figure \[fig:local\_energy\]. In consequence, although the mean energy and its variance are well defined, the variance of the variance is infinity. For other quantities the problem may be even more severe; for example, the probability distributions for the Pul
for different time steps and extrapolating to zero meter gradation or by simply choosing a small adequate time step. An example of clock time - step extrapolation is shown in human body \[fig: time\_step\_errors\ ]. 

 -    Population control bias. The $ f$ distribution is represent by a finite population of configurations which fluctuates ascribable to branching. The population may be controlled in various ways, but this introduces a population control bias which is convinced and falls off as the reciprocal of the population. In practice the population command bias is normally so small that it is difficult to detect [ @umrigar_1993; @drummond_2004 ]. 

 -    Finite size errors within periodic boundary weather calculations. It is important to correct for finite size effects carefully, as mention in section \[sec: pbc\ ]. 

 -    The pseudopotential approximation inevitably introduces errors. In DMC there is an additional error arising from the localisation [ @mitas_1991 ] or semi - localisation [ @casula_2006 ] of the non - local pseudopotential operator. The localisation error appears to be quite small in the cases for which it has been test [ @drummond_2006b ]. 

! [ DMC energy against time dance step for a 64 - electron ferromagnetic 2D hexagonal Wigner crystal at density parameter $ r_s=50 $ a.u.   with a Slater - Jastrow wave affair. The solid line is a linear fit to the data. [ ] { datum - label="fig: time_step_errors"}](ferro_rs50_N064_exp0002_E_v_dt){width=".7\textwidth " } 

 Practical methods for handling statistical errors in QMC results { # subsec: statistical errors } 
 ---------------------------------------------------------------- 

 Two main practical problems are encountered when dealing with errors in the QMC datum: the data are serially correlate and the underlie probability distributions are non - Gaussian. The probability distribution of the local energies has $ |E - E_0|^{-4}$ tails, where $ E_0 $ is a constant. These tail arise from singularities in the local energy such as the deviation at the nodal surface [ @trail_2005_1; @trail_2005_2 ], as shown in figure \[fig: local\_energy\ ]. In consequence, although the mean department of energy and its variance are well defined, the variance of the variance is infinity. For early quantities the problem may be even more severe; for example, the probability distribution for the Pul
 fog different time steps akd extrapolating to zerm time step or by simply choosing a small xnouth tine step. An example of gime-step vxtrapolarion us shown ii figure \[fig:tijc\_step\_zrcors\].

-   Populatipn control bias. The $f$ divtfiyution is represented by a finite po[ulatiom lf configuratijns erich fluctuates due to branching. The pkpulatimn may be conyrolled in various ways, buh thls introduces a pooulation cobtroj bias which ks positivt cnd falls orf as the reciprocal of the popjlatipn. In pracjnxe hve populatiin cogtrol bias is normalny so skall that it ix dmffixult to detect [@umrigac_1993; @drummond_2004].

-   Finite fize erross within periodic voyndari congitiubs zaldukafions. Lt ms importanf to correcr for finite size egfqbys carefully, as megtyoned in section \[sec:pbc\].

-   The pseudopottntiam approximation inevitavly introduces errors. In DMC trere is an additional error arising from the locanisatmov [@mnbwr_1991] og semi-localisation [@casula_2006] of the non-local psetsokotvntial operator. Tme localisation ertog skpears to be qoite smcml in the cases for ahich ij has veen testtd [@drimmond_2006b].

![DMC energy against tume step for q 64-electron ferromaynetic 2D hexcgonal Wignrr crystal at density pcrametsr $r_s=50$ a.u. witj a Slated-Gastrow wave funztipn. The solid line is a lineaw fit to vhe dcta. []{data-uabek="fig:tiie_step_errogs"}](fervm_rs50_N064_exp0002_E_v_dt){width=".7\tedtwidjh"}

Pracdical methlds for handling statistical ercprs in QMC rgsunts {#subsec:scatistlcal errors}
----------------------------------------------------------------

Two iain practical probleis ard encounteged when vealing with errors in tha QMC data: thx data arq seeialoy corrduated and the inderlyiny provability distributlons xde non-Gaussian. Cke probability dixtrkbueiln os the local etergkes nas $|E-D_0|^{-4}$ tails, wheve $D_0$ is a constant. These tahls zrise from singulatibies in tye local energy such ss the divergence wt thx nodau surgacg [@trail_2005_1; @trail_2005_2], as shown in figurs \[fig:locap\_encrgy\]. In consezuenge, ajthough thz mean energy and its variance are well vefined, the variance of the variance is inyikity. For othxr quagtities tve problem may be evwn more severe; fov example, the probabiljty divtribktions for the Pul
for different time steps and extrapolating to step by simply a small enough time-step is shown in \[fig:time\_step\_errors\]. - Population bias. The $f$ distribution is represented a finite population of configurations which fluctuates due to branching. The population may controlled in various ways, but this introduces a population control bias which is and off the of the population. In practice the population control bias is normally so small that it is to detect [@umrigar_1993; @drummond_2004]. - Finite size errors periodic boundary conditions calculations. is important to correct for size carefully, as in \[sec:pbc\]. The pseudopotential approximation introduces errors. In DMC there is an additional error arising from the localisation [@mitas_1991] or semi-localisation [@casula_2006] the non-local The localisation appears be small in the which it has been tested [@drummond_2006b]. time step for a 64-electron ferromagnetic 2D hexagonal crystal at parameter $r_s=50$ a.u. with a Slater-Jastrow function. The solid line is a linear fit the data. []{data-label="fig:time_step_errors"}](ferro_rs50_N064_exp0002_E_v_dt){width=".7\textwidth"} Practical methods for handling statistical errors in QMC results {#subsec:statistical errors} ---------------------------------------------------------------- practical problems are encountered dealing with errors the data: data serially correlated the underlying probability distributions are non-Gaussian. The probability distribution of the energies has $|E-E_0|^{-4}$ tails, where $E_0$ is a constant. These from in the local such as the divergence the surface [@trail_2005_1; @trail_2005_2], as figure In mean and variance are well defined, variance of the variance is For other quantities the severe; for example, the probability distributions for the
 for different time steps and eXtrapolatiNg to zEro TimE sTep oR by sImply choosing a SMall Enough time step. An examplE of tiMe-STep eXTrApolaTion is sHOwN IN fiGuRe \[Fig:TiME\_sTep\_erRorS\].

-   PopulaTion controL biAs. the $f$ distribuTIoN is represeNteD by a finite poPulAtion oF cOnfIGuratIonS whicH fluctUAtes duE to branchInG. the popULation mAY Be ContRolled in various waYS, bUT this introduceS a popuLaTIoN COntRol Bias which iS pOsitiVE and falLS oFF AS thE Reciprocal of tHe populatioN. in pRacticE tHe pOPulatiOn conTrOL biAs is normallY so sMall that iT is difFIcult to DEtect [@umRigar_1993; @dRumMonD_2004].

-   FinITe SiZe eRrORs wIThIn pERioDic boundArY cOnditIons CALCUlatIonS. It iS impoRtant to correcT foR finITe sIze efFects CareFuLly, as MentioNed in SeCtion \[sec:pbc\].

-   The pSeudOpotentiaL apPrOxiMaTion iNEvitabLy iNtrOduces eRrors. In dmC tHeRE IS aN additional error arIsING fRom the loCalisaTIoN [@mITas_1991] or semI-lOcaLisaTIOn [@casUla_2006] oF ThE non-locaL pseudOPoTeNtial opErAtor. ThE lOcaLisAtion ERror AppearS to be quiTe smaLL in the cases for WHich it has been TEsTED [@dRUmmoNd_2006b].

![dMC energy agAinsT Time Step FOr A 64-elECtron FerroMaGNeTIc 2D hexagonal Wigner cRyStal at DensiTy parameter $r_s=50$ A.u. with a SlaTER-jastrow wAve fUNcTIon. The solid linE is a lInear fit to THe data. []{daTa-labEl="fig:timE_step_erroRS"}](Ferro_rs50_N064_Exp0002_e_v_dT){wiDth=".7\TEXtWidth"}

PracticaL MEthoDs For handLinG statisTicAl eRroRs iN QmC results {#Subsec:stAtIsTiCaL erRors}
----------------------------------------------------------------

TWO main praCtIcaL pRobLems aRE encouNtereD wheN dEaLIng With errORs IN The QmC DaTa: thE daTa Are seRialLY coRrelateD and the unDerLYing PrObAbility Distributions ArE non-GaussiAn. the ProbabILIty distrIbution of the local energiES has $|E-E_0|^{-4}$ tAilS, wherE $E_0$ is A constant. theSe tailS arISe from SingulAritiEs In tHE Local ENErGy sUcH as the diveRGEncE at thE nOdal Surface [@Trail_2005_1; @trail_2005_2], as shown iN FigUre \[fig:local\_enErgY\]. In cONSeQueNCe, ALthOuGH thE MEan energy and its Variance arE wELl Defined, the VAriAnCe of the VariancE is inFInity. FoR other quaNtities thE pRoblEM May Be even more Severe; foR example, tHE probABiLity dIstRibutiOnS foR the PUl
 for different time stepsand extrap olati ngtoze ro t imestep or by sim p ly c hoosing a small enough  time  s t ep.A nexamp le of t i me - s tep  e xt rap ol a ti on is  sh own infigure \[f ig: ti me\_step\_er r or s\].

-  Pop ulation cont rol  bias.  T he$ f$ di str ibuti on isr eprese nted by a  f i nite p o pulatio n of  con figurations which fl u ctuates due to  branc hi n g. T hepop ulation ma ybe co n trolled in v a rio u s ways, but t his introdu c esa popu la tio n  contr ol bi as whi ch is posit iveand falls  off a s  the re c iprocal  of th e p opu lati o n.  I n p ra c tic e  t hep opu lation c on tr ol bi as i s n o rmal lyso s mallthat it is di ffi cult todetec t [@u mrig ar _1993 ; @dru mmond _2 004].

-   Fini te s ize error s w it hin  p eriod i c boun dar y c onditio ns calc u lat io n s .  I t is important toco r r ec t for fi nite s i ze  e f fects ca re ful ly,a s  ment ione d  i n sectio n \[se c :p bc \].

- The ps eu dop ote ntial appr oximat ion inev itabl y  introduces er r ors. In DMC t h er e is an a ddi tional erro r ar i sing  fro m  t hel ocali satio n[ @m i tas_1991] or semi-l oc alisat ion [ @casula_2006]  of the no n - l ocal pse udop o te n tial operator.  Thelocalisati o n errorappea rs to be  quite sm a l l in the  ca ses  fo r w h i ch  it has beent e sted  [ @drummo nd_ 2006b]. 

! [DM C e ner gy  againsttime ste pfo ra64- elect r on ferro ma gne ti c 2 D hex a gonalWigne r cr ys ta l  at  densit y  p a r amet er  $ r_s= 50$  a .u. w itha  Sl ater-Ja strow wav e f u ncti on .The sol id line is ali near fit t othe  data. [ ]{data-l abel="fig:time_step_err o rs"}](f err o_rs5 0_N0 64_exp000 2_E _v_dt) {wi d th=".7 \textw idth" }
Pr a c tical m et hod sfor handli n g  st atist ic al e rrors i n QMC results {#su b sec :statisticalerr ors} - -- --- - -- - --- -- - --- - - --------------- ---------- -- - -- ---------- - --
T wo main  practi cal p r oblemsare encou ntered wh en  dea l i ngwith error s in the  QMC data :  thed at a are  se rially  c orr elate d andt heunder lyingpr obabil ity d is tributio ns are non-Gaussian. Th e prob abili tydistribut ion ofthe local  ene rgies has$|E -E_ 0|^{- 4}$ tails , wh e re  $E _ 0$ is  a c o nstant. T h es e t a i ls  arise from s i ngu larit ies in the  loc al energy such as the divergence  att h e n oda l  sur fa ce [@trail_200 5_1 ;@ t rail_200 5_ 2], as show n in fig ur e  \[fi g:loca l\_ene rgy\].I n  c o nseque nce,  al though th e m ea n  energy  a nd its va rian ce  are w ell de f ined , the variance ofthe v a r iance isinfin it y. Foro ther  quantitie s the probl em may  beevenmore se ve re; fo r e xa mple, thep robabilit y dis tributi on s fo r t he Pul 
 for_different time_steps and extrapolating to_zero time_step_or by_simply_choosing a small_enough time step._An example of time-step_extrapolation is shown_in_figure \[fig:time\_step\_errors\].

-   Population control bias. The $f$ distribution is represented by a_finite_population of_configurations_which_fluctuates due to branching. The_population may be controlled in_various ways,_but this introduces a population control bias which_is_positive and falls_off as the reciprocal of the population. In practice_the population control bias is normally_so small that_it_is_difficult to detect [@umrigar_1993;_@drummond_2004].

-   Finite size errors_within periodic boundary conditions calculations. It_is important to correct for finite size_effects carefully, as mentioned in section_\[sec:pbc\].

-   The pseudopotential_approximation inevitably_introduces errors. In DMC there_is an additional_error arising_from the localisation_[@mitas_1991] or semi-localisation [@casula_2006] of the_non-local pseudopotential operator._The localisation error appears to be_quite_small in the_cases_for_which it_has been tested_[@drummond_2006b].

![DMC_energy against_time_step for a 64-electron ferromagnetic 2D_hexagonal_Wigner crystal at density parameter $r_s=50$ a.u. with_a Slater-Jastrow wave function._The_solid line is a_linear fit to the data._[]{data-label="fig:time_step_errors"}](ferro_rs50_N064_exp0002_E_v_dt){width=".7\textwidth"}

Practical methods for handling statistical errors_in QMC_results {#subsec:statistical_errors}
----------------------------------------------------------------

Two main practical problems are encountered when dealing with errors in_the QMC data: the data are_serially correlated and the_underlying probability_distributions_are non-Gaussian. The_probability_distribution of_the local energies has $|E-E_0|^{-4}$ tails, where_$E_0$ is_a constant. These tails arise from_singularities in the local_energy_such as the divergence at the_nodal surface [@trail_2005_1; @trail_2005_2], as shown_in figure \[fig:local\_energy\]. In consequence,_although_the_mean energy and its variance_are well defined, the variance of_the variance is_infinity. For other quantities the problem may_be_even more severe; for example, the_probability_distributions for the Pul