Hi,

May I know on what basis the values of the bath parameters (epsilon and V) are chosen?

# Levels of the bath sites
epsilon = [-1.0, 1.0]
# Hopping amplitudes
V = [0.5, 0.5]

Are these parameters optimized during the calculations or not?

Another inequality please:
Is there an available tutorial for using pomerol2triqs with DMFT calculations?

Best,
Ragab.

Dear Igor,

I have spotted an oddity in the PP-channel two particle Green's function, when benchmarking triqs/cthyb. When asking for 15x15x15 frequencies the shape of both the PP and PH gfs are:

gpp.data.shape == (29, 30, 30, 1, 1, 1, 1)

However the zeroth bosonic frequency (the first index) is located at index 13 for the PP-gf and at index 14 for the PH-gf. Please see the triqs/cthyb benchmark and the PP channel comparison and compare with the PH case:

https://github.com/TRIQS/cthyb/blob/dev_measures/benchmark/measure_g2/hubbard_atom_2bath/CompareEDandCTHYB.ipynb

Is this frequency shift intended or a bug? Or have I messed something up?

Best regards,
Hugo

I have successfully used TRIQS on Ubuntu 22.04.2, and now I am attempting to install pomerol2triqs. However, I am encountering difficulties in locating the path mentioned below:

<path_to_triqs_install_dir>/share/triqs/triqsvars.sh

I can't find triqsvars.sh anywhere.

Could you please provide some guidance or assistance?

Dear Igor,

I am not sure whether to post this question at pomerol2triqs or pomerol. I am calculating two particle Green's functions and it is very slow. I think I am using some mpi parallelism since pomerol says it is distributing work over "processors".

My question is whether I somehow can convince pomerol to compute the different two-particle Green's function elements in parallel? (I have a 4x4x4x4 = 256 elements to calculate in total, so that would be a nice speedup if done in parallel)

Is this possible to do this out of the box with pomerol? If not, what would I have to do to accomplish this?

Cheers, Hugo

Dear Igor,

I have looked through the example scripts but I can not find any example on how got compute static expectation values. Is there such functionality? or do I have to take the \beta=0 limit of a single-particle Green's function.

Cheers, Hugo

When using CMAKE_CXX_STANDARD 17, I think 3.8 is the minimum required version of cmake. There were some commits along these lines at some point, but they were reverted again, and the CXX_STANDARD was increased later leading to an inconsistency. This is what CMakeLists.txt currently looks like:

Start configuration

cmake_minimum_required(VERSION 3.1 FATAL_ERROR)
project(pomerol2triqs CXX)
set(CMAKE_CXX_STANDARD 17)

Best,

Erik

Dear Igor,

At present, only Sz is used for state classification as in the default mode of pomerol. In the Slater Hamiltonian, for example, Lz (w/o SO) or Jz (w/ SO) could be additionally used to reduce the block size.

Is there any way to pass an Operator op to call symm.checkSymmetry(op)? Or, if not, which is the easiest way to extend the code?

It seems for me that one should extend the function pomerol_ed::diagonalize to be able to pass a list of operators that is used for additional symmetry check.

Best regards,
Junya

Hi,

Is pomerol2triqs compatible with triqs 2.2? I am having some issues compiling it. For example, triqs/mpi/boost.hpp no longer exists in triqs 2.2.

Best,

Erik