text:
# Cosmological constraints on a Peccei-Quinn flatino as the lightest supersymmetric particle

## I Introduction

In supersymmetric extensions of the Standard Model, both the $`\mu `$ problem and the strong CP problem can be resolved through a simple extension of the Higgs sector, which implements spontaneously broken Peccei-Quinn (PQ) symmetry pq . Suppose that the Higgs bilinear term in the minimal supersymmetric standard model (MSSM) superpotential takes the non-renormalizable form;

$$W_{\mu \text{-term}}=h\frac{\underset{i}{}\varphi _i^{l_i}}{M_P^{p1}}H_1H_2$$
(1)
where $`M_P2.4\times 10^{18}`$ GeV is the reduced Planck mass, $`p=_il_i2`$, and $`W_{\mu \text{-term}}`$ respects a PQ symmetry. Then, nonzero vacuum expectations values (vevs) $`\varphi _i`$ can spontaneously break PQ symmetry, providing the axion solution to the strong CP problem areview , while generating the right order of magnitude for the $`\mu `$ parameter muprob ;

$$\mu =h\frac{\underset{i}{}\varphi _i^{l_i}}{M_P^{p1}}.$$
(2)
In the case of a polynomial of the second order ($`p=2`$), with $`h1`$, this gives $`\mu F_\mathrm{a}^2/M_P100`$ GeV if the axion scale $`F_\mathrm{a}\varphi _i`$ is of order $`10^{10}`$ GeV, its lowest feasible value. With a higher power $`p>2`$, one would need a larger value of $`F_\mathrm{a}`$ to get the right size of $`\mu `$.

From now on, we consider the case of supergravity with gravity-mediated supersymmetry breaking nilles . There are two ways of spontaneously breaking PQ symmetry. One is to introduce a renormalizable tree-level superpotential. In this case, the particles corresponding to the fields $`\varphi _i`$, and their spin $`1/2`$ partners, all have mass of order $`F_\mathrm{a}`$ except for a supermultiplet comprising the axion (the pseudo-Goldstone boson of the PQ symmetry), the spin zero saxion, and the spin $`1/2`$ axino. The mass of the axino in this case is generically of order $`100`$ GeV, but it can be as small as O(keV) goto ; andre . If it is the lightest supersymmetric particle (LSP), there follow some interesting consequences in cosmology and collider physics. If the axino mass is $`10`$ GeV, it can be the cold dark matter of the universe if its (non-thermal) production comes mostly from the decay of the next-to-lightest supersymmetric particle (NLSP) covi . For this to happen, of course, the primordial LSP relic density has to be diluted away, which would be the case when the reheat temperature after inflation is sufficiently low chun2 . (In contrast, the axino mass is required to be less than a few keV if its primordial abundance is not washed out raja .)

The alternative to a renormalizable superpotential is a non-renormalizable superpotential with soft supersymmetry breaking, leading to a flaton model of of PQ symmetry breaking mox denis . (A flaton yam85 ; lyth  is a scalar field with vev much bigger than its mass, corresponding to a very flat potential.) For simplicity, the superpotential is supposed to be dominated by a single term,

$$W_{PQ}=\frac{f}{M_P^{n3}}\underset{i}{}\varphi _i^{k_i}$$
(3)
where $`n=_ik_i4`$. If the soft mass-squared of one or more of the singlet fields $`\varphi _i`$ is negative, e.g., through a radiative mechanism msy , there exists always a nontrivial global minimum. The $`\varphi _i`$ with nonzero vevs are then flatons. The axion scale is $`F_\mathrm{a}\varphi _i(m_0M_P^{n3})^{1/n2}`$, which is naturally of the right size. For instance, for $`n=4`$, $`f1`$ and $`m_0100`$ GeV, the vacuum expectation value is of order $`10^{10}`$ GeV, but it can be bigger if $`n`$ is bigger.

In a flaton model, the masses of the particles corresponding to the fields $`\varphi _i`$, and their spin $`1/2`$ partners, all have mass of order $`100`$ GeV, with the sole exception of the axion. The saxion and axino, defined as the superpartners of the axion, have no special significance and are in general not mass eigenstates. Instead, with $`N2`$ fields breaking PQ symmetry, there are $`2N1`$ spin zero flaton particles, and $`N`$ spin $`1/2`$ flatinos. These play the same role as, respectively, the saxion and the axino of renormalizable models. In particular, the lightest flatino may be the LSP and hence a dark matter candidate.

With either an axino or a flatino LSP, the NLSP, which is typically the lightest neutralino in the MSSM, decays to the LSP with the rate proportional to $`1/F_\mathrm{a}^2`$. An interesting implication for collider experiments has been pointed out in Ref. martin . Although the NLSP decays are very weak, and thus the corresponding decay lengths can be much larger than a collider size, there is an opportunity to observe these decays in future colliders with a large number of supersymmetric events so as to provide direct information about physics at very high energy scale.

In this paper, we study the cosmology of models with a flatino LSP, identifying various bounds on the parameters which are required for a viable cosmology. The paper may regarded as a sequel to denis , in which all of the flatinos were supposed to be unstable. We shall focus particularly on the case of a low axion scale, $`F_\mathrm{a}10^{10}`$ GeV, which increases the possibility of collider signatures.

## II Cosmology of flatino LSP

