lime/heom/heom.py
def _calc_matsubara_params(N_exp, coup_strength, cut_freq, temperature):
    """
    Calculate the Matsubara coefficents and frequencies

    Returns
    -------
    c, nu: both list(float)

    """
    c = []
    nu = []
    lam0 = coup_strength
    gam = cut_freq
    hbar = 1.
    beta = 1.0/temperature
    N_m = N_exp

    g = 2*np.pi / (beta*hbar)
    for k in range(N_m):
        if k == 0:
            nu.append(gam)
            c.append(lam0*gam*
                (1.0/np.tan(gam*hbar*beta/2.0) - 1j) / hbar)
        else:
            nu.append(k*g)
            c.append(4*lam0*gam*nu[k] /
                  ((nu[k]**2 - gam**2)*beta*hbar**2))

#    exp_coeff = c
#    exp_freq = nu
    return c, nu

    Calculate the Matsubara coefficents and frequencies

    Returns
    -------
    c, nu: both list(float)

    
Calculate the Matsubara coefficents and frequencies
Returns


Calculate the Matsubara coefficents and frequencies Returns
def _calc_matsubara_params(N_exp, coup_strength, cut_freq, temperature):
    c = []
    nu = []
    lam0 = coup_strength
    gam = cut_freq
    hbar = 1.
    beta = 1.0/temperature
    N_m = N_exp
    g = 2*np.pi / (beta*hbar)
    for k in range(N_m):
        if k == 0:
            nu.append(gam)
            c.append(lam0*gam*
                (1.0/np.tan(gam*hbar*beta/2.0) - 1j) / hbar)
        else:
            nu.append(k*g)
            c.append(4*lam0*gam*nu[k] /
                  ((nu[k]**2 - gam**2)*beta*hbar**2))
    return c, nu
def _calc_matsubara_params ( N_exp , coup_strength , cut_freq , temperature ) : c = [ ] nu = [ ] lam0 = coup_strength gam = cut_freq hbar = 1. beta = 1.0 / temperature N_m = N_exp g = 2 * np . pi / ( beta * hbar ) for k in range ( N_m ) : if k == 0 : nu . append ( gam ) c . append ( lam0 * gam * ( 1.0 / np . tan ( gam * hbar * beta / 2.0 ) - 1j ) / hbar ) else : nu . append ( k * g ) c . append ( 4 * lam0 * gam * nu [ k ] / ( ( nu [ k ] ** 2 - gam ** 2 ) * beta * hbar ** 2 ) ) return c , nu
Calculate the Matsubara coefficents and frequencies
Returns
Calculate the Matsubara coefficents and frequencies Returns
"""
    Calculate the Matsubara coefficents and frequencies

    Returns
    -------
    c, nu: both list(float)

    """ #    exp_coeff = c #    exp_freq = nu
N_exp
coup_strength
cut_freq
temperature