Note: The method below is focused mainly on the case where $N$ is integer (which OP confirms is the case in the comment below)
You can get a symbolic function/expression if you solve for $\sigma$ with the caveat that you will have to consider a finite set of explicit integers.
Notice that for a fixed $N$, say $N=3$, you have a polynomial equation if you change variables and rearrange the equation like this :
f[n, s] - ϵ /. s -> Sqrt[2] Sqrt[-Log[z]] /. n -> 3 //
Together // Numerator // Simplify
$$-z^5-z^4+2 z^3-3 z (\epsilon +1)+3 \epsilon$$
The above equation should equate to zero.
You do not have to use a change of variables however, you can solve the equation directly.
Here is a function to generate nmax solutions:
table[nmax_] :=
Table[{n -> m,
Refine[Solve[{f[m, s] == ϵ, s > 0}, s,
Reals], ϵ > 0]}, {m, nmax}]
The three first solution for $N\in\{1,2,3\} $:
Important note ! : $\sigma$ positive below is left as a means to simplify the discussion below but you should remove it in your code. See important footnote in bold at the end.
table[4]
$$\left(
\begin{array}{cc}
n\to 1 & \left\{\left\{s\to \sqrt{2} \sqrt{\log \left(\frac{\epsilon +1}{\epsilon }\right)}\right\}\right\} \\
n\to 2 & \left\{\left\{s\to \sqrt{2} \sqrt{\log \left(\frac{\sqrt{4 \epsilon ^2+12 \epsilon +1}+2 \epsilon +1}{4 \epsilon }\right)}\right\}\right\} \\
n\to 3 & \left\{\left\{s\to \sqrt{2} \sqrt{\log \left(\text{Root}\left[3 \text{$\#$1}^5 \epsilon +\text{$\#$1}^4 (-3 \epsilon -3)+2 \text{$\#$1}^2-\text{$\#$1}-1\&,1\right]\right)}\right\}\right\} \\
\end{array}
\right)$$
The first two are explicit and the third involves the root of a quintic which can not be solved in terms of radicals but it can be solved in terms of hypergeometric functions as shown in this Wikipedia page : https://en.wikipedia.org/wiki/Bring_radical#Other_characterizations
The next root involves the solution to a polynomial of degree 10. Maybe one can find an explicit solution but for the question here you do not need to and in general functions based from Root
are full fledged symbolic functions for which you can find series expansions at small and large arguments, locate singularities symbolically with FunctionSingularities
and much more. You just might have trouble with integration which such functions but that is not the question here.
Consider now a value for $\epsilon$ such as $\epsilon=0.4$
table[4] /. ϵ -> 0.4
(* {{n -> 1, {{s -> 1.58289}}}, {n -> 2, {{s -> 1.41233}}}, {n -> 3, {{s -> 1.56111}}}, {n -> 4, {{s -> 1.58272}}}} *)
We can plot s
=$\sigma$ for 10 values of n
=$N$ with $\epsilon=0.4$
Note: Takes a while you might want to store the generated list if you plan to plot a lot with a fixed value of nmax
table[10] // Map[{n, s} /. {#[[1]], #[[2, 1, 1]]} &] //
ReplaceAll[\[Epsilon] -> 0.4] // ListPlot

We can also plot $\sigma$=s
for different $N$ as a function of $\epsilon$ where each curve is labeled by the value of $N$:
table[4] // Map[{s, n} /. {#[[1]], #[[2, 1, 1]]} &] //
MapApply[Callout] // Plot[#, {ϵ, 0, 0.5}] &

Important footnote
I imagine that s
=$\sigma$>0 imposes that the solution is positive for any $\epsilon$. In that case, one is excluding solutions that are positive or negative depending on the value of $\epsilon$. To be safe one should perhaps remove the $\sigma$>0 constraint above.
n
notN
because capitalN
is a builtin. In any case, this equation is too complex to solve usingSolve
. I'm assuming you want a fully symbolic solution and you haven't got numerical values forS
etc. ? $\endgroup$