Julia code to generate finite temperature chain coefficients (on-site energy, hopping energy and coupling constant) of a chain-mapped bosonic environment [1,2].

Parts of this code is a translation in Julia from a previous Matlab code [3,4].

[1]Chin, A.W. et al. (2010) ‘Exact mapping between system-reservoir quantum models and semi-infinite discrete chains using orthogonal polynomials’, Journal of Mathematical Physics, 51(9), p. 092109. doi:https://doi.org/10.1063/1.3490188.

[2]Tamascelli, D. et al. (2019) ‘Efficient Simulation of Finite-Temperature Open Quantum Systems’, Physical Review Letters, 123(9), p. 090402. doi:https://doi.org/10.1103/PhysRevLett.123.090402.

[3] Gautschi, W. (1994) ‘Algorithm 726: ORTHPOL–a package of routines for generating orthogonal polynomials and Gauss-type quadrature rules’, ACM Transactions on Mathematical Software, 20(1), pp. 21–62. doi:https://doi.org/10.1145/174603.174605.

[4] Gautschi, W. (2005) ‘Orthogonal polynomials (in Matlab)’, Journal of Computational and Applied Mathematics, 178(1–2), pp. 215–234.