Questaal Home
Navigation

Susceptibility calculations for Nickel

Preparing the PARAMS file

Here we provide a standard PARAMS file for CT-QMC calculation for two-particle Green’s functions for Nickel at =50 (eV) which should be run on at least 1000 cores for a day to produce high quality data over all energies.

OffDiagonal  real
Sig  Sig.out
Gf  Gf.out
mode GH
SampleVertex 100
nOm 2
Delta  Delta.inp
cix  actqmc.cix
Nmax  3000         # Maximum perturbation order allowed
nom  200          # Number of Matsubara frequency points sampled
exe  ctqmc        # Name of the executable
tsample  50       # How often to record measurements
Ed [ -26.407760, -26.407519, -27.011761, -26.407397, -27.011682, -25.662041, -25.662114, -26.272431, -25.662012, -26.272482] # Ed: Eimp-Edc for PARAMS
M  80000000     # Total number of Monte Carlo steps per core
Ncout  2000000            # How often to print out info
PChangeOrder  0.9         # Ratio between trial steps: add-remove-a-kink / move-a-kink
mu   26.407760 # mu = -Eimp[0] for PARAMS
warmup  5000000           # Warmup number of QMC steps
sderiv  0.02      # Maximum derivative mismatch accepted for tail concatenation
aom  10           # Number of frequency points used to determin the value of sigma at nom
U    4.0
J    0.3
nf0  8.0
beta 50.0

If you have run CTQMC with this PARAMS file, you have in your local directory a new file named tvertex.dat. It has stored the two-particle Green’s functions for Nickel. I shall explain to you the data format now. So now we want to have the local Green’s functions and the non-local Green’s functions. These three quantities are inputs to the BSE. So for that, we want to run lmfdmft in a speecial mode called gprt:mode=3 to produce the loca and non-local Green’s functions.

Producing local and non-local Green’s functions:

cp Sig.out sig.inp
lmfdmft ni -vnk=12 --rs=1,0 --ldadc=28.95 -job=1 --gprt:mode=3:nom=100 --fullg

This produces two files named gloc.ni and gkloc.ni. Suscept.py, the script that solved the BSE and produces the dynamic momentum resolved Susceptibilities, reads them. It would also look for a file with the k-list and a file with the Q-list. The k-list is the same k-grid information on which gkloc.ni is written. And Q-list file contains the Q-grid on which you want to find the susceptibilities .

The input files and their paths are identified in the suscept.in. You can have the input files from here this link. The suscept.in file looks like below:

suscept.klist     # filename with k-list
suscept.qlist     # filename with Qlist
rmesh.dat      # real axis mesh
./G_k1r         # file with real axis k-dependent Gren's function
./G_local1r     # file with real axis local Gren's function
./chi0_real.     # name of the Bubble on real axis
gkloc.ni         # imaginary axis k-dependent Green's function
gloc.ni          # imaginary axis local Green's function
./tvertex.dat    # ctqmc local vertex function
50            # inverse temperature for bose function in Sq(omega)

Now, copy all essential files to a directory for Susceptibility calculations and launch Suscept.py to solve BSE:

mkdir suscept
cp *.klist *.qlist tvertex.dat gkloc.ni gloc.ni suscept.in suscept
cd suscept
Suscept.py suscept.in

You can find a good vertex from here and also other input files from this link this link:

mkdir good_vertex
cp ni-sus-inp/* good_vertex
cd good_vertex
cp /home/cseg/scarf705/ni-suscept/tvertex.dat .

Here are our cuprates results published in this paper . You should be able to reproduce the results with the input files and the script. You can copy a directory with all converged results for NCO and run the script to see if the results match with what we show in our paper.

cp -r /home/cseg/scarf705/*.tar.gz .