fixed, only the Gamma_i matrices, i.e. the Yukawa couplings of the second Higgs doublet, are required. We choose as free parameters the Gamma_i matrices, while the other Yukawa couplings, the Delta_i matrices, are deduced from the matching with the observed fermion masses. Conventionally, the indices of the elements of these Yukawa matrices refer to the generations of the SU(2) doublet and singlet, respectively The 2HDM Lagrangian implemented in !FeynRules is based on the Standard Model default implementation, where the scalar potential and Yukawa interactions have been modified as explained above. An important feature of this model is the freedom to redefine the two scalar fields using arbitrary "horizontal" U(2) transformations acting on the two doublets simultaneously since this transformation leaves the gauge-covariant kinetic energy terms invariant. Since a given set of Lagrangian parameter values is only meaningful for a given basis, let us take advantage of this invariance property to select the Higgs basis (by defining the additional file HiggsBasis.fr) where only one of the two Higgs fields acquires a non-zero vev, namely H_1. Note that the Higgs basis is not univocally defined since a phase reparametrization of H_2 leaves the Higgs basis condition invariant, so that the phase of H_2 can be fixed in such a way that lambda_5 becomes real. Other basis choices can in principle be easily implemented as different extension files for the main Lagrangian file Lag.fr. There are two independent minimization conditions for general 2HDM potential, one relating m_1 to lambda_1 and one relating m_3 to lambda_6. This reduces the number of free parameters in the most general 2HDM to ten (seven real parameters, three complex ones and three minimization conditions). Besides the usual three massless would-be Goldstone bosons, the physical spectrum also contains a pair of charged Higgs with a mass directly related to lambda_3 and m_2, so that m_2 can be directly extracted from this mass, given as an external input. The symmetric matrix squared mass matrix for the three neutral Higgs field is diagonalized by an orthogonal matrix T which describe the relation between the physical scalar fields and the doublet neutral components. Even though this matrix is directly related to the potential parameter, it is still considered as an external input in the current implementation and must be provided by the user. In the current implementation of the 2HDM into FeynRules, the user has to provide numerical values for all the lambda_i parameters in the Higgs basis, together with the charged Higgs mass. The other parameters of the potential, such as the \mu_i, are then deduced from these inputs. Contrary, the T matrix must be calculated externally. This, together with the change of basis required if the user wants to provide potential parameters and Yukawa coupling values in bases different from the Higgs basis, as it is often the case, can be done using the ''TwoHiggsCalc'' calculator introduced for the original implementation of the 2HDM in MadGraph. This code has been modified to produce a parameter file compatible with the present implementation. This calculator can also be used to calculate the required Higgs boson tree-level decay widths. It can be found at http://cp3wks05.fynu.ucl.ac.be/Calculators/TwoHiggsCalc/index.html === References === * " ''CP violation'' ", G. Branco, L. Lavoura and J. P. Silva, Clarandon Press, Oxford, 1999. Chapter 22. === Model files === * [http://europa.fyma.ucl.ac.be/~duhr/twiki/pub/Main/2HDM/2HDM.tar.gz 2HDM.tar.gz]: This archive contains all the model files. Should be expanded in the FR model directory. * [http://europa.fyma.ucl.ac.be/~duhr/twiki/pub/Main/2HDM/HBreal.rst HBreal.rst]: Additionnal restriction file to define all parameters as being real, to be used, for example, with !CalcHEP? === Instructions === The model-file is loaded in the usual way. Feynman gauge is not supported, only unitary gauge is available. === Examples === * [http://europa.fyma.ucl.ac.be/~duhr/twiki/pub/Main/2HDM/2HDM.nb 2HDM.nb] : This is an example Mathematica® notebook that loads the model and calculates Feynman rules. * [http://europa.fyma.ucl.ac.be/~duhr/twiki/pub/Main/2HDM/param_card_FR.dat param_card_FR.dat] : This is an example of LHA parameter file to be used with this model. It also contains the defalut values used for validation. === Validation === * [http://europa.fyma.ucl.ac.be/~duhr/twiki/pub/Main/2HDM/PSPoints.jpg Phase space point comparison with stock version] * [http://europa.fyma.ucl.ac.be/~duhr/twiki/pub/Main/2HDM/TANU.jpg Cross section comparison, for tau/nu_tau initial state] * [http://europa.fyma.ucl.ac.be/~duhr/twiki/pub/Main/2HDM/TATA.jpg Cross section comparison, for tau/tau initial state] * [http://europa.fyma.ucl.ac.be/~duhr/twiki/pub/Main/2HDM/TATAbis.jpg Cross section comparison, for tau/tau initial state (part 2)] * [http://europa.fyma.ucl.ac.be/~duhr/twiki/pub/Main/2HDM/VV.jpg Cross section comparison, for VV initial state] * [http://europa.fyma.ucl.ac.be/~duhr/twiki/pub/Main/2HDM/VSSS.jpg Cross section comparison, for VSSS external legs]