Flaton models lead to very interesting cosmological effects. On the reasonable assumption that the flaton fields are nonzero in the early Universe, they lead to thermal inflation which can dilute away all the unwanted relics lyth ; laza2 . In a flaton model of PQ symmetry breaking, some restrictions have to be put on the model in order not to overproduce unwanted relics again from flaton decay. As a (scalar) flaton field $`\phi `$ has always the decay channel $`\phi aa`$ with a rate $`\mathrm{\Gamma }_{\phi aa}m_\phi ^3/32\pi F_\mathrm{a}^2`$ andre , flaton decay may produce too many (unthermalized) axions and upset standard nucleosynthesis. In the scenario of the flatino LSP under consideration, it is also required that the LSP is not overproduced from flaton decay. In general, a flaton $`\phi `$ can decay into ordinary particles and their superpartners $`X`$, the axion $`a`$ and the flatino LSP $`\stackrel{~}{F}_1`$. After flaton decay, the energy densities are given by $`\rho _X=B_X\rho _\phi =(\pi ^2/30)g_{RH}T_{\mathrm{RH}}^4`$, $`\rho _a=B_a\rho _\phi `$ and $`\rho _{\stackrel{~}{F}_1}=B_{\stackrel{~}{F}_1}(m_{\stackrel{~}{F}_1}/m_\phi )\rho _\phi `$ where $`\rho _\phi `$ is the energy density of the flatons before they decay and the $`B`$’s are essentially branching ratios of the flaton decay. The precise relation between $`B_a`$ and the flaton branching ratios is given in denis . As is well-known, nucleosynthesis (NS) puts a constraint on the extra amount of relativistic energy density, which is conveniently given in terms of the equivalently number of neutrino species, $`\delta N_\nu `$. At present nsbound , the bound is something like $`\delta N_\nu <0.3`$ in the favored ‘low deuterium’ scenario. This can be translated to upper limits on $`B_a`$ and $`B_{\stackrel{~}{F}_1}`$. Applying $`(\rho _a/\rho _\nu )_{NS}<\delta N_\nu `$ for the energy densities at the time of nucleosynthesis, one finds choi ; denis

$$B_a<\frac{7}{43}\left(\frac{g_{RH}}{43/4}\right)^{1/3}\delta N_\nu .$$
(4)
In the similar way, the energy density of the stable and massive LSP $`\stackrel{~}{F}_1`$ has to be constrained; $`(\rho _{\stackrel{~}{F}_1}/\rho _\nu )_{NS}<\delta N_\nu `$, which leads to

$$B_{\stackrel{~}{F}_1}<\frac{7/4}{g_{RH}}\frac{m_\phi }{m_{\stackrel{~}{F}_1}}\frac{T_{NS}}{T_{\mathrm{RH}}}\delta N_\nu .$$
(5)
To get a rough estimate of the reheat temperature, let us take the rate of the decay $`\phi X`$ to be $`\mathrm{\Gamma }_X\mathrm{\Gamma }_a/B_a`$ resulting in

$$T_{\mathrm{RH}}1.2g_{RH}^{1/4}\sqrt{M_P\mathrm{\Gamma }_X}19\mathrm{GeV}\left(\frac{m_\phi }{100\mathrm{GeV}}\right)^{3/2}\left(\frac{10^{10}\mathrm{GeV}}{F_\mathrm{a}}\right)$$
(6)
where $`g_{RH}100`$ is the effective number of particle degrees at $`T_{\mathrm{RH}}`$ and we took the numerical factor $`B_a=0.1`$.

For $`T_{\mathrm{RH}}10`$ GeV, we get $`B_{\stackrel{~}{F}_1}10^4`$ with $`B_a0.1`$. A more stringent limit on $`B_{\stackrel{~}{F}_1}`$ comes from the flatino contribution to the present energy density. That is, in order to avoid the overclosure by the LSP, $`(\rho _{\stackrel{~}{F}_1}/\rho _c)_0<1`$, one gets

$$B_{\stackrel{~}{F}_1}\frac{4}{3}\frac{m_\phi }{m_{\stackrel{~}{F}_1}}\left(\frac{3.55\mathrm{eV}}{T_{\mathrm{RH}}}\right)B_X.$$
(7)
Taking $`T_{\mathrm{RH}}10`$ GeV, one obtains $`B_{\stackrel{~}{F}_1}10^{10}`$. This is a very strong constraint. It appears impossible to get such a small number without relying on a severe fine-tuning of parameters in the model once flaton decay into flatino is allowed kinematically. Therefore, we have to forbid the decay of flatons into the flatino LSP imposing

$$m_\phi <2m_{\stackrel{~}{F}_1}.$$
(8)
Of course, the condition (7) can be invalidated in the case that the flatino is heavy enough to decay into ordinary particles, e.g., a neutralino and a light Higgs boson denis  which is opposite to our consideration.

Forbidding the decay $`\phi \stackrel{~}{F}_1`$ kinematically, we now have to fulfill the condition (4). Important thermalizable decay modes of a flaton to achieve $`B_aB_X1`$ include the decays into two top quarks, two light stops, two light Higgs bosons $`h^0`$, two $`Z`$/$`W`$ gauge bosons and a Z boson plus a Higgs boson. In this paper, we take the flatons to be as light as possible, so as to allow the lightest possible masses for the flatino LSP and for the NLSP. We therefore assume that only the last two decay modes are kinematically allowed. The decay mode into Higgs bosons comes from direct couplings of Higgs bosons and flatons given by Eq. (1) and has been considered in Ref. denis  in the context of an ordinary neutralino LSP. The other modes listed above come from mixing between Higgs bosons and flatons induced by the same $`\mu `$ term interaction (1), which is the bosonic counterparts of mixing between neutralinos and flatinos which has been worked out in Ref. martin . (The qualitative features of the decay modes into top quarks or stops have also been considered in Ref. choi .)

