NMinimize[{10^18 EtotalEM[10^-9 R, 10^-9 Δ],
   {Δ <= R/ ArcSinh[1], Δ > 0, Δ < 100, R > 5, R < 10^4}},
   {Δ, R}