harvmac.tex =[ ]{} \#1\#2[\#1\#2==]{}\#1

**Abstract**

\#1

=cmr10 =cmr7 =cmr5 =cmmi10 =cmmi7 =cmmi5 =cmsy10 =cmsy7 =cmsy5 =cmti10 =’177 =’177 =’177 =’60 =’60 =’60 \#1[\#1\^]{} \#1[[e]{}\^[\^[\#1]{}]{}]{} \#1\#2[\^[\#1]{}\_[\#2]{}]{} /\#1[/]{} \#1 \#1\#2 \#1[\_[[\#1]{}]{}]{} \#1\#2[\_[[\#1]{}]{}\_[[\#2]{}]{}]{} \#1\#2[[\#1\#2]{}]{} \#1\#2\#3[[\^2 \#1\#2 \#3]{}]{} |\#1 \#1[\#1 ]{} \#1[\#1|]{} \#1[| \#1]{} \#1[| \#1|]{} \#1 \#1[1.5ex-16.5mu \#1]{} \#1\#2 \#1[.3ex]{} \#1\#2\#3[[Nucl. Phys. ]{} B[\#1]{} (\#2) \#3]{} \#1\#2\#3[[Phys. Lett. ]{} [\#1]{}B (\#2) \#3]{} \#1\#2\#3[[Phys. Rev. Lett. ]{} [\#1]{} (\#2) \#3]{} \#1\#2\#3[[Phys. Rev. ]{} [\#1]{} (\#2) \#3]{} \#1\#2\#3[[Ann. Phys. ]{} [\#1]{} (\#2) \#3]{} \#1\#2\#3[[Phys. Rep. ]{} [\#1]{} (\#2) \#3]{} \#1\#2\#3[[Rev. Mod. Phys. ]{} [\#1]{} (\#2) \#3]{} \#1\#2\#3[[Comm. Math. Phys. ]{} [\#1]{} (\#2) \#3]{}

å[ ]{} \#1[\#1-0.5em /]{}

epsf \#1\#2 $${\left[}
\def$$[\]]{} ${\left(}
\def$[)]{} 1.in

[**Submitted to Physics Letters B**]{}

A. G. Cohen[^1][Email: cohen@andy.bu.edu, anelson@ucsd.edu]{} [^2][DOE Outstanding Junior Investigator]{}

[*Physics Department*]{}

[*Boston University*]{}

[*Boston, MA 02215*]{}

A. E. Nelson[^3][Sloan fellow, SSC fellow]{}

Despite the success of the standard model of weak and strong interactions, we still remain ignorant of the mechanism of electroweak symmetry breaking; we are just beginning to probe this sector directly through Higgs searches at LEP. We are also still ignorant of the origins of the CP violation observed in the kaon system. Unfortunately, the minimal standard model accounts for all solid experimental results observed to date, and we have few experimental constraints on the symmetry breaking or CP violating aspects of the theory. There is a strong cosmological argument that the minimal standard model cannot be the whole particle physics story. The baryon to entropy ratio of the universe is $(0.4 - 1.0)10^{-10}$ \[\], and explaining this observable requires baryon number violation coupled with out-of-equilibrium CP violation in the early universe \[\]. While the standard model does have CP violation, the effects of the CP violating phase in the Kobayashi-Maskawa matrix are too suppressed by small masses and mixing angles in order to account for the observed baryonic asymmetry of the universe (BAU) . The standard model also contains anomalous baryon number violating interactions \[\], which should be rapid enough at high temperatures to affect cosmology . Furthermore, the standard model can satisfy the out of equilibrium condition for baryogenesis if the phase transition is first order, proceeding via nucleation and expansion of bubbles of the broken phase . Unfortunately it is necessary for the vev of the Higgs field after the transition to be large in order to avoid washing out any baryon number created during the transition . In the minimal standard model this requirement cannot be satisfied unless the Higgs is lighter than $\sim35$ GeV , which conflicts with current experimental bounds. Thus we should look beyond the standard model in order to discover the origin of the BAU. However we may not need to look very far, as possible mechanisms for baryogenesis have been suggested in several reasonable extensions such as axion models with additional light scalar doublets \[\], the singlet Majoron model \[\], the two Higgs model , the supersymmetric standard model , extended supersymmetric models \[\] and left-right symmetric models \[\]. Requiring that sufficient BAU be generated during the weak transition, and that it not be washed out afterwards can give us new information about the CP violating and symmetry breaking sectors of the weak interactions, allowing us to rule out some models (such as the minimal standard model), and to constrain parameters in others. It is the aim of this letter to use baryogenesis to find new constraints on supersymmetric models. We are able to severely constrain almost all unknown parameters of the minimal model, while considerable freedom remains for models with additional gauge singlets.