Under the condition of no CP violation, scalar flatons (denoted by $`F`$) can decay into two gauge bosons, e.g., $`FWW`$, and pseudoscalars (denoted by $`F^{}`$) into a Z boson and a scalar Higgs, $`F^{}h^0Z`$, as we describe below. The former mode gives a lower limit on the scalar flaton masses, $`2m_W<m_F`$. Combining it with Eq. (8) we find

$$2m_W<m_F<2m_{\stackrel{~}{F}_1}<2m_{\stackrel{~}{\chi }_1^0}$$
(9)
where $`\stackrel{~}{\chi }_1^0`$ is the NLSP neutralino. (When $`m_F<2m_W`$ for instance, flatons can of course decay into three ordinary light fermions mediating a sfermion. But the corresponding decay rates have a large phase space suppression which makes it hard to dominate over the decay mode $`Faa`$.) Concerning pseudoscalar flatons, there are more possibilities. When $`m_F^{}>m_F`$, the decay mode $`F^{}aF`$ is open and has to be suppressed against the mode $`F^{}h^0Z`$. This requires

$$m_{h^0}+m_Z<m_F^{}<2m_{\stackrel{~}{F}_1}<2m_{\stackrel{~}{\chi }_1^0}.$$
(10)
On the other hand, if $`m_F^{}<m_F`$, $`F^{}`$ has no two body decay mode into an axion, and thus $`F^{}`$ can be so light as to have the decay modes into light fermions; $`F^{}f\overline{f}(Z)`$ or $`f\overline{f}h^0`$. In this case, the corresponding decay rate has a large Yukawa-coupling or phase-space suppression factor leading to lower reheat temperature than in Eq. (6), and only the bound (9) is applied.

Without resorting to a specific model, let us describe how the flaton couplings to gauge bosons arise. If we assume no CP violation, scalar flatons mix with the two CP-even Higgs bosons $`h^0,H^0`$ and pseudoscalar flatons mix with the CP-odd Higgs boson $`A`$, through the $`\mu `$-term potential as in Eq. (1). Quantifying small mixtures of flatons in Higgs bosons by $`\epsilon _{h^0F}`$, $`\epsilon _{H^0F}`$ and $`\epsilon _{AF^{}}`$ in a self-explaining notation, it is straightforward to write down the flaton-Higgs couplings from the MSSM Lagrangian:

$``$ $`=`$ $`gm_W\epsilon _{VVF}[FW_\mu ^+W_{}^\mu +{\displaystyle \frac{1}{c_W^2}}FZ_\mu Z^\mu ]`$ (11)
$`{\displaystyle \frac{g}{c_W}}\epsilon _{AF^{}}[\mathrm{cos}(\alpha \beta )h^0+\mathrm{sin}(\alpha \beta )H^0]\overline{_\mu }F^{}Z^\mu `$
where $`\epsilon _{VVF}\mathrm{sin}(\alpha \beta )\epsilon _{h^0F}\mathrm{cos}(\alpha \beta )\epsilon _{H^0F}`$, $`\alpha `$ is the diagonalization angle of CP-even Higgses, $`\beta `$ is defined by $`\mathrm{tan}\beta =H_2/H_1`$ and $`c_W=\mathrm{cos}\theta _W`$ with $`\theta _W`$ being the weak mixing angle.

Let us now consider the flatino–neutralino mixing martin  which arises at tree level due to the $`\mu `$ term (1). For typical parameters, this mixing will give the main contribution to the flatino interactions with ordinary particles and their superpartners, dominating the one-loop processes covi  which are otherwise responsible for these interactions. The relevant interactions for our purpose are the decay of the NLSP to the flatino LSP, and the decays of sfermions to the flatino LSP.

The term (1) leads to the mass matrix that mixes $`\stackrel{~}{F}_i`$ and ($`\stackrel{~}{B}`$, $`\stackrel{~}{W}_3`$, $`\stackrel{~}{H}_1`$, $`\stackrel{~}{H}_2`$) as follows;

$$\left(\begin{array}{cccc}0& 0& \mu s_\beta \delta _i& \mu c_\beta \delta _i\end{array}\right),$$
(12)
where $`c_\beta \mathrm{cos}\beta `$, $`s_\beta \mathrm{sin}\beta `$, and the coefficients $`\delta _iv/F_\mathrm{a}`$ depend on the specific form of the superpotential (1). Here $`v=264`$ GeV is the Higgs vev. Let $`N`$ be the diagonalization matrix of the MSSM neutralino mass matrix $`M_N`$: $`\stackrel{~}{M}_N=N^TM_NN`$ where $`\stackrel{~}{M}_N`$ is the diagonalized neutralino mass matrix for the eigenstate $`\stackrel{~}{\chi }_j^0`$. Then, further diagonalizing (12) can be done by the small mixing elements;

