Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

PBC gaussian #706

Merged
merged 31 commits into from
Mar 29, 2018
Merged
Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
31 commits
Select commit Hold shift + click to select a range
cc2db2d
Merge branch 'develop' of https://github.com/QMCPACK/qmcpack into dev…
anbenali Jan 31, 2018
a9c445f
Step 1: All piping Done + Periodic eval of Psi; Remaining: Gradiant! …
anbenali Feb 13, 2018
ecd4280
Workinggit add ../src/QMCWaveFunctions/lcao/SoaAtomicBasisSet.h ../sr…
anbenali Feb 16, 2018
f5bcdce
Working and tested. Need sposet for multidets
anbenali Feb 18, 2018
0ae8873
Tests from Molecular pyscf corrected! All good to go. Testing for sup…
anbenali Feb 20, 2018
d522994
END OF PBC Gaussians. This is the split in the PR
anbenali Feb 21, 2018
4c8d8b1
Should probably be new PR : SPOSET SOA with HDF5 handling for Multid…
anbenali Feb 21, 2018
7baf4ff
Step 2 in sposet. still debugging
anbenali Feb 22, 2018
234511e
Working VMC/DMC for Multideterminants Gaussians with PBC - with H5rea…
anbenali Feb 22, 2018
81e44e2
Fixes for Spherical harmonics
anbenali Mar 14, 2018
f0f11cf
Test case Carbon Diamond 1x1x1
anbenali Mar 22, 2018
4bb8e22
symlinks to reduce number of files
anbenali Mar 22, 2018
d10bfd9
Tests LDA Gaussian complete
anbenali Mar 23, 2018
30199e2
Forgotten file
anbenali Mar 23, 2018
3d8a1e7
Adding endpoint implementation of NbImages to be used in the evaluati…
anbenali Mar 25, 2018
903842d
Remove old commented implemnatation
anbenali Mar 25, 2018
0c34c3b
FULLY COMPLETED: Nb of Periodic images to evaluate orbitals effectiv…
anbenali Mar 25, 2018
d4e8bc6
Merge branch 'develop' into PBC-Gaussian
anbenali Mar 25, 2018
0e18560
Forgotten to re-enable when kpoints are defined as Complex in pyscf b…
anbenali Mar 25, 2018
9eb35e5
I do not know why this is here
anbenali Mar 25, 2018
152d396
Fixed qmc-ref
anbenali Mar 26, 2018
b10a92f
Removing test for kpoints since not completed yet
anbenali Mar 26, 2018
2d03c2a
bug in CMakeList
anbenali Mar 26, 2018
c263642
Fixed molecular H5 files to the latest version of builder
anbenali Mar 26, 2018
47ded4a
Fixing more failed tests in Molecules
anbenali Mar 26, 2018
c306573
Fixing all H5 files in Molecules
anbenali Mar 26, 2018
b444e24
more cleanup
anbenali Mar 26, 2018
f3878f8
Reference for diamond 1x1x1
anbenali Mar 29, 2018
e819f3a
All tests + minor fixes
anbenali Mar 29, 2018
292bb3b
Final Push. All requests answered
anbenali Mar 29, 2018
a0e8aa3
Merge branch 'develop' into PBC-Gaussian
anbenali Mar 29, 2018
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
Tests LDA Gaussian complete
  • Loading branch information
anbenali committed Mar 23, 2018
commit d10bfd964e2c812ce4d13603d706cd4454443bf1
8 changes: 4 additions & 4 deletions tests/solids/diamondC_1x1x1-Gaussian_pp/CMakeLists.txt
Original file line number Diff line number Diff line change
Expand Up @@ -13,7 +13,7 @@ IF(ENABLE_SOA)
LIST(APPEND DIAMOND_SCALARS "samples" "128000 0.0")
# LIST(APPEND DIAMOND_SCALARS "flux" "0.0 0.4")

