 signal the presence of an instability in the potential [@Weinberg:1987vp].[^7] We neglect the contributions from $h$, $s$ and the corresponding Goldstone bosons. We explicitly checked that this leads to shifts in physical observables of a few percent at most. We also do not take into account additional corrections coming from the anomalous dimensions of the Higgs and the scalar field, as they are typically only at the few percent level, as well.

All couplings as well as all logarithms in the effective potential depend on a renormalization scale. In (\[eq:Veff\_app\]) we introduced three scales $\mu_h$, $\mu_s$, and $\mu_{hs}$ that cancel separately up to terms suppressed by two loops. In our numerical analysis we set the renormalization scales to the corresponding field values $\mu_h = h$, $\mu_s = s$, and $\mu_{hs} = \sqrt{h 
s}$, which is expected to keep higher order corrections to the effective potential small.

Beta Functions \[sec:betafunctions\]
====================================

The one loop beta functions of the couplings of our framework read ($t = 
\log\mu$) $$\begin{aligned}
 \frac{d \lambda_H}{dt} = \beta_{\lambda_H} &=&  \beta_{\lambda_H}^{\rm SM} + 
\frac{1}{16\pi^2} 4 \lambda_{\Sigma H}^2  ~, 
\\[8pt]
 \frac{d \lambda_{\Sigma}}{dt} = \beta_{\lambda_\Sigma} &=& \frac{1}{16\pi^2} 
\Big( 12 \lambda_\Sigma^2 + 4 \lambda_{\Sigma H}^2 
- 9 g_X^2 \lambda_\Sigma - 3 (g_X^\prime)^2 \lambda_\Sigma + \frac{9}{4} g_X^4 
+\frac{3}{4} (g_X^\prime)^4 +\frac{3}{2} g_X^2 (g_X^\prime)^2
\nonumber \\ 
&& - 4 \sum_i (Y_{\xi_i}^4 + Y_{\chi_i}^4) + 4 \lambda_\Sigma \sum_i 
(Y_{\xi_i}^2 + Y_{\chi_i}^2) \Big) ~,  
signal the presence of an instability in the potential   [ @Weinberg:1987vp].[^7 ] We fail the contribution from $ h$, $ s$ and the corresponding Goldstone bosons. We explicitly see that this leave to shifts in physical observables of a few percentage at most. We also do not take into history additional corrections come from the anomalous dimensions of the Higgs and the scalar airfield, as they are typically merely at the few percent level, as well. 

 All couplings as well as all logarithm in the effective potential depend on a renormalization scale. In   (\[eq: Veff\_app\ ]) we introduced three scale $ \mu_h$, $ \mu_s$, and $ \mu_{hs}$ that cancel separately up to term suppress by two loops. In our numerical analysis we set the renormalization scales to the comparable field values $ \mu_h = h$, $ \mu_s = s$, and $ \mu_{hs } = \sqrt{h 
 s}$, which is expected to keep higher order correction to the effective potential small. 

 Beta Functions \[sec: betafunctions\ ] 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 The one loop beta functions of the couplings of our framework read ($ t = 
 \log\mu$) $ $ \begin{aligned } 
  \frac{d \lambda_H}{dt } = \beta_{\lambda_H } & = &   \beta_{\lambda_H}^{\rm SM } + 
 \frac{1}{16\pi^2 } 4 \lambda_{\Sigma H}^2   ~, 
 \\[8pt ] 
  \frac{d \lambda_{\Sigma}}{dt } = \beta_{\lambda_\Sigma } & = & \frac{1}{16\pi^2 } 
 \Big (12 \lambda_\Sigma^2 + 4 \lambda_{\Sigma H}^2 
 - 9 g_X^2 \lambda_\Sigma - 3 (g_X^\prime)^2 \lambda_\Sigma + \frac{9}{4 } g_X^4 
 + \frac{3}{4 } (g_X^\prime)^4 + \frac{3}{2 } g_X^2 (g_X^\prime)^2 
 \nonumber \\ 
 & & - 4 \sum_i (Y_{\xi_i}^4 + Y_{\chi_i}^4) + 4 \lambda_\Sigma \sum_i 
 (Y_{\xi_i}^2 + Y_{\chi_i}^2) \Big) ~,
 sihnal the presence of an lnstability in tkw poteitial [@Wejnberg:1987vp].[^7] We neglect the contributiond drom $y$, $s$ and the correspondkng Goldsnone bosobs. Wt explicitly checked that this mcads co shifts in phyxical obsesvables of a faw pzrcent at most. We also do not take igto accpujt additional sorrtctyons boning from the anomalous dimenaions oh the Higgs and the scalar field, as they wre hypically only at hhe few perxent oevel, as welu.