$$\epsilon _{\stackrel{~}{\chi }_j^0\stackrel{~}{F}_i}=\delta _i\frac{\mu (s_\beta N_{3j}+c_\beta N_{4j})}{m_{\stackrel{~}{\chi }_j^0}m_{\stackrel{~}{F}_i}}$$
(13)
which describe the small content of the flatino $`\stackrel{~}{F}_i`$ in the neutralino $`\stackrel{~}{\chi }_j^0`$. Note that we can have an enhancement in the mixing (13), $`\epsilon _{\stackrel{~}{\chi }_j^0\stackrel{~}{F}_i}10v/F_a`$, when $`m_{\stackrel{~}{\chi }_j^0}m_{\stackrel{~}{F}_i}=O(10)`$ GeV.

Flaton-Higgs mixing parameters and flaton/flatino masses are dependent on specific forms of the terms (1) and (3). As we discussed, in order for the flatino LSP scenario to be consistent with cosmological considerations, the model should fulfill the mass relations (8), (9) and (10), and have sufficiently large mixing elements $`\epsilon `$ in Eq. (11) to suppress the axion decay modes. However, note that mixing between the flaton sector and the MSSM sector is determined by the term (1) and flaton/flatino mass spectrum by the term (3) and its soft-breaking term. Therefore, arrangement for the mass relations (8,9,10) and large mixing elements $`\epsilon `$ can be done by independent parameters in two separate sectors, which implies that such an arrangement can be easily achieved in generic flaton models.

As the flatino LSP has a mass in the 100 GeV region, its slight regeneration after the thermal inflation can provide a sizable contribution to the matter density at present. There are two important sources of the flatino regeneration. The first is the neutralino decay covi  after its decoupling from the thermal bath. In the flaton models, the NLSP neutralino can decay into $`Z(h^0)\stackrel{~}{F}_1`$ or $`f\overline{f}\stackrel{~}{F}_1`$ within 0.01 sec martin  and thus without affecting nucleosynthesis. Then, the ratio of the flatino LSP density to the critical density is given by

$$\mathrm{\Omega }_{\stackrel{~}{F}_1}=\frac{m_{\stackrel{~}{F}_1}}{m_{\stackrel{~}{\chi }_1^0}}\mathrm{\Omega }_{\stackrel{~}{\chi }_1^0}.$$
(14)
For our case, $`m_{\stackrel{~}{F}_1}/m_{\stackrel{~}{\chi }_1^0}=𝒪(1)`$ and thus we can have $`\mathrm{\Omega }_{\stackrel{~}{F}_1}\mathrm{\Omega }_{\stackrel{~}{\chi }_1^0}0.11`$ which is valid in a wide range of the MSSM parameter space drees . For this mechanism to work, the reheat temperature (6) after thermal inflation should be larger than the neutralino decoupling temperature $`m_{\stackrel{~}{\chi }_1^0}/20`$. As can be seen from Eq. (6), this condition holds for low $`F_\mathrm{a}10^{10}10^{11}`$ GeV which increases collider signals as well.

A potentially more important source is thermal regeneration covi2 . Since our reheat temperature is typically below 100 GeV, decay processes dominate over scattering processes for the thermal regeneration. To get a qualitative calculation of the flatino population, let us take a sfermion–fermion–flatino interaction arising from the neutralino–flatino mixing,

$$=g\epsilon _{\stackrel{~}{f}f\stackrel{~}{F}_1}\stackrel{~}{f}f\stackrel{~}{F}_1+h.c.$$
(15)
where $`\epsilon _{\stackrel{~}{f}f\stackrel{~}{F}_1}v/F_\mathrm{a}`$ contains the factors from gauge quantum numbers and neutralino–flatino mixing and $`\stackrel{~}{f}`$/$`f`$ denotes a left-handed or right-handed sfermion/fermion field. The thermally regenerated flatino population is given by $`\mathrm{\Omega }_{\stackrel{~}{F}_1}h^2=m_{\stackrel{~}{F}_1}Y_{\stackrel{~}{F}_1}/3.55\mathrm{eV}`$ with the factor covi2 ; hbkim  :

$$Y_{\stackrel{~}{F}_1}2\times 10^5\frac{M_P\mathrm{\Gamma }_{\stackrel{~}{f}}}{T_{\mathrm{RH}}^2}F(x_{\stackrel{~}{f}})$$
(16)
where $`\mathrm{\Gamma }_{\stackrel{~}{f}}=g^2\epsilon _{\stackrel{~}{f}f\stackrel{~}{F}_1}^2m_{\stackrel{~}{f}}/8\pi `$ is the decay rate and $`F(x_{\stackrel{~}{f}})`$ is the Boltzmann suppression factor as a function of $`x_{\stackrel{~}{f}}=m_{\stackrel{~}{f}}/T_{\mathrm{RH}}`$ given by

$$F(x_{\stackrel{~}{f}})\frac{1}{x_{\stackrel{~}{f}}^2}_{x_{\stackrel{~}{f}}}^{\mathrm{}}\frac{dx}{e^x1}\left[\left(\frac{\pi }{2}\mathrm{tan}^1\frac{x}{\sqrt{x^2x_{\stackrel{~}{f}}^2}}\right)x^4+x_{\stackrel{~}{f}}\left(x^22x_{\stackrel{~}{f}}\sqrt{x^2x_{\stackrel{~}{f}}^2}\right)\right].$$

