This is the source code for the paper “Strong First Order Electroweak Phase Transition from a Naturally Light Singlet Scalar”. All the plots in the paper are run by the Main Mathematica file.
This section talks about some fundamental features of this new model. The problems talked here are related to new physics at zero temperature.
The scalar mass spectrum are calculated in Mathematica code mass_spectrum.nb. Masses for other relevant particles such as gauge bosons and fermions are trivial, and did not included in any code.
We need to transfer from the bare parameters to observable physical parameters. This process is performed by Mathematica in the same code.
The vacuum stability is calculated with the modified version of simplebounce code. The model file is written in benchmark.cc.
The phase transition is calculated both semi-analytically using the high-temperature expansion and in full numerical analysis.
The high-temperature expansion of the scalar effective potential offers a convenient way to perform analytical computation. The analytical calculation of high-T effective potential are calculated in Mathematica code. Plots can also be find in the same Mathematica code.
Near the boundary where the critical temperature is ill-defined, the critical temperature quickly diverges up and plots suffer from numerical fluctuation.
The Higgs and the new light scalar are mixed, leading to new signal. We calculated the Higgs exotic decay as well as the new scalar production and detection.
Scalar coupling with SM particles under our convention are calculated by hand. The collider bounds are derived from the references described in the paper. The code for the LEP bound of scalar production is in LEP.nb. The one for the Higgs exotic decay search from LHC can be found in Higgs_exotic.nb.