So in the past few months when trying to learn about the properties of the fixed points in ordinals as I move from $0$ to $\epsilon_{\epsilon_0}$ I noticed when moving from $\epsilon_n$ to the next one $\epsilon_{n+1}$, the operation that does that is a left associative operation i.e.
$${\tiny⋰}^{\epsilon_n^{\epsilon_n}}={}^{\omega}(\epsilon_n)=\epsilon_{n+1}$$
Since the exact same rule holds for all epsilons to those that are indexed by successor ordinals, unwrapping the whole thing gives $\epsilon_{n},n < \omega$ in terms of $\epsilon_0$:
$$\epsilon_n=\underbrace{{}^{\omega}(\cdots{}^{\omega}({}^{\omega}(}_{\textrm{n times}}\epsilon_0)))=\underbrace{{}^{\omega}(\cdots{}^{\omega}({}^{\omega}(}_{\textrm{n-1 times}}\epsilon_1)))=\underbrace{{}^{\omega}(\cdots{}^{\omega}({}^{\omega}(}_{\textrm{n-2 times}}\epsilon_2)))=etc.$$
That means, while $\epsilon_0$ is informally speaking the same as ${}^{\omega}\omega=\omega^{\omega^⋰}$ hence $\omega [5] 2=\omega\uparrow^3 2$, a pentation, all subsequent epsilon numbers are not tetration because the iterative operations that define them recursively is left associative but not right associative (as it would be the case for hyperoperators).
This then prompt the question: What prevent us from defining ordinal tetration, do they end up collapsing into $\epsilon_0$ or related terms thus making them unnecessary, or something else happens
To begin, since $\epsilon_0=\omega[5]2$, the next probable candidate for $\omega [5]3$ will be ${}^{{}^{\omega}}{}^{\omega}\omega=\left(\left({}^{({}^{\omega})}{}^{\omega}\right)\omega\right)$ Note the definition of $\epsilon_0$ means
$${}^{{}^{\omega}}{}^{\omega}\omega={}^{\epsilon_0}\omega$$
However, because exponentiation already lost associativity, tetration have almost no useful general identities except for finite integers $a,b,n$ (Using exponential identity $\alpha^{\beta\gamma}=(\alpha^\beta)^\gamma$ and that finite integers at the index commute)
$$(a^b)^{{}^{n-1}a}=({}^na)^b$$
Therefore, no known ordinal arithmetic can be used to simplify ${}^{\epsilon_0}\omega$. This prompt me to try a sandwiching approach as follows (where $f$ is some increasing normal function that serves as the iterative operation that climb up this sequence):
$${}^{({}^j\omega)}\omega<f({}^{({}^j\omega)}\omega)<{}^{({}^{j+1}\omega)}\omega$$
If we move through all $j < \omega$ and take the supremum, then we discover that ${}^{\epsilon_0}\omega$ is a fixed point of $f$. But nothing can be deduced further as I have no idea how to derive $f$ explicitly other than it is an operation that is similar to exponentiation, but only applied to the height i.e. $f: {}^j\omega \to \omega^{{}^j\omega}$. Either way, it had said nothing on whether ${}^{\epsilon_0}\omega$ is one of the epsilon number or where it is ordered wrt the usual ordinals.
Wikpedia also has a talk section saying that ordinal tetration is trivial, or that it will collapse into $\epsilon_0$. But as shown with ${}^{\epsilon_0}\omega$, I have so far failed to reproduce their results.
What other squeezing schemes (i.e. where should I put my $j$) I can use (or other properties of ordinals) in order to determine the value of ${}^{\epsilon_0}\omega$ ?
(NB You know my main focus is not large countable numbers since Veblen functions will eventually overshoot any hyperoperation)
So while in the main chat room, I devised a way to define tetration, and it works quite nicely.
$$\alpha\uparrow\uparrow\beta=\begin{cases}0,&\beta=-1\\1,&\beta=0\\\alpha,&\beta=1\\\alpha^{\alpha\uparrow\uparrow\zeta},&\alpha<\omega,\beta=\zeta+1>1\\(\alpha\uparrow\uparrow\zeta)^{\alpha\uparrow\uparrow\zeta},&\alpha\ge\omega,\beta=\zeta+1>1\\\sup\{\alpha\uparrow\uparrow(\beta[\eta])|\eta<\operatorname{cf}(\beta)\},&\beta\in\mathbb{Lim}\end{cases}$$
For finite $\alpha,\beta$, $\alpha\uparrow\uparrow\beta$ is what you expect it to be. Then,
$$\alpha\uparrow\uparrow\beta=\omega\forall\alpha<\omega\land\beta\ge\omega$$
Then we have,
$$\omega\uparrow\uparrow1=\omega\\\omega\uparrow\uparrow2=\omega^\omega\\\omega\uparrow\uparrow3=(\omega^\omega)^{\omega^\omega}=\omega^{\omega^\omega}$$
And so forth. Then,
$$\omega\uparrow\uparrow\omega=\varepsilon_0\\\omega\uparrow\uparrow(\omega+1)=\varepsilon_0^{\varepsilon_0}\\\vdots\\\omega\uparrow\uparrow(\omega2)=\varepsilon_1\\\vdots\\\omega\uparrow\uparrow(\omega(1+\beta))=\varepsilon_\beta\forall\beta\le\zeta_0$$
Thus, by this definition,
$$\omega\uparrow\uparrow\varepsilon_0=\varepsilon_{\varepsilon_0}$$
Edit:
The following definition may be nicer to use:
$$\alpha\uparrow^\beta\delta=\begin{cases}\alpha,&\delta=1\\\alpha^\delta,&\beta=1\\\sup\{(\alpha\uparrow^\beta\psi)\uparrow^\gamma(\alpha\uparrow^\beta\psi)|0<\gamma<\beta,0<\psi<\delta\},&\text{else}\end{cases}$$
And likewise extends to higher operations.