mirror of https://github.com/QMCPACK/qmcpack.git
434 lines
18 KiB
Plaintext
434 lines
18 KiB
Plaintext
Input file(s): vmc_b3lyp_noj.in.xml
|
|
|
|
=====================================================
|
|
QMCPACK 3.4.0
|
|
|
|
(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: ea23ca60d73178e72bd4639677282e070a3be112
|
|
Last git commit date: Tue Jun 26 14:08:17 2018 -0500
|
|
Last git commit subject: New Converged Pyscf references for FeCO6
|
|
=====================================================
|
|
Global options
|
|
|
|
MPI Nodes = 512
|
|
MPI Nodes per group = 512
|
|
MPI Group ID = 0
|
|
OMP_NUM_THREADS = 64
|
|
|
|
Precision used in this calculation, see definitions in the manual:
|
|
Base precision = double
|
|
Full precision = double
|
|
|
|
Input XML = vmc_b3lyp_noj.in.xml
|
|
|
|
Project = vmc_b3lyp_noj
|
|
date = 2018-07-03 21:34:22 UTC
|
|
host = Q02-I0-J01.vesta.itd
|
|
|
|
Random Number
|
|
-------------
|
|
Offset for the random number seeds based on time: 990
|
|
|
|
Range of prime numbers to use as seeds over processors and threads = 7853-398711
|
|
|
|
|
|
Particle Set
|
|
------------
|
|
Name: ion0
|
|
All the species have the same mass 1
|
|
Particle set size: 13
|
|
|
|
Particle Set
|
|
------------
|
|
Name: e
|
|
All the species have the same mass 1
|
|
Particle set size: 74
|
|
|
|
Wavefunction setup:
|
|
-------------------
|
|
Name: psi0
|
|
LCAO: SoaAtomicBasisSet<MultiQuintic,0>
|
|
|
|
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.
|
|
|
|
|
|
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.
|
|
|
|
Reading BasisSet from HDF5 file:../FeCO6.orbs.h5
|
|
<input node="atomicBasisSet" name="Gaussian" Morder="gaussian" angular="cartesian" elementType="Fe" normalized="no" basisType="Numerical" addSign="0" />
|
|
AO BasisSet for Fe
|
|
Angular momentum expanded in cartesian functions x^lx y^ly z^lz according to Gamess
|
|
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 0 0 0
|
|
R(n,l,m,s) 4 0 0 0
|
|
R(n,l,m,s) 5 1 0 0
|
|
R(n,l,m,s) 6 1 0 0
|
|
R(n,l,m,s) 7 1 0 0
|
|
R(n,l,m,s) 8 1 0 0
|
|
R(n,l,m,s) 9 1 0 0
|
|
R(n,l,m,s) 10 2 0 0
|
|
R(n,l,m,s) 11 2 0 0
|
|
R(n,l,m,s) 12 2 0 0
|
|
R(n,l,m,s) 13 2 0 0
|
|
R(n,l,m,s) 14 3 0 0
|
|
R(n,l,m,s) 15 3 0 0
|
|
R(n,l,m,s) 16 4 0 0
|
|
Expanding Ylm (angular function) according to Gamess using cartesian gaussians
|
|
Adding 1 cartesian gaussian orbitals for l= 0
|
|
Adding 1 cartesian gaussian orbitals for l= 0
|
|
Adding 1 cartesian gaussian orbitals for l= 0
|
|
Adding 1 cartesian gaussian orbitals for l= 0
|
|
Adding 1 cartesian gaussian orbitals for l= 0
|
|
Adding 3 cartesian gaussian orbitals for l= 1
|
|
Adding 3 cartesian gaussian orbitals for l= 1
|
|
Adding 3 cartesian gaussian orbitals for l= 1
|
|
Adding 3 cartesian gaussian orbitals for l= 1
|
|
Adding 3 cartesian gaussian orbitals for l= 1
|
|
Adding 6 cartesian gaussian orbitals for l= 2
|
|
Adding 6 cartesian gaussian orbitals for l= 2
|
|
Adding 6 cartesian gaussian orbitals for l= 2
|
|
Adding 6 cartesian gaussian orbitals for l= 2
|
|
Adding 10 cartesian gaussian orbitals for l= 3
|
|
Adding 10 cartesian gaussian orbitals for l= 3
|
|
Adding 15 cartesian gaussian orbitals for l= 4
|
|
Maximum Angular Momentum = 4
|
|
Number of Radial functors = 17
|
|
Basis size = 79
|
|
|
|
<input node="atomicBasisSet" name="Gaussian" Morder="gaussian" angular="cartesian" elementType="C" normalized="no" basisType="Numerical" addSign="0" />
|
|
AO BasisSet for C
|
|
Angular momentum expanded in cartesian functions x^lx y^ly z^lz according to Gamess
|
|
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 1 0 0
|
|
R(n,l,m,s) 6 2 0 0
|
|
R(n,l,m,s) 7 2 0 0
|
|
R(n,l,m,s) 8 3 0 0
|
|
Expanding Ylm (angular function) according to Gamess using cartesian gaussians
|
|
Adding 1 cartesian gaussian orbitals for l= 0
|
|
Adding 1 cartesian gaussian orbitals for l= 0
|
|
Adding 1 cartesian gaussian orbitals for l= 0
|
|
Adding 3 cartesian gaussian orbitals for l= 1
|
|
Adding 3 cartesian gaussian orbitals for l= 1
|
|
Adding 3 cartesian gaussian orbitals for l= 1
|
|
Adding 6 cartesian gaussian orbitals for l= 2
|
|
Adding 6 cartesian gaussian orbitals for l= 2
|
|
Adding 10 cartesian gaussian orbitals for l= 3
|
|
Maximum Angular Momentum = 3
|
|
Number of Radial functors = 9
|
|
Basis size = 34
|
|
|
|
<input node="atomicBasisSet" name="Gaussian" Morder="gaussian" angular="cartesian" elementType="O" normalized="no" basisType="Numerical" addSign="0" />
|
|
AO BasisSet for O
|
|
Angular momentum expanded in cartesian functions x^lx y^ly z^lz according to Gamess
|
|
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 1 0 0
|
|
R(n,l,m,s) 6 2 0 0
|
|
R(n,l,m,s) 7 2 0 0
|
|
R(n,l,m,s) 8 3 0 0
|
|
Expanding Ylm (angular function) according to Gamess using cartesian gaussians
|
|
Adding 1 cartesian gaussian orbitals for l= 0
|
|
Adding 1 cartesian gaussian orbitals for l= 0
|
|
Adding 1 cartesian gaussian orbitals for l= 0
|
|
Adding 3 cartesian gaussian orbitals for l= 1
|
|
Adding 3 cartesian gaussian orbitals for l= 1
|
|
Adding 3 cartesian gaussian orbitals for l= 1
|
|
Adding 6 cartesian gaussian orbitals for l= 2
|
|
Adding 6 cartesian gaussian orbitals for l= 2
|
|
Adding 10 cartesian gaussian orbitals for l= 3
|
|
Maximum Angular Momentum = 3
|
|
Number of Radial functors = 9
|
|
Basis size = 34
|
|
|
|
Created SPOSet builder named 'LCAOBSet' of type molecularorbital
|
|
Reuse SPOSetBuilder "LCAOBSet" type MolecularOrbital
|
|
Building SPOSet '' with '' basis set.
|
|
Reuse SPOSetBuilder "LCAOBSet" type MolecularOrbital
|
|
Building SPOSet '' with '' basis set.
|
|
Creating a determinant updet group=0 sposet=updet
|
|
Reusing a SPO set updet
|
|
Using DiracDeterminantBase
|
|
|
|
Creating a determinant downdet group=1 sposet=downdet
|
|
Reusing a SPO set downdet
|
|
Using DiracDeterminantBase
|
|
|
|
FermionWF = SlaterDet
|
|
QMCHamiltonian::addOperator Kinetic to H, physical Hamiltonian
|
|
QMCHamiltonian::addOperator ElecElec to H, physical Hamiltonian
|
|
QMCHamiltonian::addOperatorType added type coulomb named ElecElec
|
|
|
|
Distance table for similar particles (A-A):
|
|
source/target: ion0
|
|
Using structure-of-arrays (SoA) data layout
|
|
Distance computations use open boundary conditions in 3D.
|
|
|
|
QMCHamiltonian::addOperator IonIon to H, physical Hamiltonian
|
|
QMCHamiltonian::addOperatorType added type coulomb named IonIon
|
|
|
|
ECPotential builder for pseudopotential
|
|
|
|
Adding pseudopotential for C
|
|
Linear grid ri=0 rf=10 npts = 10001
|
|
ECPComponentBuilder::buildSemiLocalAndLocal
|
|
Assuming Hartree unit
|
|
Number of angular momentum channels 2
|
|
Maximum angular momentum channel 1
|
|
Creating a Linear Grid Rmax=1.7
|
|
Using global grid with delta = 0.001
|
|
Making L=1 a local potential with a radial cutoff of 9.999
|
|
Non-local pseudopotential parameters
|
|
Maximum angular mementum = 0
|
|
Number of non-local channels = 1
|
|
l(0)=0
|
|
Cutoff radius = 1.7
|
|
Spherical grids and weights:
|
|
1 0 0 0.08333333333
|
|
-1 1.224646799e-16 0 0.08333333333
|
|
0.4472135955 0.894427191 0 0.08333333333
|
|
-0.4472135955 0.7236067977 0.5257311121 0.08333333333
|
|
0.4472135955 0.2763932023 0.8506508084 0.08333333333
|
|
-0.4472135955 -0.2763932023 0.8506508084 0.08333333333
|
|
0.4472135955 -0.7236067977 0.5257311121 0.08333333333
|
|
-0.4472135955 -0.894427191 1.095357397e-16 0.08333333333
|
|
0.4472135955 -0.7236067977 -0.5257311121 0.08333333333
|
|
-0.4472135955 -0.2763932023 -0.8506508084 0.08333333333
|
|
0.4472135955 0.2763932023 -0.8506508084 0.08333333333
|
|
-0.4472135955 0.7236067977 -0.5257311121 0.08333333333
|
|
Maximum cutoff radius 1.7
|
|
|
|
Adding pseudopotential for Fe
|
|
Linear grid ri=0 rf=10 npts = 10001
|
|
ECPComponentBuilder::buildSemiLocalAndLocal
|
|
Assuming Hartree unit
|
|
Number of angular momentum channels 3
|
|
Maximum angular momentum channel 2
|
|
Creating a Linear Grid Rmax=1.355653
|
|
Using global grid with delta = 0.001
|
|
Making L=2 a local potential with a radial cutoff of 9.999
|
|
Non-local pseudopotential parameters
|
|
Maximum angular mementum = 1
|
|
Number of non-local channels = 2
|
|
l(0)=0
|
|
l(1)=1
|
|
Cutoff radius = 1.356
|
|
Spherical grids and weights:
|
|
1 0 0 0.08333333333
|
|
-1 1.224646799e-16 0 0.08333333333
|
|
0.4472135955 0.894427191 0 0.08333333333
|
|
-0.4472135955 0.7236067977 0.5257311121 0.08333333333
|
|
0.4472135955 0.2763932023 0.8506508084 0.08333333333
|
|
-0.4472135955 -0.2763932023 0.8506508084 0.08333333333
|
|
0.4472135955 -0.7236067977 0.5257311121 0.08333333333
|
|
-0.4472135955 -0.894427191 1.095357397e-16 0.08333333333
|
|
0.4472135955 -0.7236067977 -0.5257311121 0.08333333333
|
|
-0.4472135955 -0.2763932023 -0.8506508084 0.08333333333
|
|
0.4472135955 0.2763932023 -0.8506508084 0.08333333333
|
|
-0.4472135955 0.7236067977 -0.5257311121 0.08333333333
|
|
Maximum cutoff radius 1.356
|
|
|
|
Adding pseudopotential for O
|
|
Linear grid ri=0 rf=10 npts = 10001
|
|
ECPComponentBuilder::buildSemiLocalAndLocal
|
|
Assuming Hartree unit
|
|
Number of angular momentum channels 2
|
|
Maximum angular momentum channel 1
|
|
Creating a Linear Grid Rmax=1.31
|
|
Using global grid with delta = 0.001
|
|
Making L=1 a local potential with a radial cutoff of 9.999
|
|
Non-local pseudopotential parameters
|
|
Maximum angular mementum = 0
|
|
Number of non-local channels = 1
|
|
l(0)=0
|
|
Cutoff radius = 1.31
|
|
Spherical grids and weights:
|
|
1 0 0 0.08333333333
|
|
-1 1.224646799e-16 0 0.08333333333
|
|
0.4472135955 0.894427191 0 0.08333333333
|
|
-0.4472135955 0.7236067977 0.5257311121 0.08333333333
|
|
0.4472135955 0.2763932023 0.8506508084 0.08333333333
|
|
-0.4472135955 -0.2763932023 0.8506508084 0.08333333333
|
|
0.4472135955 -0.7236067977 0.5257311121 0.08333333333
|
|
-0.4472135955 -0.894427191 1.095357397e-16 0.08333333333
|
|
0.4472135955 -0.7236067977 -0.5257311121 0.08333333333
|
|
-0.4472135955 -0.2763932023 -0.8506508084 0.08333333333
|
|
0.4472135955 0.2763932023 -0.8506508084 0.08333333333
|
|
-0.4472135955 0.7236067977 -0.5257311121 0.08333333333
|
|
Maximum cutoff radius 1.31
|
|
QMCHamiltonian::addOperator LocalECP to H, physical Hamiltonian
|
|
|
|
Using NonLocalECP potential
|
|
Maximum grid on a sphere for NonLocalECPotential: 12
|
|
QMCHamiltonian::addOperator NonLocalECP to H, physical Hamiltonian
|
|
ECPotential::Rmax 9.999
|
|
QMCHamiltonian::addOperatorType added type pseudo named PseudoPot
|
|
|
|
QMCHamiltonian::add2WalkerProperty added
|
|
5 to P::PropertyList
|
|
0 to P::Collectables
|
|
starting Index of the observables in P::PropertyList = 9
|
|
ParticleSetPool::randomize
|
|
<init source="ion0" target="e">
|
|
</init>
|
|
Initialization Execution time = 3.887 secs
|
|
=========================================================
|
|
Summary of QMC systems
|
|
=========================================================
|
|
ParticleSetPool has:
|
|
|
|
ParticleSet e : 0 39 74
|
|
|
|
74
|
|
|
|
u 1.5258666473e+00 5.4722977408e-02 4.7827153997e-01
|
|
u -1.6235204341e-01 6.9517927260e-01 1.4319099103e+00
|
|
u -1.5850049684e+00 5.9776436312e-02 -2.1018633901e-01
|
|
u -7.7051097795e-01 -1.0741388707e+00 -9.0140749960e-01
|
|
u 8.9985419468e-01 6.9838981564e-01 -1.1236185632e+00
|
|
u 1.4108500116e+00 6.7342542969e-01 -3.4059835793e-01
|
|
u -7.8574043398e-02 1.5200140279e+00 4.9333794449e-01
|
|
u 9.3886575120e-01 -1.0568139888e+00 7.4945314950e-01
|
|
u -5.9707221260e-01 3.3931912820e-01 4.9016922561e+00
|
|
u -7.1715203329e-01 1.5738062694e-01 4.0149063627e+00
|
|
(... and 64 more particle positions ...)
|
|
|
|
ParticleSet ion0 : 0 1 7 13
|
|
|
|
13
|
|
|
|
Fe 2.4566000000e-06 -1.8900000000e-07 -1.8900000000e-07
|
|
C 6.9920000000e-06 -1.3642687800e-02 4.4230707260e+00
|
|
O -5.6502800000e-05 -7.8540602500e-02 6.5288519205e+00
|
|
C 4.3520243988e+00 -2.5889200000e-05 -2.5511300000e-05
|
|
O 6.4595484876e+00 -3.2125300000e-05 -3.1747400000e-05
|
|
C -4.3520157061e+00 2.5133400000e-05 2.4755400000e-05
|
|
O -6.4595397949e+00 3.1558400000e-05 3.1180500000e-05
|
|
C -3.5905000000e-06 1.3643065700e-02 -4.4230705371e+00
|
|
O 4.7243100000e-05 7.8541169400e-02 -6.5288515426e+00
|
|
C -5.2912000000e-06 -4.4162633661e+00 1.2432696200e-02
|
|
(... and 3 more particle positions ...)
|
|
|
|
Hamiltonian h0
|
|
Kinetic Kinetic energy
|
|
ElecElec CoulombAA source/target e
|
|
|
|
IonIon CoulombAA source/target ion0
|
|
|
|
LocalECP LocalECPotential: ion0
|
|
NonLocalECP NonLocalECPotential: ion0
|
|
|
|
=========================================================
|
|
Start VMCSingleOMP
|
|
File Root vmc_b3lyp_noj.s000 append = no
|
|
=========================================================
|
|
Adding 64 walkers to 0 existing sets
|
|
Total number of walkers: 3.2768000000e+04
|
|
Total weight: 3.2768000000e+04
|
|
Resetting Properties of the walkers 1 x 14
|
|
|
|
<vmc function="put">
|
|
qmc_counter=0 my_counter=0
|
|
time step = 1.0000000000e-01
|
|
blocks = 200
|
|
steps = 100
|
|
substeps = 3
|
|
current = 0
|
|
target samples = 1.6000000000e+05
|
|
walkers/mpi = 64
|
|
|
|
stepsbetweensamples = 4000
|
|
<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">5.0000000000e+00</parameter>
|
|
<parameter name="maxcpusecs" condition="real">3.6000000000e+05</parameter>
|
|
<parameter name="record_configs" condition="int">0</parameter>
|
|
<parameter name="record_walkers" condition="int">4000</parameter>
|
|
<parameter name="recordconfigs" condition="int">0</parameter>
|
|
<parameter name="recordwalkers" condition="int">4000</parameter>
|
|
<parameter name="rewind" condition="int">0</parameter>
|
|
<parameter name="samples" condition="real">1.6000000000e+05</parameter>
|
|
<parameter name="samplesperthread" condition="real">5.0000000000e+00</parameter>
|
|
<parameter name="steps" condition="int">100</parameter>
|
|
<parameter name="stepsbetweensamples" condition="int">4000</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">64</parameter>
|
|
<parameter name="warmup_steps" condition="int">50</parameter>
|
|
<parameter name="warmupsteps" condition="int">50</parameter>
|
|
DumpConfig==false Nothing (configurations, state) will be saved.
|
|
Walker Samples are dumped every 4000 steps.
|
|
</vmc>
|
|
Adding a default LocalEnergyEstimator for the MainEstimator
|
|
CloneManager::makeClones makes 64 clones for W/Psi/H.
|
|
Cloning methods for both Psi and H are used
|
|
===== Memory Usage before cloning =====
|
|
Allocated heap: 27.35 MB, avail. heap: 16228.63 MB
|
|
Allocated stack: 0.02 MB, avail. stack: 16228.63 MB
|
|
==================================================
|
|
===== Memory Usage after cloning =====
|
|
Allocated heap: 1040.00 MB, avail. heap: 15215.98 MB
|
|
Allocated stack: 0.02 MB, avail. stack: 15215.98 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
|
|
|
|
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
|
|
Using Locality Approximation
|
|
===== Memory Usage after the buffer registration =====
|
|
Allocated heap: 1040.00 MB, avail. heap: 15215.98 MB
|
|
Allocated stack: 0.02 MB, avail. stack: 15215.98 MB
|
|
==================================================
|
|
Anonymous Buffer size per walker : 114368 Bytes.
|
|
MEMORY increase 6 MB VMCSingleOMP::resetRun
|
|
====================================================
|
|
SimpleFixedNodeBranch::finalize after a VMC block
|
|
QMC counter = 0
|
|
time step = 0.1
|
|
reference energy = -250.477
|
|
reference variance = 38.3319
|
|
====================================================
|
|
QMC Execution time = 5.3101e+03 secs
|
|
Total Execution time = 5.3101e+03 secs
|
|
|
|
=========================================================
|
|
A new xml input file : vmc_b3lyp_noj.s000.cont.xml
|