QuaCaTeam / quaca

QuaCa is an open-source library which provides fast calculations of nonequilibrium atom-surface interactions.
https://quacateam.github.io/quaca/
MIT License
3 stars 2 forks source link

ReflectionCoefficientLocBulk with PermittivityDrude blows up at omega = 0 #49

Closed SiMaHer closed 4 years ago

SiMaHer commented 4 years ago

The reflection coefficient returns nan's for a frequency of oemga = 0. This is due the fact, that the PermittivityDrude.calculate returns a value of (-inf, nan), leading to a return value of (-nan, -nan) in ReflectionCoefficientLocBulk.calculate

myoelmy commented 4 years ago

Used scaled permittivities to fix the bug.