First let us consider baryogenesis in the minimal supersymmetric standard model (MSSM) \[\]. Several authors have claimed that this model is ruled out for baryogenesis because the bound on the lightest scalar mass in this model is the same as in the minimal standard model . However Myint  finds that the bounds on the scalar mass are relaxed somewhat, to 64 GeV, due to top quark and squark corrections to the high temperature effective potential when the top quark is heavy and the squark masses are not too heavy. The bound is relaxed because top squarks can play the role envisioned by Anderson and Hall for a gauge singlet scalar, whose coupling to the Higgs doublet increases the strength of the transition \[\]. The upper bound on the Higgs mass was computed in one loop perturbation theory and higher order corrections to the gauge propagator will reduce this number by a factor of about $\sqrt{2/3}$ , to $\sim 50$ GeV. Furthermore this bound will receive corrections proportional to $(m_H/m_Z)^2$ from two loop corrections. Thus avoiding baryon number washout is difficult in the MSSM, but not impossible. To increase the upper bound on the Higgs mass as much as possible, one should take the top quark to be heavy ( $ >  150$ GeV), the squark masses not much heavier than 150 GeV, the trilinear soft supersymmetry breaking terms (“A-terms”) small, and the parameter $\tan\beta$ (the ratio of the two Higgs vevs) less than 1.7. Then, in improved one loop perturbation theory, avoiding baryon number washout requires the lightest scalar mass to be lighter than $\sim 50$ GeV, which is not in conflict with current bounds provided that higher order corrections do not decrease the mass bound.

The MSSM has two possible sources of CP violation which are absent in the minimal standard model \[\]; a combination of these will be constrained by the BAU. The interactions in this model are given by Note that unlike the generic two Higgs doublet model discussed in refs. , the Higgs potential does not contain any CP violating phases at tree level—the phases $\phi_A$ and $\phi_B$ occur only in interactions involving the super-partners of the ordinary particles. We find that the contribution of the phase $\phi_A$ to the BAU is small, and so we focus on the effects of the phase $\phi_B$, which appears in the mass matrices of the supersymmetric partners of the gauge bosons and the Higgs scalars (the “[*inos*]{}”).

