\section{\bf Introduction}\label{sec:intro}

In this work, we are interested in studying the ramification behaviour of the $p$-torsion points of the formal group associated to an abelian variety over an unramified local field.

Let $p$ be a rational prime, let $F$ denote a finite, unramified extension of ${\mathbb Q}_p$, let $K$ be the completion of the maximal unramified extension of ${\mathbb Q}_p$, let ${\overline{K}}$ be some fixed algebraic closure of $K$, and let ${\cal O} \coloneqq \mathcal{O}_{{\overline{K}}}$.
Let $A$ be an abelian variety defined over $F$, with good reduction, let $\mathcal{A}$ denote the N\'eron model of $A$ over ${\rm Spec}({\cal O}_F)$, and let $\widehat{\mathcal{A}}$ be the formal completion of $\mathcal{A}$ along the identity of its special fiber, i.e. the formal group of $A$.

To state these results,  we need the following definitions.
In \cite{fontaine:differentials}, Fontaine studied the ${\cal O}$-module $\Omega:=\Omega^1_{{\cal O}/{\cal O}_K}\cong \Omega^1_{{\cal O}/{\cal O}_F}$ of K\"ahler differentials of ${\cal O}$ over ${\cal O}_K$, or over ${\cal O}_F$.
The ${\cal O}$-module $\Omega$ is a torsion and $p$-divisible ${\cal O}$-module, with a semi-linear action of $G_F$.
Let $d\colon{\cal O}\longrightarrow \Omega$ denote the canonical derivation, which is surjective.
We denote by ${\cal O}^{(1)}:=\ker(d)$, the kernel of $d$, which is an ${\cal O}_K$-sub-algebra of ${\cal O}$.
Membership of an element $x$ of ${\cal O}$ inside of ${\cal O}^{(1)}$ reflects ramification properties of $x$ (see e.g., \autoref{defn:delta} and \autoref{lemma:propertiesofdelta}).

By using previous work of the authors \cite{IMZ:Fontaine}, we are able to prove that the $\mathcal{O}^{(1)}$-points of the Tate module of $\mathcal{A}$ are trivial, which implies that following theorem.

 \begin{thmx}\label{xthm:main1}
Let $A$ be an abelian variety over $F$ with good reduction. Then there is $n_0\ge 1$ such that
for every $m\ge n_0$ and $0\neq P\in \widehat{\mathcal{A}}[p^m]({\cal O})\setminus \widehat{\mathcal{A}}[p^{n_0 -1 }]({\cal O})$, we have $P\notin\widehat{\mathcal{A}}({\cal O}^{(1)})$.
\end{thmx}

For a more concretely description of what this means in terms of the coordinates of the torsion point $P$, we refer the reader to  \autoref{remark:ram}.

Our second result gives conditions on $\widehat{\mathcal{A}}$ for which one may take $n_0 = 1$ in \autoref{xthm:main1}.
More precisely, we describe a condition on a formal group ${\mathscr F}$ of dimension $g$ over $\Spf({\cal O}_K)$ which implies that $0\neq P = (x_1,\dots,x_g) \in {\mathscr F}[p]({\cal O})$, the field of definition $K(P)/K$ is tamely ramified.
The condition is discussed in Section \ref{sec:ramification} and is related to a symmetric formal group law from \cite{arias}.

\begin{thmx}\label{xthm:tamelyramifiedtorsion}
Let ${\mathscr F}$ be a strict (\autoref{defn:strict}) formal group of dimension $g$ over $\Spf(\mathcal{O}_K)$.
For $0\neq P = (x_1,\dots,x_g) \in {\mathscr F}[p]({\cal O})$, the field of definition $K(P)/K$ is tamely ramified and ${\cal O}_{K(P)} \cong {\cal O}_K[x_1,\dots,x_g]$.
Moreover, $K({\mathscr F}[p])/K$ is tamely ramified.
\end{thmx}