We have defined the function $`F(x_{\stackrel{~}{f}})`$ to remove the dependence on the axion scale $`F_\mathrm{a}`$ for the remaining factor; $`M_P\mathrm{\Gamma }_{\stackrel{~}{f}}/T_{\mathrm{RH}}^2\mathrm{\Gamma }_{\stackrel{~}{f}}/\mathrm{\Gamma }_\phi `$ where $`\mathrm{\Gamma }_\phi `$ is the flaton decay rate. To a very good approximation for the range of $`x5`$, the function $`F(x)`$ can be expressed in terms of an analytic function; $`F(x)36e^{0.98x}`$ giving

$$\mathrm{\Omega }_{\stackrel{~}{F}_1}h^25\times 10^4\left(\frac{m_{\stackrel{~}{F}_1}}{m_{\stackrel{~}{f}}}\right)\left(\frac{\epsilon _{\stackrel{~}{f}f\stackrel{~}{F}_1}}{10^8}\right)^2x_{\stackrel{~}{f}}^2e^{0.98x_{\stackrel{~}{f}}}.$$
(17)
From the above equation, one finds that the requirement $`\mathrm{\Omega }_{\stackrel{~}{F}_1}h^21`$ gives $`x_{\stackrel{~}{f}}17`$ with $`m_{\stackrel{~}{F}_1}=m_{\stackrel{~}{f}}`$ and $`\epsilon _{\stackrel{~}{f}f\stackrel{~}{F}_1}=10^8`$. This consideration gives some meaningful bounds on sfermions masses, $`m_{\stackrel{~}{f}}>17T_{\mathrm{RH}}`$ given the reheat temperature (6), or vice versa. The thermally regenerated population of the flatino LSP can also be the cold dark matter which, however, needs a fine arrangement for the sfermion mass and reheat temperature as $`\mathrm{\Omega }_{\stackrel{~}{F}_1}h^2`$ is a very sensitive function of $`x_{\stackrel{~}{f}}`$. Note that the thermal regeneration becomes easily negligible as $`F_\mathrm{a}`$ becomes large. For instance, for $`F_\mathrm{a}=10^{11}`$ GeV, one has $`T_{\mathrm{RH}}2`$ GeV (6) and thus $`x_{\stackrel{~}{f}}50`$ for $`m_{\stackrel{~}{f}}>100`$ GeV.

## III A Model

To illustrate our discussion more explicitly, we consider the minimal case of two fields and $`n=4`$, which has been analyzed in Ref. denis . This model contains two scalar flatons, one pseudoscalar flaton, and two flatinos. The flaton superpotential is

$$W_{PQ}=\frac{f}{M_P}P^3Q$$
(18)
with $`U(1)_{PQ}`$ charges of $`P`$ and $`Q`$ being $`1`$ and $`3`$, respectively. Analyzing the scalar potential including soft supersymmetry breaking terms

$$V_{soft}=\frac{fA_f}{M_P}P^3Q+h.c.,$$
(19)
one finds scalar flatons $`F_{2,1}`$ with masses-squared

$$m_{F_2,F_1}^2=\frac{1}{2}f^2\mu _0^2\left(3(12\xi )+x^2(12+\xi )\pm |12\xi |\sqrt{x^4+42x^2+9}\right)$$
(20)
and a pseudoscalar flaton $`F^{}`$ with mass-squared

$$m_F^{}^2=f^2\mu _0^2\xi (x^2+9)$$
(21)
where the parameters are defined by

$$x\frac{P}{Q},\xi \frac{A_f}{f\mu _0}>0,\mu _0\frac{PQ}{2M_P}.$$

In the above, the flatons $`P,Q`$ are expanded as

$`P`$ $`=`$ $`{\displaystyle \frac{1}{\sqrt{2}}}\left(P+P^{}3i{\displaystyle \frac{Q}{F_\mathrm{a}}}F^{}\right)`$
$`Q`$ $`=`$ $`{\displaystyle \frac{1}{\sqrt{2}}}\left(Q+Q^{}i{\displaystyle \frac{P}{F_\mathrm{a}}}F^{}\right)`$ (22)
where $`F_\mathrm{a}^2P^2+9Q^2`$. The diagonalization matrix of the scalar flatons $`P^{}`$ and $`Q^{}`$ is

$$\left(\begin{array}{c}P^{}\\ Q^{}\end{array}\right)=\left(\begin{array}{cc}\mathrm{cos}\phi & \mathrm{sin}\phi \\ \mathrm{sin}\phi & \mathrm{cos}\phi \end{array}\right)\left(\begin{array}{c}F_2\\ F_1\end{array}\right)$$
(23)
where the mixing angle $`\phi `$ is determined by

$`\mathrm{cos}2\phi `$ $`=`$ $`sgn(12\xi ){\displaystyle \frac{x^2+3}{\sqrt{x^4+42x^2+9}}}`$
$`\mathrm{sin}2\phi `$ $`=`$ $`sgn(12\xi ){\displaystyle \frac{6x}{\sqrt{x^4+42x^2+9}}}.`$ (24)
From the superpotential (18), one also finds the flatino masses

