\section{Introduction}

\noindent
Since the first experimental realization of {\it Bose-Einstein (BE) condensation} in the 90's, BE condensates and cold atoms in general have been extensively studied to investigate quantum properties on almost macroscopic scales. Among the typical features of BE condensates, one of the most striking is certainly  {\it superfluidity}, which has been studied in several experiments in the last years by putting the quantum system under rotation and observing its response (see, e.g., the reviews \cite{Co,Fe1}). Because of the quantum nature of BE condensates, the only possible change to the condensate profile due to the imposed rotation is the nucleation of isolated defects, i.e., quantum {\it vortices}. The generation of vortices has been observed in various experiments as well as the growth of their number when the angular velocity increases \cite{MCWD,RAVXK,CHES}. For even larger angular velocities the number of vortices becomes so large that they fill the bulk of the system and arrange in a typical Abrikosov lattice \cite{ARVK}. In presence of harmonic trapping the rotation can not be arbitrarily fast, otherwise the centrifugal forces would break down the trapping and system would eventually fly apart. On the opposite when the trapping contains some stronger confinement, e.g., some anharmonic potential growing faster that $ |\mathbf{r}|^2$ for large $ |\mathbf{r}| $, regimes with much more rapid rotation can in principle be reached. Unfortunately so far a loss of coherence of the system has prevented the exploration of such regimes in the experiments \cite{BSSD}, although a depression at the center of the trap has been observed for large angular velocities.

However it has been predicted \cite{CD,CDY2,Fe2,FJS,FB,KTU,KB,KF,R1} that, besides the nucleation of vortices, other phase transitions should be observed in rapid rotating condensates in case of anharmonic confinement, with the occurrence of macroscopic defects or the transition to {\it giant vortex} states: when the rotational velocity gets very large, the centrifugal forces constrain the condensate in some thin annular region around a macroscopic hole and, if the rotation gets even more rapid, vortices disappear from the bulk of the system, which seems then to carry a  huge circulation centered at the origin.

