qmcpack/tests/molecules/O_ae_pyscf_UHF/qmc-ref/qmc-ref.out

260 lines
10 KiB
Plaintext

Input file(s): qmc-ref.in.xml
=====================================================
QMCPACK 3.9.2
(c) Copyright 2003- QMCPACK developers
Please cite:
J. Kim et al. J. Phys. Cond. Mat. 30 195901 (2018)
https://doi.org/10.1088/1361-648X/aab9c3
Git branch: HEAD
Last git commit: 8064ebfec9a8dd88748e8864013168c9aa7e05dd-dirty
Last git commit date: Tue Apr 28 13:44:48 2020 -0500
Last git commit subject: Merge pull request #2412 from QMCPACK/rc_392
=====================================================
Global options
Total number of MPI ranks = 128
Number of MPI groups = 1
MPI group ID = 0
Number of ranks in group = 128
MPI ranks per node = 1
OMP 1st level threads = 128
OMP nested threading disabled or only 1 thread on the 2nd level
Precision used in this calculation, see definitions in the manual:
Base precision = double
Full precision = double
Structure-of-arrays (SoA) optimization enabled
Input XML = qmc-ref.in.xml
Project = qmc-ref
date = 2020-08-25 21:24:45 UTC
host = nid02749
Random Number
-------------
Offset for the random number seeds based on time: 413
Range of prime numbers to use as seeds over processors and threads = 2857-185621
Particle Set
------------
Name: ion0
All the species have the same mass 1
Particle set size: 1
Particle Set
------------
Name: e
All the species have the same mass 1
Particle set size: 8
Wavefunction setup:
-------------------
Name: psi0
LCAO: SoaAtomicBasisSet<MultiQuintic,1>
Reading BasisSet from HDF5 file:../O-UHF-Triplet.h5
<input node="atomicBasisSet" name="cc-pcvdz" expandYlm="pyscf" angular="spherical" elementType="O" normalized="no" type="Numerical" expM="0" />
AO BasisSet for O
Spherical Harmonics contain (-1)^m factor
Angular momentum m expanded as -l, ... ,l, with the exception of L=1 (1,-1,0)
Grid is created by the input paremters in h5
Using log grid ri = 1e-06 rf = 100 npts = 1001
R(n,l,m,s) 0 0 0 0
R(n,l,m,s) 1 0 0 0
R(n,l,m,s) 2 0 0 0
R(n,l,m,s) 3 1 0 0
R(n,l,m,s) 4 1 0 0
R(n,l,m,s) 5 2 0 0
R(n,l,m,s) 6 0 0 0
R(n,l,m,s) 7 1 0 0
Expanding Ylm as L=1 as (1,-1,0) and L>1 as -l,-l+1,...,l-1,l
Adding 1 spherical orbitals
Adding 1 spherical orbitals
Adding 1 spherical orbitals
Adding 3 spherical orbitals
Adding 3 spherical orbitals
Adding 5 spherical orbitals
Adding 1 spherical orbitals
Adding 3 spherical orbitals
Setting cutoff radius 5.36056e+06
Maximum Angular Momentum = 2
Number of Radial functors = 8
Basis size = 18
Created SPOSet builder named 'LCAOBSet' of type molecularorbital
Reuse SPOSetBuilder "LCAOBSet" type MolecularOrbital
Building SPOSet '' with '' basis set.
Reading cusp info from : ../updet.cuspInfo.xml
Found precomputed cusp data for spo set: updet
Reuse SPOSetBuilder "LCAOBSet" type MolecularOrbital
Building SPOSet '' with '' basis set.
Reading cusp info from : ../downdet.cuspInfo.xml
Found precomputed cusp data for spo set: downdet
Creating a determinant updet group=0 sposet=updet
Reusing a SPO set updet
Using DiracDeterminant with DelayedUpdate engine
Setting delay_rank by default!
Using rank-1 Sherman-Morrison Fahy update
Creating a determinant downdet group=1 sposet=downdet
Reusing a SPO set downdet
Using DiracDeterminant with DelayedUpdate engine
Setting delay_rank by default!
Using rank-1 Sherman-Morrison Fahy update
Added a fermionic WaveFunctionComponent SlaterDet
QMCHamiltonian::addOperator Kinetic to H, physical Hamiltonian
QMCHamiltonian::addOperator ElecElec to H, physical Hamiltonian
QMCHamiltonian::addOperatorType added type coulomb named ElecElec
CoulombAA for ion0 is not created. Number of particles == 1 and nonPeriodic
QMCHamiltonian::addOperator IonElec to H, physical Hamiltonian
QMCHamiltonian::addOperatorType added type coulomb named IonElec
QMCHamiltonian::add2WalkerProperty added
3 to P::PropertyList
0 to P::Collectables
starting Index of the observables in P::PropertyList = 9
ParticleSetPool::randomize 1 ParticleSet.
<init source="ion0" target="e">
</init>
Initialization Execution time = 0.4639 secs
=========================================================
Summary of QMC systems
=========================================================
ParticleSetPool has:
ParticleSet 'e' contains 8 particles : u(5) d(3)
u 3.5805105668e-01 1.2738498874e+00 2.9598736260e+00
u -3.4516470307e-01 -2.6255194705e-01 9.7135235510e-01
u -9.6195779554e-01 -3.2158824604e-01 5.3304454552e-01
u 1.0513600893e+00 1.8645583639e+00 8.2647665167e-01
u -4.0333505059e-01 1.5380386945e+00 -3.4886659200e-01
d 1.2669481503e+00 1.6409343376e+00 -8.3350320616e-01
d -2.6816571456e+00 2.0146360928e+00 -5.1108979306e-01
d -1.4895366316e+00 6.8386985525e-01 1.8565503426e+00
Distance table for dissimilar particles (A-B):
source: ion0 target: e
Using structure-of-arrays (SoA) data layout
Distance computations use open boundary conditions in 3D.
Distance table for similar particles (A-A):
source/target: e
Using structure-of-arrays (SoA) data layout
Distance computations use open boundary conditions in 3D.
ParticleSet 'ion0' contains 1 particles : O(1)
O 0.0000000000e+00 0.0000000000e+00 0.0000000000e+00
Hamiltonian h0
Kinetic Kinetic energy
ElecElec CoulombAB source=e
IonElec CoulombAA source/target ion0
=========================================================
Start VMC
File Root qmc-ref.s000 append = no
=========================================================
Resetting walkers
Adding 128 walkers to 0 existing sets
Total number of walkers: 1.6384000000e+04
Total weight: 1.6384000000e+04
Resetting Properties of the walkers 1 x 12
<vmc function="put">
qmc_counter=0 my_counter=0
time step = 1.0000000000e-01
blocks = 200
steps = 30000
substeps = 3
current = 0
target samples = 1.6000000000e+05
walkers/mpi = 128
stepsbetweensamples = -2931
SpinMoves = no
<parameter name="blocks" condition="int">200</parameter>
<parameter name="blocks_between_recompute" condition="int">0</parameter>
<parameter name="check_properties" condition="int">100</parameter>
<parameter name="checkproperties" condition="int">100</parameter>
<parameter name="current" condition="int">0</parameter>
<parameter name="dmcwalkersperthread" condition="real">1.0000000000e+01</parameter>
<parameter name="maxcpusecs" condition="real">3.6000000000e+05</parameter>
<parameter name="record_configs" condition="int">0</parameter>
<parameter name="record_walkers" condition="int">-2931</parameter>
<parameter name="recordconfigs" condition="int">0</parameter>
<parameter name="recordwalkers" condition="int">-2931</parameter>
<parameter name="rewind" condition="int">0</parameter>
<parameter name="samples" condition="real">1.6000000000e+05</parameter>
<parameter name="samplesperthread" condition="real">1.0000000000e+01</parameter>
<parameter name="spinmass" condition="double">1.0000000000e+00</parameter>
<parameter name="spinmoves" condition="string">no</parameter>
<parameter name="steps" condition="int">30000</parameter>
<parameter name="stepsbetweensamples" condition="int">-2931</parameter>
<parameter name="store_configs" condition="int">0</parameter>
<parameter name="storeconfigs" condition="int">0</parameter>
<parameter name="sub_steps" condition="int">3</parameter>
<parameter name="substeps" condition="int">3</parameter>
<parameter name="tau" condition="au">1.0000000000e-01</parameter>
<parameter name="time_step" condition="au">1.0000000000e-01</parameter>
<parameter name="timestep" condition="au">1.0000000000e-01</parameter>
<parameter name="use_drift" condition="string">yes</parameter>
<parameter name="usedrift" condition="string">yes</parameter>
<parameter name="walkers" condition="int">128</parameter>
<parameter name="warmup_steps" condition="int">50</parameter>
<parameter name="warmupsteps" condition="int">50</parameter>
DumpConfig==false Nothing (configurations, state) will be saved.
</vmc>
Set drift_modifier UNR parameter a = 1.0000000000e+00
Adding a default LocalEnergyEstimator for the MainEstimator
CloneManager::makeClones makes 128 clones for W/Psi/H.
Cloning methods for both Psi and H are used
===== Memory Usage before cloning =====
Available memory on node 0, free + buffers : 189302 MB
Memory footprint by rank 0 on node 0 : 58 MB
==================================================
===== Memory Usage after cloning =====
Available memory on node 0, free + buffers : 189302 MB
Memory footprint by rank 0 on node 0 : 58 MB
==================================================
Initial partition of walkers 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128
Using Particle by Particle moves
Walker moves with drift
Total Sample Size =163840
Walker distribution on root = 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128
===== Memory Usage after the buffer registration =====
Available memory on node 0, free + buffers : 189185 MB
Memory footprint by rank 0 on node 0 : 77 MB
==================================================
Anonymous Buffer size per walker : 4256 Bytes.
MEMORY increase 0 MB VMC::resetRun
====================================================
SimpleFixedNodeBranch::finalize after a VMC block
QMC counter = 0
time step = 0.1
reference energy = -74.7993
reference variance = 16.9662
====================================================
QMC Execution time = 4.8162e+03 secs
Total Execution time = 4.8163e+03 secs
=========================================================
A new xml input file : qmc-ref.s000.cont.xml
Application 21533412 resources: utime ~73502637s, stime ~79663s, Rss ~85912, inblocks ~1263984, outblocks ~19456