$$m_{\stackrel{~}{F}_2,\stackrel{~}{F}_1}=3f\mu _0[\sqrt{x^2+1}\pm 1]$$
(25)
The light flatino $`\stackrel{~}{F}_1`$ is supposed to be the LSP. The rotation matrix from the flavor states $`\stackrel{~}{P},\stackrel{~}{Q}`$ to the mass eigenstates $`\stackrel{~}{F}_2,\stackrel{~}{F}_1`$ is given by

$$\left(\begin{array}{c}\stackrel{~}{P}\\ \stackrel{~}{Q}\end{array}\right)=\left(\begin{array}{cc}\mathrm{cos}\stackrel{~}{\phi }& \mathrm{sin}\stackrel{~}{\phi }\\ \mathrm{sin}\stackrel{~}{\phi }& \mathrm{cos}\stackrel{~}{\phi }\end{array}\right)\left(\begin{array}{c}\stackrel{~}{F}_1\\ \stackrel{~}{F}_2\end{array}\right)$$
(26)
with the mixing angle $`\stackrel{~}{\phi }`$ satisfying

$$\mathrm{cos}2\stackrel{~}{\phi }=\frac{1}{\sqrt{x^2+1}},\mathrm{sin}2\stackrel{~}{\phi }=\frac{x}{\sqrt{x^2+1}}.$$
(27)
As all the masses contain the overall factor $`|f\mu _0|`$, it is now straightforward to find the region of the parameters $`x`$ and $`\xi `$ satisfying the desired mass relations. It turns out that $`m_{F_{1,2}},m_F^{}<2m_{\stackrel{~}{F}_1}`$ requires

$`(A)`$ $`\xi <12;x>2.6,\xi _{a,b}<\xi <\xi _c`$
$`(B)`$ $`\xi >12;x>3.5,12<\xi <\xi _c`$ (28)
$`\mathrm{with}`$ $`\xi _a=12{\displaystyle \frac{\sqrt{x^4+42x^2+9}x^23}{\sqrt{x^4+42x^2+9}+x^23}}`$
$`\xi _b=12{\displaystyle \frac{\sqrt{x^4+42x^2+9}+12\sqrt{x^2+1}5x^29}{\sqrt{x^4+42x^2+9}+x^23}}`$
$`\xi _c=36{\displaystyle \frac{x^2+22\sqrt{x^2+1}}{x^2+9}}`$
We always have $`m_{F_1}<m_F^{},m_{F_2}`$ in the regions (A,B) and also $`m_{F_2}<m_F^{}`$ in the region (B). The requirement (10) has to be fulfilled in both regions.

The important interactions of flatons (flatinos) and MSSM fields comes from the mixings between flatons (flatinos) and Higgs bosons (neutralinos) due to the $`\mu `$-term superpotential of our choice

$$W_{\mu \text{-term}}=h\frac{PQ}{M_P}H_1H_2$$
(29)
with $`U(1)_{PQ}`$ charge of $`H_1H_2`$ being $`+2`$. Note that $`\mu =hPQ/2M_P=h\mu _0`$. First of all, Eq. (29) determines the flatino–Higgsino mixing elements $`\delta _i`$ defined in Eq. (12) as follows:

$$\delta _1=\frac{v}{F_\mathrm{a}}\frac{\sqrt{x^2+9}}{x}(c_{\stackrel{~}{\phi }}+xs_{\stackrel{~}{\phi }}),\delta _2=\frac{v}{F_\mathrm{a}}\frac{\sqrt{x^2+9}}{x}(s_{\stackrel{~}{\phi }}+xc_{\stackrel{~}{\phi }}),$$
(30)
with $`c_{\stackrel{~}{\phi }}\mathrm{cos}_{\stackrel{~}{\phi }}`$ and $`s_{\stackrel{~}{\phi }}\mathrm{sin}_{\stackrel{~}{\phi }}`$. Including the soft term of (29), we get the Higgs-flaton scalar potential,

$`V`$ $`=`$ $`\left|H_1\right|^2\left(m_{H_1}^2+\left|h{\displaystyle \frac{PQ}{M_P}}\right|^2\right)+\left|H_2\right|^2\left(m_{H_2}^2+\left|h{\displaystyle \frac{PQ}{M_P}}\right|^2\right)`$
$`+\left\{hH_1H_2\left(A_h{\displaystyle \frac{PQ}{M_P}}+3f^{}{\displaystyle \frac{P^2|Q|^2}{M_P^2}}+f^{}{\displaystyle \frac{P^2|P|^2}{M_P^2}}\right)+\text{c.c.}\right\}`$
$`+{\displaystyle \frac{1}{8}}(g^2+g_{}^{}{}_{}{}^{2})\left(|H_1|^2|H_2|^2\right)^2.`$
We note that the CP-odd Higgs boson has mass-squared

$$m_A^2=\frac{2}{\mathrm{sin}2\beta }\mu ^2\left(\frac{A_h}{\mu }+\frac{f}{h}(x^2+3)\right).$$
(32)
It is then straightforward to find the following flaton–Higgs mixing:

