TRIQS Hybridization Expansion Solver
  • Install CTHYB
    • Packaged Versions of CTHYB
      • Ubuntu Debian packages
      • Anaconda
      • Docker
    • Compiling CTHYB from source
      • Prerequisites
      • Installation steps
      • Version compatibility
      • Custom CMake options
  • Documentation
    • Basic notions
      • A word on the algorithm
      • QMC moves
        • Insert one pair of operators
        • Remove one pair of operators
        • Insert two pairs of operators
        • Remove two pairs of operators
        • Shift one operator
        • Global move - change of operator indices
      • Measurements: definitions
        • Green’s function
          • Imaginary time
          • Legendre polynomial basis
        • Local susceptibility and other operator pairs
        • Two-particle Green’s functions
          • Imaginary time binning
          • Matsubara frequency measurements
          • Mixed Matsubara Frequency and Legendre measurements
        • Impurity density matrix
        • Perturbation order histograms
        • Average sign
    • User guide
      • Setting the parameters
        • Step 1 - construct the solver instance
          • Examples
        • Step 2 - the Hamiltonian
          • Examples
        • Step 3 - the Monte Carlo parameters
        • Step 6 - Legendre or not?
        • Step 7 - prepare the non-interacting Green’s function
        • Step 8 - we’re ready to go!
        • Final Step - analyze the output
      • Building DMFT calculations
        • Plain-vanilla DMFT: the Bethe lattice
      • Random number generators
      • Measuring static observables / impurity density matrix
        • Measuring \(\hat{\rho}_\mathrm{imp}\) in cthyb
        • Compute expactation values of static observables
        • Measuring moments of Green’s function and self-energy
        • Evaluation of the impurity interaction energy
      • Multiplet analysis & particle number histograms
        • Analysis of results
          • plot particle number histogram
        • More complex example
      • Dynamical spin-spin susceptibility \(\chi_{S_z S_z}(\tau)\)
      • Perturbation order histograms
        • Plotting perturbation order histograms
      • Tail fitting via the CRM Dyson solver
        • Setup
        • Using the Constrained Residual Minimization Dyson solver
    • Tutorials
      • The Anderson impurity model
      • Multiorbital impurity model
      • Convergence test of CTHYB solver parameters
        • Overview
        • 1. length_cycle
        • 2. n_warmup_cycles
        • 3. n_cycles
        • 4. n_iw
        • 5. n_l
        • 6. n_tau
      • High Frequency Moments of the Green’s Functions and Self-energy
        • Derivation
        • Implentation in TRIQS/cthyb
    • Reference manual
      • triqs_cthyb.configuration
        • triqs_cthyb.configuration.Configuration
          • Configuration
          • triqs_cthyb.configuration.Configuration.beta
          • triqs_cthyb.configuration.Configuration.clear
          • triqs_cthyb.configuration.Configuration.erase
          • triqs_cthyb.configuration.Configuration.finalize
          • triqs_cthyb.configuration.Configuration.get_id
          • triqs_cthyb.configuration.Configuration.insert
          • triqs_cthyb.configuration.Configuration.replace
        • triqs_cthyb.configuration.OpDesc
          • OpDesc
          • triqs_cthyb.configuration.OpDesc.block_index
          • triqs_cthyb.configuration.OpDesc.dagger
          • triqs_cthyb.configuration.OpDesc.inner_index
          • triqs_cthyb.configuration.OpDesc.linear_index
      • triqs_cthyb.multiplet_tools
        • triqs_cthyb.multiplet_tools.multiplet_analysis
          • multiplet_analysis()
      • triqs_cthyb.solver
        • triqs_cthyb.solver.Solver
          • Solver
          • triqs_cthyb.solver.Solver.Delta_infty
          • triqs_cthyb.solver.Solver.Delta_tau
          • triqs_cthyb.solver.Solver.G0_iw
          • triqs_cthyb.solver.Solver.G2_iw
          • triqs_cthyb.solver.Solver.G2_iw_nfft
          • triqs_cthyb.solver.Solver.G2_iw_ph
          • triqs_cthyb.solver.Solver.G2_iw_ph_nfft
          • triqs_cthyb.solver.Solver.G2_iw_pp
          • triqs_cthyb.solver.Solver.G2_iw_pp_nfft
          • triqs_cthyb.solver.Solver.G2_iwll_ph
          • triqs_cthyb.solver.Solver.G2_iwll_pp
          • triqs_cthyb.solver.Solver.G2_tau
          • triqs_cthyb.solver.Solver.G_l
          • triqs_cthyb.solver.Solver.G_tau
          • triqs_cthyb.solver.Solver.G_tau_accum
          • triqs_cthyb.solver.Solver.O_tau
          • triqs_cthyb.solver.Solver.asymmetry_G_tau
          • triqs_cthyb.solver.Solver.auto_corr_time
          • triqs_cthyb.solver.Solver.auto_corr_time_converged
          • triqs_cthyb.solver.Solver.average_order
          • triqs_cthyb.solver.Solver.average_sign
          • triqs_cthyb.solver.Solver.constr_parameters
          • triqs_cthyb.solver.Solver.density_matrix
          • triqs_cthyb.solver.Solver.h_loc
          • triqs_cthyb.solver.Solver.h_loc0
          • triqs_cthyb.solver.Solver.h_loc_diagonalization
          • triqs_cthyb.solver.Solver.hybridisation_is_complex
          • triqs_cthyb.solver.Solver.last_configuration
          • triqs_cthyb.solver.Solver.last_constr_parameters
          • triqs_cthyb.solver.Solver.last_solve_parameters
          • triqs_cthyb.solver.Solver.local_hamiltonian_is_complex
          • triqs_cthyb.solver.Solver.performance_analysis
          • triqs_cthyb.solver.Solver.perturbation_order
          • triqs_cthyb.solver.Solver.perturbation_order_total
          • triqs_cthyb.solver.Solver.solve_parameters
          • triqs_cthyb.solver.Solver.solve_status
          • triqs_cthyb.solver.Solver.solve
      • triqs_cthyb.solver_core
        • triqs_cthyb.solver_core.ConstrParametersT
          • ConstrParametersT
          • triqs_cthyb.solver_core.ConstrParametersT.beta
          • triqs_cthyb.solver_core.ConstrParametersT.delta_interface
          • triqs_cthyb.solver_core.ConstrParametersT.gf_struct
          • triqs_cthyb.solver_core.ConstrParametersT.n_iw
          • triqs_cthyb.solver_core.ConstrParametersT.n_l
          • triqs_cthyb.solver_core.ConstrParametersT.n_tau
        • triqs_cthyb.solver_core.SolveParametersT
          • SolveParametersT
          • triqs_cthyb.solver_core.SolveParametersT.det_init_size
          • triqs_cthyb.solver_core.SolveParametersT.det_n_operations_before_check
          • triqs_cthyb.solver_core.SolveParametersT.det_precision_error
          • triqs_cthyb.solver_core.SolveParametersT.det_precision_warning
          • triqs_cthyb.solver_core.SolveParametersT.det_singular_threshold
          • triqs_cthyb.solver_core.SolveParametersT.h_int
          • triqs_cthyb.solver_core.SolveParametersT.h_loc0
          • triqs_cthyb.solver_core.SolveParametersT.imag_threshold
          • triqs_cthyb.solver_core.SolveParametersT.initial_configuration
          • triqs_cthyb.solver_core.SolveParametersT.length_cycle
          • triqs_cthyb.solver_core.SolveParametersT.loc_n_max
          • triqs_cthyb.solver_core.SolveParametersT.loc_n_min
          • triqs_cthyb.solver_core.SolveParametersT.max_time
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_block_order
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_blocks
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_iw
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_iw_nfft
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_iw_ph
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_iw_ph_nfft
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_iw_pp
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_iw_pp_nfft
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_iwll_nfft_buf_size
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_iwll_ph
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_iwll_pp
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_n_bosonic
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_n_fermionic
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_n_l
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_n_tau
          • triqs_cthyb.solver_core.SolveParametersT.measure_G2_tau
          • triqs_cthyb.solver_core.SolveParametersT.measure_G_l
          • triqs_cthyb.solver_core.SolveParametersT.measure_G_tau
          • triqs_cthyb.solver_core.SolveParametersT.measure_O_tau
          • triqs_cthyb.solver_core.SolveParametersT.measure_O_tau_min_ins
          • triqs_cthyb.solver_core.SolveParametersT.measure_density_matrix
          • triqs_cthyb.solver_core.SolveParametersT.measure_pert_order
          • triqs_cthyb.solver_core.SolveParametersT.move_double
          • triqs_cthyb.solver_core.SolveParametersT.move_global
          • triqs_cthyb.solver_core.SolveParametersT.move_global_prob
          • triqs_cthyb.solver_core.SolveParametersT.move_shift
          • triqs_cthyb.solver_core.SolveParametersT.n_cycles
          • triqs_cthyb.solver_core.SolveParametersT.n_warmup_cycles
          • triqs_cthyb.solver_core.SolveParametersT.nfft_buf_sizes
          • triqs_cthyb.solver_core.SolveParametersT.off_diag_threshold
          • triqs_cthyb.solver_core.SolveParametersT.partition_method
          • triqs_cthyb.solver_core.SolveParametersT.performance_analysis
          • triqs_cthyb.solver_core.SolveParametersT.proposal_prob
          • triqs_cthyb.solver_core.SolveParametersT.quantum_numbers
          • triqs_cthyb.solver_core.SolveParametersT.random_name
          • triqs_cthyb.solver_core.SolveParametersT.random_seed
          • triqs_cthyb.solver_core.SolveParametersT.use_norm_as_weight
          • triqs_cthyb.solver_core.SolveParametersT.use_trace_estimator
          • triqs_cthyb.solver_core.SolveParametersT.verbosity
        • triqs_cthyb.solver_core.SolverCore
          • SolverCore
          • triqs_cthyb.solver_core.SolverCore.Delta_infty
          • triqs_cthyb.solver_core.SolverCore.Delta_tau
          • triqs_cthyb.solver_core.SolverCore.G0_iw
          • triqs_cthyb.solver_core.SolverCore.G2_iw
          • triqs_cthyb.solver_core.SolverCore.G2_iw_nfft
          • triqs_cthyb.solver_core.SolverCore.G2_iw_ph
          • triqs_cthyb.solver_core.SolverCore.G2_iw_ph_nfft
          • triqs_cthyb.solver_core.SolverCore.G2_iw_pp
          • triqs_cthyb.solver_core.SolverCore.G2_iw_pp_nfft
          • triqs_cthyb.solver_core.SolverCore.G2_iwll_ph
          • triqs_cthyb.solver_core.SolverCore.G2_iwll_pp
          • triqs_cthyb.solver_core.SolverCore.G2_tau
          • triqs_cthyb.solver_core.SolverCore.G_l
          • triqs_cthyb.solver_core.SolverCore.G_tau
          • triqs_cthyb.solver_core.SolverCore.G_tau_accum
          • triqs_cthyb.solver_core.SolverCore.O_tau
          • triqs_cthyb.solver_core.SolverCore.asymmetry_G_tau
          • triqs_cthyb.solver_core.SolverCore.auto_corr_time
          • triqs_cthyb.solver_core.SolverCore.auto_corr_time_converged
          • triqs_cthyb.solver_core.SolverCore.average_order
          • triqs_cthyb.solver_core.SolverCore.average_sign
          • triqs_cthyb.solver_core.SolverCore.constr_parameters
          • triqs_cthyb.solver_core.SolverCore.density_matrix
          • triqs_cthyb.solver_core.SolverCore.h_loc
          • triqs_cthyb.solver_core.SolverCore.h_loc0
          • triqs_cthyb.solver_core.SolverCore.h_loc_diagonalization
          • triqs_cthyb.solver_core.SolverCore.hybridisation_is_complex
          • triqs_cthyb.solver_core.SolverCore.last_configuration
          • triqs_cthyb.solver_core.SolverCore.last_constr_parameters
          • triqs_cthyb.solver_core.SolverCore.last_solve_parameters
          • triqs_cthyb.solver_core.SolverCore.local_hamiltonian_is_complex
          • triqs_cthyb.solver_core.SolverCore.performance_analysis
          • triqs_cthyb.solver_core.SolverCore.perturbation_order
          • triqs_cthyb.solver_core.SolverCore.perturbation_order_total
          • triqs_cthyb.solver_core.SolverCore.solve_parameters
          • triqs_cthyb.solver_core.SolverCore.solve_status
          • triqs_cthyb.solver_core.SolverCore.solve
      • triqs_cthyb.tail_fit
        • triqs_cthyb.tail_fit.green_high_frequency_moments
          • green_high_frequency_moments()
        • triqs_cthyb.tail_fit.sigma_high_frequency_moments
          • sigma_high_frequency_moments()
      • triqs_cthyb.util
        • triqs_cthyb.util.block_size_from_gf_struct
          • block_size_from_gf_struct()
        • triqs_cthyb.util.estimate_nfft_buf_size
          • estimate_nfft_buf_size()
        • triqs_cthyb.util.orbital_occupations
          • orbital_occupations()
    • FAQs
      • Frequently-asked questions
        • Q: Why is my code running so slowly?
        • Q: How do I save the triqs and cthyb hash and script for debugging purposes?
        • Q: Why does my data look so noisy?
        • Q: How do I use the segment picture?
  • Reporting issues
  • Changelog
    • Version 4.0.0
      • General
      • doc
      • cmake
      • jenkins
    • Version 3.3.0
      • General
      • doc
      • jenkins
    • Version 3.2.1
      • General
      • doc
      • jenkins
    • Version 3.2.0
      • General
      • doc
      • fix
    • Version 3.1.1
      • General
      • Doc
      • Numpy
    • Version 3.1.0
      • Delta Interface
      • General
      • jenkins
      • doc
      • cmake
      • gh-actions
    • Version 3.0.0
      • Fundamental Green Function Symmetries
      • Dependency Management
      • Other Changes
    • Version 2.2.1
      • doc
      • General
    • Version 2.2.0
      • cmake
      • doc
      • General
      • h5
      • packaging
    • Version 2.1.0
      • Contributors
      • Testers
      • Changes in installation and cmake files
      • Solver Interface
      • Changes in behavior
      • Dependencies
  • Authors
    • Citation
    • License
    • Disclaimer