Our strategy for calculating the BAU in the MSSM is as follows. We will fix the top mass at 170 GeV, the lightest scalar mass at 48 GeV, the soft supersymmetry breaking top squark masses at 150 GeV, the A-terms at zero, and $\beta$ at 0.85 ($\tan\beta=1.14$), since these parameters are already constrained to be near these values , and in any case allowing them to vary will only affect the BAU by $\CO(1)$. Using these values we then calculate the critical temperature, the effective potential at the critical temperature, and the shape of the bubble walls. (Note that these phase transition quantities are not very sensitive to other supersymmetric parameters such as masses.) We use the improved one loop approximation for the effective potential, including the order $T^2$ corrections to the gauge propagators, and neglecting the contribution of Higgs doublet loops. When the mass of the lightest scalar is far below the gauge boson masses this approximation is reliable at the critical temperature around the symmetry breaking minimum and in the vicinity of the symmetric minimum for scalar field expectation values larger than $\CO(g_{\rm wk}T)$. Although perturbation theory is not valid for a calculation of the effective potential between the two minima, which will affect the width and shape of the phase boundary, fortunately our calculation of the BAU will turn out to be insensitive to the detailed shape of this boundary, provided it is much thicker than $\CO(1)/T$. The one loop estimate gives the width of the phase boundary to be $11/T$ at the temperature where the two minima are degenerate, so we will assume the phase boundary is thick. The one loop effective potential indicates that the ratio of the expectation values of the two Higgs remains constant during the transition, so we take where $\vev{H'}$ and $\vev{H}$ are the expectation values at the critical temperature in the symmetry breaking minimum. The temperature of the transition is 59 GeV, and at this temperature the minimum of the effective potential occurs at $H=63$ GeV and $H'=53$ GeV.

For thick bubble walls the relevant baryogenesis mechanism is known as “spontaneous baryogenesis”, reviewed below. This mechanism, first introduced for baryogenesis during a second order phase transition \[\], involves a space-time dependent field, which evolves coherently during the transition. This time evolution produces a CPT-violating term in the effective Hamiltonian called a “charge potential”, which resembles a chemical potential. A charge potential will cause the free energy density inside the bubble walls to be minimized for nonzero baryon number, and hence the production of a net baryon number via anomalous weak interactions. Spontaneous electroweak baryogenesis during the weak phase transition has been suggested for the two Higgs model in refs.  and as a baryogenesis mechanism for supersymmetric models in refs. .

In the MSSM the mass matrices are space-time varying during the transition, and contain an irremovable CP violating phase, which leads to a charge potential for baryon number in the effective fermion Hamiltonian. If the are not too heavy and the phase is not too small, this charge potential is large enough to result in generation of an acceptable baryon number during the weak transition. We will find that the resultant BAU depends sensitively on the CP violating phase $\phi_B$, and on the masses of the . We are able to use the BAU to place upper bounds on masses and lower bounds on $\phi_B$ and the electric dipole moment of the neutron.

Before launching into the specifics of the MSSM calculation, we show how to calculate the charge potentials resulting from space-time varying fermion mass matrices.

Consider a fermion mass term of the form where we take all fermions to be left-handed, and $C$ is the charge conjugation matrix . We can make a space-time dependent unitary change of basis on the fermions: in order to make the fermion mass terms everywhere real, positive and diagonal; however the space-time dependence of $U$ requires that we replace the kinetic energy terms in the Lagrangian by Note that since $U$ is a unitary matrix, $U^\dagger\partial_\mu U$ may be written where for $n$ fermions the $T_a$s are generators of $U(n)$, and the functions $\alpha(x_\mu)$ are defined by eq. . Thus the Lagrangian with mass term  is equivalent to a Lagrangian with a real diagonal mass term but also containing a term If the transformation  has a gauge anomaly there will also be a modification of the Lagrangian where the $t_\beta$s are gauge generators in the left-handed fermion representation and the $F_\beta$s are the gauge field strengths.

The presence of these anomalous terms complicates the discussion of the charge potentials, and consequently we will choose our unitary transformation to have no gauge anomaly; this is the strategy followed in .

Finally, if the transformation  does not correspond to a symmetry of the interactions the coupling constants will be affected by the change of basis; however the effect of these coupling constants on energy levels is higher order in perturbation theory and will not concern us here.

What effect does a charge potential have on a thermal system? There are two possibilities:

[1)]{} If there is a charge potential for an exactly conserved charge, (e.g. electric charge or B-L) we can ignore it. Although it looks like the system could lower its free energy by producing a net charge, charge conservation imposes a zero charge constraint. It can easily be seen by integration by parts that a charge potential for an exactly conserved charge has no physical effect; equivalently, we can always redefine fields so such charge potentials never arise.

