diff --git a/src/run.jl b/src/run.jl
index 324ae12..32cdc32 100644
--- a/src/run.jl
+++ b/src/run.jl
@@ -81,7 +81,7 @@ function run_sim(
         n_snapshots=n_snapshots,
         c₁=4 * sim_consts.ϵ * 6 * sim_consts.σ^6 * sim_consts.δt * sim_consts.μ,
         c₂=2 * sim_consts.σ^6,
-        c₃=sqrt(2 * sim_consts.D₀ * sim_consts.δt),
+        c₃=sqrt(2 * sim_consts.Dₜ * sim_consts.δt),
         c₄=sqrt(2 * sim_consts.Dᵣ * sim_consts.δt),
         v₀δt=sim_consts.v₀ * sim_consts.δt,
         μ=sim_consts.μ,
diff --git a/src/setup.jl b/src/setup.jl
index 818b11f..6cf7bdd 100644
--- a/src/setup.jl
+++ b/src/setup.jl
@@ -55,10 +55,10 @@ function gen_sim_consts(
     @assert packing_ratio > 0
 
     μ = 1.0
-    D₀ = 1.0
+    Dₜ = 1.0
     particle_radius = 0.5
-    Dᵣ = 3 * D₀ / ((2 * particle_radius)^2)
-    σ = 1.0
+    Dᵣ = 3 * Dₜ / ((2 * particle_radius)^2)
+    σ = 2 * particle_radius
     ϵ = 100.0
     interaction_r = 2^(1 / 6) * σ
 
@@ -92,7 +92,7 @@ function gen_sim_consts(
         packing_ratio,
         # Internal
         μ,
-        D₀,
+        Dₜ,
         particle_radius,
         Dᵣ,
         σ,