Hello Prof. Gull,


> Moreover, I tried to follow the instruction in doc/hybdoc.pdf to run hybridization from scratch, but there's a file "delta.dat" containing hybridization functions which I don't know how to produce. Did I miss something in the tutorial?

I can help you with that. For your application, the python scripts may not be the best way to proceed.

My intention was simply using the python scripts to check if I could reproduce Figs 1 & 2 in Comp. Phy. Comm. 184, 1280 (2013) for the single-orbital Anderson model as the first step, so tutorial hybridization-02 is essential to me. But I agree python scripts may not be the best way for my project (see below).
 

> ps. I'm aware that CT-HYB is implemented as an impurity solver for the DMFT applications. But my project does not need DMFT calculations and I intend to simply use hybridization as a stand-alone module.

In order to simulate an Anderson impurity model you will need to have a hybridization function, parametrized in imaginary time. We can set one up if you’re interested, but ideally it comes from your physics setup.

I was looking for the format of this dat file but it seems not mentioned in the documentation. However, I just noticed that a delta.dat file with constant hybridization function could be generated from the script of tutorial hybridization-01 (after skipping the cthyb and boost.mpi errors). I assume the first column is for the imaginary time index and the second and the third for spin up and down, respectively. Is it correct? By feeding this dat to the hybridization executable, I successfully finished a calculation with a bunch of dat file and a h5 file generated. 

My problem at hand involves one quantum dot coupled to two metal leads, so basically it's kind of like a single-orbital Anderson impurity model. Later I may need to increase the number of impurities and to include interactions with other physical degree of freedoms, which may complicate the update scheme. I guess eventually I may have to tweak the hybridization executable to fit my need.

Thank you.


Sincerely,
Leo (Fang Yao-Lung)
Duke Physics