\subsection{Related results}
In \cite[Section 1]{Serre:PropertiesGaloisEllipticCurves}, Serre showed that for $E/\mathbb{Q}_p$ an elliptic curve with good supersingular reduction, the field extension $\mathbb{Q}_p(E[p])/\mathbb{Q}_p$ is tamely ramified, and his proof relies on a detailed study of the formal group attached to $E$. In particular, Serre explicitly determined the $p$-adic valuation of the points on $E[p]$ (note that when $E/\mathbb{Q}_p$ has good supersingular reduction we have that $E[p]\cong \widehat{{\mathcal E}}[p]$ where $\widehat{{\mathcal E}}$ is the formal group of the N\'eron model of $E$), which allowed him to embed $E[p]$ into a certain vector space on which the wild inertia group acts trivially.

In her thesis, Arias-de-Reyna generalized this approach, and in \cite[Theorem 3.3]{arias}, she showed that if there exists a positive rational number such that for all $0\neq P = (x_1,\dots,x_g) \in \widehat{\mathcal{A}}[p]({\cal O})$, the minimum of $p$-adic valuation of the coordinates of $P$ equals $\alpha$, then the action of wild inertia on $\widehat{\mathcal{A}}[p]$ is trivial, and so $K(\widehat{\mathcal{A}}[p])/K$ is tamely ramified.
She goes on to define the notion of a symmetric formal group law on a formal group of dimension 2, and then proves \cite[Theorem 4.15]{arias} that if formal group of dimension 2 has a symmetric formal group law and height 4, then the above statement about the $p$-adic valuation of the $p$-torsion points holds.
Later in \cite[Theorem 5.9]{arias}, she identifies a family of genus 2 curves whose Jacobians have associated formal groups with a symmetric formal group law and height 4, and hence their $p$-torsion defines a tamely ramified extension.  We also mention work of Rosen and Zimmerman \cite{RosenZimmerman:TorsionFormal}, in which the authors study the Galois group of $K({\mathscr F}[p^n])$ where ${\mathscr F}$ is a generic commutative formal group of dimension 1 and height $h$.

In the global setting i.e., when working over a number field $F/\mathbb{Q}$, Coleman \cite{coleman1987ramified} studied the ramification properties of torsion points on abelian varieties in relation to the Manin--Mumford conjecture.
More precisely, he conjectured (\textit{loc.~cit.~}Conjecture B) that for a smooth, projective, geometrically integral curve $C/F$ and any Galois stable torsion packet $T$ in $C(\overline{F})$, the field $F(T)/F$ is unramified at a certain prime $\mathfrak{p}$ of $F$.
Coleman proved this conjecture when $\mathfrak{p}$ is large enough, and using work of Bogomolov, he provided a new proof of the Manin--Mumford conjecture.

We conclude by noting that our \autoref{xthm:tamelyramifiedtorsion} generalizes \cite[Theorem 4.15]{arias} and theoretically provides examples of abelian varieties of arbitrary dimension for which the extension $K(\widehat{\mathcal{A}}[p])/K$ is tamely ramified.

\subsection{Outline of paper}
In Section \ref{sec:Fontaine}, we recall the definition of the Fontaine integral, our previous work on the kernel of the Fontaine integral \cite{IMZ:Fontaine}, and a different perspective on the Fontaine integral via work of Wintenberger.
In Section \ref{sec:differentperspective}, we prove \autoref{xthm:main1}.
We conclude in Section \ref{sec:ramification} with the definition of a strict formal group and our proof of \autoref{xthm:tamelyramifiedtorsion}.

 \subsection{Conventions}\label{subsec:conventions}
We establish the following notations and conventions throughout the paper.

\subsubsection*{Fields}
Fix a rational prime $p>2$.
Let $K$ denote the completion of maximal unramified extension of ${\mathbb Q}_p$, let $\overline{K}$ be a fixed algebraic closure of $K$, and let ${\mathbb C}_p$ denote the completion of $\overline{K}$ with respect to the unique extension $v$ of the $p$-adic valuation on ${\mathbb Q}_p$ (normalized such that $v(p) = 1$).
For a tower of field extensions ${\mathbb Q}_p\subset F\subset K$, we denote by $G_K$ and respectively $G_F$ the absolute Galois groups of $K$ and $F$ respectively.
We denote ${\cal O}:={\cal O}_{{\overline{K}}}$.

\subsubsection*{Abelian varieties}
We will consider an abelian variety $A$ defined over some  subfield $F\subset K$ such that $[F:{\mathbb Q}_p]<\infty$, with good reduction over
$F$.
Let $\mathcal{A}$ denote the N\'eron model of $A$ over ${\rm Spec}({\cal O}_F)$ and also denote by $\widehat{\mathcal{A}}$ the formal completion of $\mathcal{A}$ along the identity of its special fiber, i.e. the formal group of $A$.
We note that the formation of N\'eron models commutes with unramified base change.
We will denote the Tate module of $A$ (resp.~the N\'eron model $\mathcal{A}$ of $A$) by $T_p(A)$ (resp.~$T_p(\mathcal{A})$).
We note that $T_p(A)\cong T_p(\mathcal{A})$ as $G_F$-modules.

\subsubsection*{Formal groups}
We will let ${\mathscr F}$ denote a formal group over $\Spf({\cal O}_K)$.
Recall that $\widehat{\mathcal{A}}$ is a formal group of dimension $\dim (A)$ and of height $h$ which satisfies $\dim(A) \leq h \leq 2\dim(A)$.
We refer the reader to \cite{Hazewinkel:FormalGroups} for an extensive treatment of formal groups and to \cite[Chapter 4.2 and Chapter 5]{AriasDeReyna:Thesis} for a more concise treatment.

 \section{\bf Fontaine integration for abelian varieties with good reduction}\label{sec:Fontaine}

In this section, we recall the construction of the Fontaine integration as well as our previous work concerning the kernel of the Fontaine integral.

\subsection*{The differentials of the algebraic integers}
First, we recall for the reader's convenience the notation established above.
Let $K$ denote the maximal unramified extension of ${\mathbb Q}_p$, let $\overline{K}$ be a fixed algebraic closure of $K$, and let ${\mathbb C}_p$ denote the completion of $\overline{K}$.
Let $G_K$ denote the absolute Galois group of $K$.
We denote ${\cal O}:={\cal O}_{{\overline{K}}}$. Fix a finite extension $F$ of ${\mathbb Q}_p$ in $K$.
For a $G_K$-representation $V$, the $n$-th Tate twist of $V$ is denoted by $V(n)$, which is just the tensor product of $V$ with the $n$-fold product of the $p$-adic cylcotomic character $\mathbb{Q}_p(1)$.

In \cite{fontaine:differentials}, Fontaine studied a fundamental object related to these choices, namely the ${\cal O}$-module $\Omega:=\Omega^1_{{\cal O}/{\cal O}_K}\cong \Omega^1_{{\cal O}/{\cal O}_F}$ of K\"ahler differentials of ${\cal O}$ over ${\cal O}_K$, or over ${\cal O}_F$.
The ${\cal O}$-module $\Omega$ is a torsion and $p$-divisible ${\cal O}$-module, with a semi-linear action of $G_F$.
Let $d\colon{\cal O}\longrightarrow \Omega$ denote the canonical derivation, which is surjective.

Important examples of algebraic differentials arise as follows:~Let $(\varepsilon_n)$ denote a compatible sequence of primitive $p$th roots of unity in $\overline{K}$.
Then
\[
\frac{d\varepsilon_n}{\varepsilon_n} = d (\log \varepsilon_n) \in \Omega \quad \mbox{ and } \quad  p\left(\frac{d\varepsilon_{n+1}}{\varepsilon_{n+I}}\right)=\frac{d\varepsilon_n}{\varepsilon_n}.
\]

Next, we recall a theorem of Fontaine.

\begin{theorem}[\protect{\cite[Th\'eor\`eme 1']{fontaine:differentials}}]\label{thm:Fontainedifferentials}
Let $(\varepsilon_n)$ denote a compatible sequence of primitive $p$th roots of unity in $\overline{K}$.
For $\alpha \in K$, write $\alpha = a/p^r$ for some $a\in \mathcal{O}$.
The morphism $\xi\colon \overline{K}(1) \longrightarrow \Omega$ defined by
\[
\xi(\alpha \otimes (\varepsilon_n)_n) = a \frac{d\varepsilon_r}{\varepsilon_r}
\]
is surjective and $G_K$-equivariant with kernel
\[
\ker(\xi) = \underline{a}_{K} \coloneqq \brk{x\in \overline{K} : v(x) \geq - \frac{1}{p-1}}.
\]
Moreover, $\Omega\cong{\overline{K}(1)}/{\underline{a}_K(1)} \cong (\overline{K}/\underline{a}_K)(1)$ and $V_p(\Omega) = \Hom_{\mathbb{Z}_p}(\mathbb{Q}_p,\Omega) \cong \mathbb{C}_p(1)$.
\end{theorem}

\autoref{thm:Fontainedifferentials} implies the following:
\[
T_p(\Omega)\otimes_{{\mathbb Z}_p}{\mathbb Q}_p:=\left(\varprojlim_n \Omega[p^n]\right) \otimes_{{\mathbb Z}_p}{\mathbb Q}_p\cong \left(\varprojlim\left(\Omega\stackrel{p}{\leftarrow}\Omega\stackrel{p}{\leftarrow}\cdots\stackrel{p}{\leftarrow}\Omega\cdots\right)\right)\otimes_{{\mathbb Z}_p}{\mathbb Q}_p\cong{\mathbb C}_p(1)
\]
as $G_F$-modules.

We denote by ${\cal O}^{(1)}:=\ker(d)$, the kernel of $d$, which is an ${\cal O}_K$-sub-algebra of ${\cal O}$.
Indeed, if $a,b\in {\cal O}^{(1)}$, then $d(ab)=ad(b)+bd(a)=0$, and so $ab\in {\cal O}^{(1)}$.
In order to better understand ${\cal O}^{(1)}$, we recall a construction from the first and last author \cite{iovita_zaharescu}.

\begin{definition}\label{defn:delta}
Let $a\in {\cal O}$. Let $L/K$ be a finite extension which contains $a$, let $\pi$ be a uniformizer of $L$, and let $f\in {\cal O}_K[x]$ be such that $a = f(\pi)$. Then, define
\[
\delta(a) \coloneqq \min \left( v\left( \frac{f'(\pi)}{\mathcal{D}_{L/K}}\right),0 \right)
\]
where $\mathcal{D}_{L/K}$ denotes the different ideal of $L/K$.
Note that $\delta$ does not depend on $\pi$, $f$, or $F$, and so it defines a function $\delta\colon {\cal O} \longrightarrow (-\infty,0]$.
\end{definition}

\begin{lemma}[Properties of $\delta$]\label{lemma:propertiesofdelta}
The function $\delta$ from \autoref{defn:delta} satisfies the following properties.
\begin{enumerate}
\item If $a,b\in {\cal O}$, then $\delta(a + b) \geq \min(\delta(a),\delta(b))$, and if $\delta(a) \neq \delta(b)$, then we have equality.
\item If $a,b\in {\cal O}$, then $\delta(ab) \geq \min(\delta(a) + v(b),\delta(b) + v(a))$.
\item If $f\in {\cal O}_K[x]$ and $\alpha \in {\cal O}$, then $\delta(f(\alpha)) = \min(v(f'(\theta)) + \delta(\theta) ,0)$.
\item If $x,y \in $, then $xdy = 0$ if and only if $v(x) + \delta(y) \geq 0$.
\item For $a\in {\cal O}$, $\delta(a) = 0$ if and only if $a\in {\cal O}^{(1)}$.
\item The formula $\delta(a db) \coloneqq \min(v(a) + \delta(b),0)$ is well-defined and give a map $\delta\colon \Omega \longrightarrow (-\infty,0]$, which makes the obvious diagram commutative.
\end{enumerate}
\end{lemma}

We will use the follow properties of $\delta$ in our study of the Fontaine integral.

\begin{lemma}[\protect{\cite[Lemma 2.2]{iovita_zaharescu}}]\label{lemma:relateddelta}
Let $a,b \in {\cal O}$ be such that $\delta(a)\leq \delta(b)$. Then there exists $c\in {\cal O}_{K[a,b]}$ such that $c d a = d b$.
\end{lemma}

\begin{proposition}[\protect{\cite[Theorem 2.2]{iovita_zaharescu}}]\label{prop:deeplyramified}
Let $L/K$ be an algebraic extension. Then $L$ is deeply ramified  (\textit{loc.~cit.~}Definition 1.1) if and only if $\delta({\cal O}_L)$ is unbounded.
\end{proposition}

\subsection*{The definition of Fontaine's integration}\label{sec:Fontaineintegral}
We are now ready to define Fontaine's integration.
Let $ H^0(\mathcal{A}, \Omega^1_{\mathcal{A}/{\cal O}_F})$ and respectively $\Lie(\mathcal{A})({\cal O}_F)$ denote the ${\cal O}_F$-modules of invariant differentials on $\mathcal{A}$ and respectively its Lie algebra.
Note that $\omega\in  H^0(\mathcal{A}, \Omega^1_{\mathcal{A}/{\cal O}_F})$ being invariant implies that $(x\oplus_{
\mathcal{A}} y)^*(\omega) = x^*(\omega) + y^*(\omega)$ and $[p]^*(\omega) = p\omega$ where $\oplus_{
\mathcal{A}}$ is the group law in $A({\overline{K}})$.

\begin{definition}\label{defn:Fontaineintegral}
Let $\underline{u}=(u_n)_{n\in {\mathbb N}}\in T_p(A)$ and $\omega\in H^0(\mathcal{A}, \Omega^1_{\mathcal{A}/{\cal O}_F})$.
Each $u_n\in \mathcal{A}({\cal O})$ corresponds to a morphism $u_n\colon {\rm Spec}({\cal O})\longrightarrow \mathcal{A}$, and hence we can pullback $\omega$ along this map giving us a K\"ahler differential $u_n^*(\omega)\in \Omega$.
The sequence $\left(u_n^*(\omega)\right)_{n\geq 0}$ is a sequence of differentials in $\Omega$ satisfying
$pu_{n+1}^\ast(\omega)=u_n^\ast(\omega)$, and hence defines an element in $V_p(\Omega)\cong {\mathbb C}_p(1)$.

The \cdef{Fontaine integration map}
\[
\varphi_\mathcal{A}\colon T_p(A)\longrightarrow \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F}{\mathbb C}_p(1)
\]
is a non-zero, $G_F$-equivariant map defined by
\[
\varphi_\mathcal{A}(\underline{u})(\omega):=\left(u_n^\ast(\omega)\right)_{n\ge 0}\in V_p(\Omega)\cong {\mathbb C}_p(1).
\]
\end{definition}

\begin{remark}\label{rem:differentdefn}
Using \autoref{thm:Fontainedifferentials} and the function $\delta$ from \autoref{defn:delta}, we can give an alternative description of the Fontaine integration map.
Let $\underline{u}=(u_n)_{n\geq 0}\in T_p(A)$ and $\omega\in H^0(\mathcal{A}, \Omega^1_{\mathcal{A}/{\cal O}_F})$.
Each $u_n\in \mathcal{A}({\cal O})$ corresponds to a morphism $u_n\colon {\rm Spec}({\cal O})\longrightarrow \mathcal{A}$, and hence we can pullback $\omega$ along this map giving us a K\"ahler differential $u_n^\ast(\omega)\in \Omega$.

For every $n \geq 0$, there is a maximal $m(n) \geq 0$ such that $u_n^*(\omega) = \alpha_n (d\varepsilon_{m(n)}/\varepsilon_{m(n)})$ with $\alpha_n \in {\cal O}$ where $\varepsilon_{m(n)}$ is some primitive $p^{m(n)}$-th root of unity.
To see this, we first note that
\[
\delta\left(\frac{d\varepsilon_r}{\varepsilon_r}\right) = -r - \frac{1}{p^r(p-1)}
\]
for any primitive $p^r$-th root of unity. This result follow from the definition of $\delta$ and a result of Tate \cite[Proposition 5]{tate:pdiv} on the valuation of the different ideal of $K(\varepsilon_r)/K$.
By taking $m(n) = -[\delta(u_n^\ast(\omega))]$ where $[x]$ denotes the greatest integer of the real number $x$, we can use  \autoref{lemma:propertiesofdelta}.(6) and \autoref{lemma:relateddelta} to deduce the above equality.

Now using \autoref{thm:Fontainedifferentials}, we have that
\[
\varphi_{\mathcal{A}}(\underline{u})(\omega) = \lim_{n\longrightarrow \infty}p^{n-m(n)}\alpha_n \in {\mathbb C}_p.
\]
Moreover, using the definition of $\delta$ and this above interpretation, we can see that if $\underline{u} \in T_p(A)^{G_K}$ (i.e., if $\underline{u}$ is an unramified path), then $\varphi_{\mathcal{A}}(\underline{u})(\omega) = 0$. Indeed, it is clear from the definition of $\delta$ that $m(n) = 0$.
\end{remark}

\subsection*{The kernel of the Fontaine integral}
In \cite{IMZ:Fontaine}, we studied the kernel of $\varphi_A$.
As noted in \autoref{rem:differentdefn},  we have that $T_p(A)^{G_K}$ lies in $\ker(\varphi_A)$, and in \cite[Theorem 4.5, Theorem A.4]{IMZ:Fontaine}, we showed that $T_p(A)^{G_K} = \ker(\varphi_A)$.
In proving these results, we determined the kernel of the Fontaine integral when restricted to the Tate module of the formal group of $A$. This result will play a role later on, and so we present it below.

\begin{theorem}[\protect{\cite[Theorem 5.5]{IMZ:Fontaine}}]\label{thm:kernelformal}
Let $A$ be an abelian variety over $F$ with good reduction, let $\mathcal{A}$ denote its N\'eron model, and let $\widehat{\mathcal{A}}$ be the formal group of $A$.
The Fontaine integral restricted to the Tate module of $\widehat{\mathcal{A}}$ is injective i.e., $\ker((\varphi_A)_{|T_p(\widehat{\mathcal{A}})}) = 0$.
\end{theorem}

 \subsection*{\bf Another point of view on the Fontaine integration map}\label{sec:another}

 In this subsection, we give another perspective on the Fontaine integration map, which will naturally lead us towards an application of \autoref{thm:kernelformal}.

We keep all the notations from the previous sections and Subsection \ref{subsec:conventions}.
Recall that we let ${\cal O}:={\cal O}_{{\overline{K}}}$ and we have the ${\cal O}$-module $\Omega:=\Omega^1_{{\cal O}/{\cal O}_K}$ with its canonical derivation $d\colon {\cal O}\longrightarrow\Omega$.
Note that $d$ is surjective and $\Omega$ is $p$-divisible, and let us denote ${\cal O}^{(1)}:=\ker(d)$.

\begin{lemma}\label{lemma:squarezero}
Let $A_{\rm inf}^{(1)}$ denote the $p$-adic completion of ${\cal O}^{(1)}$.
Then, the exact sequence of $G_F$-modules
\[
0\longrightarrow {\cal O}^{(1)}\longrightarrow {\cal O}\stackrel{d}{\longrightarrow}\Omega\longrightarrow 0
\]
induces another exact sequence:
\[
0\longrightarrow T_p(\Omega)\longrightarrow A_{\rm inf}^{(1)}\stackrel{\gamma}{\longrightarrow}{\cal O}_{{\mathbb C}_p}\longrightarrow 0,
\]
where $\gamma$ is an ${\cal O}_F$-algebra homomorphism and $T_p(\Omega)$ is seen as an ideal of
$A_{\rm inf}^{(1)}$ of square $0$.
\end{lemma}

\begin{proof}
The statement follows from \cite[Lemme 3.8]{Colmez:BdR} and also from \cite[Corollary 1.1]{iovita_zaharescu}, but we present another proof below.

We consider the diagram
\[
\begin{tikzcd}[row sep = 1.2em]
0 \arrow{r} & {\cal O}^{(1)} \arrow{r} \arrow{d}{p^n}&{\cal O} \arrow{r}{d}\arrow{d}{p^n} &\Omega \arrow{r}\arrow{d}{p^n} &0\\
0 \arrow{r} & {\cal O}^{(1)} \arrow{r} &{\cal O} \arrow{r}{d} &\Omega \arrow{r} &0.
\end{tikzcd}
\]
The snake lemma gives the exact sequence of $G_F$-modules:
\[
0\longrightarrow \Omega[p^n]\longrightarrow {\cal O}^{(1)}/p^n{\cal O}^{(1)}\longrightarrow {\cal O}/p^n{\cal O}\longrightarrow 0.
\]
By taking the projective limit with respect to $n$ of this exact sequence, we obtain the claim.
\end{proof}

Recall that we have the isomorphism $\Lie(\mathcal{A})({\cal O}_F)\cong H^0(\mathcal{A}, \Omega^1_{\mathcal{A}/{\cal O}_F})^\vee$.
By \autoref{lemma:squarezero}, we have the short exact sequence
\[
0\longrightarrow  T_p(\Omega) \longrightarrow A_{\rm inf}^{(1)}\longrightarrow {\cal O}_{{\mathbb C}_p}\longrightarrow 0,
\]
where $T_p(\Omega)$ is an ideal of $A_{\rm inf}^{(1)}$ such that $(T_p(\Omega))^2=0$.

By definition, we have
\[
\Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F}T_p(\Omega) \cong \ker\left(\mathcal{A}(A_{\rm inf}^{(1)})\longrightarrow \mathcal{A}({\cal O}_{{\mathbb C}_p})\right),
\]
and hence we have the following short exact sequence of abelian groups with $G_F$-action
\[
0\longrightarrow \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F}T_p(\Omega)\longrightarrow \mathcal{A}(A_{\rm inf}^{(1)})\longrightarrow \mathcal{A}({\cal O}_{{\mathbb C}_p})\longrightarrow 0.
\]
Consider the following commutative diagram with exact rows
\[
\begin{tikzcd}[row sep = 1.2em]
0 \arrow{r} & \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F}T_p(\Omega)\arrow{r} \arrow{d}{p^n}&\mathcal{A}(A_{\rm inf}^{(1)}) \arrow{r}{d}\arrow{d}{p^n} &\mathcal{A}({\cal O}_{{\mathbb C}_p}) \arrow{r}\arrow{d}{p^n} &0\\
0 \arrow{r} & \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F}T_p(\Omega)\arrow{r} &\mathcal{A}(A_{\rm inf}^{(1)}) \arrow{r}{d}&\mathcal{A}({\cal O}_{{\mathbb C}_p}) \arrow{r} &0.
\end{tikzcd}
\]
The snake lemma gives a $G_K$-equivariant map
\[
\nu_n\colon \mathcal{A}({\cal O}_{{\mathbb C}_p})[p^n]\cong A({\overline{K}})[p^n]\longrightarrow \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F} \Omega[p^n]
\]
and by taking the projective limit over $n$'s, we obtain a map
\[
\nu\colon T_p(A) \longrightarrow \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F}T_p(\Omega).
\]
\begin{proposition}\label{prop:anotherperspective}
The map obtained above
\[
\nu: T_p(A) \longrightarrow \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F}T_p(\Omega) \subset {\rm Lie}(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F}{\mathbb C}_p(1)
\]
coincides with Fontaine's integral, i.e.~we have $\nu=(\varphi_\mathcal{A})$.
\end{proposition}

\begin{proof}
In \cite[Section 4, page 394]{wintenberger:comparision}, Wintenberger used a generalization of the above construction to obtain an integration pairing which coincides with the Colmez integration pairing $\langle \cdot ,\cdot \  \rangle_{\Cz}$.
The result now follows from \cite[Proposition 6.1]{Colmez1992}.
\end{proof}

 \section{\bf Consequences of \autoref{thm:kernelformal}: ramification of $p$-power torsion points on $\widehat{\mathcal{A}}$}
 \label{sec:differentperspective}

In this section, we use the interpretation of the Fontaine integral from \autoref{prop:anotherperspective} and \autoref{thm:kernelformal} to deduce properties concerning the ramification of $p$-power torsion points on the formal group of $A$.