$`\epsilon _{h^0F_1}`$ $`=`$ $`{\displaystyle \frac{v}{F_\mathrm{a}}}{\displaystyle \frac{\sqrt{x^2+9}}{x}}{\displaystyle \frac{\mu ^2}{m_{h^0}^2m_{F_1}^2}}\{\mathrm{sin}(\alpha \beta )[2s_\phi 2xc_\phi ]`$
$`+\mathrm{cos}(\alpha \beta )[{\displaystyle \frac{f}{h}}(4x^2s_\phi +6s_\phi 6xc_\phi )+{\displaystyle \frac{A_h}{\mu }}(s_\phi xc_\phi )]\}`$
$`\epsilon _{H^0F_1}`$ $`=`$ $`{\displaystyle \frac{v}{F_\mathrm{a}}}{\displaystyle \frac{\sqrt{x^2+9}}{x}}{\displaystyle \frac{\mu ^2}{m_{H^0}^2m_{F_1}^2}}\{\mathrm{cos}(\alpha +\beta )[2s_\phi 2xc_\phi ]`$
$`+\mathrm{sin}(\alpha +\beta )[{\displaystyle \frac{f}{h}}(4x^2s_\phi +6s_\phi 6xc_\phi )+{\displaystyle \frac{A_h}{\mu }}(s_\phi xc_\phi )]\}`$
$`\epsilon _{h^0F_2}`$ $`=`$ $`{\displaystyle \frac{v}{F_\mathrm{a}}}{\displaystyle \frac{\sqrt{x^2+9}}{x}}{\displaystyle \frac{\mu ^2}{m_{h^0}^2m_{F_2}^2}}\{\mathrm{sin}(\alpha \beta )[2c_\phi +2xs_\phi ]`$
$`\mathrm{cos}(\alpha \beta )[{\displaystyle \frac{f}{h}}(4x^2c_\phi +6c_\phi +6xs_\phi )+{\displaystyle \frac{A_h}{\mu }}(c_\phi +xs_\phi )]\}`$
$`\epsilon _{H^0F_2}`$ $`=`$ $`{\displaystyle \frac{v}{F_\mathrm{a}}}{\displaystyle \frac{\sqrt{x^2+9}}{x}}{\displaystyle \frac{\mu ^2}{m_{H^0}^2m_{F_2}^2}}\{\mathrm{cos}(\alpha +\beta )[2c_\phi +2xs_\phi ]`$
$`\mathrm{sin}(\alpha +\beta )[{\displaystyle \frac{f}{h}}(4x^2c_\phi +6c_\phi +6xs_\phi )+{\displaystyle \frac{A_h}{\mu }}(c_\phi +xs_\phi )]\}`$
$`\epsilon _{AF^{}}`$ $`=`$ $`{\displaystyle \frac{v}{F_\mathrm{a}}}{\displaystyle \frac{x^2+3}{x}}{\displaystyle \frac{\mu ^2}{m_A^2m_F^{}^2}}\{{\displaystyle \frac{A_h}{\mu }}6{\displaystyle \frac{f}{h}}\}`$ (37)
where $`v\sqrt{H_1^2+H_2^2}`$. The flaton decay rates for the processes $`FWW`$ and $`F^{}h^0Z`$ of interest are given by

$`\mathrm{\Gamma }(F_iWW)`$ $`=`$ $`\alpha _2|\epsilon _{VVF_i}|^2{\displaystyle \frac{m_W^2}{m_{F_i}}}\sqrt{14{\displaystyle \frac{m_W^2}{m_{F_i}^2}}}`$ (38)
$`\mathrm{\Gamma }(F^{}h^0Z)`$ $``$ $`{\displaystyle \frac{\alpha _2}{8c_W^2}}\mathrm{cos}(\alpha \beta )^2|\epsilon _{AF^{}}|^2m_F^{}(1+{\displaystyle \frac{m_{h^0}^2\frac{1}{2}m_Z^2}{m_F^{}^2}})`$
where $`\epsilon _{VVF_i}\mathrm{sin}(\alpha \beta )\epsilon _{h^0F_i}\mathrm{cos}(\alpha \beta )\epsilon _{H^0F_i}`$. The above rates have to be compared with the rates for the decay processes $`Faa`$ and $`F^{}aF_i`$ \[restricting ourselves to the region (B) in Eq. (III)\] which are given by denis

$`\mathrm{\Gamma }(F_1aa)`$ $`=`$ $`{\displaystyle \frac{1}{32\pi }}{\displaystyle \frac{m_{F_1}^3}{F_\mathrm{a}^2}}{\displaystyle \frac{(xs_\phi +9c_\phi )^2}{(x^2+9)}}`$
$`\mathrm{\Gamma }(F_2aa)`$ $`=`$ $`{\displaystyle \frac{1}{32\pi }}{\displaystyle \frac{m_{F_2}^3}{F_\mathrm{a}^2}}{\displaystyle \frac{(xc_\phi +9s_\phi )^2}{(x^2+9)}}`$ (39)
$`\mathrm{\Gamma }(F^{}aF_i)`$ $`=`$ $`{\displaystyle \frac{1}{16\pi }}{\displaystyle \frac{m_F^{}^3}{F_\mathrm{a}^2}}\left(1{\displaystyle \frac{m_{F_i}^2}{m_F^{}^2}}\right)^3{\displaystyle \frac{(3c_\phi 3xs_\phi )^2}{(x^2+9)}}.`$
Let us define $`R(F_i)\mathrm{\Gamma }(F_iaa)/\mathrm{\Gamma }(F_iWW)`$ and $`R(F^{})\mathrm{\Gamma }(F^{}aF_i)/\mathrm{\Gamma }(F^{}h^0Z)`$ which behave like

