Work Flow

1. Construct a spectral density

As an examples, we consider a Ohmic form with exponential cutoff

\[J(\omega) = \pi\alpha\]

PowerLawExpSD(s::Float64, γ::Float64; alpha::Union{Float64,Nothing}=nothing, reorgene::Union{Float64,Nothing}=nothing) 

We need to specify either the coupling strength alpha or the reorganization energy reorgene. Note that in the Ohmic case ($s=1$), the reorganization energy equals to alpha.

    julia> using QFiND
    julia> sdens = PowLowExpSD(1.0, 50.0; alpha=35.0)
    (::PowerLawExpSD) (generic function with 1 method)

This sdens works as a function.

    julia> Er = reorganization_energy(sdens)
    35.0

1. Construct a QNSD

    julia> sbeta = BosonicQNSD(sdens, 300.0)
    (::BosonicQNSD) (generic function with 1 method)

1. Construct a bath correlation function

    julia> bcf = BosonicBCF(sbeta, 300.0)
    (::BosonicBCF) (generic function with 1 method)

1. Create an initial dataset