All couplpugs as well as all logarithms in the effecgive 'otential dg'wnd mn a renormelizatpon scale. In (\[eq:Veff\_app\]) we inttoduced three xcanes $\mu_h$, $\mu_s$, and $\mu_{hs}$ thav cancel separately op to termv auppressed by two liops. Hn ogr njnerkcam enamysis ae aet the rehormalizatiin scales to the cotrqwponding fiels valuqs $\mu_h = h$, $\mu_s = s$, and $\mu_{hs} = \sqrt{h 
s}$, which ps espected to keep higher irder corrections to jhe effectyve potential small.

Beta Functions \[sec:betafunctionv\]
====================================

The kve oojo bfta functions of the couplings of our framewowi tesd ($t = 
\log\mu$) $$\benin{aligned}
 \frac{d \lsmhds_R}{dt} = \beta_{\lambaa_H} &=&  \ystz_{\lambda_H}^{\rm SM} + 
\fraf{1}{16\pi^2} 4 \laibda_{\Sugma H}^2  ~, 
\\[8kt]
 \frsc{d \lambda_{\Sigma}}{dt} = \beta_{\lambea_\Sigma} &=& \frab{1}{16\pi^2} 
\Vig( 12 \lambda_\Sigma^2 + 4 \lambda_{\Sigmc H}^2 
- 9 b_X^2 \lakbda_\Sigma - 3 (g_X^\prime)^2 \lamyda_\Sigja + \frac{9}{4} g_X^4 
+\frac{3}{4} (g_X^\pdkme)^4 +\frac{3}{2} g_X^2 (g_X^\prkme)^2
\kongmber \\ 
&& - 4 \sum_i (Y_{\xi_i}^4 + Y_{\chi_i}^4) + 4 \lambda_\Wigmc \sum_i 
(Y_{\bi_i}^2 + Y_{\chi_i}^2) \Big) ~,  
signal the presence of an instability in [@Weinberg:1987vp].[^7] neglect the from $h$, $s$ We checked that this to shifts in observables of a few percent at We also do not take into account additional corrections coming from the anomalous of the Higgs and the scalar field, as they are typically only at few level, well. couplings as well as all logarithms in the effective potential depend on a renormalization scale. In we introduced three scales $\mu_h$, $\mu_s$, and $\mu_{hs}$ cancel separately up to suppressed by two loops. In numerical we set renormalization to corresponding field values = h$, $\mu_s = s$, and $\mu_{hs} = \sqrt{h s}$, which is expected to keep higher order to the small. Beta \[sec:betafunctions\] The loop beta functions couplings of our framework read ($t \frac{d \lambda_H}{dt} = \beta_{\lambda_H} &=& \beta_{\lambda_H}^{\rm SM} + 4 \lambda_{\Sigma ~, \\[8pt] \frac{d \lambda_{\Sigma}}{dt} = \beta_{\lambda_\Sigma} \frac{1}{16\pi^2} \Big( 12 \lambda_\Sigma^2 + 4 \lambda_{\Sigma H}^2 9 g_X^2 \lambda_\Sigma - 3 (g_X^\prime)^2 \lambda_\Sigma + \frac{9}{4} g_X^4 +\frac{3}{4} (g_X^\prime)^4 +\frac{3}{2} g_X^2 (g_X^\prime)^2 && - 4 \sum_i + Y_{\chi_i}^4) + \lambda_\Sigma (Y_{\xi_i}^2 Y_{\chi_i}^2) ~,
 signal the presence of an instAbility in tHe potEntIal [@weInbeRg:1987vp].[^7] we neglect the coNTribUtions from $h$, $s$ and the corrEsponDiNG GolDStOne boSons. We eXPlICItlY cHeCkeD tHAt This lEadS to shifTs in physicAl oBsErvables of a fEW pErcent at moSt. WE also do not taKe iNto accOuNt aDDitioNal CorreCtions COming fRom the anoMaLOus dimENsions oF THe higgS and the scalar fielD, As THey are typicallY only aT tHE fEW PerCenT level, as weLl.

all coUPlings aS WeLL AS alL Logarithms in tHe effective POteNtial dEpEnd ON a renoRmaliZaTIon Scale. In (\[eq:VeFf\_apP\]) we introdUced thREe scaleS $\Mu_h$, $\mu_s$, aNd $\mu_{hs}$ ThaT caNcel SEpArAteLy UP to TErMs sUPprEssed by tWo LoOps. In Our nUMERIcal AnaLysiS we seT the renormaliZatIon sCAleS to thE corrEspoNdIng fiEld valUes $\mu_H = h$, $\Mu_s = s$, and $\mu_{hs} = \sqrt{H 
s}$, whIch is expeCteD tO keEp HigheR Order cOrrEctIons to tHe effecTIve PoTENTiAl small.

Beta FunctioNs \[SEC:bEtafunctIons\]
====================================

ThE OnE lOOp beta fuNcTioNs of THE coupLingS Of Our frameWork reAD ($t = 
\LoG\mu$) $$\begiN{aLigned}
 \FrAc{d \LamBda_H}{dT} = \Beta_{\Lambda_h} &=&  \beta_{\lamBda_H}^{\rM sM} + 
\frac{1}{16\pi^2} 4 \lambda_{\sIgma H}^2  ~, 
\\[8pt]
 \frac{d \lAMbDA_{\siGMa}}{dt} = \BetA_{\lambda_\SigmA} &=& \fraC{1}{16\Pi^2} 
\BiG( 12 \lamBDa_\sigMA^2 + 4 \lambDa_{\SigMa h}^2 
- 9 G_X^2 \LAmbda_\Sigma - 3 (g_X^\prime)^2 \laMbDa_\SigmA + \frac{9}{4} G_X^4 
+\frac{3}{4} (g_X^\prime)^4 +\Frac{3}{2} g_X^2 (g_X^\prIME)^2
\Nonumber \\ 
&& - 4 \Sum_i (y_{\Xi_I}^4 + y_{\chi_i}^4) + 4 \lambda_\SigMa \sum_I 
(Y_{\xi_i}^2 + Y_{\chi_i}^2) \bIg) ~,  
 signal the presence of an  instabili ty in  th e p ot enti al [ @Weinberg:1987 v p].[ ^7] We neglect the con tribu ti o ns f r om  $h$,  $s$ an d  t h e  co rr es pon di n gGolds ton e boson s. We expl ici tl y checked th a tthis leads  to  shifts in p hys ical o bs erv a blesofa few  perce n t at m ost. We a ls o  do no t  take i n t oacco unt additional co r re c tions coming f rom th ea no m a lou s d imensionsof  theH iggs an d  t h e sca l ar field, asthey are ty p ica lly on ly  at the fe w per ce n t l evel, as we ll.
All coup lingsa s wella s all l ogarit hms  in  the ef fe cti ve pot e nt ial dep end on a  r en ormal izat i o n scal e.In ( \[eq: Veff\_app\])weintr o duc ed th ree s cale s$\mu_ h$, $\ mu_s$ ,and $\mu_{hs}$that  cancel s epa ra tel yup to termssup pre ssed by  two lo o ps.  I n o ur  numerical analysi sw e  s et the r enorma l iz at i on scale stothec o rresp ondi n gfield va lues $ \ mu _h  = h$,$\ mu_s =  s $,and  $\mu _ {hs}  = \sq rt{h 
s} $, wh i ch is expected to keep highe r  o r d er corr ect ions to the  eff e ctiv e po t en tia l  smal l.

B et a  F u nctions \[sec:betaf un ctions \]
== ============= ========== = = = ======== 

Th e  o n e loop beta fu nctio ns of thec ouplings  of o ur frame work read ( $t = 
\l og\ mu$ ) $ $\b e g in {aligned}
 \f r a c{d\l ambda_H }{d t} = \b eta _{\ lam bda _H } &=&  \b eta_{\la mb da _H }^ {\r m SM} + 
\frac {1 }{1 6\ pi^ 2} 4\ lambda _{\Si gmaH} ^2  ~,  
\\[8p t ]\ frac {d  \ lamb da_ {\ Sigma }}{d t } =  \beta_ {\lambda_ \Si g ma}&= &\frac{1 }{16\pi^2} 
\ Bi g( 12 \lam bd a_\ Sigma^ 2 + 4 \lam bda_{\Sigma H}^2 
- 9 g _ X^2 \la mbd a_\Si gma- 3 (g_X^ \pr ime)^2  \l a mbda_\ Sigma+ \fr ac {9} { 4 } g_X ^ 4  +\f ra c{3}{4} (g _ X ^\p rime) ^4  +\f rac{3}{ 2} g_X^2 (g_X^\pri m e)^ 2
\nonumber \ \ && - 4  \ sum _ i( Y_{ \x i _i} ^ 4  + Y_{\chi_i}^4 ) + 4 \lam bd a _\ Sigma \sum _ i (Y _{\xi_i }^2 + Y _{\ch i _i}^2)\Big) ~, 
 signal_the presence_of an instability in_the potential [@Weinberg:1987vp].[^7]_We_neglect the_contributions_from $h$, $s$_and the corresponding_Goldstone bosons. We explicitly_checked that this_leads_to shifts in physical observables of a few percent at most. We also do_not_take into_account_additional_corrections coming from the anomalous_dimensions of the Higgs and_the scalar_field, as they are typically only at the_few_percent level, as_well.

All couplings as well as all logarithms in the_effective potential depend on a renormalization_scale. In (\[eq:Veff\_app\]) we_introduced_three_scales $\mu_h$, $\mu_s$, and_$\mu_{hs}$ that cancel separately up to_terms suppressed by two loops. In_our numerical analysis we set the renormalization_scales to the corresponding field values_$\mu_h = h$, $\mu_s =_s$, and_$\mu_{hs} = \sqrt{h 
s}$, which_is expected to_keep higher_order corrections to_the effective potential small.

Beta Functions \[sec:betafunctions\]
====================================

The_one loop beta_functions of the couplings of our_framework_read ($t =_
\log\mu$)_$$\begin{aligned}
_\frac{d \lambda_H}{dt}_= \beta_{\lambda_H} &=&__\beta_{\lambda_H}^{\rm SM}_+_
\frac{1}{16\pi^2} 4 \lambda_{\Sigma H}^2  ~,_
\\[8pt]
_\frac{d \lambda_{\Sigma}}{dt} = \beta_{\lambda_\Sigma} &=& \frac{1}{16\pi^2} 
\Big(_12 \lambda_\Sigma^2 + 4_\lambda_{\Sigma_H}^2 
- 9 g_X^2_\lambda_\Sigma - 3 (g_X^\prime)^2 \lambda_\Sigma_+ \frac{9}{4} g_X^4 
+\frac{3}{4} (g_X^\prime)^4 +\frac{3}{2}_g_X^2 (g_X^\prime)^2
\nonumber_\\ 
&&_- 4 \sum_i (Y_{\xi_i}^4 + Y_{\chi_i}^4) + 4 \lambda_\Sigma \sum_i 
(Y_{\xi_i}^2_+ Y_{\chi_i}^2) \Big) ~,  