The following limit is left unevaluated (Edit: added the assumption that $\epsilon$ is real thanks to the comment below):
Limit[ 1/(1 + Exp[ϵ/T]), T->0, Direction->-1, Assumptions->ϵ ∈ Reals]
However, the limit is easy to calculate (it is the Fermi-Dirac distribution), and it gives $\theta(-\epsilon)$. It seems that MMA has a hard time when the result is discontinuous in the parameters.
Nevertheless, MMA can actually calculate the limit, if we help it by explicitly splitting up the parameter space:
In[2]:= Limit[1/(1 + Exp[ϵ/T]), T->0, Direction->-1, Assumptions->#] & /@ {ϵ > 0, ϵ < 0}
Out[2]= {0, 1}
Is there a way to make MMA automatically calculate such limits, without explicitly splitting up the parameter space?
Assuming[]
or theAssumptions
option)... $\endgroup$HeavisideTheta[]
andUnitStep[]
have a different definition from your $\theta(u)$; your limit evaluates to $\dfrac12$ for $\varepsilon=0$, for instance. $\endgroup$Limit[]
-related questions cropping up, I don't see the need); at the moment, calculus-and-analysis is sufficient to cover your question. $\endgroup$Limit[1/(1 + Exp[x/T]), T -> 0, Direction -> #, Assumptions -> x > 0] & /@ {1, -1}
$\endgroup$