To begin, we recall the diagram
\[
\begin{tikzcd}[row sep = 1.2em]
0 \arrow{r} & \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F}T_p(\Omega)\arrow{r} \arrow{d}{p^n}&\mathcal{A}(A_{\rm inf}^{(1)}) \arrow{r}{d}\arrow{d}{p^n} &\mathcal{A}({\cal O}_{{\mathbb C}_p}) \arrow{r}\arrow{d}{p^n} &0\\
0 \arrow{r} & \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F}T_p(\Omega)\arrow{r} &\mathcal{A}(A_{\rm inf}^{(1)}) \arrow{r}{d}&\mathcal{A}({\cal O}_{{\mathbb C}_p}) \arrow{r} &0.
\end{tikzcd}
\]
Above, we only wrote a piece of the snake lemma, and by writing more of it, we have an exact sequence of $G_K$-modules
\[
0\longrightarrow \mathcal{A}(A_{\rm inf}^{(1)})[p^n]\longrightarrow \mathcal{A}({\cal O})[p^n]\longrightarrow \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F} \Omega[p^n].
\]
By taking projective limits, we have the \textit{exact} sequence
\begin{equation}\label{eqn:reducingconj}
0\longrightarrow T_p(\mathcal{A}(A_{\rm inf}^{(1)}))\longrightarrow T_p(A)\stackrel{\varphi_\mathcal{A}}{\longrightarrow}\Lie(\mathcal{A})({\cal O}_F)\otimes _{{\cal O}_F} T_p(\Omega)\subset \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F} {\mathbb C}_p(1).
\end{equation}
Therefore, \autoref{thm:kernelformal} implies that $T_p(\widehat{\mathcal{A}}(A_{\rm inf}^{(1)}))=0$.

To study consequences of this property, we will use another ring instead of $A_{\rm inf}^{(1)}$.

\begin{definition}[\protect{\cite{fontaine}}]\label{defn:Df}
Let $\theta\colon A_{\rm inf}^{(1)}\longrightarrow {\cal O}_{{\mathbb C}_p}$ denote the projection map. Then, we define $D_f:=\theta^{-1}({\cal O})$.
In \cite[Remark 1.4.7]{fontaine}, Fontaine gives the following construction of $D_f$.
Let us recall that
\[
V_p(\Omega)=T_p(\Omega)\otimes_{{\mathbb Z}_p} {\mathbb Q}_p=\varprojlim \left(\Omega\stackrel{p}{\leftarrow} \Omega\stackrel{p}{\leftarrow}\cdots\stackrel{p}{\leftarrow}\Omega \cdots  \right)
\]
and that $\Omega$ and $V_p(\Omega)$ are ${\cal O}$-modules. We make $R:=V_p(\Omega)\oplus {\cal O}$ into a commutative ring by defining multiplication as follows:~$(u, \alpha)(v, \beta)=(\beta u+\alpha v, \alpha\beta)$ for $(u, \alpha), (v, \beta)\in R$, i.e. we require that $V_p(\Omega)$ is an ideal of $R$ of square $0$. Then we have
\[
D_f=\{\left(u=(u_n)_{n\ge 0},\alpha\right)\in R\ |\ d(\alpha)=u_0\}.
\]
\end{definition}

By \autoref{defn:Df}, we have an exact sequence of $G_K$-modules
\[
0\longrightarrow T_p(\Omega)\longrightarrow D_f\stackrel{\theta}{\longrightarrow} {\cal O}\longrightarrow 0,
\]
where $\theta(u,\alpha)=\alpha$, and the $p$-adic completion of $D_f$ is $A_{\rm inf}^{(1)}$.
We note that we may construct the diagram above in the same way using $D_f$ instead of $A_{\rm inf}^{(1)}$, which produces the exact sequence \eqref{eqn:SESDf} with $D_f$ instead of $A_{\rm inf}^{(1)}$.
Instead of the exact sequence \eqref{eqn:reducingconj} above, we will have the following \textit{exact} sequence
\[
0\longrightarrow T_p(\mathcal{A}(D_f))\longrightarrow T_p(A)\stackrel{\varphi_\mathcal{A}}{\longrightarrow}\Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F} T_p(\Omega)\subset \Lie(\mathcal{A})({\cal O}_F)\otimes_{{\cal O}_F} {\mathbb C}_p(1).
\]
Again, the \autoref{thm:kernelformal} implies that $T_p\bigl(\widehat{\mathcal{A}}(D_f)\bigr)=0$.

We will use this observation to deduce that $T_p\bigl(\widehat{\mathcal{A}}({\cal O}^{(1)})\bigr)=0$.
In order to do so, we need to show that $\widehat{\mathcal{A}}[p^n](D_f)\cong \widehat{\mathcal{A}}[p^n]({\cal O}^{(1)})$, for all $n\ge 1$, which is accomplished through the following two lemmas.

\begin{lemma}\label{lemma:O^1torsionDftorsion}
Let $x\in \mathcal{A}[p^n]({\cal O}^{(1)})$, then there is $x'\in \mathcal{A}(D_f)$ with $\theta(x')=x$ and such that $[p^n](x')=0$.
\end{lemma}

\begin{proof}
Recall that $D_f:=\left\lbrace\bigl((x_n)_n, y)\in V_p(\Omega)\times {\cal O}\ |\ x_0=dy\right\rbrace$, i.e. we have an exact sequence
\[
0\longrightarrow T_p(\Omega)\longrightarrow D_f\stackrel{\theta}{\longrightarrow}{\cal O}\longrightarrow 0
\]
and this exact sequence splits over ${\cal O}^{(1)}\subset {\cal O}$, i.e. the following diagram is cartesian and has exact rows
\[
\begin{array}{cccccccccc}
0&\longrightarrow &T_p(\Omega)&\longrightarrow& D_f&\stackrel{\theta}{\longrightarrow}&{\cal O}&\longrightarrow& 0\\
&&||&&\cup&&\cup\\
0&\longrightarrow &T_p(\Omega)&\longrightarrow &T_p(\Omega)\oplus {\cal O}^{(1)}&\stackrel{\theta}{\longrightarrow}&{\cal O}^{(1)}&\longrightarrow& 0.
\end{array}
\]
In particular, the section $s\colon {\cal O}^{(1)} \longrightarrow D_f$ is defined by $s(x):=(0,x)$. Then $s$ defines a morphism $s\colon \mathcal{A}({\cal O}^{(1)})\longrightarrow \mathcal{A}(D_f)$, and if $x\in \mathcal{A}[p^n]({\cal O}^{(1)})$, then $s(x)\in \mathcal{A}[p^n](D_f)$.
\end{proof}

In the next lemma, we show a converse of \autoref{lemma:O^1torsionDftorsion} at least for the formal group $\widehat{\mathcal{A}}$ of $A$.

\begin{lemma}
\label{prop:dftorsion}
Let $\widehat{\mathcal{A}}$ denote the formal group of the abelian scheme $\mathcal{A}$ and fix $n\ge 1$ an integer.
Let
\[
0\neq P\in \widehat{\mathcal{A}}({\cal O})[p^n]\backslash \widehat{\mathcal{A}}({\cal O}^{(1)})
\]
and let $Q\in \widehat{\mathcal{A}}(D_f)$ be a point such that $\theta(Q)=P$. Then $[p^m]Q\neq 0$ for all $m\ge n$.
\end{lemma}

\begin{proof}
Let $m\ge n$ and denote by $$[p^m](X_1,\dots,X_g):=\left(f_1(X_1,\dots,X_g), f_2(X_1,\dots,X_g),\dots, f_g(X_1,\dots,X_g)\right)$$ the multiplication by $p^m$ on $\widehat{\mathcal{A}}$.
Let  $S:={\cal O}_F[[X_1,\dots,X_g)]]/I$, where $I$ is the ideal generated by $f_1, f_2,\dots,f_g$, then we know $S$ is a finite flat ${\cal O}_F$-algebra,  in particular $S$ is a free ${\cal O}_F$-module and $\widehat{\mathcal{A}}[p^m]:={\rm Spec}(S)$ with the co-multiplication of $\widehat{\mathcal{A}}$, is a finite flat group-scheme, so $\widehat{\mathcal{A}}\times_{{\cal O}_F}{\rm Spec}(F)$ is an \'etale, therefore smooth, group-scheme over $F$. This implies that the image in $S\otimes_{{\cal O}_F}F$ of the determinant of the matrix:
$$
\left( \begin{array}{cccccc} \frac{\partial(f_1)}{\partial(X_1)} & \frac{\partial(f_1)}{\partial(X_2)}&\dots&\frac{\partial(f_1)}{\partial(X_g)} \\
\vdots&\vdots&\ddots&\vdots\\
\frac{\partial(f_g)}{\partial(X_1)} & \frac{\partial(f_g)}{\partial(X_2)}&\dots&\frac{\partial(f_g)}{\partial(X_g)}\end{array}\right)
$$
is a unit.

