add atom O open BC noj test.

This commit is contained in:
Ye Luo 2017-01-18 00:11:20 -06:00
parent e5885a261a
commit 010731a7e9
3 changed files with 100 additions and 13 deletions

View File

@ -229,17 +229,17 @@ else()
# test case for O atom
# from a long run
# totenergy
# LocalEnergy = -15.870084 +/- 0.000055
# Kinetic = 11.85262 +/- 0.00052
# LocalPotential = -27.72271 +/- 0.00051
# ElecElec = 10.36870 +/- 0.00018
# LocalECP = -40.15674 +/- 0.00071
# NonLocalECP = 2.06533 +/- 0.00030
LIST(APPEND O_ATOM_SCALARS "totenergy" "-15.870084 0.0045")
LIST(APPEND O_ATOM_SCALARS "kinetic" "11.85262 0.04")
LIST(APPEND O_ATOM_SCALARS "localecp" "-40.15674 0.057")
LIST(APPEND O_ATOM_SCALARS "nonlocalecp" "2.06533 0.020")
LIST(APPEND O_ATOM_SCALARS "eeenergy" "10.36870 0.014")
# LocalEnergy = -15.870035 +/- 0.000016
# Kinetic = 11.85233 +/- 0.00014
# ElecElec = 10.368642 +/- 0.000050
# LocalECP = -40.15639 +/- 0.00020
# NonLocalECP = 2.065389 +/- 0.000083
LIST(APPEND O_ATOM_SCALARS "totenergy" "-15.870035 0.0045")
LIST(APPEND O_ATOM_SCALARS "kinetic" "11.85233 0.04")
LIST(APPEND O_ATOM_SCALARS "localecp" "-40.15639 0.057")
LIST(APPEND O_ATOM_SCALARS "nonlocalecp" "2.065389 0.020")
LIST(APPEND O_ATOM_SCALARS "eeenergy" "10.368642 0.014")
LIST(APPEND O_ATOM_SCALARS "samples" "64000 0.0")
LIST(APPEND O_ATOM_SCALARS "flux" "0.0 0.25")
@ -261,4 +261,23 @@ else()
SET_TESTS_PROPERTIES( ${QE_TEST_NAME}-atomO_pp-vmc_sdj-1-16
PROPERTIES DEPENDS ${LAST_TEST_NAME} )
# pw2casino Kinetic 11.64963538 affected by smearing.
# measued Kinetic 11.64798 +/- 0.00017
LIST(APPEND O_ATOM_NOJ_SCALARS "kinetic" "11.64798 0.04")
#LIST(APPEND O_ATOM_NOJ_SCALARS "flux" "0.0 0.3")
SOFTLINK_H5( ${QE_TEST_NAME} atomO_pp-vmc_noj-1-16 atomO atomO.pwscf.h5 LAST_TEST_NAME )
QMC_RUN_AND_CHECK(${QE_TEST_NAME}-atomO_pp-vmc_noj
"${CMAKE_SOURCE_DIR}/tests/solids/atomO_pp"
qmc_short_noj
qmc_short_noj.in.xml
1 16
O_ATOM_NOJ_SCALARS
0 # VMC
TRUE)
SET_TESTS_PROPERTIES( ${QE_TEST_NAME}-atomO_pp-vmc_noj-1-16
PROPERTIES DEPENDS ${LAST_TEST_NAME} )
endif()

View File

@ -12,7 +12,7 @@
ecutwfc = 200,
occupations = 'smearing'
smearing = 'gauss'
degauss = 0.02
degauss = 0.001
nspin = 2,
starting_magnetization(1) = 1.0
/

View File

@ -0,0 +1,68 @@
<?xml version="1.0"?>
<simulation>
<project id="qmc_short_noj" series="0">
<application name="qmcapp" role="molecu" class="serial" version="1.0"/>
</project>
<qmcsystem>
<simulationcell>
<parameter name="lattice" units="bohr">
20.0 0.0 0.0
0.0 20.0 0.0
0.0 0.0 20.0
</parameter>
<parameter name="bconds">
n n n
</parameter>
<parameter name="LR_dim_cutoff" > 15 </parameter>
</simulationcell>
<particleset name="ion0">
<group name="O" size="1" mass="29164.3928678">
<parameter name="charge" > 6 </parameter>
<parameter name="valence" > 6 </parameter>
<parameter name="atomicnumber" > 8 </parameter>
<parameter name="mass" > 29164.3928678 </parameter>
<attrib name="position" datatype="posArray" condition="1">
0.50000000 0.50000000 0.50000000
</attrib>
</group>
</particleset>
<particleset name="e" random="yes" randomsrc="ion0">>
<group name="u" size="4" mass="1.0">
<parameter name="charge" > -1 </parameter>
<parameter name="mass" > 1.0 </parameter>
</group>
<group name="d" size="2" mass="1.0">
<parameter name="charge" > -1 </parameter>
<parameter name="mass" > 1.0 </parameter>
</group>
</particleset>
<wavefunction name="psi0" target="e">
<determinantset type="einspline" href="atomO.pwscf.h5" tilematrix="1 0 0 0 1 0 0 0 1" twistnum="0" source="ion0" meshfactor="1.0" precision="double">
<slaterdeterminant>
<determinant id="updet" size="4">
<occupation mode="ground" spindataset="0"/>
</determinant>
<determinant id="downdet" size="2">
<occupation mode="ground" spindataset="1"/>
</determinant>
</slaterdeterminant>
</determinantset>
</wavefunction>
<hamiltonian name="h0" type="generic" target="e">
<pairpot type="coulomb" name="ElecElec" source="e" target="e"/>
<pairpot type="coulomb" name="IonIon" source="ion0" target="ion0"/>
<pairpot type="pseudo" name="PseudoPot" source="ion0" wavefunction="psi0" format="xml">
<pseudo elementType="O" href="O.BFD.xml"/>
</pairpot>
<estimator type="flux" name="Flux"/>
</hamiltonian>
</qmcsystem>
<qmc method="vmc" move="pbyp">
<parameter name="walkers" > 16 </parameter>
<parameter name="blocks" > 100 </parameter>
<parameter name="steps" > 40 </parameter>
<parameter name="subSteps" > 4 </parameter>
<parameter name="timestep" > 0.3 </parameter>
<parameter name="warmupSteps" > 100 </parameter>
</qmc>
</simulation>