Show Posts

This section allows you to view all posts made by this member. Note that you can only see posts made in areas you currently have access to.


Topics - Leonardo

Pages: [1]
1
Dear Colleagues,
I need your help in understanding an observation: Take any kind of small system (e.g. a small cluster of atoms or molecule). You can calculate the Hesienberg exchange coupling Jij(R) in two different ways: either as a "molecular configuration" or as a "bulk configuration" (using a large enough computational cell to avoid image interactions). For the latter, QATK together with your script "plot_Jij_vs_r.py" provides the coupling values Jij for R values corresponding to the interatomic distances. Doing exactly the same in the "molecular configuration" gives exactly the same Jij(R>0) values, BUT it gives an additional value for R=0! It turns out that Jij(R=0) is in general much stronger than all the other Jij values.
Could you please explain me the physical meaning for Jij(R=0)? Is it a measure for the coupling strength of two spins within the same atom? But why does this term not show up in the "bulk configuration" calculation?
With many sincere thanks for any hint in advance!

PS: The script "plot_Jij_vs_r.py" proposed in the QATK documentation is the same for QATK 2019 and the latest versions of QATK. However, this script provides much more convincing results when used to analyze the results obtained with the old QATK 2019 then with those given by QATK 2021. Is there a chance that the script "plot_Jij_vs_r.py" was not properly updated to reliably analyze the results given by the newer versions of QATK? Who is the author of "plot_Jij_vs_r.py"?

2
Dear Colleagues,

I'm trying to use QATK for the first time and I would like to calculate the absorption spectrum of the Si20H20 fullerene (see attached Si20H20_OPT.py).

The calculation runs without any problem, but to my big surprise the off-diagonal elements of the absorption tensor turn out to be about 4 times stronger than the diagonal ones.

This looks quite strange to me, but I have not found a way to diagonalize the tensor.

Any help or hint would be highly appreciated.
With many thanks in advance!

Leonardo

Pages: [1]