05Diamond with k twist average

00 Introduction

From this tutorial, you can learn how to calculate Diamond (with k twist average) with JDFT ansatz. You can download all the input and output files from here.

01 DFT

cd 01trial_wavefunction/00makefort10/
turbogenius makefort10 -g -str 2101499.cif -detbasis cc-pVTZ -jasbasis cc-pVDZ -detcutbasis -jascutbasis -pp ccECP -complex
turbogenius makefort10 -r -post

cp fort.10 fort.10_in
turbogenius convertfort10mol -g -r -post

cd ../01DFT
cp ../00makefort10/fort.10 ./
cp ../00makefort10/pseudo.dat .

turbogenius prep -g -grid 0.10 0.10 0.10 -kpts 2 2 2 0 0 0

# on a local machine (parallel version)
mpirun -np 8x prep-mpi.x < prep.input > out_prep
# on a cluster machine (PBS)
qsub submit.sh
# on a cluster machine (Slurm)
sbatch submit.sh

turbogenius prep -post

02 Jastrow optimization

cd ../../02optimization/
cp ../01trial_wavefunction/01DFT/fort.10 fort.10
cp -r ../01trial_wavefunction/01DFT/turborvb.scratch turborvb.scratch
cp ../01trial_wavefunction/01DFT/pseudo.dat ./
turbogenius vmcopt -g -opt_onebody -opt_twobody -opt_jas_mat -optimizer lr -vmcoptsteps 100 -steps 200 -twist -kpts 2 2 2 0 0 0

# on a local machine (parallel version)
mpirun -np 8x turborvb-mpi.x < datasmin.input > out_min
# on a cluster machine (PBS)
qsub submit.sh
# on a cluster machine (Slurm)
sbatch submit.sh

turbogenius vmcopt -post -optwarmup 50 -plot

03 JDFT ansatz - VMC

cd ../03vmc/
cp ../02optimization/fort.10 fort.10
cp -r ../02optimization/turborvb.scratch turborvb.scratch
cp ../02optimization/pseudo.dat .
turbogenius vmc -g -twist -kpts 2 2 2 0 0 0

# on a local machine (parallel version)
mpirun -np 8x turborvb-mpi.x < datasvmc.input > out_vmc
# on a cluster machine (PBS)
qsub submit.sh
# on a cluster machine (Slurm)
sbatch submit.sh

turbogenius vmc -post -bin 10 -warmup 5