\section{Introduction}

The idea of using gravitational wave as a complementary approach to explore particle physics started some time ago e.g.~\cite{Apreda:2001us,Grojean:2006bp,Jarvinen:2009mh}. However, the bulk of the research concentrated, so far, on the electroweak phase transition which is typically in the detection region of the LISA gravitational wave detector  as nicely summarised in \cite{Caprini:2015zlo,Caprini:2019egz}. Of special interest is the possible detection of gravity waves originated in Grand Unified Theories (GUT). The interest arises also because typically the new physics energy scale of GUTs is beyond the reach of the existing and even future $100\,\rm{TeV}$ colliders.

A prerequisite to start even discussing gravitational wave detection is that the underlying theory must undergo a strong first order phase transition at some point during the evolution of the Universe. Additionally the higher is the energy scale of the first order phase transition the higher will be the peak frequency of the gravitational wave that needs to be detected. Inevitably, the upper frequency limit of the existing and planned gravitational wave detectors (roughly at order of $10^3\,\rm{Hz}$) provides an upper bound on the detectable energy scale (roughly at $10^4\--10^5\,\rm{TeV}$). In this sense, among the different types of GUTs, only the two semi-simple GUTs: Pati-Salam model \cite{Pati:1974yy} and Trinification model \cite{Babu:1985gi} satisfy this criterion. In this work, we will focus mainly on the gravitational wave signatures of the minimal Pati-Salam model. Our investigation differs from the one in  \cite{Croon:2018kqn} in which an alternative model of Pati-Salam model was considered. In that work the authors employed a rather involved matter content that featured, however, a simpler first order phase transition structure \footnote{In the work of \cite{Croon:2018kqn}, the authors try to realize the gauge coupling unification and symmetry breaking to an intermediate step (left-right model) first and thus their scalar sectors are overall more complicated. However, their first order phase transition occurs only when $SU(4)$ is breaking while in our case both $SU(4)$ and $SU(2)_R$ breaks and thus their analysis of the first order phase transition is simpler and fewer couplings are involved.}.

The Pati-Salam model of matter field unification \cite{Pati:1974yy} is a time-honoured example in which one can address the hypercharge triviality issue by embedding it in an asymptotically free theory. From a phenomenological standpoint it can be commended because it does not induce fast proton decay, and it can even be extended to provide a stable proton \cite{FileviezPerez:2016laj} while automatically providing a rationale for the existence of right handed neutrinos (see more details in a recent nice review \cite{Pati:2017ysg}).

So far, asymptotic freedom has been the well traveled route to resolve the triviality problem. An alternative route is that in which the UV theory acquires an interacting fixed point, before gravity sets in, de facto {\it saving} itself from the presence of a cutoff.  This unexplored route was opened when the first safe gauge-Yukawa theory was discovered in \cite{Litim:2014uca}.

To achieve a safe theory with a small number of colours we employ large number of matter fields techniques \cite{PalanquesMestre:1983zy,Gracey:1996he}. The first phenomenological applications of the large $N_f$ limit appeared in \cite{Mann:2017wzh} where it was first explored whether the SM augmented by a large number of vector-like fermions can have an ultra-violet fixed-point in all couplings. The full treatment appeared in \cite{Pelaggi:2017abg} and further generalized in \cite{Antipin:2018zdg}.
It was found in \cite{Pelaggi:2017abg} and later on proved in \cite{Antipin:2018zdg} that while the non-abelian gauge couplings, Higgs quartic and Yukawa coupling can exhibit a safe fixed point, the hypercharge remains troublesome. In fact, for abelian  theories  the fermion mass anomalous dimension diverges at the alleged fixed point \cite{Antipin:2017ebo} suggesting that a safe extension of the SM, like the asymptotically free counterpart, is best obtained by embedding the SM in a non-abelian gauge structure.
The first non-abelian safe PS and Trinification embeddings were put forward in  \cite{Molinaro:2018kjz, Wang:2018yer}.  However, in the minimal models, only one generation of SM fermions can be modelled, since all the Yukawa couplings are determined by the same UV fixed point value with no resulting hierarchy at low energy.  Yukawa hierarchies among three generations of SM fermions are discussed in \cite{Sannino:2019sch}.

In this work, we will start by  investigating  gravitational wave signatures emerging in   Pati-Salam extensions of the Standard Model  embedded in an asymptotically safe scenario. We use these predictions as an initial seed value to study the first order phase transition and gravitational wave signatures. Later, we will depart from the safety scenario and will  explore a more general parameter space.
Therefore, our work of studying the phase transition and gravitational wave generation is very general and applies to
both safe and non-safe embeddings of the Pati-Salam model.

We discover that the next-generation gravity waves detectors are  able to explore time honoured extensions of the Standard Model occurring at energy scales inaccessible by present and future particle physics colliders. More precisely we show that the peak frequency of the expected gravitational wave signals ranges within $0.1-10$ Hz.

The paper is organised as follows. In Section~\ref{PS-model} we introduce the Pati-Salam model while in Section~\ref{Temperature} we compute the finite temperature corrections to the relevant part of the potential of the theory. The order of the phase transition as well as gravitational waves generation and detection are studied in Section~\ref{PhaseTransition}. The predictions for the gravity waves signals stemming from the model parameters are presented in Section~\ref{GravityWaves}. We conclude in Section~\ref{Conclusions}. In the appendix we provide some detailed computations.

\section{Introducing the Pati-Salam model}
\label{PS-model}

We first briefly review the Pati-Salam embedding of the SM suggested in \cite{Molinaro:2018kjz}.

Consider the time-honored PS gauge symmetry group $G_\text{PS}$ \cite{Pati:1974yy}
\begin{equation}
	G_\text{PS} = SU(4)\otimes SU(2)_{L} \otimes SU(2)_{R}\,,\label{Pati_Group}
\end{equation}
with gauge couplings $g_4$, $g_L$ and $g_R$, respectively. Here the gauge group $SU(4)\supset SU(3)_C \otimes U(1)_{B-L}$, where $SU(3)_C$ denotes the SM QCD gauge group.
The SM quark and lepton fields are unified into the $G_{\rm PS}$ irreducible representations
\begin{eqnarray}
\begin{split}\label{fermionsLR}
	\psi_{L i} &= \left(\begin{array}{cccc} u_L^1  & u_L^2 & u_L^3 & \nu_L\\ d_L^1  & d_L^2 & d_L^3 & e_L\end{array}\right)_i \sim  (4,2,1)_i \,, \\
	\psi_{R i} &= \left(\begin{array}{cccc} u_R^1  & u_R^2 & u_R^3 & \nu_R\\ d_R^1  & d_R^2 & d_R^3 & e_R\end{array}\right)_i \sim (4,1,2)_i \,,
\end{split}
\end{eqnarray}
where $i=1,2,3$ is a flavor index.
In order to induce the breaking of $G_{\rm PS}$ to the SM gauge group, we introduce a scalar field $\phi_R$ which transforms as the fermion multiplet $\psi_R$, that is, $\phi_R\sim (4,1,2)$:
\begin{equation}
	\phi_R \; = \; \left(\begin{array}{cccc}  \phi_R^{u_1} & \phi_R^{u_2} &\phi_R^{u_3} & \phi_R^0 \\ \phi_R^{d_1}  & \phi_R^{d_2}  &\phi_R^{d_3}  & \phi_R^- \end{array} \right)\,,
\end{equation}
where the neutral component $\phi_R^0$ takes a non-zero vev, $\langle\phi_R^0 \rangle \equiv v_R$, such that $G_{\rm PS} \vRto SU(3)_C \otimes SU(2)_L \otimes U(1)_Y$.
We also introduce an additional (complex) scalar field $\Phi\sim (1,2,2)$, with
\begin{eqnarray}
	\Phi & = & \left(\begin{array}{cc} \phi_1^0 & \phi_2^+ \\ \phi_1^- & \phi_2^0 \end{array} \right) \equiv \left(\begin{array}{cc} \Phi_1 & \Phi_2 \end{array}\right)\,,
\end{eqnarray}
which is responsible of the breaking of the EW symmetry.

The  most general Yukawa Lagrangian for the matter fields $\psi_{L/R}$ is:
\begin{eqnarray}
\mathcal{L}_{\rm Yuk}^\psi & = & y \, \text{Tr}\left[\overline{\psi_L}\, \Phi\, \psi_R \right] \,+\, y_c \, \text{Tr}\left[\overline{\psi_L}\, \Phi^c\, \psi_R \right]\,+\, \text{h.c.}\,,\label{LYuk1}
\end{eqnarray}
where $y$ and $y_c$ are the Yukawa couplings for the third generation only. Note that the Yukawa couplings of the first two generations can be generated through the clockwork mechanism \cite{Sannino:2019sch}.

In the case of a self-conjugate bi-doublet field $\Phi \equiv \Phi^c$,  one obtains degenerate masses at tree-level, namely
\begin{equation}
	m_t = m_b = m_\tau = m_{\nu_\tau}\,.\label{spectrum2}
\end{equation}
 In order to separate the neutrino and top masses in Eq.~(\ref{spectrum2}), we implement the seesaw mechanism
\cite{Minkowski:1977sc,Yanagida:1979as,GellMann:1980vs,Mohapatra:1979ia} by adding a new chiral fermion singlet  $N_L\sim (1,1,1)$, which has Yukawa interaction (see e.g.~\cite{Volkas:1995yn,Molinaro:2018kjz} for more details)
\begin{equation}
	\mathcal{L}_{\rm Yuk}^N \; = \; - y_\nu \,\overline{N_L} {\rm Tr}\left[ \phi_R^\dagger \,\psi_R \right]\,+\,{\rm h.c.}
\end{equation}
In order to split the mass of top, bottom and tau lepton in Eq.~\eqref{spectrum2}, we introduce a new vector-like fermion $F\sim (10, 1, 1)$ with mass $M_F$ and Yukawa interactions (see e.g.~\cite{Volkas:1995yn,Molinaro:2018kjz} for more details):
\begin{eqnarray}\label{YukF}
	\mathcal{L}_{\rm Yuk}^F & = & y_F\, {\rm Tr}\left( \overline{F_L} \, \phi_R^T \,i \tau_2 \,\psi_R \right) \,+\, {\rm h.c.}	\label{Yuk_F}
\end{eqnarray}
All the field contents and couplings are summarized in Tab.~\ref{couplings}.
\begin{table}[t!]
\centering
  \begin{tabular}{|| l | l | l ||}
    	\hline
Gauge & Yukawa & Scalar \\ \hline
$SU(4):\,g_4$ & $\psi_{L/R}: \,y,\,y_c$ & $\phi_R:\,\lambda_{R1},\,\lambda_{R2}$\\\hline
$SU(2)_L:\,g_L$ & $N_L:\,y_{\nu}$ & portal: $\lambda_{R\Phi_1},\,\lambda_{R\Phi_2},\,\lambda_{R\Phi_3}$\\ \hline
$SU(2)_R:\,g_R$ & $F :\,y_F$ & $\Phi:\,\lambda_1,\,\lambda_2,\,\lambda_3,\,\lambda_4$\\ \hline
\end{tabular}
\caption{\small Gauge, Yukawa and scalar quartic couplings of the PS model.}
\label{couplings}
\end{table}