QMC_RUN_AND_CHECK(short-diamondC_1x1x1_pp-vmc_sdj
QMC_RUN_AND_CHECK(short-diamondC_1x1x1_pp-vmc_gausian_sdj
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

fix the test name typo gausian to gaussian

"${CMAKE_SOURCE_DIR}/tests/solids/diamondC_1x1x1-Gaussian_pp"
qmc_short
qmc_short.in.xml
Expand All @@ -25,7 +25,7 @@ IF(ENABLE_SOA)

# Reference DMC run in qmc-ref "-10.541 0.031"
LIST(APPEND DIAMOND_DMC_SCALARS "totenergy" "-10.541 0.031")
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Please put consistent energy value.

QMC_RUN_AND_CHECK(short-diamondC_1x1x1_pp-dmc_sdj
QMC_RUN_AND_CHECK(short-diamondC_1x1x1_pp-dmc_gaussian_sdj
"${CMAKE_SOURCE_DIR}/tests/solids/diamondC_1x1x1-Gaussian_pp"
qmc_short_vmc_dmc
qmc_short_vmc_dmc.in.xml
Expand All @@ -42,7 +42,7 @@ IF(ENABLE_SOA)
LIST(APPEND LONG_DIAMOND_SCALARS "totenergy" "-10.4951 0.0017")
LIST(APPEND LONG_DIAMOND_SCALARS "samples" "192640 0.0")

QMC_RUN_AND_CHECK(long-diamondC_1x1x1_pp-vmc_sdj
QMC_RUN_AND_CHECK(long-diamondC_1x1x1_pp-vmc_gaussian_sdj
"${CMAKE_SOURCE_DIR}/tests/solids/diamondC_1x1x1-Gaussian_pp"
qmc_long
qmc_long.in.xml
Expand All @@ -54,7 +54,7 @@ IF(ENABLE_SOA)

# Reference DMC run in qmc-ref "-10.544774 0.014660"
LIST(APPEND LONG_DIAMOND_DMC_SCALARS "totenergy" "-10.544774 0.014408")
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Put consistent value.

QMC_RUN_AND_CHECK(long-diamondC_1x1x1_pp-dmc_sdj
QMC_RUN_AND_CHECK(long-diamondC_1x1x1_pp-dmc_gaussian_sdj
"${CMAKE_SOURCE_DIR}/tests/solids/diamondC_1x1x1-Gaussian_pp"
qmc_long_vmc_dmc
qmc_long_vmc_dmc.in.xml
Expand Down
Binary file not shown.
1 change: 1 addition & 0 deletions tests/solids/diamondC_1x1x1-Gaussian_pp/C_Diamond.h5
Original file line number Diff line number Diff line change
@@ -1,13 +1,13 @@
<?xml version="1.0"?>
<simulation>
<project id="qmc_long" series="0"/>
<include href="../C_Diamond.structure.xml"/>
<include href="../C_Diamond.wfj-Twist0.xml"/>
<include href="C_Diamond.structure.xml"/>
<include href="C_Diamond.wfj-Twist0.xml"/>
<hamiltonian name="h0" type="generic" target="e">
<pairpot name="ElecElec" type="coulomb" source="e" target="e" physical="true"/>
<pairpot name="IonIon" type="coulomb" source="ion0" target="ion0"/>
<pairpot name="PseudoPot" type="pseudo" source="ion0" wavefunction="psi0" format="xml">
<pseudo elementType="C" href="../C.BFD.xml"/>
<pseudo elementType="C" href="C.BFD.xml"/>
</pairpot>
</hamiltonian>
<qmc method="vmc" move="pbyp">
Expand Down
Original file line number Diff line number Diff line change
@@ -1,13 +1,13 @@
<?xml version="1.0"?>
<simulation>
<project id="qmc_long_vmc_dmc" series="0"/>
<include href="../C_Diamond.structure.xml"/>
<include href="../C_Diamond.wfj-Twist0.xml"/>
<include href="C_Diamond.structure.xml"/>
<include href="C_Diamond.wfj-Twist0.xml"/>
<hamiltonian name="h0" type="generic" target="e">
<pairpot name="ElecElec" type="coulomb" source="e" target="e" physical="true"/>
<pairpot name="IonIon" type="coulomb" source="ion0" target="ion0"/>
<pairpot name="PseudoPot" type="pseudo" source="ion0" wavefunction="psi0" format="xml">
<pseudo elementType="C" href="../C.BFD.xml"/>
<pseudo elementType="C" href="C.BFD.xml"/>
</pairpot>
</hamiltonian>
<qmc method="vmc" move="pbyp">
Expand Down
1 change: 1 addition & 0 deletions tests/solids/diamondC_2x1x1-Gaussian_pp/C.BFD.xml
66 changes: 66 additions & 0 deletions tests/solids/diamondC_2x1x1-Gaussian_pp/CMakeLists.txt
Original file line number Diff line number Diff line change
@@ -0,0 +1,66 @@


IF(ENABLE_SOA)
LIST(APPEND DIAMOND_SCALARS "totenergy" "-21.702829 0.017")
LIST(APPEND DIAMOND_SCALARS "kinetic" "20.706 0.067")
LIST(APPEND DIAMOND_SCALARS "potential" "-42.409 0.067")
LIST(APPEND DIAMOND_SCALARS "eeenergy" "-5.361 0.021")
LIST(APPEND DIAMOND_SCALARS "ionion" "-25.55 0.00")
LIST(APPEND DIAMOND_SCALARS "localecp" "-13.244 0.079")
LIST(APPEND DIAMOND_SCALARS "nonlocalecp" "1.748 0.045")
LIST(APPEND DIAMOND_SCALARS "samples" "6880 0.0")
# LIST(APPEND DIAMOND_SCALARS "flux" "0.0 0.4")

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Update all the values according to your ref run.

QMC_RUN_AND_CHECK(short-diamondC_2x1x1_pp-vmc_gaussian_sdj
"${CMAKE_SOURCE_DIR}/tests/solids/diamondC_2x1x1-Gaussian_pp"
qmc_short
qmc_short.in.xml
1 16
TRUE
0 DIAMOND_SCALARS # VMC
)


# Reference DMC run in qmc-ref "-21.807657 0.095185"
LIST(APPEND DIAMOND_DMC_SCALARS "totenergy" "-21.807657 0.095185")
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Put consistent value

QMC_RUN_AND_CHECK(short-diamondC_2x1x1_pp-dmc_gaussian_sdj
"${CMAKE_SOURCE_DIR}/tests/solids/diamondC_2x1x1-Gaussian_pp"
qmc_short_vmc_dmc
qmc_short_vmc_dmc.in.xml
1 16
TRUE
1 DIAMOND_DMC_SCALARS # DMC
)


#
# Long tests
#

LIST(APPEND LONG_DIAMOND_SCALARS "totenergy" "-21.709193 0.006340")
LIST(APPEND LONG_DIAMOND_SCALARS "samples" "49600 0.0")

QMC_RUN_AND_CHECK(long-diamondC_2x1x1_pp-vmc_gaussian_sdj
"${CMAKE_SOURCE_DIR}/tests/solids/diamondC_2x1x1-Gaussian_pp"
qmc_long
qmc_long.in.xml
1 16
TRUE
0 LONG_DIAMOND_SCALARS # VMC
)


# Reference DMC run in qmc-ref "-21.840087 +/- 0.026449"
LIST(APPEND LONG_DIAMOND_DMC_SCALARS "totenergy" "-21.840087 0.026449")
QMC_RUN_AND_CHECK(long-diamondC_2x1x1_pp-dmc_gaussian_sdj
"${CMAKE_SOURCE_DIR}/tests/solids/diamondC_2x1x1-Gaussian_pp"
qmc_long_vmc_dmc
qmc_long_vmc_dmc.in.xml
1 16
TRUE
1 LONG_DIAMOND_DMC_SCALARS # DMC
)

ELSE()
MESSAGE("Skipping all tests since PBC gaussians only enabled with SOA")
ENDIF()
1 change: 1 addition & 0 deletions tests/solids/diamondC_2x1x1-Gaussian_pp/C_Diamond.h5
36 changes: 36 additions & 0 deletions tests/solids/diamondC_2x1x1-Gaussian_pp/C_Diamond.structure.xml
Original file line number Diff line number Diff line change
@@ -0,0 +1,36 @@
<?xml version="1.0"?>
<qmcsystem>
<simulationcell>
<parameter name="lattice">
6.74632230000000e+00 6.74632230000000e+00 0.00000000000000e+00
0.00000000000000e+00 3.37316115000000e+00 3.37316115000000e+00
3.37316115000000e+00 0.00000000000000e+00 3.37316115000000e+00
</parameter>
<parameter name="bconds">p p p</parameter>
<parameter name="LR_dim_cutoff">15</parameter>
</simulationcell>
<particleset name="ion0" size="4">
<group name="C">
<parameter name="charge">4</parameter>
<parameter name="valence">4</parameter>
<parameter name="atomicnumber">6</parameter>
</group>
<attrib name="position" datatype="posArray">
0.0000000000e+00 0.0000000000e+00 0.0000000000e+00
1.6865805750e+00 1.6865805750e+00 1.6865805750e+00
3.3731611500e+00 3.3731611500e+00 0.0000000000e+00
5.0597417250e+00 5.0597417250e+00 1.6865805750e+00
</attrib>
<attrib name="ionid" datatype="stringArray">
C C C C
</attrib>
</particleset>
<particleset name="e" random="yes" randomsrc="ion0">
<group name="u" size="8">
<parameter name="charge">-1</parameter>
</group>
<group name="d" size="8">
<parameter name="charge">-1</parameter>
</group>
</particleset>
</qmcsystem>
30 changes: 30 additions & 0 deletions tests/solids/diamondC_2x1x1-Gaussian_pp/C_Diamond.wfj-Twist0.xml
Original file line number Diff line number Diff line change
@@ -0,0 +1,30 @@
<?xml version="1.0"?>
<qmcsystem>
<wavefunction name="psi0" target="e">
<determinantset type="MolecularOrbital" name="LCAOBSet" source="ion0" transform="yes" twist="0 0 0" href="../C_Diamond.h5" PBCimages="5 5 5">
<slaterdeterminant>
<determinant id="updet" size="8">
<occupation mode="ground"/>
<coefficient size="116" spindataset="0"/>
</determinant>
<determinant id="downdet" size="8">
<occupation mode="ground"/>
<coefficient size="116" spindataset="0"/>
</determinant>
</slaterdeterminant>
</determinantset>
<jastrow name="J2" type="Two-Body" function="Bspline" print="yes">
<correlation size="10" speciesA="u" speciesB="u">
<coefficients id="uu" type="Array"> 0.3019254237 0.2517951909 0.2061570882 0.1641736643 0.1269421094 0.09561604135 0.06839479985 0.04632683097 0.02721411268 0.01465252</coefficients>
</correlation>
<correlation size="10" speciesA="u" speciesB="d">
<coefficients id="ud" type="Array"> 0.5022561801 0.4048443902 0.3216075783 0.2492321374 0.1884420469 0.1383293395 0.09666256943 0.0639067195 0.03677653401 0.01946279858</coefficients>
</correlation>
</jastrow>
<jastrow name="J1" type="One-Body" function="Bspline" source="ion0" print="yes">
<correlation size="10" cusp="0" elementType="C">
<coefficients id="eC" type="Array"> -0.2418983031 -0.2330046528 -0.2188560081 -0.1998074922 -0.1762538476 -0.1457518169 -0.1132881432 -0.08069267659 -0.04840654078 -0.0254674871</coefficients>
</correlation>
</jastrow>
</wavefunction>
</qmcsystem>
Binary file not shown.
Loading