TRIQS Hybridization Expansion Solver
  • Documentation
  • View page source

Documentation

Basic notions

  • A word on the algorithm
  • QMC moves
  • Measurements: definitions

User guide

  • Setting the parameters
  • Building DMFT calculations
  • Random number generators
  • Measuring static observables / impurity density matrix
  • Multiplet analysis & particle number histograms
  • Dynamical spin-spin susceptibility \(\chi_{S_z S_z}(\tau)\)
  • Perturbation order histograms
  • Tail fitting via the CRM Dyson solver

Tutorials

  • The Anderson impurity model
  • Multiorbital impurity model
  • Convergence test of CTHYB solver parameters
  • High Frequency Moments of the Green’s Functions and Self-energy

Reference manual

triqs_cthyb.configuration

MC configuration for triqs_cthyb.

triqs_cthyb.multiplet_tools

functions for analyzing the multiplet structure in cthyb

triqs_cthyb.solver

the triqs_cthyb solver class

triqs_cthyb.solver_core

The TRIQS cthyb solver

triqs_cthyb.tail_fit

tail fitting and high frequency moments

triqs_cthyb.util

utility functions

Links to all relevant solver parameters:

  • Construction parameters

  • Solve parameters

The C++ reference manual can be found here.

FAQs

  • Frequently-asked questions
    • Q: Why is my code running so slowly?
    • Q: How do I save the triqs and cthyb hash and script for debugging purposes?
    • Q: Why does my data look so noisy?
    • Q: How do I use the segment picture?
Previous Next

© Copyright 2014-2025, The Simons Foundation, P. Seth, I. Krivenko, H. Strand, M. Ferrero, O. Parcollet.