\subsection{Renormalization}

One result that is rather generic is that for a highly irregular chemical potential, we find a very smooth dependence
of the condensate on the coordinate $x$. 
A typical form of $\mu(x)$ and its corresponding ${\cal O}(x)$ are represented in Fig. (\ref{Smoothing}).
\begin{figure}[htp]
\begin{center}
\includegraphics[width=1.68in]{muvxfr.pdf}
\includegraphics[width=1.68in]{Condensatevxfr.pdf}
\caption{\label{Smoothing} Initial chemical potential profile $\mu(x)=6.0+0.1 \sum\limits_{n=1}^{100} \frac{1}{2\pi n}\cos (2\pi n x +\delta_n)$ (left panel) and the corresponding condensate profile (right panel).}
\end{center}
\end{figure}
On the other hand, the opposite happens for the charge density. A noisy chemical potential will translate into an
even noisier charge density. A typical form of $\mu(x)$ and its corresponding $\rho(x)$ are represented in
Fig. (\ref{Roughening}).
\begin{figure}[htp]
\begin{center}
\includegraphics[width=1.68in]{muvxfr.pdf}
\includegraphics[width=1.68in]{rhovxfr.pdf}
\caption{\label{Roughening} Initial chemical potential profile $\mu(x)=6.0+0.1 \sum\limits_{n=1}^{100} \frac{1}{2\pi n}\cos (2\pi n x +\delta_n)$ (left panel) and the corresponding charge density profile (right panel).}
\end{center}
\end{figure}

This kind of smoothing/roughening points to a renormalization of sorts, where higher harmonics in ${\cal O}$ are 
suppressed with respect to their spectral weight in $\mu$. Here we pursue this idea further. To characterize this 
renormalization quantitatively we consider a boundary chemical potential of the form presented in equation 
(\ref{noisefunc}) but considering now different values for $\alpha$.
The choice of this parameter $\alpha$ determines the degree of differentiability (smoothness) of the initial 
profile. To make the concept of renormalization more precise we consider the power spectrum of the signal $\mu(x)$ 
which is essentially proportional to $k^{-2\alpha}$; we also consider the power spectrum of the condensate 
${\cal O}(x)$ which we find numerically well approximated by $k^{-2\Delta}$. We find that 
$\Delta\simeq 1.9+1.0\alpha$ is therefore larger than $\alpha$ signaling that the weight of the high-$k$ harmonics 
is smaller in ${\cal O}$ than in $\mu$. Similarly, we approximate the power spectrum of the charge density 
$\rho(x)$ as $k^{-2\Gamma(\alpha)}$, but this time we find $\Gamma<\alpha$. Below we present a plot of $\Delta$ 
and $\Gamma$ versus $\alpha$ for a wide range of values. The error bars were computed based on considering many 
realizations of the given noise.

\begin{figure}[htp]
\begin{center}
\includegraphics[width=1.68in]{Deltavalphafr.pdf}
\includegraphics[width=1.68in]{Gammavalphafr.pdf}
\caption{\label{RGToy} Renormalization of the disorder: Condensate $\Delta=1.9+1.0\,\alpha$ (left panel) and charge
density $\Gamma= -0.88+ 1.0\,\alpha $ (right panel). This plot was made considering $k_0=2 \pi$, $\mu_0=6$ and
$\epsilon=1$. We can see that the flater the spectrum, the larger the error. This might be associated with the fact
that large momenta are more sensitive to the discretization. }
\end{center}
\end{figure}

We have gathered evidence that this behavior is also rather independent of the value of the mass of the scalar field, or, in the field theory language, it is independent of the conformal dimension of the order parameter. For example, we have confirmed a similar behavior for $m^2=0$. This universality of RG is one of the main observations of our work and its origin seems to be in the strongly coupled nature of the problem. The weak field theory intuition would dictate that $\Delta$ should be well approximated by the conformal dimension associated with the order parameter and here we verify that it is not.

