holstein mobility is 1 / lattice_constant ^2

jarvist · Apr 12, 2024 · b8ecc52 · b8ecc52
1 parent 4ef156c
commit b8ecc52
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/PolaronUnits.jl b/src/PolaronUnits.jl
@@ -230,7 +230,7 @@ function addunits!(polaron::HolsteinPolaron; unit="pu")
     polaron.Ma = pustrip.(polaron.Ma .* polaron.mb) .* punit(1Unitful.me)
     polaron.Mr = pustrip.(polaron.Mr .* polaron.mb) .* punit(1Unitful.me)
     polaron.R = pustrip.(polaron.R ./ sqrt(2 * polaron.ωeff * polaron.mb)) .* punit(1Unitful.Å)
-    polaron.μ = pustrip.(polaron.μ ./ polaron.mb) .* punit(1u"cm^2/V/s")
+    polaron.μ = pustrip.(polaron.μ ./ polaron.a.^2) .* punit(1u"cm^2/V/s")
     polaron.χ = pustrip.(polaron.χ .* polaron.mb) .* punit(u"Ω") .* punit(1Unitful.THz2π)
     polaron.z = pustrip.(polaron.z .* polaron.mb) .* punit(u"Ω")
     polaron.σ = pustrip.(polaron.σ ./ polaron.mb) .* punit(u"S")