I'll have a try on #4, although I admit that my solution hides AC in a general result, which I won't prove here. And my first idea hasn't been that far from what I have (not exactly a derivation, but something linear instead). The whole problem is to connect the algebraic property of linearity with the topological property of the ##\lim \inf## and ##\lim \sup##.
I'll start with a description of
Hewitt, Stromberg, Real and Abstract Analysis, Springer in which the reals are defined as the cosets of a Cauchy sequence ##\mathfrak{C}## modulo null sequences ##\mathfrak{N}##. Both are subsets of all bounded real sequences ##\mathfrak{B}##.
With ##\mathfrak{N} \subset \mathfrak{C} \subset \mathfrak{B}## we get ##\mathbb{R} \cong \mathfrak{C} / \mathfrak{N} \subset \mathfrak{B} / \mathfrak{N} =: \mathfrak{A}## a real (topological, infinite dimensional, commutative) ##\mathbb{R}-##algebra ##\mathfrak{A}## (with ##1##).
The linear function ##L## can now be viewed as a linear functional on ##\mathfrak{A}##.
I will use the sub-additivity ##\overline{\lim}((x_n)+(y_n)) \leq \overline{\lim} (x_n)+\overline{\lim} (y_n)## and ##\underline{\lim}(x_n) = - \overline{\lim}(-x_n)##
Now here's the trick and the point where the real work is done:
If we define ##L(x_n)\vert_{\mathfrak{C} / \mathfrak{N}} := \lim(x_n)## then ##L(x_n) = \lim(x_n) \leq \overline{\lim}(x_n)## holds on ##\mathfrak{C} / \mathfrak{N} \subset \mathfrak{A} \;##, is linear and we may apply the theorem of Hahn-Banach. (Its proof uses Zorn's Lemma and therefore AC.)
This theorem guarantees us the existence of a linear functional ##L## on ##\mathfrak{A}## which extends our ##L## defined on ##\mathfrak{C} / \mathfrak{N} \; ##.
Linearity is already given by Hahn-Banach.
With that we have ##L(x_n) \leq \overline{\lim}(x_n)## by the theorem and
$$L(x_n)=-L(-x_n) \geq -\overline{\lim}(-x_n) = \underline{\lim}(x_n)$$
For the cutting process in point 5) we define a sequence ##z_n## by ##z_1=x_1## and ##z_n=x_{n} - y_{n-1}=0## for ##n > 1##.
Then ##L(z_n)=L(x_1,0, \dots) = 0 = L(x_n) - L(0,y_1,y_2, \dots) = L(x_n) - L(y_n)## and 5) is equivalent to ##L(0,y_1,y_2, \dots) = L(y_1,y_2, \dots)##. This is true for elements of ##\mathfrak{C} / \mathfrak{N}## which are dense in ##\mathfrak{A}## according to ##\Vert \, . \Vert_{\infty}##.
For ##0 \leq x_n## we have ##0 \leq \underline{\lim} (x_n) \leq L(x_n)##; and for ##\lim (x_n) = x## we have ##x = \underline{\lim} (x_n) \leq L(x_n) \leq \overline{\lim} (x_n) = x## and get equality everywhere. Therefore 4) and 6) hold.