Let now $P=(x_1,x_2,\dots,x_g)\in \mathfrak{m}_{{\cal O}}^g\in \widehat{\mathcal{A}}[p^n]({\cal O})\backslash \widehat{\mathcal{A}}[p^n]({\cal O}^{(1)})$, i.e. there is $1\le i\le g$ such that
$x_i$ is not in ${\cal O}^{(1)}$. Let $P'=(y_1,y_2,\dots,y_g)\in \widehat{\mathcal{A}}(D_f)$ such that $\theta(P')=P$, i.e. $y_j=\alpha_j+x_j$, with
$\alpha_j\in V_p(\Omega)$ and $d(x_j)=\alpha_{j,0}$, for all $1\le j\le g$. By the above assumption $\alpha_i\neq 0$.
As in $D_f$ we have $\alpha_j\alpha_k=0$ for all $1\le j,k\le g$, the Taylor formula implies that if $[p^m](P')=0$ we must have:
\[
  f_s(x_1,\dots,x_g)+ \sum_{j=1}^g\frac{\partial(f_s)}{\partial(X_j)}(x_1,\dots,x_g)\alpha_j=\sum_{j=1}^g\frac{\partial(f_s)}{\partial(X_j)}(x_1,\dots,x_g)\alpha_j =0
\]
for every $1\le s\le g$.
But the determinant of the matrix
\[
\left( \begin{array}{cccccc} \frac{\partial(f_1)}{\partial(X_1)}(x_1,\dots,x_g) & \frac{\partial(f_1)}{\partial(X_2)}(x_1,\dots,x_g)&\dots&\frac{\partial(f_1)}{\partial(X_g)}(x_1,\dots,x_g) \\
\vdots&\vdots&\ddots&\vdots\\
\frac{\partial(f_g)}{\partial(X_1)}(x_1,\dots,x_g) & \frac{\partial(f_g)}{\partial(X_2)}(x_1,\dots,x_g)&\dots&\frac{\partial(f_g)}{\partial(X_g)}(x_1,\dots,x_g)\end{array}\right)
\]
is a unit in $\overline{K}$, i.e. it is non-zero and $\alpha_i\neq 0$. This is a contradiction.
\end{proof}

\begin{remark}
We note that the group-scheme $\widehat{\mathcal{A}}[p^m]$ is not smooth over ${\cal O}_F$ (for example $S/pS$ could have nilpotents).
We also remark that as $\widehat{\mathcal{A}}[p^m]\times_{{\cal O}_F}{\rm Spec}(F)$ is smooth, the map $\theta\otimes 1\colon \widehat{\mathcal{A}}[p^m](D_f\otimes_{{\cal O}_F}K)\longrightarrow \widehat{\mathcal{A}}[p^m](\overline{F})$ is surjective, but this is clear as $D_f\otimes_{{\cal O}_F}F=V_p(\Omega)\oplus \overline{F}$.
\end{remark}

\autoref{lemma:O^1torsionDftorsion} and \autoref{prop:dftorsion} imply that the map $\theta$ gives an isomorphism $\widehat{\mathcal{A}}[p^n](D_f)\cong
\widehat{\mathcal{A}}[p^n]({\cal O}^{(1)})$, for all $n\ge 1$.
Combining this with the fact that $T_p\bigl(\widehat{\mathcal{A}}(D_f)\bigr)=0$, we have the following result.

\begin{theorem}[$=$\autoref{xthm:main1}]
\label{lemma:ram}
Let $A$ be an abelian variety over $F$ with good reduction. Then there is $n_0\ge 1$ such that
for every $m\ge n_0$ and $0\neq P\in \widehat{\mathcal{A}}[p^m]({\cal O})\setminus \widehat{\mathcal{A}}[p^{n_0 -1 }]({\cal O})$, we have $P\notin\widehat{\mathcal{A}}({\cal O}^{(1)})$
\end{theorem}

\begin{remark} \label{remark:ram}
We observe that the above is a result regarding ramification properties of the $p
$-power torsion points of the formal group of our abelian  variety with good reduction over $F$
(c.f.~\autoref{prop:deeplyramified}). More precisely, let $m\ge n_0$ and $0\neq P\in
\widehat{\mathcal{A}}[p^m]({\cal O})\setminus \widehat{\mathcal{A}}[p^{n_0 -1 }]({\cal O})$ be as in the theorem above.
Let $P=(x_1,x_2,...,x_g)$ with $x_i\in {\cal O}$ for $1\le i\le g$.  \autoref{lemma:ram} says the following: let
$L=K[P]:=K[x_1,x_2,...,x_g]$, let $\pi$ denote a uniformizer of $L$ and let $\mathcal{D}_{L/K}$ denote
the different ideal of $L/K$. For every $1\le i\le g$ let $f_i(X)\in {\cal O}_K[X]$
be polynomials such that $f_i(\pi)=x_i$, for every $i$. Then there is $1\le j\le g$ such that
$v(f_j'(\pi))< v(\mathcal{D}_{L/K})$ (i.e. $x_j\notin {\cal O}^{(1)}$).
\end{remark}

\section{\bf A theorem on the ramification type of the field obtained by adjoining a $p$-torsion point of a formal group}
\label{sec:ramification}
In this subsection we study the ramification properties of the extension $K[P]/K$, where $P$ is a non-zero $p$-power torsion point of the formal group of $A$.

We continue to denote by $K$ the completion of the maximal unramified extension of ${\mathbb Q}_p$ in an algebraic closure of ${\mathbb Q}_p$, which we denote ${\overline{K}}$.
Let ${\mathscr F}$ denote a formal group of dimension $g$ over ${\rm Spf}({\cal O}_K)$.
For example,  ${\mathscr F}$ can be the formal group of the N\'eron model of $A$ over ${\rm Spf}({\cal O}_K)$.

To begin, we define the notion of a strict formal group.

\begin{definition}\label{defn:strict}
Consider the multiplication-by-$p$ map \[[p](X_1,\dots.,X_g) = (f_1(X_1,\dots.,X_g),f_2(X_1,\dots.,X_g),\dots,f_g(X_1,\dots.,X_g))\] on ${\mathscr F}$ where each $f_i(X_1,\dots.,X_g)$ is a power series in with coefficients in ${\cal O}_K$.
For each $1\leq i\leq g$, define $F_i(X_1,\dots,X_g)$ to be the form comprised of monomials of $f_i$ which have unit coefficient and minimal degree, where we consider each monomial $X_1,\dots,X_g$ to be of degree 1.

Let $d_1,\dots,d_g$ denote the degree of these forms $F_1(X_1,\dots,X_g),\dots,F_g(X_1,\dots,X_g)$, respectively, which we note are (possibly distinct) powers of $p$.
Let $G_1(X_1,\dots,X_g),\dots,G_g(X_1,\dots,X_g)$ denote the reductions modulo $p$ of the forms $F_1(X_1,\dots,X_g),\dots,F_g(X_1,\dots,X_g)$.

Consider the system of equations
\begin{equation}\label{eqn:strict}
G_1(X_1,\dots,X_g) = G_2(X_1,\dots,X_g) = \cdots = G_g(X_1,\dots,X_g) = 0,
\end{equation}
We say that the formal group ${\mathscr F}$ is \cdef{strict} if $d_1 = d_2 = \cdots = d_g$ and the only solution to \eqref{eqn:strict} is $(0,0,\dots,0) \in (\overline{\mathbb{F}_p})^g$.
\end{definition}

\begin{remark}\label{rem:strict1}
If ${\mathscr F}$ is a formal group of dimension $1$, then it is clear that ${\mathscr F}$ is strict since we will have that $F_1(X_1) = uX_1^{p^h}$, where $h$ is the height of ${\mathscr F}$.
Moreover, if ${\mathscr F}$  is the product of $1$-dimensional formal groups, then again ${\mathscr F}$ is strict.
\end{remark}

\begin{remark}\label{rem:strict2}
We can given an equivalent characterization of strict as follows.
Consider the $g\times g$ matrix $M = (a_{ij})$ where the entry $a_{ij}$ consists of the coefficient of $X_i$  in the linear form $G_j$ for each $1\leq i,j \leq g$.
The condition that ${\mathscr F}$ be strict is equivalent to the determinant of $M$ being non-zero.

We refer the reader to \cite[Remark 4.14]{arias} for an example of a $2$-dimensional formal group of height $4$ where this condition holds, and here we note that the above degrees all equal $p^2$. Moreover, we remark that the proof from \textit{loc.~cit.~}holds for any $g$-dimensional formal group of height $2g$ where the degrees $d_1 = d_2 = \cdots = d_g$ all equal $p^2$.
\end{remark}

With this definition, we can state our main result.

\begin{theorem}\label{thm:tamelyramifiedtorsion}
Let ${\mathscr F}$ be a strict formal group of dimension $g$.
For $0\neq P = (x_1,\dots,x_g) \in {\mathscr F}[p]({\cal O})$, the field of definition $K(P)/K$ is tamely ramified and ${\cal O}_{K(P)} \cong {\cal O}_K[x_1,\dots,x_g]$.
Moreover, $K({\mathscr F}[p])/K$ is tamely ramified.
\end{theorem}

\begin{proposition}
\label{prop:torsionram}
Let ${\mathscr F}$ be a strict the formal group of dimension $g$.
For every $0\neq P\in{\mathscr F}[p]({\cal O})$, the coordinates of $P$ are not all in ${\cal O}^{(1)}$.
\end{proposition}

\begin{proof}
Let $0\neq P = (x_1,\dots,x_g) \in {\mathscr F}[p]({\cal O})$ be a non-zero $p$-torsion point.
By \autoref{thm:tamelyramifiedtorsion}, we know that the extension $K(P)/K$ is tamely ramified and that there exists some coordinate $x_i$ which is a uniformizer for $K(P)/K$.
By \cite[Proposition III.6.13]{Serre}, we have that $v(\mathcal{D}_{K(P)/K}) > 0 $ where $\mathcal{D}_{K(P)/K}$ is the different ideal of $K(P)/K$.
Now since $x_i$ is a uniformizer for $K(P)/K$, we have that
\[
\delta(x_i) = -v(\mathcal{D}_{K(P)/K})  < 0,
\]
where $\delta$ is the function defined in \autoref{defn:delta}.
By \autoref{lemma:propertiesofdelta}.(5), $x_i \notin {\cal O}^{(1)}$ as desired.
\end{proof}

For the remainder of this section, we focus on proving \autoref{thm:tamelyramifiedtorsion}. The proof can be broken down into three steps.
\begin{enumerate}
\item Given a non-zero $p$-torsion point $P = (x_1,\dots,x_g) \in {\mathscr F}[p]({\cal O})$, we will carefully construct linear combinations $z_i^*$ of the $x_1,\dots,x_n$ which satisfy nice properties in terms of their valuations and distances between their $\overline{K}$-conjugates. See \autoref{lemma:niceisom}.
\item Next, we consider the change of variables (i.e., the isomorphism of formal groups) which sends the coordinate $X_i$ to the linear combination $Z_i^*$ described above. We use the properties of the $z_i^*$ and the strictness of ${\mathscr F}$ to precisely determine the valuation of $z_i^*$ and to estimate the valuation of the difference between them and their $\overline{K}$-conjugates.
\item Finally, we use Krasner's lemma to deduce that one of the original coordinates $x_1,\dots,x_g$ of $P$ must be a uniformizer for the maximal order of $K(P)$, from which \autoref{thm:tamelyramifiedtorsion} follows.
\end{enumerate}

\begin{lemma}\label{lemma:niceisom}
Let ${\mathscr F}$ be a formal group of dimension $g$ over $\Spf({\cal O}_K)$.
Let $0\neq P= (x_1,\dots,x_g) \in {\mathscr F}[p]({\cal O})$.
There exist linear combinations $z_1^*,\dots,z_g^*$ of $x_1,\dots,x_g$ with coefficients in $({\cal O}_K)^{\times} \cup \{0\}$ which satisfy:
\begin{enumerate}
\item $K(z_{i}^*) \cong K(P)$,
\item $v(z_{i}^*) =  \min\{ v(x_{1}),\dots, v(x_{g})\}$,
\item $v(z_{i}^* - \sigma(z_{i}^*)) = \min\{ v(x_{1} - \sigma(x_{1})),\dots, v(x_{g} - \sigma(x_{g}))\}$ for all $\sigma \in \gal(\widetilde{K(P)}/K)$ where $\widetilde{K(P)}$ is the Galois closure of $K(P)$,
\end{enumerate}
and such the matrix $M$ representing the change of coordinates $(z_1^*,\dots,z_g^*)^t= M (x_1,\dots,x_g)^t$ is invertible. Here
the exponent $t$ indicates the transpose of a matrix.
\end{lemma}

\begin{proof}
Let $e := [K(P) : K]$.
Our proof will involve making a series of linear combinations.
To begin, we will construct the element $z_1^*$.
First, consider all the linear combinations of the form
\begin{equation}\label{eqn:1s}
\mathcal{B}_1:=\{z = u_1x_1 + \cdots + u_gx_g \text{ where } u_i\in({\cal O}_K)^{\times} \cup \{0\} \text{ and }u_1 \neq 0 \}.
\end{equation}

By our assumptions on $K$, the set of $u_i (\mbox{mod }p)$ is infinite, with $u_i$ as in the above formula, and hence we may find one linear combination, call it $z_1$, satisfying the following two conditions:
\begin{enumerate}
\item $v(z) \geq v(z_1)$ for all other linear combinations $z$ from $\mathcal{B}_1$,
\item $K(z_1) \cong K(P)$.
\end{enumerate}
To show that condition (2) holds, consider the following.
There are exactly $e$ embeddings of $K(x_1, \dots ,x_g)$ into the fixed algebraic closure ${\overline{K}}$, call them $\sigma_1, \dots, \sigma_e$. Note that the vectors ${\textbf{w}}_j  := (\sigma_j(x_1), \dots, \sigma_j(x_g))$, $ 1\le j \le e$, are distinct. Indeed, if for some $ i \ne j$ the vectors ${\textbf{w}}_i$ and ${\textbf{w}}_j$ coincide, then $\sigma_i$ and $\sigma_j$ will coincide at $x_1,\dots,x_g$ and so they will coincide on $K(x_1,\dots,x_g)$, which is not the case.

Consider now, for each pair $(i,j)$ with $1\le i, j \le e$ and $i\ne j$, the hyperplane ${\mathcal H}_{i,j}$ given by
\[
{\mathcal H}_{i,j} = \brk{ (c_1,c_2,\dots,c_g) :  c_1, \dots, c_g \in K , \sum_{l = 2}^g  c_l ( \sigma_i(x_l) - \sigma_j(x_l)) = 0 }.
\]
Since the vectors ${\textbf{w}}_j$, $1 \le j \le e$ are distinct, none of ${\mathcal H}_{i,j}$ covers the full space.
Denote by ${\mathcal H}$ the union of these finitely many hyperplanes.
Choose now any
$c_1,\dots,c_g \in K$ such that the point $(c_1,c_2,\dots,c_g)$ lies outside ${\mathcal H}$. Then we claim that the element
\[
z := c_1x_1+ \cdots + c_g x_g
\]
satisfies $K(z) = K(x_1,\dots,x_g)$. Indeed, $\sigma_1(z),\dots,\sigma_e(z)$ are distinct. For, if two of them are equal, say
$\sigma_i(z) = \sigma_j(z)$ with $i\ne j$, then $(c_1,c_2\dots,c_g) $ is forced to lie in ${\mathcal H}_{i,j}$. Thus
$\sigma_1(z),\dots,\sigma_e(z)$ are distinct, so $z$ has at least $e$ distinct conjugates over $K$. Hence $[K(z):K] \ge e$ and in conclusion $K(z) = K(x_1,\dots,x_g)$ i.e., $z_1$ has degree $e$ over $K$.

We pause to note that the matrix $M$ representing the change of coordinates $(z_1,x_2,\dots,x_g)^t= M (x_1,x_2,\dots,x_g)^t$ is invertible.
Indeed, the matrix $M$ has units along the diagonal, the coefficients of the linear combination $z_1$ in the first row, and zeros elsewhere, hence the determinant is a unit.

We now look at the distances between these linear combinations and various of their conjugates over $K$.
Fix $\sigma \in \text{Gal}(\overline{{\mathbb Q}_p}/K)$ and consider the infimum of the values $v(z - \sigma(z))$ for the linear combinations $z$ as in \eqref{eqn:1s}, i.e. we look at $\inf\{ v\left(z-\sigma(z)\right)\ |\ z\in \mathcal{B}_1\}$.
We first show that this infimum exists and is attained by some linear combination.
Note that if $\sigma_{|K(z_1)} = \text{id}$, then since all linear combinations belong to $K(x_1,\dots,x_g)$, we have that $z - \sigma(z) = 0$ for all linear combinations $z$, i.e. $\inf\{ v\left(z-\sigma(z)\right)\ |\ z\in \mathcal{B}_1\}=\infty=v(z_1-\sigma(z_1))$.
If we consider $\sigma$ such that $\sigma_{|K(z_1)} \neq \text{id}$, then we at least have that $v(z - \sigma(z)) < \infty$ for some linear combinations $z$, for example for $z = z_1$.
We have that the set of possible valuations is discrete, and the set of values $v(z - \sigma(z))$ has a lower bound, namely $\min\{v(x_1 - \sigma(x_1)),\dots,v(x_n - \sigma(x_g))\}$.
Thus, the infimum of the set of values $v(z - \sigma(z))$ is attained by some linear combination $z$ from \eqref{eqn:1s}, but note it need not be obtained by $z_1$.

Let $G$ denote $\gal(\widetilde{K(P)}/K)$ where $\widetilde{K(P)}$ is the Galois closure of $K(P)/K$.
For a fixed  $\sigma\in G$, let $z_{\sigma}\in \mathcal{B}_1$ denote one such linear combination attaining the minimum $v(z_\sigma - \sigma(z_\sigma))=\min\{ v\left(z-\sigma(z)\right)\ |\ z\in \mathcal{B}_1\}$ .
We claim that we can find a linear combination of $z_1$ and of all of these $z_{\sigma}$ where $\sigma \in G$ which simultaneously achieves these minima.
To do this, consider all linear combinations
\begin{equation}\label{eqn:2}
z^* = z_1 + \sum_{\sigma \in G}u_{\sigma}z_{\sigma} \quad  \text{ where } u_{\sigma} \in ({\cal O}_K)^{\times}\cup \{0\}.
\end{equation}
We will choose the $u_\sigma$'s such that each such $z^*$ will live in $\mathcal{B}_1$.

To achieve our desired simultaneous minima, we start with one $\sigma \in G$, call it $\sigma_1$.
First, we let $z^*= z_1$.
This linear combination might work in that it already attains the minimum at $\sigma_1$; by this we mean that $v(z - \sigma_1(z)) \ge v(y_1 - \sigma_1(y_1))$ holds for all linear combinations $z\in \mathcal{B}_1$ from \eqref{eqn:1s}.
If this is the case, then we set $y_1:=z_1$.
Now suppose that $z_1$ does not attain the minimum at $\sigma$.
In this case, we may use any unit $u\in ({\cal O}_K)^{\times}$ and set $z^*:=  z_1 + uz_{\sigma_1}$.
Indeed, for any unit $u \in ({\cal O}_K)^{\times}$ and $z^*$ above, we have that
$$v(z^* - \sigma_1(z^*)) =v\left(z_1-\sigma_1(z_1)+u(z_{\sigma_1}-\sigma_1(z_{\sigma_1})\right)=
 v(z_{\sigma_1} - \sigma_1(z_{\sigma_1})),
 $$
 because $v(z_1-\sigma_1(z_1))>v(z_{\sigma_1}-\sigma_1(z_{\sigma_1})$.
Let $y_1:=z_1+uz_{\sigma_1}$ with unit $u \in ({\cal O}_K)^{\times}$ such that $y_1\in \mathcal{B}_1$. We have that for such
 $y_1\in \mathcal{B}_1$, $v(y_1-\sigma_1(y_1))\le v(z-\sigma_1(z)$ for all $z\in \mathcal{B}_1$ and that the $u's$ with $y_1\in \mathcal{B}_1$
 have the property that $u \pmod p$ avoids a finite number of elements in $\overline{{\mathbb F}_p}$.

For another automorphism $\sigma_2\in G$,  we proceed along the same lines, that is: if $y_1$ has the property
that $v(y_1-\sigma_2(y_1))\le v(z-\sigma_2(y_1))$ for all $z\in \mathcal{B}_1$ we set $y_2:=y_1$.
If the above is not true, let $z^*:=y_1+ux_{\sigma_2}$, for some $u\in {\cal O}_K^\times$. Then as above we have:
$v(z^*-\sigma_2(z^*))=v\left(x_{\sigma_2}-\sigma_2(x_{\sigma_2})\right)$ therefore $z^*$ realizes the minimum for $\sigma_2$, for all
$u's$ for which $z^*\in \mathcal{B}_1$. What about for $\sigma_1$. The worst that can happen is that $v(y_1-\sigma_1(y_1))=v\left(x_{\sigma_2}-\sigma_1(x_{\sigma_2})\right)$, i.e. ~if we denote by $\pi$ a uniformizer of $K(P)$, we have $y_1-\sigma_1(y_1)=a\pi^\alpha$
$x_{\sigma_2}-\sigma_1(x_{\sigma_2})=b\pi^\alpha$, with $a,b\in {\cal O}_{K(P)}^\times$. Therefore $z^*-\sigma_1(z^*)=(a+ub)\pi^a.$
Now the residue field of ${\cal O}_{K(P)}$ is $k$, therefore by choosing $u\in {\cal O}_K^{\times}$ such that
$a+ub(\mbox{ mod }\pi)\neq 0$ we have $v(z^*-\sigma_1(z^*))=v(y_1-\sigma_1(y_1))$ and therefore $y_2:=z^*$ realizes the minima for
both $\sigma_1$ and $\sigma_2$.

Continuing in this fashion, we arrive at the conclusion that there exist linear combinations of the form
\begin{equation}
z^*_1 = z_1 + \sum_{\sigma \in G}u_{\sigma}z_{\sigma}
\end{equation}
where $u_{\sigma} \in ({\cal O}_K)^{\times}\cup \{0\}$ which satisfy the following four conditions:
\begin{enumerate}
\item $z_{1}^* \in \mathcal{B}_1$,
\item $K(z^*_1) \cong K(P)$,
\item $v(z^*_1) = \min\{ v(x_1),\dots, v(x_n)\}$,
\item $v(z^{*}_1 - \sigma(z^*_1)) = \min\{ v(x_1 - \sigma(x_1)),\dots, v(x_n - \sigma(x_g))\}$ for all $\sigma \in G$,
\end{enumerate}
as desired.
Again, we pause to note that the matrix $M$ representing the change of coordinates $(z_1^*,x_2,\dots,x_g)^t= M (x_1,x_2,\dots,x_g)^t$ is invertible.
Indeed, the matrix $M$ has units along the diagonal, the coefficients of the linear combination $z_1^*$ in the first row, and zeros elsewhere, hence the determinant is clearly a unit.

We now wish to iterate this construction as follows.
First, consider the set of all linear combinations
\begin{equation}\label{eqn:firstlinearcombo}
\mathcal{B}_2:=\{z' = u_1z_1^* + u_2x_2 + \cdots + u_gx_g \text{ where } u_i\in({\cal O}_K)^{\times} \cup \{0\} \text{ and }u_2 \neq 0 \}.
\end{equation}
Then, we can repeat the above construction to arrive at a linear combination $z_2 \in \mathcal{B}_2$ satisfing:
\begin{enumerate}
\item $v(z') \geq v(z_2)$ for all other linear combinations $z$ from $\mathcal{B}_2$,
\item $K(z_2) \cong K(P)$.
\end{enumerate}
Furthermore, we can follow the above construction to say that there exist linear combinations of the form
\begin{equation}
z^*_2 = z_2 + \sum_{\sigma \in G}u_{\sigma}z_{\sigma}
\end{equation}
where $u_{\sigma} \in ({\cal O}_K)^{\times}\cup \{0\}$ which satisfy the following four conditions:
\begin{enumerate}
\item $z_{2}^* \in \mathcal{B}_2$,
\item $K(z^*_2) \cong K(P)$,
\item $v(z^*_2) = \min\{ v(x_1),\dots, v(x_n)\}$,
\item $v(z^{*}_2 - \sigma(z^*_2)) = \min\{ v(x_1 - \sigma(x_1)),\dots, v(x_g - \sigma(x_g))\}$ for all $\sigma \in G$,
\end{enumerate}
Note that the matrix $M'$ representing the change of coordinates $(z_1^*,z_2^*,\dots,x_g)^t = M'(z_1^*,x_2,\dots,x_g)^t$ is invertible.
Indeed, $M'$ has units on the diagonal, the coefficients of the linear combination $z_2^*$ in the second row, and has zeros everywhere else.
Although this matrix is not triangular, we can make it so by switching the second row with the first and interchanging the first and second columns; these operations will not change the determinant. After these operations, the matrix becomes triangular with units on the diagonal, and hence the will be invertible.
Moreover, we see that the matrix $M''$ representing the change of coordinates $(z_1^*,z_2^*,\dots,x_g)^t = M''(x_1,x_2,\dots,x_g)^t$ is invertible since $M'' = M' \cdot M$.

We continue in this fashion for all of the remaining coordinates $x_3,\dots,x_g$ and arrive at our desired claim, namely that there exists linear combinations $z_1^*,\dots,z_g^*$ of $x_1,\dots,x_g$ with coefficients in $({\cal O}_K)^{\times} \cup \{0\}$ which satisfy:
\begin{enumerate}
\item $K(z_{i}^*) \cong K(P)$,
\item $v(z_{i}^*) =  \min\{ v(x_{1}),\dots, v(x_{g})\}$,
\item $v(z_{i}^* - \sigma(z_{i}^*)) = \min\{ v(x_{1} - \sigma(x_{1})),\dots, v(x_{g} - \sigma(x_{g}))\}$ for all $\sigma G$
\end{enumerate}
and such the matrix $M$ representing the change of coordinates $(z_1^*,\dots,z_g^*)^t= M (x_1,\dots,x_g)^t$ is invertible.
\end{proof}

We now complete the proof of  \autoref{thm:tamelyramifiedtorsion}.
\begin{proof}[Proof of \autoref{thm:tamelyramifiedtorsion}]
Fix $0\neq P = (x_1,\dots,x_g) \in {\mathscr F}[p]({\cal O})$ and let $e = [K(P) : K]$.
We first remark that since $K$ is completion of the maximal unramified extension, we have that $e > 1$, and hence the extension $[K(P) : K]$ is totally ramified. Indeed, this follows from that fact that the group-scheme ${\mathscr F}[p]$ is connected.
In \autoref{lemma:niceisom}, we constructed linear combinations $z_i^*$ of $x_1,\dots,x_g$ satisfying certain properties.
Let $Z_i^*$ denote the same linear combinations of the coordinates $X_1,\dots,X_g$ (so if we were to evaluate $Z_i^*$ at $(x_1,\dots,x_g)$ we would recover $z_i^*$).
The last condition from \autoref{lemma:niceisom} implies that the change of variables $(X_1,\dots,X_g) \mapsto (Z_1^*,\dots,Z_g^*)$ is an isomorphism of formal groups.

We claim that the isomorphism of formal groups $(X_1,\dots,X_g) \mapsto (Z_1^*,\dots,Z_g^*)$ will preserve strictness.
As each of the $Z_i^*$ are linear combinations of $X_1,\dots,X_g$ with coefficients in $({\cal O}_K)^{\times} \cup \{0\}$, this isomorphism of formal groups will act linearly on terms of minimal degree, and hence it changes $F_1,\dots,F_g$ by a linear transformation, which is invertible.
We note that it also does the same to $G_1,\dots,G_g$, therefore, it transforms the set of solutions of the system by an invertible transformation. Moreover, the system having or not having a single solution $(0,\dots,0)$ is the same before or after an isomorphism.
We pause to note that it is crucial that the degrees $d_1, d_2, \dots, d_g$ from \autoref{defn:strict} are all equal

For the remainder of the proof, we work with this isomorphic formal group with coordinates $(Z_1^*,\dots,Z_g^*)$.
We note that the vector $(z_1^*,\dots,z_g^*)$ will reduce mod $p$ to the point $(0,\dots,0) \in k^g$ because all $z_i^*$ have valuations strictly positive.
But the $z_i^*$ have the same valuation so we can divide all of them by one of them, and consider the vector $(z_1^*/z_g^*,\dots,z_{g-1}^*/z_g^*,1)$, which will not reduce the zero vector over the residue field.
Since ${\mathscr F}$ was assumed to be strict, the reduction of $(z_1^*/z_g^*,\dots,z_{g-1}^*/z_g^*,1)$ cannot be a common root of all of $G_1,\dots,G_g$. Therefore, there exists an index $j$ for which $G_j(z_1^*/z_g^*,\dots,1)$ is not zero in the residue field $k$, and hence the valuation of  $F_j(z_1^*,\dots,z_n^*)$ equals the valuation of each of its individual monomials.

We now want to determine the valuation of $z_j^*$, and hence the valuation of every other $z_i^*$ as they have the same valuation.
By considering the equation $f_j(z_1^*,\dots,z_g^*)$, we have
\[
0 = pz_j^* + p(\text{terms of degree between $2$ and $d_j - 1$}) + F_j(z_1^*,\dots,z_n^*) + (\text{higher degree terms}).
\]

We claim that $v(z_j^*) = 1/(d_j-1)$ where $d_j$ is the degree of $F_j$.
To see this choose a unit $u_1$ in ${\cal O}_K$ that is a representative for the element in the residue field corresponding to $z_1^*/z_g^*$, and similarly choose $u_2,\dots,u_{g-1}$. Then $F_j(u_1,\dots,u_{g-1},1)$ is a unit in ${\cal O}_K$, because its image in the residue field is nonzero, by the above choice of $j$.
But this is a form (of degree $d_j$) so we can divide by $u_j$ inside $F_j$, and re-denoting the $u_j$'s in consideration, we have that $F_j( u_1,\dots,u_{j-1}, 1, u_{j+1} ,\dots,u_g)$ is a unit, and moreover,
\[
F_j(z_1^*,\dots,z_g^*)  =  F_j(u_1,...1,...u_g) z_j^{*^{d_j}}  +  ( \text{terms of strictly larger valuation}).
\]
Plugging this into the the above equation, we arrive at the equation
\begin{equation}\label{eqn:firstpowerseries}
0 = pz_j^* + F_j(u_1,\dots,1,\dots,u_g) z_j^{*^{d_j}}  + (\text{terms of strictly larger valuation}).
\end{equation}
The minimum valuation in the equality of \eqref{eqn:firstpowerseries} must be attained in at least two terms, and these terms are forced to be $pz_j^*$ and $F_j(u_1,\dots,1,\dots,u_g) z_j^{*^{d_j}}$. Our claim now follows since $F_j(u_1,\dots,1,\dots,u_g) $ is a unit in ${\cal O}_K$.

We now want to study the relationship between the valuations $v(z_j^* - \sigma(z_j^*))$ where $\sigma \in G$.
Let $u := F_j(u_1,\dots,1,...u_g)$ which is a unit in ${\cal O}_K$.
For each $\sigma \in G$, we will consider \eqref{eqn:firstpowerseries} and
\begin{equation}\label{eqn:secondpowerseries}
0 = p\sigma(z_j^*) + u \sigma(z_j^{*})^{d_j} + (***)
\end{equation}
where $(***)$ corresponds to terms strictly larger valuation.
If we subtract equality \eqref{eqn:secondpowerseries} from \eqref{eqn:firstpowerseries}, we arrive at the following:
\begin{equation}\label{eqn:finalpoweseries}
0 = p(z_j^* - \sigma(z_j^*)) + u(z_j^{*^{d_j}} - \sigma(z_j^{*})^{d_j}) + (\text{interesting terms}).
\end{equation}
By condition (3) of \autoref{lemma:niceisom}, the valuation of the ``interesting terms" from \eqref{eqn:finalpoweseries} will be larger than $v(z_j^* - \sigma(z_j^*))$.
Indeed, these interesting terms are in fact monomials of the form $z_1^{*^{m_1}} z_2^{*^{m_2}} \dots z_n^{*^{m_g}}$ where some of the $m_i$ could be zero and the total degree is strictly greater than $d_j$.
In any case, we may deal with the difference of such monomials and their conjugates as follows. Suppose for example, that we have a term of the for $z_2^{*^{m_2}}z_3^{*^{m_3}} - \sigma ( z_2^{*^{m_2}}z_3^{*^{m_3}})$. Then by adding and subtracting the term  $z_2^{*^{m_3}} \sigma( z_3^{*^{m_3}})$, the difference we need to deal with can then be written as $(z_2^* - \sigma(z_2^*))$ times something of positive valuation, plus $(z_3^* - \sigma(z_3^*)$ times something of positive valuation, and so we can use property (3) of \autoref{lemma:niceisom} for this particular $\sigma$ to get our desired claim.

We now arrive at the crucial claim of the proof.
Recall that $G$ denotes the Galois group of the Galois closure of $K(P)/K$.
We claim that $v(z_j^* - \sigma(z_j^*)) = v(z_j^*)$ for all $\sigma \in G$.
Assume that $v(z_j^* - \sigma(z_j^*)) = v(z_j^*) + t$  where $t > 0$.
By the above discussion and condition (3) of \autoref{lemma:niceisom}, we can save the value $t$ from each term, and since there must be at least two terms of equal valuation in \eqref{eqn:finalpoweseries}, we have that
\[
v(p(z_j^* - \sigma(z_j^*))) \geq v(u(z_j^{*^{d_j}} - \sigma(z_j^{*})^{d_j})).
\]
Note that we cannot guarantee the equality of these valuations because there could be other terms of total minimal degree other than $u(z_j^{*^{d_j}} - \sigma(z_j^{*})^{d_j})$, but the above inequality will suffice.
Using the above inequality and previous equality $v(pz_j^*) = v(uz_j^{*^{d_j}})$, and letting $w = \sigma(z_j^*)/z_j^*$, we have that
\[
t = v(z_j^* - \sigma(z_j^*))  - v(z_j^*) = v((z_j^* - \sigma(z_j^*))/z_j^*) = v(1 - w) \geq v(1 - w^{d_j}).
\]
Let $y = 1 - w$. We have that $v(y) = t$, which by assumption is strictly greater than 0.
We now arrive at a contradiction by considering the above inequality and the equation
\[
1 - w^{d_j} = 1 - (1 - y)^{d_j} = d_jy + (\text{terms times }y^2) + y^{d_j},
\]
and noting that all terms in the above have valuation strictly greater than $v(y) = t$.
Therefore, we have that $t= 0 $, and hence $v(z_j^* - \sigma(z_j^*)) = v(z_j^*)$ for all $\sigma \in G$.

To conclude our proof, we use Krasner's lemma \cite[Exercise II.2.1]{Serre} to explicitly describe the extension $K(z_j^*)/K$ and show that $[K(z_j^*) : K] = d_j - 1$.
Recall that our $z_j^*$ satisfies the equation \eqref{eqn:firstpowerseries}.
Consider the polynomial $P(Z) = p + uZ^{d_j - 1}$ where $u = F_j(u_1,...1,...u_n)$ as above.
Note that $z_j^*$ is not a root of $P(Z)$, but it satisfies the following inequality $v(P(z_j^*)) > v(pz_j^*) $ where the right side here is also equal to $v( u z_j^{*^{d_j}})$.

On the other hand, the roots $\theta_1,\dots,\theta_{d_j - 1}$ of $P(Z)$ each have valuation exactly $1/(d_j - 1)$ since $P(Z)$ is Eisenstein at $p$.
We now compute the valuation of the derivative of $P(Z)$ evaluated at a root in two different ways.
First, $P'(\theta_l) = u \prod_{i\neq l} (\theta_l - \theta_i)$ and hence
\[
v(P'(\theta_l)) = \sum_{i\neq l} v(\theta_l - \theta_i) \geq \frac{d_j - 2}{d_j - 1}.
\]
Second, we directly compute that $P'(\theta_l) = (d_j - 1)u\theta_l^{d_j - 2},$ which yields
\[
v(P'(\theta_l)) = \frac{d_j - 2}{d_j - 1}.
\]
The first inequality and the second equality imply that $v(\theta_l - \theta_i) = 1/(d_j - 1)$ for all $1 \leq i \neq j \leq d_j - 1$.
We also have that $P(z_j^*) = u(z_j^* - \theta_1)\cdots (z_j^* - \theta_{d_j - 1})$ and hence
\[
v(P(z_j^*)) = \sum_{i = 1}^{d_j -1} v(z_j^* - \theta_i).
\]
There are exactly $d_j - 1$ terms here and since $v(P(z^*)) > 1$, it follows that at least one term must be strictly larger than $1/(d_j - 1)$.
Without loss of generality, we may assume that $v(z_j^* - \theta_1) > 1/(d_j - 1).$
Now we have the strict inequality
\[
v(z_j^* - \theta_1) > 1/(d_j - 1) = v(\theta_1 - \theta_i)
\]
for all $1 < i \leq d_j - 1$, and hence Krasner's lemma implies that $K(\theta_1) \subseteq K(z_j^{*})$.
Recall that we have just shown that $v(z_j^* - \sigma(z_j^*)) = v(z_j^*) = 1/(d_j - 1)$ for all $\sigma \in G$.
Moreover, we can just switch $z_j^*$ and $\theta_1$ to arrive at the inequality
\[
v(\theta_1 - z_j^*) > v(z_j^* - \sigma(z_j^*)) = 1/(d_j - 1),
\]
and so we can apply Krasner's lemma again to deduce that $K(z_j^*) \subseteq K(\theta_1)$.
Therefore, we have shown that $K(P)\cong K(z_j^*)\cong K(\theta_1)$ and this extension is totally ramified of degree $d_j- 1$, hence tamely ramified. We also have that $z_j^*$, which is a linear combination of $x_1,\dots, x_g$ with unit coefficients, is a uniformizer for $K(P)$ and hence ${\cal O}_{K(P)} \cong {\cal O}_K[x_1,\dots,x_g]$. Finally, we note that there exists some coordinate $x_i$ of $P$ which has valuation $v(x_i) = v(z_j^*)$ by condition (3) of \autoref{lemma:niceisom} and hence $x_i$ is a uniformizer for $K(P)$ as well.

The second statement of \autoref{thm:tamelyramifiedtorsion} follows because $K({\mathscr F}[p])$ is the compositum of the fields $K(P)/K$ as $P$ varies over points in ${\mathscr F}[p]$ and the compositum of tamely ramified extensions is again tamely ramified.
\end{proof}


\section{Introduction}

Gaseous systems seen in absorption along quasar sightlines are excellent probes of astrophysics throughout cosmic time, from measuring the physical and chemical evolution of gaseous reservoirs of galaxies \citep{Berg15,Neeleman15,Noterdaeme21} to imposing constraints on dark matter properties \citep{Irsic17} and fundamental constant evolution in our Universe \citep[e.g.,][]{Murphy17,Schmidt21}. Accurately quantifying these probes requires a detailed kinematic decomposition of all absorption lines along a quasar sightline. For distinct QSO absorption line systems  (QALs),  the metal absorption profiles are often quite complex, consisting of several sub-components (with a typical width between $\approx5-100$~km~s$^{-1}${}) spread out over a few hundred km~s$^{-1}${}. In addition, the neutral hydrogen absorption of individual absorbers are sometimes significantly blended within the Ly$\alpha$ forest absorption from the intervening intergalactic medium. While low resolution spectra (with resolving power $R<1000$) can easily identify the strongest Ly$\alpha$ QALs such as damped Lyman $\alpha$ systems \citep[DLAs; H\ion{i} column densities of $\geq 2\times10^{20}~{\rm atom~cm^{-2}}$;][]{Wolfe05}, spectrographs like UVES \citep[][]{Dekker00} are frequently used to  identify the weaker H\ion{i} absorbers, which are required to probe the intergalactic and circumgalactic media. Typical studies of faint quasars with UVES use a $1''$ slit ($R\approx48~000$) to maximize $R$ while minimizing light loss based on the median seeing of Paranal ($\approx0.8''$).

Although UVES has proved to be an effective instrument for the study of QALs over the past 20 years, the QAL community are now exploring higher spectral resolving power ($R\gtrsim100~000$; i.e., the smallest cloud widths of $\gtrsim3$~km~s$^{-1}${}) to pinpoint the detailed structure of the absorption line profiles, and search for rarely detected species, such as:  (i) deuterium and lithium for constraining Big Bang nucleosynthesis \citep{Cooke18}, (ii) molecular hydrogen to constrain the gas and dust properties of high redshift galactic gas reservoir \citep{Tchernyshyov15, deCia18, Krogager18}, and (iii) exotic nucleosynthetic species to place constraints on chemical evolution processes \citep{Ellison02,Berg13}. To do these experiments with UVES would require the smallest slit widths (0.3$''$; thus $R\approx100~000$), leading to large amounts of light loss at the slit interface unless observations are completed in the most pristine atmospheric seeing conditions. In essence, it is not efficient for UVES to continue to explore the emerging questions of the QAL community.

The ESPRESSO spectrograph \citep{Pepe10,Pepe21} provides an excellent alternative to UVES for observing QALs at high spectral resolution.  In its medium-resolution mode, the 4 Unit Telescope (UT) MULTIMR42 mode ($R=75~000$ with 4x2 binning), ESPRESSO has already proven to be effective at using QALs to constrain our understanding of the Universe \citep{Cooke20,Welsh20}. However, these observations require all 4 UTs simultaneously, and the collection of light from MULTIMR42 mode means a reduction of $R$ from $\approx 137\,000$ to $75\,000$; a resolution UVES can obtain without excessive slit losses in  average Paranal atmospheric conditions. As of Period P108, ESPRESSO supports a single UT mode with $4\times2$ binning (HR42 mode). This mode is designed to reduce read-out noise while maintaining a high spectral resolving power ($R\approx137~000$) and is ideal for observing faint targets such as QSOs.

In this paper, we present the results from a pilot program to test the ability of ESPRESSO's HR42 mode to analyze the chemical contents of a DLA seen along a quasar sightline, and compare the performance of ESPRESSO HR42 mode to the typical observing set-up used for UVES.

\section{Target selection, observations and data reduction}

The target, QSO~J0003$-$2603 (right ascension 00~hr $03'$ $22.953''$, declination $-26^\circ$ $03'$ $18.307''$, and redshift $z_{\rm em}=4.01$), was chosen as it is: ($i$) sufficiently bright for the ESPRESSO autoguiding system (AB magnitude $R = 18.5$~mag),  ($ii$) known to host an intervening metal-poor ([M/H]~$\approx-2$) DLA absorber along the line of sight (redshift $z_{\rm abs}=3.39$, and ($iii$) has archival UVES observations allowing a direct comparison between instruments \citep[][]{Molaro00,Levshakov00}.

The ESPRESSO observations were carried out on October 2$^{\rm nd}$ 2021 on UT2 at Paranal (Program ID 60.A-9801(W)). Four sequential 60-minute on-target exposures were taken. The observations were made in average seeing conditions ($\approx0.77''$; image quality at 5000~\AA{} of $\approx0.85''$) across the four hours of observations at an average airmass of 1.05. The data were reduced with the default settings of the ESPRESSO ESO-reflex pipeline \citep[version 2.3.3;][]{Pepe21,Freudling13} and post-processed with the ESPRESSO data analysis software \citep[hereafter DAS; version 1.3.3;][]{Cupani19,Pepe21}. While we use the sky-subtracted data in this paper to avoid skyline contribution, we point out that the S/N is marginally better in the spectrum without sky subtraction. The S/N (in units of pixel$^{-1}$; $\approx1$~km~s$^{-1}${} per pixel) is $\approx20$ within the Ly$\alpha$ forest, and $\approx30$ beyond the Ly$\alpha$ emission line of the quasar.

The quoted resolving power of the HR42 mode, based on the ThAr arc calibrations taken by ESO, is $R\approx131~000$. To compare what was achieved by our data, we quantified $R$ using the telluric lines from two different methods. We first assumed that the measured full width half maximum (FWHM) of a weak, unresolved telluric line is equal to the FWHM of the resolution element (${\rm FWHM_{res}}${}); the measured ${\rm FWHM_{res}}${} of the telluric line at $\approx7685$~\AA{} within the median-combined spectrum is $2.269\pm0.292$~km~s$^{-1}${} (corresponding to $R\approx132~100$). However, this method assumes that telluric lines are unresolved by ESPRESSO, which is not necessarily true given the high spectral resolution of the HR42 mode. To confirm the specific telluric line was unresolved and to obtain a second, independent measurement of the resolving power, we fitted the telluric lines using using \textsc{molecfit} \citep{Smette15}, which provides both an estimate of ${\rm FWHM_{res}}${} as well as the intrinsic FWHM of the telluric line based on a model atmosphere at the time of observations. Taking the mean value from each individual exposure, the \textsc{molecfit}-estimated intrinsic FWHM of the telluric line at $\approx7685$~\AA{} is $1.480\pm0.050$~km~s$^{-1}${} and ${\rm FWHM_{res}}${}~$=2.193\pm0.050$~km~s$^{-1}${} ($R\approx136~700$). It is clear from the results of \textsc{molecfit} that the telluric line at $\approx7685$~\AA{} is marginally unresolved within our dataset. However we caution that depending on the airmass and atmospheric conditions at time of observations, the telluric lines may be resolvable by ESPRESSO's HR mode and will not provide an accurate estimate of the resolving power of the instrument. We adopt $R\approx136~700$ from \textsc{molecfit} as the best estimate of the resolving power of the HR42 mode.

The comparison archival UVES data used 4.5 hours of on-target exposures, and was combined and continuum normalized by the UVES Spectral Quasar Absorption Database team \citep[SQUAD;][]{Murphy19}. Observations were executed during the commissioning of UVES (Program ID 60.A-9022(A); slit width 0.9$''$ with the DIC2 437+860 instrument setup), and were observed at similar airmass and at various seeing measurements (from 0.34$''$ to 2.4$''$; mean of $1.0''$ at 5000~\AA{}) across two nights. Using unresolved telluric lines within the data, the obtained FWHM of the spectral resolution element is ${\rm FWHM_{res}}${}~$=5.981\pm0.126$~km~s$^{-1}${} ($R\approx50~100$; archival $R$ measurements from ThAr arc lines obtained $R\approx47~900$). The UVES pixel size is $\approx2.5$~km~s$^{-1}${}. The S/N  per pixel is $\approx2\times$ (alternatively, $\approx 1.33\times$~per \AA{}) larger within the UVES spectrum relative to that of ESPRESSO (Figure \ref{fig:snr}). Once accounting for the difference in pixel size between the two instruments, the S/N of UVES becomes significantly better than ESPRESSO at $\lesssim5000$~\AA{}.

\section{Results}
\subsection{Comparison with the ESO Exposure Time Calculator}

Figure \ref{fig:snr} shows the observed S/N (red line) across the entire spectrum compared to what was estimated using version P109 of the ESPRESSO exposure time calculator (ETC; black points). The ETC S/N estimates have been computed assuming a $0.84''$ image quality, the average airmass through the observations (1.05), and fractional lunar illumination FLI=0.0 as the moon was not visible during observations. The S/N estimated by the ETC is consistently $\approx1.3\times$ higher across the entire wavelength range. However, we point out that the ETC solely includes emission from the QSO continuum and not absorption lines observed towards QSO J0003$-$2603 (particularly those of the Ly$\alpha$ forest; $<6202$~\AA{}), thus overestimating the estimated S/N. This is particularly demonstrated by the presence of the strong Ly$\alpha$ absorption from the DLA at $z_{\rm abs}=3.39$ ($\approx5300$~\AA{} in red line of Figure~\ref{fig:snr}). Using only pixels of the observed ESPRESSO spectrum with predominantly quasar continuum flux near the centers of ESPRESSO's spectral orders, the observed S/N is only $\approx1.15\times$ smaller than the ETC predicts.

To obtain a near-equivalent spectrum with UVES (i.e., $R\approx110~000$) and wavelength coverage from 3800~\AA{} to 7880~\AA{}), one would need to observe with two UVES settings (e.g., DIC1 390+580 and DIC2 437+760) and with slit widths of 0.3$''$. The cyan points in Figure \ref{fig:snr} show the ETC-predicted S/N for the UVES spectrum using these two settings and 0.3$''$ slit width, with an exposure time of four hours \emph{per setting} under the same observing conditions as the ESPRESSO  observations. The S/N of the cyan points is $\approx 2\times$ higher per pixel than what is observed with ESPRESSO for two reasons. First, as there is overlap in wavelength between the two UVES settings, there is an increase in the S/N in the wavelength ranges $\approx$3760-4480~\AA{},  $\approx$4770-4960~\AA{}, and $\approx$5700-6780~\AA{}. Secondly, the UVES unbinned pixel sampling of the resolution element of this setup is 1.39 pixels per element in the blue and 1.94 pixels in the red. Not only does this pixel size lead to under-sampling in the spectral direction, but also it results in an increase in S/N per pixel of 1.18--1.39$\times$ compared to ESPRESSO. Thus the expected ratio of the S/N between UVES and ESPRESSO -- after accounting for the doubled exposure time in overlapping regions of the two UVES settings, and the resolution element size -- is between 1.18 and 1.97, and explains the majority of the discrepancy between the red curve and cyan points of Figure \ref{fig:snr}. While ESPRESSO is expected to be as efficient as UVES with a 0.3" slit to obtain a $R>100~000$ spectrum, ESPRESSO is able to obtain full wavelength coverage (3800\AA{} to 7880\AA{}) in a single setup whilst Nyquist sampling the spectral resolution element, making the HR42 mode more efficient in terms of observing time required.

\begin{figure*}
    \centering
    \includegraphics[width=\textwidth]{SNR.pdf}
    \caption{The S/N (per pixel; $\approx 1$~km~s$^{-1}${} per pixel) of the observed ESPRESSO spectrum (solid red line) and ETC-estimated values based on a QSO emission template (black circles). The observed spectrum S/N is estimated as the measured flux divided by the standard deviation from the error spectrum. Note that the red line is smoothed by taking the median S/N within a wavelength bin of the same size as the average space between ETC points. For reference, the vertical grey line denotes the Ly$\alpha$ emission line of the observed QSO. The S/N from the ETC consistently $\lesssim20$~percent larger than measured in the observations. For reference, the smoothed S/N from the comparison UVES data ($\approx2.5$~km~s$^{-1}${} per pixel) is shown as the blue dotted line. Due to the gaps in the UVES data wavelength coverage, the observed UVES data does not cover the Ly$\alpha$ emission of the QSO. Additionally, the small cyan dots represent the combined ETC estimate for UVES data taken with two UVES settings, with each setting using a 0.3$''$ slit (near equivalent $R$) with the same exposure time \emph{per setting} (i.e., doubled observing time) and observing conditions as the ESPRESSO data.}
    \label{fig:snr}
\end{figure*}

\subsection{Characterization of the intervening absorber with ESPRESSO}

By simultaneously fitting the continuum and a Voigt profile to both the Ly$\alpha$ and Ly$\beta$ lines, the measured H\ion{i} column density of the  intervening absorber towards QSO~J0003$-$2603 is logN(H\ion{i})~$=21.35\pm 0.12$ at $z_{\rm abs}=3.3901$. The measured H\ion{i} column density is consistent with previous measurements \citep{Lu96}. The Voigt profile fit to the Ly$\alpha$ absorption line is shown in Figure \ref{fig:HI}.

\begin{figure}
    \centering
    \includegraphics[width=0.5\textwidth]{HI.pdf}
    \caption{The best-fit Voigt profile of Ly$\alpha$ (cyan line) to the continuum-normalized ESPRESSO data (black line) for the absorption seen at $z_{\rm abs}=3.3901$. The best fitting H\ion{i} column density is logN(H\ion{i})~$=21.35\pm 0.12$. The shaded cyan region denotes the range of profiles from the error on logN(H\ion{i}). }
    \label{fig:HI}
\end{figure}

Additionally, Voigt profile fitting was done for absorption lines of key metal species (C\ion{iv}, N\ion{i}, Si\ion{ii}, Al\ion{ii}, Fe\ion{ii}, Ni\ion{ii}). We note that in the analysis of the UVES data done by \cite{Molaro00}, there was sufficient S/N in the Ly$\alpha$ forest ($\ge10$~pixel$^{-1}$) to detect the O\ion{i} lines $\lambda$ 925~\AA{} and 950~\AA{} ($\approx 4115$~\AA{} in the observed frame). However, the efficiency of ESPRESSO at blue wavelengths is too poor to be able to accurately model the continuum within the Ly$\alpha$ forest, although it is possible to see evidence of saturation of the O\ion{i} $\lambda$ 929~\AA{} and 950~\AA{} lines. Unfortunately, the wavelength coverage of ESPRESSO does not include coverage of  Cr\ion{ii} and Zn\ion{ii} at the redshift of the DLA.

The fitting was done using the code \textsc{ALIS}\footnote{\url{https://github.com/rcooke-ast/ALIS}}, which uses a $\chi^{2}$ minimization to simultaneously fit the continuum and profile components from a list of starting redshifts ($z$), the turbulent component  of the Doppler broadening parameters ($b_{\rm turb}${}), gas temperature (${\rm T_{gas}}${}), and column densities (N; in units atoms cm${-2}$) for each ionic species. A Gaussian FWHM of 2.193 km~s$^{-1}${} is used as the ESPRESSO line spread function in the fitting procedure. As the HR42 mode provides sufficient resolution to possibly constrain gas temperature \citep{Narayanan06}, we attempted including ${\rm T_{gas}}${} as a free parameter per absorption component in the fitting process. Despite having a range of ions of varying atomic mass, many lines in this system  are blended or saturated preventing a realistic constraint on ${\rm T_{gas}}${}; including ${\rm T_{gas}}${} as a free parameter in the model suggests a purely turbulent medium. In order to obtain an accurate measure of the column densities for all components of each ion, we elect to assume ${\rm T_{gas}}${}~$=10~000$~K, a typical temperature of metal-poor DLAs \citep{Cooke15} rather than separately fitting a single Doppler parameter to each ions' component. While assuming a single ${\rm T_{gas}}${} may not provide accurate measurements of $b_{\rm turb}${} for each component, the Voigt profile fitting procedure will be more precise as the fit is constrained by more data (i.e.~multiple ions' probing the same $b_{\rm turb}${}) and effects from atomic mass are properly accounted for. As a result, this will reduce the overall errors in both the Doppler parameter and logN. Furthermore, by deriving $b_{\rm turb}${} for each component, the fitting also breaks the degeneracy between the Doppler parameter and column density for saturated lines, enabling a measurement of the column density for Si\ion{ii} and Al\ion{ii}. We note that repeating the fitting procedure with a ${\rm T_{gas}}${} within the range of DLA values (i.e.~between 1000K -- 100~000K) has little effect on the derived column densities, which remain consistent with the assumed ${\rm T_{gas}}${}$=10~000$~K column densities. By assuming a fixed ${\rm T_{gas}}${}, the errors derived from fitting $b_{\rm turb}${} encapsulate the errors on the total Doppler parameter. The high ions (i.e., C\ion{iv}) and low ions (N\ion{i}, Si\ion{ii}, Al\ion{ii}, Fe\ion{ii}, and Ni\ion{ii}) were fitted separately, and assuming that the $z$ of each component of the low ions' velocity profile are the same across all species. Given the level of blending and saturation for many of these lines which introduces model degeneracies while fitting all lines simultaneously, the velocity profiles were fitted in several iterations to constrain individual components free of these degeneracies. First, the strongest components at $\approx0$~km~s$^{-1}${} were fitted using the unsaturated Fe\ion{ii} $\lambda$ 1611~\AA{} and the two redmost lines of the N\ion{i} $\lambda$ 1134~\AA{} triplet  ($\lambda$ 1134.4 and 1135.0~\AA{}) to fix $b_{\rm turb}${}, $z$ and N(Fe\ion{ii}) for these two components that cannot be easily determined from the saturated lines. The redshifts of the rest of the profile (five components total) were then fixed using the unblended Fe\ion{ii} $\lambda$ 1608~\AA{}, before running the fit on the remaining lines to determine the remaining parameters. These ALIS models were subsequently cross-checked with the ESPRESSO DAS; the results of the two methods are in agreement within the uncertainties.  The best-fit absorption profile parameters are presented in the top of Table \ref{tab:loComp}, while the velocity profiles and their corresponding fits are shown in the two leftmost columns of Figure \ref{fig:lowions}. For reference, the  velocity of each component relative to $z_{\rm abs}$ ($\Delta v$) is tabulated in Table \ref{tab:loComp}.

\begin{figure*}
\centering
\includegraphics[width=\textwidth]{J0003-2603_lowions_velplots.pdf}
\caption{Velocity profiles of key low-ionization absorption lines (labelled in the bottom right corner of each panel) detected within the continuum-normalized data (black lines) of QSO~J0003$-$2603. The panels of the left two columns of the figure denote data from ESPRESSO, while the rightmost column shows the absorption lines of the UVES data. The solid magenta curve denotes the complete Voigt profile fitted to the data, while the red dashed lines denote the components that constitute the fit. Overplotted in cyan are the blended components. Although the profile for N\ion{i} $\lambda$ 1134.1~\AA{} was not fit due uncertain continuum blueward of the line, the expected profile shape based on the fits of the other two triplet lines ($\lambda$ 1134.4~\AA{} [1134b] and $\lambda$ 1135.0~\AA{} [1134c]) is shown in the N\ion{i} 1134b panel at $\approx-70$~km~s$^{-1}${}.}
\label{fig:lowions}
\end{figure*}

 Similarly, we applied the same profile fitting method to the low ions in the UVES data (assuming an instrumental resolution element of 5.981 km~s$^{-1}${}; right column of Figure \ref{fig:lowions}; bottom of Table \ref{tab:loComp}), and detected only four components. Components 2, 3 and 4 of the ESPRESSO fit (Table \ref{tab:loComp}) are separated by $\approx20$~km~s$^{-1}${}, and all three components have a $b_{\rm turb}${} below the resolving power of the 0.9$''$ UVES set-up (FWHM~$\approx6.4$~km~s$^{-1}${}). As component 2 is weaker than the other two components, the asymmetry in the line profile can be detected within the UVES data, while components 3 and 4 cannot be distinguished. Thus, as a result of the $\approx3\times$ larger spectral resolving power of the HR42 mode, more components with small $b_{\rm turb}${} can be detected with ESPRESSO despite the modest S/N.

 Given the strength of the Si\ion{ii} $\lambda$ 1526 line in the ESPRESSO data, we point out that there appears to be an artefact in the UVES data near the Si\ion{ii} $\lambda$ 1808~$\AA$ line. The profile shape appears shifted blueward in velocity space ($\approx -5$~km~s$^{-1}${}) with reference to the other lines (right column of Figure \ref{fig:lowions}), leading to a poor fit using the fixed redshift and $b_{\rm turb}${} from the Al\ion{ii} and Fe\ion{ii} lines. The resulting logN(Si\ion{ii})  is severely inconsistent with the ESPRESSO value. Fitting a separate, independent two-component profile to the Si\ion{ii} 1808 line results in a total column density estimate of logN(Si\ion{ii})~$=15.65$, which is inconsistent with the sum of the components 2, 3 and 4 of the equivalent absorption feature in the ESPRESSO data (logN(Si\ion{ii})~$=15.07\pm0.16$).

 \begin{table*}
 \small
\begin{center}
\caption{Component Voigt profile fits for low ions}
\label{tab:loComp}
\begin{tabular}{lcccccccc}
\hline
Comp.& $z$& $\Delta v$ (km s$^{-1}$)& $b_{\rm turb}${} (km s$^{-1}$)& logN(N\ion{i})& logN(Si\ion{ii})& logN(Al\ion{ii})& logN(Fe\ion{ii})& logN(Ni\ion{ii})\\
\hline
\multicolumn{9}{c}{\textbf{ESPRESSO}}\\
1& $3.389512$& $-42$& $8.1\pm0.27$& $13.25\pm0.175$& . . .{}& $12.14\pm0.017$& $13.38\pm0.058$& $12.59\pm0.271$\\
2& $3.389892$& $-16$& $4.3\pm0.09$& $13.30\pm0.082$& $13.71\pm0.068$& $12.30\pm0.017$& $13.59\pm0.015$& . . .{}\\
3& $3.390056$& $-5$& $5.1\pm0.12$& $14.26\pm0.021$& $14.56\pm0.198$& $13.17\pm0.068$& $14.57\pm0.026$& $13.14\pm0.045$\\
4& $3.390188$& $4$& $5.1\pm0.09$& $14.32\pm0.019$& $14.73\pm0.053$& $12.86\pm0.030$& $14.07\pm0.015$& $12.92\pm0.052$\\
5& $3.390481$& $24$& $13.7\pm0.68$& . . .{}& $13.29\pm0.028$& $12.06\pm0.024$& $12.97\pm0.062$& $11.96\pm0.911$\\
\hline
\multicolumn{9}{c}{\textbf{UVES}}\\
1& $3.389511$& $-42$& $8.4\pm0.39$& . . .{}& . . .{}& $12.09\pm0.054$& $13.40\pm0.019$& . . .{}\\
2& $3.389949$& $-12$& $7.0\pm0.11$& . . .{}& $14.96\pm0.009$& $12.79\pm0.012$& $13.99\pm0.008$& . . .{}\\
3& $3.390130$& $0$& $6.5\pm0.09$& . . .{}& $14.59\pm0.018$& $13.29\pm0.029$& $14.63\pm0.018$& . . .{}\\
4& $3.390505$& $26$& $15.4\pm1.17$& . . .{}& . . .{}& $11.92\pm0.035$& $12.82\pm0.086$& . . .{}\\

\end{tabular}
\end{center}
\end{table*}

The same fitting method was also applied to the C\ion{iv} doublet velocity profiles  (Figure \ref{fig:CIV}). The complexity of the C\ion{iv} absorption profile combined with either the modest S/N of the ESPRESSO or the resolution of the UVES data leads to a degenerate parameter space, requiring anywhere between 10 and 16 components to provide a satisfactory fit. Many of the additional components beyond the initial 10 are much weaker, with either the column density being several orders of magnitude smaller than the total value and thus consistent with the noise, or the components are shallow and/or broad in shape and could be accounted for by modifying the continuum fit. In many cases, these additional components contribute minimally to the total column density. For the UVES data, the minimum 10 components were sufficient to reproduce the absorption profile (right column of Figure \ref{fig:CIV}). As with the low ion profile, the ESPRESSO data is able to resolve an additional two components (at $\Delta v = 22$ and 57~km~s$^{-1}${}; left column of Figure \ref{fig:CIV}).  The properties of each component for both datasets are tabulated in Table \ref{tab:hiComp}. The recovered total column density of C\ion{iv} from UVES and ESPRESSO are equivalent.

\begin{figure*}
\centering
\includegraphics[width=\textwidth]{J0003-2603_CIV_velplots.pdf}
\caption{ Velocity profiles of the C\ion{iv} doublet absorption seen in the ESPRESSO (left column) and UVES (right column) data of QSO~J0003$-$2603. Notation is the same as in Figure \ref{fig:lowions}.}
\label{fig:CIV}
\end{figure*}

\begin{table*}
\begin{center}
\caption{Component Voigt profile fits for high ions}
\label{tab:hiComp}
\begin{tabular}{lcccc}
\hline
Comp.& $z$& $\Delta v$ (km s$^{-1}$)& $b$ (km s$^{-1}$)& logN(C\ion{iv})\\
\hline
\multicolumn{5}{c}{\textbf{ESPRESSO}}\\
1& $3.388100$& $-139$& $11.5\pm0.17$& $13.96\pm0.017$\\
2& $3.388645$& $-101$& $41.4\pm1.63$& $13.76\pm0.023$\\
3& $3.388729$& $-96$& $7.5\pm0.55$& $13.06\pm0.039$\\
4& $3.389040$& $-74$& $12.8\pm0.78$& $13.31\pm0.036$\\
5& $3.389703$& $-29$& $5.9\pm0.18$& $13.22\pm0.012$\\
6& $3.389918$& $-15$& $5.4\pm0.44$& $13.13\pm0.067$\\
7& $3.390172$& $3$& $11.9\pm1.04$& $13.91\pm0.035$\\
8& $3.390453$& $22$& $8.5\pm1.42$& $13.56\pm0.122$\\
9& $3.390717$& $40$& $10.3\pm1.78$& $13.87\pm0.082$\\
10& $3.390960$& $57$& $8.3\pm1.68$& $13.55\pm0.156$\\
11& $3.391181$& $72$& $9.4\pm1.09$& $13.52\pm0.084$\\
12& $3.391497$& $93$& $16.2\pm0.87$& $13.48\pm0.031$\\
\hline
\multicolumn{5}{c}{\textbf{UVES}}\\
1& $3.388094$& $-139$& $11.2\pm0.11$& $13.95\pm0.005$\\
2& $3.388660$& $-100$& $43.3\pm1.05$& $13.84\pm0.015$\\
3& $3.388702$& $-98$& $4.3\pm0.42$& $12.86\pm0.028$\\
4& $3.389022$& $-76$& $14.3\pm0.68$& $13.31\pm0.029$\\
5& $3.389681$& $-31$& $4.0\pm0.22$& $13.06\pm0.015$\\
6& $3.390156$& $2$& $1.6\pm0.66$& $12.89\pm0.076$\\
7& $3.390182$& $4$& $22.7\pm0.64$& $14.09\pm0.013$\\
8& $3.390781$& $44$& $19.7\pm0.83$& $14.10\pm0.019$\\
9& $3.391200$& $73$& $11.5\pm0.68$& $13.54\pm0.043$\\
10& $3.391548$& $97$& $12.8\pm0.49$& $13.38\pm0.021$\\
\end{tabular}
\end{center}
\end{table*}

\subsection{Comparison of DLA properties between ESPRESSO and UVES}

Table \ref{tab:Ntot} provides the total column densities derived from the ESPRESSO and UVES Voigt profile fitting analysis. Errors on the total column densities include continuum fitting errors. Comparing the total column densities derived from ESPRESSO with those from UVES \citep[both from the analysis in this paper and ][]{Molaro00}, all values are consistent. We note that \cite{Molaro00} only fitted the low ions absorption profile as a single component profile, thus excluding the weak components 1, 2, and 5 from the ESPRESSO fit (Table \ref{tab:loComp}; components 1 and 4 from UVES). However, we emphasize that these three weak components contribute minimally to the total column densities, and our derived total column densities are consistent with \cite{Molaro00}. Thus, while in Table \ref{tab:Ntot} we conservatively assume logN(Si\ion{ii}) is a lower limit as components 1 and 2 are severely blended, it is likely that they have little contribution to the total column density of Si\ion{ii}. Ignoring these missing components would imply a total logN(Si\ion{ii})~$=14.99\pm0.217$ from ESPRESSO.

Despite the $\approx2\times$ lower S/N in the ESPRESSO data, the errors on the ESPRESSO total ion column densities are only $\approx2\times$ higher than those of UVES. However, we emphasize that the number of components used in the fitting of the low ions and C\ion{iv} profiles are not identical, preventing a fair comparison of the errors due to the degeneracy in component structure. Taking the redshifts of the fitted components for the ions Al\ion{ii}, Fe\ion{ii}, and C\ion{iv} from ESPRESSO and fixing them in the UVES profile fitting results in both $b_{\rm turb}${} and logN being generally consistent when the same number of components are used in the fitting (Tables \ref{tab:NloFixEsp} and \ref{tab:NhiFixEsp}; for low ions and C\ion{iv} respectively). We exclude Si\ion{ii} from the comparison due to the artefact seen in the UVES data. Figure \ref{fig:errs} shows the distribution of errors on $b_{\rm turb}${} (top panel) and logN (bottom panel) for this alternative fit. While the errors on the UVES fits typically vary  between $\approx 0.25-4.0\times$ that of ESPRESSO, the distribution of errors are comparable in size between the two datasets. Therefore, the combination of lower S/N but higher resolution does not significantly change the precision of fitting the same absorption line profile model. However, the combination of complex absorption profiles with either the lower S/N of the ESPRESSO data or lower resolution of UVES implies that the parameter space of the fitting process becomes highly degenerate and thus $\chi^2$ minimization techniques may not converge on similar profiles fitted.

\begin{figure}
    \includegraphics[width=0.5\textwidth]{err_dist.pdf}
    \caption{Distributions of the errors on the fitted components' $b_{\rm turb}${} (top panel) and logN (bottom panel) for ions Al\ion{ii}, Fe\ion{ii}, and C\ion{iv} from ESPRESSO (solid red bars) and UVES (hollow black bars) fits assuming the same number of components as the ESPRESSO fits for the respective ions (Tables \ref{tab:loComp} and \ref{tab:hiComp}). The UVES distributions show the errors derived from fitting the same redshift components as used in the ESPRESSO fits in order to provide an equal comparison. Despite the lower S/N in the ESPRESSO data, the errors on the component fits for ESPRESSO are similar in magnitude as those of UVES.}
    \label{fig:errs}
\end{figure}

\begin{table*}
\small
\begin{center}
\caption{Total column densities}
\label{tab:Ntot}
\begin{tabular}{lcccccc}
\hline
Instrument& logN(C\ion{iv})& logN(N\ion{i})& logN(Si\ion{ii})& logN(Al\ion{ii})& logN(Fe\ion{ii})& logN(Ni\ion{ii})\\
\hline
ESPRESSO& $14.70\pm0.021$& $14.63\pm0.015$& $>15.01$& $13.44\pm0.040$& $14.75\pm0.018$& $13.42\pm0.060$\\
UVES& $14.71\pm0.007$& . . .{}& $>15.12$& $13.47\pm0.022$& $14.76\pm0.015$& . . .{}\\

\end{tabular}
\end{center}
\end{table*}

The two sets of data span two decades and we checked for a temporal evolution in the column density of the individual profiles' components by refitting the UVES Fe\ion{ii} and Al\ion{ii} lines with the same five components used to fit the ESPRESSO data. While keeping $z$ and $b_{\rm turb}${} fixed, logN was refit for the UVES data (Tables \ref{tab:Nevolve}). No significant evolution in column densities was detected. Within $1\sigma$, all logN(Fe\ion{ii}) for each fitted component in the UVES data were consistent with the equivalent component in ESPRESSO (Table \ref{tab:loComp}), while only three components of the saturated Al\ion{ii} line were consistent within $1\sigma$ (the other two were within $2\sigma$). We therefore cannot confirm potential changes in the column densities over the past two decades.

\begin{table*}
\small
\begin{center}
\caption{Fitted UVES Voigt profile parameters for low ions using the same five components' redshifts of the ESPRESSO fit from Table \ref{tab:loComp}}
\label{tab:NloFixEsp}
\begin{tabular}{lccccc}
\hline
Comp.& $z$& $\Delta v$ (km s$^{-1}$)& $b_{\rm turb}${} (km s$^{-1}$)& logN(Al\ion{ii})& logN(Fe\ion{ii})\\
\hline
1& $3.389512$& $-42$& $8.5\pm0.27$&  $12.12\pm0.017$& $13.39\pm0.017$\\
2& $3.389892$& $-16$& $4.3\pm0.12$&  $12.49\pm0.018$& $13.60\pm0.019$\\
3& $3.390056$& $-5$& $6.0\pm0.16$&  $12.97\pm0.035$& $14.42\pm0.022$\\
4& $3.390188$& $4$& $4.6\pm0.11$&  $12.89\pm0.036$& $14.19\pm0.021$\\
5& $3.390481$& $24$& $15.2\pm1.13$&  $12.05\pm0.033$& $12.81\pm0.086$\\
\end{tabular}
\end{center}
\end{table*}

\begin{table}
\small
\begin{center}
\caption{Fitted UVES Voigt profile parameters for C\ion{iv} using the same twelve components' redshifts of the ESPRESSO fit from Table \ref{tab:hiComp}}
\label{tab:NhiFixEsp}
\begin{tabular}{lcccc}
\hline
Comp.& $z$& $\Delta v$ (km s$^{-1}$)& $b_{\rm turb}${} (km s$^{-1}$)& logN(C\ion{iv})\\
\hline
1& $3.388100$& $-139$& $11.2\pm0.11$& $13.95\pm0.004$\\
2& $3.388645$& $-101$& $43.8\pm1.05$& $13.84\pm0.014$\\
3& $3.388729$& $-96$& $4.3\pm0.43$& $12.86\pm0.028$\\
4& $3.389040$& $-74$& $14.1\pm0.68$& $13.29\pm0.029$\\
5& $3.389703$& $-29$& $5.7\pm0.22$& $13.25\pm0.012$\\
6& $3.389918$& $-15$& $5.4\pm0.64$& $13.17\pm0.083$\\
7& $3.390172$& $3$& $11.5\pm1.28$& $13.90\pm0.046$\\
8& $3.390453$& $22$& $8.2\pm1.61$& $13.41\pm0.184$\\
9& $3.390717$& $40$& $16.2\pm1.84$& $14.06\pm0.044$\\
10& $3.390960$& $57$& $4.7\pm1.22$& $12.98\pm0.186$\\
11& $3.391181$& $72$& $10.3\pm0.69$& $13.59\pm0.037$\\
12& $3.391497$& $93$& $13.8\pm0.52$& $13.43\pm0.020$\\
\end{tabular}
\end{center}
\end{table}

\begin{table*}
\small
\begin{center}
\caption{Fitted UVES component column densities using fixed ESPRESSO $b_{\rm turb}${} and $z$ from Table \ref{tab:loComp}}
\label{tab:Nevolve}
\begin{tabular}{lccccc}
\hline
Ion & logN$_{\rm Comp.~1}$ & logN$_{\rm Comp.~2}$ & logN$_{\rm Comp.~3}$ & logN$_{\rm Comp.~4}$ & logN$_{\rm Comp.~5}$ \\
\hline
Al\ion{ii} & $12.15\pm0.013$ &$12.31\pm0.018$ &$13.18\pm0.068$ &$12.90\pm0.026$ &$12.06\pm0.018$ \\
Fe\ion{ii} & $13.37\pm0.018$ &$13.60\pm0.015$ &$14.59\pm0.024$ &$14.08\pm0.017$ &$12.99\pm0.060$ \\
\end{tabular}
\end{center}
\end{table*}

\section{Conclusions}

This paper presents an analysis of the performance of ESPRESSO's high spectral resolution, $4\times2$ binning mode (HR42) for observing faint quasars. Using four hours of on-target integration, we observed the quasar QSO J0003$-$2603 (AB magnitude $R = 18.5$~mag). An archival spectrum of the same quasar obtained with UVES for approximately the same exposure time allows one to compare ESPRESSO's performance at extracting properties of an intervening DLA absorber. Despite having $\approx2\times$ smaller S/N (per pixel) than the UVES spectra (Figure \ref{fig:snr}), the ESPRESSO HR42 mode observations were able to recover column densities and absorption profiles to nearly the same precision as the UVES data. For complex absorption profiles with many overlapping narrow features (such as the typical C\ion{iv} profiles seen in strong Ly$\alpha$ absorbers), fitting Voigt profiles becomes degenerate with the modest S/N of our ESPRESSO data. The low ESPRESSO throughput below $\lambda \approx5000$~\AA{} prevents accurate continuum fitting and modelling of lines. The higher resolving power of ESPRESSO ($R\approx136\,700$) allowed to identify additional narrow (i.e., small turbulent component of the Doppler broadening parameter) substructure within the metal absorption profiles that were missed by the typical UVES setup used for observing DLAs ($R\approx48~000$) despite the $\approx2\times$ lower S/N (per pixel). Whilst UVES can achieve a resolving power beyond 100~000, in order to match the same S/N depth  as ESPRESSO (after accounting for differing pixel size), observations with UVES would require double the exposure time. Furthermore, the spectral resolution element of UVES would not be properly sampled at this resolution. Thus ESPRESSO's HR42 mode provides an optimal combination of extended wavelength coverage, resolving power, and observing efficiency above 5000~\AA{} that UVES cannot obtain.

In summary, ESPRESSO's HR42 mode can provide additional stability and $\approx3\times$ the spectral resolving power with a comparable performance to UVES beyond $\lambda \geq 5000$ \AA{} for analysing quasar absorption line systems. As the equivalent width of a line scales as FWHM divided by S/N, the HR42 mode of ESPRESSO would be an ideal avenue for searching for weaker lines, such as C\ion{i}, Si\ion{i}, and rare elements, typically not easily identified with UVES. The increased resolution and sensitivity of the HR42 mode combined with the stability and precise wavelength calibration of ESPRESSO enables precision physics and cosmological experiments \citep{Irsic17, Cooke20, Schmidt21}, and observations quantifying isotopic ratios and gas temperatures of individual absorbers \citep{Narayanan06, Cooke15, Welsh20, Noterdaeme21} whilst minimizing observing time requirements using a single UT.

\begin{acknowledgements}
We are grateful for assistance, discussions and comments provided by Ryan Cooke, Valentina D'Odorico, Alain Smette, and Louise Welsh that significantly improved this manuscript, and for the comments from the anonymous referee that helped improve the clarity of the text. We are thankful for the ESPRESSO observations that were supported by and collected at the European Southern Observatory under ESO programme 60.A-9801(W).
\end{acknowledgements}

\bibliographystyle{aa}