We present the effects of vector-like fermions (VLF) on the stability of the Higgs electroweak vacuum, using the renormalization group improved Higgs effective potential. We review the calculation of the one-loop beta-functions of the standard model couplings paying particular attention to the fermion contributions. From this, we derive the VLF contributions to the betafunctions. We also include the significant two-loop contributions to the beta-functions. Using these beta-functions, we determine the scale at which the effective Higgs quartic-coupling becomes zero and goes negative, signaling vacuum instability. We find that for certain VLF masses and Yukawa couplings, the Higgs quartic stays positive for field values all the way up to the Planck scale, implying that the meta-stable vacuum of the standard model can be rendered absolutely stable if VLFs are present with certain parameters. For other values of VLF parameters, the Higgs vacuum is metastable as in the standard model. For cases where the vacuum is metastable, we compute the probability of quantum tunneling from the false electroweak vacuum into a deeper true vacuum in our Hubble volume by numerically solving for the bounce configuration in Euclidean space-time and computing the bounce action for it. We compare our numerical solution with the analytical approximation for the bounce action commonly used in the literature, and comment on when the latter may be used. * shri@imsc.res.in † arunprasath@imsc.res.in arXiv:1812.11303v3 [hep-ph] 13 Jul 2019Denoting the field value as h ≡ µ and denoting λ eff (µ) as just λ(µ), it can be shown (see for example Ref. [41]) that λ(µ) obeys a renormalization group equation (RGE) of the formThe RGE is interpreted now as an evolution with field value h = µ, and the β-function β λ is the usual β-function for the coupling λ, governed by the RGE. The λ(µ) obtained by integrating the RGE has the leading logs of the form log n (µ/M ) resummed. β λ is shown as a function of λ itself, and also of the other couplings that contribute significantly, which, in the SM, are the top Yukawa coupling y t and the SU (3), SU (2), U (1) gauge couplings g a = {g 3 , g 2 , g 1 }. All these couplings also evolve with µ via analogous RGE equations with their corresponding β-functions β yt , β ga . We neglect the contributions of the other SM couplings to the β-functions as they contribute insignificantly. From Eq.(3) we see that for h m h , the instability is signalled by the Higgs quartic effective coupling λ(µ) becoming negative.As we show explicitly later, β λ obtains a negative contribution from y t , while it obtains a positive contribution from λ(µ) itself and from gauge couplings. Thus, the top quark has the important effect of decreasing λ(µ), and for y t as large as in the SM, for the observed m h , it drives λ(µ) negative at higher energies, signaling vacuum instability. The effect of fermions coupled to the Higgs is generally to destabilize the electroweak vacuum, although in this work we show that this statement is not so definit...