[2)]{} Charge potentials for non-conserved charges will lead to an asymmetry in the rates between processes which create and destroy the charges, until the system reaches thermal equilibrium. The thermal occupation numbers in general will be different for particles and their CP conjugates, due to CP violation and the dynamical CPT violation from the space-time varying scalar fields. For a charge which is approximately (but not exactly) conserved the system will take a long time to reach equilibrium. If the system is near thermal equilibrium the net rate of charge production can be computed using thermodynamic arguments. For instance if there is a small charge potential $\dot\alpha_B$ for the baryon number current and no other charge potentials, for a system starting with no net quantum numbers the constraint of zero net baryon number can be implemented by introducing a baryon chemical potential The chemical potential $\mu_B$ is just the force of constraint on the system, [*i.e.*]{} the derivative of the free energy with respect to baryon number. Now the difference between the rates of anomalous processes which create and destroy baryon number is just proportional to the difference in the change in the free energy per event, leading to an anomalous baryon creation rate of where $\Gamma_B$ is the rate of anomalous baryon violating events per unit volume. The factor of 9 comes about because each anomalous event changes the free energy by $3 \mu_B$ and changes the baryon number by 3 units.

We can now calculate the BAU produced in the MSSM during the weak phase transition. We first compute the charged and neutral mass matrices, as functions of the scalar field expectation values, assuming equ. , and find the transformation on the fields $U_I(  H )$ which renders their masses real, positive, and diagonal. $U_I$ will in general be anomalous and would, by itself, give rise to an effective $W\tilde{W}$ operator, as well as giving rise to charge potentials for the via equ. . By further making a space-time dependent baryon number rotation of the light fermions we can remove this anomalous operator, at the cost of introducing a charge potential for baryon number: We then need to know how anomalous processes inside the bubble walls are affected by CP violating terms. Since the can quickly come into equilibrium with the charge potentials via ordinary [*non-anomalous*]{} interactions the charge potentials will not affect anomalous processes. The main effect on anomalous baryogenesis will come from the term . If the system is near thermal equilibrium then eqs. , and  can be used to find the total baryon number density produced during the transition We still need to know the rate of anomalous baryon density production $\Gamma_B$, which has been estimated in the symmetric phase to be \[\] while in the broken phase it is vanishingly small. Unfortunately, there is currently no reliable way of computing $\Gamma_B$ inside the wall where the scalar expectation values are changing. Furthermore it has been claimed, based on some $1+1$ dimensional simulations \[\] that anomalous baryon production inside the wall is an inherently non-equilibrium process. (Dine suggests this is unlikely to be the case for a transition with thick walls and small latent heat such as occurs in the MSSM \[\].) We will use a conservative estimate for the maximum baryon number produced by simply computing the integral  using equ.  for $\Gamma_B$. We think it is more likely that the baryon production inside the wall will be suppressed when the Higgs vevs become large. For instance McLerran \[\] estimates that baryon production is shut off when the Higgs vev reaches about (1/2) its value inside the wall, while Dine estimate that this shutoff occurs for vevs of order $g_{\rm wk} T/(4 \pi)$. In Figure 1 we plot the baryon number produced as a function of the value of the Higgs vev where this shutoff occurs for a typical choice of parameters; note that the dependence is approximately quadratic. We conclude that it is most likely that the actual baryon number produced during the transition will be between $\sim 1/4$ and $\sim 10^{-2}$ times smaller than our most favorable estimate.

In order to find the allowed range of parameters in the MSSM, for each value of the gaugino and higgsino masses we adjust the CP violating phase $\phi_B$ to be as large as is compatible with current constraints on the EDMN \[\]. We then compute the BAU, from equs.  and . The allowed range of masses are those for which the upper bound on the baryon to entropy ratio is greater than $0.4\times 10^{-10}.$ For masses in the allowed range, we compute the minimum value of $\phi_B$ which is consistent with the BAU, and use the calculation of Kizukuri and Oshimo in ref.  to find the lower bound on the EDMN.

