hamiltonians.hpp

Go to the documentation of this file.

// Copyright (c) 2025--present, The Simons Foundation
// This file is part of TRIQS/modest and is licensed under the terms of GPLv3 or later.
// SPDX-License-Identifier: GPL-3.0-or-later
// See LICENSE in the root of this distribution for details.
 
#pragma once
#include <triqs/gfs.hpp>
#include <triqs/operators.hpp>
#include <triqs/operators/many_body_operator.hpp>
 
using namespace triqs::gfs;
 
namespace triqs {
 
  std::pair<nda::matrix<double>, nda::matrix<double>> U_matrix_kanamori(long n_orb, double U_int, double U_prime, double J_hund);
 
  nda::array<double, 4> U_matrix_slater_spherical(long l, double U_int, double J_hund);
 
  nda::array<dcomplex, 4> U_matrix_slater_local(long l, nda::matrix<dcomplex> s2l, double U_int, double J_hund);
 
 
  using op_name = std::pair<std::string, int>;
 
  operators::many_body_operator rename_op(operators::many_body_operator const &op, std::map<op_name, op_name> const &map);
 
  struct kanamori_params {
    bool spin_flip    = true;
    bool pair_hopping = true;
  };
 
  operators::many_body_operator h_int_kanamori(nda::matrix<double> const &U, nda::matrix<double> const &Uprime, double const &J_hund,
                                               int const &n_orb, std::vector<std::string> const &spin_names, kanamori_params const &params = {});
 
  operators::many_body_operator h_int_density(nda::matrix<double> const &U, nda::matrix<double> const &Uprime, double const &J_hund, int const &n_orb,
                                              std::vector<std::string> const &spin_names);
 
  operators::many_body_operator h_int_slater(nda::array<dcomplex, 4> const &U, int const &n_orb, std::vector<std::string> const &spin_names);
 
  operators::many_body_operator make_density_density(std::vector<std::string> const &tau_names, std::vector<long> const &dim_gamma, double U_int,
                                                     double U_prime, double J_hund);
 
  operators::many_body_operator make_kanamori(std::vector<std::string> const &tau_names, std::vector<long> const &dim_gamma, double U_int,
                                              double U_prime, double J_hund, bool spin_flip = true, bool pair_hopping = true);
  operators::many_body_operator make_slater(std::vector<std::string> const &tau_names, std::vector<long> const &dim_gamma, double U_int,
                                            double J_hund, nda::matrix<dcomplex> const &spherical_to_dft,
                                            std::optional<nda::matrix<dcomplex>> const &dft_to_local);
 
} // namespace triqs