mirror of https://github.com/QMCPACK/qmcpack.git
437 lines
17 KiB
Plaintext
437 lines
17 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-dirty
|
|
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 = 1
|
|
MPI Nodes per group = 1
|
|
MPI Group ID = 0
|
|
OMP_NUM_THREADS = 16
|
|
|
|
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-05 08:47:03 CDT
|
|
host = beboplogin2
|
|
|
|
Random Number
|
|
-------------
|
|
Offset for the random number seeds based on time: 343
|
|
|
|
Range of prime numbers to use as seeds over processors and threads = 2333-2437
|
|
|
|
|
|
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,1>
|
|
|
|
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.h5
|
|
<input node="atomicBasisSet" name="bfd-vtz" Morder="pyscf" angular="spherical" elementType="Fe" normalized="no" basisType="Numerical" addSign="0" />
|
|
AO BasisSet for Fe
|
|
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 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 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 1 spherical orbitals
|
|
Adding 1 spherical orbitals
|
|
Adding 3 spherical orbitals
|
|
Adding 3 spherical orbitals
|
|
Adding 3 spherical orbitals
|
|
Adding 3 spherical orbitals
|
|
Adding 3 spherical orbitals
|
|
Adding 5 spherical orbitals
|
|
Adding 5 spherical orbitals
|
|
Adding 5 spherical orbitals
|
|
Adding 5 spherical orbitals
|
|
Adding 7 spherical orbitals
|
|
Adding 7 spherical orbitals
|
|
Adding 9 spherical orbitals
|
|
Maximum Angular Momentum = 4
|
|
Number of Radial functors = 17
|
|
Basis size = 63
|
|
|
|
<input node="atomicBasisSet" name="bfd-vtz" Morder="pyscf" angular="spherical" elementType="C" normalized="no" basisType="Numerical" addSign="0" />
|
|
AO BasisSet for C
|
|
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 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 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 3 spherical orbitals
|
|
Adding 5 spherical orbitals
|
|
Adding 5 spherical orbitals
|
|
Adding 7 spherical orbitals
|
|
Maximum Angular Momentum = 3
|
|
Number of Radial functors = 9
|
|
Basis size = 29
|
|
|
|
<input node="atomicBasisSet" name="bfd-vtz" Morder="pyscf" angular="spherical" elementType="O" normalized="no" basisType="Numerical" addSign="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 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 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 3 spherical orbitals
|
|
Adding 5 spherical orbitals
|
|
Adding 5 spherical orbitals
|
|
Adding 7 spherical orbitals
|
|
Maximum Angular Momentum = 3
|
|
Number of Radial functors = 9
|
|
Basis size = 29
|
|
|
|
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 = 0.1882 secs
|
|
=========================================================
|
|
Summary of QMC systems
|
|
=========================================================
|
|
ParticleSetPool has:
|
|
|
|
ParticleSet e : 0 39 74
|
|
|
|
74
|
|
|
|
u -1.1939420404e-01 2.0937909212e-01 -1.5817412405e+00
|
|
u -1.4419240406e+00 2.3217184298e-01 6.5340942535e-01
|
|
u -9.6175855203e-01 -1.2766240414e+00 -7.2438450046e-02
|
|
u 8.0908192819e-01 8.8438116371e-01 1.0598394365e+00
|
|
u -1.1801661212e+00 1.0777500871e+00 7.5209320504e-02
|
|
u 1.3155823127e+00 -7.0622586347e-03 9.1060387949e-01
|
|
u -4.0209130721e-01 1.4084367563e+00 -6.4391476117e-01
|
|
u -1.1863485716e+00 9.1731615179e-01 5.5776496650e-01
|
|
u 2.3372381527e-01 -5.1136850995e-01 5.0616708784e+00
|
|
u -1.6568707504e-01 -1.3454882728e-01 3.6057024184e+00
|
|
(... and 64 more particle positions ...)
|
|
|
|
ParticleSet ion0 : 0 1 7 13
|
|
|
|
13
|
|
|
|
Fe 2.4566439619e-06 -1.8897261246e-07 -1.8897261246e-07
|
|
C 6.9919866609e-06 -1.3642688784e-02 4.4230710463e+00
|
|
O -5.6502811124e-05 -7.8540608217e-02 6.5288523933e+00
|
|
C 4.3520247140e+00 -2.5889247907e-05 -2.5511302682e-05
|
|
O 6.4595489554e+00 -3.2125344118e-05 -3.1747398893e-05
|
|
C -4.3520160212e+00 2.5133357457e-05 2.4755412232e-05
|
|
O -6.4595402626e+00 3.1558426280e-05 3.1180481055e-05
|
|
C -3.5904796367e-06 1.3643066729e-02 -4.4230708574e+00
|
|
O 4.7243153114e-05 7.8541175134e-02 -6.5288520154e+00
|
|
C -5.2912331488e-06 -4.4162636859e+00 1.2432697146e-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 16 walkers to 0 existing sets
|
|
Total number of walkers: 1.6000000000e+01
|
|
Total weight: 1.6000000000e+01
|
|
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 = 15
|
|
substeps = 3
|
|
current = 0
|
|
target samples = 1.6000000000e+04
|
|
walkers/mpi = 16
|
|
|
|
stepsbetweensamples = 3
|
|
<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+03</parameter>
|
|
<parameter name="maxcpusecs" condition="real">3.6000000000e+05</parameter>
|
|
<parameter name="record_configs" condition="int">0</parameter>
|
|
<parameter name="record_walkers" condition="int">3</parameter>
|
|
<parameter name="recordconfigs" condition="int">0</parameter>
|
|
<parameter name="recordwalkers" condition="int">3</parameter>
|
|
<parameter name="rewind" condition="int">0</parameter>
|
|
<parameter name="samples" condition="real">1.6000000000e+04</parameter>
|
|
<parameter name="samplesperthread" condition="real">1.0000000000e+03</parameter>
|
|
<parameter name="steps" condition="int">15</parameter>
|
|
<parameter name="stepsbetweensamples" condition="int">3</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">16</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 3 steps.
|
|
</vmc>
|
|
Adding a default LocalEnergyEstimator for the MainEstimator
|
|
CloneManager::makeClones makes 16 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 : 119363 MB
|
|
Memory footprint by rank 0 on node 0 : 68 MB
|
|
==================================================
|
|
===== Memory Usage after cloning =====
|
|
Available memory on node 0, free + buffers : 119345 MB
|
|
Memory footprint by rank 0 on node 0 : 87 MB
|
|
==================================================
|
|
Initial partition of walkers 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
|
|
|
|
Using Particle by Particle moves
|
|
Walker moves with drift
|
|
Total Sample Size =16000
|
|
Walker distribution on root = 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
|
|
Using Locality Approximation
|
|
===== Memory Usage after the buffer registration =====
|
|
Available memory on node 0, free + buffers : 119277 MB
|
|
Memory footprint by rank 0 on node 0 : 157 MB
|
|
==================================================
|
|
Anonymous Buffer size per walker : 114368 Bytes.
|
|
MEMORY increase 1 MB VMCSingleOMP::resetRun
|
|
====================================================
|
|
SimpleFixedNodeBranch::finalize after a VMC block
|
|
QMC counter = 0
|
|
time step = 0.1
|
|
reference energy = -250.208
|
|
reference variance = 16.7191
|
|
====================================================
|
|
QMC Execution time = 3.3151e+01 secs
|
|
Total Execution time = 3.3154e+01 secs
|
|
|
|
=========================================================
|
|
A new xml input file : vmc_b3lyp_noj.s000.cont.xml
|