In Figures 2a and 2b we plot the lower bound on the EDMN, as a function of $\mu$ and the wino mass parameter $m_2$. The first figure has a phase $\phi_B$ near $0$, while the second figure has the phase near $\pi$. The bino mass parameter $m_1$ is assumed to satisfy the GUT relation $m_1=m_2(5/3)\tan^2\theta_W$. The signs in fig. 2b indicate the sign of the EDMN, while in fig. 2a the EDMN is negative. The black excluded region corresponds to an EDMN greater than the current bound of $10^{-25}$ e-cm \[\], the dark grey region is an EDMN greater than $10^{-26}$ e-cm, and the light grey region greater than $10^{-27}$ e-cm. Thus for any portion of the allowed parameter space the EDMN is greater than $10^{-27}$ e-cm. (Similarly, if we assume the selectron mass is equal to the squark mass, the electric dipole moment of the electron is greater than $\CO(3) \times 10^{-29}$ e-cm.) Throughout all of the region where the EDMN could be less than $10^{-26}$ e-cm the lightest chargino mass is lighter than 88 GeV and the lightest neutralino mass is lighter than 44 GeV. Also, if the EDMN is less than $10^{-26}$ e-cm the phase $\phi_B$ must be near $\pi$ in order to avoid a chargino mass lighter than $M_Z/2$.

These figures have assumed that the rate of baryon violation is given by equ. , and that this violation is occurring throughout the bubble wall. If this is not the case, the baryon number produced will be reduced and the resulting phase must then be larger to compensate, giving rise to larger dipole moments.

If the experimental bound on the EDMN is pushed down by an order of magnitude, or if numerical calculation finds the rate of anomalous baryon creation inside the phase boundary to be much less than it is in the symmetric phase, then there are strict upper bounds on masses. If the MSSM is responsible for baryogenesis the prospects for discovering electric dipole moments and supersymmetric particles in the next few years are excellent. Should we fail to make these discoveries, there are several possible conclusions. The BAU may come from physics above the weak scale, from weak scale physics other than supersymmetry, or from a more complicated supersymmetric model. The simplest extension of the MSSM would be to add a gauge singlet superfield, which changes the allowed form of the Higgs potential and greatly relaxes all the constraints on scalar masses \[\]. The tree level potential will generically contain large cubic terms, and gives a very strongly first order transition with a thin boundary between the two phases. Furthermore, with gauge singlets there is the possibility for CP violating phases in the Higgs potential. With thin bubble walls we expect baryogenesis to be dominated by the mechanism of refs. , in which CP violating particle scattering processes from the phase boundary leads to a transport of particle quantum numbers, biasing anomalous baryon production throughout the symmetric phase. The charge transport mechanism could produce the BAU for phases as small as $10^{-4}$ . Unfortunately such extended supersymmetric models currently have too many free parameters to allow for calculation of the BAU.

In summary, we have shown that unlike the minimal standard model, the MSSM is still viable, but soon either it should be ruled out (for baryogenesis), or new experimental discoveries such as electric dipole moments will give us an important clue towards the explanation of why there are more baryons than anti-baryons. We are still a long way from testable predictions from baryogenesis in extensions of the MSSM.

**Acknowledgments**

We would like to thank David Kaplan and Stanley Myint for useful discussions. A. C. was supported in part by DOE contracts \#DE-AC02-89ER40509, \#DE-FG02-91ER40676, by the Texas National Research Laboratory Commission grant \#RGFY91B6, and NSF contract \#PHY-9057173; A. N. was supported in part by DOE contract \#DE-FG03-90ER40546, by a fellowship from the Alfred P. Sloan Foundation, and by an SSC Fellowship from the Texas National Research Laboratory Commission.

[^1]: 

[^2]: $^a$

[^3]: $^b$