It is also interesting to point out that this behavior does not depend on any of the parameters of our theory,
i.e. $k_0$, $\mu_0$ or $\epsilon$. This means that we can redo Fig. (\ref{RGToy}) for the charge density in the
normal phase. This particular case is interesting, since the theory becomes linear and we can therefore separate
variables. Being that the case, we can recompute the power spectrum solving the equations of motion using a simple
$Mathematica$'s NDSolve command  and we get $\Gamma= -1.0+ 1.0\,\alpha$.
The slope of this fit agrees with that found both for the condensate and the charge density in the broken phase
(Fig. (\ref{RGToy}))
\footnote{A similar behavior was observed in \cite{Hartnoll:2014cua} for the spectrum of a scalar operator.}.
As for the ordinate we expect the aparent discrepancy to vanish for numerics with thinner grids.



{\it Conclusions.}-- In this paper we report two interesting findings: (i) The critical temperature of holographic superconductors increases with the increase in the strength of disorder and (ii) The power spectrum of the condensate and charge density are governed by fairly universal relations depending on the power spectrum of the original random signal.

Let us cautiously compare our results with the situation in the condensed matter literature. In the condensed 
matter literature about superconductors some results point to a degrading of $T_c$ with the strength of the disorder
\cite{Kapitulnik}, 
\cite{s-wave}. Other results, however, point to an enhancement of $T_c$ \cite{Fisher-Fisher-Huse}. The precise role 
of the interactions in these studies is hard to gauge. 
Our results clearly point to an enhancement of $T_c$ but we 
should warn the reader that the role of interactions in our context is central and direct comparison with previous 
studies in the condensed matter literature will require a considerable amount of work likely at the level of 
\cite{BAA} for dirty superconductors, that is, an analysis able to sum the electron-electron perturbation theory to 
all orders.


It might be more pertinent to compare our results with the literature for strongly coupled disordered superfluids.
Indeed the numerical simulations of \cite{Dang, Chandra} have shown
that disorder can triger an insulator to superfluid phase transition in systems that can be of relevance both for
superfluids and high-Tc superconductors.







It is worth pointing out that other holographic discussions, which could be considered as technical precursors to our work in that they solved simplified versions of our system, seem to also point to an enhancement of $T_c$. For example,  \cite{Bao:2011pa} considered a time-dependent stimulation that lead, for a range of frequencies to increasing $T_c$, see however \cite{Natsuume:2013lfa} (also \cite{Li:2013fhw}).  Spatially modulated chemical potentials considered, for example, in \cite{pajer}, \cite{Erdmenger:2013zaa} (see also \cite{Domokos:2013kha}) point instead to a reduction of $T_c$.


The universality of the result $k^{-\alpha} \mapsto (k^{-\Delta(\alpha)}, k^{-\Gamma(\alpha)})$  for the order parameter and the charge density seems to be a general property of the gravity equations of motion. We have also provided supporting evidence that this result is largely independent of the value of the mass of the scalar field. An immediate and intriguing conclusion is that the operator mostly responsible for the RG is not necessarily the one to which the scalar field couples. This ``universality'' is interesting in the framework of the AdS/CFT and deserves further investigation.


The enhancement of  $T_c$ with the strength of the noise in holographic superconductors deserves more scrutiny.  It could be an important prediction of holographic superconductivity.
It would be particularly interesting to consider other types of holographic superconductors (like, for instance, $p$-wave \cite{Gubser:2008wv}), and carry out a similar analysis there. We leave this for future investigation.

In this manuscript we focused largely on the behavior of the condensate and the charge density averaged over the $x$ direction. There seems to be a rich structure in the $x$-dependence of such quantities (see Figs. (\ref{Smoothing}) and (\ref{Roughening})). In particular, the condensate seems to show potential islands of superconductivity. It would be interesting to pursue the appearance of islands of superconductivity and the effect of different noises on the conductivity (see \cite{Hartnoll:2012rj}). We hope to address such questions in a separate publication.




{\bf Acknowledgments} We are thankful to C. Herzog, J. Maldacena, D. Musso, J. Sonner, D. Tielas, D. Vaman and C. Wu.  D.A.
and L.PZ.  acknowledge hospitality of the  MCTP  and ICTP, respectively. L.PZ. and A.S. acknowledge the hospitality
of GGI during the gestation of the project. D.A. thanks the FRoGS for unconditional support. This work was supported
by the DoE Grant DE-SC0007859.


