diff --git a/Source/Particles/PulsarParameters.cpp b/Source/Particles/PulsarParameters.cpp
index d8d943f678e..e9b13f41253 100644
--- a/Source/Particles/PulsarParameters.cpp
+++ b/Source/Particles/PulsarParameters.cpp
@@ -396,8 +396,8 @@ Pulsar::ReadParameters () {
         pp.get("beta_rec_RR",m_beta_rec_RR); // reconnection rate
         amrex::Real S = m_RLC/m_R_star; // Scale separation ratio
         m_BLC = m_B_star/(S*S*S); // Bfield at LC
-        m_re_scaled = m_beta_rec_RR * (3./2.) * PhysConst::m_e * PhysConst::c
-                    / (m_gammarad_RR * m_gammarad_RR * PhysConst::q_e * m_BLC);
+        m_re_scaled = m_beta_rec_RR * (3./2.) * PhysConst::q_e
+                    / ( 4 * MatchConst::pi * PhysConst::ep0 * m_gammarad_RR * m_gammarad_RR * PhysConst::c * m_BLC);
         amrex::Print() << " re_scaled at gammarad : " << m_gammarad_RR << " is : " << m_re_scaled << "\n";
         m_re_scaledratio = m_re_scaled / PhysConst::r_e;
         amrex::Print() << " ratio : re_scaled/re : " << m_re_scaledratio << "\n";