mirror of https://github.com/abinit/abinit.git
129 lines
4.3 KiB
Plaintext
129 lines
4.3 KiB
Plaintext
#Testing CTQMC local susceptibility
|
|
fftalg 112 #only for test
|
|
iomode 3 #only for test
|
|
|
|
#Multi_dataset parameters
|
|
ndtset 2 # Two datasets
|
|
jdtset 1 2 # one followed by two
|
|
getwfk -1 # Second one takes the density of first
|
|
|
|
pawprtvol 3 # Printing additional information
|
|
prtvol 4 #
|
|
prtden -1 #
|
|
|
|
############################
|
|
# Structural Parameters
|
|
############################
|
|
#Definition of unit cell
|
|
acell 3*5.44761 # lattice parameter
|
|
rprim -0.5 0.5 0.5 # Primitive cell
|
|
0.5 -0.5 0.5
|
|
0.5 0.5 -0.5
|
|
|
|
#Definition atom
|
|
ntypat 1 # number of different types of atoms
|
|
znucl 26 # z number of atom
|
|
|
|
#Definition of the atoms
|
|
natom 1 # number of atom
|
|
typat 1 # type of atom
|
|
xred 0.0 0.0 0.0 # reduced coordinates for primitive cell
|
|
|
|
#Definition pseudopotentials
|
|
pp_dirpath "$ABI_PSPDIR" # repository with pseudopotentials
|
|
pseudos "26fe.lda2.paw"
|
|
|
|
#Definition of the k-point grid
|
|
kptopt 4
|
|
ngkpt 4 4 4 # Reciprocal space vectors are built from
|
|
# the rprim parameters. This is the size of the
|
|
# reciprocal k-points.
|
|
nshiftk 2 # Convergence of the density with regular shifts.
|
|
shiftk 0.25 0.25 0.25
|
|
-0.25 -0.25 -0.25
|
|
|
|
############################
|
|
# Energy Parameters
|
|
############################
|
|
|
|
ecut 20 # Maximal plane-wave kinetic energy cut-off, in Hartree
|
|
pawecutdg 50 # PAW: Energy Cutoff for the Double Grid
|
|
nband 40 # Number of bands
|
|
occopt 3 # Smearing type
|
|
tsmear 3000 K # Electronic smearing temperature
|
|
|
|
|
|
############################
|
|
# DFT Calculation
|
|
############################
|
|
usedmft1 0 # DFT alone
|
|
nstep1 30 # Number of iterations for the DFT convergence
|
|
tolvrs 1.0d-10 # Tolerance of the density for DFT convergence
|
|
|
|
#Definition of the spin properties
|
|
nsppol 1 # No spin-polarization
|
|
nspden 4 # No Scalar magnetism
|
|
nspinor 2 # No Spin-orbit
|
|
spinat 0 0 0 # Magnetization vector
|
|
spnorbscl 10 # SOC scaling
|
|
|
|
############################
|
|
# DMFT Calculation
|
|
############################
|
|
usedmft2 1 # DMFT calculation
|
|
nstep2 1 # Number of iterations for the DFT+DMFT convergence
|
|
nline2 5 # Number of line minimisations
|
|
nnsclo2 5 # Number of non-self consistent loops
|
|
|
|
#Definition for DMFT loop
|
|
dmftbandi 9 # First KS band included in the projection for Wannier orbitals
|
|
dmftbandf 26 # Last KS band included in the projection for Wannier orbitals
|
|
dmft_nwlo 500 # Logarythmic frequency mesh
|
|
dmft_nwli 100000 # Linear frequency mesh
|
|
dmft_iter 1 # Number of iterations of the DMFT part.
|
|
# We often use single-shot, since anyway the charge density
|
|
# changes through the DFT+DMFT anyway.
|
|
dmftcheck 0
|
|
dmft_rslf 1 # Read self-energy, if nothing (like here) initialize.
|
|
dmft_mxsf 0.7 # Mixing of the old and new self-energy at every iterations.
|
|
dmft_dc 5 # Double counting type. 1 is Fully Localized Limit (FLL)
|
|
dmftqmc_seed 2
|
|
|
|
#Definition of impurity solver
|
|
dmft_solv = 8 # CT-QMC solver in density-density approximation
|
|
dmftqmc_l 200 # Number of time slices for G(tau).
|
|
dmftqmc_n 1.d6 # Number of QMC sweeps
|
|
dmftqmc_therm 10000 # Thermalization
|
|
dmftctqmc_gmove 0 # Global move occurence in QMC
|
|
dmftctqmc_order 50 # Perturbation order
|
|
dmftctqmc_config 2 # Activate Histogram and spin susceptibility
|
|
|
|
#Definition of DFT+U
|
|
usepawu1 4 #
|
|
usepawu2 14 #
|
|
dmatpuopt 2 #
|
|
lpawu 2 # Angular momentum for the projected Hamiltonian
|
|
upawu1 0.00 eV
|
|
upawu2 5.5 eV # Values of U for each angular momentum
|
|
jpawu1 0.00 eV
|
|
jpawu2 0.84 eV # Values of J
|
|
|
|
|
|
|
|
#%%<BEGIN TEST_INFO>
|
|
#%% [setup]
|
|
#%% executable = abinit
|
|
#%% exclude_builders = eos_gnu_13.2_mpich
|
|
#%% [paral_info]
|
|
#%% nprocs_to_test = 10
|
|
#%% max_nprocs = 10
|
|
#%% [NCPU_10]
|
|
#%% files_to_test =
|
|
#%% t101_MPI10o_DS2_LocalMagnSuscept_atom_01.dat, tolnlines =0, tolabs =0.001, tolrel = 0.001, fld_options=-medium
|
|
#%% [extra_info]
|
|
#%% authors = F. Gendron
|
|
#%% keywords = DMFT, CTQMC
|
|
#%% description = Test local susceptibility from CTQMC at spin-orbit level. Not realistic
|
|
#%% topics = DMFT
|
|
#%%<END TEST_INFO>
|