\section{Finite Temperature Effective Potential}
\label{Temperature}
\subsection{Tree Level Effective Potential of Pati-Salam Model}
The relevant terms in the tree level effective potential can be written as:
\begin{equation}
V_{\rm{tree}}\left(\phi_R\right)=\lambda_{R1} \trace^2\left(\phi_R^\dagger\phi_R \right)
				+ \lambda_{R2}\,\trace\left( \phi_R^\dagger\phi_R \phi_R^\dagger\phi_R \right) \,.\label{potential tree}
\end{equation}
It is important to note that we do not include any explicit mass terms in the tree level potential. The symmetry breaking in this work is induced by Coleman-Weinberg mechanism.

If we write out $\phi_R$ explictly as :
\begin{equation}
\frac{1}{\sqrt2}\left(
\begin{array}{cccc}
 \phi_{R1}+i\phi_{R2} & \cdots &\cdots & \phi_{R7}+i\phi_{R8} \\
 \phi_{R9}+i\phi_{R10}& \cdots &\cdots & v+\phi_{R15}+i\phi_{R16} \\
\end{array}
\right) \, ,
\end{equation}
where we choose the symmetry breaking direction of $\phi_R$ and thus all field components except the $\phi_{R15}$ direction are   zero.  As mentioned above, $\langle \phi_R \rangle $ triggers breaking of $G_{\rm PS} \vRto SU(3)_C \otimes SU(2)_L \otimes U(1)_Y$.
Out of sixteen scalar fields, there are nine Goldstone bosons and seven physical bosons.
Therefore, eight gauge bosons of $SU(4)$ (corresponding to QCD gluons) and one gauge field from $SU(4) \otimes SU(2)_R$
(which is simply $U(1)_Y$, a linear combination of the $U(1)_{B-L}$ from $SU(4)$ and $U(1)_R$ from $SU(2)_R$, with $Y=2I_R+B-L$) remain massless. The other nine gauge bosons of $SU(4) \otimes SU(2)_R$ (six lepto-quark, two right boson $W_R^{\pm}$ and one $Z'$) become massive.

With Eq.\eqref{potential tree}, we can construct the mass matrix of the scalar fields and obtain sixteen tree level mass eigenvalues.
These mass eigenvalues can be divided into nine Goldstone bosons with a mass $M_{\rm{Gold}}^2=v^2 \left(\lambda _{\text{R1}}+\lambda _{\text{R2}}\right)$ and seven physical Higgses, one out of which has a mass of $M_{\rm{Higgs1}}^2=3 v^2 \left(\lambda _{\text{R1}}+\lambda _{\text{R2}}\right)$ and six other Higgses with a  mass $M_{\rm{Higgs2}}^2=v^2 \lambda _{\text{R1}}$.

\subsection{Loop Level Effective Potential of Pati-Salam Model}
In this section, we will discuss the one loop contributions to the effective potential from scalar, gauge fields and fermions. The general formula is well known and can be written as:
\begin{equation}
V_{\rm{1loop}}=\sum_{i}\pm n_i\frac{m_i^4}{64\pi^2}\left(\log\left[\frac{m_i^2}{\mu^2}\right]-C_i\right)
\end{equation}
where the sum runs over the bosons $\left(+\right)$ and fermions $\left(-\right)$ and $n_i$ counts the internal degrees of freedom (d.o.f.) of each species $i$.
The symbols $m_i$, $\mu$ and $C_i$ correspond respectively to the tree level mass terms, renormalization scale and constant (equal to $5/6$ for gauge bosons and $3/2$ for scalars and fermions in Minimal Subtraction Scheme).
We define the background field as $\rho$. In the following, we write out the scalars, gauge fields and fermions contribution explicitly.

The Higgs fields contributions (7 d.o.f.) to the one loop effective potential $V_{\text{Higgs}}$ are:
\begin{equation}
\begin{split}
&\frac{1}{64 \pi ^2}\left(3 \rho ^2 \left(\lambda _{R1}+\lambda _{R2}\right)\right)^2 \left(\log \left(\frac{3 \rho ^2 \left( \lambda _{R1}+\lambda _{R2}\right) }{\mu ^2}\right)-\frac{3}{2}\right)\\
&+\frac{6}{64 \pi ^2} \left(\rho ^2 \lambda _{R1}\right){}^2 \left(\log \left(\frac{\rho ^2 \lambda _{R1}}{\mu ^2}\right)-\frac{3}{2}\right)\,.
\end{split}
\end{equation}

The Goldstone contributions (9 d.o.f.) to the one-loop effective potential $V_{\text{Gold}}$ are:
\begin{equation}
\frac{9}{64 \pi ^2} \left(\rho ^2 \left(\lambda _{R1}+\lambda _{R2}\right)\right)^2 \left(\log \left(\frac{\rho ^2 \left( \lambda _{\text{R1}}+\lambda _{\text{R2}}\right) }{\mu ^2}\right)-\frac{3}{2}\right)
\end{equation}

The lepto-quark contributions from $SU(4)$ gauge fields (6 lepto-quark $\times3$ polarization=$18$ d.o.f.) to the one-loop effective potential are:
\begin{equation}
V_{\text{lepto}}=\frac{18}{64 \pi ^2} \left(\frac{1}{4} g_4^2 \rho ^2\right)^2 \left(\log \left(\frac{g_4^2 \rho ^2}{4 \mu ^2}\right)-\frac{5}{6}\right)\,,
\end{equation}
where the tree level lepto-quark mass is given by $M_{\rm{lepto}}^2=\frac{1}{4}g_4^2v^2$,
and $g_4~(g_R) $ is the $SU(4)~(SU(2)_R)$ gauge coupling.
The gauge boson $W_R^{\pm}$ contributions (2 $W_R$ $\times3$ polarization=$6$ d.o.f.) to the one loop effective potential are:
\begin{equation}
V_{W_R^{\pm}}=\frac{6}{64 \pi ^2} \left(\frac{1}{4} g_R^2 \rho ^2\right)^2 \left(\log \left(\frac{g_R^2 \rho ^2}{4 \mu ^2}\right)-\frac{5}{6}\right)\,,
\end{equation}
where the tree level $W_R$ mass is given by $M_{W_R^{\pm}}^2=\frac{1}{4}g_R^2v^2$.

The $Z'$ boson contribution (1 $Z'$ $\times3$ polarization=$3$ d.o.f.) to the one loop effective potential $V_{Z'}$ is:
\begin{equation}
\frac{3}{64 \pi ^2} \left(\frac{1}{8} \left(2g_R^2 +3g_4^2\right)\rho ^2\right)^2 \left(\log \left(\frac{\left(2g_R^2 +3g_4^2\right) \rho ^2}{8 \mu ^2}\right)-\frac{5}{6}\right)\,,
\end{equation}
where the tree level $Z'$ mass is given by $M_{Z'}^2=\frac{1}{8} \left(2g_R^2 +3g_4^2\right)v^2$.

The neutrino singlet contribution ($4$ d.o.f.~of Dirac Fermion) to the one loop effective potential is:
\begin{equation}
V_{\rm{\nu}}=-\frac{4}{64 \pi ^2} \left(\frac{1}{2}y_{\nu}^2\rho ^2\right)^2 \left(\log \left(\frac{y_{\nu}^2\rho ^2}{2 \mu ^2}\right)-\frac{3}{2}\right)\,,
\end{equation}
where the tree level neutrino singlet mass is given by $M_{\rm{\nu}}^2=\frac{1}{2}y_{\nu}^2v^2$