Although BE condensates are many-body quantum system composed of a number of atoms ranging from few thousands to many millions, all the physical prediction about them are made by using an effective theory, the {\it Gross-Pitaevskii (GP) theory}, namely a one-particle approximation in which the energy of the system is given by a suitable nonlinear functional (see below). In spite of its simplicity the agreement with experimental observations is quite good, specially in the so called Thomas-Fermi regime, i.e., when the effective coupling becomes large. One of the major advantages of GP theory is the possibility of run very sophisticated and accurate numerical simulations \cite{Dan,FJS,KTU}. See also the webpage \texttt{http://gpelab.math.cnrs.fr/}, where one can find an efficient free code for simulations of the GP energy or dynamics developed by \textsc{X. Antoine} and \textsc{R. Duboscq} \cite{AD1,AD2}.

In the framework of the GP theory the energy of a two-dimensional rotating BE condensate in physical units on the plane orthogonal to the rotational axis reads
\begin{equation}
\label{eq: gGPf}
	\mathcal{E}^{\mathrm{GP}}_{\mathrm{phys}}[\Psi] = \int_{\mathbb R^2} {\mathrm{d}} \mathbf{r} \: \left\{ \tx\frac{1}{2}\left|\left(\nabla - i \mathbf{A}_{\mathrm{rot}}\right) \Psi \right|^2 +
\left( V(r) - \tx\frac{1}{2} \Omega_{\rm rot}^2 r^2 \right) |\Psi|^2 + \frac{|\Psi|^4}{\varepsilon^2} \right\},
\end{equation}
where $ \Omega_{\mathrm{rot}} $ is the angular velocity, $ \mathbf{A}_{\mathrm{rot}} : = \Omega_{\mathrm{rot}} r \mathbf{e}_{\vartheta} $,  $r=|\mathbf r|$, with $\mathbf r = (x,y) \in\mathbb R^2$, and $ \mathbf{e}_{\vartheta} = (-y,x)/|\mathbf{r}| $ is the unit vector in the transverse direction. The trapping potential is assumed to be of the form
\begin{equation}
	\label{eq: potential}
	V(r) : = k  r^s + \tx\frac{1}{2} \Omega_{\mathrm{osc}}^2 r^2
\end{equation}
with $k>0$ and
\begin{equation}
	2 < s <\infty,
\end{equation}
i.e., the harmonic trapping is corrected by some anharmonic perturbation. Finally, we will focus on the study of the Thomas-Fermi (TF) regime $\varepsilon \to 0$. The ground state energy of the system is thus obtained by minimizing the functional \eqref{eq: gGPf} under the normalization constraint $ \left\| \Psi \right\|_2 = 1 $, which amounts to require conservation of the particle number. Any minimizer is called condensate wave function and its modulus square, i.e., the associated probability distribution, is what can be observed experimentally.

The range of validity of the GP description as well the derivation of the GP effective theory from the quantum mechanical description of a condensed Bose gas is an interesting topic on its own, which has been completely solved in the non-rotating case \cite{LSY1,LSSY}. In presence of rotation on the other hand \cite{LS} contains a derivation of the GP functional, which is however restricted to bounded angular velocities and therefore not directly applicable to the case under discussion (see also \cite{BCPY} for further results). However we will not investigate further such questions and take as a starting point the GP theory.

The mathematical physics literature contains now a large number of works studying the behavior of the GP minimization problem in different asymptotic regimes of the angular velocity. If we restrict the discussion to trapping potentials of the form \eqref{eq: potential}, three {\it phase transitions} have been identified (see \cite{CPRY3} for an extensive discussion or \cite{CPRY4} for a more concise exposition), corresponding to three critical values of the rotational velocity. Here we briefly sum up the most relevant features of the physics of rotating condensates in anharmonic traps:
\begin{itemize}
	\item for small angular velocities $ \Omega_{\mathrm{rot}} $, the rotation has no effect on the condensate wave function, i.e., the minimizer of $ \mathcal{E}^{\mathrm{GP}}_{\mathrm{phys}} $ coincides with the one in absence of rotation \cite{AJR};
	\item when the first critical speed $ \Omega_{\mathrm{c_1}} \propto \varepsilon^{\frac{4}{s+2}}|\log\varepsilon| $ is crossed, one observes the nucleation of quantum vortices, i.e., isolated zero of the condensate wave function \cite{CR1};
	\item if $ \Omega_{\mathrm{rot}} $ stays far from a second critical speed $ \Omega_{\mathrm{c_2}} \propto \varepsilon^{-\frac{s-2}{s+2}} $, the number of vortices might increase but the profile of the condensate wave function is still close to the non-rotating one. Close to $ \Omega_{\mathrm{c_1}} $ it is possible to derive the explicit distribution of vortices \cite{CR1}, which eventually cover the whole bulk of the condensate. In this regime one expects that they arrange in a regular (Abrikosov) lattice to minimize the interaction energy. This remains an open question although it has been proven that the vorticity is uniformly distributed \cite{CPRY3};
	\item for $ \Omega_{\mathrm{rot}} \propto \varepsilon^{-\frac{s-2}{s+2}} $ a first change of the macroscopic profile of the condensate is observed, due to the effect of the centrifugal forces. When the second critical speed $ \Omega_{\mathrm{c_2}} $ is crossed this change has a dramatic effect since a macroscopic hole is created at the center of the trap. However the vorticity remains uniform in the bulk of the system \cite{CPRY3};
	\item for very rapid rotations above $ \Omega_{\mathrm{c_2}} $, the bulk of the condensate becomes essentially annular and its width shrinks as $  \varepsilon \to 0 $. No further changes are however observed until a third critical speed $ \Omega_{\mathrm{c_3}} \propto \varepsilon^{-\frac{4(s-2)}{s+2}} $ is crossed. Then vortices are expelled from the bulk and the condensate behaves as if the whole vorticity was concentrated at the origin of the trap \cite{CPRY3,CPRY5}. This is the giant vortex state that we plan to study in this paper.
\end{itemize}

So far we have only discussed condensates in anharmonic traps of the type \eqref{eq: potential} but a lot of results are also available for other classes of trapping potentials. First of all the harmonic case has been extensively studied both in the physics and mathematical literature and, while there exists a first critical value of the angular velocity \cite{IM1,IM2} corresponding to the occurrence of vortices and the behavior of the condensate for not too rapid rotation is similar to the one described above (vortex lattice, uniform distribution of vorticity, etc.), when the angular velocity approaches the harmonic frequency of the trap, some new physical features come into play and fractional quantum Hall states emerge \cite{ABD,ABN,LSY2}. As we have already mentioned larger angular velocity are not allowed because the system would otherwise be no longer trapped. See however \cite{Ka} for an alternative setting in which the trapping potential is suitably rescaled to reach fast rotation regimes.

Even if we restrict to the anharmonic traps \eqref{eq: potential} there is an extreme case which is of certain interest, namely $ s = \infty $. Formally this corresponds to a confinement of the system to a two-dimensional disc of unit radius. Naturally one has then to provide suitable boundary conditions and both the Neumann \cite{CDY1,CY,CRY} and Dirichlet \cite{CPRY1} cases have been deeply studied. Indeed phase transitions analogous to the one described above has been found out even in this extreme case, although the nature of the third one is much more subtle.

Let us now go back to the functional \eqref{eq: gGPf} and introduce more convenient parameters: if we set
\begin{equation}
	\Omega_{\mathrm{phys}} : = \sqrt{\Omega_{\mathrm{rot}}^2 - \Omega_{\mathrm{osc}}^2},
\end{equation}
and obviously assume that $ \Omega_{\mathrm{osc}} < \Omega_{\mathrm{rot}} $, the trapping potential can be cast in the form
\begin{equation}
	\label{eq: eff potential}
	V(r) = k r^s - \textstyle\frac{1}{2} \Omega_{\mathrm{phys}}^2 r^2.
\end{equation}
Since we are interested in exploring a regime in which both $ \varepsilon \to 0 $ and $ \Omega_{\mathrm{phys}} \to \infty $ (or, equivalently,  $ \Omega_{\mathrm{rot}} \to \infty $), it is convenient to rescale units in the GP functional, in order to observe a non-trivial behavior \cite[Sect. I.A]{CPRY3}: if one would naively minimize $ \mathcal{E}^{\mathrm{GP}}_{\mathrm{phys}} $ under the mass constraint $ \| \Psi \| = 1$, one would get trivially that the ground state energy diverges and the corresponding minimizer tends to $ 0 $ pointwise. The appropriate rescaling depends however on the asymptotics of $  \Omega_{\mathrm{phys}} $ and, in the regime we want to explore (very fast rotation),
\begin{equation}
	\label{eq: lower bound speed}
	\Omega_{\mathrm{phys}} \gg \varepsilon^{-\frac{s-2}{s+2}},
\end{equation}
which leads to the rescaling (see \cite[Sect. I.A]{CPRY3})
\begin{equation}
 	\label{eq: rescaling}
	\mathbf{r} =  R_{\mathrm{m}} \mathbf{x},	\qquad	\Psi(\mathbf{r}) = R_{\mathrm{m}}^{-1} \psi(\mathbf{x}),	\qquad \Omega_{\mathrm{phys}} =  R_{\mathrm{m}}^{-2} \Omega,	\qquad	\mathbf{A}_{\Omega}  = \Omega x \mathbf{e}_{\vartheta},
\end{equation}
where $ R_{\mathrm{m}} $ stands for the unique minimum point of the potential \eqref{eq: eff potential}, i.e., explicitly
\begin{equation}
	\label{eq: rm}
	R_{\mathrm{m}} : = \left( \frac{\Omega_{\mathrm{phys}}^2}{s k} \right)^{\frac{1}{s-2}}.
\end{equation}
Under the scaling \eqref{eq: rescaling}, the GP functional \eqref{eq: gGPf} becomes
\begin{equation}
	\mathcal{E}^{\mathrm{GP}}_{\mathrm{phys}}[\Psi] = R_{\mathrm{m}}^{-2} \left[ \mathcal{E}^{\mathrm{GP}}[\psi]  + \left( \textstyle\frac{s}{2} - \tx\frac{1}{2} \right) \Omega^2 \right],
\end{equation}
with
\begin{equation}
	\label{eq: gpf}
	\framebox{$\mathcal{E}^{\mathrm{GP}}[\psi] : = \displaystyle\int_{\mathbb{R}^2} \mathrm{d} \mathbf{x} \left\{ \tx\frac{1}{2} \left| \left( \nabla - i \mathbf{A}_{\Omega} \right) \psi \right|^2 + \Omega^2 W(x) |\psi|^2 + \varepsilon^{-2} |\psi|^{4} \right\}$}.
\end{equation}
The rescaled potential
\begin{equation}
	\label{eq: W}
	W(x) : = \textstyle\frac{1}{s} \left( x^s - 1 \right) - \frac{1}{2} \left( x^2 - 1 \right)
\end{equation}
is positive and has a unique minimum at $ x = 1 $, i.e., $ \inf_{x \in \mathbb{R}^+} W(x) = W(1) = 0 $. The rescaled angular velocity $ \Omega $ is related to the original physical quantities via
\begin{equation}
	 \Omega = (sk)^{-\frac{2}{s-2}} \Omega_{\rm phys}^{\frac{s+2}{s-2}},
\end{equation}
and condition \eqref{eq: lower bound speed} becomes
\begin{equation}
	\Omega \gg \varepsilon^{-1}.
\end{equation}

From now on we will focus on the analysis of the minimization of the functional \eqref{eq: gpf} on the domain
\begin{equation}
	\label{eq: domain}
	{\mathscr{D}^{\mathrm{GP}} : = \left\{ \psi \in H^1(\mathbb{R}^2)\: \big| x^{s/2} \psi \in L^2(\mathbb{R}^2), \left\| \psi \right\|_2 = 1 \right\}}.
\end{equation}
We also set
\begin{equation}
	\label{eq: ground state}
	E^{\mathrm{GP}} : = \inf_{\psi \in \mathscr{D}^{\mathrm{GP}}} \mathcal{E}^{\mathrm{GP}}[\psi],
\end{equation}
and denote by $ \psi^{\mathrm{GP}} $ any minimizer, which is known to exist by standard arguments. In addition any $ \psi^{\mathrm{GP}} $, which might be non-unique due to a breaking of the rotational symmetry and the occurrence of isolated vortices, solves the variational equation
\begin{equation}
	\label{eq: GP variational}
	- \tx\frac{1}{2} \left(\nabla - i \mathbf{A}_{\Omega} \right)^2 \psi^{\mathrm{GP}}  + \Omega^2 W(x) \psi^{\mathrm{GP}} + 2 \varepsilon^{-2} \left| \psi^{\mathrm{GP}} \right|^2 \psi^{\mathrm{GP}} = \mu^{\mathrm{GP}} \psi^{\mathrm{GP}},
\end{equation}
where the chemical potential (Lagrange multiplier) is fixed by imposing the $L^2-$normalization of $ \psi^{\mathrm{GP}} $:
\begin{equation}
	\label{eq: chem}
	\mu^{\mathrm{GP}} = E^{\mathrm{GP}} + \varepsilon^{-2} \int_{\mathbb{R}^2} \mathrm{d} \mathbf{x} \: \left| \psi^{\mathrm{GP}} \right|^4.
\end{equation}

As discussed in details in \cite[Sect. I.B]{CPRY3}, when $ \Omega \gg \varepsilon^{-1} $ the condensate has already crossed the second critical speed, i.e., its profile approaches a density supported on an annulus centered in the origin, whose inner and outer radii tend to $ 1 $ as $ \varepsilon \to 0 $. More precisely $ |\psi^{\mathrm{GP}}|^2 $ is close {in $ L^p $, $ p < \infty $,} to the TF profile
\begin{equation}
	\label{eq: tfm}
	{\rho^{\mathrm{TF}}}(x) = \textstyle\frac{1}{2} \left[ \mu^{\mathrm{TF}} - \varepsilon^2 \Omega^2 W(x) \right]_+,
\end{equation}
with $  \mu^{\mathrm{TF}} $ the chemical potential fixed by the $ L^1$-normalization of the function. A straightforward analysis shows indeed that $ {\rho^{\mathrm{TF}}} $ is compactly supported and $ \mathrm{supp}({\rho^{\mathrm{TF}}}) = \left[ x_{\mathrm{in}}, x_{\mathrm{out}} \right] $ with \cite[Eq. (2.7)]{CPRY3}
\begin{equation}
	x_{\mathrm{out}} - x_{\mathrm{in}} = C \left(\varepsilon \Omega\right)^{-2/3} {\ll 1},		\qquad		{x_{\mathrm{in}/\mathrm{out}} = 1 + \mathcal{O}\big(\left(\varepsilon \Omega\right)^{-2/3}\big)},
\end{equation}
{as it can be proven by taking a Taylor expansion of $ W $ around $ x = 1 $ in \eqref{eq: tfm} and imposing the $ L^1 $ normalization.}

The vortex structure of $ \psi^{\mathrm{GP}} $ is richer: being well above the first critical speed $ \Omega_{\mathrm{c_1}} \sim |\log\varepsilon| $ for the nucleation of vortices, the GP minimizer contains a very large number of vortices distributed all over its support. More precisely one can prove that the vorticity is uniformly distributed in the bulk of the condensate. {As in \cite[Eq. (1.42)]{CPRY3}, we denote by $ {\mathcal R}_{\rm bulk} \subset \mathrm{supp}({\rho^{\mathrm{TF}}}) $ a suitable annulus $ \left\{ \mathbf{x} \: | \: x_< \leq x \leq x_> \right\} $ with $ x_{>/<} = x_{\mathrm{out}/\mathrm{in}} + o((\varepsilon \Omega)^{-2/3}) $}.

	\begin{teo}[\textbf{\cite[Theorem 1.2]{CPRY3}}]
		\label{teo: vortex distribution CPRY3}
		\mbox{}	\\
		If $ \varepsilon^{-1} \lesssim \Omega \ll \varepsilon^{-4} $ as $ \varepsilon \to 0 $, there exists a finite family of disjoint balls $ \{ \mathcal{B}_i \} : = \{ \mathcal{B}(\mathbf{x}_i, \varrho_i) \} \subset {\mathcal R}_{\rm bulk} $, $ i = 1, \ldots, N $, such that 		\begin{enumerate}
			\item
$ \varrho_i \leq \mathcal{O}\left(\Omega^{-1/2}\right) $,
$ \sum \varrho_i^2 \leq (1 + (\varepsilon\Omega)^{2/3})^{-1} $;
			\item $ \left|\psi^{\mathrm{GP}}\right| > 0 $ on $ \partial \mathcal{B}_i $, $ i = 1, \ldots, N $.
		\end{enumerate}
		Moreover, setting $ d_i : = \deg\{\psi^{\mathrm{GP}}, \partial \mathcal{B}_i\} $ and defining the vorticity measure as  $ \mu : =  2\pi\sum_{i = 1}^N d_i \delta(\mathbf{x}- \mathbf{x}_i) $, then, for any set $ \mathcal{S} \subset {\mathcal R}_{\rm bulk} $ such that {$ |\partial \mathcal{S}| = 0 $ and} $ |\mathcal{S}| \gg \Omega^{-1} |\log(\varepsilon^4\Omega)|^2 $ as $ \varepsilon \to 0 $,
		\begin{equation}
			\frac{\mu(\mathcal{S})}{\Omega |\mathcal{S}|} \underset{\varepsilon \to 0}{\longrightarrow} 1.
		\end{equation}
	\end{teo}

	\noindent
	The inner region $ \left\{ \mathbf{x} \in \mathbb{R}^2 \: | \: x \leq x_{\mathrm{in}} \right\} $ is presumably also filled with vortices but, because of the exponential smallness of $ \psi^{\mathrm{GP}} $ there, the vortex structure in that region is practically inaccessible. An important condition contained in Theorem \ref{teo: vortex distribution CPRY3} is the request
	\begin{displaymath}
		\Omega \ll \varepsilon^{-4}.
	\end{displaymath}
	The reason is that at angular velocities of that order the proof of Theorem \ref{teo: vortex distribution CPRY3} might fail due to the occurrence of a further phase transition, i.e., the transition to a giant vortex state. This paper is precisely devoted to the investigation of such a transition.

	From the heuristic point of view it is quite simple to explain why one should expect a change in the vortex structure when $ \Omega \sim \varepsilon^{-4} $: from energy considerations it is easy to see that the average size of the vortex core, i.e., the radius of the region around a vortex when $ |\psi^{\mathrm{GP}}|^2 $ is substantially far from $ {\rho^{\mathrm{TF}}} $, is of order $ \varepsilon^{2/3} \Omega^{-1/3} $. The width of $ \mathrm{supp}({\rho^{\mathrm{TF}}}) $ is on the other hand of order $ (\varepsilon\Omega)^{-2/3} $ and the two quantities are clearly of the same order when $ \Omega \sim \varepsilon^{-4} $. Hence it must happen that for $ \Omega = \Omega_0 \varepsilon^{-4} $ with $ \Omega_0 $ a large enough constant, the vortex core becomes larger that the bulk of the condensate, i.e., vortices can no longer be accommodated in $ \mathrm{supp}({\rho^{\mathrm{TF}}}) $. A non trivial phase factor of $ \psi^{\mathrm{GP}} $ is however needed in order to compensate the effect of the rotation but, because no vortex can occur in the bulk of the condensate, all the vorticity should get concentrated in the inner region where $ \psi^{\mathrm{GP}} $ is exponentially small. In fact when this occurs it is impossible to distinguish from the energetic point of view such a state with vortices distributed in the inner hole from a giant vortex state of the form $ f(x) e^{i n \vartheta} $, $ n \in \mathbb{Z} $.

	Notice that although this might seem to suggest that the rotational symmetry is restored, such a phenomenon never occurs as proven in \cite[Theorem 1.6]{CPRY3}. However the GP energy is expected to be well approximated above the critical speed for the transition to a giant vortex state by a one-dimensional energy functional obtained by evaluating $ \mathcal{E}^{\mathrm{GP}} $ on functions of the form $ f(x) e^{i n \vartheta} $. In fact by some very simple observations one can show that $ n = \lfloor \Omega \rfloor (1 + o(1)) $, where $ \lfloor \: \cdot \: \rfloor $ stands for the integer part. Let us now fix the angular velocity to be
	\begin{equation}
		\label{eq: speed}
		\Omega = \frac{\Omega_0}{\varepsilon^4},
	\end{equation}
	with $ \Omega_0 $ a positive constant. Concerning the giant vortex regime, the main results proven in \cite{CPRY3} are stated below. {We denote by $ \mathcal{A}_{\mathrm{bulk}} $ a suitable annular layer around $ x = 1 $ containing the bulk of the condensate (see next \eqref{eq: bulk} for a precise definition).}

	\begin{teo}[\textbf{\cite[Theorem 1.3]{CPRY3}}]
		\label{teo: giant vortex CPRY3}
		\mbox{}	\\
		If $ \Omega $ is given by \eqref{eq: speed}, there exists a finite constant $ \bar\Omega_0 $ such that for any $ \Omega_0 > \bar\Omega_0 $, no minimizer $ \psi^{\mathrm{GP}} $ has a zero inside $ \mathcal{A}_{\mathrm{bulk}} $ if $ \varepsilon $ is sufficiently small.
	\end{teo}

	\begin{teo}[\textbf{\cite[Theorem 1.4]{CPRY3}}]
		\label{teo: gv energy CPRY3}
		\mbox{}	\\
		If $ \Omega $ is given by \eqref{eq: speed} with $ \Omega_0 > \bar\Omega_0 $  as in Theorem \ref{teo: giant vortex CPRY3}, then as $ \varepsilon \to 0 $\footnote{{We use here polar coordinates $ (x,\vartheta) \in \mathbb{R}^+ \times [0,\pi) $ on the plane.}}
		\begin{equation}
			\label{teo: gpe asympt CPRY3}
			E^{\mathrm{GP}} = \min_{\left\| f \right\|_2 = 1} \mathcal{E}^{\mathrm{GP}}\left[f(x) e^{i \lfloor \Omega \rfloor \vartheta}\right] + \mathcal{O}(|\log \varepsilon| ^{9/2}).
		\end{equation}
	\end{teo}

	\noindent
	The first result, although being a consequence of the energy asymptotics \eqref{teo: gpe asympt CPRY3}, is the most relevant one, since it shows the occurrence of the giant vortex transition for angular velocities of order $ \varepsilon^{-4} $. The precise mathematical statement is a pointwise estimate in the bulk region of $ |\psi^{\mathrm{GP}}| $ in terms of a strictly positive function, i.e., the minimizer of the functional appearing on the r.h.s. of \eqref{teo: gpe asympt CPRY3}: since the latter is bounded from below by a positive constant in the bulk and the difference is pointwise small in $ \varepsilon $, also $ \psi^{\mathrm{GP}} $ can not vanish there.

	For the analysis of the present paper it is very important to remark that both results hold true if the angular velocity is expressed by \eqref{eq: speed} with $ \Omega_0 $ {\it large enough}, namely no precise estimate is {derived there} on the sharp value for the transition {(see also Remark \ref{rem: comparison})}. We indeed expect that the giant vortex structure appears as soon as $ \Omega $ becomes (asymptotically) larger than
	\begin{equation}
		\label{eq: othird}
		\Omega_{\mathrm{c_3}} = \frac{\Omega_{\mathrm{c}}}{\varepsilon^4},
	\end{equation}
	for some explicit value $ \Omega_{\mathrm{c}} $. In this paper we will indeed investigate such a question and exhibit a finite value $ \Omega_{\mathrm{c}} $ which is a good candidate for the sharp constant. Actually we are going to see that such a constant is a solution of some algebraic equation (see \eqref{eq: Omegac}) involving quantities relative to a limit problem independent of $ \varepsilon $. Although we have not proven it yet, we do expect that next \eqref{eq: Omegac} has a unique solution, thus providing the sharp value of the critical velocity\footnote{Strictly speaking in order to show that $ \Omega_{\mathrm{c}} $ is the sharp value for the transition one should also prove that, below $ \Omega_{\mathrm{c}} $, vortices are still present in the bulk of the condensate, as done in \cite{R2} for hard anharmonic traps. We will come back to this question later.}.

	We outline here the structure of the paper. Next Section contains the main results, i.e., the identification of the explicit value of the angular velocity for the transition to the giant vortex state, together with an asymptotic expansion of the GP ground state energy which is actually on the main ingredients of the proof of the above mentioned result. We also show that, as in \cite[Theorem 1.5]{CPRY3}, one can deduce a (better) estimate of the total winding number of any GP minimizer.

	Sections \ref{sec: preliminary} contains some preliminary estimates and a detailed analysis of the effective functionals that will play a significant role throughout the proofs. In Section \ref{sec: positivity} we prove the main properties of the cost function and in particular its positivity, which is the main mathematical tool used in the proof of the giant vortex transition as in several other works \cite{CPRY1, CPRY2, CPRY3, CR1, CRY}.

	Sections \ref{sec: energy} and \ref{sec: gv transition} are devoted to the proofs of the main results: we first (Section \ref{sec: energy}) obtain the asymptotic expansion of the GP energy by comparing suitable upper and lower bounds and then (Section \ref{sec: gv transition}) use such a result to deduce the pointwise estimate of $ |\psi^{\mathrm{GP}}| $ showing the absence of vortices in the bulk.

\medskip
\noindent
{{\bf Notation:} In the asymptotic analysis $ \varepsilon \to 0  $ we will often use the Landau symbols: given a positive function $ g $, we say that
	\begin{itemize}
		\item 	$ f = \mathcal{O}(g) $ (resp. $ = o(g) $), if $ \lim_{\varepsilon \to 0} |f|/g \leq C < \infty $ (resp. $ = 0 $);
		\item 	$ f \propto g $, whenever $ \lim_{\varepsilon \to 0} |f|/g = C $, with $ 0 < C < \infty $;
		\item 	if $ f \geq 0 $, $ f \ll g $ is synonimous of $ f = o(g) $ and $ f \gg g $ simply means that $ g \ll f $.
	\end{itemize}
	Sometimes we will use the notation $ \mathcal{O}(|\log\varepsilon|^{\infty}) $ to indicate a quantity of order $ |\log\varepsilon|^a $ for some finite but possibly large $ a $. Since such a quantity will typically appear multiplied by powers of $ \varepsilon $, the explicit value of $ a $ will be irrelevant. 	\\
	We denote by $ \mathcal{B}_{\varrho}(\mathbf{x}) $ any two-dimensional ball centered in $ \mathbf{x} $ and with radius $ \varrho $ and by $ \lfloor x \rfloor $ the integer part of the real number $ x $. The symbol $ C $ will stand for a finite constant independent of $ \varepsilon $, whose value might change from line to line.
}

\medskip
\noindent
{\bf Acknowledgements:} The authors acknowledge the support of MIUR through the FIR grant 2013 ``Condensed Matter in Mathematical Physics (COND-MATH)'' (code RBFR13WAET).

\section{Main Results}
\label{sec: main results}

The first non trivial observation to improve the results proven in \cite{CPRY3} is that instead of making a special choice of the giant vortex winding number ($ \lfloor \Omega \rfloor $ in \cite{CPRY3}), one might try and optimize w.r.t. such a parameter, so obtaining a better candidate for the giant vortex state. This leads to consider the functional obtained evaluating $ \mathcal{E}^{\mathrm{GP}} $ on a giant vortex ansatz $ f(x) e^{i n \vartheta}$ and minimize w.r.t. both $ f $ and $ n $ to find out the optimal giant vortex phase, i.e., explicitly
\begin{equation}
	\label{eq: gvbf}
	\mathcal{E}_\beta^\mathrm{gv}[g] = \int _{-\eta}^{\eta} \mathrm{d} y \: (1 + \varepsilon^2 y) \left\{ \tx\frac{1}{2}
		\left|
			\nabla g(y)
		\right| ^2
		+U _\beta (y) g ^2 (y)
		+\varepsilon ^2 y ^3 v(y) g ^2(y)
		+\frac{1}{2\pi}
		g ^4(y)
	\right\},
\end{equation}
where we have set for convenience $ n = \Omega + \beta $ and exploited the exponential fall off of $ \psi^{\mathrm{GP}} $ to cut the tails $ |y| \geq \eta \propto |\log\varepsilon| $. The spatial coordinate has also been rescaled around $ |\mathbf{x}| = 1 $ by setting $ x = 1 + \varepsilon^2 y $. The potentials $ U_\beta $ and $ v $ are obtained via a Taylor expansion of $ W(x) $ around $ x = 1 $ and to the leading order in $ \varepsilon $ are simply given by a shifted quadratic potential (see \eqref{eq: pot U} and \eqref{eq: pot v} for their explicit expressions). We remark however that in $ U_\beta $, the parameter $ \beta $ always appears multiplied (at least) by $ \varepsilon^2 $, so showing that the correction is only lower order.

Setting $ E_\beta^\mathrm{gv} : = \inf_{\| f \| = 1} \mathcal{E}_\beta^\mathrm{gv}[f] $ and denoting by $ g_\beta $ the corresponding minimizer, which can be proven to exist and be unique (up to multiplication by a phase factor) (see Proposition \ref{pro: min gvbf}), one can subsequently minimize w.r.t. $ \beta \in \mathbb{R} $, obtaining the energy $ E_{\star}^\mathrm{gv} $, an optimal phase $ {\beta_{\star}} $ and a density $ g_{\star} $, i.e.,
\begin{equation}
	\label{eq: gcre}
	E_{\star}^\mathrm{gv} := \min_{\beta \in \mathbb{R}} E_\beta^\mathrm{gv} = E^{\mathrm{gv}}_{{\beta_{\star}}} = \mathcal{E}^{\mathrm{gv}}_{{\beta_{\star}}}[g_{\star}].
\end{equation}

In Subsection \ref{sec: betas} we will prove that $ {\beta_{\star}} = \mathcal{O}(1) $, so that, by the above argument, one expects the functional $ \mathcal{E}_\beta^\mathrm{gv} $ to be close in the limit $ \varepsilon \to 0 $ to the following simplified giant vortex functional
\begin{equation}
	\label{eq: gvf}
	\mathcal{E}^{\mathrm{gv}}[g] = \int _{\mathbb{R}} \mathrm{d} y
	\left\{
		\textstyle\frac{1}{2}\left(
			g '
		\right) ^2
		+ \frac{\alpha ^2}{2} y^2 g ^2
		+\frac{1}{2\pi} g ^4
	\right\},
\end{equation}
with ground state energy $ E^{\mathrm{gv}} $ and minimizer $ g_{\mathrm{gv}} $, i.e.,
\begin{equation}
	E^{\mathrm{gv}} : = \inf_{g \in \mathscr{D}^{\mathrm{gv}}} \mathcal{E}^{\mathrm{gv}}[g] = \mathcal{E}^{\mathrm{gv}}[g_{\mathrm{gv}}],
\end{equation}
where
\begin{displaymath}
	\mathscr{D}^{\mathrm{gv}} : = \left\{ g \in H^1(\mathbb{R}) \: \big| \: y g \in L^2(\mathbb{R}), \left\| g \right\|_{L^2(\mathbb{R})}  = 1 \right\}.
\end{displaymath}
Here we have denoted for short
\begin{equation}
	\label{eq: alpha}
	\alpha : = \Omega_0 \sqrt{s+2}.
\end{equation}
The minimizer $ g_{\mathrm{gv}} $ solves the variational equation
\begin{equation}
	- \textstyle\frac{1}{2} g^{\prime\prime} + \frac{1}{2} \alpha^2 y^2 g + \frac{1}{\pi} g^3 = \mu^{\mathrm{gv}} g,
\end{equation}
where $ \mu^{\mathrm{gv}} = E^{\mathrm{gv}} + \frac{1}{2\pi} \left\| g_{\mathrm{gv}} \right\|_4^4 $.

We are now in position to introduce the explicit value of the constant $ \Omega_{\mathrm{c}} $ appearing in the critical value of the angular velocity $ \Omega_{\mathrm{c_3}} $, which can be expressed in terms of the critical quantities associated with the effective one-dimensional functional $ \mathcal{E}^{\mathrm{gv}} $ and, specifically, $ g_{\mathrm{gv}} $ and $ \mu^{\mathrm{gv}} $: we denote by $ \Omega_{\mathrm{c}} $ the {\it largest} solution of the equation
\begin{equation}
	\label{eq: Omegac}
 	\framebox{$\Omega_0 = \displaystyle\frac{4}{s+2} \left[ \mu^{\mathrm{gv}} - \frac{1}{2\pi} g_{\mathrm{gv}}^2(0) \right],$}
\end{equation}
where the r.h.s. depends on $ \Omega_0 $ through $ \mu^{\mathrm{gv}} $ and $ g_{\mathrm{gv}} $. The existence of such a solution is proven in Proposition \ref{pro: Omegac}.
Note that thanks to the estimate $ \left\| g_{\mathrm{gv}} \right\|^2_{\infty} \leq  \pi \mu^{\mathrm{gv}} $ (see \eqref{eq: gvm infty bound}), $ \Omega_{\mathrm{c}} > 0 $.

Before stating the main result of this paper, we have to define more precisely the region we identify with the bulk of the condensate: we set for any $ a > 0 $
\begin{equation}
	\label{eq: bulk}
	\mathcal{A}_{\rm bulk} : = \left\{ \mathbf{x} \in \mathbb{R}^2 \: | \: g_{\mathrm{gv}}\left(\textstyle\frac{x-1}{\varepsilon^2}\right) \geq |\log\varepsilon|^{-a} \right\},
\end{equation}
and observe that by the exponential decay proven in Proposition \ref{pro: gpm exp small}, $ \left\| \psi^{\mathrm{GP}} \right\|_{L^2(\mathcal{A}_{\rm bulk})} = 1 + o(1) $, i.e., it certainly contains the bulk of the system.

\begin{teo}[Absence of vortices in $ \mathcal{A}_{\rm bulk} $]
	\label{teo: no vortices}
	\mbox{}	\\
	If $ \Omega = \Omega_0/ \varepsilon^{4} $ with $ \Omega_0 > \Omega_{\mathrm{c}} $ as $ \varepsilon \to 0 $, then no GP minimizer $ \psi^{\mathrm{GP}} $ contains vortices in $ \mathcal{A}_{\rm bulk} $. More precisel
	for any $ \mathbf{x} \in \mathcal{A}_{\rm bulk} $
	\begin{equation}
		\label{eq: no vortices}
		\framebox{$\left| \psi^{\mathrm{GP}}(\mathbf{x}) \right|  = \displaystyle\frac{1}{\sqrt{2\pi} \varepsilon} \: g_{\mathrm{gv}}\left(\textstyle\frac{x-1}{\varepsilon^2}\right)  \left(1 + \mathcal{O}(\varepsilon^{1/2} |\log\varepsilon|^{\infty}) \right)$.}
	\end{equation}
\end{teo}

\begin{rem}[Giant vortex structure]
	\mbox{}	\\
	The pointwise estimate \eqref{eq: no vortices} suggests that $ |\psi^{\mathrm{GP}}| $ is approximately radial within $ \mathcal{A}_{\rm bulk} $. As already mentioned, this does not imply that the rotational symmetry is restored, since one expects that $ |\psi^{\mathrm{GP}}| $ is far from being radial in the inner region $ x \leq x_{\mathrm{in}} $, where several vortices should presumably be distributed more or less uniformly. In any case no GP minimizer is invariant under rotations if $ \Omega $ is large enough \cite[Theorem 1.6]{CPRY3}, e.g., in the giant vortex regime.
\end{rem}

\begin{rem}[Third critical velocity]
	\label{rem: sharp}
	\mbox{}	\\
	In Proposition \ref{pro: Omegac} we will prove that the equation \eqref{eq: Omegac} has a solution. Although we do not prove it, we strongly believe that such a solution is in fact unique and identifies the sharp constant in the value of the third critical speed.
	\newline
	More precisely Theorem \ref{teo: no vortices} indicates that above $ \Omega_{\mathrm{c}}/\varepsilon^{4} $ the system undergoes the phase transition to the giant vortex state and the bulk of the condensate becomes vortex free. Hence
	\begin{equation}
		\Omega_{\mathrm{c_3}} \leq \frac{\Omega_{\mathrm{c}}}{\varepsilon^4}.
	\end{equation}
	We actually expect that
	$
		\Omega_{\mathrm{c_3}} =\Omega_{\mathrm{c}}/\varepsilon^4
	$,
	which obviously requires to prove that the solution to \eqref{eq: Omegac} is unique. In addition one should also prove that for slower rotations vortices are still present in the bulk of the system. We plan to attack such a problem in a future work, but here we want to stress that the negativity of the cost function (see next Section \ref{sec: heuristics}) for $ \Omega_0 < \Omega_{\mathrm{c}} $ is a very strong indication that vortices are indeed convenient in this case and thus the sharp value of the critical speed is precisely $ \Omega_{\mathrm{c}}/\varepsilon^4 $.
\end{rem}

{\begin{rem}[Comparison with \cite{CPRY3}]
	\label{rem: comparison}
	\mbox{}	\\
	We want here to discuss in more details the comparison between Theorem \ref{teo: no vortices} and the analogous result proven in \cite[Theorem 1.3]{CPRY3}: in principle, one could indeed derive an estimate of the threshold $ \bar{\Omega}_0 $ for the transition to the giant vortex state there and then it would be natural to compare it with the explicit value found here. However we provide here some heuristic arguments showing that such a comparison is actually not needed (see however next Remark \ref{rem: cost functions} for futher details).	\\
	First of all an explicit estimate of $ \bar{\Omega}_0 $ is not an easy task to achieve, due to the proof structure in \cite{CPRY3}: the result proven there is indeed obtained through an asyptotic analysis as $ \Omega_0 \to \infty $ and one should then estimate all the coefficients of the error terms appearing in the formulae. Such quantities ultimately depends on the pointwise estimate of the difference between the giant vortex profile and the ground state of the harmonic oscillator given in \cite[Proposition 3.5]{CPRY3}, which is not explicit at all. 	\\
	However, even assuming that one could obtain a sharp value $ \bar{\Omega}_0 $, there are strong reasons to believe that, unlike $ \Omega_{\mathrm{c}} $ (see also the previous Remark \ref{rem: sharp}), it can not be the coefficient of the critical speed. First of all the condition $ \Omega_0 > \bar{\Omega}_0 $ guarantees the positivity of the vortex energy cost in \cite{CPRY3} (Remark \ref{rem: sharp}) and therefore $ \bar{\Omega}_0 > \Omega_{\mathrm{c}} $. Moreover, as explained in \cite{CPRY3} (see also \cite{CPRY2}), when $ \Omega_0 \to \infty $, another transition takes place, i.e., the condensate density profile goes from a TF-like shape \eqref{eq: tfm} to a gaussian function minimizing some suitable harmonic energy. The key fact is that such a transition is expected to take place \underline{after} the giant vortex one. Indeed here we show that, for finite $ \Omega_0 $, when the profile change has not yet occurred, the condensate is already in a giant vortex state. On the opposite, a quick inspection to the proof in \cite{CPRY3} reveals that the transition to the giant vortex is proven there by imposing that the profile is already gaussian. Hence any so obtained threshold value can not be meaningful.
\end{rem}}

\begin{rem}[Giant vortex density]
	\mbox{}	\\
	We have formulated the pointwise estimate \eqref{eq: no vortices} with $ g_{\mathrm{gv}} $, but an analogous statement holds true with $ g_{\mathrm{gv}} $ replaced with $ g_{\star} $. The error in \eqref{eq: no vortices} is indeed so large that one can not appreciate the difference between the two reference profiles (see Proposition \ref{pro: infty est diff gvm}). Let us stress however that the use of $ g_{\star} $ as a reference profile in the proof is on the opposite crucial to obtain the result (compare, e.g., the asymptotics \eqref{eq: energy asymptotics} and \eqref{eq: energy asympt gv}).
\end{rem}

	\noindent
	The absence of vortices proven in Theorem \ref{teo: no vortices} and the pointwise estimate of $  \psi^{\mathrm{GP}}$ follows from a refined result about the energy asymptotics in the same regime, that we state in the following

\begin{teo}[Energy asymptotics]
	\label{teo: energy asymptotics}
	\mbox{}	\\
	If $ \Omega = \Omega_0 \varepsilon^{-4} $ with $ \Omega_0 > \Omega_{\mathrm{c}} $ as $ \varepsilon \to  0 $, then
	\begin{equation}
		\label{eq: energy asymptotics}
		\framebox{$E^{\mathrm{GP}} = \displaystyle\frac{E^{\mathrm{gv}}_{\star}}{\varepsilon^4} + \mathcal{O}(1).$}
	\end{equation}
\end{teo}

\begin{rem}[Energy expansion]
	\mbox{}	\\
	The leading term $ E^{\mathrm{gv}}_{\star} /\varepsilon^{4} $ contains the main energy contribution due to the inhomogeneity of the GP profile together with the subleading kinetic energy of $ |\psi^{\mathrm{GP}}| $. The absence of vortices in $ \mathcal{A}_{\rm bulk} $ can be read in the very small remainder term $ \mathcal{O}(1) $. It is indeed interesting to compare \eqref{eq: energy asymptotics} with the analogous result \cite[Theorem 1.4]{CPRY3}, where the error term is much larger, i.e., $ \mathcal{O}(|\log\varepsilon|^{9/2}) $, in addition to the fact that the result proven there holds true only for $ \Omega_0 $ large enough.

	Notice however that the coefficient of the leading term $ E^{\mathrm{gv}}_{\star} $ still depends on $ \varepsilon $, through the boundaries of the integration domain as well as the optimal phase $ {\beta_{\star}} $ and the potential $ U_{{\beta_{\star}}} $. If one wanted to extract a proper asymptotic expansion then the natural statement would be
	\begin{equation}
		\label{eq: energy asympt gv}
		E^{\mathrm{GP}} = \frac{E^{\mathrm{gv}}}{\varepsilon^4} + \mathcal{O}(|\log\varepsilon|^7),
	\end{equation}
	with a much worse error term.
\end{rem}

	\noindent
	Thanks to the pointwise statement \eqref{eq: no vortices}, one can deduce that $ \psi^{\mathrm{GP}} $ does not vanish on $ \mathcal{A}_{\rm bulk} $. In particular for any $ R = 1  + \mathcal{O}(\varepsilon^2) $, $ |\psi^{\mathrm{GP}}| > 0 $ on $ \partial \mathcal{B}_R $. Hence it is possible to define the winding number of $ \psi^{\mathrm{GP}} $ on $ \partial \mathcal{B}_R $ for any such $ R $. A consequence of the energy asymptotics and the estimate \eqref{eq: no vortices} is thus the following

\begin{teo}[Winding number]
	\label{teo: winding number}
	\mbox{}	\\
	Let $ \Omega = \Omega_0 \varepsilon^{-4} $ with $ \Omega_0 > \Omega_{\mathrm{c}} $ and $ R $ be any radius such that $ R = 1 + \mathcal{O}(\varepsilon^2) $ as $ \varepsilon \to  0 $, then
	\begin{equation}
		\framebox{$ \deg \left(\psi^{\mathrm{GP}}, \partial \mathcal{B}_R\right) = \displaystyle\frac{\Omega_0}{\varepsilon^4} + \mathcal{O}(1). $}
	\end{equation}
\end{teo}

\noindent
{Note that the combination of the above result with the proof of the rotational symmetry breaking given in \cite[Theorem 1.6]{CPRY3} implies the presence of vortices in the inner hole region where $ \psi^{\mathrm{GP}} $ is exponentially small.}

\subsection{Heuristics}
\label{sec: heuristics}

Before discussing the proofs of the main results, we briefly expose the proof strategy from a heuristic point of view, i.e., not tracking down the error terms and neglecting most technical points. As usual the main result about the behavior of the condensate wave function is deduced from the energy asymptotics \eqref{eq: energy asymptotics}. We thus focus on such a proof.

Most of the relevant features of a fast rotating Bose-Einstein condensate were already discussed in details in \cite{CPRY3} and recalled in the Introduction. Here we take as a starting point the effective functional \eqref{eq: gvbf} which is expected to provide the leading order term in the energy asymptotics in units $ \varepsilon^{-4} $. Note that the ground state energy of $ \mathcal{E}_\beta^\mathrm{gv} $ always provides an upper bound to $ E^{\mathrm{GP}} $ for any {\it integer} phase, i.e., whenever $ \Omega + \beta \in \mathbb{Z} $. Actually the same upper bound can be proven to hold true up to some small error term even if $ \Omega + \beta $ is not an integer (see Section \ref{sec: energy ub}). Hence we can neglect the upper bound part of the proof and discuss only the lower estimate to $ E^{\mathrm{GP}} $.

A preliminary step which is already described in details in \cite{CPRY3} is the restriction of the integration in $ \mathcal{E}^{\mathrm{GP}} $ to the bulk of the condensate, i.e., to an annulus centered in the origin with radius $ \simeq 1 $ and width $ \mathcal{O}(\varepsilon^2) $. This can be done by exploiting the exponential decay of $ \psi^{\mathrm{GP}} $ outside. From now on we will then assume that the integration in $ \mathbf{x} $ is restricted to the annulus $ \left| 1 - x \right| \leq \mathcal{O}(\varepsilon^2|\log\varepsilon|) $.

The main steps in the energy lower bound are then the following:
\begin{enumerate}
	\item \underline{optimal giant vortex phase and profile}: we minimize $ E_\beta^\mathrm{gv} $ w.r.t. to $ \beta \in \mathbb{R} $ and obtain a minimizing $ {\beta_{\star}} $ and an associated density $ g_{\star} $. It is crucial to observe that such a minimization yields an additional equation involving $ g_{\star} $, which is in fact nothing but the vanishing of the first derivative of $ E_\beta^\mathrm{gv} $ w.r.t. $ \beta $. Such an equation will play a crucial role at point 4 below;
	\item \underline{splitting of the energy}: using a technique introduced in \cite{LM}, which is now rather standard, we decouple $ \psi^{\mathrm{GP}} = \frac{1}{\sqrt{2\pi} \varepsilon}g_{\star}\left(\frac{x-1}{\varepsilon^2}\right) u(\mathbf{x}) $ and, exploiting the variational equation satisfied by $ g_{\star} $, we obtain
		\begin{equation}
			E^{\mathrm{GP}} = \frac{E^{\mathrm{gv}}_{\star}}{\varepsilon^{4}} + \frac{\mathcal{E}[u]}{2\pi \varepsilon^2},
		\end{equation}
		with $ u $ essentially minimizing the reduced energy functional
		\begin{equation}
			\label{eq: reduced energy}
			\mathcal{E}[u] = \int \mathrm{d}\mathbf{x} \: g_{\star}^2 \left\{
			\textstyle\frac{1}{2} {|\nabla u|} ^2 + {\bf a}(\mathbf{x}) \cdot {\bf j}_u(\mathbf{x})
			+\frac{1}{2\pi\varepsilon ^4} g_{\star} ^2{(1-|u| ^2)} ^2
		\right\},
		\end{equation}
		where the ``magnetic potential'' $ \mathbf{a} $ depends on $ \Omega $ and $ {\beta_{\star}} $ and $ {\bf j}_u $ is the {\it superconducting current}
		\begin{equation}
			\label{eq: current}
			{\bf j}_u(\mathbf{x}) = \textstyle\frac{i}{2} \left( u \nabla u^* - u^* \nabla u \right).
		\end{equation}
		Completing the lower bound means to show that $ \mathcal{E}[u] $ is positive;

	\item \underline{hydrodynamic estimate}: we note that the ``magnetic potential'' is divergence free and therefore it exists a {\it potential function} $ F(x) $ such that $ 2g_{\star}^2(x) \mathbf{a}(\mathbf{x}) = - \nabla^{\perp} F(x) $. This trick was first used in \cite{CRY} in the context of the GP theory for rotating condensates. For later applications to the GL function see also \cite{CR2,CR3}. We can thus integrate by parts the second term in \eqref{eq: reduced energy} obtaining
	\begin{equation}
		\label{eq: angular momentum term}
		\int \mathrm{d} \mathbf{x} \: F(x) \: \mbox{curl} \left( {\bf j}_u \right).
	\end{equation}
	At this stage we observe that since $ {\beta_{\star}} = \mathcal{O}(1) $ and it appears in \eqref{eq: gvbf} always multiplied by $ \varepsilon^2 $, a good approximation of the functional $ \mathcal{E}^{\mathrm{gv}}_{{\beta_{\star}}} $ can be obtained by taking the limit $ \varepsilon \to 0 $, which yields the functional \eqref{eq: gvf}, with ground state energy $ E^{\mathrm{gv}} $ and minimizer $ g_{\mathrm{gv}} $. We can also replace $ F(x) $ with its limiting counterpart $ F^{\mathrm{gv}}(x) $, which is in fact a negative function. The last step to estimate \eqref{eq: angular momentum term} is to use the trivial inequality $ |\mbox{curl} \left( {\bf j}_u \right)| \leq \left| \nabla u \right|^2 $ and the negativity of $ F^{\mathrm{gv}} $ to get the lower bound
	\begin{equation}
		\mathcal{E}[u] \geq \int \mathrm{d}\mathbf{x} \:
			\left(\textstyle\frac{1}{2} g_{\mathrm{gv}}^2 + F^{\mathrm{gv}} \right) {|\nabla u|} ^2,
	\end{equation}
	where we have also dropped the last positive term in \eqref{eq: reduced energy};

	\item \underline{positivity of the cost function}: the above lower bound suggests that any topological defect of $ u $ should carry an energy cost given by the {\it cost function}
		\begin{equation}
			\label{eq: kgv}
			K^{\mathrm{gv}} = \tx\frac{1}{2} g_{\mathrm{gv}}^2 + F^{\mathrm{gv}}.
		\end{equation}
		Positivity of such a function in the bulk would then imply that vortices are not energetically favorable anywhere in the condensate. This is turn can be proven by direct inspection of the function itself. First we observe that both $ g_{\mathrm{gv}} $ and $ F^{\mathrm{gv}} $ are radial functions and we therefore change coordinates $ x = 1 + \varepsilon^2 y $, so that in the new variable $ y $ the bulk of the condensate is basically the whole real line. In the new variable the explicit expression of $ F^{\mathrm{gv}} $ (that we still denote by $ F^{\mathrm{gv}} $) is
		\begin{equation}
			\label{eq: fgv}
			F^{\mathrm{gv}}(y) = -2\Omega_0 \int_y^{\infty} \mathrm{d} t \: t \: g_{\mathrm{gv}}^2(t).
		\end{equation}
		Notice that by symmetry\footnote{Unlike $ g_{\mathrm{gv}} $, the profile $ g_{\star} $ is not exactly symmetric, but $ F $ satisfies analogous properties thanks to the optimality condition of $ {\beta_{\star}} $, i.e., the additional equation involving $ g_{\star} $ and $ {\beta_{\star}} $ which was mentioned at point 1.}  of $ g_{\mathrm{gv}} $, $ F^{\mathrm{gv}}(-\infty) = 0 $ and $ F^{\mathrm{gv}}(y) \leq 0 $ for any $ y \in \mathbb{R} $. The cost function $ K^{\mathrm{gv}} $ is therefore smooth and $ K^{\mathrm{gv}}(\pm \infty) = 0 $, so that, if it becomes negative, it must have a minimum. The derivative of $ K^{\mathrm{gv}} $ can be easily computed
		\begin{equation}
			{K^{\mathrm{gv}}}^{\prime}(y) = g_{\mathrm{gv}}(y) g_{\mathrm{gv}}^{\prime}(y) + 2\Omega_0 y g_{\mathrm{gv}}^2(y),
		\end{equation}
		so that, by strict positivity of $ g_{\mathrm{gv}} $, at any critical point $ y_0 $ for $ K^{\mathrm{gv}} $, one has
		\begin{equation}
			\label{eq: critical point}
			g_{\mathrm{gv}}^{\prime}(y_0) = - 2\Omega_0 y_0 g_{\mathrm{gv}}(y_0).
		\end{equation}
		Now using the variational equation for $ g_{\mathrm{gv}} $ and manipulating the expression \eqref{eq: fgv} of the potential function, it is possible to show that the cost function can be equivalently rewritten as
		\begin{equation}
			K^{\mathrm{gv}} = \left[\frac{1}{2}  + \Omega _0 y ^2 + \frac{\Omega _0}{\pi \alpha ^2} g_{\mathrm{gv}}^2 - \frac{2 \Omega _0 \mu^{\mathrm{gv}}}{\alpha ^2} \right] g_{\mathrm{gv}} ^2 - \frac{\Omega _0}{\alpha ^2}{g_{\mathrm{gv}}'}^2
		\end{equation}
		and, inserting the condition \eqref{eq: critical point} satisfied at any minimum point $y_0 $ of $ K^{\mathrm{gv}} $, we get
		\begin{equation}
			K^{\mathrm{gv}}(y_0) = \left[\frac{1}{2}  + \frac{\Omega _0(s+1)}{s+2} y_0 ^2 + \frac{\Omega _0}{\pi \alpha ^2} g_{\mathrm{gv}}^2(y_0) - \frac{2 \Omega _0 \mu^{\mathrm{gv}}}{\alpha ^2} \right] g_{\mathrm{gv}} ^2(y_0).
		\end{equation}
		Using the parity of $ g_{\mathrm{gv}} $ as well as the variational equation, one can prove that the quantity between brackets on the r.h.s. of the expression above is positive if and only if it is positive at the origin (see Proposition \ref{pro: kgv positive})
		\begin{equation}
			\frac{1}{2}  +  \frac{\Omega _0}{\pi \alpha ^2} g_{\mathrm{gv}}^2(0) - \frac{2 \Omega _0 \mu^{\mathrm{gv}}}{\alpha ^2} = \frac{1}{2}  +  \frac{2}{\Omega_0 (s+2)} \left[ \frac{1}{2\pi} \left\| g_{\mathrm{gv}} \right\|_{\infty}^2 -  \mu^{\mathrm{gv}} \right] \geq 0	\quad	\Longleftarrow \quad \Omega_0 \geq \Omega_{\mathrm{c}}.
		\end{equation}
\end{enumerate}

Once the energy asymptotics is proven, the pointwise estimate of $ |\psi^{\mathrm{GP}}|$, which allows to exclude the presence of vortices in the bulk for $ \Omega_0 > \Omega_{\mathrm{c}} $, is a simple consequence: putting back the positive term we have dropped in the lower bound, one first obtains an estimate of the region where $ |u| $ can differ from $ 1 $. Then combining this with an $ L^{\infty} $ estimate of the gradient of $ u $, one gets the result.

	It is worth mentioning at this stage a technical difference with previous approaches. Indeed in \cite{CPRY3} two potential functions were actually used instead of one, in order to get rid of boundary terms coming from the integration by parts described at  step 3 (see the discussion in \cite[Sect. C]{CPRY3}). Here on the opposite we are able to use only one potential function by estimating in a more refined way the boundary terms (compare, e.g., with next \eqref{eq: int by parts}). As in \cite{CPRY3} we also exploit the symmetry properties of the profile $ g_{\star} $, which is to a very good approximation invariant under reflections w.r.t. the origin.

\section{Preliminary Estimates}
\label{sec: preliminary}

Here we collect some useful technical results as well as the main properties of the effective functionals involved in the analysis. An important piece of information is contained in Section \ref{sec: positivity} where we prove the positivity of the cost function.

\subsection{Giant Vortex Functionals}
\label{sec: gv profiles}

We start by describing the derivation of the functional \eqref{eq: gvbf} from the GP energy. As anticipated in Section \ref{sec: main results} the idea is to evaluate the energy of a trial state of the form $ f(x) e^{i n \vartheta} $ in polar coordinates $ \mathbf{x} = (x, \vartheta) $ and with $ n = \Omega + \beta $. In addition we assume that $ f $ is real as it will be for any giant vortex profile. The result of a rather simple computation is
\begin{equation}
	\label{eq: proto gvf}
	\mathcal{E}^{\mathrm{GP}}[f e^{i n \vartheta}]
	= 2\pi \int_0^{\infty} \mathrm{d} x \: x  \left\{
		\tx\frac{1}{2} \left|
			\nabla f
		\right| ^2
		+ \Omega^2 \left[
			\textstyle\frac{1}{2}\left(
				x ^2
				-\frac{\Omega + \beta}{\Omega}
			\right) ^2
			\frac{1}{x ^2}
			+W(x)
		\right] f ^2
		+ \frac{1}{\varepsilon^2} f ^4
	\right\}.
\end{equation}
Exploiting the exponential smallness of $ \psi^{\mathrm{GP}} $ outside of the bulk of the condensate proven in \cite{CPRY3} and recalled in next Proposition \ref{pro: gpm exp small}, we can restrict the integration domain to the annulus
\begin{equation}
	\label{eq: ring}
	{\mathcal{A}_{\eta}}:
	=\left\{
		\mathbf{x} \in {\mathbb{R}^2} :\:
		\left|
			1
			-\mathbf{x}
		\right|
		\le \varepsilon ^2 \eta
	\right\},\qquad
	\eta:
	= \textstyle\frac{\eta _0}{2\sqrt{\Omega _0}}|\log\varepsilon|,
\end{equation}
where $ \eta_0 >  0 $ is an arbitrary finite constant and the prefactor in the definition of $ \eta $ has been chosen of that form for further convenience. Thanks to \eqref{eq: gpm exp small}
\begin{equation}
	\label{eq: change of coordinates}
	\psi^{\mathrm{GP}}(\mathbf{x}) = \mathcal{O}(\varepsilon^{\infty}),	\qquad \mbox{for any } \mathbf{x} \notin {\mathcal{A}_{\eta}},
\end{equation}
and the restriction is thus well motivated. In addition we will also see that a similar estimate holds true for any giant vortex profile. In terms of the one-dimensional functional \eqref{eq: proto gvf} we are then integrating in the interval $ [1-\varepsilon^2 \eta,1+ \varepsilon^2 \eta] $ and a change of variable is called for: setting
\begin{equation}
	x = 1 + \varepsilon ^2 y,	\qquad		g(y) = \sqrt{2\pi}\varepsilon \: f(1 + \varepsilon^2y)
\end{equation}
so that $ g $ is normalized in\footnote{We set in fact $ L^p_{\eta} : = L^p([-\eta,\eta],(1 + \varepsilon^2 y) \mathrm{d} y) $ for any $ 1 \leq p \leq \infty $.} $ L^2_{\eta} : = L^2([-\eta,\eta], (1 + \varepsilon^2y) \mathrm{d} y) $, we obtain the energy
\bml{
	\label{eq: quasi gvf}
	\tilde{\mathcal{E}}^{\mathrm{GP}}[g]
	= \frac{1}{\varepsilon ^4}
	\int _{-\eta}^{\eta}
	\mathrm{d} y \: (1 + \varepsilon^2 y)
	\left\{
		\textstyle\frac{1}{2}
		\left(
			g^{\prime}
		\right)^2 +
	\right.
	\\
	\left.
		+\varepsilon ^4 \Omega ^2
		\left[
			\textstyle\frac{1}{2}\left(
				1 + \varepsilon^2 y
				-\frac{\Omega + \beta}{\Omega\lf(1+\eps^2y\ri)}
			\right) ^2
			+W(1+\eps^2y)
		\right]g ^2
		+ \textstyle\frac{1}{2\pi}
		g ^4
	\right\}.
}
Now we expand $ W(1 + \varepsilon^2 y) $ in Taylor series around $ y = 0 $ to get
\begin{equation}
	\label{eq: taylor}
	W(1 + \varepsilon^2 y)
	= \textstyle\frac{s-2}{2} \varepsilon^4 y^2	+ \textstyle\frac{(s-2)(s-1)}{6} \varepsilon^6 y^3 + \varepsilon^8 \varphi(y)
\end{equation}
where $ \varphi (y) = \mathcal{O}(y^4) $.
Using this fact we can rewrite the potential in \eqref{eq: quasi gvf} as (recall that $ \alpha^2 = \Omega_0^2 (s + 2) $)
\begin{equation}
	\varepsilon ^4 \Omega ^2
		\left[
			\textstyle\frac{\left( 2 \varepsilon^2 y - \varepsilon^4 \beta/\Omega_0 + \varepsilon^4 y^2\right)^2}{2(1 + \varepsilon^2 y)^{2}}
			+W(1+\eps^2y)
		\right] = U_{\beta}(y) + \varepsilon^2 y^3 v(y),
\end{equation}
with $ v $ independent of $ \beta $ and of lower order w.r.t. to $ U_{\beta} $. Explicitly
\begin{equation}
	\label{eq: pot U}
	U _\beta(y):
	= \frac{1}{\left(1+\eps^2y\right) ^2}
	\left(
		\frac{\alpha ^2}{2} y ^2
		-2 \Omega _0 \varepsilon ^2 \beta y
		-\Omega _0 \varepsilon ^4 \beta y ^2
		+\frac{1}{2}\varepsilon ^4 \beta ^2
	\right),
\end{equation}
\begin{equation}
	\label{eq: pot v}
	v(y):
	= \frac{\Omega _0 ^2(s + (s-1)\varepsilon ^2 y)}{\left(1+\eps^2y\right) ^2}
	+\frac{(s-1)(s-2) \Omega_0^2}{6}
	+\frac{\varepsilon^2 \Omega _0 ^2}{y ^3} \varphi\left(y\right).
\end{equation}
Some trivial estimate using the Taylor expansion \eqref{eq: taylor} implies that for $ y \in [-\eta,\eta] $
\begin{equation}
	\label{eq: taylor U}
	U _\beta(y) = \textstyle\frac{1}{2} \alpha^2 y^2 + \mathcal{O}\left(\varepsilon^2 (1 + |\beta|) \eta + \varepsilon^4 \beta^2\right),
\end{equation}
which shows that, if, e.g., $ \beta $ is uniformly bounded in $ \varepsilon $, the potential $ U_{\beta}(y) $ is harmonic up to corrections of higher order in $ \varepsilon $. Alternatively one can think of $ U_{\beta} $ as a shifted harmonic oscillator by writing
\begin{equation}
	U_{\beta}(y) = \textstyle\frac{1}{2} \alpha^2 \left( y - \textstyle\frac{2 \Omega_0 \varepsilon^2 \beta}{\alpha^2} \right)^2 + \mathcal{O}(\varepsilon^4 |\beta| \eta^2 + \varepsilon^4 \beta^2).
\end{equation}
In fact, since the optimal value of $ \beta $ we are going to choose is $ \mathcal{O}(1) $, both representations are equivalent since the shift will be $ \mathcal{O}(\varepsilon^2) $.
Concerning the rest $ v(y) $ one trivially has the upper bound
\begin{equation}
	\label{eq: taylor v}
	\left| v(y) \right| \leq C_{\Omega_0} + \mathcal{O}(\varepsilon^2 \eta),
\end{equation}
for $ y \in [-\eta,\eta] $ and with a finite constant $ C_{\Omega_0} $. The rest in the above expression is a consequence of the bound $ |\varphi(y)| \leq C |y|^4 $, which follows from the Taylor expansion \eqref{eq: taylor}.

In conclusion we have recovered the expression \eqref{eq: gvbf}, i.e.,
\begin{displaymath}
	\mathcal{E}_\beta^\mathrm{gv}[g] = \int _{-\eta}^{\eta} \mathrm{d} y \: (1 + \varepsilon^2 y) \left\{ \tx\frac{1}{2}
		\left(
			g^{\prime}
		\right)^2
		+U _\beta(y) g ^2
		+\varepsilon ^2 y ^3 v(y) g ^2
		+\frac{1}{2\pi}
		g ^4
	\right\}.
\end{displaymath}
We now discuss the minimization of such a function w.r.t. $ g $ and for that purpose we have to identify the proper minimization domain, i.e.,
\begin{equation}
	\label{eq: gvbeta domain}
	\mathscr{D}_\beta^\mathrm{gv}:
	=\left\{
		g \in H ^1 (-\eta,\eta)\: \big|\:
		g = g ^*,\:
		\left\|
			g
		\right\| _{L ^2_{\eta}}
		= 1
	\right\}.
\end{equation}
The ground state energy of $ \mathcal{E}_\beta^\mathrm{gv} $ is defined as
\begin{equation}
	\label{eq: gvbeta energy}
	E_\beta^\mathrm{gv} =\inf _{g \in \mathscr{D}_\beta^\mathrm{gv}} \mathcal{E}_\beta^\mathrm{gv}[g].
\end{equation}
Notice that the assumption $ g = g^* $, i.e., reality of the argument, does not imply any loss of generality because the ground state can always be chosen real (see next Proposition).

\begin{pro}[Minimization of $ \mathcal{E}_\beta^\mathrm{gv} $]
	\label{pro: min gvbf}
	\mbox{}	\\
	There exists a minimizer $ g_\beta \in \mathscr{D}_\beta^\mathrm{gv} $ of \eqref{eq: gvbf} that is unique up to a sign, radial
	and can be chosen strictly positive. 	In addition $ g_\beta \in C^{\infty}(-\eta,\eta) $ and it solves the variational equation
	\begin{equation}
		\label{eq: gvbm variational}
		-\textstyle\frac{1}{2}g_\beta''
		-\frac{\varepsilon ^2}{2\lf(1+\eps^2y\ri)}g_\beta'
		+U _\beta (y) g_\beta
		+\varepsilon ^2 y ^3 v(y)g_\beta
		+\frac{1}{\pi} g_\beta ^3
		=\mu_\beta g_\beta
	\end{equation}
	with Neumann boundary conditions $ g_\beta^{\prime}(\pm \eta) = 0 $ and $ \mu_\beta = E_\beta^\mathrm{gv} + \frac{1}{2\pi} \|g_\beta\| _{L ^4_{\eta}} ^4 $. 	\\
	Finally $ g_\beta $ has a unique maximum point at $ {y_\beta}  $ and it decreases monotonically anywhere else.

\end{pro}

\begin{proof}
	Existence and uniqueness of the minimizer follow from strict convexity of the functional $ \mathcal{E}_\beta^\mathrm{gv} [\sqrt{\rho}] $ with respect to the density $ \rho = g ^2 $. The variational equation \eqref{eq: gvbm variational} is satisfied at least in weak sense. Then one deduces the strict positivity of $ g_\beta $ noticing that it is actually a ground state of a suitable one-dimensional Schr\"odinger operator. The equality for $ \mu_\beta $ follows integrating the \eqref{eq: gvbm variational} and recalling the fact that $ g_\beta $ has $ L ^2 $-norm equal to one. Finally a trivial bootstrap argument allows to deduce smoothness of $ g_\beta $ and therefore that \eqref{eq: gvbm variational} is solved in a classical sense.

	The only non-trivial result is the one about the existence of the a single maximum point for $ g_\beta $. However it follows from the property of the potential $ U_{\beta}(y) + \varepsilon^2 y^3 v(y) $: going back to the expression of the potential in \eqref{eq: proto gvf}, one can easily compute, with $ x = 1  + \varepsilon^2 y $,
	\begin{displaymath}
		\frac{\partial \left[ U_{\beta}(y) + \varepsilon^2 y^3 v(y) \right]}{\partial x} = \frac{1}{x^3} \left[ x^{s+2} - \left(1 + \textstyle\frac{\varepsilon^4 \beta}{\Omega_0} \right)^2 \right],
	\end{displaymath}
	which vanishes at a single point $ y_{\mathrm{pot}} $, i.e., where
	\begin{displaymath}
		1 + \varepsilon^2 y_{\mathrm{pot}} = \left(1 + \textstyle\frac{\varepsilon^4 \beta}{\Omega_0} \right)^{\frac{2}{s+2}} = 1 + \textstyle\frac{2 \beta \varepsilon^4}{(s+2) \Omega_0} + \mathcal{O}(\varepsilon^8 \beta^2).
	\end{displaymath}
	The Taylor expansion also shows that
	\begin{displaymath}
		y_{\mathrm{pot}} = \textstyle\frac{2 \beta \varepsilon^2}{(s+2) \Omega_0} (1 + \mathcal{O}(\varepsilon^4 \beta)).
	\end{displaymath}
	The monotonicity property of $ g_\beta $ can then be obtained by a simple rearrangement argument (see, e.g., \cite[Proposition 2.2]{CPRY1}): since the potential has a single maximum point, if $ g_\beta $ had more than one maximum besides $ {y_\beta} $, one could move mass from the further maximum to the minimum in between and lower the energy. Since $ g_\beta $ is a minimizer one gets a contradiction.
\end{proof}

\noindent
Another effective one-dimensional functional which is going to play an important role in the analysis is \eqref{eq: gvf}, i.e., the formal limit $ \varepsilon \to 0 $ of $ \mathcal{E}_\beta^\mathrm{gv} $, assuming that $ \beta = o(\varepsilon^{-2}) $:
\begin{displaymath}
	\mathcal{E}^{\mathrm{gv}}[g]
	=\int _{\mathbb{R}}
	\mathrm{d} y
	\left\{
		\textstyle\frac{1}{2}\left(
			g'
		\right)^2
		+\frac{\alpha ^2}{2} y ^2 g ^2
		+\frac{1}{2\pi} g ^4
	\right\}.
\end{displaymath}
The minimization domain is in this case given by
\begin{equation}
	\mathscr{D}^{\mathrm{gv}}:
	=\left\{
		g \in H ^1(\mathbb{R})\:\big|\:
			g
			= g ^*,
			\left\| g \right\| _{L ^2(\mathbb{R})}
			= 1
		\right\},
\end{equation}
and the ground state energy will be denoted by $ E^{\mathrm{gv}} = \inf _{g \in \mathscr{D}^{\mathrm{gv}}} \mathcal{E}^{\mathrm{gv}}[g] $.

\begin{pro}[Minimization of $ \mathcal{E}^{\mathrm{gv}} $]
	\label{pro: min gvf}
	\mbox{}	\\
	There exists a minimizer $ g_{\mathrm{gv}} \in \mathscr{D}^{\mathrm{gv}} $ of \eqref{eq: gvf} that is unique up to a sign, radial
	and can be choose strictly positive. In addition $ g_{\mathrm{gv}} \in C^{\infty}(\mathbb{R}) $ and it solves the variational equation
	\begin{equation}
		\label{eq: gvm variational}
		-\textstyle\frac{1}{2}g_{\mathrm{gv}}''
		+\frac{\alpha ^2}{2} y^2 g_{\mathrm{gv}}
		+\frac{1}{\pi} g_{\mathrm{gv}} ^3
		=\mu^{\mathrm{gv}} g_{\mathrm{gv}}
	\end{equation}
	with $ \mu^{\mathrm{gv}} = E^{\mathrm{gv}} + \frac{1}{2\pi} \|g_{\mathrm{gv}}\| _4 ^4 $.	\\
	Finally $ g_{\mathrm{gv}} $ is even w.r.t. the origin and has only one maximum at $ y = 0 $, which fulfills the inequality
	\begin{equation}
		\label{eq: gvm infty bound}
		g_{\mathrm{gv}}^2(0) = \left\|g_{\mathrm{gv}}\right\| _\infty ^2 \leq \pi \mu^{\mathrm{gv}}.
	\end{equation}
\end{pro}

\begin{proof}
	See the proof of Proposition \ref{pro: min gvbf}. Parity of $ g_{\mathrm{gv}}$  is a trivial consequence of the parity of the potential. The inequality \eqref{eq: gvm infty bound} follows from direct inspection of the variational equation \eqref{eq: gvm variational}: at any maximum point $ g_{\mathrm{gv}}'' \leq 0 $, which immediately implies the result.
\end{proof}

\subsection{Estimates of the Gross-Pitaevskii and Giant Vortex Profiles}

In this Section we collect several technical estimates of the profiles involved in the discussion. Such estimates will play a key role in the proofs but can be typically obtained by standard techniques in functional analysis.

We start by recalling a result which was in fact proven in \cite[Propositions 3.1 and 3.2]{CPRY3}: {let $ \eta_0 $ be the parameter appearing in the definition \eqref{eq: ring} of $ {\mathcal{A}_{\eta}} $, then}

	\begin{pro}[Exponential decay of $ \psi^{\mathrm{GP}} $]
		\label{pro: gpm exp small}
		\mbox{}	\\
		If $ \Omega = \Omega_0 / \varepsilon^4 $, there exists two finite constants $ c, C > 0 $ (independent of $ \eta_0 $) such that, for any $ \mathbf{x} \notin {\mathcal{A}_{\eta}} $,
		\begin{equation}
			\label{eq: gpm exp small}
			\left| \psi^{\mathrm{GP}}(\mathbf{x}) \right|^2 \leq \frac{C}{\varepsilon^2} \max\left[ \varepsilon^{\frac{c \eta_0^2}{4}}, \exp\left\{ - \textstyle\frac{\sqrt{\Omega_0}}{\varepsilon^2} |1 - x|\right\} \right].
		\end{equation}
	\end{pro}

	\noindent
	In particular the above result implies that by taking $ \eta_0 $ large enough we can make $ \psi^{\mathrm{GP}} $ arbitrarily small outside $ {\mathcal{A}_{\eta}} $. This fact will be crucial in restricting the computation of the GP energy within $ {\mathcal{A}_{\eta}} $. Notice also that as soon as $ |1 - x| \gg \varepsilon^2|\log\varepsilon| $, $ \psi^{\mathrm{GP}} = \mathcal{O}(\varepsilon^{\infty}) $.

	Let us now focus on the giant vortex profiles. Before stating the main technical estimates we first formulate a simple preliminary bound on the giant vortex energy $ E_\beta^\mathrm{gv} $:

	\begin{pro}[Preliminary bound on $ E_\beta^\mathrm{gv} $]
		\label{pro: estimates on energy}
		\mbox{}	\\
		If $ \beta = \bigO{\varepsilon ^{-2}} $ as $ \varepsilon \to 0$, then
		\begin{equation}
			\label{eq: en preliminary est}
			E_\beta^\mathrm{gv} = \mathcal{O}(1),	\qquad
			\mu_\beta
			= \mathcal{O}(1).
		\end{equation}
	\end{pro}

\begin{proof}
	Since $ E_\beta^\mathrm{gv} $ is positive (compare with \eqref{eq: quasi gvf}), it suffices to prove a suitable upper bound: to that purpose one can simply evaluate the functional $ \mathcal{E}_\beta^\mathrm{gv} $ on the ground state of the one-dimensional harmonic oscillator with frequency $ \alpha $. The result on $ \mu_\beta $ follows from the trivial estimates $ E_\beta^\mathrm{gv} \leq \mu_\beta \leq 2 E_\beta^\mathrm{gv} $.
\end{proof}

\noindent
The giant vortex profile $ g_\beta $ decays exponentially for large $ |y|$ and one can actually show that this decay captures the correct asymptotics of $ g_\beta $:

	\begin{pro}[Pointwise estimates of $ g_\beta $]
		\label{pro: gvbm point est}
		\mbox{}	\\
		If $ \beta = \bigO{\varepsilon ^{-2}} $ as $ \varepsilon \to 0 $, then there exists a finite constant $ C $ such that
		\begin{equation}
			\label{eq: estimate on gvbm}
			g_\beta(y)
			\le C e ^{-2\sqrt{\Omega _0}|y|},	\qquad		\mbox{so that } g_\beta(\pm \eta) =\bigO{\varepsilon ^{\eta _0}}.
		\end{equation}
		If  $ \beta = \bigO{1} $ then there exist two finite constants $ C _1, C _2 > 0 $ such that the following inequalities hold true
		\begin{equation}
			\label{eq: better estimate on gvbm}
			C_1 \| g_\beta \| _{L _\eta ^4} ^2
			\exp\left\{
				- \textstyle\frac{\alpha}{2} y ^2
			\right\}
			\le g_\beta(y)
			\le C_2
			\exp\left\{
				- \textstyle\frac{\alpha}{4} y ^2
			\right\}.
		\end{equation}
	\end{pro}

	\begin{proof}
		The results are proven by means of standard super- and sub-solution techniques. We spell however the proofs in full details for the sake of clarity.

		To prove \eqref{eq: estimate on gvbm} it is somehow more convenient to go back to the variational equation satisfied by $ f_{\beta}(x) = (\sqrt{2\pi}\varepsilon)^{-1} \: g_\beta((x-1)/\varepsilon^2) $, i.e.,
		\begin{equation}
			\label{eq: vareq for f}
			-\textstyle\frac{1}{2} f^{\prime\prime}_\beta - \frac{1}{2r} f^{\prime}_\beta + \textstyle\frac{\Omega ^2}{2x ^2} \left (
			x ^2
			- \frac{\Omega + \beta}{\Omega}
			\right) ^2
			f _\beta
			+\Omega ^2 W (x) f _\beta
			+\frac{1}{\varepsilon ^2} f _\beta ^3
			=\frac{1}{\varepsilon ^4} \mu_\beta f _\beta.
		\end{equation}
		The first simple observation is that by positivity of $ f_{\beta} $ and $ W(x) $, we get
		\[
			-\textstyle\frac{1}{2} f^{\prime\prime}_\beta - \frac{1}{2r} f^{\prime}_\beta
			\le \frac{1}{\varepsilon ^4}
			\left(
			\mu _\beta
			-\varepsilon ^2 f _\beta ^2
			\right) f _\beta,
		\]
		which, by negativity of the second derivative of $ f_\beta $ at any maximum point of $ f_{\beta} $, immediately implies the upper bound
		\begin{displaymath}
			\left\| f_\beta \right\|_{L^\infty({\mathcal{A}_{\eta}})}^2 \leq \textstyle\frac{1}{\varepsilon^2} \mu_\beta,
		\end{displaymath}
		which in terms of $ g_\beta $ becomes, via \eqref{eq: en preliminary est} (here we are assuming that $ \beta = \mathcal{O}(\varepsilon^{-2}) $),
		\begin{equation}
			\label{eq: gvbm infty est}
			\left\|g_\beta \right\|_{L^\infty_{\eta}} = \mathcal{O}(1).
		\end{equation}

		In order to prove \eqref{eq: estimate on gvbm} we will provide an explicit supersolution to the equation \eqref{eq: vareq for f}. Notice that the first two terms of the equation form the two-dimensional Laplacian, i.e., for any radial function $ f $, $ - \Delta f = - \frac{1}{r} \partial_r ( r f^{\prime}) $. We will use this fact to construct a supersolution in dimension two. Let then $ a > 0 $ be a parameter independent of $ \varepsilon $ that is going to be chosen later and consider the two-dimensional region
		\begin{equation}
			\mathcal{A}:
			= \mathcal{B} _{1 - a \varepsilon ^2}(0)
			\cap {\mathcal{A}_{\eta}}
			=\left\{
			\mathbf{x} \in {\mathbb{R}^2}
			\;\big|\;
			1 - \eta \varepsilon ^2
			\leq x
			\leq 1 - a \varepsilon ^2
			\right\}.
		\end{equation}
		Inside $ \mathcal{A} $ one has the lower bound
		\[
			W (x)
			\ge \textstyle\frac{s-2}{2}\left(
			x-1
			\right) ^2
			+\bigO{\left|x-1\right| ^3}
			\ge \frac{\left(s-2\right)}{2} a ^2 \varepsilon ^4
			+\bigO{\eta ^3 \varepsilon ^6}
			\ge C _0 a ^2 \varepsilon ^4
		\]
		with $ C _0 > 0 $, so that \eqref{eq: vareq for f} and \eqref{eq: en preliminary est} yield
		\[
			- \textstyle\frac{1}{2} \Delta f _\beta
			\le \frac{1}{\varepsilon ^4} \mu_\beta f _\beta
			-\Omega ^2 W(x) f _\beta
			\le \frac{1}{\varepsilon ^4}
			\left(
			C _1
			-\Omega _0 ^2 C _0 a ^2
			\right)f _\beta
		\]
		where also $ C_1 > 0 $. If now we pick $ a ^2 \ge \frac{2 C _1}{C _0 \Omega _0 ^2} $, we get that $ f _\beta $ is a subsolution of the following differential problem
		\begin{equation}
			\label{eq: f subsol}
			-\textstyle\frac{1}{2} \Delta f
			+\frac{1}{2} C _0 a ^2 \Omega ^2 \varepsilon ^4 f
			=0.
		\end{equation}
		To get rid of the inner boundary we now extend $ f_\beta $ to the whole ball $ \mathcal{B}_{1 - a\varepsilon^2}(0) $ in a smooth (in fact at least $ C^2 $) way. We denote by $ \tilde{f} $ such a new function and we require that $ \tilde{f}(x) = 0 $ for $ x \leq 1 - 2 \varepsilon^2 \eta $ and
		\begin{equation}
			 -\textstyle\frac{1}{2} \Delta f
			+\frac{1}{2} C _0 a ^2 \Omega ^2 \varepsilon ^4 f
			\leq 0,
		\end{equation}
		for any $ \mathbf{x} \in \mathcal{B}_{1 - a\varepsilon^2}(0) $. We omit the explicit details of such a construction for the sake of brevity.

		A supersolution to the same problem can be constructed by taking $$ f _\mathrm{sup}(x):=C _a \| f _\beta \| _\infty e ^{-\sqrt{\Omega}\left(1-x ^2\right)}$$ with $ C _a $ a constant to be suitably chosen:
		\[
			-\textstyle\frac{1}{2} \Delta f _\mathrm{sup}
			+\frac{1}{2} C _0 a ^2 \Omega ^2 \varepsilon ^4 f _\mathrm{sup}
			=\frac{1}{2}\left(
			-2\sqrt{\Omega}
			-2 \Omega x ^2
			+ C _0 \Omega _0 a ^2 \Omega
			\right)  f _\mathrm{sup} >0,
		\]
		if we choose $ a ^2 > \frac{4}{C _0 \Omega _0} $. The constant $ C _a $ is then used to guarantees that $ f_{\mathrm{sup}} $ satisfies the proper boundary conditions. In order to apply the maximum principle (see, e.g., \cite[\S~6.4.1, Theorem 2]{E}), we need that $ \tilde{f}(x) \le f _\mathrm{sup} (x) $ on $ \partial \mathcal{B}_{1 - a\varepsilon^2}(0) $, which holds true if  $ C _a \geq e ^{2\sqrt{\Omega _0}a} $:
		\[
			\left.
				f _\mathrm{sup}
			\right| _{\partial \mathcal{B}_{1 - a\varepsilon^2}(0)}
			= C _a \| f _\beta \| _\infty e^{-\sqrt{\Omega _0}a \left(2-a\varepsilon ^2\right)}
			\ge \left.
				f_{\beta}
			\right| _{\partial \mathcal{B}_{1 - a\varepsilon^2}(0)}.
		\]
		Hence we conclude that $ \tilde{f} \le f _\mathrm{sup} $ in the whole $ \mathcal{B}_{1 - a\varepsilon^2}(0) $, and therefore, using the monotonicity of $ f _\mathrm{sup} $, $ f_\beta \leq  f _\mathrm{sup} $ in the whole region $ \mathcal{B} _1 (0) \cap {\mathcal{A}_{\eta}} $. Going back to $ g_\beta $ and using \eqref{eq: gvbm infty est}, we obtain \eqref{eq: estimate on gvbm} in $ \mathcal{B} _1 (0) \cap {\mathcal{A}_{\eta}} $. To extend the result to the complementary region, one can use a very similar argument with the trivial change $ x^2 - 1 \to 1 - x^2 $ in the supersolution.

		For the refined estimates \eqref{eq: better estimate on gvbm}, we consider the variational equation \eqref{eq: gvbm variational} for $ a \leq |y| \leq \eta $, with $ a > 0 $ such that $ a ^2 > \frac{8\mu_\beta}{3\alpha ^2} $ and $ \varepsilon $ small enough, which imply

	%

		\[
			U _\beta (y)
			+\varepsilon ^2 y ^3 v(y)
			-\mu_\beta
			= \textstyle\frac{\alpha ^2}{2} y ^2
			-\mu_\beta
			+\bigO{\varepsilon ^2 \eta ^3}
			\ge \frac{\alpha ^2}{8} y ^2
		\]
		and therefore in that region $ g_\beta $ is a subsolution of the equation
		\begin{equation}
			\label{eq: subequation}
			-\textstyle\frac{1}{2} g''
			-\frac{\varepsilon ^2}{2\lf(1+\eps^2y\ri)} g'
			+\frac{\alpha ^2}{8} y ^2 g = 0.
		\end{equation}
		As before we extend $ g_\beta $ to the whole region $ |y| \geq a $ in a $C^2$ way and preserving the differential inequality satisfied in $ a \leq |y| \leq \eta $, i.e.,
		 \begin{displaymath}
			-\textstyle\frac{1}{2} g''
			-\frac{\varepsilon ^2}{2\lf(1+\eps^2y\ri)} g'
			+\frac{\alpha ^2}{8} y ^2 g \leq 0.
		\end{displaymath}
		Again we skip the details for brevity.

		Now for some $ C > 0 $ to be fixed later the following function
		\[
			g _\mathrm{sup} (y):
			=C e ^{-\frac{\alpha}{4}y ^2}
		\]
		is a supersolution to \eqref{eq: subequation}: for $ \varepsilon $ small enough
		\[
			-\textstyle\frac{1}{2} g _\mathrm{sup}''
			-\frac{\varepsilon}{2\lf(1+\eps^2y\ri)} g _\mathrm{sup}'
			+\frac{\alpha ^2}{8} y ^2g _\mathrm{sup}
			=\left(
			\frac{\alpha}{4}
			+\frac{\alpha y \varepsilon ^2}{4\lf(1+\eps^2y\ri)}
			\right)g _\mathrm{sup}
			\geq 0.
		\]
		Choosing the $ C \geq \|g_\beta\| _\infty e ^{\frac{\alpha a ^2}{4}} $ to ensure that $ g_\beta(\pm a) \leq g_{\mathrm{sup}}(a) $, we get the upper estimate.
	%

	%

	%

		Analogously we can choose $ C > 0 $ in such a way that
		\[
			g_{\mathrm{sub}}:
			=C e ^{-\frac{\alpha}{2}y ^2}
		\]
		is a subsolution to \eqref{eq: gvbm variational}: first one notes that
		\[
			\textstyle\frac{\alpha}{2}
			\le \mu_\beta
			-\frac{1}{\pi} \| g_\beta \| _{L _\eta ^4} ^4
			+\bigO{\varepsilon ^{2} \eta},
		\]
		which follows from the fact that the harmonic oscillator on the real line is bounded from below by $ \alpha/2 $; then using this inequality in \eqref{eq: gvbm variational}, we obtain
		\bmln{
			 -\textstyle\frac{1}{2} g_{\mathrm{sub}}''
	 		-\frac{\varepsilon ^2}{2\lf(1+\eps^2y\ri)} g_{\mathrm{sub}}'
			 +U _\beta (y) g_{\mathrm{sub}}
			 +\varepsilon ^2 y ^3 v(y) g_{\mathrm{sub}}
			 +\frac{1}{\pi} g_{\mathrm{sub}}^3
			 - \mu_\beta g_{\mathrm{sub}}	\\
			=\left[
			\textstyle\frac{\alpha}{2}
			-\mu_\beta
			+\frac{1}{\pi} g_{\mathrm{sub}}^2
			+U _\beta (y)
			-\frac{\alpha ^2}{2}y ^2
			+\frac{\alpha y \varepsilon ^2}{2\lf(1+\eps^2y\ri)}
			+\varepsilon ^2 y ^3 v(y)
			\right] g_{\mathrm{sub}}	\\
			\le \left[
			\textstyle\frac{1}{\pi}\left(
				 g _\mathrm{sub} ^2 (y)
				 -\|g_\beta\| _{L _\eta ^4} ^4
			\right)
			+\bigO{\varepsilon ^2 \eta ^3}
			\right] g _\mathrm{sub}
			< 0,
		}
		if we pick $ C < \|g_\beta\| _{L _\eta ^4} ^2 $. To conclude we use the fact that $  g_{\mathrm{sub}} $ goes to 0 as $ |y| $ goes to infinity: indeed it is sufficient to observe that there certainly exists a point $ \bar{y} >0  $ such that $  g_{\mathrm{sub}} (\pm\bar{y}) = \min \left\{g_\beta (\eta), g_\beta (-\eta) \right\} $ and
		\[
			\widetilde{g}(y):
			=\left\{
				\begin{array}{ll}
					g_\beta(y) & |y| \le \eta, \\
					g_\beta(\eta) & \eta \le y \le \bar{y}, \\
					g_\beta(-\eta) & -\bar{y} \le y \le -\eta,
				\end{array}
			\right.
		\]
		 is a supersolution to \eqref{eq: gvbm variational}, satisfying $ \tilde{g}(\pm \bar{y}) \geq g_{\mathrm{sub}}(\pm \bar{y}) $. Hence $ g_{\mathrm{sub}} \leq \tilde{g} $ for any $ |y| \leq \bar{y} $, which implies the lower estimate \eqref{eq: better estimate on gvbm} for $ |y| \leq \eta $.
	\end{proof}

\noindent
We conclude this Section by stating analogous pointwise estimate for the limiting profile $ g_{\mathrm{gv}} $:

	\begin{pro}[Pointwise estimates of $ g_{\mathrm{gv}} $]
		\label{pro: point est gvm}
		\mbox{}	\\
		There exists a finite constant $ C >0 $ such that
		\begin{equation}
			\label{eq: estimate on gvm}
			\|
				g_{\mathrm{gv}}
			\| _4 ^2 \exp\left\{-\textstyle\frac{\alpha}{2}y ^2\right\}
			\le g_{\mathrm{gv}}(y)
			\le C \exp\left\{-\textstyle\frac{\alpha}{4}y ^2\right\}.

		\end{equation}
	\end{pro}

	\begin{proof}
		The estimate can be proven exactly as \eqref{eq: better estimate on gvbm} in Proposition \ref{pro: gvbm point est} and we skip the details.

	%

	%

	%

	%

	%
	\end{proof}

\subsection{Optimal Giant Vortex Phase and Profile}
\label{sec: betas}

In this Section we investigate the minimization of $ E_\beta^\mathrm{gv} $ w.r.t. $ \beta \in \mathbb{R} $. The main result is the following

	\begin{pro}[Optimal phase]
		\label{pro: betas}
		\mbox{}	\\
		For $ \varepsilon $ small enough there exists a unique minimizer $ {\beta_{\star}} \in \mathbb{R} $ such that
		\begin{equation}
			\label{eq: minimal energy}
			E_{\star}^\mathrm{gv}
			: =\inf _{\beta \in \mathbb{R}} E_\beta^\mathrm{gv} = E^{\mathrm{gv}}_{{\beta_{\star}}}.
		\end{equation}
	 	Such an optimal phase is explicitly given by
		\begin{equation}
			\label{eq: betas}
			{\beta_{\star}}
			= - \frac{2}{\Omega _0 \left(s-2\right)}\left[
			(s-2) V - Q
			+\bigO{\varepsilon ^2}
			\right],
		\end{equation}
		where we set $ g_{\star} : = g_{{\beta_{\star}}} $ and
		\begin{equation}
			\label{eq: defining V Q}
			V:
			=\frac{\alpha ^2}{2}
			\int _{-\eta} ^\eta
			\mathrm{d} y \:
			y ^2 g_{\star}^2,\qquad
			Q:
			=\frac{1}{2\pi}
			\int _{-\eta} ^\eta
			\mathrm{d} y \:
			g_{\star} ^4.
		\end{equation}
	\end{pro}

	\begin{proof}
		The existence of a minimizer $ {\beta_{\star}} $ is guaranteed from the fact that
		\[
			U _\beta (y)
			\ge \textstyle\frac{1}{\lf(1+\eps^2y\ri) ^2}\left[
			\frac{s-2}{2(s+2)}\varepsilon ^4 \beta ^2
			-\Omega _0 \varepsilon ^4 \beta \eta
			\right]
		\]
		which implies that $ \lim_{|\beta| \to \infty} E_\beta^\mathrm{gv} = + \infty $ (recall that $ s > 2 $).  By the same lower bound on the potential together with the trivial bound $ E_{\star}^\mathrm{gv} \le E _0 ^\mathrm{gv} = \mathcal{O}(1) $, we also deduce that $ {\beta_{\star}} = \bigO{\varepsilon ^{-2}} $.

		In order to find the explicit expression of $ {\beta_{\star}} $, we first observe that by standard arguments $ E_\beta^\mathrm{gv} $ is a smooth function of $ \beta $ and therefore by the Feynman-Hellmann principle\footnote{The notation $ \braket{\cdot}{\cdot}_{\eta} $ stands for the scalar product in $ L^2([-\eta,\eta], (1 + \varepsilon^2 y) \mathrm{d} y) $.}
		\begin{equation}
			\partial _\beta E_\beta^\mathrm{gv}
			=\left\langle
				g_\beta\left|
				\partial _\beta U _\beta
			\right|g_\beta
			\right\rangle _\eta = \mean{g_\beta}{\textstyle\frac{\varepsilon ^2}{\lf(1+\eps^2y\ri) ^2}\left(
				-2\Omega _0 y
				-\Omega _0 \varepsilon ^2 y ^2
				+\varepsilon ^2 \beta
			\right)}{g_\beta}_{\eta}.
		\end{equation}
		Since $ {\beta_{\star}} $ is a minimizer, we must have $ \left. \partial_\beta E_\beta^\mathrm{gv}\right|_{{\beta_{\star}}} = 0 $, i.e.,
		\begin{equation}
			\label{eq: solved by bcr}
			\varepsilon ^2 {\beta_{\star}}
			\inpr{g_{\star}}{\textstyle\frac{1}{\lf(1+\eps^2y\ri) ^2}}{g_{\star}}
			-2\Omega _0
			\inpr{g_{\star}}{\textstyle\frac{y}{\lf(1+\eps^2y\ri) ^2}}{g_{\star}}
			-\Omega _0 \varepsilon ^2
			\inpr{g_{\star}}{\textstyle\frac{y ^2}{\lf(1+\eps^2y\ri) ^2}}{g_{\star}}
			=0.
		\end{equation}
		We compute the first and last terms of the expression above:
		\begin{equation}
			\label{eq: exp val betas}
			\inpr{g_{\star}}{\textstyle\frac{1}{\lf(1+\eps^2y\ri) ^2}}{g_{\star}}
			=1
			+\bigO{\varepsilon ^2},\qquad
			\inpr{g_{\star}}{\textstyle\frac{1}{\lf(1+\eps^2y\ri) ^2}y ^2}{g_{\star}}
			= \frac{2V}{\alpha ^2}
			+\bigO{\varepsilon ^2}.
		\end{equation}
		Indeed thanks to the exponential decay proven in \eqref{eq: estimate on gvbm}, one can easily realize that
		\begin{equation}
			\label{eq: int powers}
			\int_{-\eta}^{\eta} \mathrm{d} y \: |y|^k \: g_{\star}^2 = \mathcal{O}(1),	\qquad		\mbox{for any } k < \infty.
		\end{equation}
		In fact an analogous estimate holds true if $ g_{\star} $ is replaced with $ g_{\star}^{\prime} $, in particular
		\begin{equation}
			\label{eq: int powers derivative}
			\int_{-\eta}^{\eta} \mathrm{d} y \: y \: \left( g_{\star}'\right)^2 = \mathcal{O}(1).
		\end{equation}
		To see this it suffices to integrate by parts and use the variational equation \eqref{eq: gvbm variational} to go back to an expression involving only $ g_{\star} $ and there one can use the above estimate. We omit the details for the sake of brevity.	Notice that at this stage we are implicitly exploiting the bound $ {\beta_{\star}} = \mathcal{O}(\varepsilon^{-2}) $, which is among the hypothesis of Proposition \ref{pro: gvbm point est}. Next we integrate by parts the second term in \eqref{eq: solved by bcr} to get
		\bmln{
			\inpr{g_{\star}}{\textstyle\frac{1}{\lf(1+\eps^2y\ri) ^2}y}{g_{\star}}
			=\int _{-\eta} ^\eta
			\mathrm{d} y\
			\textstyle\frac{1}{1+\eps^2y} \: y \: g_{\star} ^2 \\
			= \left[
			\textstyle\frac{1}{2\lf(1+\eps^2y\ri)} \:y ^2 \: g_{\star} ^2
			\right] _{-\eta} ^\eta
			- \displaystyle\int _{-\eta} ^\eta
			\mathrm{d} y\
			\textstyle\frac{1}{1+\eps^2y}\: y ^2 \: g_{\star} g_{\star}'
			+ \displaystyle\frac{\varepsilon ^2V}{\alpha ^2}
			\left(
			1
			+\bigO{\varepsilon ^2}
			\right)

		}
		where the boundary terms (first term on the r.h.s. of the expression above) can be included in the remainder $ \mathcal{O}(\varepsilon^2) $ if we choose $ \eta _0 > 2 $ (see again \eqref{eq: estimate on gvbm}). For the rest we can compute
		\bmln{
			-\int _{-\eta} ^\eta
			\mathrm{d} y \:
			\textstyle\frac{1}{1+\eps^2y} \: y ^2 \: g_{\star} g_{\star} ' = -
			\displaystyle\frac{2}{\alpha ^2}\displaystyle\int _{-\eta} ^\eta
			\mathrm{d} y \: (1 + \varepsilon^2 y) \:
			U _{\beta_{\star}}(y)
			g_{\star}
			g_{\star}'	\\
			- \frac{2}{\alpha ^2}
			\displaystyle\int _{-\eta} ^\eta
			\mathrm{d} y \:
			\frac{
				2\Omega _0
				\varepsilon ^2 {\beta_{\star}} y
				+\Omega _0
				\varepsilon ^4 {\beta_{\star}} y ^2
				-\textstyle\frac{1}{2}\varepsilon ^4 {\beta_{\star}} ^2
			}{1 + \varepsilon^2 y} \:
			g_{\star}\gcrm'	\\
			=\frac{1}{\alpha ^2}
			\int _{-\eta} ^\eta
			\diff y \: \xinyp \:
			\partial_y \left[
				-\textstyle\frac{1}{2} \left(
				g_{\star}'
				\right) ^2
				+\frac{1}{2\pi} g_{\star} ^4
				-\mu_{\star} g_{\star} ^2
			\right]
			+\frac{\varepsilon ^2}{\alpha ^2}
			\int _{-\eta} ^\eta
			\diff y \: \xinyp \:
			y ^3
			v(y)
			\partial_y g_{\star} ^2 	\\
			-\frac{\varepsilon ^2}{\alpha ^2}
			\int _{-\eta} ^\eta
			\mathrm{d} y \:
			\left(
				g_{\star}'
				\right) ^2
			-\frac{\varepsilon ^2 {\beta_{\star}}}{\alpha ^2}
			\int _{-\eta} ^\eta \mathrm{d} y \: \frac{
				2\Omega _0
				y
				+\Omega _0
				\varepsilon ^2  y ^2
				-\textstyle\frac{1}{2}\varepsilon ^2{\beta_{\star}}
			}{1 + \varepsilon^2 y} \:
			\partial_y g_{\star}^2=	\\
			=\frac{\varepsilon ^2}{\alpha ^2}
			\int _{-\eta} ^\eta
			\mathrm{d} y\
			\left\{
				-\textstyle\frac{1}{2} \left(
				g_{\star}'
				\right) ^2
				-\frac{1}{2\pi} g_{\star} ^4
				+\mu_{\star} g_{\star} ^2
			\right\}
			 - \frac{\varepsilon^2(s+1) V}{\alpha^2}+ \frac{2 \varepsilon^2 \Omega_0 {\beta_{\star}}}{\alpha^2} + \mathcal{O}(\varepsilon ^4 {\beta_{\star}})
			+\bigO{\varepsilon ^{4}}=	 	\\
			=\frac{\varepsilon ^2}{\alpha^2}
			\left[
				-K
				-Q
				+\mu_{\star}
				- (s+1) V + 2 \Omega_0 {\beta_{\star}} +\bigO{\varepsilon ^2 {\beta_{\star}}}
				+\bigO{\varepsilon ^2}
			\right],
			}
			where we have made use repeatedly of \eqref{eq: int powers} and exploited the identity
			\begin{displaymath}
				\left( y^3 v(y) \right)^{\prime} = \textstyle\frac{1}{2} \alpha^2 (s+1) y^2 + \mathcal{O}(\varepsilon^2 |y|^3).
			\end{displaymath}
			We have also set
			\begin{equation}
				T :
			= \frac{1}{2}
		\int _{-\eta} ^\eta
		\mathrm{d} y \: \left(
			g_{\star}'
		\right)^2.
			\end{equation}
			Hence
			\begin{equation}
				\label{eq: implicit per bcr}
				\inpr{g_{\star}}{\textstyle\frac{1}{\lf(1+\eps^2y\ri) ^2}y}{g_{\star}}
				=\frac{\varepsilon ^2}{\alpha^2}
				\left[
				-T - s V
				-Q
				+\mu_{\star} + 2 \Omega_0 {\beta_{\star}}
				+\bigO{\varepsilon ^2 {\beta_{\star}}}
				+\bigO{\varepsilon ^2}
				\right]
			\end{equation}
			and plugging this together with \eqref{eq: exp val betas} into \eqref{eq: solved by bcr}, we obtain
			\[
				\frac{s-2}{s+2} {\beta_{\star}} (1 + \mathcal{O}(\varepsilon^2))
				+ \frac{2\Omega _0}{\alpha ^2}
				\left[
				(s-2) V - Q
				\right]
				+ \mathcal{O}(\varepsilon ^2)
				=0,
			\]
		since $ \mu_{\star} = T + V + 2Q + \mathcal{O}(\varepsilon^2) $ (see \eqref{eq: int powers} and \eqref{eq: int powers derivative}). The expression \eqref{eq: betas} is then recovered.
	\end{proof}

\noindent
Along the proof we have also proven in \eqref{eq: solved by bcr} that
\begin{equation}
	\label{eq: optimality}
	\int_{-\eta}^{\eta} \mathrm{d} y \: \textstyle\frac{1}{1 + \varepsilon^2 y} \left(y + \frac{1}{2} \varepsilon^2 y^2 - \frac{1}{2 \Omega_0} \varepsilon^2 {\beta_{\star}} \right) g_{\star}^2 = 0,
\end{equation}
which, thanks to the result about $ {\beta_{\star}} $, also implies that
\begin{equation}
	\label{eq: gcrm almost symm}
	\inpr{g_{\star}}{y}{g_{\star}}
	= \mathcal{O}(\varepsilon^2),
\end{equation}
i.e., the profile $ g_{\star} $ is almost symmetric w.r.t. the origin.

In fact this latter information can be deduced also by looking at the relation between the functional $ \mathcal{E}^{\mathrm{gv}}_{{\beta_{\star}}} $ and its minimization and the limiting model $ \mathcal{E}^{\mathrm{gv}} $. From now on we fix $ \beta $ equal to the optimal value $ {\beta_{\star}} $.

Before discussing this question further we have however to state an useful estimate on $ g_{\star} $.

	\begin{lem}
		\mbox{}	\\
		There exists a finite constant $ C > 0 $ such that for any $ y \in [-\eta,\eta] $
		\begin{equation} \label{eq: derivative gcrm}
			\left|
				g_{\star} ' (y)
			\right|
			\le C \eta ^3 g_{\star} (y).
		\end{equation}
	\end{lem}

	\begin{proof}
		It suffices to integrate the variational equation \eqref{eq: gvbm variational} between $ y \geq y_{{\beta_{\star}}} $ and $ \eta $ or $ -\eta $ and $ y \leq y_{{\beta_{\star}}} $ (recall that $ y_{{\beta_{\star}}} $ stands for the unique maximum point of $ g_{\star} $): let us assume that $ y \geq {y_\beta} $, then by positivity of $ g_{\star} $
		\bml{
			\textstyle\frac{1}{2} \left|
			g_{\star}'(y)
			\right|
			=-\frac{1}{2} g_{\star}' (y)
			=\displaystyle\int _{y} ^{\eta}
			\mathrm{d} t\:
			\left\{ - \textstyle\frac{\varepsilon^2}{2(1 + \varepsilon^2 t)} g_{\star}' +
			\left(\textstyle\frac{1}{\pi} g_{\star} ^2
			+ U_{{\beta_{\star}}}(t) + \varepsilon^2 t^3 v(t)
			- \mu^{\mathrm{gv}}_{{\beta_{\star}}}
			\right) g_{\star} \right\}		\\
			\leq \displaystyle\int _{y} ^{\eta}
			\mathrm{d} t\:
			\left\{ - \textstyle\frac{\varepsilon^4}{2(1 + \varepsilon^2 t)^2} +
			\textstyle\frac{1}{\pi} g_{\star} ^2
			+ U_{{\beta_{\star}}}(t) + \varepsilon^2 t^3 v(t)
			 \right\}g_{\star} +  \textstyle\frac{\varepsilon^2}{2(1 + \varepsilon^2 y)} g_{\star}(y).
		}
		Now given that the quantity between brackets can be easily bounded from above by $ C \eta^2 $, it only remains to use the monotonicity of $ g_{\star} $ to conclude the proof.
	\end{proof}

	\noindent
	We are now in position to prove the first result about the energy difference $ E_{\star}^\mathrm{gv} - E^{\mathrm{gv}} $. As a matter of fact this will involve a corresponding statement about the closeness of $ g_{\star}^2 $ to $ g_{\mathrm{gv}}^2 $ in $ L^2 $. {We recall the expressions of the energy functionals
	\begin{displaymath}
		\mathcal{E}_{\star}^{\mathrm{gv}}[g] = \int _{-\eta}^{\eta} \mathrm{d} y \: (1 + \varepsilon^2 y) \left\{ \tx\frac{1}{2}
		\left(
			g^{\prime}
		\right)^2
		+U _{{\beta_{\star}}}(y) g ^2
		+\varepsilon ^2 y ^3 v(y) g ^2
		+\frac{1}{2\pi}
		g ^4
	\right\},
	\end{displaymath}
	\begin{displaymath}
		\mathcal{E}^{\mathrm{gv}}[g] = \int _{\mathbb{R}} \mathrm{d} y
		\left\{
		\textstyle\frac{1}{2}\left(
			g '
		\right) ^2
		+ \frac{\alpha ^2}{2} y^2 g ^2
		+\frac{1}{2\pi} g ^4
		\right\}.
	\end{displaymath}}

	\begin{pro}[Estimate of $ E_{\star}^\mathrm{gv} - E^{\mathrm{gv}} $]
		\label{pro: gv energy diff}
		\mbox{}	\\
		As $ \varepsilon \to 0 $
		\begin{equation}
			\label{eq: gv energy diff}
			E_{\star}^\mathrm{gv}
			=E^{\mathrm{gv}}
			+\bigO{\varepsilon ^4 \eta ^7},	\qquad			\left\|
			g_{\star} ^2
			-g_{\mathrm{gv}} ^2
			\right\| _{L ^2 _\eta}^2
			=\bigO{\varepsilon ^4 \eta ^{7}}.
		\end{equation}
	\end{pro}

	\begin{proof}
		We test the two functionals $ \mathcal{E}^{\mathrm{gv}}_{{\beta_{\star}}} $ and $ \mathcal{E}^{\mathrm{gv}} $ on suitable test functions. Let us first regularize $ g_{\star} $ outside $ [-\eta,\eta] $ to make it an admissible test function for $ \mathcal{E}^{\mathrm{gv}} $: we define
		\[
			g _\mathrm{trial}(y):
			= c_{\varepsilon}
			\begin{cases}
				g_\beta(y), & |y|\le\eta, \\
				r_{1}\left(y\right), & \eta\le y\le2\eta, \\
				r_{2}\left(y\right), & -2\eta\le y\le\eta, \\
				0, & |y|\ge2\eta,
			\end{cases}
		\]
		with $ r_{1,2} $ positive smooth functions chosen in such a way that $ g_{\mathrm{trial}} $ is at least $ C^2 $. We also assume that both functions $ r_{1,2} $ are also monotonically decreasing. The normalization constant $ c_{\varepsilon}$, which ensures that $ \left\| g_{\mathrm{trial}} \right\| _{L^2(\mathbb{R})} = 1 $, can be easily estimated: assuming that $ \eta_0 > 2 $, we have
		\begin{equation}
			c_{\varepsilon} = 1 + \mathcal{O}(\varepsilon^4),
		\end{equation}
		since, e.g.,
		\begin{displaymath}
			\int_{\eta}^{2\eta} \mathrm{d} y \: r_1^2(y) \leq \eta g_{\star}^2(\eta) = \mathcal{O}(\eta\varepsilon^{2\eta_0}).
		\end{displaymath}
		Notice also that we need to use \eqref{eq: gcrm almost symm} to reconstruct the norm of $ g_{\star} $:
		\begin{displaymath}
			\int_{-\eta}^{\eta} \mathrm{d} y \: g_{\mathrm{trial}}^2 = c_{\varepsilon}^2 \int_{-\eta}^{\eta} \mathrm{d} y \: (1 + \varepsilon^2 y) \: g_{\star}^2 + c_{\varepsilon}^2 \varepsilon^2 \int_{-\eta}^{\eta} \mathrm{d} y \: y \: g_{\star}^2 = c_{\varepsilon}^2 + \mathcal{O}(\varepsilon^4).
		\end{displaymath}
		Now we estimate
		\[
			E^{\mathrm{gv}}
			\le\mathcal{E}^{\mathrm{gv}}[g_{\mathrm{trial}}]
			=c_{\varepsilon} ^2
			\int _{-\eta} ^\eta
			\mathrm{d} y\
			\left\{
				\textstyle\frac{1}{2} (
				g_\beta '
			) ^2
			+\textstyle\frac{\alpha ^2}{2} y ^2 g_\beta ^2
			+\textstyle\frac{1}{2\pi} g_\beta ^4
			\right\}
			+\mathcal{O}(\varepsilon^4).
		\]
		Thanks to \eqref{eq: gcrm almost symm} we can easily estimate the error we make by replacing $ \frac{\alpha ^2}{2} y ^2 $ with $ U _\beta + \varepsilon ^2 y ^3 v(y) $: denoting for short $ \langle f \rangle : = \mean{g_{\star}}{f}{g_{\star}} $, we have
		\bmln{
			\left\langle
			g_{\star}\left|
				\textstyle\frac{\alpha ^2}{2}y ^2
				-U_\beta(y)
				-\varepsilon ^2 y ^3v(y)
			\right|g_{\star}
			\right\rangle_{L^2(-\eta,\eta)}		\\
			=\left\langle
			g_{\star}\left|
				\textstyle\frac{\alpha ^2 \varepsilon ^2 y ^3\left(
					2
					+\varepsilon ^2 y
				\right)}{2\lf(1+\eps^2y\ri) ^2}
				+\frac{\varepsilon ^2 {\beta_{\star}}}{\lf(1+\eps^2y\ri) ^2}\left(
					2\Omega _0 y
					+\Omega _0 \varepsilon ^2 y
					-\frac{1}{2}\varepsilon ^2 {\beta_{\star}}
				\right)
				-\varepsilon ^2 y ^3 v(y)
			\right|g_{\star}
			\right\rangle
			=\bigO{
			\varepsilon ^4
			+\varepsilon ^2 \left\langle	y ^3 \right\rangle
		},}
		so that
		\begin{equation}
			\label{eq: gve ub}
			E^{\mathrm{gv}}
			\le
			E_{\star}^\mathrm{gv}
			+\bigO{
			\varepsilon ^4
			+\varepsilon ^2 \left\langle	y ^3 \right\rangle
		},
		\end{equation}
		where we have also used \eqref{eq: en preliminary est}, which in turn requires $ {\beta_{\star}} = \mathcal{O}(\varepsilon^{-2}) $.

		The trial state for the functional $ \mathcal{E}^{\mathrm{gv}}_{{\beta_{\star}}}$ is simply the truncation of $ g_{\mathrm{gv}} $, i.e., $ c_{\varepsilon} g_{\mathrm{gv}} $, where now the normalization factor can be estimated in this case as
		\begin{equation}
			c_{\varepsilon} = 1 + \mathcal{O}(\varepsilon^\infty),
		\end{equation}
		since by symmetry of $ g_{\mathrm{gv}} $
		\begin{displaymath}
			\int_{-\eta}^{\eta} \mathrm{d} y \: (1 + \varepsilon^2 y) \: g_{\mathrm{gv}}^2 = 1 - 2 \int_{\eta}^{\infty} \mathrm{d} y \: g_{\mathrm{gv}}^2 = 1 + \mathcal{O}(\varepsilon^{\infty}),
		\end{displaymath}
		where we have used the pointwise estimate \eqref{eq: estimate on gvm} on $ g_{\mathrm{gv}} $ and the fact that the integral of a gaussian, i.e., the error function, is bounded by the value of the gaussian at the boundary, i.e., $ \exp\{- c |\log\varepsilon|^2 \} = \mathcal{O}(\varepsilon^{\infty}) $ (see, e.g., \cite[Eq. (5.1.19)]{AS}). Then we have
		\bml{
			\label{eq: gcre ub}
			E_{\star}^\mathrm{gv}
			\le \mathcal{E}^{\mathrm{gv}}_{{\beta_{\star}}}[c_{\varepsilon} g_{\mathrm{gv}}]
			= (1 + \mathcal{O}(\varepsilon^{\infty})) \int _{-\eta} ^\eta
			\diff y \: \xinyp \:
			\left\{
			\textstyle\frac{1}{2} (g_{\mathrm{gv}}^{\prime})^2+U _{\beta_{\star}}(y) g_{\mathrm{gv}}^2
			+\varepsilon ^2 y ^3 v(y) g_{\mathrm{gv}}^2
			+\frac{1}{2\pi}g_{\mathrm{gv}}^4
			\right\}
			\\
			=
			\int _{-\eta} ^\eta
			\mathrm{d} y
			\left\{
			\textstyle\frac{1}{2} (g_{\mathrm{gv}}^{\prime})^2+ \textstyle\frac{\alpha^2}{2} g_{\mathrm{gv}}^2
			+\frac{1}{2\pi}g_{\mathrm{gv}}^4
			\right\}
			+ \mathcal{O}(\varepsilon^4)
			= E^{\mathrm{gv}}
			+ \mathcal{O}(\varepsilon^4).
		}
		Putting together \eqref{eq: gve ub} with \eqref{eq: gcre ub}, we obtain
		\begin{equation}
			\label{eq: en diff est 1}
			E_{\star}^\mathrm{gv}
			= E^{\mathrm{gv}}
			+\bigO{
			\varepsilon ^4
			+\varepsilon ^2 \left\langle	y ^3 \right\rangle
			}.
		\end{equation}

		Now we decouple the energy $ E^{\mathrm{gv}} $: first we bound from below $ E^{\mathrm{gv}} $ as
		\begin{displaymath}
			E^{\mathrm{gv}} \geq \int_{-\eta}^{\eta} \mathrm{d} y \: (1 + \varepsilon^2 y) \left\{ \textstyle\frac{1}{2} (g_{\mathrm{gv}}^{\prime})^2+ \textstyle\frac{\alpha^2}{2} y^2 g_{\mathrm{gv}}^2
			+\frac{1}{2\pi}g_{\mathrm{gv}}^4 \right\},
		\end{displaymath}
		where we have just dropped some positive quantities and used the symmetry of $ g_{\mathrm{gv}} $. Then we set $ g_{\mathrm{gv}} = u g_{\star} $ for some unknown smooth function $ u $ (recall that $ g_{\star} $ never vanishes in $ [-\eta,\eta] $) and using the variational equation for $ g_{\star} $ as well as Neumann boundary conditions, we obtain
		\bmln{
			E^{\mathrm{gv}} \geq E_{\star}^\mathrm{gv} + \int_{-\eta}^{\eta} \mathrm{d} y \: (1 + \varepsilon^2 y) \: g_{\star}^2 \left\{\textstyle\frac{1}{2} (u')^2 + \left( \frac{\alpha^2}{2} y^2 - U_{{\beta_{\star}}}(y) - \varepsilon^2 y^3 v(y) \right) u^2 + \frac{1}{2\pi} g_{\star}^2 \left( 1 - u^2\right)^2 \right\} 	\\
		- \textstyle\frac{1}{2}\varepsilon^2 \displaystyle\int_{-\eta}^{\eta} \mathrm{d} y \: u^2 g_{\star} g_{\star}' + \mathcal{O}(\varepsilon^{\infty}).
		}
		Then we estimate
		\bmln{
			 \left|\int_{-\eta}^{\eta} \mathrm{d} y \: (1 + \varepsilon^2 y) \: \left( \textstyle\frac{\alpha^2}{2} y^2 - U_{{\beta_{\star}}}(y) - \varepsilon^2 y^3 v(y) \right) g_{\star}^2 u^2 \right| \leq C \varepsilon^2 \int_{-\eta}^{\eta} \mathrm{d} y \: |y| g_{\star}^2 \left|1 - u^2 \right| \\
			 +	 \left| \int_{-\eta}^{\eta} \mathrm{d} y \: (1 + \varepsilon^2 y) \: \left( \textstyle\frac{\alpha^2}{2} y^2 - U_{{\beta_{\star}}}(y) - \varepsilon^2 y^3 v(y) \right) g_{\star}^2 \right| \\
			 \leq C \varepsilon^2 \eta^{3/2} \left\| g_{\star}^2(1 - u^2) \right\|_{L^2_{\eta}} +  \mathcal{O}\left(\varepsilon^2 \left\langle y^3 \right\rangle +\varepsilon^4\right),
		}
		and by \eqref{eq: derivative gcrm} (notice that the factor $1 + \varepsilon^2 y $ is uniformly bounded from above and below by a constant)
		\bmln{
			\left| \displaystyle\int_{-\eta}^{\eta} \mathrm{d} y \: u^2 g_{\star} g_{\star}' \right| = \displaystyle\int_{-\eta}^{\eta} \mathrm{d} y \: \left|1 - u^2 \right| |g_{\star}| |g_{\star}'| + \mathcal{O}(\varepsilon^4) \leq C \eta^3 \displaystyle\int_{-\eta}^{\eta} \mathrm{d} y \: \left|1 - u^2 \right| g_{\star}^2 + \mathcal{O}(\varepsilon^4)	\\
			 \leq C \eta^{7/2} \left\| g_{\star}^2(1 - u^2) \right\|_{L^2_{\eta}} + \mathcal{O}(\varepsilon^4),
		}
		which imply by dropping the kinetic term
		\bml{
			\label{eq: en diff est 2}
			E^{\mathrm{gv}} \geq E_{\star}^\mathrm{gv} + \textstyle\frac{1}{2} \left\| g_{\star} u^{\prime} \right\|_{L^2_{\eta}}^2 + \textstyle\frac{1}{\pi} \left\| g_{\star}^2(1 - u^2) \right\|^2_{L^2_{\eta}}- C \varepsilon^2 \eta^{7/2} \left\| g_{\star}^2(1 - u^2) \right\|_{L^2_{\eta}} +  \mathcal{O}\left(\varepsilon^2 \left\langle y^3 \right\rangle +\varepsilon^4\right)	\\
			\geq E_{\star}^\mathrm{gv} + \textstyle\frac{1}{\pi} \left(\left\| g_{\star}^2(1 - u^2) \right\|_{L^2_{\eta}}- C \varepsilon^2 \eta^{7/2} \right)^2 +  \mathcal{O}\left(\varepsilon^2 \left\langle y^3 \right\rangle +\varepsilon^4 \eta^7\right).
		}
		If we compare what we have obtained with \eqref{eq: en diff est 1}, we conclude that
		\begin{equation}
			\left\| g_{\star}^2 - g_{\mathrm{gv}}^2\right\|^2_{L^2_{\eta}} = \left\| g_{\star}^2(1 - u^2) \right\|_{L^2_{\eta}}^2 =  \mathcal{O}\left(\varepsilon^2 \left\langle y^3 \right\rangle +\varepsilon^4 \eta^7\right),
		\end{equation}
		but on the other hand
		\begin{equation}
			\label{eq: y3 proof}
			\left| \int_{-\eta}^{\eta} \mathrm{d} y \: y^3 \: g_{\star}^2 \right| = \left| \int_{-\eta}^{\eta} \mathrm{d} y \: y^3 \: \left( g_{\star}^2 - g_{\mathrm{gv}}^2 \right) \right| \leq C \eta^{7/2} \left\| g_{\star}^2 - g_{\mathrm{gv}}^2\right\|_{L^2_{\eta}},
		\end{equation}
		so that finally
		\begin{displaymath}
			\left\| g_{\star}^2 - g_{\mathrm{gv}}^2\right\|^2_{L^2_{\eta}} = \mathcal{O}\left(\varepsilon^4 \eta^7\right).
		\end{displaymath}
		This proves the second inequality in \eqref{eq: gv energy diff} but the first is obtained by replacing the above estimate into \eqref{eq: en diff est 1}.
	\end{proof}

	\noindent
	It is interesting to remark that a by-product of the proof of Proposition \ref{pro: gv energy diff} is that (see \eqref{eq: y3 proof})
	\begin{equation}
		\label{eq: y3}
		\mean{g_{\star}}{y^3}{g_{\star}} = \mathcal{O}(\varepsilon^2 \eta^7),
	\end{equation}
	which in combination with \eqref{eq: gcrm almost symm} is a strong indication of $ g_{\star} $ being symmetric w.r.t. the origin with a very high precision. In fact this is also made apparent by the estimate of the difference $ g_{\star} - g_{\mathrm{gv}} $.

	The $L^2$-statement in \eqref{eq: gv energy diff} can indeed be improved to an $ L^{\infty}$-one, showing that $ g_{\star} $ and $ g_{\mathrm{gv}} $ are pointwise close. The price to pay to have a result in a stronger norm is the restriction of the region under consideration to the annulus $ \widetilde{\mathcal{A}}_{\eta} \subset {\mathcal{A}_{\eta}} $ defined as
	\begin{equation}
		\widetilde{\mathcal{A}}_{\eta} :
		=\left\{
			y \in \mathbb{R}\;\big|\;
			g_{\mathrm{gv}}(y) \geq
			 \textstyle\frac{1}{\eta ^\nu}
		\right\},
	\end{equation}
	for some $ \nu > 0 $ independent of $ \varepsilon $. Note that thanks to the pointwise estimates \eqref{eq: estimate on gvm} and the monotonicity of $ g_{\mathrm{gv}} $,
	\begin{equation}
		\widetilde{\mathcal{A}}_{\eta} = [-y_{\eta},y_{\eta}],	\qquad		\mbox{with } y_{\eta} \gg 1.
	\end{equation}

	\begin{pro}[Pointwise estimate of $ g_{\star} - g_{\mathrm{gv}} $]
		\label{pro: infty est diff gvm}
		\mbox{}	\\
		As $ \varepsilon \to 0 $ and for any $ \nu > 0 $
		\begin{equation}
			\label{eq: infty est diff gvm}
			\left\|
				g_{\star}
				-g_{\mathrm{gv}}
				\right\| _{L ^\infty \left(\widetilde{\mathcal{A}}_{\eta}\right)}
			=\bigO{\varepsilon ^2 \eta ^{\frac{7+4\nu}{2}}}.
		\end{equation}
	\end{pro}

	\begin{proof}
		Going back to \eqref{eq: en diff est 2} and retaining the kinetic term, we see that we obtain via \eqref{eq: gv energy diff} the upper bound
		\begin{equation}
			\left\| g_{\star} u^{\prime} \right\|_{L^2_{\eta}}^2 = \mathcal{O}(\varepsilon^4 \eta^7),
		\end{equation}
		where we recall that $ u = g_{\mathrm{gv}}/g_{\star} $. Now let us introduce the set
		\begin{displaymath}
			\widehat{\mathcal{A}}_{\eta} : = \left\{
			y \in \mathbb{R}\;\big|\;
			g_{\star}(y) \geq
			 \textstyle\frac{1}{2\eta ^\nu}
		\right\},
		\end{displaymath}
		so that above inequality  together with \eqref{eq: gv energy diff} and the bound $ |1 - u| \leq |1 - u^2| $ imply
		\begin{displaymath}
			\left\| u^{\prime} \right\|_{L^2(\widehat{\mathcal{A}}_{\eta})}^2 = \mathcal{O}(\varepsilon^4 \eta^{7+2\nu}),	\qquad \left\| 1 - u \right\|^2_{L^2(\widehat{\mathcal{A}}_{\eta})} = \mathcal{O}(\varepsilon^4 \eta^{7+4\nu}).
		\end{displaymath}
		Then it suffices to use Sobolev inequality in one-dimension:
		\begin{equation}
			\left\| 1 - u \right\|_{L^{\infty}(\widehat{\mathcal{A}}_{\eta})}^2 \leq C \left( \left\| u' \right\|_{L^{2}(\widehat{\mathcal{A}}_{\eta})}^2 + \left\| 1 - u \right\|_{L^{2}(\widehat{\mathcal{A}}_{\eta})}^2 \right) = \mathcal{O}\left(\varepsilon^4 \eta^{7+4\nu} \right).
		\end{equation}
		Finally to obtain the result it remains to observe that $ \widetilde{\mathcal{A}}_{\eta} \subset \widehat{\mathcal{A}}_{\eta} $, because in the region where $ g_{\star} \geq 1/(2\eta^\nu) $, by the pointwise estimate, $ g_{\mathrm{gv}} $ is larger than $ (1 + o(1))/(2\eta^\nu) $, which is obviously satisfied if $ g_{\mathrm{gv}} \geq 1/\eta^\nu  $.
	\end{proof}

	\noindent
	The above bound shows that inside $ \widetilde{\mathcal{A}}_{\eta} $ one can estimate the distance of $ g_{\star} $ from a perfectly even function: for any $ y \in \widetilde{\mathcal{A}}_{\eta} $
	\begin{displaymath}
		g_{\star}(-y)
		= g_{\star}(y)
		+\bigO{\varepsilon ^2 \eta ^{\frac{7+4\nu}{2}}},
	\end{displaymath}
	which is perfectly compatible with the estimates \eqref{eq: gcrm almost symm} and \eqref{eq: y3}.

	Another useful consequence of the above pointwise statement is the following

	\begin{cor}[Maximum point of $ g_{\star} $]
		\label{cor: max gcrm}
		\mbox{}	\\
		Let $ y_{{\beta_{\star}}} $ be the unique maximum point of $ g_{\star} $, then as $ \varepsilon \to 0 $
		\begin{equation}
			y_{{\beta_{\star}}} = \bigO{\varepsilon ^2 \eta ^{\frac{7+4\nu}{2}}}.
		\end{equation}
	\end{cor}

	\begin{proof}
		The result is a straightforward consequence of the pointwise estimate \eqref{eq: infty est diff gvm} and the properties of $ g_{\mathrm{gv}} $ (see Proposition \ref{pro: min gvf}).
	\end{proof}

\subsection{Critical Velocity and Positivity of the Cost Function}
\label{sec: positivity}

From now we fix the phase to be optimal one, i.e., $ \beta = {\beta_{\star}} $. The {\it potential function} is defined as
\bml{
	\label{eq: potential f}
	F(y):
	=- \frac{1}{\varepsilon^2} \int _{-\eta} ^y
	\mathrm{d} t \: (1 + \varepsilon^2 t)
	\left. \partial _\beta U _\beta (t) \right|_{\beta = {\beta_{\star}}}
	g_{\star} ^2	\\
	=2 \Omega_0 \int_{-\eta}^{y} \mathrm{d} t \: \frac{1}{1 + \varepsilon^2 t} \left(t + \frac{1}{2} \varepsilon^2 t^2 - \frac{\varepsilon^2 {\beta_{\star}}}{2 \Omega_0}  \right) g_{\star}^2.
}
The main object under investigation is the {\it cost function}
\begin{equation}
	\label{eq: cost f}
	K(y) : = \textstyle\frac{1}{2} g_{\star}^2(y) + F(y),
\end{equation}
and our main goal in this Section is to prove that it is positive in the bulk of the condensate when $ \Omega_0 \geq \Omega_{\mathrm{c}} $. To this purpose we will clearly have to investigate the equation \eqref{eq: Omegac} and prove at least that there exists a positive solution to it. Notice the equation \eqref{eq: Omegac} involves only quantities relative to the limiting functional $ \mathcal{E}^{\mathrm{gv}} $ and is independent of $ \varepsilon $. We thus introduce the analogue of \eqref{eq: cost f} for the limiting case, i.e., the function \eqref{eq: kgv}
\begin{displaymath}
	K^{\mathrm{gv}} = \tx\frac{1}{2} g_{\mathrm{gv}}^2 + F^{\mathrm{gv}},
\end{displaymath}
where $ F^{\mathrm{gv}} $ is defined in \eqref{eq: fgv}:
\begin{displaymath}
	F^{\mathrm{gv}}(y) = - 2 \Omega_0 \int_y^{\infty} \mathrm{d} t \: t \: g_{\mathrm{gv}}^2(t).
\end{displaymath}

We will start by studying the positivity of \eqref{eq: kgv} and show that the condition $ \Omega_0 > \Omega_{\mathrm{c}} $, where the latter is defined as the biggest solution to \eqref{eq: Omegac}, is sufficient to deduce that $ K^{\mathrm{gv}}(y) \geq 0 $ for any $ y \in \mathbb{R} $. In the second part of the Section we will turn our attention to the cost function \eqref{eq: cost f} and prove that the same condition on $ \Omega_0 $ guarantees positivity of $ K $ as well.

We first observe that $ F^{\mathrm{gv}} $ is a negative function vanishing at $ \pm \infty $: at $ y = +\infty $ it is obvious, at $ -\infty $ it is a consequence of parity of $ g_{\mathrm{gv}} $. By this property one can rewrite
\begin{displaymath}
	F^{\mathrm{gv}}(y) = 2 \Omega_0 \int_{-\infty}^y \mathrm{d} t\: t \: g_{\mathrm{gv}}^2(t).
\end{displaymath}
In fact there is another explicit expression of $ F^{\mathrm{gv}} $, which can be obtained by using the variational equation \eqref{eq: gvm variational}:
\bml{
	\label{eq: fgv alternative}
	F^{\mathrm{gv}}(y) =  - \Omega_0 \int_y^{\infty} \mathrm{d} t \: \partial_t (t^2) \: g_{\mathrm{gv}}^2 = \Omega_0 y^2 g_{\mathrm{gv}}^2(y) + 2 \Omega_0 \int_y^{\infty} \mathrm{d} t \: t^2 \: g_{\mathrm{gv}} g_{\mathrm{gv}}^{\prime}	\\
	 = \Omega_0 y^2 g_{\mathrm{gv}}^2(y) + \frac{4 \Omega_0}{\alpha^2} \int_y^{\infty} \mathrm{d} t \: g_{\mathrm{gv}}^{\prime} \left[ \textstyle\frac{1}{2} g_{\mathrm{gv}}^{\prime\prime} - \frac{1}{\pi} g_{\mathrm{gv}}^3 + \mu^{\mathrm{gv}} g_{\mathrm{gv}} \right] 	\\
	 =
	- \textstyle\frac{1}{\Omega_0(s+2)} \left( g_{\mathrm{gv}}^{\prime}(y) \right)^2 + \left[ \Omega_0 y^2 + \frac{1}{\pi \Omega_0 (s+2)} g_{\mathrm{gv}}^2(y) - \frac{2\mu^{\mathrm{gv}}}{\Omega_0(s+2)} \right] g_{\mathrm{gv}}^2(y),
	}
where we have used the exponential decay at $ \infty $ of $ g_{\mathrm{gv}} $ to cancel the missing boundary terms. Consequently we can rewrite $ K^{\mathrm{gv}} $ as
\begin{equation}
	\label{eq: kgv alternative}
	K^{\mathrm{gv}}(y) = - \textstyle\frac{1}{\Omega_0(s+2)} \left( g_{\mathrm{gv}}^{\prime}(y) \right)^2 + \left[ \textstyle\frac{1}{2} + \Omega_0 y^2 + \frac{1}{\pi \Omega_0 (s+2)} g_{\mathrm{gv}}^2(y) - \frac{2\mu^{\mathrm{gv}}}{\Omega_0(s+2)} \right] g_{\mathrm{gv}}^2(y).
\end{equation}
The main result about $ K^{\mathrm{gv}} $ is the following

	\begin{pro}[Positivity of $ K^{\mathrm{gv}} $]
		\label{pro: kgv positive}
		\mbox{}	\\
		Let $ \Omega_0 > 0 $, then
		\begin{equation}
			\label{eq: kgv positive}
			K^{\mathrm{gv}}(y) \geq 0 \mbox{ for any } y \in \mathbb{R} 	\Longleftrightarrow \Omega_0 \geq \textstyle\frac{4}{s+2} \left[ \mu^{\mathrm{gv}} - \frac{1}{2\pi} g_{\mathrm{gv}}^2(0) \right].
		\end{equation}
		Moreover if the strict inequality is verified on the r.h.s., $ K^{\mathrm{gv}}(y) > 0 $ for any $ y $ finite.
	\end{pro}

	{
	\begin{rem}[Comparison with \cite{CPRY3}]
		\label{rem: cost functions}
		\mbox{}	\\
		Despite the use of two different potential functions $ F_1 $ and $ F_2 $ in \cite{CPRY3}, one should realize that \cite[Lemma 3.3]{CPRY3} yields the poitwise positivity of a cost function, which is the analogue of $ K^{\mathrm{gv}} $ \underline{in the asymptotic regime} $ \Omega_0 \gg 1 $. In fact it can be easily seen that the cost function in \cite{CPRY3} is bounded from below by $ K^{\mathrm{gv}} $ and therefore the positivity of the latter implies the positivity of the first. Hence any threshold $ \bar\Omega_0 $ one might deduce there must be larger than $ \Omega_{\mathrm{c}} $ by definition.
	\end{rem}}

	\begin{proof}
		One side of the statement, i.e., the fact that the condition $ K^{\mathrm{gv}}(0) \geq 0 $ is necessary for the positivity of $ K^{\mathrm{gv}} $ everywhere, is obviously trivial, so we focus on the other side of the implication, namely that $ K^{\mathrm{gv}}(0) \geq 0 $ is also sufficient.

		 The core of the proof is to show that
		 \begin{equation}
		 	\label{eq: condition kgv}
		 	K^{\mathrm{gv}}(y) \geq 0 \mbox{ for any } y \in \mathbb{R} 	\Longleftrightarrow K^{\mathrm{gv}}(0) \geq 0.
		\end{equation}
		Indeed if we assume that this double implication is true, a straightforward computation yields
		\begin{equation}
			\label{eq: kgv0}
			K^{\mathrm{gv}}(0) = \left[ \frac{1}{2} + \frac{1}{\pi \Omega_0 (s+2)} g_{\mathrm{gv}}^2(0) - \frac{2\mu^{\mathrm{gv}}}{\Omega_0(s+2)} \right] g_{\mathrm{gv}}^2(0),
		\end{equation}
		since $ g_{\mathrm{gv}} $ is symmetric w.r.t. the origin and has a maximum at $ y = 0  $ (see Proposition \ref{pro: min gvf}). The result is then a trivial consequence of strict positivity of $ g_{\mathrm{gv}}(0) $.

		In order to prove \eqref{eq: condition kgv}, we first observe that $ K^{\mathrm{gv}}(\pm \infty) = 0 $ and $ K^{\mathrm{gv}} $ is smooth, so, if there was a point $ y_0 $ where $ K^{\mathrm{gv}} $ becomes negative, it must be $ |y_0| < +\infty $. Moreover as $ g_{\mathrm{gv}} $, $ K^{\mathrm{gv}} $ is symmetric w.r.t. to the origin, so it suffices to consider $ y \in \mathbb{R}^+ $. The derivative of $ K^{\mathrm{gv}} $ is easily computed from the expression \eqref{eq: kgv}:
		\begin{equation}
			{K^{\mathrm{gv}}}^{\prime}(y) = g_{\mathrm{gv}} g_{\mathrm{gv}}^{\prime} + 2 \Omega_0 y g_{\mathrm{gv}}^2,
		\end{equation}
		and one immediately has that $ {K^{\mathrm{gv}}}^{\prime}(0) = 0 $, i.e., $ K^{\mathrm{gv}} $ has a critical point there. Whether it is a minimum or a maximum depends on $ s $ and $ \Omega_0 $, but as we are going to see this does not matter. We can in any case compute easily the second derivative of $ K^{\mathrm{gv}} $ exploiting once more \eqref{eq: gvm variational}:
		\begin{equation}
			\label{eq: kgv second derivative}
			{K^{\mathrm{gv}}}^{\prime\prime}(y) = \left(g_{\mathrm{gv}}^{\prime}\right)^2 + 4 \Omega_0 y g_{\mathrm{gv}} g_{\mathrm{gv}}^{\prime} + 2 \left[ \textstyle\frac{1}{2} \alpha^2 y^2 + \Omega_0  - \mu^{\mathrm{gv}}+ \textstyle\frac{1}{\pi} g_{\mathrm{gv}}^2 \right] g_{\mathrm{gv}}^2.
		\end{equation}

		Then we prove the crucial property of $ K^{\mathrm{gv}} $: suppose that $ K^{\mathrm{gv}} $ has a maximum at $ y_1  \geq 0 $ and then a minimum at $ y_2 > y_1 $, then
		\begin{equation}
			\label{eq: kgv crucial pro}
			\frac{K^{\mathrm{gv}}(y_2)}{g_{\mathrm{gv}}^2(y_2)} \geq \frac{K^{\mathrm{gv}}(y_1)}{g_{\mathrm{gv}}^2(y_1)},
		\end{equation}
		and in particular $ K^{\mathrm{gv}}(y_2) \geq 0 $ if $ K^{\mathrm{gv}}(y_1) \geq 0 $.

		To conclude the argument once \eqref{eq: kgv crucial pro} is proven, it is sufficient to observe that $ K^{\mathrm{gv}} $ has a critical point in $ y = 0 $, which by parity must be either a maximum or a minimum: if it is a maximum, then \eqref{eq: kgv crucial pro} shows that at any minimum point $ y_2 > 0 $, $ K^{\mathrm{gv}}(y_2) \geq 0 $. Notice that it does not matter whether $ K^{\mathrm{gv}} $ has a single or multiple minima, because any minimum after the first requires the presence of a preceding maximum point, where $ K^{\mathrm{gv}} $ is larger than its first minimum and therefore positive. If on the opposite $ K^{\mathrm{gv}} $ has a minimum at the origin, then it means that there must be a maximum at some $ y_1 > 0 $, where obviously $ K^{\mathrm{gv}}(y_1) \geq K^{\mathrm{gv}}(0) \geq 0 $ and we can repeat the argument for any minimum after $ y_1 $.

		Let us now prove \eqref{eq: kgv crucial pro}: we assume again that $ K^{\mathrm{gv}} $ has a maximum in $ y_1 \geq 0 $ and a minimum in $ y_2 > y_1 $. Then it must be $ {K^{\mathrm{gv}}}^{\prime}(y_{1,2}) = 0 $, i.e.,
		\begin{equation}
			g_{\mathrm{gv}}^{\prime}(y_{1,2}) = - 2 \Omega_0 y_{1,2} g_{\mathrm{gv}}^2(y_{1,2}).
		\end{equation}
		Moreover replacing this condition in \eqref{eq: kgv} and \eqref{eq: kgv second derivative}, we get
		\begin{equation}
			\label{eq: kgv min}
			K^{\mathrm{gv}}(y_{1,2}) = \left[ \Omega_0 \textstyle\frac{s-2}{s+2} y_{1,2}^2 + \frac{1}{2}  - \frac{2}{\Omega_0(s+2)} \left( \mu^{\mathrm{gv}} - \frac{1}{2\pi} g_{\mathrm{gv}}^2(y_{1,2}) \right) \right] g_{\mathrm{gv}}^2(y_{1,2}),
		\end{equation}
		\begin{equation}
			\label{eq: kgv second derivative min}
			{K^{\mathrm{gv}}}^{\prime\prime}(y_{1,2}) = \left[ \Omega_0^2 (s-2) y_{1,2}^2 + 2\Omega_0 - 2 \mu^{\mathrm{gv}} + \textstyle\frac{2}{\pi} g_{\mathrm{gv}}^2(y_{1,2}) \right] g_{\mathrm{gv}}^2(y_{1,2}).
		\end{equation}
		Moreover
		\begin{equation}
			{K^{\mathrm{gv}}}^{\prime\prime}(y_1) \leq 0 \leq {K^{\mathrm{gv}}}^{\prime\prime}(y_2),
		\end{equation}
		which implies
		\begin{equation}
			\Omega_0^2(s-2) y_1^2 + 2 \Omega_0 - 2 \mu + \textstyle\frac{2}{\pi} g_{\mathrm{gv}}^2(y_{1}) \leq \Omega_0^2(s-2) y_2^2 + 2 \Omega_0 - 2 \mu + \textstyle\frac{2}{\pi} g_{\mathrm{gv}}^2(y_{2}),
		\end{equation}
		and using this inequality in the expression of $ K^{\mathrm{gv}}(y_2) $, we obtain
		\bml{
			\frac{K^{\mathrm{gv}}(y_2)}{g_{\mathrm{gv}}^2(y_2)} \geq \Omega_0 \frac{s-2}{s+2} y_{1}^2 - \frac{2}{\pi\Omega_0(s+2)} \left( g_{\mathrm{gv}}^2(y_2) - g_{\mathrm{gv}}^2(y_1) \right) + \frac{1}{2}  - \frac{2}{\Omega_0(s+2)} \left( \mu^{\mathrm{gv}} - \frac{1}{2\pi} g_{\mathrm{gv}}^2(y_{2}) \right) \\
			= \frac{K^{\mathrm{gv}}(y_1)}{g_{\mathrm{gv}}^2(y_1)} +  \frac{1}{\pi\Omega_0(s+2)} \left( g_{\mathrm{gv}}^2(y_1) - g_{\mathrm{gv}}^2(y_2) \right) \geq \frac{K^{\mathrm{gv}}(y_1)}{g_{\mathrm{gv}}^2(y_1)},
		}
		because by hypothesis $ y_2 > y_1 $ and $ g_{\mathrm{gv}} $ is decreasing.

		Notice that as a by-product of our analysis we found out that can have no global minima, since $ \inf_{y \in \mathbb{R}} K^{\mathrm{gv}}(y) = 0 $ and therefore at any such minimum point $ y_0 $, we would have $ K^{\mathrm{gv}}(y_0) = 0 $, but this clearly contradicts \eqref{eq: kgv crucial pro}. Hence if the inequality on r.h.s. of \eqref{eq: kgv positive} is strict then $ K^{\mathrm{gv}}(y) > 0 $ for any finite $ y $.
	\end{proof}

\noindent
Proposition \ref{pro: kgv positive} introduces the equation \eqref{eq: Omegac}. The next step is obviously to prove that such an equation as at least one solution:

	\begin{pro}[Equation \eqref{eq: Omegac}]
		\label{pro: Omegac}
		\mbox{}	\\
		The equation \eqref{eq: Omegac}
		\begin{displaymath}
			 \Omega_0 = \frac{4}{s+2} \left[ \mu^{\mathrm{gv}} - \frac{1}{2\pi} g_{\mathrm{gv}}^2(0) \right].
		\end{displaymath}
		has at least one solution $ \Omega_0 > 0 $.
	\end{pro}

	\begin{proof}
		Let us first set
		\begin{displaymath}
			G(\Omega_0) : = \frac{4}{s+2} \left[ \mu^{\mathrm{gv}} - \frac{1}{2\pi} g_{\mathrm{gv}}^2(0) \right],
		\end{displaymath}
		so that \eqref{eq: Omegac} reads $ \Omega_0 = G(\Omega_0) $.	We will show that $ G(\Omega_0) $ is asymptotically smaller than $ \Omega_0 $ (resp. larger) $ \Omega_0 $ for large (resp. small) $ \Omega_0 $.

		Let us first consider $ \Omega_0 \gg 1 $: using the trivial bound $ \left\| g_{\mathrm{gv}} \right\|_4^4 \leq \left\| g_{\mathrm{gv}} \right\|_{\infty}^2 $ and the definition of $ \mu^{\mathrm{gv}} $, we get
		\begin{displaymath}
			G(\Omega_0) \leq \textstyle\frac{4}{s+2} E^{\mathrm{gv}}.
		\end{displaymath}
		If now we plug into $ \mathcal{E}^{\mathrm{gv}} $ as a trial state the ground state of the harmonic oscillator $ h_{\mathrm{osc}} = - \frac{1}{2} \Delta + \frac{1}{2} \alpha^2 y^2 $, we easily obtain
		\begin{displaymath}
			E^{\mathrm{gv}} \leq \textstyle\frac{1}{2} \Omega_0 \sqrt{s+2} \left(1 + \mathcal{O}(\Omega_0^{-1/2}) \right),
		\end{displaymath}
		so that
		\begin{displaymath}
			G(\Omega_0) \leq \textstyle\frac{2}{\sqrt{s+2}} \Omega_0  \left(1 + \mathcal{O}(\Omega_0^{-1/2}) \right) < \Omega_0,
		\end{displaymath}
		if $ \Omega_0 \gg 1 $, because $ s > 2 $.

		On the other hand for small $ \Omega_0 $, thanks to the estimate \eqref{eq: gvm infty bound} and again the definition of $ \mu^{\mathrm{gv}} $, we have
		\begin{displaymath}
			G(\Omega_0) \geq \textstyle\frac{2}{s+2} E^{\mathrm{gv}}.
		\end{displaymath}
		To bound from below $ E^{\mathrm{gv}} $ for small $ \Omega_0 $, we can simply drop the kinetic term to get
		\begin{displaymath}
			E^{\mathrm{gv}} \geq \inf_{\left\| \rho \right\|_1 = 1} \int_{\mathbb{R}} \mathrm{d} y \: \left\{ \textstyle\frac{1}{2} \alpha^2 y^2 \rho + \frac{1}{2\pi} \rho^2 \right\},
		\end{displaymath}
		i.e., a TF-like functional. By scaling we immediately obtain
		\begin{displaymath}
			\inf_{\left\| \rho \right\|_1 = 1} \int_{\mathbb{R}} \mathrm{d} y \: \left\{ \textstyle\frac{1}{2} \alpha^2 y^2 \rho + \frac{1}{2\pi} \rho^2 \right\} = C \Omega_0^{2/3},
		\end{displaymath}
		as $ \Omega_0 \to 0 $, so that
		\begin{displaymath}
			G(\Omega_0) \geq C \Omega_0^{2/3} > \Omega_0,
		\end{displaymath}
		for $ \Omega_0 $ small enough.
	\end{proof}

	In the above Proposition we have not investigated the uniqueness of the solution. We indeed expect that such a solution is in fact unique, but without a proof of this fact, we have to choose $ \Omega_{\mathrm{c}} $ equal to the {\it largest} possible solution.

	Now we turn our attention back to the cost function $ K $: as $ K^{\mathrm{gv}} $ it is given by the sum of a the positive density $ \frac{1}{2} g_{\star}^2 $ and the negative potential function $ F $ (see next Proposition \ref{pro: pro F}). In addition $ g_{\star} $ is monotonically decreasing for $ y \geq y_{{\beta_{\star}}} = o(1) $ and, like $ F $, almost symmetric. Close the boundary of the interval $ [-\eta,\eta] $, $ g_{\star} $ gets extremely small (in fact exponentially small in $ \varepsilon $) but $ F $ vanishes identically at $ \pm \eta $. In conclusion it is clear that the overall positivity of $ K $ should then emerge from a very delicate balance between the two opposite contributions.

	We first state some simple properties of $ F $ collected in the following

	\begin{pro}[Properties of $ F $]
		\label{pro: pro F}
		\mbox{}	\\
		The potential function defined in \eqref{eq: potential f} is such that
		\begin{equation}
			\label{eq: negativity F}
			F(y) \leq 0,	\qquad 	\mbox{for any } y \in [-\eta,\eta],
		\end{equation}
		\begin{equation}
			\label{eq: boundary F}
			F(\pm \eta) = 0.
		\end{equation}
	\end{pro}

	\begin{proof}
		One of the identities \eqref{eq: boundary F} is trivial, the other is a direct consequence of \eqref{eq: optimality}. In order to show that $ F $ is negative everywhere we compute the derivative
		\begin{equation}
			\label{eq: derivative F}
			F'(y)
			= \textstyle\frac{1}{1+\eps^2y}
			\left(
			2\Omega _0 y
			+\Omega _0 \varepsilon ^2 y ^2
			- {\beta_{\star}} \varepsilon ^2
			\right)g_{\star} ^2(y),
		\end{equation}
		and it is easy to verify that because of the first term $ F'(-\eta) < 0 $ while $ F'(\eta) > 0 $. Moreover $ F' $ vanishes at a single point $ y_F = \mathcal{O}(\varepsilon^2) $, where $ F $ has a global minimum. Hence it is negative everywhere in $ [-\eta,\eta] $.
	\end{proof}

	A very crucial piece of information about the potential function formulated in the next Proposition is an alternative expression of it, which relies on the variational equation \eqref{eq: gvbm variational} and is the analogue of \eqref{eq: fgv alternative} for $ F^{\mathrm{gv}} $.

	\begin{pro}[Alternative expressions of $ F $]
		\label{pro: alt expression F}
		\mbox{}	\\
		For any $ y \in [-\eta,\eta] $ the potential function $ F $ admits the following alternative expressions
		\begin{equation}
			F(y)
			= - \frac{\Omega _0}{\alpha ^2} \left(
				g_{\star}'
			\right) ^2
			+ \frac{2\Omega _0}{\alpha ^2}\left[ \frac{1}{2} \alpha^2 y^2
			+\frac{1}{2\pi} g_{\star}^2
			-\mu_{\star}
			\right] g_{\star} ^2 +
			\begin{cases}
				R_{+}(y) + R_{+},	&	\mbox{if } y \geq y_{\beta_{\star}},	\\
				R_{-}(y) + R_{-},	&	\mbox{if } y \leq y_{\beta_{\star}},
			\end{cases}
		\end{equation}
		where
		\begin{equation}
			\label{eq: tails}
			R_{\pm}(y) = \mathcal{O}(\varepsilon ^2 \eta ^7) g_{\star}^2(y), 	\qquad
			R_{\pm} = -\Omega _0 \eta ^2 g_{\star} ^2(\pm\eta)
			\left(
			1
			+ o(1) \right).
		\end{equation}
	\end{pro}

	\begin{proof}
		We consider only the case $ y \ge {y_\beta} $, since the other one is analogous. The key ingredient of the proof is an integration by parts, exactly as for \eqref{eq: fgv alternative}. We spell all the details nevertheless for the sake of clarity. Thanks to the vanishing of $ F $ at $ \eta $, we have
		\bml{
			\label{eq: F computed 0}
			F(y)
			=\int _y ^\eta
			\mathrm{d} t\
			\textstyle\frac{1}{1+\eps^2t}
			\left(
			-2\Omega _0 t
			-\Omega _0 \varepsilon ^2 t ^2
			+{\beta_{\star}} \varepsilon ^2
			\right)
			g_{\star} ^2 = \displaystyle\int _y ^\eta
			\mathrm{d} t\
			\textstyle\frac{1}{1+\eps^2t}
			\: g_{\star} ^2	\: \partial_t \left(
			-\Omega_0 t ^2
			-\frac{1}{3} \Omega _0 \varepsilon ^2 t ^3
			+{\beta_{\star}} \varepsilon ^2 t
			\right)
			\\
			= \textstyle\frac{1}{1+\varepsilon ^2 \eta}
			\left(
			-\Omega _0\eta ^2
			-\frac{1}{3} \Omega _0 \varepsilon ^2 \eta ^3
			+{\beta_{\star}} \varepsilon ^2 \eta
			\right)
			g_{\star} ^2 (\eta)
			-\frac{1}{1+\eps^2y}
			\left(
			-\Omega _0y ^2
			-\frac{1}{3}\Omega _0 \varepsilon ^2 y ^3
			+{\beta_{\star}} \varepsilon ^2 y
			\right)
			g_{\star} ^2 (y)+
			\\
			+\int _y ^\eta
			\mathrm{d} t\
			\left(
			-\Omega _0 t ^2
			-\textstyle\frac{1}{3}\Omega _0 \varepsilon ^2 t ^3
			+{\beta_{\star}} \varepsilon ^2 t
			\right)
			\left[ \textstyle\frac{\varepsilon ^2}{\lf(1+\eps^2y\ri) ^2} g_{\star} ^2 - \frac{2}{\lf(1+\eps^2y\ri)} g_{\star}  g_{\star}' \right]  \\
			=
			-\Omega _0 \eta ^2 g_{\star} ^2 (\eta) (1 + o(1))
			+ \left( \Omega _0 y ^2 + \mathcal{O}(\varepsilon^2\eta^4) \right) g_{\star}^2(y) 	\\
			- 2 \displaystyle\int _y ^\eta
			\mathrm{d} t\: \textstyle\frac{1}{\lf(1+\eps^2y\ri)}
			\left(
			-\Omega _0 t ^2
			-\textstyle\frac{1}{3}\Omega _0 \varepsilon ^2 t ^3
			+{\beta_{\star}} \varepsilon ^2 t
			\right)
			g_{\star}  g_{\star}'.
		}
		We now rewrite the last term by reconstructing the potential $ U_{{\beta_{\star}}} $ and using the variational equation \eqref{eq: gvbm variational}: since
		\bmln{
			-\textstyle\frac{1}{1+\eps^2y}
			\left(
			-\Omega _0 t ^2
			-\frac{1}{3} \Omega _0\varepsilon ^2 t ^3
			+{\beta_{\star}} \varepsilon ^2 t
			\right) = \frac{2\Omega _0}{\alpha ^2} (1 + \varepsilon^2 t) U _{\beta_{\star}}(t) \\
			+ \textstyle\frac{\varepsilon ^2}{\lf(1+\eps^2t\ri)}
			\left(
			\frac{1}{6} \alpha^2 t ^3
			- \frac{1}{2} \Omega_0 (s-2) {\beta_{\star}} t
			+\Omega _0 \varepsilon ^2 {\beta_{\star}} t ^2-\frac{1}{2}\varepsilon ^2 {\beta_{\star}}^2
			\right),
		}
		we obtain
		\bml{
			\label{eq: F computed 1}
			- 2 \int _y ^\eta
			\mathrm{d} t\	\frac{1}{\lf(1+\eps^2y\ri)}
			\left(
			-\Omega _0 t ^2
			-\textstyle\frac{1}{3}\Omega _0 \varepsilon ^2 t ^3
			+{\beta_{\star}} \varepsilon ^2 t
			\right) g_{\star}  g_{\star}' =\frac{4\Omega _0}{\alpha ^2}
			\int _y ^\eta
			\mathrm{d} t \: (1 + \varepsilon^2 t) \:
			U _{\beta_{\star}}(t)
			g_{\star}
			g_{\star}'
			\\
			+  2\varepsilon ^2\int _y ^\eta
			\mathrm{d} t\
			\textstyle\frac{1}{\lf(1+\eps^2t\ri)}
			\left(
			\frac{1}{6} \alpha^2 t ^3
			- \frac{1}{2} \Omega_0 (s-2) {\beta_{\star}} t
			+\Omega _0 \varepsilon ^2 {\beta_{\star}} t ^2-\frac{1}{2}\varepsilon ^2 {\beta_{\star}}^2
			\right)
			g_{\star}\gcrm' \\
			= \frac{4\Omega _0}{\alpha ^2}
			\int _y ^\eta
			\mathrm{d} t \: (1 + \varepsilon^2 t) \:
			U _{\beta_{\star}}(t)
			g_{\star}
			g_{\star}' + \mathcal{O}(\varepsilon^2 \eta^7) g_\beta^2(y),
		}
		where we have used the bound \eqref{eq: derivative gcrm} and the monotonicity of $ g_{\star} $ for $ y \geq y_{{\beta_{\star}}} $. The first term on the r.h.s. can be rewritten by means of \eqref{eq: gvbm variational}:
		\bmln{
			\frac{4\Omega _0}{\alpha ^2}
			\int _y ^\eta
			\mathrm{d} t \: (1 + \varepsilon^2 t) \:
			U _{\beta_{\star}}(t)
			g_{\star}
			g_{\star}'
			=\frac{2\Omega _0}{\alpha ^2}
			\int _y ^\eta
			\mathrm{d} t \: (1 + \varepsilon^2 t) \:
			\left[
			\textstyle\frac{1}{2}
			\left(
				g_{\star}'
			\right) ^2
			-\frac{1}{2\pi} g_{\star} ^4
			+\mu_{\star} g_{\star} ^2
			\right]'
			\\
			+\frac{2\Omega _0 \varepsilon ^2}{\alpha ^2}
			\int _y ^\eta
			\mathrm{d} t\
			\left(
			g_{\star} '
			\right) ^2
			-\frac{4\Omega _0 \varepsilon ^2}{\alpha ^2}
			\int _y ^\eta
			\mathrm{d} t \: (1 + \varepsilon^2 t) \:
			t ^3 v(t) g_{\star} g_{\star}'=
			\\
			=\frac{2\Omega _0}{\alpha ^2}
			\left(
			1
			+\varepsilon ^2 \eta
			\right)
			\left[
			-\textstyle\frac{1}{2\pi}g_{\star} ^4(\eta)
			+\mu_{\star} g_{\star} ^2(\eta)
			\right]
			-\frac{2\Omega _0}{\alpha ^2}\lf(1+\eps^2y\ri)
			\left[
			\textstyle\frac{1}{2}
			\left(
				g_{\star}'(y)
			\right) ^2
			-\frac{1}{2\pi}
			g_{\star} ^4(y)
			+\mu_{\star} g_{\star} ^2(y)
			\right]-
			\\
			+\frac{2\Omega _0 \varepsilon ^2}{\alpha ^2}
			\int _y ^\eta
			\mathrm{d} t\
			\left[
			\textstyle\frac{1}{2}
			\left(
				g_{\star}'
			\right) ^2
			+\frac{1}{2\pi}g_{\star} ^4
			-\mu_{\star} g_{\star} ^2
			\right]
			-\frac{4\Omega _0 \varepsilon ^2}{\alpha ^2}
			\int _y ^\eta
			\mathrm{d} t \: (1 + \varepsilon^2 t) \:
			t ^3 v(t) g_{\star}  g_{\star}'	\\
			= \frac{2\Omega _0}{\alpha ^2}
			\mu_{\star} g_{\star} ^2(\eta)
			-\frac{2\Omega _0}{\alpha ^2}
			\left[
			\textstyle\frac{1}{2}
			\left(
				g_{\star}'(y)
			\right) ^2
			-\frac{1}{2\pi}
			g_{\star} ^4(y)
			+\mu_{\star} g_{\star} ^2(y)
			\right] + \mathcal{O}(\varepsilon^2\eta^7) g_{\star}^2(y).
		}
		Putting together the above estimate with \eqref{eq: F computed 0} and \eqref{eq: F computed 1}, we obtain the result.
	\end{proof}

\noindent
Thanks to Proposition \ref{pro: alt expression F}, the cost function $ K$ can also be expressed as
\begin{equation}
	\label{eq: K alternative}
	K(y) = - \frac{\Omega _0}{\alpha ^2} \left(
				g_{\star}'
			\right) ^2
			+ \frac{2\Omega _0}{\alpha ^2}\left[ \frac{\alpha^2}{4 \Omega_0} + \frac{1}{2} \alpha^2 y^2
			-\frac{1}{2\pi} g_{\star}^4
			+\mu_{\star}
			\right] g_{\star} ^2
			+
			\begin{cases}
				R_{+}(y) + R_{+},	&	\mbox{if } y \geq y_{\beta_{\star}},	\\
				R_{-}(y) + R_{-},	&	\mbox{if } y \leq y_{\beta_{\star}}.
			\end{cases}
\end{equation}
This alternative expression will play an important role in the proof of its positivity, exactly as for $ K^{\mathrm{gv}} $. Another important ingredient of the proof is also the closeness of $ K $ to $ K^{\mathrm{gv}} $ as $ \varepsilon \to 0 $:

	\begin{lem}
		\label{lem: K close kgv}
		\mbox{}	\\
		For any $ \Omega_0 > 0 $ and $ y \in \widetilde{\mathcal{A}}_{\eta} $
		\begin{equation}
			\label{eq: K close kgv}
			K(y) - K^{\mathrm{gv}}(y) = \mathcal{O}(\varepsilon^2|\log\varepsilon|^{\infty}).
		\end{equation}
	\end{lem}

	\begin{proof}
		The result is a direct consequence of the pointwise estimate \eqref{eq: infty est diff gvm}.
	\end{proof}

\noindent
The above Lemma in combination with Proposition \ref{pro: kgv positive} might seem to give also the positivity of $ K $ inside $ \widetilde{\mathcal{A}}_{\eta} $. However this is not the case because, although we proved that $ K^{\mathrm{gv}} $ is positive on the whole real line, we did not provide any lower bound to it. In fact even by just looking at its minima, one could conclude from \eqref{eq: kgv crucial pro} that $ K^{\mathrm{gv}}(y_2) \geq K^{\mathrm{gv}}(y_1) g_{\mathrm{gv}}^2(y_2)/g_{\mathrm{gv}}^2(y_1) $, where $ y_{2}, y_1 $ are the positions of the minimum point and the preceding maximum point (consider for simplicity the half-line $ \mathbb{R}^+ $). Now even if $ K^{\mathrm{gv}}(y_1) > C > 0 $ as it occurs for instance at the origin, the ratio between the densities can become extremely small in $ {\mathcal{A}_{\eta}} $. In addition to that the inequality holds true only for the minima of $ K^{\mathrm{gv}} $ and it might be that it has no minimum inside $ \widetilde{\mathcal{A}}_{\eta} $ or $ {\mathcal{A}_{\eta}} $, in which case we only know that it is positive there, but without any meaningful lower bound.

In fact we will be able to prove positivity of $ K $ only in domain strictly smaller than $ {\mathcal{A}_{\eta}} $, because of the additional constant terms $ R_{\pm} $ in \eqref{eq: K alternative}, irrespective of their smallness. We thus set
\begin{equation}
	\label{eq: annm}
	\mathcal{A}_{>} : =
	\left\{
		y \in [-\eta,\eta]\: \big|\:
		g_{\star} ^2 (y)
		\ge \eta ^6 \max\left\{
			g_{\star} ^2 (\eta),
			g_{\star} ^2 (-\eta)
		\right\}
	\right\}.
\end{equation}
By monotonicity of $ g_{\star} $ for large $ y $ is easy to see that $ \mathcal{A}_{>} = [- y_{-}, y_{+}] $ with $ y_{\pm} \to \infty $ as $ \varepsilon \to 0 $. Notice also that $ g_{\star} $ is very small at the boundary of $ \mathcal{A}_{>} $, although not as small as $ g_{\star}(\eta) $.

We can now state the main result of this Section:

	\begin{pro}[Positivity of $ K $]
		\label{pro: K positive}
		\mbox{}	\\
		If $ \Omega_0 > \Omega_{\mathrm{c}} $ as $ \varepsilon \to 0 $,
		\begin{equation}
			K(y) > 0,	\qquad		\mbox{for any } y \in \mathcal{A}_{>}.
		\end{equation}

	\end{pro}

	\begin{proof}
		As in the proof of positivity of $ K^{\mathrm{gv}} $ in Proposition \ref{pro: kgv positive} the key idea is to show that positivity at the origin is indeed sufficient to get the result. This in turn is easily inherited from positivity of $ K^{\mathrm{gv}} $ whenever $ \Omega_0 > \Omega_{\mathrm{c}} $, via \eqref{eq: K close kgv}: by \eqref{eq: kgv0}
		\begin{displaymath}
			K^{\mathrm{gv}}(0) > C > 0,
		\end{displaymath}
		but
		\begin{displaymath}
			K(0) - K^{\mathrm{gv}}(0) = \mathcal{O}(\varepsilon^2|\log\varepsilon|^{\infty}),
		\end{displaymath}
		and thus $ K(0) > C > 0 $ for a possibly different constant $ C $.

		The rest of the proof follows the same line of reasoning of the proof of Proposition \ref{pro: kgv positive}. There are however two complications: first we have two alternative expressions of $ K $ in $ [y_{{\beta_{\star}}},\eta] $ and $ [-\eta,y_{{\beta_{\star}}}] $ respectively. Recall that $ y_{{\beta_{\star}}} = o(1) $ denotes the unique maximum point of $ g_{\star} $. Second the presence of the constant terms $ R_{\pm} $ in \eqref{eq: K alternative} is very annoying and in fact it is responsible of the restriction to $ \mathcal{A}_{>} $.

		In order to handle the first issue it is sufficient to take into account the two intervals $ [y_{{\beta_{\star}}},\eta] $ and $ [-\eta,y_{{\beta_{\star}}}] $ separately and use a different expressions for $ K $ (see \eqref{eq: K alternative}).

		The second issue on the other hand leads to the introduction of the modified cost function
		\begin{equation}
			\label{eq: tilde K}
			\widetilde{K}(y) : = K(y)  - \delta_{\varepsilon} g_{\star}^2(y) - R_{\pm}
		\end{equation}
		for some
		\begin{equation}
			\label{eq: deps}
			 0 < \delta_{\varepsilon} \ll \eta^{-2} \ll 1
		\end{equation}
		to be chosen later. Here we have used a compact notation to mean that we subtract $ R_+ $ (resp. $ R_- $) in $ [y_{{\beta_{\star}}},\eta] $ (resp. $ [-\eta,y_{{\beta_{\star}}}] $).

		Now we observe that if $ \Omega_0 > \Omega_{\mathrm{c}} $
		\begin{equation}
			\widetilde{K}(y_{{\beta_{\star}}}) = \widetilde{K}(0) + o(1) = K(0) + o(1) > 0,
		\end{equation}
		thanks to the pointwise estimate \eqref{eq: K close kgv} and since $ \left\| g_{\star} \right\|_{\infty} \leq C $. It is interesting to remark that this is the only point in the proof where we use the condition $ \Omega_0 > \Omega_{\mathrm{c}} $, although several later estimates are affected by this one.
		Moreover at the boundary of the domain we have
		\begin{equation}
			\widetilde{K}(\pm\eta) = \left( \textstyle\frac{1}{2} - \delta_{\varepsilon} + \mathcal{O}(\varepsilon^3 \eta^7) \right) g_{\star}^2(\pm\eta) > 0.
		\end{equation}
		Therefore in order to exclude that $ \widetilde{K} $ becomes negative, it suffices to prove that it is positive at any possible {\it global} minimum point $ - \eta < y_m < \eta $.

		We claim that, for $ \Omega_0 > \Omega_{\mathrm{c}} $, any global minimum point $ y_m $ of $ \widetilde{K} $ must satisfy the condition
		\begin{equation}
			\label{eq: ym condition}
			|y_m| \gg 1.
		\end{equation}
		The reason is the pointwise estimate \eqref{eq: K close kgv} and the observation contained in Proposition \ref{pro: kgv positive}: $ \widetilde{K} $ and $ K^{\mathrm{gv}} $ are pointwise close and therefore we would have
		\begin{displaymath}
			K^{\mathrm{gv}}(y_m) \leq \min [K(-\eta),K(\eta)] + o(1) \leq \tx\frac{1}{2} \min [g_{\star}^2(-\eta),g_{\star}^2(\eta)]+ o(1) = o(1),
		\end{displaymath}
		 which in turn implies $ K^{\mathrm{gv}}(y_m) \leq 0 $ since $ K^{\mathrm{gv}} $ is independent of $ \varepsilon $. For $ \Omega_0 > \Omega_{\mathrm{c}} $ this contradicts the statement of Proposition \ref{pro: kgv positive}.

		The key point in the proof is the following property: let $ y_m $ be a point where $ \widetilde{K} $ reaches its global minimum $ \widetilde{K}(y_m) < \widetilde{K}(y_{{\beta_{\star}}}) $ (otherwise there would be nothing to prove) and $ y_M $ any maximum point of $ \widetilde{K} $ such that $ y_M < y_m $, if $ y_m > y_{{\beta_{\star}}} $, or $ y_M > y_m $ in the opposite case $ y_m < y_{{\beta_{\star}}} $. Notice that such a maximum needs not to be the global maximum but its existence is a consequence of smoothness of $ \widetilde{K} $ and the inequalities
		\begin{displaymath}
			\widetilde{K}(\pm\eta) < \widetilde{K}(y_{{\beta_{\star}}}), \qquad	\widetilde{K}(y_m) < \widetilde{K}(y_{{\beta_{\star}}}).
		\end{displaymath}
		Then we are going to prove that
		\begin{equation}
			\label{eq: K crucial pro 1}
			\frac{\widetilde{K}(y_m)}{g_{\star}^2(y_m)} \geq \frac{\widetilde{K}(y_M)}{g_{\star}^2(y_M)} + o(1).
		\end{equation}
		Now suppose that this is true, then we can pick a maximum point $ y_M $ of $ \widetilde{K} $ where $ \widetilde{K}(y_M) \geq \widetilde{K}(y_{{\beta_{\star}}}) > C > 0 $. In addition it must be
		\begin{equation}
			y_M = \mathcal{O}(1),
		\end{equation}
		because $ \widetilde{K}(y) \leq C g^2(y) $ and the decay estimate \eqref{eq: better estimate on gvbm} or the pointwise estimate \eqref{eq: infty est diff gvm} imply that $ \widetilde{K}(y) = o(1) $, if $ |y| \gg 1 $. Hence by the lower bound \eqref{eq: better estimate on gvbm} $ g_{\star}(y_M) \geq C > 0 $, \eqref{eq: K crucial pro 1} yields
		\begin{equation}
			\label{eq: K crucial pro 2}
			\widetilde{K}(y_m) \geq g_{\star}^2(y_m) \left( C^{-2} \widetilde{K}(y_{{\beta_{\star}}}) + o(1) \right) \geq C_0 g_{\star}^2(y_m) > 0
		\end{equation}
		for some $ C_0 > 0 $. In fact we have obtained something more: for any $ y \in {\mathcal{A}_{\eta}} $ either $ \widetilde{K}(y) \geq \min \{ \widetilde{K}(-\eta), \widetilde{K}(\eta) \} > 0  $ or
		\begin{displaymath}
			\frac{\widetilde{K}(y)}{g_{\star}^2(y)} \geq \frac{\widetilde{K}(y_m)}{g_{\star}^2(y)} \geq C \frac{g_{\star}^2(y_m)}{g_{\star}^2(y)} > 0.
		\end{displaymath}
		Either way $ \widetilde{K} $ is positive everywhere in $ [-\eta,\eta] $. Moreover the positivity of $ \widetilde{K} $ implies that
		\begin{displaymath}
			K(y) > R_{\pm} + \delta_{\varepsilon} g_{\star}^2(y) \geq 0,	\qquad \mbox{if } g_{\star}^2(y) \geq \delta_{\varepsilon}^{-1} |R_+|,
		\end{displaymath}
		for any $ y \in {\mathcal{A}_{\eta}} $. If now we restrict the inequality to $ \mathcal{A}_{>} $ and we choose, e.g., $ \delta_{\varepsilon} = \eta^{-3} $, the estimates \eqref{eq: tails} imply that inside $ \mathcal{A}_{>} $
		\begin{displaymath}
			g_{\star}^2(y) \geq \eta^6 \max\left\{g_{\star}^2(-\eta), g_{\star}^2(\eta) \right\} \gg C\eta^{5} \max\left\{g_{\star}^2(-\eta), g_{\star}^2(\eta) \right\} \geq \delta_{\varepsilon}^{-1} |R_{\pm}|,
		\end{displaymath}
		so that $ K(y) $ is strictly positive for any $ y \in \mathcal{A}_{>} $. In fact a closer look to the chain of inequalities reveals that we have proven something more, i.e., for $ \varepsilon $ small enough
		\begin{equation}
			\label{eq: lb K}
			K(y) \geq \eta^{-3} g_{\star}^2(y),	\qquad	\mbox{for any } y \in \mathcal{A}_{>}.
		\end{equation}

		We have now to prove \eqref{eq: K crucial pro 1}. Recall the assumption: we have a global minimum of $ \widetilde{K} $ at $ y_m $ and a maximum at $ y_M $, which is on the left (resp. right) of $ y_m $, if $ y_m > y_{{\beta_{\star}}} $ (resp. $ y_m < y_{{\beta_{\star}}} $). The idea is the same used in the proof of Proposition \ref{pro: kgv positive}: the derivative of $ \widetilde{K} $, i.e.,
		\begin{equation}
			\label{eq: kt derivative}
			\widetilde{K}^{\prime}(y) =\left[
			(1 - \delta_{\varepsilon}) g_{\star}'
			+ \textstyle\frac{1}{1+\eps^2y}
			\left(
				2\Omega _0 y
				+\Omega _0 \varepsilon ^2 y ^2
				- {\beta_{\star}} \varepsilon ^2
			\right)g_{\star}
		\right] g_{\star},
		\end{equation}
		must vanish both at $ y_m $ and $ y_M $ and therefore
		\begin{equation}
			\label{eq: extreme K}
			g_{\star}^{\prime}(y_{m,M}) = - \frac{1}{(1 - \delta_{\varepsilon})(1 + \varepsilon^2 y_{m,M})}\left(
				2\Omega _0 y_{m,M}
				+\Omega _0 \varepsilon ^2 y_{m,M} ^2
				- {\beta_{\star}} \varepsilon ^2
			\right)g_{\star}(y_{m,M}).
		\end{equation}
		The second derivative of $K $ can be computed as well:
		\begin{equation}
			\widetilde{K}''(y) =
			(1 - \delta_{\varepsilon}) \left(g_{\star}'\right)^2 + 4 \Omega_0 y g_{\star} g_{\star}' + \left[ (1 - \delta_{\varepsilon})\left( \alpha^2 y^2 + \textstyle\frac{2}{\pi} g_{\star}^2 - 2 \mu_{\star} \right) + 2 \Omega_0 + \mathcal{O}(\varepsilon^2 \eta^5) \right] g_{\star}^2,
		\end{equation}
		so that at any extreme point of $ \widetilde{K} $, one has
		\begin{equation}
			\label{eq: kt primeprime ext}
			\widetilde{K}''(y_{m,M}) = \left[ \Omega_0^2 (s-2) y_{m,M}^2 + 2 \Omega_0 + \textstyle\frac{2}{\pi} g_{\star}^2(y_{m,M}) - 2 \mu_{\star} + o(1) \right] g_{\star}^2(y_{m,M}),
		\end{equation}
		where we have exploited the condition \eqref{eq: deps}. Similarly by \eqref{eq: K alternative} we get
		\begin{equation}
			\label{eq: kt ext}
			\widetilde{K}(y_{m,M}) = \left[ \Omega_0 \textstyle\frac{s-2}{s+2} y_{m,M}^2  + \textstyle\frac{1}{2} - \textstyle\frac{1}{\pi \Omega_0(s+2)} g_{\star}^2(y_{m,M}) - \frac{2}{\Omega_0(s+2)} \mu_{\star} + o(1) \right] g_{\star}^2(y_{m,M}),
		\end{equation}
		and a direct comparison between \eqref{eq: kt primeprime ext} and \eqref{eq: kt ext} yields
		\begin{equation}
			\label{eq: key identity}
			\frac{\widetilde{K}(y_{m,M})}{g_{\star}^2(y_{m,M})} = \frac{s-2}{s+2} -\frac{g_{\star}^2(y_{m,M})}{\pi \Omega_0 (s+2)} + \frac{\widetilde{K}''(y_{m,M}) }{\Omega_0 (s+2) g_{\star}^2(y_{m,M})} + o(1).
		\end{equation}
		Now this is the key identity because by assumption (recall also \eqref{eq: ym condition})
		\begin{displaymath}
			\widetilde{K}''(y_{M}) \leq 0 \leq \widetilde{K}''(y_{m}),	\qquad		g_{\star}(y_m) < g_{\star}(y_M),
		\end{displaymath}
		so that
		\begin{displaymath}
				\frac{\widetilde{K}(y_{m})}{g_{\star}^2(y_{m})} \geq \frac{\widetilde{K}(y_{M})}{g_{\star}^2(y_{M})} + o(1),
		\end{displaymath}
		i.e., \eqref{eq: K crucial pro 1} is proven. Note that the fact that we have two different explicit expressions of $ \widetilde{K} $ for $ y > y_{{\beta_{\star}}} $ and $ y < y_{{\beta_{\star}}} $ did not affect the proof, because the difference between the two expressions is $ o(1) $ and therefore can be included in the error term.
	\end{proof}

\section{Energy Asymptotics}
\label{sec: energy}

We attack in this Section the proof of Theorem \ref{teo: energy asymptotics}, which will imply the main result of the paper. The result is obtained by combining upper (Proposition \ref{pro: upper bound}) and lower (Proposition \ref{pro: lower bound}) bounds on $ E^{\mathrm{GP}} $.

\subsection{Upper Bound}
\label{sec: energy ub}

The upper bound on $ E^{\mathrm{GP}} $ is stated in next

	\begin{pro}[GP energy upper bound]
		\label{pro: upper bound}
		\mbox{}	\\
		As $ \varepsilon \to 0 $,
		\begin{equation}
			\label{eq: upper bound}
			E^{\mathrm{GP}} \leq \frac{E_{\star}^\mathrm{gv}}{\varepsilon^4} + \mathcal{O}(1).
		\end{equation}
	\end{pro}

	\begin{proof}
		The proof is rather simple because it is sufficient to test $ \mathcal{E}^{\mathrm{GP}} $ on suitable trial function of the form
		\begin{equation}
			\Psi_{\mathrm{trial}}(\mathbf{x}) :
			=\frac{1}{\sqrt{2\pi}\varepsilon}
			g_{\mathrm{trial}}\left(
			\textstyle\frac{|\mathbf{x}|-1}{\varepsilon ^2}
			\right)
			\exp \left\{i\left\lfloor\Omega +{\beta_{\star}}\right\rfloor\theta \right\},
		\end{equation}
		where $ g_{\mathrm{trial}} $ coincides up to a normalization constant with $ g_{\star} $ within $ {\mathcal{A}_{\eta}} $ and is suitably regularized outside. The calculation is rather straightforward and we omit it for the sake of brevity. Note that the remainder $ \mathcal{O}(1) $ is entirely due to the fact that the phase $ \Omega + {\beta_{\star}} $ might not be an integer number. Otherwise one would obtain a much better error term $ \mathcal{O}(\varepsilon^4) $.
	\end{proof}

\subsection{Lower Bound}
\label{sec: energy lb}

A lower bound for $ E^{\mathrm{GP}} $ matching the upper bound of Proposition \ref{pro: upper bound} is formulated in next

	\begin{pro}[GP energy lower bound]
		\label{pro: lower bound}
		\mbox{}	\\
		If $ \Omega_0 > \Omega_{\mathrm{c}} $, as $ \varepsilon \to 0 $,
		\begin{equation}
			\label{eq: lower bound}
			E^{\mathrm{GP}} \geq \frac{E_{\star}^\mathrm{gv}}{\varepsilon^4} + \mathcal{O}(\varepsilon^{\infty}).
		\end{equation}
	\end{pro}

	\begin{proof}
		We first restrict the integration in the GP energy functional to the domain $ {\mathcal{A}_{\eta}} $ (recall its definition in \eqref{eq: ring}: to this purpose we just have to observe that all the three terms in the GP energy functional are pointwise positive and thus we can simply drop their integrals outside $ {\mathcal{A}_{\eta}} $. Of course $ \psi^{\mathrm{GP}} $ is not normalized in $ L^2({\mathcal{A}_{\eta}}) $ but the exponential decay proven in Proposition \ref{pro: gpm exp small} guarantees that
		\begin{equation}
			\left\| \psi^{\mathrm{GP}} \right\|_{L^2({\mathcal{A}_{\eta}})} = 1 + \mathcal{O}(\varepsilon^{\infty}),
		\end{equation}
		by taking $ \eta_0 $ large enough.

		The first step in the proof is a splitting of the energy, in order to extract the leading order term $ E_{\star}^\mathrm{gv}/\varepsilon^4 $. This is now rather standard and we do not spell all the details of the computation. We just note that one sets
		\begin{equation}
			\label{def: u}
			\psi^{\mathrm{GP}}(\mathbf{x})
			=: \frac{1}{\sqrt{2\pi}\varepsilon}
			u\left(
			x,\vartheta
			\right)
				g_{\star} \left(
			\textstyle\frac{x-1}{\varepsilon ^2}
			\right)
		e^{i\left(\Omega +{\beta_{\star}}\right)\theta}.
		\end{equation}
		Since $ \Omega + {\beta_{\star}} $ needs not to be an integer, $ u $ is not single-valued in general, but
		\begin{equation}
			\label{eq: u semiperiodic}
			u(x,\vartheta+2 k \pi) = e^{-i 2\pi k \left(\Omega +{\beta_{\star}}\right)} u(x,\vartheta),
		\end{equation}
		for any $ k \in \mathbb{Z} $. A part from that $ u $ is finite for any $ \mathbf{x} \in {\mathcal{A}_{\eta}} $ thanks to the strict positivity of $ g_{\star} $. A long but simple computation using the variational equation for $ g_{\star} $ gives
		\begin{equation}
			\label{eq: splitting}
			E^{\mathrm{GP}} \geq \frac{E_{\star}^\mathrm{gv}}{\varepsilon^4} + \frac{\mathcal{E}[u]}{2 \pi \varepsilon^2} + \mathcal{O}(\varepsilon^{\infty}),
		\end{equation}
		where the inequality is mainly due to the restriction of the integration domain and setting $ y = 1 + \varepsilon^2 x $ for short
		\begin{equation}
			\label{eq: eu}
			\mathcal{E}[u] = \int _{\mathcal{A}_{\eta}}
			\mathrm{d} \mathbf{x}\
			g_{\star} ^2(y)
			\left\{
			\textstyle\frac{1}{2}
			\left|
				\nabla u
			\right| ^2
			+{\bf a}
			\cdot {\bf j} _u
			+\frac{1}{2\pi \varepsilon ^4} g_{\star} ^2(y)
			(
			1
			-|
				u
			| ^2
			)^2
			\right\},
		\end{equation}
		\begin{equation}
			{\bf a}(x):=\left(
			\frac{\Omega + {\beta_{\star}}}{x}
			-\Omega x
		\right){\bf e_\vartheta},
		\end{equation}
		and the superfluid current is defined in \eqref{eq: current}. The rest of the proof is devoted to prove that
		\begin{equation}
			\mathcal{E}[u] \geq \mathcal{O}(\varepsilon^{\infty}).
		\end{equation}

		In order to exploit the cost function trick mentioned in Section \ref{sec: heuristics} and the positivity of $ K $ proven in Proposition \ref{pro: K positive}, we need to restrict again the integration domain in $ \mathcal{E}[u] $ to $ \mathcal{A}_{>}^{\mathrm{2D}} \subset {\mathcal{A}_{\eta}} $, where
		\begin{displaymath}
			\mathcal{A}_{>}^{\mathrm{2D}} = \left\{ \mathbf{x} \in \mathbb{R}^2 \: \big| \: \left|1 - |\mathbf{x}| \right|/\varepsilon^2 \in \mathcal{A}_{>} \right\}.
		\end{displaymath}
		The only annoying term is the only one which is not positive, i.e., the second one in \eqref{eq: eu}:
		\bml{
			\label{eq: ang momentum layer}
			\left| \int_{{\mathcal{A}_{\eta}} \setminus \mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \: g_{\star}^2(y) \: {\bf a} \cdot {\bf j} _u \right| \leq \left\| {\bf a} \right\|_{L^{\infty}({\mathcal{A}_{\eta}})} \int_{{\mathcal{A}_{\eta}} \setminus \mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \: g_{\star}^2(y) \left|u\right| \left| \nabla_{\vartheta} u \right| \\
			\leq C \eta^2 \int_{{\mathcal{A}_{\eta}} \setminus \mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \: \left| \psi^{\mathrm{GP}} \right| \left| \nabla_{\vartheta} \psi^{\mathrm{GP}} \right|
			\leq C \varepsilon^2 \eta^4 \left\| \psi^{\mathrm{GP}} \right\|_{L^{\infty}({\mathcal{A}_{\eta}}\setminus \mathcal{A}_{>}^{\mathrm{2D}})} \left\| \nabla \psi^{\mathrm{GP}} \right\|_{L^{\infty}({\mathcal{A}_{\eta}})}	\\
			 \leq C \varepsilon^{-4} \eta^4 \left\| \psi^{\mathrm{GP}} \right\|_{L^{\infty}({\mathcal{A}_{\eta}}\setminus \mathcal{A}_{>}^{\mathrm{2D}})},
			}
			where we have used the bound $ \left\| \nabla \psi^{\mathrm{GP}} \right\|_{L^{\infty}({\mathcal{A}_{\eta}})} \leq C \varepsilon^{-6} $, following from
			\begin{equation}
				\label{eq: GN}
				\left\| \nabla \psi \right\|_{\infty} \leq C \left( \left\| \Delta \psi \right\|_{\infty}^{1/2} \left\| \psi \right\|_{\infty}^{1/2} + \left\| \psi \right\|_{\infty} \right),
			\end{equation}
			which can be proven from Gagliardo-Nirenberg inequalities exactly as in \cite[Lemma 5.1]{CRY}. However the lower bound \eqref{eq: better estimate on gvbm} easily implies that if we set $ \mathcal{A}_{>} =: [-y_-,y_+] $, then
			\begin{displaymath}
				y_{\pm} = \eta(1 + o(1)),
			\end{displaymath}
			so that
			\begin{equation}
				\left| \psi^{\mathrm{GP}} \right|_{\partial \mathcal{A}_{>}^{\mathrm{2D}}} \leq \mathcal{O}(\varepsilon^{\infty}),
			\end{equation}
			again by \eqref{eq: gpm exp small} and the arbitrariness in the choice of $ \eta_0 $. Hence \eqref{eq: ang momentum layer} yields an error which can be made smaller than any power of $ \varepsilon $ and we get the lower bound
			\begin{equation}
				\mathcal{E}[u] \geq  \int _{\mathcal{A}_{>}^{\mathrm{2D}}}
				\mathrm{d} \mathbf{x}\
				g_{\star} ^2(y)
				\left\{
				\textstyle\frac{1}{2}
				\left|
				\nabla u
				\right| ^2
				+{\bf a}
				\cdot {\bf j} _u
				+\frac{1}{2\pi \varepsilon ^4} g_{\star} ^2(y)
				(
				1
				-|
				u
				| ^2
				)^2
				\right\} + \mathcal{O}(\varepsilon^{\infty}).
			\end{equation}

			We can now finally integrate by the angular momentum term by using the potential function $ F $ defined in \eqref{eq: potential f}: it is trivial to verify that
			\begin{equation}
				2 g_{\star}^2\left(\textstyle\frac{x - 1}{\varepsilon^2}\right) {\bf a}(x) = - \partial_x F \left( \textstyle\frac{x - 1}{\varepsilon^2} \right) {\bf e}_{\vartheta},
			\end{equation}
			so that
			\bml{
				\label{eq: int by parts}
				 \int_{\mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \: g_{\star}^2(y) \: {\bf a} \cdot {\bf j} _u = - \frac{1}{2} \int_0^{2\pi} \mathrm{d} \vartheta \int_{1 - \varepsilon^2 y_-}^{1 + \varepsilon^2 y_+} \mathrm{d} x \: \partial_x F \left( \textstyle\frac{x - 1}{\varepsilon^2} \right)  \Re \left[ i u(x,\vartheta) \partial_{\vartheta} u^*(x,\vartheta) \right] 	\\
				 = \frac{1}{2} \int_0^{2\pi} \mathrm{d} \vartheta \int_{1 - \varepsilon^2 y_-}^{1 + \varepsilon^2 y_+} \mathrm{d} x \: F\left( \textstyle\frac{x - 1}{\varepsilon^2} \right) \Re \left[ i \partial_x u(x,\vartheta) \partial_{\vartheta} u^*(x,\vartheta) + i u(x,\vartheta) \partial^2_{x,\vartheta} u^*(x,\vartheta) \right]	\\
				  - \frac{1}{2} \int_0^{2\pi} \mathrm{d} \vartheta \: \Big| F\left( \textstyle\frac{x - 1}{\varepsilon^2} \right) \Re \left[ i u(x,\vartheta) \partial_{\vartheta} u^*(x,\vartheta) \right] \Big|_{1-\varepsilon^2 y_-}^{1 + \varepsilon^2 y_+}.
			}
			The boundary term can be easily proven to provide an exponentially small correction: consider, e.g., the term at $ 1 + \varepsilon^2 y_+ $, since $ \left| F(y_{\pm}) \right| \leq C \eta^8 g_{\star}^2(\pm\eta) $, one can reconstruct a term, which can be bounded exactly as \eqref{eq: ang momentum layer}. The result is an error $ \mathcal{O}(\varepsilon^{\infty}) $. The rest is integrated by parts once more but this time w.r.t. $ \vartheta $:
			\bmln{
				\frac{1}{2} \int_0^{2\pi} \mathrm{d} \vartheta \int_{1 - \varepsilon^2 y_-}^{1 + \varepsilon^2 y_+} \mathrm{d} x \: F\left( \textstyle\frac{x - 1}{\varepsilon^2} \right) \Re \left[ i u(x,\vartheta) \partial^2_{x,\vartheta} u^*(x,\vartheta) \right]		\\
				 = - \frac{1}{2} \int_0^{2\pi} \mathrm{d} \vartheta \int_{1 - \varepsilon^2 y_-}^{1 + \varepsilon^2 y_+} \mathrm{d} x \: F\left( \textstyle\frac{x - 1}{\varepsilon^2} \right) \Re \left[ i \partial_{\vartheta} u(x,\vartheta) \partial{x} u^*(x,\vartheta) \right],
			}
			where the vanishing of boundary terms is due to the periodicity of $ u^* \partial_x u $ and  its complex conjugate (compare with \eqref{eq: u semiperiodic}). Altogether we have thus obtained that
			\bml{
				 \int_{\mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \: g_{\star}^2(y) \: {\bf a} \cdot {\bf j} _u = \int_{\mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \: F\left( \textstyle\frac{x - 1}{\varepsilon^2} \right) \Re \left[ i \nabla_x u(x,\vartheta) \nabla_{\vartheta} u^*(x,\vartheta) \right]	 + \mathcal{O}(\varepsilon^{\infty})\\
				  \geq - \int_{\mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \: \left| F\left( \textstyle\frac{x - 1}{\varepsilon^2} \right) \right| \left| \nabla u \right|^2 + \mathcal{O}(\varepsilon^{\infty}) = - \int_{\mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \: F\left( \textstyle\frac{x - 1}{\varepsilon^2} \right) \left| \nabla u \right|^2  + \mathcal{O}(\varepsilon^{\infty}),
			}
			and therefore
			\bml{
				\label{eq: lb final}
				\mathcal{E}[u] \geq \int_{\mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \: \left\{ K\left( \textstyle\frac{x - 1}{\varepsilon^2} \right) \left| \nabla u \right|^2 + \frac{1}{2\pi \varepsilon ^4} g_{\star} ^2(y)
				(
				1
				-|
				u
				| ^2
				)^2
				\right\} + \mathcal{O}(\varepsilon^{\infty})	\\
				\geq \eta^{-3} \int_{\mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \: g_{\star}^2(y) \left| \nabla u \right|^2 + \mathcal{O}(\varepsilon^{\infty}) \geq \mathcal{O}(\varepsilon^{\infty}),
			}
			thanks to Proposition \ref{pro: K positive} and in particular \eqref{eq: lb K}.
	\end{proof}

\section{Giant Vortex Transition}
\label{sec: gv transition}

In this Section we prove the results regarding absence of vortices and total vorticity of the condensate.

	\begin{proof}[Proof of Theorem \ref{teo: no vortices}]
		Combining \eqref{eq: lb final} with \eqref{eq: splitting} and the upper bound proven in Proposition \ref{pro: upper bound}, we get
		\begin{equation}
			\int_{\mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \: g_{\star} ^4(y)
				(
				1
				-|
				u
				| ^2
				)^2 = \mathcal{O}(\varepsilon^6),
		\end{equation}
		which already means that $ |u| $ can not differ too much from $ 1 $. To deduce the pointwise estimate of Theorem \ref{teo: no vortices}, we need to combine this with an estimate of $ \left\| \nabla  u \right\|_{\infty} $.

		As in \cite[Lemma 4.3]{CRY} we obtain from \eqref{eq: GP variational} and \eqref{eq: gvbm variational} the following variational equation for $ u $:
		\begin{equation}
			\label{eq: u variational}
			-\tx\frac{1}{2} g_\beta \Delta u
			-\frac{1}{\varepsilon ^2} g_{\star}' \partial _x u
			-i g_{\star} {\bf a} \cdot \nabla u
			+\frac{1}{\pi \varepsilon ^4} g_{\star} ^3
			\left(
			\left|
				u
			\right| ^2
			-1
			\right) u
			=\left(
			\mu^{\mathrm{GP}}
			-\frac{1}{\varepsilon ^4} \mu_{\star}
			\right) g_{\star} u,
		\end{equation}
		which yields (recall the definition of $ \mathcal{A}_{\mathrm{bulk}} \subset \mathcal{A}_{>}^{\mathrm{2D}} $ in \eqref{eq: bulk})
		\begin{displaymath}
			\left\| \Delta u \right\|_{L^{\infty}(\mathcal{A}_{\mathrm{bulk}})} \leq  C\left[
			\varepsilon ^{-2} \eta
			\left\|
				\nabla u
			\right\| _{L ^\infty(\mathcal{A}_{\mathrm{bulk}})}
			+\varepsilon ^{-4} \eta^{3a}
			\right].
		\end{displaymath}
		Now using the elliptic estimate \eqref{eq: GN} we conclude that
		\begin{equation}
			\label{eq: estimate on grad u}
			\left\|
				\nabla u
			\right\| _{L ^\infty(\mathcal{A}_{\mathrm{bulk}})}
			= \mathcal{O}\left(\varepsilon ^{-2} \eta ^{1+\frac{3a}{2}}\right).
		\end{equation}

		Suppose now that it exists $ \mathbf{x} _0 \in \mathcal{A}_{\mathrm{bulk}} $ such that $	\left|
			u(\mathbf{x} _0)
			-1
		\right|
		\geq \varepsilon^{1/2} |\log\varepsilon|^b
		$,
		for some $ b > 0  $ to be chosen later. Then from \eqref{eq: estimate on grad u} we get that
		\[
			\textstyle
			\left|
			|u|
			-1
			\right|
			\ge \tx\frac{1}{2} \varepsilon ^{1/2} |\log\varepsilon|^b,
			\qquad \mbox{for }
			\mathbf{x} \in \mathcal{B}_{\varrho}\left(
			\mathbf{x} _0\right)	\cap \mathcal{A}_a,
		\]
		with $ \varrho = \varepsilon^{5/2} |\log\varepsilon|^{b - 1 - \frac{3a}{2}} $, and
		\begin{displaymath}
			\mathcal{O}(\varepsilon^6) =  \int _{\mathcal{A}_{\mathrm{bulk}} \cap  \mathcal{B}_{\varrho}\left(
			\mathbf{x} _0\right)}
			\mathrm{d} \mathbf{x}\
			g_{\star} ^4(y) \left(
			1 - |u| ^2
			\right) ^2
			\ge C \varepsilon^{6} |\log\varepsilon|^{2b - 7a - 2},
		\end{displaymath}
		which is a contradiction for all $ b \geq 4a - 1 $.
	\end{proof}

\noindent
We now focus on the proof of Theorem \ref{teo: winding number} and for later purposes we state a useful Lemma, which is the analogue of \cite[Lemma 3.5]{CPRY3}:

	\begin{lem}
		\label{lem: deg u estimate}
		\mbox{}	\\
		Let $ \Omega_0 > \Omega_{\mathrm{c}} $ and $ R $ be a radius satisfying $ R = 1 +\bigO{\varepsilon ^2} $, then
		\begin{equation}
			\label{eq: winding u}
			\left|
			\deg\left(
				u,
				\partial \mathcal{B} _R
			\right)
			\right|	= \mathcal{O}(1)
		\end{equation}
	\end{lem}

	\begin{proof}
		We use a smooth radial cut-off function $ \chi $ with support in $ [\widetilde{R},R] $ such that $ \chi(\widetilde{R}) =0 $ and $ \chi(R) = 1 $, for some radius
	\begin{displaymath}
		\widetilde{R}
		=R
		- c \varepsilon ^2
	\end{displaymath}
	with $ c>0 $. We also require that $ \left|
			\chi
		\right|
		\le 1 $ and $ \left|
			\nabla \chi
		\right| = \mathcal{O}(\varepsilon ^{-2}) $. Then by Stokes formula
		\bml{
			\deg\left(
			u,
			\partial \mathcal{B} _R
			\right)
			=\frac{1}{\pi}
			\int _{\partial \mathcal{B} _R}
			\mathrm{d} \sigma \: \Im
			 \left(
			\frac{\nabla_{\vartheta} u}{u}
			\right) = \frac{1}{\pi}
			\int _{\partial \mathcal{B} _R}
			\mathrm{d} \sigma \: \chi(R) \: \Im
			 \left(
			\frac{\nabla_{\vartheta} u}{u}
			\right) \\
			=
			\frac{1}{\pi} \int _{\mathcal{B}_R\setminus \mathcal{B}_{\widetilde{R}} }
			\mathrm{d} \mathbf{x} \:	\nabla ^\perp \chi \cdot
			\Im \left(
			\frac{\nabla u}{u}
			\right).
		}
		Therefore
		\begin{equation}
			\label{eq: deg est}
			\left|
			\deg\left(
				u,
				\partial \mathcal{B} _R
			\right)
			\right|
			\le \frac{C}{\varepsilon ^2}
			\int _{\mathcal{B}_R\setminus \mathcal{B}_{\widetilde{R}} }
			\mathrm{d} \mathbf{x}\
			\frac{\left|\nabla u\right|}{|u|}
			\leq \frac{C}{\varepsilon^2} \left| \mathcal{B}_R\setminus \mathcal{B}_{\widetilde{R}} \right|^{1/2} \left\| \nabla u \right\|_{L^2(\mathcal{B}_R\setminus \mathcal{B}_{\widetilde{R}})},
		\end{equation}
		where we used that $ \left\| 1 - |u| \right\|_{L^{\infty}(\mathcal{B}_R\setminus \mathcal{B}_{\widetilde{R}})} = o(1) $. Now the result proven in Proposition \ref{pro: kgv positive} in fact says that for any $ \mathbf{x} \in \mathcal{B}_R\setminus \mathcal{B}_{\widetilde{R}} $ and for $ \Omega_0 > \Omega_{\mathrm{c}} $, there exists a constant $ C > 0  $ such that
		\begin{displaymath}
			K^{\mathrm{gv}}\left(\textstyle\frac{x-1}{\varepsilon^2}\right) \geq C,
		\end{displaymath}
		and thanks to \eqref{eq: K close kgv} the same inequality holds true for $ K $, i.e.,
		\begin{equation}
			K\left(\textstyle\frac{x-1}{\varepsilon^2}\right) \geq C > 0, 	\qquad		\mbox{for any } \mathbf{x} \in \mathcal{B}_R\setminus \mathcal{B}_{\widetilde{R}}.
		\end{equation}
		Going back to \eqref{eq: lb final} this yields
		\begin{displaymath}
			\mathcal{O}(\varepsilon^2) \geq \int_{\mathcal{A}_{>}^{\mathrm{2D}}} \mathrm{d} \mathbf{x} \:  K\left( \textstyle\frac{x - 1}{\varepsilon^2} \right) \left| \nabla u \right|^2 \geq  \int_{\mathcal{B}_R\setminus \mathcal{B}_{\widetilde{R}}} \mathrm{d} \mathbf{x} \:  K\left( \textstyle\frac{x - 1}{\varepsilon^2} \right) \left| \nabla u \right|^2 \geq C \left\| \nabla u \right\|_{L^2(\mathcal{B}_R\setminus \mathcal{B}_{\widetilde{R}})},
		\end{displaymath}
		which gives the result once plugged into \eqref{eq: deg est}.
	\end{proof}

	\noindent
	We are now in position to complete the estimate of the winding number of $ \psi^{\mathrm{GP}} $:

	\begin{proof}[Proof of Theorem \ref{teo: winding number}]
		We follow \cite[proof of Theorem 1.5]{CPRY3}. The positivity of $ |\psi^{\mathrm{GP}}| $ on $ \partial \mathcal{B} _R $ is guaranteed for any radius $ R = 1 +\bigO{\varepsilon ^2} $ thanks to \eqref{eq: no vortices}. A simple computation shows that
		\[
			\deg \left(
			\psi^{\mathrm{GP}},
			\partial \mathcal{B}_R
			\right)
			=\Omega
			+{\beta_{\star}}
			+\deg \left(
			u,
			\partial \mathcal{B}_R
			\right),
		\]
		which yields the result in combination with \eqref{eq: winding u}.
	\end{proof}