, if it is not, the energy is always raised. It follows that the DMC energy is always less than or equal to the VMC energy with the same trial wave function, and always greater than or equal to the exact ground-state energy.

The fixed-node DMC algorithm described above is extremely inefficient and a vastly superior algorithm can be obtained by introducing an importance sampling transformation [@grimm_1971; @kalos_1974]. Consider the mixed distribution, $$\begin{aligned}
\label{eq:f}
f({\bf R},t) = \Psi_{\rm T}({\bf R}) \Phi({\bf R},t) \;,\end{aligned}$$ which has the same sign everywhere if and only if the nodal surface of $\Phi({\bf R},t)$ equals that of $\Psi_{\rm T}({\bf R})$. Substituting in equation (\[eq:imaginary\_time\_se\]) for $\Phi$ we obtain $$\begin{aligned}
\label{eq:importance_sampled_imaginary_time_se}
-\frac{\partial f}{\partial t} = -\frac{1}{2} \nabla_{\bf R}^2 f + \nabla_{\bf
 R} \cdot [{\bf v}f] + [E_{\rm L}-E_{\rm T}]f \;,\end{aligned}$$ where the $3N$-dimensional drift velocity is defined as $$\begin{aligned}
\label{eq:drift_velocity}
{\bf v}({\bf R}) = \Psi_{\rm T}^{-1}({\bf R}) \nabla_{\bf R} \Psi_{\rm T}({\bf
R}) \;.\end{aligned}$$ The three terms on the right-hand side of equation (\[eq:importance\_sampled\_imaginary\_time\_se\]) correspond to diffusion, drift and branching processes, respectively. The importance sampling transformation has several consequences. First, the density of configurations is increased where $|\Psi_{\rm T}|$ is large, so that the more important parts of the wave function are sampled more often. Second, the rate of branching is now controlled by the local energy which is normally a much smoother function than the potential energy. This is particularly important for the Coulomb interaction, which diverges when particles are coincident. The importance sampling transformation, together with an algorithm that imposes $f({\bf R},t)
\ge 0$, ensures that $\Psi_{\rm T}$ and $\Phi({\bf
, if it is not, the energy is always raised. It follows that the DMC department of energy is constantly less than or equal to the VMC energy with the like trial wave routine, and always greater than or adequate to the accurate ground - state energy. 

 The fixed - node DMC algorithm identify above is extremely inefficient and a vastly ranking algorithm can be prevail by introducing an importance sampling transformation [ @grimm_1971; @kalos_1974 ]. Consider the mixed distribution, $ $ \begin{aligned } 
 \label{eq: f } 
 f({\bf R},t) = \Psi_{\rm T}({\bf R }) \Phi({\bf R},t) \;,\end{aligned}$$ which has the like sign everywhere if and only if the nodal surface of $ \Phi({\bf R},t)$ equal that of $ \Psi_{\rm T}({\bf R})$. Substituting in equation (\[eq: imaginary\_time\_se\ ]) for $ \Phi$ we receive $ $ \begin{aligned } 
 \label{eq: importance_sampled_imaginary_time_se } 
 -\frac{\partial f}{\partial t } = -\frac{1}{2 } \nabla_{\bf R}^2 f + \nabla_{\bf 
  roentgen } \cdot [ { \bf v}f ] + [ E_{\rm L}-E_{\rm T}]f \;,\end{aligned}$$ where the $ 3N$-dimensional drift velocity is defined as $ $ \begin{aligned } 
 \label{eq: drift_velocity } 
 { \bf v}({\bf R }) = \Psi_{\rm T}^{-1}({\bf R }) \nabla_{\bf R } \Psi_{\rm T}({\bf 
 R }) \;.\end{aligned}$$ The three terms on the right - bridge player side of equation (\[eq: importance\_sampled\_imaginary\_time\_se\ ]) correspond to diffusion, drift and branching processes, respectively. The importance sampling transformation has several consequences. First, the density of configurations is increased where $ |\Psi_{\rm T}|$ is large, so that the more important parts of the wave function are sampled more frequently. Second, the rate of branching is now master by the local department of energy which is normally a a lot smoother function than the potential energy. This is particularly important for the Coulomb interaction, which diverges when particles are coincident. The importance sample transformation, together with an algorithm that imposes $ f({\bf R},t) 
 \ge 0 $, ensures that $ \Psi_{\rm T}$ and $ \Phi({\bf
, if it is not, the energy is always raised. Nr follmws thzt the DOC energy is always less thai or equao to the VMC energy wigh the sale trial wavt function, and alxzys greater tgwn ox xqual to the exsct ground-vtate energy.

Tha wired-node DMC algorithm described abovq is exyrfmely inefficignt amq a bastly superior algorithm can be ogtained by introducimg an importance sampling hrandformation [@grimm_1971; @kwlos_1974]. Considgd trw mixed distfibution, $$\btgnn{aligned}
\lagel{eq:f}
f({\bf R},t) = \Psi_{\rm T}({\bf R}) \Phi({\bf R},t) \;,\eud{aligned}$$ wyixh jds the same sign everywhere lg and mnly if the nodal surnace mf $\Phi({\bf R},t)$ equals that mf $\Psi_{\rm T}({\bf R})$. Subftituting iu equation (\[eq:imaginart\_tume\_se\]) for $\Phi$ qe ubtzii $$\bsgin{allgnxd}
\label{eq:imlortance_sampled_imaginary_time_se}
-\grwb{\lartial f}{\parfial t} = -\srac{1}{2} \nabla_{\bf R}^2 f + \nabla_{\bf
 R} \cdot [{\bf v}f] + [E_{\rj L}-E_{\rm T}]f \;,\end{aligned}$$ whwre the $3N$-dimensional frift veljcity is defined as $$\begin{aligned}
\label{eq:drift_velochty}
{\bf x}({\bf R}) = \Psl_{\rm T}^{-1}({\bf R}) \nabla_{\bf R} \Psi_{\rm T}({\bf
R}) \;.\end{aligned}$$ The fhtev terms on the rinht-hand side of eqiahipg (\[eq:importancg\_sampleb\_jmzginary\_time\_se\]) corrfspond jo difdusion, dryft snd branching processes, respectively. Thv importance sampling cransformatiun hss seferal consequences. Firsc, the sensity of fonfigurafkons is increasea wmera $|\Psi_{\rm U}|$ is large, so that the more impprtant oartx of tre wave fujctiok are sampled more lften. Sacond, the gate of branching is now controlled by the lovan etergy whnch is normally a mtch smoother fonction tkan thd potentiam energb. This is pawticularly im[lrtant for tie Coulomf inreraxtion, wfkch diverges wnen partiblts are coibcident. The importence aampling transfirmqtion, together eitf ag wljoritrk that imposas $f({\cf F},y)
\ge 0$, dnsures thab $\Pri_{\rm T}$ and $\Phi({\bf
, if it is not, the energy raised. follows that DMC energy is to VMC energy with same trial wave and always greater than or equal the exact ground-state energy. The fixed-node DMC algorithm described above is extremely inefficient a vastly superior algorithm can be obtained by introducing an importance sampling transformation @kalos_1974]. the distribution, \label{eq:f} f({\bf R},t) = \Psi_{\rm T}({\bf R}) \Phi({\bf R},t) \;,\end{aligned}$$ which has the same sign everywhere and only if the nodal surface of $\Phi({\bf equals that of $\Psi_{\rm R})$. Substituting in equation (\[eq:imaginary\_time\_se\]) $\Phi$ obtain $$\begin{aligned} -\frac{\partial t} -\frac{1}{2} \nabla_{\bf R}^2 + \nabla_{\bf R} \cdot [{\bf v}f] + [E_{\rm L}-E_{\rm T}]f \;,\end{aligned}$$ where the $3N$-dimensional drift velocity is as $$\begin{aligned} v}({\bf R}) \Psi_{\rm R}) R} \Psi_{\rm T}({\bf The three terms on the right-hand (\[eq:importance\_sampled\_imaginary\_time\_se\]) correspond to diffusion, drift and branching processes, The importance transformation has several consequences. First, the of configurations is increased where $|\Psi_{\rm T}|$ is so that the more important parts of the wave function are sampled more often. Second, of branching is now by the local which normally much function than potential energy. This is particularly important for the Coulomb interaction, which when particles are coincident. The importance sampling transformation, together with that $f({\bf R},t) \ge ensures that $\Psi_{\rm T}$ $\Phi({\bf
, if it is not, the energy is alwayS raised. It fOllowS thAt tHe dMC eNergY is always less tHAn or Equal to the VMC energy witH the sAmE TriaL WaVe funCtion, anD AlWAYs gReAtEr tHaN Or Equal To tHe exact Ground-statE enErGy.

The fixed-noDE DmC algorithM deScribed above Is eXtremeLy IneFFicieNt aNd a vaStly suPErior aLgorithm cAn BE obtaiNEd by intRODuCing An importance samplINg TRansformation [@gRimm_1971; @kaLoS_1974]. coNSIdeR thE mixed distRiButioN, $$\Begin{alIGnED}
\LAbeL{Eq:f}
f({\bf R},t) = \Psi_{\rm t}({\bf R}) \Phi({\bf R},t) \;,\ENd{aLigned}$$ WhIch HAs the sAme siGn EVerYwhere if and Only If the nodaL surfaCE of $\Phi({\bF r},t)$ equalS that oF $\PsI_{\rm t}({\bf R})$. sUbStItuTiNG in EQuAtiON (\[eq:ImaginarY\_tImE\_se\]) foR $\Phi$ WE OBTain $$\BegIn{alIgned}
\Label{eq:importAncE_samPLed_ImagiNary_tIme_sE}
-\fRac{\paRtial f}{\PartiAl T} = -\frac{1}{2} \nabla_{\bf R}^2 f + \nAbla_{\Bf
 R} \cdot [{\bf V}f] + [E_{\Rm l}-E_{\rM T}]F \;,\end{aLIgned}$$ wHerE thE $3N$-dimenSional dRIft VeLOCItY is defined as $$\begin{aLiGNEd}
\Label{eq:dRift_veLOcItY}
{\Bf v}({\bf R}) = \PsI_{\rM T}^{-1}({\bF R}) \naBLA_{\bf R} \PSi_{\rm t}({\Bf
r}) \;.\end{aligNed}$$ The THrEe Terms on ThE right-HaNd sIde Of equATion (\[Eq:impoRtance\_saMpled\_IMaginary\_time\_se\]) COrrespond to diFFuSIOn, DRift And Branching prOcesSEs, reSpecTIvEly. tHe impOrtanCe SAmPLing transformation hAs SeveraL consEquences. First, The density OF COnfiguraTionS Is INcreased where $|\PSi_{\rm T}|$ Is large, so tHAt the morE impoRtant parTs of the waVE Function Are SamPleD moRE OfTen. Second, the rATE of bRaNching iS noW controLleD by The LocAl Energy whiCh is normAlLy A mUcH smOotheR Function ThAn tHe PotEntiaL Energy. this iS parTiCuLArlY importANt FOR the coUlOmb iNteRaCtion, WhicH DivErges whEn particlEs aRE coiNcIdEnt. The iMportance sampLiNg transforMaTioN, togetHER with an aLgorithm that imposes $f({\bf R},T)
\Ge 0$, ensurEs tHat $\PsI_{\rm T}$ And $\Phi({\bf
, if it is not, the energy  is always  rais ed.  It  f ollo ws t hat the DMC en e rgyis always less than or  equa lt o th e  V MC en ergy wi t ht h e s am etri al wa ve fu nct ion, an d always g rea te r than or eq u al  to the ex act  ground-stat e e nergy. 
The fixed -no de DM C algo r ithm d escribedab o ve ise xtremel y in effi cient and a vastl y  s u perior algorit hm can  b e  o b t ain edby introdu ci ng an importa n ce s a mpl i ng transforma tion [@grim m _19 71; @k al os_ 1 974].Consi de r  th e mixed dis trib ution, $$ \begin { aligned } 
\label {eq:f} 
f( {\b f R} , t)  =  \P si _ {\r m  T }({ \ bfR}) \Phi ({ \b f R}, t) \ ; , \ e nd{a lig ned} $$ wh ich has the s ame  sig n  ev erywh ere i f an donlyif the  noda lsurface of $\Ph i({\ bf R},t)$  eq ua lsth at of $\Psi_ {\r m T }({\bfR})$. S u bst it u t i ng  in equation (\[eq :i m a gi nary\_ti me\_se \ ])  f o r $\Phi$  w e o btai n $$\be gin{ a li gned}
\l abel{e q :i mp ortance _s ampled _i mag ina ry_ti m e_se }
-\fr ac{\part ial f } {\partial t} = -\frac{1}{2}\ na b l a_ { \bfR}^ 2 f + \nabl a_{\ b f
 R } \c d ot  [{ \ bf v} f] +[E _ {\ r m L}-E_{\rm T}]f \; ,\ end{al igned }$$ where the  $3N$-dime n s i onal dri ft v e lo c ity is defined  as $ $\begin{al i gned}
\l abel{ eq:drift _velocity } {\bf v}( {\b f R })= \ P s i_ {\rm T}^{-1}( { \ bf R })  \nabla _{\ bf R} \ Psi _{\ rmT}( {\ bf
R}) \; .\end{al ig ne d} $$  Th e thr e e termson  th erig ht-ha n d side  of e quat io n( \[e q:impor t an c e \_sa mp le d\_i mag in ary\_ time \ _se \]) cor respond t o d i ffus io n,  driftand branching  p rocesses,re spe ctivel y .  The imp ortance sampling transf o rmation  ha s sev eral  conseque nce s. Fir st, the de nsityof co nf igu r a tions i sinc re ased where $ |\P si_{\ rm  T}| $ is la rge, so that the m o reimportant par tsof t h e  w ave fu n cti on are s ampled more oft en. Second ,t he  rate of b r anc hi ng is n ow cont rolle d  by the  local en ergy whic his n o r mal ly a muchsmoother  function thant he  pote nti al ene rg y.Thisis par t icu larly  impor ta nt for  theCo ulomb in teraction, which diverg es whe n par tic les are c oin c ide nt. The i mpor tance samp lin g t ransf orm a tion,  tog e th erw ith a n al g orithm th a timp o s es  $f({\bf R} , t ) 
\g e 0$,  en s ures t hat$\Psi_{\rm T}$ an d  $\Phi({\bf
, if_it is_not, the energy is_always raised._It_follows that_the_DMC energy is_always less than_or equal to the_VMC energy with_the_same trial wave function, and always greater than or equal to the exact ground-state_energy.

The_fixed-node DMC_algorithm_described_above is extremely inefficient and_a vastly superior algorithm can_be obtained_by introducing an importance sampling transformation [@grimm_1971; @kalos_1974]._Consider_the mixed distribution,_$$\begin{aligned}
\label{eq:f}
f({\bf R},t) = \Psi_{\rm T}({\bf R}) \Phi({\bf R},t) \;,\end{aligned}$$_which has the same sign everywhere_if and only_if_the_nodal surface of $\Phi({\bf_R},t)$ equals that of $\Psi_{\rm T}({\bf_R})$. Substituting in equation (\[eq:imaginary\_time\_se\]) for_$\Phi$ we obtain $$\begin{aligned}
\label{eq:importance_sampled_imaginary_time_se}
-\frac{\partial f}{\partial t} =_-\frac{1}{2} \nabla_{\bf R}^2 f + \nabla_{\bf
_R} \cdot [{\bf v}f] +_[E_{\rm L}-E_{\rm_T}]f \;,\end{aligned}$$ where the $3N$-dimensional_drift velocity is_defined as_$$\begin{aligned}
\label{eq:drift_velocity}
{\bf v}({\bf R})_= \Psi_{\rm T}^{-1}({\bf R}) \nabla_{\bf R}_\Psi_{\rm T}({\bf
R}) \;.\end{aligned}$$_The three terms on the right-hand_side_of equation (\[eq:importance\_sampled\_imaginary\_time\_se\])_correspond_to_diffusion, drift_and branching processes,_respectively._The importance_sampling_transformation has several consequences. First, the_density_of configurations is increased where $|\Psi_{\rm T}|$_is large, so that_the_more important parts of_the wave function are sampled_more often. Second, the rate of_branching is_now controlled_by the local energy which is normally a much smoother function_than the potential energy. This is_particularly important for the_Coulomb interaction,_which_diverges when particles_are_coincident. The_importance sampling transformation, together with an algorithm_that imposes_$f({\bf R},t)
\ge 0$, ensures that $\Psi_{\rm_T}$ and $\Phi({\bf