$$R(F_i)\frac{m_{F_i}^4}{v^2m_W^2}\frac{m_{F_i}^4}{\mu ^4},R(F^{})\frac{m_F^{}^2}{v^2}\frac{(m_A^2m_F^{}^2)^2}{\mu ^4}$$
(40)
apart from the other numerical factors. Eq. (40) shows that small $`R`$ can be easily obtained with small flaton masses and large $`\mu `$. In the case of small $`\mu `$, one would need to have small $`m_A`$ as well as $`m_F^{}<v`$.

Having calculated all of the relevant quantities, we could in principle identify the region of parameter space which gives a viable cosmology as well as viable particle physics. Here, we present instead two sets of parameters, with (i) a bino-like and (ii) Higgsino-like NLSP, which satisfy all the requirements. The sets give an NLSP mass around the lowest possible value $`(m_{h^0}+m_Z)/2`$. The flatino LSP is only a little lighter, while the other flatino as well as the flatons have masses around $`200\text{GeV}`$, features which appear to be typical for the region of parameter space corresponding to a light NLSP.

In making the choice of parameters, we focussed on the region $`2m_W<m_{F_1},m_{F_2}<2m_Z`$, which can be arranged for $`\xi 12`$ \[see Eq. (20)\]. Both parameter sets have

$$\mathrm{tan}\beta =3,m_{h^0}=110\mathrm{G}\mathrm{e}\mathrm{V},M_2=2M_1,$$

and the other parameters are listed in in Table I along with the important output quantities. In both cases, we have $`R(F_i,F^{})<O(0.01)`$ and the reheat temperature $`T_{\mathrm{RH}}=16`$ GeV using the decay rates (38). The neutralino NLSP decouples after reheating, and can decay into the flatino LSP to provide the cold dark matter of the universe. To check that potential overproduction of the flatino LSP from the decay of thermal sfermions can be avoided, let us take as an example the right-handed stau. From $`\epsilon _{\chi ^0\stackrel{~}{F}_1}`$ in Table I, the coupling in Eq. (17) is given by $`\epsilon _{\stackrel{~}{\tau }_R\tau _R\stackrel{~}{F}_1}=4.7\times 10^8`$ which gives $`\mathrm{\Omega }_{\stackrel{~}{F}_1}h^21`$ for $`m_{\stackrel{~}{\tau }_R}300`$ GeV for the case (i) and a similar figure is obtained for the case (ii).

## IV Conclusion

In an interesting class of models, a non-renormalizable term of the superpotential is responsible for the spontaneous breaking of PQ symmetry, while another such term generates the $`\mu `$ term of the MSSM. The flaton fields which break PQ symmetry are accompanied by flatinos, and the lightest flatino can be the LSP. Through the $`\mu `$ term of the superpotential, the NLSP neutralino decay might then be visible at colliders.

In this paper, we have examined the cosmology of this kind of model, following our earlier work denis  which explored the same flaton models on the assumption that all of the flatinos were unstable. We make the reasonable assumption that the flaton fields are nonzero in the early Universe, giving thermal inflation which eliminates pre-existing relics and leads to a rather well-defined cosmology. The decay of flatons into relativistic axions will interfere with nucleosynthesis unless the branching ratio is small, while the decay of flatons into flatinos should be forbidden altogether so as to avoid overproduction of the flatino LSP. We have argued that these combined requirements lead to the bounds $`m_{\stackrel{~}{\chi }_1^0}>m_{\stackrel{~}{F}_1}>m_W`$ or $`(m_Z+m_{h^0})/2`$ on the masses of the neutralino NLSP $`\stackrel{~}{\chi }_1^0`$ and the flatino LSP $`\stackrel{~}{F}_1`$. In order to have the lightest possible NLSP, one can focus on the case of light scalar (pseudoscalar) flatons, which can decay only into two $`W`$ bosons (a light Higgs and a $`Z`$ boson). These decays, as well as the neutralino NLSP decay into the flatino LSP, come from the mixing between flatons (flatinos) and Higgses (Higgsinos) due to the $`\mu `$ term of the superpotential.

We focus on the case of a low axion scale, $`F_\mathrm{a}10^{10}`$ GeV, which increases the feasibility of observing NLSP decays in future colliders. While the axion can hardly be the dark matter with such a low scale, the flatino LSP becomes a good dark matter candidate, because the reheat temperature can be high enough to generate the NLSP in thermal equilibrium, but low enough to suppress flatino production from the decay of thermally produced sfermions.

To make some definite statements, we studied an explicit model with $`F_\mathrm{a}10^{10}`$ GeV. We verified that flaton decays into $`WW`$ or $`h^0Z`$ can be large enough compared to their decays into axions in a wide range of parameter space, giving rise to the reheat temperature $`T_{\mathrm{RH}}O(10)`$ GeV. In this model, the NLSP can be as light as $`120\text{GeV}`$ with $`m_{h^0}=110`$ GeV, while the other flatino and the flatons have masses $`\text{ }>200\text{GeV}`$. We further notice that the parameter space admits the interesting mass region, $`m_{\stackrel{~}{\chi }_1^0}m_{\stackrel{~}{F}_1}=O(10)`$ GeV, where an enhancement of the flaton–neutralino mixing, $`\epsilon _{\stackrel{~}{\chi }_1^0\stackrel{~}{F}_1}10^7`$, results and thus the NLSP decays becomes faster.