On the other hand, the Yukawa coupling in Eq.~\eqref{YukF} also contributes to the potential
as (4 colours $\times4$ d.o.f. of Dirac fermion=$16$ d.o.f.):
\begin{equation}
V_{\rm{F}}=-\frac{16}{64 \pi ^2} \left(\frac{1}{2}y_{F}^2\rho ^2\right)^2 \left(\log \left(\frac{y_{F}^2\rho ^2}{2 \mu ^2}\right)-\frac{3}{2}\right)\,,
\end{equation}
with a mass term $M_{F}^2=\frac{1}{2}y_{F}^2v^2$.
All in all, the total one-loop effective potential is:
\begin{equation}
V_{\rm{1loop}}=V_{\rm{Higgs}}+V_{\rm{Gold}}+V_{\rm{lepto}}+V_{W_R^{\pm}}+V_{Z'}+V_{\rm{\nu}}+V_F\,.
\end{equation}

\subsection{Finite Temperature Effective Potential of Pati-Salam Model}
The one loop finite temperature effective potential has the following general form
\begin{equation}
V_T=\sum_i\pm n_i \frac{T^4}{2\pi^2}\int_0^\infty dy y^2\log\left[1\mp e^{-\sqrt{y^2+m_i^2/T^2}}\right] \, ,\label{finite_temperature_full}
\end{equation}
where $+n_i~( - n_i)$ corresponds to bosons~(fermions).
We can further write the thermal integral in the form of the polynomials which can significantly simplify the calculations. We focus on the integral part of Eq.~\eqref{finite_temperature_full} and define:
\begin{equation}
I_{B,F}\left(a\right)=\pm\int_o^\infty dy y^2\log\left[1\mp e^{-\sqrt{y^2+a}}\right]\,,\label{thermal_integral}
\end{equation}
where we have used $a\equiv m_i^2/T^2$. For high temperature expansions ($m_i/T\ll1$), the thermal integral can be expanded respectively for bosons and fermions as:
\begin{equation}
\begin{split}
I_B^H\left(a\right)&=-\frac{\pi^4}{45}+\frac{\pi^2}{12}a-\frac{\pi}{6}a^{\frac{3}{2}}-\frac{a^2}{32}\left(\log\left(a\right)-c_B\right)\\
I_F^H\left(a\right)&=-\frac{7\pi^4}{360}+\frac{\pi^2}{24}a+\frac{a^2}{32}\left(\log\left(a\right)-c_F\right)\,,\label{high}
\end{split}
\end{equation}
where $c_B$ and $c_F$ are respectively $c_B=3/2-2\gamma_E+2\log\left(4\pi\right)$ and $c_F=3/2-2\gamma_E+2\log\left(\pi\right)$ and $\gamma_E\approx0.5772$. For low temperature expansions ($m_i/T\gg1$), the thermal integral for both bosons and fermions can be expanded as \footnote{Note that there are typos in the expressions of low energy expansion in \cite{Li:2014wia}.}:
\begin{equation}
I_{B,F}^L\left(a\right)=- \sqrt{\frac{\pi}{2}}a^{\frac{3}{4}}e^{- \sqrt{a}}
\left(1 +  \frac{15}{8} a^{- \frac{1}{2}} + \frac{105}{128} a^{-1}\right)\,.\label{low}
\end{equation}
To include the information for both high temperature and low temperature, we need to have an expression to connect the above two expressions Eq.~\eqref{high} and Eq.~\eqref{low}. We find:
\begin{equation}
\begin{split}
I_B\left(a\right)&=e^{-\left(\frac{a}{6.3}\right)^4}I_B^H\left(a\right)+\left(1-e^{-\left(\frac{a}{6.3}\right)^4}\right) I_{B}^L\\
I_F\left(a\right)&=e^{-\left(\frac{a}{3.25}\right)^4}I_F^H\left(a\right)+\left(1-e^{-\left(\frac{a}{3.25}\right)^4}\right) I_{F}^L\\
\end{split}
\end{equation}
Thus, we have the finite temperature effective potential (without ring contributions so far) as:
\begin{equation}
\begin{split}
V_T^{\rm{tot}}=&\frac{T^4}{2\pi^2}\left(I_B\left[\frac{M_{\rm{Higgs1}}^2}{T^2}\right]+6I_B\left[\frac{M_{\rm{Higgs2}}^2}{T^2}\right]\right.\\
+&9I_B\left[\frac{M_{\rm{Gold}}^2}{T^2}\right]+6I_B\left[\frac{M_{W_R^{\pm}}^2}{T^2}\right]+3I_B\left[\frac{M_{Z'}^2}{T^2}\right]\\
+&\left.18I_B\left[\frac{M_{\rm{lepto}}^2}{T^2}\right]+4I_F\left[\frac{M_\nu^2}{T^2}\right]+16I_F\left[\frac{M_{F}^2}{T^2}\right]\right)\,.
\end{split}
\end{equation}

\subsection{Ring Contribution to the Effective Potential of Pati-Salam Model}
The general formula for the ring contributions can be written as:
\begin{equation}
V_{\rm{ring}}^i=-\frac{T}{12\pi}\left(\left[m_i^2\left(\rho\right)+\sum_{\text{bosons}~j }
\pi_i^j\left(0\right)\right]^{3/2}-m_i^3\left(\rho\right)\right)\,,\label{V_ring}
\end{equation}
where $\pi_i\left(0\right)$ denotes the corresponding thermal mass contributions to the species $i$ from the relevant
bosonic d.o.f.~$j$ (in the outside rings of the daisy diagram). To consider the ring diagram contributions to the Higgs field, for example, $\pi_{\rm{Higgs}}$ should include all the scalar field (thermal mass) contributions denoted as $\pi_{\rm{Higgs}}^{\rm{Higgs1}},\,\pi_{\rm{Higgs}}^{\rm{Higgs2}},\,\pi_{\rm{Higgs}}^{\rm{Gold}}$
as well as the gauge field contributions.
For thermal mass contributions to the scalar field from the gauge and scalar fields (i.e.~scalar field in the big central ring of the Daisy diagram), we have the following general formula for the contributions of different species $j$ in the outside ring of the daisy diagram i.e.
\begin{equation}
\label{ring_general}
\pi_{\rm{scalar}}^j\left(0\right)=\frac{1}{12}\frac{m_j^2\left(v\right)}{v^2}T^2\,.
\end{equation}
Thus, we obtain the thermal mass from the two Higgs fields and Goldstone fields respectively as:
\begin{equation}
\begin{split}
&\pi_{\rm{scalar}}^{\rm{Higgs1}}\left(0\right)=\frac{1}{4}\left(\lambda_{R1}+\lambda_{R2}\right)T^2,\quad \pi_{\rm{scalar}}^{\rm{Higgs2}}\left(0\right)=\frac{1}{12}\lambda_{R1}T^2\\
&\pi_{\rm{scalar}}^{\rm{Gold}}\left(0\right)=\frac{1}{12}\left(\lambda_{R1}+\lambda_{R2}\right)T^2\,.\label{scalar_scalar_Ring}
\end{split}
\end{equation}
Similarly, the scalar thermal mass contributions from the gauge fields are obtained in the following:
\begin{equation}
\begin{split}
&\pi_{\rm{scalar}}^{\rm{lepto}}\left(0\right)=\frac{1}{48}g_4^2T^2,\qquad \pi_{\rm{scalar}}^{W_R^{\pm}}\left(0\right)=\frac{1}{48}g_R^2T^2\\
&\pi_{\rm{scalar}}^{Z'}\left(0\right)=\frac{1}{96}\left(2g_R^2+3g_4^2\right)T^2\,.\label{gauge_scalar_Ring}
\end{split}
\end{equation}
To obtain the total thermal mass contributions to the Higgs field, we need to include all the above thermal masses i.e.~Eq.\eqref{scalar_scalar_Ring}, Eq.~\eqref{gauge_scalar_Ring} and
we have:
\begin{equation}
\begin{split}
\sum_j\pi_{\rm{scalar}}^j\left(0\right)&=\pi_{\rm{scalar}}^{\rm{Higgs1}}\left(0\right)+6\pi_{\rm{scalar}}^{\rm{Higgs2}}\left(0\right)+9\pi_{\rm{scalar}}^{\rm{Gold}}\left(0\right)\\
&+18\pi_{\rm{scalar}}^{\rm{lepto}}\left(0\right)+6\pi_{\rm{scalar}}^{W_R^{\pm}}\left(0\right)+3\pi_{\rm{scalar}}^{Z'}\left(0\right)\,.
\label{scalar_TOT}
\end{split}
\end{equation}
Note that for each scalar field d.o.f.~(either the Higgs or Goldstone bosons), it receives the same ring diagram contributions $\sum_j \pi^j_i$.
Thus, by using Eq.~\eqref{V_ring} and Eq.~\eqref{scalar_TOT}, we obtain the total ring contributions to the scalar fields in the Pati-Salam model are:
\begin{equation}
V_{\rm{ring}}^{\rm{scalar,tot}}=V_{\rm{ring}}^{\rm{Higgs1}}+6V_{\rm{ring}}^{\rm{Higgs2}}+9V_{\rm{ring}}^{\rm{Gold}}\,.
\end{equation}

Now we consider the case where the gauge fields are in the central ring of the Daisy diagram. We have the following general formulas to calculate the gauge, scalar and fermion fields contributions to the gauge thermal masses for both abelian and Non-abelian cases. For abelian case, we have:
\begin{equation}
U\left(1\right):\quad \pi^{L,S}_{\rm{gauge}}=\frac{g'^2T^2}{3}\sum_S Y_S^2,\quad \pi^{L,F}_{\rm{gauge}}=\frac{g'^2T^2}{6}\sum_F Y_F^2\,,
\end{equation}
where $L$ denotes the longitudinal thermal mass since it can be shown that that the transverse thermal mass is suppressed and $Y_S,\,Y_F$ correspond respectively to the hypercharge of relevant scalar and fermion fields. For non-abelian case, we have:
\begin{equation}
\begin{split}
SU\left(N\right):\quad \pi^{L,S}_{\rm{gauge}}&=\frac{g^2T^2}{3}\sum_S t_2\left(R_S\right),\\
\pi^{L,F}_{\rm{gauge}}&=\frac{g^2T^2}{6}\sum_Ft_2\left(R_F\right),\\
\pi^{L,V}_{\rm{gauge}}&=\frac{N}{3}g^2T^2\,,
\end{split}
\end{equation}
where $t_2\left(R_S\right),\,t_2\left(R_F\right)$ corresponds respectively to the Dynkin indices of the scalar and fermion representations, $\text{Tr}[T_R^a T_R^b] =t_2(R) \delta^{ab}$.
We obtain the total thermal mass contributions to the lepto-quark, $W_R^{\pm}$, and $Z'$ are:
\begin{equation}
\pi^{L,\rm{Tot}}_{\rm{lepto}}=\frac{5}{3}g_4^2T^2,~~\pi^{L,\rm{Tot}}_{W_R^{\pm}}=\frac{4}{3}g_R^2T^2,~~\pi^{L,\rm{Tot}}_{Z'}=\frac{4}{3}g_4^2T^2\,.
\end{equation}
When computing ring contributions for gauge fields, we use the original basis instead of the mass eigenstates.
 Thus, both $m_i^2\left(\rho\right)$ and $\sum_i\pi_i^j\left(0\right)$ are rewritten as matrices $\mathbf{M^2\left(\rho\right)}$ and $\mathbf{\Pi\left(0\right)}$ respectively rather than eigenvalues as in the above scalar case. Eq.~\eqref{ring_general} can be correspondingly modified as:
\begin{equation}
V^{\text{gauge,tot}}_{\rm{ring}}=-\frac{T}{12\pi}\Tr\left(\left[\mathbf{M^2\left(\rho\right)}+\mathbf{\Pi\left(0\right)}\right]^{3/2}-\mathbf{M^3\left(\rho\right)}\right)\,.
\end{equation}
where we include all contributions to the gauge rings and take into account only the massive gauge bosons
for the big rings.
The $\mathbf{\Pi\left(0\right)}$ is a diagonal 10-by-10 matrix with the entries of $(i,i)$ being $5 g_4^2 T^2/3$
and entries of $(j,j)$ being $4 g_R^2 T^2/3$ for $i=(1, \dots,7)$ and $j = (8,9,10)$.

In contrast, $\mathbf{M^2\left(\rho\right)}$ is a nearly-diagonal symmetric 10-by-10 matrix with
the first six diagonal elements being $g^2_4 \rho^2/4$, the seventh being $3g^2_4 \rho^2/8$, and the last three diagonal
being $g^2_R \rho^2/4$, plus two off-diagonal elements:
$ [\mathbf{M^2\left(\rho\right)} ]_{7,\,10}$$=$$[\mathbf{M^2\left(\rho\right)}]_{10,\, 7}$$= \sqrt{3/32} g_R g_4 \rho^2$.

\subsection{Complete Finite Temperature  Potential}

Now we are ready to write out the total finite temperature effective potential of the Pati-Salam model. It can be written as:
\begin{equation}
V_{\rm{tree}}+V_{\rm{1loop}}+V_T^{\rm{tot}}+V_{\rm{ring}}^{\rm{scalar,tot}}+V_{\rm{ring}}^{\rm{gauge,tot}}\,.\label{total_potential}
\end{equation}

\section{First Order Phase Transition and Gravitational Wave}
\label{PhaseTransition}

In this section, we will discuss the order of the possible early time Pati-Salam phase transition and the impact on possible gravitational wave signals.

\subsection{Strong First Order Phase Transition}
Here we focus on showing  that a strong first order phase transition can occur at around the Pati-Salam symmetry breaking scale with a sample coupling solutions shown in Tab.~\ref{sample solution_1}.
\begin{table}[t!]
\centering
  \begin{tabular}{|| l | l | l | l | l | l | l | l | l | l | l ||}
    	\hline
  $\alpha_L$ & $\alpha_R$ & $\alpha_4$  &  $\lambda_{R1}$ & $\lambda_{R2}$ & $y_F$ & $y_\nu$ \\ \hline
   0.0038 & 0.0015 & 0.0109 & 0.291 & -0.291 & 0.004 & 0.645 \\ \hline
\end{tabular}
\caption{\small This table summarizes the sample coupling solutions at the Pati-Salam symmetry breaking scale. We did not include $\lambda_1,\,\lambda_2,\,\lambda_3,\,\lambda_4,\,$$\lambda_{{R\Phi_1}},\,$$\lambda_{{R\Phi_2}},$ $\,\lambda_{{R\Phi_3}},\,y,\,y_c$ since they are irrelevant in studying the finite temperature effective potential. Note that this set of solutions is obtained from a safe UV fixed point.}
\label{sample solution_1}
\end{table}
We did not include all the couplings in  the table since the remaining couplings are irrelevant in the analysis of our effective potential. We further note that the sample solutions in Tab.~\ref{sample solution_1}  are the ones leading to an asymptotically safe extension of the Pati-Salam model.   However, we will show that the occurrence of a first order phase transition is not limited to this set of
specific values of the couplings.

The finite temperature effective potential Eq.~\eqref{total_potential} is shown in Fig.~\ref{effective_potential}. Here we have set the renormalization scale $\mu$ at $5000\,\rm{TeV}$ that is  reasonable as the lower bound on the Pati-Salam physics scale is at  $2000\,\rm{TeV}$ or so, derived from the upper limit ${\rm Br}\left(K_L\rightarrow\mu^{\pm}e^{\mp}\right)<4.7\times10^{-12}$ \cite{Ambrose:1998us}. We have also chosen the temperature $T$ to match the critical temperature i.e.~$T=T_c=2680\,\rm{TeV}$ at which the potential has degenerate minima.

A positive non-trivial (away from the origin) minimum  occurs for  $\phi_R\sim8400\,\rm{TeV}$ and it is denoted as $\phi_{Rc}$ and thus $\phi_{Rc}/T_c\sim 3.13>1$. This  shows  that the associated phase transition is a strong first order one.

\begin{figure}[t!]
\centering
\includegraphics[width=0.8\columnwidth]{finite_temperature_Effective_Potential_2}\hspace{0.07\columnwidth}
\caption{\small We plot the finite temperature effective potential by using the set of the couplings in Tab.~\ref{sample solution_1}. The renormalization scale $\mu$ is set at $5000\,\rm{TeV}$ while the temperature is chosen at $T=T_c=2680\,\rm{TeV}$ which is the critical temperature.}
\label{effective_potential}
\end{figure}

\subsection{Connection between First Order Phase Transition and Coleman-Weinberg Symmetry Breaking}

We noticed that a strong first order phase transition occurs
when spontaneous symmetry breaking happens via the Coleman-Weinberg mechanism. This is in line with the results and expectations of \cite{Sannino:2015wka}.  Of course, in other models first order phase transitions can still occur when symmetry breaking is generated via a hard negative mass square in the potential  \cite{Cline:1996mga}.

Around the finite temperature transition the Coleman-Weinberg values of the couplings reported in Tab.~\ref{sample solution_1} are such that  $\lambda_{R1} \simeq - \lambda_{R2}$ canceling each other. From the Renormalization Group (RG) flow point of view, Coleman-Weinberg symmetry breaking occurs when the RG flows of $\lambda_{R1}\left(\mu\right)+\lambda_{R2}\left(\mu\right)$ run from positive to negative flowing from the UV to the IR. The transition point (the scale $\lambda_{R1}\left(\mu\right)+\lambda_{R2}\left(\mu\right)=0$) defines the dynamical symmetry breaking scale of the Pati-Salam model which is  below $10000\,\rm{TeV}$.

\begin{figure}[t!]
\centering
\includegraphics[width=0.8\columnwidth]{Coleman_Weinberg}\hspace{0.07\columnwidth}
\caption{\small We plot the RG running of $\lambda_{R1}\left(t\right)+\lambda_{R2}\left(t\right)$ from UV to IR. The transition point (the scale $\lambda_{R1}\left(t\right)+\lambda_{R2}\left(t\right)=0$) defines the Coleman-Weinberg symmetry breaking scale of the Pati-Salam model.}
\label{Coleman_Weinberg}
\end{figure}

\begin{figure}[t!]
\centering
\includegraphics[width=0.8\columnwidth]{stream_plot_2}\hspace{0.07\columnwidth}
\caption{\small We show the stream plot of $\lambda_{R1},\,\lambda_{R2}$ where the flow direction is defined from UV to IR. The red and black plots are both the fixed point. The two green lines are the symmetry breaking lines which are defined as $\lambda_{R1}+\lambda_{R2}=0$ for $\lambda_{R2}<0$ and $\lambda_{R2}/2+\lambda_{R1}=0$ for $\lambda_{R2}>0$. The purple line is the particular RG flow corresponding to the sample solution in Tab.~\ref{sample solution_1}.}
\label{stream_plot}
\end{figure}

To gain insight it is interesting to show the symmetry breaking phenomenon via the stream plot provided in Fig.~\ref{stream_plot}.
The green line consisting of two symmetry breaking lines~($\lambda_{R1}+\lambda_{R2}=0$ for $\lambda_{R2}<0$ and $\lambda_{R2}/2+\lambda_{R1}=0$ for $\lambda_{R2}>0$) divides the plot into two phases. The right hand side of the green line corresponding to the vacuum stable phase while the left side is related to the symmetry breaking phase. In our convention the arrows point towards the infrared. The two dots correspond respectively to a saddle point (the red one)  and to an UV fixed point in both couplings. The bare couplings are meant to be fixed at some high energy scale on the right hand side of the plot. A glance at the plot shows that the only consistent way to radiatively cross the green line is by initiating the flow in the bottom right corner of the plot. One might be tempted to cross it from left to right by starting near the black dot. However this scenario would lead to an unstable potential at high energies and therefore is discarded.

Focussing on the bottom right corner there is a special asymptotically safe trajectory emanating from the red dot.  On that trajectory the theory will avoid a Landau pole and can be considered fundamental (up to gravity) in the deep ultraviolet. Another point is that the trajectory leads to a predictive infrared physics. We are also pleased to see that there is a wider region of UV bare couplings values that lead to a Coleman-Weinberg phenomenon beyond the asymptotically safe limit.

\subsection{Bubble nucleation}

The time is ripe to discuss bubble nucleation within our model. We will provide a brief review of the method and apply it to our case.

The general picture is that as the universe cools down, a second minimum, away from the origin, develops below a critical temperature. This triggers the tunnelling from the false vacuum, at the origin,  to the stable vacuum below the critical temperature. Assuming the transition to be first order, the tunnelling rate per unit volume $\Gamma\left(T\right)$ from the metastable (false) vacuum to the stable one is suppressed by the three dimensional Euclidean action $S_3\left(T\right)$ and we have~\cite{Kobakhidze:2017mru}:
\begin{equation}
\Gamma\left(T\right)=\left(\frac{S_3\left(T\right)}{2\pi T}\right)^{3/2}T^4 e^{-S_3\left(T\right)/T}\label{decay_rate}
\end{equation}
The Euclidean action has the form:
\begin{multline}
S_3\left(\rho,T\right)=4\pi\int_0^\infty dr r^2\left[\frac{1}{2}\left(\frac{d\rho}{dr}\right)^2  +V\left(\rho,T\right) \right.
\\
\left. -V\left(0,T\right)\right]\,,
\label{Euclidean_Action}
\end{multline}
where we use the difference of the potential $F\left(\rho,T\right)\equiv V\left(\rho,T\right)-V\left(0,T\right)$ to adjust the ``datum point" of the potential at zero. The bubble configuration (instanton solution) is give by solving the following equation of motion of  the action in Eq.~\eqref{Euclidean_Action}:
\begin{equation}
\frac{d^2\rho}{dr^2}+\frac{2}{r}\frac{d\rho}{dr}-\frac{\partial F}{\partial\rho}\left(\rho,T\right)=0\,,
\end{equation}
with the associated boundary conditions:
\begin{equation}
\frac{d\rho}{dr}\left(0,T\right)=0,\qquad \lim_{r\rightarrow\infty} \rho\left(r,T\right)=0\,.
\end{equation}
To find the solutions we use the so called overshooting and under shooting method. We also used the numerical package, CosmoTransitions~\cite{Wainwright:2011kj} to cross-check our results.
 For $T=2200\,\rm{TeV}$ the bubble profile is shown in  Fig.~\ref{instanton_solution}.
\begin{figure}[t!]
\centering
\includegraphics[width=0.8\columnwidth]{instanton_solution}\hspace{0.07\columnwidth}
\caption{\small We plot the bubble profile $\rho\left(r\right)$, where $T$ is chosen at $T=2200\,\rm{TeV}$ which is slightly lower than the critical temperature at $T_c=2680\,\rm{TeV}$.}
\label{instanton_solution}
\end{figure}
We can insert the bubble profile $\rho\left(r,T\right)$ into the Euclidean action Eq.~\eqref{Euclidean_Action} and thus $S_3$ will be  dependent on $T$ only.

The next step is to obtain the nucleation temperature which is defined as the temperature at which the rate of bubble nucleation per Hubble volume and time is approximately one. This means:
\begin{equation}
\Gamma\left(T\right)\sim H^4\,,
\end{equation}
where $H$ is the Hubble constant. By using Eq.\eqref{decay_rate}, we obtain:
\begin{equation}
T\ln\frac{T}{m_{pl}}\simeq -\frac{S_3\left(T\right)}{4}\,,\label{nucleation}
\end{equation}
where $m_{pl}$ is the Planck mass.
By solving Eq.~\eqref{nucleation} numerically, we find the nucleation temperature $T_n$ is around $1260\,\rm{TeV}$.
The inverse duration of the phase transition $\beta$ relative to the Hubble rate $H_{*}$ at the nucleation temperature $T_n$ is given by:
\begin{equation}
\frac{\beta}{H_*}=\left[T\frac{d}{dT}\left(\frac{S_3\left(T\right)}{T}\right)\right]\bigg\vert_{T=T_n}\,.\label{beta}
\end{equation}
We numerically obtain $\beta/H_* \simeq 183$.

Next, we will calculate another important parameter $\alpha$ which is the ratio of the latent heat released by the phase transition normalized against the radiation density:
\begin{equation}
\begin{split}
\alpha&=\frac{\epsilon}{\rho_{\rm{rad}}}=\frac{1}{\frac{\pi^2}{30}g_*T_n^4}\left(-\Delta V+T_n\Delta s\right)\\
\Delta V&=V\left(v_{T_n},T_n\right)-V\left(0,T_n\right)\\
\Delta s&=\frac{\partial V}{\partial T}\left(v_{T_n},T_n\right)-\frac{\partial V}{\partial T}\left(0,T_n\right)\,,\label{alpha}
\end{split}
\end{equation}
where $v_{T_n}$ is the vacuum expectation value of the finite temperature effective potential at the nucleation temperature, and $g_*$~(=150) is the relativistic d.o.f. in the universe.
We find $\alpha_{T_n}\equiv
\alpha(T=T_n) =0.217$.

\subsection{Gravitational Waves}

We are now have all the instruments to address the generation and potential observation of gravitational waves stemming from the Pati-Salam early times phase transition.

For the reader's benefit we provide a brief review of the ingredients needed to discuss the acoustic gravitational waves signals by following Ref.~\cite{Weir:2017wfa}. The  discussion about collision dynamics of scalar field shells and turbulence can be found in ~\cite{Weir:2017wfa} and their effects can be safely neglected in light  of being sub-leading.

 The power spectrum of the acoustic gravitational wave is given by:
\begin{equation}
h^2\Omega_{sw}\left(f\right)=8.5*10^{-6}\left(\frac{100}{g_*}\right)^{\frac{1}{3}}\Gamma_{AI}^2 \overline{U}_f^4\left(\frac{H_*}{\beta}\right)v_w S_{sw}\left(f\right)\,,\label{power spectrum}
\end{equation}
where the adiabatic index $\Gamma_{AI}=\overline{\omega}/\overline{\epsilon}\simeq 4/3$. $\overline{\omega}$ and $\overline{\epsilon}$ denote respectively the volume-averaged enthalpy and energy density respectively. $\overline{U}_f$ is a measure of the root-mean-square (rms) fluid velocity and is given by:
\begin{equation}
\overline{U}_f^2\simeq\frac{3}{4}\kappa_f \alpha_{T_n}\,,
\end{equation}
where $\kappa_f$ is the efficiency parameter and it is well approximated by
\begin{equation}
\kappa_f\sim\frac{\alpha}{0.73+0.083\sqrt\alpha+\alpha}
\end{equation}
when $v_\omega~(\text{wall speed}) \rightarrow 1$. The spectral shape $S_{sw}\left(f\right)$ is given by:
\begin{equation}
S_{sw}\left(f\right)=\left(\frac{f}{f_{sw}}\right)^3\left(\frac{7}{4+3\left(f/f_{sw}\right)^2}\right)^{\frac{7}{2}}
\end{equation}
with peak frequency $f_{sw}$ approximated by:
\begin{equation}
f_{sw}=8.9\mu\text{Hz}\,\frac{1}{v_\omega}\left(\frac{\beta}{H_*}\right)\left(\frac{z_p}{10}\right)\left(\frac{T_n}{100\,\rm{GeV}}\right)\left(\frac{g_*}{100}\right)^{\frac{1}{6}}
\end{equation}
with $z_p$ a simulation-derived factor that is of order $10$, and following \cite{Hindmarsh:2017gnf} we take it to be $6.9$.

By substituting $\alpha_{T_n}$ and $\beta/H_*$ from Eq.~\eqref{beta} and Eq.~\eqref{alpha} into the above power spectrum formula for the acoustic gravitational wave Eq.~\eqref{power spectrum}, we plot the curves of energy density against frequency (solid lines and the sample solution in Tab.~\ref{sample solution_1} is in red) in Fig.~\ref{voyager} where the coupling solutions in Tab.~\ref{parameter_space_1} are used.
We have also included the future bounds~(dashed lines) coming from planned gravitational wave detection experiments such as LIGO Voyager \cite{Evans:2016mbw,Yagi:2017zhb}, LISA \cite{Caprini:2015zlo}, TianQing\footnote{The project's name consists of two Chinese words: "Tian", meaning sky or heavens, and "Qin", meaning the stringed instrument.}\cite{Luo:2015ght}, BBO \cite{Yagi:2011wg,Thrane:2013oya}, Einstein Telescope (ET) \cite{Punturo:2010zz,Hild:2010id} and Cosmic Explorer (CE) \cite{Evans:2016mbw}. They are shown respectively in blue, cyan, orange, purple, green and magenta in Fig.~\ref{voyager}. Interestingly, we find the predicted acoustic gravitational wave signal predicted  to be within the detection region of LIGO Voyager which is planned to be operational around 2027-2028.

\begin{figure}[t!]
\centering
\includegraphics[width=0.9\columnwidth]{Bound_1}\hspace{0.09\columnwidth}
\caption{\small In this diagram, we show the bounds of the relevant future gravitational wave detections in the plot of dimensionless energy density in GWs against frequency. The bounds of LIGO Voyager, LISA, TianQing, BBO, Einstein Telescope (ET), Cosmic Explorer (CE) are shown respectively in blue, cyan, orange, purple, green and magenta.
The acoustic gravitational wave signals predicted in our Pati-Salam models are shown in black, red and brown. The red one is predicted in the asymptotically safe scenario while the black and brown ones are with the Yukawa coupling ($y_F, y_\nu$) values deviated from the safe scenario without modifying the IR SM physics. }
\label{voyager}
\end{figure}

\section{Pati-Salam driven Gravity Waves}
\label{GravityWaves}

We are now in a position to analyse in more detail the parameter space of bare couplings leading to observable gravitational waves within the Pati-Salam grand unified framework.

For convenience we start with the asymptotically safe Pati-Salam scenario that has helped us quickly identify the relevant parameter space for the occurrence of a strong first order phase transition.

\subsection{Asymptotically Safe Case}
 In this section, we discuss an asymptotically safe embedding of the Pati-Salam framework by adding a large  number of vector-like fields into the theory.  In this limit we will argue for the existence of an UV fixed point  which solves the triviality problem while yielding a highly predictive theory at lower energies.

Without further ado we introduce $N_F$ pairs of vector-like fermions charged under the fundamental representation of the Pati-Salam gauge group Eq.~\eqref{Pati_Group} with the following charge assignments:
\begin{equation}
N_F: \left(4,1,2\right)\oplus\left(4,2,1\right)\,.
\end{equation}
For simplicity, we assume that these new vector-like fermions appear at the Pati-Salam symmetry breaking scale.

\begin{table}[t!]
\centering
  \begin{tabular}{|| l | l | l | l | l | l | l | l | l | l | l | l ||}
    	\hline
   	$\lambda_1$ & $\lambda_2$ & $\lambda_3$ & $\lambda_4$ & $\lambda_{R\Phi_1}$ & $\lambda_{R\Phi_2}$ & $\lambda_{R1}$ & $\lambda_{R2}$ & $y,\,y_c$ & $y_F$ & $y_\nu$\\ \hline
   	0.13 & 0.01 & 0.03 & 0.05 & 0.10 & 0.01 & 0.34 & -0.29 & 0.53 & 0 & 0.67 \\ \hline
\end{tabular}
\caption{\small This table summarizes the UV fixed point solution for $N_F=13$
involving the bubble diagram contributions in the Yukawa and quartic RG beta functions. $y_{F}$ is asymptotically free and thus is zero at the fixed point.
}
\label{shifting UV fixed point_1}
\end{table}

Employing the large $N_F$ beta functions reported in appendix \ref{RG_beta} we can compute the RG flow connecting the UV fixed point (red dot in Fig.~\ref{stream_plot}) and the the SM in the infrared. For each $N_F\gg 1$ input, we obtain a set of UV fixed point solutions. Follow the RG flow starting from the determined UV fixed point to the electroweak scale, we can check whether it matches onto the SM.

At the PS symmetry breaking scale, we need to use matching conditions for both the gauge couplings and scalar quartic couplings. In particular, after PS symmetry breaking, the scalar bi-doublet should match the conventional two Higgs doublet model (we implement the beta functions of the two Higgs doublet model provided in \cite{Branco:2011iw}). We have searched the full parameter space in the range of $N_F\in\left(10,\,200\right)$ and find that $N_F=13$ with the UV fixed point solutions shown in Tab.~\ref{shifting UV fixed point_1} agree best with the low energy data (both the Higgs mass and the top Yukawa coupling at the electroweak scale).
We note that  $y_F$ is  asymptotically free for all viable solutions. We have therefore provided a UV safe completion of the SM \footnote{We note that even if the fixed point is not entirely established this analysis is still valid because the associated trajectories are valid for any energy scale sufficiently close to the would-be UV fixed point due to the nature of the precise results of the large $N_f$ expansion away from the fixed point.}.

 The sample solutions in Tab.~\ref{sample solution_1} are already the asymptotically safe solutions corresponding to $N_F=13$. This set is particularly interesting because:
\begin{itemize}
\item It corresponds to a possible UV safe fixed point rendering (up to gravity) our Pati-Salam model UV complete.
\item The Pati-Salam symmetry is dynamically broken through the Coleman-Weinberg mechanism below $10000\,\rm{TeV}$ (see Fig.~\ref{Coleman_Weinberg}) without adding any mass terms\footnote{This result does not depend on the existence of the fixed point but it is a welcome prediction.}.
\item Below $2680\,\rm{TeV}$
a strong first order phase transition occurs and at the nucleation temperature $T_n=1260\,\rm{TeV}$ gravitational wave signals can be generated. These are within the reach of the planned LIGO Voyager experiment detection region (see Fig.~\ref{voyager}) as well as the detection regions envisioned for the Einstein Telescope (ET), Cosmic Explorer (CE) and Big Bang Observer (BBO).
 \end{itemize}

We show the results as the red solid curve in both Fig.~\ref{voyager} and Fig.~\ref{bound_2}).

\subsection{Beyond the safe scenario}
Here, we will go beyond the safe scenario by exploring a more general parameter space able to generate testable gravitational wave signals.

We observe that the gauge couplings $g_4,\,g_R,\,g_L$ are  fixed by the Standard Model once the Pati-Salam symmetry breaking scale is chosen. In addition, when varying the quartic couplings we must ensure the presence of the Standard Model Higgs with its $125\,\rm{GeV}$  mass at the electroweak scale.
We therefore vary only the Yukawa couplings $y_F,\,y_\nu$ and the two quartic couplings $\lambda_{R1},\,\lambda_{R2}$ to satisfy this constraint.

Scanning the Yukawa coupling parameter space, we discover that  when increasing either $y_F$ or $y_\nu$ (see black row of Tab.~\ref{parameter_space_1}), the dimensionless energy density of the gravitational wave signal  increases accordingly and the peak frequency will shift slightly to the left. This is clear when comparing the black curve with the red~(safe) curve in Fig.~\ref{voyager}.

When scanning the quartic couplings parameter space, we find that the gravitational waves signal also depends on $\lambda_{R1}+\lambda_{R2}$. Varying $\lambda_{R1},\,\lambda_{R2}$ with fixed $\lambda_{R1}+\lambda_{R2}$, the dimensionless energy density and the peak of the frequency of the gravitational wave signals are roughly fixed. When increasing $\lambda_{R1}+\lambda_{R2}$ (see Brown and Grey row of Tab.~\ref{parameter_space_2}) the dimensionless energy density of the gravitational wave signal decreases accordingly and the peak frequency shifts significantly to the left with respect to the safe scenario. This can be seen  from Fig.~\ref{bound_2}.

Thus, differently from the safe scenario where the peak frequency is roughly  around $10$Hz,  going beyond the safe scenario allows for a peak of frequency ranging between $0.1$ and $10$Hz.

\begin{table}[t!]
\centering
  \begin{tabular}{|| l | l | l | l | l | l | l | l | l | l | l | l ||}
    	\hline
  & $\alpha_L$ & $\alpha_R$ & $\alpha_4$  &  $\lambda_{R1}$ & $\lambda_{R2}$ & $y_F$ & $y_\nu$ \\ \hline
   Safe &0.0038 & 0.0015 & 0.0109 & 0.291 & -0.291 & 0.004 & 0.645 \\ \hline
 Black &0.0038 & 0.0015 & 0.0109 & 0.291 & -0.291 & 0.5119 & 0.645 \\ \hline
 Brown &0.0038 & 0.0015 & 0.0109 & 0.291 & -0.291 & 0.001 & 0.001 \\ \hline
\end{tabular}
\caption{\small This table summarizes the sample coupling solutions at the Pati-Salam symmetry breaking scale. The Safe, Black, Brown represent respectively the gravitational wave energy density curves with Red, Black, Brown colours in Fig.~\ref{voyager}.}
\label{parameter_space_1}
\end{table}

\begin{table}[t!]
\centering
  \begin{tabular}{|| l | l | l | l | l | l | l | l | l | l | l | l ||}
    	\hline
  & $\alpha_L$ & $\alpha_R$ & $\alpha_4$  &  $\lambda_{R1}$ & $\lambda_{R2}$ & $y_F$ & $y_\nu$ \\ \hline
   Safe &0.0038 & 0.0015 & 0.0109 & 0.291 & -0.291 & 0.004 & 0.645 \\ \hline
 Black &0.0038 & 0.0015 & 0.0109 & 0.291 & -0.291 & 0.5119 & 0.645 \\ \hline
 Brown &0.0038 & 0.0015 & 0.0109 & 0.291 & -0.001 & 0.5119 & 0.645 \\ \hline
Grey &0.0038 & 0.0015 & 0.0109 & 0.291 & -0.093 & 0.5119 & 0.645 \\ \hline
\end{tabular}
\caption{\small This table summarizes the sample coupling solutions at the Pati-Salam symmetry breaking scale. The Safe, Black, Brown, Grey represent respectively the gravitational wave energy density curves with Red, Black, Brown, Grey colours in Fig.~\ref{bound_2}.}
\label{parameter_space_2}
\end{table}

\begin{figure}[t!]
\centering
\includegraphics[width=0.9\columnwidth]{Bound_2}\hspace{0.09\columnwidth}
\caption{\small In this diagram, we show the bounds of the relevant future gravitational wave detections in the plot of dimensionless energy density in GWs against frequency. The bounds of LIGO Voyager, LISA, TianQing, BBO, Einstein Telescope (ET), Cosmic Explorer (CE) are shown respectively in blue, cyan, orange, purple, green and magenta.
The acoustic gravitational wave signals predicted in our Pati-Salam models are shown in red, black, grey and brown. The red one is predicted in the asymptotically safe scenario while the black, grey and brown ones are with the Yukawa and Quartic coupling ($\lambda_{R1}, \lambda_{R2}$) values deviated from the safe scenario without touching the IR SM physics. }
\label{bound_2}
\end{figure}

\section{Conclusions}
\label{Conclusions}
We investigated the gravitational wave signatures stemming from the Pati-Salam model by identifying the parameter space of its couplings supporting a strong first order phase transition.

We started the analysis by employing a safe version of the Pati-Salam extension of the Standard Model and then quickly generalised to more generic situations. We find that a Coleman-Weinberg spontaneous breaking of the symmetry triggers a first order phase transition that can be observed via the next generation of gravitational wave detectors such as LIGO Voyager, the Einstein Telescope (ET) and the Cosmic Explorer (CE).

Beyond the safe scenario we notice that the Yukawa couplings $y_\nu,\,y_F$ affect mostly the gravitational wave energy density while the combination of quartic couplings $\lambda_{R1}+\lambda_{R2}$ shifts its peak frequency.

Concluding, we discover that the peak frequency of the gravitational wave signals stemming from the Pati-Salam model ranges within $0.1-10$ Hz.  Our results lead to the exciting news that the next generation of gravity waves detectors will be able to explore important extensions of the Standard Model appearing not at the electroweak scale but at much higher energy scales not accessible through present and future particle physics accelerators.

\begin{acknowledgments}
The work is partially supported by the Danish National Research Foundation under the grant DNRF:90.
WCH was supported by the Independent Research Fund Denmark, grant number DFF 6108-00623.
Zhi-Wei Wang thanks Tom Steele, Robert Mann and Steve Abel for helpful discussions and Huan Yang for recommending the relevant references.
\end{acknowledgments}


\section{Introduction}\label{Intro}

Some 220 molecular species, not counting isotopologues, have been identified in the
galactic interstellar medium (ISM). Among them are ions, neutral radicals, metallic
compounds, and complex organic molecules (COMs), whose relative abundances and isotopic ratios
vary drastically due to the type of source and environment
(see CDMS\footnote{\texttt{https://cdms.astro.uni-koeln.de/}}). Conversely,
the observation of molecular abundances offers a powerful way to
characterise the gas properties and past history. This is particularly
true for distant galaxies where we lack spatial resolution.

Distant galaxies tell us about the earlier stages of the ISM.
Heavy elements and metals appeared early on in the ISM
after a first generation of massive stars released their
nucleosynthesis products into space. Secondary elements, such as
nitrogen, mostly appeared several
billions of years later; after a second generation of stars,
low-mass longer-lived stars released their processed material
\citep{Johnson2019}. The chemical and isotopic
composition of the ISM in a galaxy at redshift $z=1$, when the Universe
was half its present age, should therefore be markedly different from
that in the Galaxy and be dominated by the products of massive stars
\citep{Muller2011}. An analysis of the molecular content and isotopic
ratios in distant galaxies, thus, offers us an opportunity to characterise
these types of sources and to check current models of galaxy chemical
evolution.

A powerful method to study the ISM of distant galaxies is through the
detection of molecular lines in absorption.  The strength of the
absorption depends on the flux of the background source and the line
opacity.
The latter can be directly obtained and, if the background
quasar image is small enough to be fully covered by the absorbing gas,
this strength is decoupled from the distance to the absorber. Therefore,
it is not diluted by the telescope beam,
allowing for the detection of low abundance species
in distant sources.

One of the most studied molecular absorbers is the system located
towards the blazar PKS\,1830$-$211 at redshift $z=2.5$
(\citealt{Lidman1999}; see \citealt{Muller2006} for a detailed
description of the system). The blazar image is gravitationally lensed
by a foreground, nearly face-on spiral galaxy at $z=0.89$ that
intercepts the line of sight \citep{Wiklind1996,Winn2002}. At radio wavelengths, lensing
gives rise to two point-like images of the blazar,
embedded in a faint Einstein ring of 1$''$ in diameter (8\,kpc at the
distance of the galaxy, \citealt{Jauncey1991}). The two bright images,
that is, only the ones that are still present at millimetre wavelengths, are
located 2\,kpc SW and 6\,kpc NE of the galaxy nucleus, respectively. The redshift of
$z$\,=\,0.89 corresponds to an age of 6.4\,Gyr (see \citealt{Wright2006}
for the assumed cosmological parameters) and hence to a lookback time
close to half the present age of the Universe.

The molecular absorption arises in two spiral arms symmetrically
unrolling on either side of the nucleus \citep{Wiklind1996}.  Several
studies using millimetre interferometers (NOEMA, ATCA, and ALMA) have
revealed the molecular richness of the SW arm \citep{Muller2006,
Muller2011, Muller2013, Muller2014a, Muller2014b, Muller2016a,
Muller2016b, Muller2017, Muller2015}. More than 45
molecular species, plus more than 15 rare isotopologues, have been
reported in those studies. Among others, they consist of COMs
(such as CH$_3$OH, CH$_3$CN, NH$_2$CHO, and
CH$_3$CHO), of hydrocarbons (l-C$_3$H, l-C$_3$H$_2$, and C$_4$H), and of
light hydrides and cations (H$_2$Cl$^+$, ArH$^+$, CF$^+$, OH$^+$,
H$_2$O$^+$, CH$^+$, and SH$^+$). \citet{Muller2011, Muller2014a}
suggest that this denotes a chemical signature similar to that of the
galactic diffuse and
translucent clouds. Furthermore, the analysis of the
light hydrides detected in this source points out a multi-phase
composition of the absorbing gas \citep{Muller2016a, Muller2016b}.

PKS\,1830$-$211 is also a strong gamma-ray emitter, with several flaring
events in recent years (see e.g. \citealt{Abdo2015}).
During 2019, flares of intense gamma-ray activity were reported
using the Large Area Telescope (LAT), one of the two instruments on board the
Fermi Gamma-ray Space Telescope \citep{Buson2019}. The preliminary analysis conducted by
the Fermi-LAT collaboration indicates that this source
has been undergoing a long-term brightening since October 2018.
Following this gamma-ray outburst, a radio-wave monitoring campaign with the 32\,m Medicina and 64\,m Sardinia
radio telescopes, at 8.3\,GHz and 25.4\,GHz, has revealed a parallel radio emission
enhancement \citep{Iacolina2019}.

We took advantage of this event to search for new extragalactic
molecular species with the Yebes 40\,m radio telescope, whose 7\,mm
receiver was recently updated with the help of the European Research
Council NANOCOSMOS project\footnote{\texttt{https://nanocosmos.iff.csic.es/}}.
In this Letter,
we report the detection of six new extragalactic molecular species towards
the SW image of PKS\,1830$-$211.  We derived their molecular column densities
and rotational temperatures, as well as those of three more
species that were detected for the first time in this source
(Section\,\ref{results}). In light of the new
identifications, we discuss the implications for the prevailing
chemistry in the source (Section\,\ref{dis}).

\section{Observations and data reduction}
\label{obs}

The observations began in April 2019 with the 40\,m diameter Yebes
radio telescope\footnote{\texttt{http://rt40m.oan.es/rt40m$\_$en.php}}, just after
the recent enhancement of the blazar radio flux \citep{Buson2019,Iacolina2019}.
They consisted in monitoring the continuum
flux and absorption line profiles in the course of commissioning
the new 7\,mm receiver (F. Tercero et al. in preparation). The
telescope \citep{deVicente2016} is located at 990\,m of altitude near
Guadalajara, Spain.

The Q band (7\,mm) HEMT receiver allows for simultaneous broad-band
observations in two linear polarisations. It is connected to 16 fast
Fourier transform spectrometers (FFTS), each covering 2.5\,GHz with a
38\,kHz channel separation. In each polarisation, this system provides an instantaneous band of 18\,GHz between 31.5\,GHz and 50\,GHz.

The 7\,mm spectra shown in this paper are the average of 12 observing
sessions, which were carried out between April 2019 and September 2019.  Due to
the low declination of the source, it could only be observed for 5 hours per day
above 15$^{\circ}$ elevation (the telescope latitude is
$+$40$^{\circ}$\,31$'$\,29.6$''$). The intensity scale was calibrated
using two absorbers at different temperatures and the atmospheric
transmission model (ATM, \citealt{Cernicharo1985, Pardo2001}).

In addition, we carried out three observing sessions at 1.3\,cm
with the K band receiver on
May 29 and 30 as well as June 2.  The goal of those observations was to confirm the
identification of C$_3$H$^+$ by observing a third rotational
transition. The covered band extended from 23.7\,GHz to 24.2\,GHz,
and the channel spacing was 30\,kHz.
For these data, the intensity
scale was calibrated using a noise diode and the ATM model.

The observational procedure was position switching with the
reference position located $-$240$''$ away in azimuth. The telescope
pointing and focus were checked every one to two hours through
pseudo-continuum observations of VX\,Sgr, a red hypergiant star close
to the target source. VX\,Sgr shows strong SiO $v$\,=\,1
$J$\,=\,1$-$0 (at 43.122\,GHz) and H$_2$O
$J_{K_a,K_c}$\,=\,6$_{1,6}$$-$5$_{2,3}$ (at 22.235\,GHz) maser
emission.
Pseudo-continuum observations consist in subtracting
the summed spectral emission from the masers from the rest
of the spectrum, while performing pointing azimuth-elevation cross drifts or a focus
scan. This is a useful technique since atmospheric emissivity fluctuations during the scans
have a small influence on the results and provide very flat continuum baselines.
In addition, we regularly checked the pointing towards the blazar
PKS\,1830$-$211 itself (coordinates \mbox{$\upalpha_{\rm
J2000}$\,=\,18$^{\rm h}$\,33$^{\rm m}$\,39.9$^{\rm s}$},
\mbox{$\updelta_{\rm J2000}$\,=\,$-$21$^{\circ}$\,03$'$\,39.6$''$}),
whose flux (NE\,+\,SW images) was $\simeq$\,10\,Jy at intermediate frequencies of the 7\,mm band.
The pointing errors were always found $<5''$ on both axes.  In order
to confirm that no spurious signals were contaminating our spectra,
the frequency of the local oscillator was changed by 20\,MHz
from one observing session to the next.

The data were reduced using the \texttt{CLASS} software of the
\texttt{GILDAS}
package\footnote{\texttt{http://www.iram.fr/IRAMFR/GILDAS/}}.  A
polynomial baseline of low order (typically second or third order) was
removed from the spectrum obtained for each session. All sessions were then
averaged together. In this paper, we adopt $z$\,=\,0.885875
($v_0$\,=\,0\,km\,s$^{-1}$ in the local standard of rest, LSR) to
convert the observed frequencies to rest-frame frequencies.
For each observed frequency, Table\,\ref{table_telescope} lists the
telescope aperture efficiency ($\eta_{\rm A}$), the conversion factor
between flux ($S$) and antenna temperature ($T_{\rm A}^*$), the half power
beam width (HPBW), the system temperature ($T_{\rm sys}$), the integration
time, and the noise root mean square (rms).

\section{Continuum emission}
\label{cont}

The intensity scales of the spectra shown in this Letter were
normalised to the total (NE\,+\,SW) continuum intensity of the source. The latter
varies with frequency as indicated (in $T_A^*$ scale) in
Table\,\ref{table_telescope}. The molecular lines presented in this
article all appear in absorption against the SW image (the NE image
absorption, which is much weaker, appears at markedly different
velocities).  In order to derive the line opacities, we had to
determine the contribution of the SW image to the total continuum flux
detected by the telescope. This can be done by assuming that the
HCO$^+$ and HCN ($J=1-0$) line absorption (see Fig.\,\ref{fig_hcn_hco})
is at its total at velocities near $v_0=0$\,km\,s$^{-1}$, which corresponds to
the arm located in front of the SW source. This assumption, which was already
discussed by \citet{Muller2008} for the $J=2-1$ lines, is supported by
the obvious saturation of the flattened absorption profiles and by the
near-stability of the saturation level despite strong blazar intensity
variations.

\citet{Muller2008} measured a roughly constant flux ratio of the NE to
SW components, $R$\,$\sim$\,1.7, over several observations between
1995 and 2007 at 3\,mm. This value is similar to that obtained by
\citet{Subrahmanyan1990} at the centimetre wavelengths. Using ATCA data from 2011, \citet{Muller2013} derived an $R$ of $\sim$1.5 at 3\,mm and
7\,mm.  Following the strategy of \citet{Muller2008}, we estimated the
NE/SW flux ratio using the  $J=1-0$ HCO$^+$ and HCN
lines. As these saturated lines
absorb 45\,\% of the total flux, the NE component (with the
residual emission of the faint Einstein ring) has to contribute 55\,\%
to the total flux. We obtain a rather constant value of
$R$\,$\sim$\,1.2 over the period from April to September 2019 (see
Fig.\,\ref{fig_hcn_hco}).  This value is close to that measured by
\citet{Muller2014a} at 250\,GHz and 290\,GHz from ALMA data from April 9 and 11, 2012, respectively, just before a
gamma-ray flare. It is, however, significantly lower than in
most of the other measurements at millimetre and submillimetre
wavelengths (see also \citealt{Marti2019}).\ This is perhaps due to a conjunction
between the blazar flux variations and the four-week delay between the
NE and SW light paths. These flux variations are thought to be
associated with plasmon ejections along the precessing jet of the
blazar \citep{Nair2005,Muller2008}.

\begin{figure*
\includegraphics[scale=0.65,angle=0]{fig_CH3SH.eps}

\includegraphics[scale=0.65,angle=0]{fig_C3H+.eps}

\includegraphics[scale=0.65,angle=0]{fig_C3N.eps}

\includegraphics[scale=0.65,angle=0]{fig_HCOOH.eps}

\includegraphics[scale=0.65,angle=0]{fig_CH2CHCN.eps}

\includegraphics[scale=0.65,angle=0]{fig_o-H2CN_4.eps}

\includegraphics[scale=0.65,angle=0]{fig_OCS.eps}

\includegraphics[scale=0.65,angle=0]{fig_SO2.eps}

\includegraphics[scale=0.65,angle=0]{fig_NH2CN.eps}
\caption{Absorption lines towards PKS\,1830$-$211 observed with the Yebes 40\,m telescope.
The species identified for the first time in an extragalactic source are highlighted in red.
The molecules highlighted in blue are those that were identified for the first time in this source, but observed previously in other extragalactic sources.
The $y$-axis is the normalised intensity to the total (NE\,+\,SW) continuum
level (see Sect.\,\ref{cont}).
We assume a $z$\,=\,0.885875 ($v_{\rm 0}$\,=\,0\,km\,s$^{-1}$ in the local standard of rest).}
\label{fig_lines}
\end{figure*}

\section{Results}
\label{results}

\begin{table*}[!ht]
\centering
\caption{Dipole moments ($\mu$), number of detected lines, table number
listing the spectroscopic and observational line parameters, rotational temperatures ($T_{\rm rot}$) and column densities ($N$) obtained with the rotational diagram (RD) analysis, and
physical parameters assumed for the \texttt{MADEX} LTE model of the new molecules detected in this work.}
\label{Table_dipole_moments}
\begin{tabular}{l c c c c c c c c c c c@{\vrule height 8pt depth 4pt width 0pt}}
\hline
\hline
 & & &  & \multicolumn{2}{c}{RD analysis} \rule[0.2cm]{0cm}{0.15cm}\ & &  \multicolumn{4}{c}{\texttt{MADEX} LTE model} \rule[0.2cm]{0cm}{0.15cm}\ & \\ \cline{5-6} \cline{8-11}

Molecule  & $\mu$   &  No. Lines & Table & $T_{\rm rot}$ & $N$\,$\times$\,10$^{12}$ & & $v_{\rm 0}$ &  $v_{\rm FWHM}$ & $T_{\rm rot}$ & $N$\,$\times$\,10$^{12}$ \rule[0.2cm]{0cm}{0.15cm}\ & Notes \\
  & [Debye]  &   &  & [K] & [cm$^{-2}$] & & [km\,s$^{-1}$] &  [km\,s$^{-1}$] & [K] & [cm$^{-2}$]  & \\
\hline
A/E-CH$_{3}$SH      &   $\mu_{a}$=1.3$^{a}$     & 3   & \ref{Table_CH3SH}    & ...                & ...                  &   & $-$1.0 &   7.0 & 5.0   &  20 &  \\
l-C$_{3}$H$^{+}$    &   $\mu_{a}$=3.0$^{b}$     & 3   & \ref{Table_l-C3H+}   & 5\,$\pm$\,1        & 1.2\,$\pm$\,0.2      &   & $-$0.6 &   6.0 & 5.0   &  1.0 &  \\
C$_{3}$N            &   $\mu_{a}$=2.8$^{c}$     & 4   &  \ref{Table_C3N}     & 5\,$\pm$\,1        & 4\,$\pm$\,1       &   &   $-$0.5   &   5.0 & 5.0   &  \,3.0 & $\dagger$   \\
t-HCOOH             &   $\mu_{a}$=1.4$^{d}$     & 5   & \ref{Table_HCOOH}    & 5\,$\pm$\,1        & 2.3\,$\pm$\,0.5   &   & $-$0.7 &   6.0 & 5.0   &  2.3 & * \\
                    &                           &     &                      & 6\,$\pm$\,1        & 3.3\,$\pm$\,0.4   &   & $-$0.7 &   6.0 & 6.0     &  3.3 & **    \\
                    &                           &     &                      &                    & 6\,$\pm$\,2       &   &  &    &    &   &***  \\
CH$_{2}$CHCN        &   $\mu_{a}$=3.8$^{e}$     & 7   & \ref{Table_CH2CHCN}  & 10\,$\pm$\,4\,\,\, & 5\,$\pm$\,2          &   & $-$0.5 &   5.0 & 10    &  5.5 &  \\
o-H$_{2}$CN         &   $\mu_{a}$=2.5$^{f}$     & 4   & \ref{Table_o-H2CN}   & ...                & ...                  &   & $-$1.0 &   5.0 & 5.0   &  4.0  & \\
\hline
OCS                 &   $\mu_{a}$=0.7$^{g}$     &  3  &  \ref{Table_OCS}     & 8\,$\pm$\,3        & 43\,$\pm$\,15        &   & $-$0.7 &   5.5 & 8.0     &  43  & \\
SO$_{2}$            &   $\mu_{b}$=1.6$^{h}$     &  3  &  \ref{Table_SO2}     & 8\,$\pm$\,1        & 25\,$\pm$\,4\,\,\,   &   & $-$0.9 &   7.0 & 8.0     &  25  &\\
O$^+$-NH$_{2}$CN          &   $\mu_{a}$=4.3$^{i}$     &  6  & \ref{Table_NH2CN}    & 30\,$\pm$\,5\,\,\, & 11\,$\pm$\,2\,\,\,   &   & $-$0.9 &   7.0 & 30    &  11 & \\
\hline

         \end{tabular}
         \tablefoot{Refs.
         (a)~\citet{Tsunekawa1989},
         (b)~\citet{Pety_2012},
         (c)~\citet{McCarthy1995},
         (d)~\citet{Kuze1982},
         (e)~\citet{Krasnicki2011},
         (f)~\citet{Cowles1991},
         (g)~\citet{Tanaka_1985},
         (h)~\citet{Patel_1979},
         (i)~\citet{Read1986}.\\
         ($\dagger$) This slight difference in $N$ (within the
uncertainty) between the two methods is mainly caused by adapting the C$_3$N spectroscopic
values in the rotational diagram in order to address the blending of
the hyperfine components (see\,Sect.\,\ref{appendix_dr}). (*) For the $K_{\rm a}=0$ ladder (see\,Sect.\,\ref{appendix_dr}). (**) For the $K_{\rm a}=1$ ladder (see\,Sect.\,\ref{appendix_dr}).
(***) Total column density (see\,Sect.\,\ref{appendix_dr}).
}
   \end{table*}

For the first time, we detected
CH$_3$SH, C$_3$H$^+$, C$_3$N, HCOOH, CH$_2$CHCN, and H$_2$CN in an extragalactic source and identified the following three more
species, also for the first time, towards PKS\,1830$-$211: OCS, SO$_2$, and NH$_2$CN.
All of those species have been detected in absorption towards the SW spiral arm of the $z$\,=\,0.89
galaxy in front of the blazar. Figure\,\ref{fig_lines} shows these detections. For each species, Table\,\ref{Table_dipole_moments}
lists its dipole moment and derived properties.

A single Gaussian profile was fitted to the strongest
component of all detected lines to obtain the observational line parameters: radial velocities ($v_{\rm 0}$), full width at half maximum (FWHM) line widths ($v_{\rm FWHM}$),
and absorption intensities ($I_{\rm \nu}$) measured from the unity on the normalised spectrum to the total (SW\,+\,NW) continuum level (see Sect.\,\ref{cont}). Then,
optical depths ($\tau_{\rm \nu}$) towards the SW line of sight were directly derived according to

\begin{equation}
\tau_{\rm \nu} = -{\rm ln} \left(1 - \frac{I_{\rm \nu}}{0.45}\right)
\label{eq_tau}
,\end{equation}
assuming an SW source covering factor of unity (see e.g. \citealt{Muller2008} and \citealt{Muller2011}; these references
also address a more general discussion regarding how to obtain optical depths of absorbing lines in this source).
This assumption does not introduce major uncertainties since \citet{Muller2014a},
using ALMA observations, conclude that the size of the SW continuum image is only 5\,$-$\,10\,\%
larger than the size of the absorbing gas at millimetre and submillimetre wavelengths.
These results, together with the spectroscopic parameters of the involved transitions are shown in Tables\,\ref{Table_CH3SH}\,$-$\,\ref{Table_NH2CN}.
We note that while the derived radial velocities are in agreement with
previous observations at 7\,mm towards the SW line of sight \citep{Muller2011,Muller2013},
most of the line widths of
the molecular lines in our data are in a range
between 5\,$-$10\,km\,s$^{-1}$. Our absorbing molecular lines
are narrower than those of these previous works by a factor of 2.
Narrower lines
may indicate either lower optical depths or different clouds intercepting the
line of sight. This variation could result from
the appearance or disappearance of discrete components in the
background continuum source \citep{Muller2008}.

As we detected several rotational transitions for many species,
we estimated their rotational temperatures ($T_{\rm rot}$)
and molecular column densities ($N$) by drawing
rotational diagrams.
The analysis assumes the Rayleigh-Jeans approximation, optically thin lines,
and local thermodynamic equilibrium (LTE). In Appendix\,\ref{appendix_dr},
we describe this method and the diagrams relative to this work (Fig.\,\ref{fig_DR}).
The resulting $T_{\rm rot}$ and $N$ are given in Table\,\ref{Table_dipole_moments}.

We note that the $T_{\rm rot}$ obtained for C$_3$H$^+$, C$_3$N, and HCOOH is close
to the temperature of the cosmic microwave background (CMB) at
$z$\,=\,0.89 ($T_{\rm CMB}$\,=\,5.08\,$\pm$\,0.10\,K,
\citealt{Muller2013}).  This result was already noted in previous
surveys
\citep{Combes1999,Menten1999,Henkel2009,Muller2011,Muller2013}.  It
implies an excitation dominated by radiative coupling with the CMB
since the gas kinetic temperature in the SW source is found from
NH$_3$ lines that are much higher ($\sim$80\,K, according to
\citealt{Henkel2008}).  It also means that the bulk of the gas has a
moderate or low density: $\leq 10^3$\,cm$^{-3}$ (e.g. \citealt{Henkel2009,Muller2013}).

It is important to note that CH$_2$CHCN, OCS, and SO$_2$ have a somewhat higher $T_{\rm rot}$
(8\,$-$\,10\,K).\ This suggests that either the observed lines arise in the denser parts
of the absorbing clouds or these molecules present
lower critical densities than those of C$_3$H$^+$, C$_3$N, and HCOOH, which are, therefore, more sensitive to collisions.

The case is more extreme for NH$_2$CN for which we find $T_{\rm
rot}=30$\,K.  Interestingly enough, the lines pertaining to the two
different $K_{\rm a}$ ladders of that latter species merge into a
single straight line in the rotational diagram after correcting for the
different statistical weights. This is typical of asymmetric top
molecules in dense environments (see Appendix\,\ref{appendix_dr}).
As we mention above, these results may indicate that either
NH$_2$CN traces a denser and hotter component of the absorbing gas
or a surprisingly low critical density for this species.

To confirm the line identifications, we produced synthetic spectra of
the putative species using the \texttt{MADEX} tool
\citep{Cernicharo2012}.  For this, we assumed LTE approximation and
adopted the $T_{\rm rot}$ and $N$ values derived from the rotational
diagrams.  In the cases of CH$_3$SH and H$_2$CN, for which we do not have enough lines to draw a rotational diagram, we fixed $T_{\rm rot}$ to
5\,K and varied $N$ until we obtained a reasonable fit to the line profiles. The
$T_{\rm rot}$ and $N$ values as well as the fitted $v_{\rm 0}$ and $v_{\rm FWHM}$ for all
species are shown in Table\,\ref{Table_dipole_moments}.  Since
\texttt{MADEX} computes the line opacities
in the absence
of a continuum background source ($\tau_{\rm \nu}$), we derived the line absorption profile
against the SW continuum image, normalised to the total (NE\,+\,SW) continuum emission by
re-arranging Eq.\,\ref{eq_tau} to obtain $I_{\rm \nu}$.
We note that the absorption by the SW spiral arm occurs at velocities of $0\pm 50$ km\,s$^{-1}$, which are quite different from those
caused by the NE arm around $-$140 km\,s$^{-1}$, and that no absorption is detected at the latter velocity for the nine
species considered here.
The resulting synthetic spectrum is overlaid over the observed lines in Fig.\,\ref{fig_lines}.
This simple model reproduces the observed lines.
The model also allows us to distinguish an additional component for C$_3$H$^{+}$ and
H$_2$CN between 0\,km\,s$^{-1}$ and 10\,km\,s$^{-1}$ that was not taken into account in our model. This component is also
seen in the lowest energy lines of CH$_3$SH and CH$_2$CHCN.
However, the weakness of this component in the mentioned lines
prevents us from further evaluating its properties.

\section{Discussion}
\label{dis}
This work
expands the inventory of extragalactic species, particularly
in the SW line of sight towards PKS\,1830$-$211
with the detection of nine new species.
With these new species, the number of molecular species towards
the SW line of sight of PKS\,1830$-$211 raises to more than 50.

The sample of molecules detected here
points to a large diversity in molecular environments:
\begin{enumerate}
\item
Carbon
chains (C$_3$N) and ion-molecule gas phase-produced $-$CN species
(H$_2$CN and CH$_2$CHCN) have been detected in dense
($\sim$10$^4$\,cm$^{-3}$) and cold ($T_{\rm K}$\,$\sim$\,10\,K)
prestellar cores \citep{Ohishi1994,Kaifu2004} of our Galaxy.
\item
CH$_3$SH and NH$_2$CN have only been detected towards hot cores
in high-mass star-forming regions
\citep{Turner1975,Linke1979,Kolesnikova2014} and low-mass protostars
\citep{Cernicharo2012b,Majumdar2016,Coutens2018}.
NH$_2$CN has also been detected in the starburst galaxies M32 and NGC253 \citep{Martin2006,Aladro2011}.
Besides prestellar cores, CH$_2$CHCN is an abundant species in hot cores \citep{Belloche2013,Lopez2014}.
\item
Although OCS, SO$_2$, and HCOOH have been detected in a large variety of environments, these molecules usually trace
hot and shocked gas in star forming regions \citep{Tercero2010,Esplugues2013,Tercero2018} .
\item
C$_3$H$^+$ traces
the edge of dense molecular clouds illuminated by UV radiation. In our Galaxy,
this species has only been unambiguously detected towards the Horsehead and the Orion Bar \citep{Pety_2012,McGuire2014,Cuadrado2015},
which are two well known photodissociation regions (PDRs).
After the first detection of C$_3$H$^+$ towards the Horsehead \citep{Pety_2012},
\citet{McGuire2013} tentatively detected C$_3$H$^+$ via absorption lines ($J$\,=1$-$0, 2$-$1) in diffuse,
spiral arm clouds along the line of sight to Sgr\,B2(N).
These authors \citep{McGuire2014}, using the CSO telescope
operating at 1\,mm, conducted a further search for C$_3$H$^+$ towards 39 galactic sources including
hot cores, evolved stars, dark clouds, class 0 objects, PDRs, HII regions, outflows, and translucent clouds.
Interestingly, they only detected this species towards the Orion Bar, a prototypical
high-UV flux, hot PDR, with a far-UV radiation field of a few 10$^4$ times the mean interstellar
field (see \citealt{Cuadrado2015} and references therein). It is possible that
the negative results of \citet{McGuire2014} towards translucent clouds were due to a bias introduced by
the high energies of the searched transitions ($E_{\rm u}$/$k_{\rm B}$\,$>$\,50\,K). Although the gas in PDRs is also
subthermally excited, these regions are usually hotter and denser ($T_{\rm K}$\,$\sim$\,150\,K and $n$(H$_{2}$)\,$\sim$\,10$^5$\,cm$^{-3}$
for the Orion Bar, see \citealt{Cuadrado2015}) than translucent clouds,
allowing for the excitation of more energetic transitions.

\end{enumerate}

A large variety of conditions coexist in sources such as SgrB2, and absorption
measurements against a bright and hot continuum source are a powerful way of revealing all types of clouds.
As we mention in Sect.\,\ref{results}, the gas intercepting the SW line of sight at $v_{\rm 0}$\,$\sim$\,0\,km\,s$^{-1}$ mostly
consists of translucent clouds with a moderate density (10$^3$\,cm\,$^{-3}$) and
with a kinetic temperature of $\sim$80\,K \citep{Henkel2009,Muller2013}.
The wide detection of cations towards the SW line of sight \citep{Muller2014b,Muller2015,Muller2016a,Muller2016b,Muller2017} indicates
that non-thermal mechanisms, such as high UV or X-ray illumination of the gas or shocks, dominate
the chemistry of this region.
The previous tentative detection of C$_3$H$^+$ in diffuse spiral arm clouds towards Sgr\,B2(N)
and its detection towards PKS\,1830$-$211
suggests that translucent clouds may produce this cation in significant abundances.
The gas-phase production of C$_3$H$^+$ via C$_2$H$_2$\,+\,C$^+$
competes with its destruction pathway via H$_2$ leading to small hydrocarbons \citep{Pety_2012}.
To observe detectable abundances of C$_3$H$^+$, the UV-radiation field should be greatly enhanced relative
to the mean interstellar value to produce sufficient C$^+$.
The detection of C$_3$H$^+$ towards the translucent clouds intercepting the SW line of sight of PKS\,1830$-$211
indicates the large average UV-field in this medium.

It is worth noting that many COMs in our Galaxy have been detected, for the first time,
in absorption and in the centimetre domain
pointing towards the hot cores of Sgr\,B2 (see e.g. CH$_2$CHCHO, CH$_3$CH$_2$CHO, CH$_2$OHCHO,
c-H$_2$C$_3$O, CH$_3$CONH$_2$, CH$_3$CHNHO, HNCHCN, and CH$_3$CHCH$_2$O;
\citealt{Hollis2004a,Hollis2004b,Hollis2006a,Hollis2006b,Loomis2013,Zaleski2013,McGuire2016}).
These detections via absorption lines indicate that these species are associated with
the moderate dense ($n$(H$_2$)\,$\sim$\,10$^3$\,$-$\,10$^4$\,cm$^{-3}$)
and warm ($T_{\rm K}$\,$\simeq$\,100\,$-$\,300\,K) envelope of Sgr\,B2.
This envelope around the hot cores near the galactic centre also consists of a subthermally excited
molecular gas \citep{Huettemeister1995,deVicente1997,Jones2008,Etxaluze2013}.
This suggests that both the heating of this envelope and the production of these complex species
might be mainly related to
non-thermal processes, such as shocks
or an enhanced UV or X-ray flux in the surrounding medium.
Thus, we may consider that similar mechanisms may produce the diversity of
molecular species detected towards the translucent clouds in the SW line of
sight towards PKS\,1830$-$211.

\begin{acknowledgements}
We would like to thank the anonymous referee for a helpful
report that led to improvements in the paper.
We thank the ERC for support under grant ERC-2013-Syg-610256- NANOCOSMOS. We also thank
the Spanish MINECO for funding support under grants AYA2012-32032 and FIS2014-52172-C2,
and  the CONSOLIDER-Ingenio programme $``$ASTROMOL$"$ CSD 2009-00038.
\end{acknowledgements}

\bibliographystyle{aa}