TREX/qmc-lttc
1
0
Fork 0
mirror of https://github.com/TREX-CoE/qmc-lttc.git synced 2025-01-03 01:56:20 +01:00
Code Issues Projects Releases Wiki Activity

deploy: a7953fe007

Browse Source
This commit is contained in:
scemama 2021-01-13 17:26:25 +00:00
parent 9e9b0a5e50
commit d9c4ffc855
3 changed files with 84 additions and 84 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

168
index.html
View File

@ -3,7 +3,7 @@
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"> "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" lang="en" xml:lang="en"> <html xmlns="http://www.w3.org/1999/xhtml" lang="en" xml:lang="en">
<head> <head>
<!-- 2021-01-13 Wed 17:25 --> <!-- 2021-01-13 Wed 17:26 -->
<meta http-equiv="Content-Type" content="text/html;charset=utf-8" /> <meta http-equiv="Content-Type" content="text/html;charset=utf-8" />
<meta name="viewport" content="width=device-width, initial-scale=1" /> <meta name="viewport" content="width=device-width, initial-scale=1" />
<title>Quantum Monte Carlo</title> <title>Quantum Monte Carlo</title>
@ -257,63 +257,63 @@ for the JavaScript code in this tag.
<h2>Table of Contents</h2> <h2>Table of Contents</h2>
<div id="text-table-of-contents"> <div id="text-table-of-contents">
<ul> <ul>
<li><a href="#org3f6ec9c">1. Introduction</a></li> <li><a href="#org42ed47a">1. Introduction</a></li>
<li><a href="#orgdbe3e57">2. Numerical evaluation of the energy</a> <li><a href="#org53e00dc">2. Numerical evaluation of the energy</a>
<ul> <ul>
<li><a href="#org38b45b3">2.1. Local energy</a> <li><a href="#orgcb06e2a">2.1. Local energy</a>
<ul> <ul>
<li><a href="#org53887c4">2.1.1. Exercise 1</a></li> <li><a href="#org0e5d0fc">2.1.1. Exercise 1</a></li>
<li><a href="#org72355fc">2.1.2. Exercise 2</a></li> <li><a href="#org9b94717">2.1.2. Exercise 2</a></li>
<li><a href="#org4ba063d">2.1.3. Exercise 3</a></li> <li><a href="#org2bcd29b">2.1.3. Exercise 3</a></li>
<li><a href="#org0a0f1ab">2.1.4. Exercise 4</a></li> <li><a href="#org1203f81">2.1.4. Exercise 4</a></li>
</ul> </ul>
</li> </li>
<li><a href="#org3e462eb">2.2. Plot of the local energy along the \(x\) axis</a> <li><a href="#org5d34261">2.2. Plot of the local energy along the \(x\) axis</a>
<ul> <ul>
<li><a href="#orga0b3c58">2.2.1. Exercise</a></li> <li><a href="#org3cb394d">2.2.1. Exercise</a></li>
</ul> </ul>
</li> </li>
<li><a href="#org4b33aa5">2.3. Numerical estimation of the energy</a> <li><a href="#org35b332e">2.3. Numerical estimation of the energy</a>
<ul> <ul>
<li><a href="#org91bfa39">2.3.1. Exercise</a></li> <li><a href="#org7927f35">2.3.1. Exercise</a></li>
</ul> </ul>
</li> </li>
<li><a href="#orga5c38b2">2.4. Variance of the local energy</a> <li><a href="#org565a2b8">2.4. Variance of the local energy</a>
<ul> <ul>
<li><a href="#orgf94cf34">2.4.1. Exercise</a></li> <li><a href="#org773f3b3">2.4.1. Exercise</a></li>
</ul> </ul>
</li> </li>
</ul> </ul>
</li> </li>
<li><a href="#orgebe9860">3. Variational Monte Carlo</a> <li><a href="#org55d88f8">3. Variational Monte Carlo</a>
<ul> <ul>
<li><a href="#org96b3bb6">3.1. Computation of the statistical error</a> <li><a href="#orgd8d3851">3.1. Computation of the statistical error</a>
<ul> <ul>
<li><a href="#orgbcb7a27">3.1.1. Exercise</a></li> <li><a href="#org96e37d8">3.1.1. Exercise</a></li>
</ul> </ul>
</li> </li>
<li><a href="#org86c0ce5">3.2. Uniform sampling in the box</a> <li><a href="#org406e2b4">3.2. Uniform sampling in the box</a>
<ul> <ul>
<li><a href="#org35ee5a9">3.2.1. Exercise</a></li> <li><a href="#org15e2594">3.2.1. Exercise</a></li>
</ul> </ul>
</li> </li>
<li><a href="#org85bd312">3.3. Gaussian sampling</a> <li><a href="#org9158965">3.3. Gaussian sampling</a>
<ul> <ul>
<li><a href="#orgaae11fe">3.3.1. Exercise</a></li> <li><a href="#org4488162">3.3.1. Exercise</a></li>
</ul> </ul>
</li> </li>
<li><a href="#org3452f49">3.4. Sampling with \(\Psi^2\)</a> <li><a href="#org451cd23">3.4. Sampling with \(\Psi^2\)</a>
<ul> <ul>
<li><a href="#org1afdd0f">3.4.1. Importance sampling</a></li> <li><a href="#orgb72825c">3.4.1. Importance sampling</a></li>
<li><a href="#org52c8984">3.4.2. Metropolis algorithm</a></li> <li><a href="#orgff63624">3.4.2. Metropolis algorithm</a></li>
</ul> </ul>
</li> </li>
</ul> </ul>
</li> </li>
<li><a href="#org8958dc7">4. <span class="todo TODO">TODO</span> Diffusion Monte Carlo</a> <li><a href="#org075fffc">4. <span class="todo TODO">TODO</span> Diffusion Monte Carlo</a>
<ul> <ul>
<li><a href="#orgea4f6d6">4.1. Hydrogen atom</a></li> <li><a href="#orgd2b43b9">4.1. Hydrogen atom</a></li>
<li><a href="#org2f5f17e">4.2. Dihydrogen</a></li> <li><a href="#orgc22643f">4.2. Dihydrogen</a></li>
</ul> </ul>
</li> </li>
</ul> </ul>
@ -321,8 +321,8 @@ for the JavaScript code in this tag.
</div> </div>
<div id="outline-container-org3f6ec9c" class="outline-2"> <div id="outline-container-org42ed47a" class="outline-2">
<h2 id="org3f6ec9c"><span class="section-number-2">1</span> Introduction</h2> <h2 id="org42ed47a"><span class="section-number-2">1</span> Introduction</h2>
<div class="outline-text-2" id="text-1"> <div class="outline-text-2" id="text-1">
<p> <p>
We propose different exercises to understand quantum Monte Carlo (QMC) We propose different exercises to understand quantum Monte Carlo (QMC)
@ -364,8 +364,8 @@ interpreted as a single precision value
</div> </div>
<div id="outline-container-orgdbe3e57" class="outline-2"> <div id="outline-container-org53e00dc" class="outline-2">
<h2 id="orgdbe3e57"><span class="section-number-2">2</span> Numerical evaluation of the energy</h2> <h2 id="org53e00dc"><span class="section-number-2">2</span> Numerical evaluation of the energy</h2>
<div class="outline-text-2" id="text-2"> <div class="outline-text-2" id="text-2">
<p> <p>
In this section we consider the Hydrogen atom with the following In this section we consider the Hydrogen atom with the following
@ -439,13 +439,13 @@ E & = & \frac{\langle \Psi| \hat{H} | \Psi\rangle}{\langle \Psi |\Psi \rangle}
\end{eqnarray*} \end{eqnarray*}
</div> </div>
<div id="outline-container-org38b45b3" class="outline-3"> <div id="outline-container-orgcb06e2a" class="outline-3">
<h3 id="org38b45b3"><span class="section-number-3">2.1</span> Local energy</h3> <h3 id="orgcb06e2a"><span class="section-number-3">2.1</span> Local energy</h3>
<div class="outline-text-3" id="text-2-1"> <div class="outline-text-3" id="text-2-1">
</div> </div>
<div id="outline-container-org53887c4" class="outline-4"> <div id="outline-container-org0e5d0fc" class="outline-4">
<h4 id="org53887c4"><span class="section-number-4">2.1.1</span> Exercise 1</h4> <h4 id="org0e5d0fc"><span class="section-number-4">2.1.1</span> Exercise 1</h4>
<div class="outline-text-4" id="text-2-1-1"> <div class="outline-text-4" id="text-2-1-1">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -489,8 +489,8 @@ and returns the potential.
</div> </div>
</div> </div>
<div id="outline-container-org72355fc" class="outline-4"> <div id="outline-container-org9b94717" class="outline-4">
<h4 id="org72355fc"><span class="section-number-4">2.1.2</span> Exercise 2</h4> <h4 id="org9b94717"><span class="section-number-4">2.1.2</span> Exercise 2</h4>
<div class="outline-text-4" id="text-2-1-2"> <div class="outline-text-4" id="text-2-1-2">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -525,8 +525,8 @@ input arguments, and returns a scalar.
</div> </div>
</div> </div>
<div id="outline-container-org4ba063d" class="outline-4"> <div id="outline-container-org2bcd29b" class="outline-4">
<h4 id="org4ba063d"><span class="section-number-4">2.1.3</span> Exercise 3</h4> <h4 id="org2bcd29b"><span class="section-number-4">2.1.3</span> Exercise 3</h4>
<div class="outline-text-4" id="text-2-1-3"> <div class="outline-text-4" id="text-2-1-3">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -607,8 +607,8 @@ So the local kinetic energy is
</div> </div>
</div> </div>
<div id="outline-container-org0a0f1ab" class="outline-4"> <div id="outline-container-org1203f81" class="outline-4">
<h4 id="org0a0f1ab"><span class="section-number-4">2.1.4</span> Exercise 4</h4> <h4 id="org1203f81"><span class="section-number-4">2.1.4</span> Exercise 4</h4>
<div class="outline-text-4" id="text-2-1-4"> <div class="outline-text-4" id="text-2-1-4">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -651,14 +651,14 @@ local energy.
</div> </div>
</div> </div>
<div id="outline-container-org3e462eb" class="outline-3"> <div id="outline-container-org5d34261" class="outline-3">
<h3 id="org3e462eb"><span class="section-number-3">2.2</span> Plot of the local energy along the \(x\) axis</h3> <h3 id="org5d34261"><span class="section-number-3">2.2</span> Plot of the local energy along the \(x\) axis</h3>
<div class="outline-text-3" id="text-2-2"> <div class="outline-text-3" id="text-2-2">
</div> </div>
<div id="outline-container-orga0b3c58" class="outline-4"> <div id="outline-container-org3cb394d" class="outline-4">
<h4 id="orga0b3c58"><span class="section-number-4">2.2.1</span> Exercise</h4> <h4 id="org3cb394d"><span class="section-number-4">2.2.1</span> Exercise</h4>
<div class="outline-text-4" id="text-2-2-1"> <div class="outline-text-4" id="text-2-2-1">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -775,8 +775,8 @@ plot './data' index 0 using 1:2 with lines title 'a=0.1', \
</div> </div>
</div> </div>
<div id="outline-container-org4b33aa5" class="outline-3"> <div id="outline-container-org35b332e" class="outline-3">
<h3 id="org4b33aa5"><span class="section-number-3">2.3</span> Numerical estimation of the energy</h3> <h3 id="org35b332e"><span class="section-number-3">2.3</span> Numerical estimation of the energy</h3>
<div class="outline-text-3" id="text-2-3"> <div class="outline-text-3" id="text-2-3">
<p> <p>
If the space is discretized in small volume elements \(\mathbf{r}_i\) If the space is discretized in small volume elements \(\mathbf{r}_i\)
@ -806,8 +806,8 @@ The energy is biased because:
</div> </div>
<div id="outline-container-org91bfa39" class="outline-4"> <div id="outline-container-org7927f35" class="outline-4">
<h4 id="org91bfa39"><span class="section-number-4">2.3.1</span> Exercise</h4> <h4 id="org7927f35"><span class="section-number-4">2.3.1</span> Exercise</h4>
<div class="outline-text-4" id="text-2-3-1"> <div class="outline-text-4" id="text-2-3-1">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -917,8 +917,8 @@ a = 2.0000000000000000 E = -8.0869806678448772E-002
</div> </div>
</div> </div>
<div id="outline-container-orga5c38b2" class="outline-3"> <div id="outline-container-org565a2b8" class="outline-3">
<h3 id="orga5c38b2"><span class="section-number-3">2.4</span> Variance of the local energy</h3> <h3 id="org565a2b8"><span class="section-number-3">2.4</span> Variance of the local energy</h3>
<div class="outline-text-3" id="text-2-4"> <div class="outline-text-3" id="text-2-4">
<p> <p>
The variance of the local energy is a functional of \(\Psi\) The variance of the local energy is a functional of \(\Psi\)
@ -940,8 +940,8 @@ energy can be used as a measure of the quality of a wave function.
</p> </p>
</div> </div>
<div id="outline-container-orgf94cf34" class="outline-4"> <div id="outline-container-org773f3b3" class="outline-4">
<h4 id="orgf94cf34"><span class="section-number-4">2.4.1</span> Exercise</h4> <h4 id="org773f3b3"><span class="section-number-4">2.4.1</span> Exercise</h4>
<div class="outline-text-4" id="text-2-4-1"> <div class="outline-text-4" id="text-2-4-1">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -1083,8 +1083,8 @@ a = 2.0000000000000000 E = -8.0869806678448772E-002 s2 = 1.806881
</div> </div>
<div id="outline-container-orgebe9860" class="outline-2"> <div id="outline-container-org55d88f8" class="outline-2">
<h2 id="orgebe9860"><span class="section-number-2">3</span> Variational Monte Carlo</h2> <h2 id="org55d88f8"><span class="section-number-2">3</span> Variational Monte Carlo</h2>
<div class="outline-text-2" id="text-3"> <div class="outline-text-2" id="text-3">
<p> <p>
Numerical integration with deterministic methods is very efficient Numerical integration with deterministic methods is very efficient
@ -1100,8 +1100,8 @@ interval.
</p> </p>
</div> </div>
<div id="outline-container-org96b3bb6" class="outline-3"> <div id="outline-container-orgd8d3851" class="outline-3">
<h3 id="org96b3bb6"><span class="section-number-3">3.1</span> Computation of the statistical error</h3> <h3 id="orgd8d3851"><span class="section-number-3">3.1</span> Computation of the statistical error</h3>
<div class="outline-text-3" id="text-3-1"> <div class="outline-text-3" id="text-3-1">
<p> <p>
To compute the statistical error, you need to perform \(M\) To compute the statistical error, you need to perform \(M\)
@ -1141,8 +1141,8 @@ And the confidence interval is given by
</p> </p>
</div> </div>
<div id="outline-container-orgbcb7a27" class="outline-4"> <div id="outline-container-org96e37d8" class="outline-4">
<h4 id="orgbcb7a27"><span class="section-number-4">3.1.1</span> Exercise</h4> <h4 id="org96e37d8"><span class="section-number-4">3.1.1</span> Exercise</h4>
<div class="outline-text-4" id="text-3-1-1"> <div class="outline-text-4" id="text-3-1-1">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -1191,8 +1191,8 @@ input array.
</div> </div>
</div> </div>
<div id="outline-container-org86c0ce5" class="outline-3"> <div id="outline-container-org406e2b4" class="outline-3">
<h3 id="org86c0ce5"><span class="section-number-3">3.2</span> Uniform sampling in the box</h3> <h3 id="org406e2b4"><span class="section-number-3">3.2</span> Uniform sampling in the box</h3>
<div class="outline-text-3" id="text-3-2"> <div class="outline-text-3" id="text-3-2">
<p> <p>
We will now do our first Monte Carlo calculation to compute the We will now do our first Monte Carlo calculation to compute the
@ -1226,8 +1226,8 @@ statistical error.
</p> </p>
</div> </div>
<div id="outline-container-org35ee5a9" class="outline-4"> <div id="outline-container-org15e2594" class="outline-4">
<h4 id="org35ee5a9"><span class="section-number-4">3.2.1</span> Exercise</h4> <h4 id="org15e2594"><span class="section-number-4">3.2.1</span> Exercise</h4>
<div class="outline-text-4" id="text-3-2-1"> <div class="outline-text-4" id="text-3-2-1">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -1337,8 +1337,8 @@ E = -0.49588321986667677 +/- 7.1758863546737969E-004
</div> </div>
</div> </div>
<div id="outline-container-org85bd312" class="outline-3"> <div id="outline-container-org9158965" class="outline-3">
<h3 id="org85bd312"><span class="section-number-3">3.3</span> Gaussian sampling</h3> <h3 id="org9158965"><span class="section-number-3">3.3</span> Gaussian sampling</h3>
<div class="outline-text-3" id="text-3-3"> <div class="outline-text-3" id="text-3-3">
<p> <p>
We will now improve the sampling and allow to sample in the whole We will now improve the sampling and allow to sample in the whole
@ -1434,8 +1434,8 @@ average energy can be computed as
</div> </div>
<div id="outline-container-orgaae11fe" class="outline-4"> <div id="outline-container-org4488162" class="outline-4">
<h4 id="orgaae11fe"><span class="section-number-4">3.3.1</span> Exercise</h4> <h4 id="org4488162"><span class="section-number-4">3.3.1</span> Exercise</h4>
<div class="outline-text-4" id="text-3-3-1"> <div class="outline-text-4" id="text-3-3-1">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -1546,8 +1546,8 @@ E = -0.49517104619091717 +/- 1.0685523607878961E-004
</div> </div>
</div> </div>
<div id="outline-container-org3452f49" class="outline-3"> <div id="outline-container-org451cd23" class="outline-3">
<h3 id="org3452f49"><span class="section-number-3">3.4</span> Sampling with \(\Psi^2\)</h3> <h3 id="org451cd23"><span class="section-number-3">3.4</span> Sampling with \(\Psi^2\)</h3>
<div class="outline-text-3" id="text-3-4"> <div class="outline-text-3" id="text-3-4">
<p> <p>
We will now use the square of the wave function to make the sampling: We will now use the square of the wave function to make the sampling:
@ -1572,8 +1572,8 @@ the local energies, each with a weight of 1.
</div> </div>
<div id="outline-container-org1afdd0f" class="outline-4"> <div id="outline-container-orgb72825c" class="outline-4">
<h4 id="org1afdd0f"><span class="section-number-4">3.4.1</span> Importance sampling</h4> <h4 id="orgb72825c"><span class="section-number-4">3.4.1</span> Importance sampling</h4>
<div class="outline-text-4" id="text-3-4-1"> <div class="outline-text-4" id="text-3-4-1">
<p> <p>
To generate the probability density \(\Psi^2\), we consider a To generate the probability density \(\Psi^2\), we consider a
@ -1686,7 +1686,7 @@ variance \(\tau\,2D\).
</div> </div>
<ol class="org-ol"> <ol class="org-ol">
<li><a id="org579ab16"></a>Exercise 1<br /> <li><a id="org8061073"></a>Exercise 1<br />
<div class="outline-text-5" id="text-3-4-1-1"> <div class="outline-text-5" id="text-3-4-1-1">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -1722,7 +1722,7 @@ Write a function to compute the drift vector \(\frac{\nabla \Psi(\mathbf{r})}{\P
</div> </div>
</li> </li>
<li><a id="org329a908"></a><span class="todo TODO">TODO</span> Exercise 2<br /> <li><a id="org22449f0"></a><span class="todo TODO">TODO</span> Exercise 2<br />
<div class="outline-text-5" id="text-3-4-1-2"> <div class="outline-text-5" id="text-3-4-1-2">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -1834,8 +1834,8 @@ E = -0.48584030499187431 +/- 1.0411743995438257E-004
</ol> </ol>
</div> </div>
<div id="outline-container-org52c8984" class="outline-4"> <div id="outline-container-orgff63624" class="outline-4">
<h4 id="org52c8984"><span class="section-number-4">3.4.2</span> Metropolis algorithm</h4> <h4 id="orgff63624"><span class="section-number-4">3.4.2</span> Metropolis algorithm</h4>
<div class="outline-text-4" id="text-3-4-2"> <div class="outline-text-4" id="text-3-4-2">
<p> <p>
Discretizing the differential equation to generate the desired Discretizing the differential equation to generate the desired
@ -1896,7 +1896,7 @@ the simulation.
</div> </div>
<ol class="org-ol"> <ol class="org-ol">
<li><a id="org9e1ca89"></a>Exercise<br /> <li><a id="orgca456a4"></a>Exercise<br />
<div class="outline-text-5" id="text-3-4-2-1"> <div class="outline-text-5" id="text-3-4-2-1">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -2052,17 +2052,17 @@ A = 0.78861366666666655 +/- 3.5096729498002445E-004
</div> </div>
<div id="outline-container-org8958dc7" class="outline-2"> <div id="outline-container-org075fffc" class="outline-2">
<h2 id="org8958dc7"><span class="section-number-2">4</span> <span class="todo TODO">TODO</span> Diffusion Monte Carlo</h2> <h2 id="org075fffc"><span class="section-number-2">4</span> <span class="todo TODO">TODO</span> Diffusion Monte Carlo</h2>
<div class="outline-text-2" id="text-4"> <div class="outline-text-2" id="text-4">
</div> </div>
<div id="outline-container-orgea4f6d6" class="outline-3"> <div id="outline-container-orgd2b43b9" class="outline-3">
<h3 id="orgea4f6d6"><span class="section-number-3">4.1</span> Hydrogen atom</h3> <h3 id="orgd2b43b9"><span class="section-number-3">4.1</span> Hydrogen atom</h3>
<div class="outline-text-3" id="text-4-1"> <div class="outline-text-3" id="text-4-1">
</div> </div>
<ol class="org-ol"> <ol class="org-ol">
<li><a id="orgdfab00a"></a>Exercise<br /> <li><a id="orgf0b78c7"></a>Exercise<br />
<div class="outline-text-5" id="text-4-1-0-1"> <div class="outline-text-5" id="text-4-1-0-1">
<div class="exercise"> <div class="exercise">
<p> <p>
@ -2221,8 +2221,8 @@ A = 0.78861366666666655 +/- 3.5096729498002445E-004
</div> </div>
<div id="outline-container-org2f5f17e" class="outline-3"> <div id="outline-container-orgc22643f" class="outline-3">
<h3 id="org2f5f17e"><span class="section-number-3">4.2</span> Dihydrogen</h3> <h3 id="orgc22643f"><span class="section-number-3">4.2</span> Dihydrogen</h3>
<div class="outline-text-3" id="text-4-2"> <div class="outline-text-3" id="text-4-2">
<p> <p>
We will now consider the H<sub>2</sub> molecule in a minimal basis composed of the We will now consider the H<sub>2</sub> molecule in a minimal basis composed of the
@ -2244,7 +2244,7 @@ the nuclei.
</div> </div>
<div id="postamble" class="status"> <div id="postamble" class="status">
<p class="author">Author: Anthony Scemama, Claudia Filippi</p> <p class="author">Author: Anthony Scemama, Claudia Filippi</p>
<p class="date">Created: 2021-01-13 Wed 17:25</p> <p class="date">Created: 2021-01-13 Wed 17:26</p>
<p class="validation"><a href="http://validator.w3.org/check?uri=referer">Validate</a></p> <p class="validation"><a href="http://validator.w3.org/check?uri=referer">Validate</a></p>
</div> </div>
</body> </body>

BIN
plot.png Normal file
View File

Binary file not shown.
Side by Side

After

Width:  |  Height:  |  Size: 7.8 KiB

BIN
plot_py.png Normal file
View File

Binary file not shown.
Side by Side

After

Width:  |  Height:  |  Size: 40 KiB